 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1908 | 回復(fù): 5 | ||||
[求助]
關(guān)于Ewald sum算法的相關(guān)問(wèn)題 已有1人參與
|
|
自己作死地用C編寫Monte Carlo的代碼，模擬的具體問(wèn)題是粘土礦物對(duì)Ar的吸附。 然后需要加上Ewald Sum算法的問(wèn)題算礦物羥基的長(zhǎng)程力。 但是這個(gè)礦物是三斜晶系的，書上給的Ewald Sum是立方格子。 我查過(guò)LAMMPS里面，前幾年有個(gè)Ewald/n方法專門算三斜格子，現(xiàn)在和Ewald本身合并了。 我決定用三斜的格子算，因?yàn)榈V物的周期性是三斜的，用立方重復(fù)不了。但是實(shí)驗(yàn)室其他人不搞算法，他們說(shuō)關(guān)于三斜粘土礦物的文獻(xiàn)，沒(méi)有聽過(guò)把格子劃分為三斜格子的說(shuō)法。 我認(rèn)為算短程力的格子可以是立方的，但長(zhǎng)程力不行，那些文獻(xiàn)里沒(méi)有寫清楚。 不知道我的觀點(diǎn)正不正確，因?yàn)槲抑荒苷业搅⒎礁褡拥脑创a和偽代碼，沒(méi)見(jiàn)過(guò)三斜的代碼。 還有，如果有編程簡(jiǎn)單一點(diǎn)的，被大家比較公認(rèn)的算長(zhǎng)程力的算法，求推薦。我已經(jīng)要被Ewald sum坑死了。 |
回復(fù)此樓» 猜你喜歡
293求調(diào)劑
已經(jīng)有14人回復(fù)
-
0703化學(xué)336分求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有9人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請(qǐng)
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于pwscf進(jìn)行優(yōu)化的問(wèn)題
已經(jīng)有5人回復(fù)
- gromacs運(yùn)行問(wèn)題
已經(jīng)有9人回復(fù)
- LDA+U遇到的一個(gè)奇怪的問(wèn)題？？？
已經(jīng)有14人回復(fù)
- [求助成功]VASP 計(jì)算 graphene 的錯(cuò)誤
已經(jīng)有15人回復(fù)
- 分子動(dòng)力學(xué)模擬長(zhǎng)程庫(kù)侖勢(shì)的處理問(wèn)題
已經(jīng)有8人回復(fù)
- 【分享】分子模擬中作圖和看軌跡的軟件vmd windows版下載及使用手冊(cè)
已經(jīng)有13人回復(fù)
lsloneil
專家顧問(wèn) (正式寫手)
-
專家經(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)行，理論上來(lái)講跟你用什么晶系關(guān)系不大。也就是說(shuō)你依然可以把盒子里的電荷密度做一個(gè)傅立葉變換然后求和。 不過(guò)有一點(diǎn)技術(shù)上我沒(méi)太想明白的是在三斜晶系里做傅立葉變換，倒格子的格矢如何選取。立方盒子格矢很容易選取，直接2pi/L就可以了。三斜晶系下可能要做點(diǎn)額外處理來(lái)獲得一組正交的倒格子格矢？ 實(shí)空間部分求和跟你的晶系有一定關(guān)系，但這一部分很容易處理，因?yàn)橹皇莻�(gè)短程相互作用＋考慮minimum image convention。 |
2樓2015-01-11 03:35:12
3樓2015-01-11 18:31:08
lsloneil
專家顧問(wèn) (正式寫手)
-
專家經(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)助回帖
|
看不懂你在說(shuō)什么。 我又看了一下DL_POLY中Ewald sum的代碼，對(duì)于所有晶格，它們k空間部分的處理是這樣的。 1.先從晶格的格矢得到倒格子(reciprocal lattice）的格矢 2. 由倒格子的格矢再構(gòu)建一組正交的k空間矢量 3. 有了這組正交的k空間矢量，就可以獲得kmax并進(jìn)行下一步k空間的計(jì)算了。 對(duì)于步驟2，如果是立方盒子，這一步相當(dāng)于什么都沒(méi)做，因?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ā)表 | |
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123！ 2026-03-11 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|