大家好,我是小马老师。
在lammps模拟聚合物时,经常出现的一个错误是“bond missing”,也就是“键丢失”错误。
出现这种错误的直接原因是某个键两端的原子距离过大,使键被过度的拉长,或者其中一个原子直接飞出box,lammps已经找不到这个键的拓扑关系,此时就会出现“bond missing”错误提示,模拟中断。
导致键丢失的原因由很多,本文介绍其中一个原因:建模方法不合理导致键丢失。
聚合物结构比较复杂,在lammps中很难直接建模,多采用MS、packmol等软件建模。
MS的AC模块可以方便的进行聚合物的建模,但AC模块建模有一个问题,建模完成后,部分原子会处于盒子的外面,如下图所示:

转换为data文件后,盒子外面的原子根据周期性转换原理被放入到盒子中:

盒子上侧的原子被放入到盒子的下侧。
如果使用ppp边界模拟,不会出现问题。
如果改变周期性条件,如改为ppf,或者改变盒子的尺寸,就会出现键丢失。
例如,当盒子尺寸扩大5埃米之后,在ovito中也会发现原子键连接情况,很明显有些键已经被拉长:

因为盒子上下层原子已没有镜像关系,上层的原子与下层原子之间的键连接就会被取消,从而造成键丢失。
解决的办法是在ms中取消默认的盒子或者扩大盒子尺寸,使所有的原子都处于盒子内部。
转换为data文件后,所有的原子都在盒子内部,上下左右前后层之间没有键连接,此时,改变周期性或者改变盒子大小就不会出现键丢失错误。

这种方法也有个缺点,因为盒子扩大,聚合物的密度变小,需要在npt下进行驰豫已达到合理的密度。
在lammps中模拟聚合物体系特别容易出错,本文仅仅介绍了键丢失的一个原因,后续文章将会继续介绍键丢失错误原因及解决办法,敬请关注。

公众号会陆续推出势函数设置方面的文章,敬请关注。
如果在使用lammps过程中遇到什么问题,可到公众号或者QQ群留言。
公众号案例代码以及学习交流,请加QQ群:754749935

扫描关注微信公众号:lammps加油站。

