计算化学公社

标题: 求助：lammps模拟固态-固态相变时在弛豫或升温过程中出现结构不合理的地方 [打印本页]

作者
Author: XWH 时间: 2021-6-3 11:11
标题: 求助：lammps模拟固态-固态相变时在弛豫或升温过程中出现结构不合理的地方
本帖最后由 XWH 于 2021-6-3 11:11 编辑

用lammps模拟的体系是固态相变，金属氧化物于70摄氏度发生固态相变，由单斜变为四方，体系目前仅考虑升温过程，目前找到的势函数有三种，1)库仑力+莫尔斯势（玻璃体系中找到的）；2）白金汉势+莫尔斯势（gulp中找的代入lammps的）；3)反应力场（含有相同元素的金属氧化物，其他软件模拟单斜到四方的相变）；250K-300K温度下都试着试着仅用nvt,npt进行弛豫，后面一用npt升温，试过iso,aniso,tri控压，但是前两种势函数无论用nvt，还是npt，只要一弛豫，连续3个蓝色的原子间的两个键长趋于相等，本该扭转的弛豫不到几步（100万步）就几乎变直了，而只有反应力场跑nvt弛豫时才会保证初始结构基本不变，但只要一跑npt，但弛豫后一升温，没跑几步就跟之前两个势函数一样的情况。有办法可以控制这种波动吗？（该体系参考相变参考资料只有第一性原理的）一直测试的是758或者1536或者2274左右的原子，看书说系综下原子数目多会提高精确度，要是在大幅度增大原子数目的情况下，可以控制这种情况吗？还是势函数不合适？（找到相关固态-固态相变也有用过混合势模拟成功的）3种模拟命令写法、结构变化如下

作者
Author: 李江山 时间: 2023-3-21 23:03
你好，我遇到同样的问题，在NPT系综下升温观察不到相变，由于测试的时候用得是500个原子左右，是不是扩大体系可以看到明显相变呢？

作者
Author: sychang1 时间: 2023-6-7 22:12

李江山发表于 2023-3-21 23:03
你好，我遇到同样的问题，在NPT系综下升温观察不到相变，由于测试的时候用得是500个原子左右，是不是扩大体 ...

扩大体系试试，小体系参数变化不明显

欢迎光临计算化学公社 (http://bbs.keinsci.com/)