 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1919 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
關(guān)于Ewald sum算法的相關(guān)問題 已有1人參與
|
||||
|
自己作死地用C編寫Monte Carlo的代碼，模擬的具體問題是粘土礦物對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ù)此樓» 猜你喜歡
材料專業(yè)求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有8人回復(fù)
-
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑
已經(jīng)有6人回復(fù)
-
請教下大家 2026年國家基金申請是雙盲審嗎？
已經(jīng)有4人回復(fù)
-
一志愿南航材料專317分求調(diào)劑
已經(jīng)有3人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271
已經(jīng)有4人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有8人回復(fù)
-
化學(xué)308分求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有11人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 關(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é)模擬長程庫侖勢的處理問題
已經(jīng)有8人回復(fù)
- 【分享】分子模擬中作圖和看軌跡的軟件vmd windows版下載及使用手冊
已經(jīng)有13人回復(fù)
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í)
- 蟲號: 556471
- 注冊: 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
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲號: 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
|
看不懂你在說什么。 我又看了一下DL_POLY中Ewald sum的代碼，對于所有晶格，它們k空間部分的處理是這樣的。 1.先從晶格的格矢得到倒格子(reciprocal lattice）的格矢 2. 由倒格子的格矢再構(gòu)建一組正交的k空間矢量 3. 有了這組正交的k空間矢量，就可以獲得kmax并進(jìn)行下一步k空間的計(jì)算了。 對于步驟2，如果是立方盒子，這一步相當(dāng)于什么都沒做，因?yàn)樵鹊垢褡拥母袷妇褪钦坏摹５珜τ谌钡母褡�，因�(yàn)槠涞垢褡痈袷覆皇钦坏模@一步就可以獲得一組正交的k空間矢量。步驟2的算法的大概意思是，對于三斜晶系，總可以找到一組正交的格矢，它們形成的晶胞體積跟原晶胞相等。 如果你想看dlpoly的代碼，可以在網(wǎng)上搜DL_POLY classic的代碼（下載需要注冊，但免費(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í)
- 蟲號: 1535805
- 注冊: 2011-12-13
- 性別: MM
- 專業(yè): 理論和計(jì)算化學(xué)
5樓2015-01-13 11:48:52
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 07化學(xué)280分求調(diào)劑 +3 | 722865 2026-03-23 | 3/150 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+4
|
孅華 2026-03-22 | 4/200 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +11 | 梨花珞晚風(fēng) 2026-03-17 | 11/550 |
|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|