lammps混合势函数设置,常见的一个错误解决方案
大家好,我是小马老师!
最近在咨询过程中,发现不少朋友对势函数混合势写法不是很清楚,以致运行出错。
比较常见的一个错误是:All pair coeffs are not set,如下图所示:
错误提示为“不是所有的原子都设置了势函数”,表示有的原子可能没有设置势参数。
出现这种错误主要有两种原因:
第一,因为原子数量多,确实少写了某种原子的势参数,这种情况仔细检查,找出未设置势的原子,补充设置即可。
第二,设置了全部原子的势函数,但是设置方法出错,这种情况再怎么检查原子势也不能解决问题。
因为表面上看,所有的势已经全部被设置。
下面以一个简单的Cu-Al界面模型为例,给出这种错误的解决方法。
界面模型如下图所示,下层为Cu,原子类别1,上层为Al,原子类型2,Cu和Al均使用meam势,cu-Al之间的势使用lj势。
势函数设置如下:
pair_style hybrid meam/c lj/cut 10
pair_coeff * * meam/c library.Cu.meam Cu Cu.meam Cu NULL #Cu
pair_coeff * * meam/c library-Al.meam Al Al.meam NULL Al #Al
pair_coeff 1 2 lj/cut 0.4 2.47 #Cu-Al
如果按照以上代码设置,就会提示“All pair coeffs are not set”,可代码中已经设置了Cu、Al以及Cu-Al的所有势。
出现这种错误,主要是因为在同一in文件使用了两个meam势,第二个Al的meam势会覆盖掉Cu的meam势,导致Cu的meam势参数丢失,所以程序给出势参数设置不全的提示。
解决的办法:
如果在同一个in文件中使用多个同种类型的势,为防止覆盖,需要对相同类型的势进行编号区分,所以正确的写法为:
#类型设置中把所有的meam/c全部列出,本文需写2个meam/c
pair_style hybrid meam/c meam/c lj/cut 10
#在meam/c后面增加编号1 2
pair_coeff * * meam/c 1 library.Cu.meam Cu Cu.meam Cu NULL
pair_coeff * * meam/c 2 library-Al.meam Al Al.meam NULL Al
pair_coeff 1 2 lj/cut 0.4 2.47 #Cu-Al
这时候再进行模拟,就没有错误提示了。
下面给出全部in文件代码以供参考:
units metal
atom_style atomicregion box block 0 80 0 10 0 40 units box
create_box 2 boxregion cu block INF INF INF INF INF 20 units box
lattice fcc 3.61
create_atoms 1 region curegion al block INF INF INF INF 21 INF units box
lattice fcc 4.05
create_atoms 2 region almass 1 64
mass 2 27pair_style hybrid meam/c meam/c lj/cut 10
pair_coeff * * meam/c 1 library.Cu.meam Cu Cu.meam Cu NULL
pair_coeff * * meam/c 2 library-Al.meam Al Al.meam NULL Al
pair_coeff 1 2 lj/cut 0.4 2.47 #Cu-Alwrite_data CuAl.data
复合势函数的设置相对比较复杂,也是lammps学习过程中的一个难点,公众号后陆续推出势函数设置技巧,敬请关注!
如果在使用lammps过程中遇到什么问题,可到公众号或者QQ群留言。
公众号案例代码以及学习交流,请加QQ群:754749935
扫描关注微信公众号:lammps加油站。
lammps混合势函数设置,常见的一个错误解决方案相关推荐
- 盘点一个网络爬虫中常见的一个错误
点击上方"Python爬虫与数据挖掘",进行关注 回复"书籍"即可获赠Python从入门到进阶共10本电子书 今 日 鸡 汤 圣朝无阙事,自觉谏书稀. 大家好, ...
- Dnn中常见的一个错误及可能原因:Could not find stored procedure 'dbo.GetPortals'
Could not find stored procedure 'dbo.GetPortals' 错误的 {objectqualifier}. 数据库的版本保存在dnn.config 中,assemb ...
- mujoco强化学习环境配置,常见OSError: __glewBindBuffer错误解决方案
文章目录 一.mujoco强化学习安装 二.遇到错误 __glewBindBuffer 一.mujoco强化学习安装 参考链接1 ubantu18参考链接 注意的是,现在mujoco开源了,直接官网下 ...
- 在lammps模拟过程中的常用势函数设置
文章目录 1 lj/cut 1.1 lj/cut在in文件中使用方法 1.2 lj/cut在data文件中使用方法 1.3 lj/cut参数查询方法 1.4 lj/cut参数单位转换方法 1.5 lj ...
- lammps模拟技巧:高熵合金势函数设置三种方法
大家好,我是小马老师. 本文继续介绍lammps模拟中的势函数设置方法. 在前面的文章中,介绍过高熵合金的建模:高熵合金建模 高熵合金包含的原子数较多,势函数的设置相对复杂,在平时一对一的咨询中,这种 ...
- co01设置错误消息_win10提示“调制解调器报告了一个错误”怎么办
在使用windows10系统过程中,一位用户反馈电脑无法连接宽带连接,而且出现提示:无法连接到宽带链接 调制解调器(或其他连接设备)报告了一个错误. 这是怎么回事呢?遇到这样的情况是由于调制解调器设置 ...
- 常见的 Git 错误解决方法
常见的 Git 错误解决方法 从索引中取消暂存文件/目录 编辑上次提交的邮件 忘记了上次提交中的一些更改 放弃本地更改 已将个人数据提交到本地存储库 password files 将最新的提交替换为新 ...
- 6种常见的Git错误以及解决的办法
我们都会犯错误,尤其是在使用像Git这样复杂的东西时.如果你是Git的新手,可以学习如何在命令行上开始使用Git.下面介绍如何解决六个最常见的Git错误. Photo by Pawel Janiak ...
- 最常见的HTTP错误
1. HTTP 500错误(内部服务器错误)对对HTTP 500错误的定义已经充分证明了这是一个最常见的HTTP错误. 一般来说,HTTP 500 错误就是web服务器发生内部错误时返回的信息. 例如 ...
最新文章
- 计算机技能需求新排名:Python 仅排第 3,第 1 你可能猜不到哦
- RAID0、RAID1、RAID0+1模式实战评测
- 联络中心的发展方向是SOA
- python如何让用户输入文件名并打开文件_(Python)如何让用户打开文本文件然后更改整数/数字...
- 第九期: 阿里巴巴程序员常用的15款开发者工具
- nrf51822添加UUID
- ReactNative学习笔记(二)Flex布局
- HDU-3460 Ancient Printer 字典树
- hooks组件封装 react_名符其实的react下一代状态管理器hox
- 唐山师范学院计算机宿舍,唐山师范学院宿舍条件怎么样宿舍图片内景
- 鹅厂优文 | ReactJS一点通
- Java 1.2.3 文件输入与输出
- 学python能做什么-学了Python都能干什么,哪个最赚钱?
- java语言有很多的优点和特点_Java语言具有许多优点和特点,下列选项中()反映了Java程序并行机制的特点 (5.0分)_学小易找答案...
- java 分词获取词性_jie分词之分词、词性分析、关键词抽取
- 2020上半年十大典型勒索软件大盘点
- Top计算机博士去BATH还是985非升即走?
- 智能窗帘控制(语音、时间、光、红外)
- 学习日志day45(2021-09-09)(1、有道翻译API使用json格式数据 2、JSONP 3、Ajax提交form表单 4、Web Uploader)
- Oracle DBlink创建访问远程数据库
热门文章
- 【Leetcode】748. Shortest Completing Word
- RTL8188CUS WIFI模块的使用方式
- cad菜单栏快捷键_天正CAD界面“菜单栏”不见了怎么办?教你3招秒解决,实用神技巧...
- Metasploit入门手册(msfconsole)
- Metasploit基本用法
- 那一年我学习软件测试,现在没有后悔
- 网易云信 QUIC 应用优化实践
- Pytorch 安装
- 官方发布历年考研平均分!你今年考多少分才能上岸
- 【MATLAB appdesigner】22_界面美化_如何添加“图标icon“?