24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4796  |  回复: 8

xueyiran2007

铜虫 (小有名气)

[求助] MS生成Lammps的data文件,文件中原子的类型和顺序并不是自己想要的,请问如何修改? 已有2人参与

我利用MS建立了水在铁表面的模型,并导出data文件,但发现文件里H, O的原子类型分别有两种,而且H.O的atom ID也不符合Lammps的要求,
Lammps要求1个O后面跟2个H,但我的data文件中H,O的atom ID并不是按照这个规律来的,请问大神,这种情况怎么处理啊?万分感激!
请各位大神帮小弟一把啊

Lammps中的H,O规定如下:
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的要求,请问各位大神,如何修改H,O的排列顺序?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

LiveandLearn
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ysshj

金虫 (小有名气)


月只蓝: 金币+1, 鼓励交流 2014-01-10 19:35:54
看见回/看低的帖子就忍不住留个名。
我也是新手,见闻有限,求指点。

H,O的atom type分别有2种,这个跟具体模型有关。有的时候为了分group或是设势函数方便,不同位置同种原子可以设不一样的atom type。

改排列或许可以试个土办法,用matlab读取文本,再编程调整,例如对每个O搜索距离最近的两个H,再把信息挨着放。
一个人要像一支队伍。
2楼2014-01-10 17:00:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茹此无声

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
月只蓝: 金币+1, 鼓励交流 2014-01-11 07:40:50
xueyiran2007: 金币+5, ★★★很有帮助 2014-01-12 03:43:29
1.ms文件导入到ME里,bdl转成data.
2.fortran编程实现。

[ 发自小木虫客户端 ]
人类一思考,上帝就发笑。
3楼2014-01-11 00:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueyiran2007

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by ysshj at 2014-01-10 17:00:14
看见回/看低的帖子就忍不住留个名。
我也是新手,见闻有限,求指点。

H,O的atom type分别有2种,这个跟具体模型有关。有的时候为了分group或是设势函数方便,不同位置同种原子可以设不一样的atom type。

改排 ...

多谢!我试试
LiveandLearn
4楼2014-01-11 03:40:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueyiran2007

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by 茹此无声 at 2014-01-11 00:23:25
1.ms文件导入到ME里,bdl转成data.
2.fortran编程实现。

多谢!大神,您可以说的详细点吗?bdl是什么?我已经用msi2lmp转换成data文件了,您的意思是用bdl重新转换吗?另外如何利用fortran编程实现?是通过编程对data文件进行排序吗?可以给个例子吗?
LiveandLearn
5楼2014-01-12 08:36:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茹此无声

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
xueyiran2007: 金币+5, ★★★很有帮助 2014-01-13 04:23:31
月只蓝: 金币+2, 感谢指导! 2014-03-12 19:21:39
引用回帖:
5楼: Originally posted by xueyiran2007 at 2014-01-12 08:36:55
多谢!大神,您可以说的详细点吗?bdl是什么?我已经用msi2lmp转换成data文件了,您的意思是用bdl重新转换吗?另外如何利用fortran编程实现?是通过编程对data文件进行排序吗?可以给个例子吗?...

ms里导出来的文件是pdb格式,用me里的转化工具可以变为bdl格式。bdl格式再转化为data格式就好了。了解每个格式的特定的规则编程。我也新手,只知道这么多了。

[ 发自小木虫客户端 ]
人类一思考,上帝就发笑。
6楼2014-01-12 11:20:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueyiran2007

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by 茹此无声 at 2014-01-12 11:20:36
ms里导出来的文件是pdb格式,用me里的转化工具可以变为bdl格式。bdl格式再转化为data格式就好了。了解每个格式的特定的规则编程。我也新手,只知道这么多了。
...

多谢!
LiveandLearn
7楼2014-01-13 04:23:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huachao90

禁虫 (正式写手)


月只蓝: 金币+1, 鼓励交流! 2014-03-12 19:21:44
本帖内容被屏蔽

8楼2014-03-12 18:22:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alwynwen

新虫 (初入文坛)

请问楼主最后是怎么解决这个问题的?谢谢
9楼2014-09-16 19:30:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xueyiran2007 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂 +16 艾志恒 2026-04-03 17/850 2026-04-06 00:27 by 永字号
[考研] 086000生物与医药 初试274求调剂 +6 小叮当来了 2026-03-30 7/350 2026-04-05 20:30 by lys0704
[考研] 材料调剂 +13 一样YWY 2026-04-03 14/700 2026-04-05 18:20 by 蓝云思雨
[考研] 求调剂 +4 wos666 2026-04-03 4/200 2026-04-05 11:48 by arrow8852
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 一志愿北京2,材料与化工308求调剂 +10 熊二想上岸 2026-04-04 10/500 2026-04-05 05:20 by houyaoxu
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 考研调剂 +6 15615482637 2026-04-04 6/300 2026-04-04 22:43 by yu221
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 331求调剂 +3 niby 2026-04-02 3/150 2026-04-04 19:56 by 蓝云思雨
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 285求调剂 +7 AZMK 2026-04-02 9/450 2026-04-03 11:12 by wanwan00
[考研] 311求调剂一志愿合肥工业大学 +15 秋二十二 2026-03-30 15/750 2026-04-03 10:19 by linyelide
[考研] 085801 总分275 本科新能源 求调剂 +19 bradoner 2026-04-01 23/1150 2026-04-03 10:07 by linyelide
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 13:03 by yulian1987
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 一志愿北交大材料工程,总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 一志愿中海洋材料357 +4 麦恩莉. 2026-03-30 4/200 2026-03-31 14:35 by 记事本2026
信息提示
请填处理意见