| 5 | 1/1 | 返回列表 |
| 查看: 4789 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
xueyiran2007銅蟲 (小有名氣)
|
[求助]
MS生成Lammps的data文件,文件中原子的類型和順序并不是自己想要的,請(qǐng)問如何修改? 已有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ī)律來的,請(qǐng)問大神,這種情況怎么處理。咳f分感激! 請(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 ################################ 問題1:通過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* ################################################### 問題2:這里H,O的排列并不符合lammps的要求,請(qǐng)問各位大神,如何修改H,O的排列順序? |

木蟲 (小有名氣)

|
看見回/看低的帖子就忍不住留個(gè)名。 我也是新手,見聞?dòng)邢蓿笾更c(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,再把信息挨著放。 |

木蟲 (小有名氣)

銅蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 324分 085600材料與化工 +7 | 呆鵝oor 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿北林食品與營養(yǎng)095500,301分,已過六級(jí),有科研經(jīng)歷 +3 | 快樂儲(chǔ)蓄罐 2026-03-31 | 3/150 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +11 | Nadiums 2026-03-25 | 12/600 |
|
|
[考研] 308求調(diào)劑 +11 | 墨墨漠 2026-03-25 | 11/550 |
|
|
[考研] 一志愿中國科學(xué)院大學(xué)265求調(diào)劑 +6 | 恬淡ye 2026-03-31 | 7/350 |
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
[考研] 英一數(shù)一408,總分284,二戰(zhàn)真誠求調(diào)劑 +3 | 12.27 2026-03-30 | 5/250 |
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|