| 查看: 3731 | 回复: 11 | ||
[求助]
求助LAMMPS两次运算结果不一致问题！ 已有1人参与
|
|
本人打算采用LAMMPS计算高聚物的导热系数，由于是新手，首先按照LAMMPS例子里面的IN文件进行了练习（文件路径：LAMMPSExampleKAPPAin.heat），运算结果与所提供的参照结果比较吻合。而考虑到将来我的研究对象是复杂的高聚物体系，应该要用到MS来建模，为了先探探路，于是先将上述in.heat文件内描述的模型应用MS进行建模，建好模型后在Discover模块进行设置（energy里面选择compass），随后typing，并顺利导出、转换成LAMMPS所能识别的data文件（只有选择了compass才能顺利导出，后来在data文件中修改了pair_coeff）。 我碰到的问题是：从MS-->LAMMPS的方法计算结果仅为参照结果的一半（即按照LAMMPS例子中提供的in文件的计算结果）。我仔细对比过MS产生data文件与LAMMPS提供的in文件的设置，关于模型的参数（如原子数目、pair_style、pair_coeff、Mass等参数均相同），以及后续对物理过程的处理代码都是相同的，可是就是结果差别太大。（附上两次运算的LOG文件，其中Neighs数字有较大区别） 特请这方面经验丰富的朋友帮忙出出主意，或者提点宝贵意见，先感谢啦！@oxox6085@smutao |
回复此楼» 猜你喜欢
厉娜
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有15人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 硅纳米管在MS中的discover模块进行不同温度下的动力学分析 和 在lammps下的计算
已经有3人回复
- 你会发表数据不完美的结果吗？
已经有20人回复
- MS生成Lammps的data文件，文件中原子的类型和顺序并不是自己想要的，请问如何修改？
已经有8人回复
- compute 命令出错，急求
已经有5人回复
- 用lammps做分子动力学新手求助
已经有6人回复
- 经过MS模型导入到LAMMPS，在运算时出现的错误求助
已经有22人回复
- 求助LAMMPS关于MSD计算的问题
已经有17人回复
- lammps计算中出现原子距离太近现象
已经有11人回复
- lammps 离子液体，离子中相邻隔的挺远的原子成键了，不知什么原因
已经有8人回复
- LAMMPS 的坐标文件怎么构建呢？不可能一个一个原子的输入吧？
已经有8人回复
- lammps计算结果xyz文件中每个数字代表的字符长短。
已经有10人回复
- 求助数值计算结果和仿真结果有什么不一样的
已经有6人回复
- LAMMPS并行计算的问题（cpu——time关系）
已经有17人回复
- 特定分子振动的不可约表示
已经有5人回复
- 用MS建模导入LAMMPS计算，是不是不能用compass力场？
已经有9人回复
- lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
已经有7人回复
- 【讨论】lammps 用sw势计算硅的熔点怎么这么高（2400k),附in文件
已经有13人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
2楼2014-04-02 17:10:32
|
您好！我重新上传附件（其中log.MS为利用MS创建模型产生的数据文件的计算结果），请帮我分析分析，谢谢！ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Logfiles.rar
2014-04-04 09:21:30, 5.67 K
3楼2014-04-04 09:22:42
4楼2014-04-04 18:06:54
|
请看我上传到Datafiles，其中in.heat是LAMMPS例子里面的文件，用这个文件直接运行得到 log.LAMMPS；运行in.hydrogen（其中会调用hydrogen.data）产生 log.MS。请帮忙检查一下，谢谢！ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Datafiles.rar
2014-04-04 19:47:01, 20.23 K
5楼2014-04-04 19:49:52
【答案】应助回帖
★ ★ ★ ★ ★
月只蓝: 金币+5, 感谢指导！ 2014-04-07 10:27:25
月只蓝: 金币+5, 感谢指导！ 2014-04-07 10:27:25
|
1、data文件里#慎用，例如“Atoms # full”。 2、region hot block INF INF INF INF -0.4705175*v_rho 0.5294825*v_rho units box region cold block INF INF INF INF 9.5294825*v_rho 10.5294825*v_rho units box 这里设置可能得改改，否则报错： Use of region with undefined lattice 3、data文件里力场参数设置改一改： Pair Coeffs 1 1.0 1.0 这样好像可以得到类似结果。按原来的设置好像多设定了个cutoff1。 参考手册里read_data和pair_style lj/cut部分。 |
6楼2014-04-06 17:04:11
7楼2014-04-06 19:21:17
★ ★ ★ ★ ★
月只蓝: 金币+5, 感谢指导！ 2014-04-07 10:27:33
月只蓝: 金币+5, 感谢指导！ 2014-04-07 10:27:33
|
试过了。你试一下这两个文件，和你的对比一下。 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : hydrogen.data
- 附件 2 : in.hydrogen
2014-04-06 20:09:09, 353.65 K
2014-04-06 20:09:10, 1.84 K
8楼2014-04-06 20:09:12
9楼2014-04-08 08:35:51
★ ★ ★
月只蓝: 金币+3, 鼓励交流 2014-04-13 09:29:06
月只蓝: 金币+3, 鼓励交流 2014-04-13 09:29:06
|
我用的windows版本的lammps，暂时只保留了弛豫部分的run。 例子的in文件不变。自己的data文件里改一下。 你把出问题的log文件发来看看也行。 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : log.lammps
- 附件 2 : logMS.lammps
2014-04-08 09:32:28, 3.78 K
2014-04-08 09:32:36, 2.34 K
10楼2014-04-08 09:38:21
