| 9 | 1/1 | 返回列表 |
| 查看: 4742 | 回復(fù): 8 | |||
xueyiran2007銅蟲 (小有名氣)
|
[求助]
MS生成Lammps的data文件,文件中原子的類型和順序并不是自己想要的,請(qǐng)問(wèn)如何修改? 已有2人參與
|
|
我利用MS建立了水在鐵表面的模型,并導(dǎo)出data文件,但發(fā)現(xiàn)文件里H, O的原子類型分別有兩種,而且H.O的atom ID也不符合Lammps的要求, Lammps要求1個(gè)O后面跟2個(gè)H,但我的data文件中H,O的atom ID并不是按照這個(gè)規(guī)律來(lái)的,請(qǐng)問(wèn)大神,這種情況怎么處理。咳f(wàn)分感激! 請(qǐng)各位大神幫小弟一把啊 Lammps中的H,O規(guī)定如下: the atom IDs for the O and 2 H atoms must be consecutive, with the O atom first. This is to enable LAMMPS to "find" the 2 H atoms associated with each O atom. For example, if the atom ID of an O atom in a TIP4P water molecule is 500, then its 2 H atoms must have IDs 501 and 502. 我的data文件如下, 5437 atoms 3136 bonds 1571 angles 6 dihedrals 0 impropers 5 atom types 3 bond types 2 angle types 1 dihedral types 6.514843568 92.404843568 xlo xhi -0.483828184 5.242171816 ylo yhi 1.713492553 58.973492553 zlo zhi Masses 1 55.847000 2 15.999400 3 15.999400 4 1.007970 5 1.007970 ################################ 問(wèn)題1:通過(guò)Masses可以看到,這里H,O的atom type分別有2種,這是怎么回事?如何修改? Atoms # full 2290 1 2 -0.820000 68.712000000 2.863000000 54.397000000 0 0 0 # o* 2291 1 2 -0.820000 70.143500000 4.294500000 55.828500000 0 0 0 # o* 2292 1 2 -0.820000 71.575000000 2.863000000 54.397000000 0 0 0 # o* 2293 1 2 -0.820000 73.006500000 4.294500000 55.828500000 0 0 0 # o* 2294 1 2 -0.820000 74.438000000 2.863000000 54.397000000 0 0 0 # o* 2295 1 2 -0.820000 75.869500000 4.294500000 55.828500000 0 0 0 # o* 2296 1 2 -0.820000 77.301000000 2.863000000 54.397000000 0 0 0 # o* 2297 1 2 -0.820000 78.732500000 4.294500000 55.828500000 0 0 0 # o* 2298 1 2 -0.820000 80.164000000 2.863000000 54.397000000 0 0 0 # o* 2299 1 2 -0.820000 81.595500000 4.294500000 55.828500000 0 0 0 # o* 2300 1 2 -0.820000 83.027000000 2.863000000 54.397000000 0 0 0 # o* 2301 1 2 -0.820000 84.458500000 4.294500000 55.828500000 0 0 0 # o* 2302 1 4 0.410000 52.595500374 2.171500101 19.349500101 0 0 0 # h* 2303 1 4 0.410000 54.075500010 1.431500000 18.609500000 0 0 0 # h* 2304 1 4 0.410000 55.458500374 2.171500101 19.349500101 0 0 0 # h* 2305 1 4 0.410000 56.938500010 1.431500000 18.609500000 0 0 0 # h* 2306 1 4 0.410000 58.321500374 2.171500101 19.349500101 0 0 0 # h* 2307 1 4 0.410000 59.801500010 1.431500000 18.609500000 0 0 0 # h* 2308 1 4 0.410000 61.184500374 2.171500101 19.349500101 0 0 0 # h* 2309 1 4 0.410000 62.664500010 1.431500000 18.609500000 0 0 0 # h* 2310 1 4 0.410000 64.047500374 2.171500101 19.349500101 0 0 0 # h* 2311 1 4 0.410000 65.527500010 1.431500000 18.609500000 0 0 0 # h* ################################################### 問(wèn)題2:這里H,O的排列并不符合lammps的要求,請(qǐng)問(wèn)各位大神,如何修改H,O的排列順序? |

|
看見(jiàn)回/看低的帖子就忍不住留個(gè)名。 我也是新手,見(jiàn)聞?dòng)邢,求指點(diǎn)。 H,O的atom type分別有2種,這個(gè)跟具體模型有關(guān)。有的時(shí)候?yàn)榱朔謌roup或是設(shè)勢(shì)函數(shù)方便,不同位置同種原子可以設(shè)不一樣的atom type。 改排列或許可以試個(gè)土辦法,用matlab讀取文本,再編程調(diào)整,例如對(duì)每個(gè)O搜索距離最近的兩個(gè)H,再把信息挨著放。 |

木蟲 (小有名氣)

銅蟲 (小有名氣)

銅蟲 (小有名氣)

木蟲 (小有名氣)

銅蟲 (小有名氣)

|
本帖內(nèi)容被屏蔽 |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +5 | yy要上岸呀 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評(píng)估A類院;瘜W(xué)工程與工藝,誠(chéng)求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 求調(diào)劑,總分315,考的生物醫(yī)藥,一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 材料371求調(diào)劑 +9 | 鱷魚? 2026-03-11 | 11/550 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|