| 3 | 1/1 | 返回列表 |
| 查看: 1749 | 回復(fù): 2 | ||
胡博洋的洋金蟲 (小有名氣)
|
[求助]
求助各位前輩關(guān)于在LAMMPS中angle_style的設(shè)置問題。
|
|
大家好,我需要進(jìn)行LiFePO4的分子動(dòng)力學(xué)模擬,因?yàn)槭堑谝淮问褂肔AMMPS,很多地方不懂,請(qǐng)大家理解。 我在angle_style命令可以選擇harmonic,對(duì)應(yīng)的angle_coeff通過查文獻(xiàn)也可以找到。但是我有一個(gè)疑問,除了這兩個(gè)命令以外,沒有感覺到有定義說這個(gè)角是由哪三個(gè)原子組成,比如我要定義O-P-O的三體彎曲勢(shì)(Here it is defined for each O-P-O bond and takes the form of a harmonic angle-bending potential about the central P ion) 我的命令是 angle_style harmonic angle_coeff 1 0.6613 109.47 結(jié)果錯(cuò)誤提示說我定義的原子類型(我定義的是charge)不可以定義角。(@#¥……&*真的是原子類型的問題嗎),然后我又看manual,發(fā)現(xiàn)在data文件里有這樣說: Angles section: • one line per angle • line syntax: ID type atom1 atom2 atom3 ID = number of angle (1-Nangles) type = angle type (1-Nangletype) atom1,atom2,atom3 = IDs of 1st,2nd,3rd atoms in angle 于是我覺得有希望了,但是有一點(diǎn)還是不明白,就是這三個(gè)原子的ID,它的意思是指原子的ID吧,也就是說要一步一步的輸入原子的ID,也就是需要自己判斷是哪三個(gè)原子所成的角。自己判斷的話我起初是打算把與中間的那個(gè)P相距最近的四個(gè)O算出來,后來發(fā)現(xiàn)只有三個(gè)距離特別近且?guī)缀跻粯,第四個(gè)最近的O至少也是剛才那三個(gè)距離的5倍。因?yàn)檎G闆r應(yīng)該是P位于四個(gè)O組成的四面體的中心,所以不應(yīng)該出現(xiàn)上述情況, 1、想請(qǐng)問下大家出現(xiàn)類似情況是什么原因呢?原子坐標(biāo)啥的計(jì)算我可以確定是沒問題的。 2、在考慮是不是周期性邊界條件啥的,真正的那個(gè)O在另一個(gè)區(qū)域? 3、推到重來,應(yīng)該怎么樣定義三體彎曲勢(shì)? 4、原子類型為charge真的不可以定義鍵角嗎? 5、請(qǐng)求各位前輩指點(diǎn)~~ |

木蟲 (正式寫手)
金蟲 (小有名氣)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +6 | 上岸6666@ 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +3 | jiajunX 2026-03-22 | 3/150 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 17/850 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗(yàn)) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|