| 查看: 1384 | 回复: 0 | |||
[交流]
求助：lammps模拟固态-固态相变在弛豫和升温过程结构结构变化不理想问题
|
用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种模拟命令写法、结构变化如下      发自小木虫Android客户端 |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有282人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
