lammps问题-关于100*dt的使用

托马斯922

已于 2022-06-26 09:32:26 修改

阅读量770

点赞数

文章标签： java 开发语言

于 2022-06-25 16:23:25 首次发布

本文链接：https://blog.csdn.net/weixin_43494235/article/details/125460595

版权

在LAMMPS的NVT ensemble模拟中，使用`$(100.0*dt)`作为`Tdamp`参数时，系统在弛豫后出现回缩现象。对比实验发现，将`Tdamp`改为固定值10.0后，问题得到改善。文章探讨了`Tdamp`设置对模拟结果的影响，并分享了完整的模拟配置和步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

在lammps官网介绍fix nvt使用时有以下介绍在这里插入图片描述关于Tdamp的使用

在偶尔试一下的时候发现了一下这个问题
1.100dt
fix 1 all nvt temp 300 300 $(100.0dt)
在这里插入图片描述弛豫后出现了往中间回缩的情况。

在不改变其他条件的情况下，变为
fix 1 all nvt temp 300 300 10.0
在这里插入图片描述
源代码如下：
shell mkdir initial-comp
shell mkdir config-comp
shell mkdir restart-comp
#--------------------------DEFINITION-------------------
units metal
dimension 3
boundary s s p
atom_style atomic
read_data Ti_bulk.lmp
region 1 block INF INF INF INF INF 5 units box
region 2 block INF INF INF INF 221 INF units b