lammps模拟“lost atoms”原子原因及解决方法

2023-11-04 01:31

本文主要是介绍lammps模拟“lost atoms”原子原因及解决方法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

lammps模拟过程中,最怕的是模拟过程中出错,其中,比较常见的一个错误是“lost atoms”,也就是常说的原子丢失。

正常情况下,分子动力学模拟要保证原子数目保持不变。

nve、nvt、npt系综中的“n”就是原子数目,这些系综已经明确了模拟过程中要保证原子数目“n”不变,原子数目减少则会中止模拟,并给出“lost atoms”错误提示。

“lost atoms”错误出现的最主要的原因是某些原子在积分过程中求解得到的速度过大,根据该速度计算出的下一步位移过大,已经超出了体系盒子的范围,换句话说,速度过快的原子已经跑出了盒子,造成“lost atoms”。

导致原子速度过快的原因有以下几个方面:

(1)模型设置不合理

一般出现在多原子建模中,比如一层石墨烯原子插入到Al原子之间,由于建模时石墨烯C原子位置设置不合理,导致与Al原子重叠。
在这里插入图片描述

在进行积分时,位置重叠或者距离过近的原子之间的排斥力过大,因此原子获得的速度过大,容易飞出box。

解决方法:
建模时注意原子排列位置,宁可离的稍远一些,也不要靠的过近。

建模之后使用minimize命令对体系进行能量最小化,minimize命令不用设置力场,是拓扑结构上对体系进行最小化,能够解决大部分的原子重叠问题。

(2)边界条件设置不合理

如果建模过程没有问题,也没出现明显的原子重叠现象,考虑一下是否是边界条件原因。

例如,边界条件为固定f条件,由于建模原因,模型的初始内应力又比较大,正常弛豫下,盒子应该会扩大以释放掉内应力从而达到系统平衡,而固定边界条件又不允许盒子尺寸发生变化,此时,边界上的原子会被挤出盒子,造成原子丢失。
在这里插入图片描述

解决方法:

将边界条件改为s,查看驰豫结果,如果盒子只是稍微扩大,没有出现盒子爆炸性的扩大,基本上属于这种问题,将初始盒子尺寸扩大可以解决此类问题。
(3)力场设置不合理

如果是发生了盒子“爆炸”,原子完全散开,单纯的扩大盒子尺寸没有意义,需要考虑力场设置是否合理。

力场描述原子之间的受力情况,如果力场设置不合理,可能会出现原子之间受力计算错误,原子之间的力远超正常值,也会造成原子丢失。

解决方法:更换力场。

(4)模拟步长设置不合理

当以上原因都排除后,如果还是存在“lost atoms”,考虑是否是timestep设置过大,设置更小的步长试试。

当出现丢失原子错误提示时,添加thermo_modify lost ignore 命令可以关掉错误提示,使程序能够运行下去。

但是,不建议这么做,这个命令仅仅是强制程序继续运行,不能解决原子丢失问题,程序“带病”运行,有可能会得到不正确的运行结果。

更多lammps教程请关注微信公众号:lammps加油站
扫码关注微信公众号
lammps交流请加qq群:
在这里插入图片描述

这篇关于lammps模拟“lost atoms”原子原因及解决方法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/342294

相关文章

使用Python进行文件读写操作的基本方法

《使用Python进行文件读写操作的基本方法》今天的内容来介绍Python中进行文件读写操作的方法,这在学习Python时是必不可少的技术点,希望可以帮助到正在学习python的小伙伴,以下是Pyth... 目录一、文件读取:二、文件写入:三、文件追加:四、文件读写的二进制模式:五、使用 json 模块读写

解决systemctl reload nginx重启Nginx服务报错:Job for nginx.service invalid问题

《解决systemctlreloadnginx重启Nginx服务报错:Jobfornginx.serviceinvalid问题》文章描述了通过`systemctlstatusnginx.se... 目录systemctl reload nginx重启Nginx服务报错:Job for nginx.javas

Oracle数据库使用 listagg去重删除重复数据的方法汇总

《Oracle数据库使用listagg去重删除重复数据的方法汇总》文章介绍了在Oracle数据库中使用LISTAGG和XMLAGG函数进行字符串聚合并去重的方法,包括去重聚合、使用XML解析和CLO... 目录案例表第一种:使用wm_concat() + distinct去重聚合第二种:使用listagg,

Java后端接口中提取请求头中的Cookie和Token的方法

《Java后端接口中提取请求头中的Cookie和Token的方法》在现代Web开发中,HTTP请求头(Header)是客户端与服务器之间传递信息的重要方式之一,本文将详细介绍如何在Java后端(以Sp... 目录引言1. 背景1.1 什么是 HTTP 请求头?1.2 为什么需要提取请求头?2. 使用 Spr

Java如何通过反射机制获取数据类对象的属性及方法

《Java如何通过反射机制获取数据类对象的属性及方法》文章介绍了如何使用Java反射机制获取类对象的所有属性及其对应的get、set方法,以及如何通过反射机制实现类对象的实例化,感兴趣的朋友跟随小编一... 目录一、通过反射机制获取类对象的所有属性以及相应的get、set方法1.遍历类对象的所有属性2.获取

Java中的Opencv简介与开发环境部署方法

《Java中的Opencv简介与开发环境部署方法》OpenCV是一个开源的计算机视觉和图像处理库,提供了丰富的图像处理算法和工具,它支持多种图像处理和计算机视觉算法,可以用于物体识别与跟踪、图像分割与... 目录1.Opencv简介Opencv的应用2.Java使用OpenCV进行图像操作opencv安装j

Ubuntu系统怎么安装Warp? 新一代AI 终端神器安装使用方法

《Ubuntu系统怎么安装Warp?新一代AI终端神器安装使用方法》Warp是一款使用Rust开发的现代化AI终端工具,该怎么再Ubuntu系统中安装使用呢?下面我们就来看看详细教程... Warp Terminal 是一款使用 Rust 开发的现代化「AI 终端」工具。最初它只支持 MACOS,但在 20

Python实现数据清洗的18种方法

《Python实现数据清洗的18种方法》本文主要介绍了Python实现数据清洗的18种方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学... 目录1. 去除字符串两边空格2. 转换数据类型3. 大小写转换4. 移除列表中的重复元素5. 快速统

Mysql DATETIME 毫秒坑的解决

《MysqlDATETIME毫秒坑的解决》本文主要介绍了MysqlDATETIME毫秒坑的解决,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着... 今天写代码突发一个诡异的 bug,代码逻辑大概如下。1. 新增退款单记录boolean save = s

Debian如何查看系统版本? 7种轻松查看Debian版本信息的实用方法

《Debian如何查看系统版本?7种轻松查看Debian版本信息的实用方法》Debian是一个广泛使用的Linux发行版,用户有时需要查看其版本信息以进行系统管理、故障排除或兼容性检查,在Debia... 作为最受欢迎的 linux 发行版之一,Debian 的版本信息在日常使用和系统维护中起着至关重要的作