lammps聚合物模拟,“bond missiong”键丢失的一个原因及解决技巧相关推荐

  1. 模拟grub损坏及丢失

    ps: 红字字体为重要部分, 仔细看 此次实验目的:  在工作中, 难免会遇到grub损坏, 现在来模拟grub损坏后修复过程. 一. 第一阶段修复MBR扇区故障; 1. 我们利用一条命令来把MBR的 ...

  2. 如何熟练掌握分子动力学LAMMPS软件模拟

    随着理论化学及计算机技术的发展,分子模拟技术被越来越多地用来加快科学研究与开发过程.分子动力学模拟方法以统计热力学.分子力学及牛顿力学为基础,对纳观体系的状态与过程进行精细研究,已被广泛用于化学化工. ...

  3. android辅助功能模拟home键,Android launcher中模拟按home键的实现

    Android launcher中模拟按home键的实现 Intent mHomeIntent = new Intent(Intent.ACTION_MAIN); mHomeIntent.addCat ...

  4. mouseup 左键_javascript中mouseup事件丢失的原因与解决办法

    这篇文章主要跟大家介绍了关于Javascript中mouseup事件丢失的原因与解决办法的相关资料,文中给出详细的示例代码供大家参考学习,需要的朋友们下面跟着小编一起来学习学习吧. 前言 当实现类似E ...

  5. android 模拟点击home,Android launcher中模拟按home键的实现

    Android launcher中模拟按home键的实现 Intent mHomeIntent = new Intent(Intent.ACTION_MAIN); mHomeIntent.addCat ...

  6. VB用API模拟截屏键PrintScreen

    很多人用 SendKeys "{PRTSC}" 模拟截屏键 PrintScreen 的时候提示<错误:'70' 拒绝的权限>,于是经常遇到人问...干脆写下来 '声明 ...

  7. php控制鼠标,HTML_用VBS控制鼠标的实现代码(获取鼠标坐标、鼠标移动、鼠标单击、鼠标双击、鼠标右击),怎么用VBS模拟鼠标左键单击、 - phpStudy...

    用VBS控制鼠标的实现代码(获取鼠标坐标.鼠标移动.鼠标单击.鼠标双击.鼠标右击) 怎么用VBS模拟鼠标左键单击.左键双击.右键单击?-- 网上搜到的答案普遍是VBS无法实现,或者是要用第三方COM( ...

  8. php左键点改成左键点,window_win10怎么将鼠标左键设置为主键?win10鼠标左右键互换的技巧,当我们的鼠标使用一段时间之 - phpStudy...

    win10怎么将鼠标左键设置为主键?win10鼠标左右键互换的技巧 当我们的鼠标使用一段时间之后,会由于左键经常按压而导致老化或者损坏,而对于突然间的损坏,又没有其他可以替代的鼠标时候,就需要对鼠标左 ...

  9. Hyper-V下的Linux虚拟机网卡丢失问题原因及解决办法

    Hyper-V下的Linux虚拟机网卡丢失问题原因及解决办法 虚拟化大势所趋 公司推行了虚拟化,全部用的是Microsoft Windows 2008 R2 Enterprise with Hyper ...

  10. 计算机中丢失setupxml.dll,Win7电脑安装VideoStudio Pro X6显示丢失SetupXML.dll文件怎么解决...

    最近有win7系统用户在电脑安装VideoStudio Pro X6软件的时候,突然出现错误的提示,显示无法启动此程序,因为计算机中丢失SetupXML.dll.尝试重新安装该程序来解决此问题,要怎么 ...

最新文章

  1. Linux下的USB总线驱动 3
  2. 简单python画圣诞树图片-python圣诞树编写实例详解
  3. OpenCASCADE绘制测试线束:OCAF 命令之树属性命令
  4. 【Windows CMD】设置编码格式
  5. oracle电梯案例,Oracle技术嘉年华的一个案例,redo的那些事,连载一
  6. python随机生成验证码_Python生成随机验证码
  7. Netty内存池泄漏问题
  8. Atitit 互联网2018风口与趋势大盘点 目录 1.1. 短视频 1 1.2. 互联网金融 1 1.3. 区块链 1 1.4. 新零售 1 2. 2018热度榜丨你能抓住的风口有几个?.html
  9. (丘维声)高等代数课程笔记:商空间
  10. (4.2.48)MVPArms源码分析
  11. centos7离线安装wget
  12. 1972 年 11 月 29 日:雅达利推出投币式街机游戏《乓》
  13. colormap保存 matlab_matlab中colormap用法
  14. 高德离线地图开发API
  15. 关于oracle预言机
  16. 时间序列基础操作:使用python与eviews对AR与ARMA模型进行定阶与预报
  17. 删除Karabiner-Elements 和 Karabiner-EventViewer提示项目被锁定
  18. 跑步感到痛苦,试试超慢跑
  19. route----添加路由的方法
  20. lammps 案例in文件汇总

热门文章

  1. 俯瞰开源工作流引擎Activiti
  2. nvme linux读写测试工具,不同硬件平台(intel/AMD)和不同OS/FS/测试工具下NVMe SSD性能简测...
  3. 网吧用电影服务器系统,网吧流媒体电影服务器搭建的解决方案
  4. breadweb控制台下载_路由器刷breed web控制台助手通用版下载
  5. Win10系统优化工具
  6. 一个小时,零基础入门,看完这篇30行代码 教你实现百度换肤!
  7. 八大常用渗透测试工具
  8. coreos 安装 mysql_手把手教你在 CoreOS 上构建你的第一个应用
  9. WPS加载项深入开发代码ribbon解析
  10. Windows 无法访问指定设备、路径或文件