 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1921 | 回復(fù): 5 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
[求助]
關(guān)于Ewald sum算法的相關(guān)問(wèn)題 已有1人參與
|
|||||
|
自己作死地用C編寫(xiě)Monte Carlo的代碼，模擬的具體問(wèn)題是粘土礦物對(duì)Ar的吸附。 然后需要加上Ewald Sum算法的問(wèn)題算礦物羥基的長(zhǎng)程力。 但是這個(gè)礦物是三斜晶系的，書(shū)上給的Ewald Sum是立方格子。 我查過(guò)LAMMPS里面，前幾年有個(gè)Ewald/n方法專(zhuān)門(mén)算三斜格子，現(xiàn)在和Ewald本身合并了。 我決定用三斜的格子算，因?yàn)榈V物的周期性是三斜的，用立方重復(fù)不了。但是實(shí)驗(yàn)室其他人不搞算法，他們說(shuō)關(guān)于三斜粘土礦物的文獻(xiàn)，沒(méi)有聽(tīng)過(guò)把格子劃分為三斜格子的說(shuō)法。 我認(rèn)為算短程力的格子可以是立方的，但長(zhǎng)程力不行，那些文獻(xiàn)里沒(méi)有寫(xiě)清楚。 不知道我的觀點(diǎn)正不正確，因?yàn)槲抑荒苷业搅⒎礁褡拥脑创a和偽代碼，沒(méi)見(jiàn)過(guò)三斜的代碼。 還有，如果有編程簡(jiǎn)單一點(diǎn)的，被大家比較公認(rèn)的算長(zhǎng)程力的算法，求推薦。我已經(jīng)要被Ewald sum坑死了。 |
回復(fù)此樓» 猜你喜歡
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專(zhuān)業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有8人回復(fù)
-
341求調(diào)劑(一志愿湖南大學(xué)070300)
已經(jīng)有6人回復(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
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲(chóng)號(hào): 556471
- 注冊(cè): 2008-05-11
- 性別: GG
- 專(zhuān)業(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
lsloneil
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲(chóng)號(hào): 556471
- 注冊(cè): 2008-05-11
- 性別: GG
- 專(zhuān)業(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
王慧1021
金蟲(chóng) (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1018.5
- 帖子: 140
- 在線: 144.6小時(shí)
- 蟲(chóng)號(hào): 1535805
- 注冊(cè): 2011-12-13
- 性別: MM
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
5樓2015-01-13 11:48:52
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +10 | 梁梁校校 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料專(zhuān)碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
|
[基金申請(qǐng)] 請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎？ +3 | lishucheng1 2026-03-22 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚(yú).. 2026-03-21 | 4/200 |
|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|