 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1909 | 回復(fù): 5 | |||||
[求助]
關(guān)于Ewald sum算法的相關(guān)問題 已有1人參與
|
|
自己作死地用C編寫Monte Carlo的代碼，模擬的具體問題是粘土礦物對(duì)Ar的吸附。 然后需要加上Ewald Sum算法的問題算礦物羥基的長程力。 但是這個(gè)礦物是三斜晶系的，書上給的Ewald Sum是立方格子。 我查過LAMMPS里面，前幾年有個(gè)Ewald/n方法專門算三斜格子，現(xiàn)在和Ewald本身合并了。 我決定用三斜的格子算，因?yàn)榈V物的周期性是三斜的，用立方重復(fù)不了。但是實(shí)驗(yàn)室其他人不搞算法，他們說關(guān)于三斜粘土礦物的文獻(xiàn)，沒有聽過把格子劃分為三斜格子的說法。 我認(rèn)為算短程力的格子可以是立方的，但長程力不行，那些文獻(xiàn)里沒有寫清楚。 不知道我的觀點(diǎn)正不正確，因?yàn)槲抑荒苷业搅⒎礁褡拥脑创a和偽代碼，沒見過三斜的代碼。 還有，如果有編程簡單一點(diǎn)的，被大家比較公認(rèn)的算長程力的算法，求推薦。我已經(jīng)要被Ewald sum坑死了。 |
回復(fù)此樓» 猜你喜歡
312求調(diào)劑
已經(jīng)有7人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有16人回復(fù)
-
280求調(diào)劑
已經(jīng)有3人回復(fù)
-
266求調(diào)劑
已經(jīng)有7人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境工程調(diào)劑
已經(jīng)有8人回復(fù)
-
302求調(diào)劑
已經(jīng)有9人回復(fù)
-
297求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有4人回復(fù)
-
焦慮
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于pwscf進(jìn)行優(yōu)化的問題
已經(jīng)有5人回復(fù)
- gromacs運(yùn)行問題
已經(jīng)有9人回復(fù)
- LDA+U遇到的一個(gè)奇怪的問題？？？
已經(jīng)有14人回復(fù)
- [求助成功]VASP 計(jì)算 graphene 的錯(cuò)誤
已經(jīng)有15人回復(fù)
- 分子動(dòng)力學(xué)模擬長程庫侖勢(shì)的處理問題
已經(jīng)有8人回復(fù)
- 【分享】分子模擬中作圖和看軌跡的軟件vmd windows版下載及使用手冊(cè)
已經(jīng)有13人回復(fù)
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲號(hào): 556471
- 注冊(cè): 2008-05-11
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
Waxwane(月只藍(lán)代發(fā)): 金幣+1, 鼓勵(lì)交流 2015-01-11 12:37:23
感謝參與，應(yīng)助指數(shù) +1
Waxwane(月只藍(lán)代發(fā)): 金幣+1, 鼓勵(lì)交流 2015-01-11 12:37:23
|
Ewald的k空間求和里的傅立葉變換還是在直角坐標(biāo)系下進(jìn)行，理論上來講跟你用什么晶系關(guān)系不大。也就是說你依然可以把盒子里的電荷密度做一個(gè)傅立葉變換然后求和。 不過有一點(diǎn)技術(shù)上我沒太想明白的是在三斜晶系里做傅立葉變換，倒格子的格矢如何選取。立方盒子格矢很容易選取，直接2pi/L就可以了。三斜晶系下可能要做點(diǎn)額外處理來獲得一組正交的倒格子格矢？ 實(shí)空間部分求和跟你的晶系有一定關(guān)系，但這一部分很容易處理，因?yàn)橹皇莻�(gè)短程相互作用＋考慮minimum image convention。 |
2樓2015-01-11 03:35:12
3樓2015-01-11 18:31:08
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲號(hào): 556471
- 注冊(cè): 2008-05-11
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
|
看不懂你在說什么。 我又看了一下DL_POLY中Ewald sum的代碼，對(duì)于所有晶格，它們k空間部分的處理是這樣的。 1.先從晶格的格矢得到倒格子(reciprocal lattice）的格矢 2. 由倒格子的格矢再構(gòu)建一組正交的k空間矢量 3. 有了這組正交的k空間矢量，就可以獲得kmax并進(jìn)行下一步k空間的計(jì)算了。 對(duì)于步驟2，如果是立方盒子，這一步相當(dāng)于什么都沒做，因?yàn)樵鹊垢褡拥母袷妇褪钦坏�。但�?duì)于三斜的格子，因?yàn)槠涞垢褡痈袷覆皇钦坏模@一步就可以獲得一組正交的k空間矢量。步驟2的算法的大概意思是，對(duì)于三斜晶系，總可以找到一組正交的格矢，它們形成的晶胞體積跟原晶胞相等。 如果你想看dlpoly的代碼，可以在網(wǎng)上搜DL_POLY classic的代碼（下載需要注冊(cè)，但免費(fèi)）。Ewald的代碼在文件ewald_module.f里，參看subroutine ewald1. |
4樓2015-01-12 01:20:43
王慧1021
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1018.5
- 帖子: 140
- 在線: 144.6小時(shí)
- 蟲號(hào): 1535805
- 注冊(cè): 2011-12-13
- 性別: MM
- 專業(yè): 理論和計(jì)算化學(xué)
5樓2015-01-13 11:48:52
DDXZS
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 553.8
- 帖子: 164
- 在線: 26.3小時(shí)
- 蟲號(hào): 9850412
- 注冊(cè): 2018-08-12
- 性別: MM
- 專業(yè): 環(huán)境污染化學(xué)
6樓2020-05-25 20:55:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Waxwane 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +3 | 煙汐憶海 2026-03-16 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過，有科研經(jīng)歷 +8 | 曦熙兮 2026-03-15 | 8/400 |
|
|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123！ 2026-03-11 | 4/200 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|