对于lammps新手来说,系综的选择令人困惑,到底是应该用nvt系综还是npt系综?

我在初学lammps时就被这个问题困扰了很久。如果去问别人,可能会告诉你一个答案:根据模拟需要确定。

其实和没说一样。

到底如何确定系综,确实不是一句两句能够说清楚的,不同的体系、不同的模拟目的要选择不同的系综。

前几天遇到的一个合金的模拟案例,比较有代表性,分享出来供大家参考。

建模方法和之前介绍过的高熵合金建模方法类似,先建立单晶模型,然后使用替换原子的命令获得合金结构。

建模完成后,按照常规,先上nvt系综进行驰豫。

驰豫完成把dump文件下载下来放到ovito里查看,发现结构出现问题。

合金内部出现了很多空洞,在正常情况下,合金不可能出现空洞。

一开始猜测可能是势函数问题,重新下载了一个势函数,上机运行,症状依旧。

既然不是势函数问题,那还是要分析一下原因。

在nvt系综下驰豫,box体积保持不变,原子在力的作用下抱团收缩,原子之间的距离变小,原子总体积小于box体积,从而产生空洞。

这种情况下就不能再用nvt进行驰豫。
npt系综根据原子总体积的变化自动调整box的尺寸,比较适合当前的体系,于是把nvt改为npt。

在npt下驰豫后发现合金在开始阶段有比较大的收缩,收缩之后空洞消失,得到比较致密的合金组织,完美解决问题。

多数情况下,用nvt系综进行驰豫即可,但部分体积变化比较大的体系,用npt系综进行驰豫效果会更好。

问题的答案又回到原点,假如别人问你如何选择系综,优雅地告诉他“根据具体情况而定”,这答案准没错。

更多lammps教程,扫码关注公众号:

lammps教程:实例讲解npt、nvt系综的选择相关推荐

  1. 工具类软件操作手册_全套广联达软件学习资料合集:教程+实例讲解+操作手册,一文搞定...

    要学习并掌握好算量软件,不仅仅需要软件知识,还需要各种造价相关知识,只有精通软件及造价知识,才能大大的提高算量的速度及准确度. 一个算量的高手掌握的不单是软件操作技术,还要掌握的是造价的基础知识,你要 ...

  2. CSS教程:实例讲解定位Position

    http://www.missyuan.com/thread-395406-1-1.html 1. position:static 所有元素的默认定位都是:position:static,这意味着元素 ...

  3. php线程教程,实例讲解php实现多线程

    我们首先来看一个示例 class execMulti { function MultithreadTest { exec("php static_massive.php > /dev/ ...

  4. lammps计算聚合物例子_LAMMPS中的系综(NPT/NVT)命令

    官方手册的翻译及补充,关于系综命令可见: fix NVT/NPT/NPH​lammps.sandia.gov 补充部分在段落结尾以"添注"形式补充. 语法: fix ID grou ...

  5. 讯搜 配置mysql_实例讲解帝国cms整合xunsearch(迅搜)教程

    实例讲解帝国cms整合xunsearch(迅搜)教程 来源:51酷酷 编辑:广告联盟评测网 时间:2017-11-28 16:50:18 帝国CMS安装迅搜_linux系统下 安装迅搜记录 目的:由于 ...

  6. 【编程实践】Git命令基础教程和代码实例讲解

    Git命令基础教程和代码实例讲解 Git是一个开源的分布式版本控制系统,用于敏捷高效地处理任何或小或大的项目.Git与常用的版本控制工具CVS.Subversion等不同,它采用了分布式版本库的方式. ...

  7. 编程python爬取网页数据教程_实例讲解Python爬取网页数据

    一.利用webbrowser.open()打开一个网站: >>> import webbrowser >>> webbrowser.open('http://i.f ...

  8. python饼状图教程_Python数据可视化:饼状图的实例讲解

    使用python实现论文里面的饼状图: 原图: python代码实现: # # 饼状图 # plot.figure(figsize=(8,8)) labels = [u'Canteen', u'Sup ...

  9. pythonappend教程_对python append 与浅拷贝的实例讲解

    在做Leetcode的第39题的时候,看到网上一个用递归的解法,很简洁.于是重写了一遍. class Solution(object): def combinationSum(self, candid ...

最新文章

  1. 思考: 什么时候需要disable MMU/i-cache/d-cache?
  2. mysql sqlite 分页查询_php基于SQLite实现的分页功能示例
  3. C# Lambda 和 匿名函数的GC总结
  4. 吴恩达《机器学习》学习笔记十二——机器学习系统
  5. 【HTML5】Web存储、通信、地理位置
  6. Java高级进阶学习资料!Java虚拟机的垃圾回收机制
  7. mysql 5.7.9 winx64_Windows 7 安装配置 mysql-5.7.17-winx64 方法-Fun言
  8. AD 10分钟画一块最简PCB
  9. GHOSTXPSP3电脑爱好者V9.9美化版
  10. 如何在你的Android手机上配置 Python 环境?
  11. python生成树状图_使用python的分层聚类树状图
  12. 计算机重启 ie 被改,ie被修改怎么办 ie被修改的解决方法【详解】
  13. BIM模型文件下载——江湖别墅Revit模型
  14. Android Studio 使用友盟进行多渠道打包
  15. Tesra超算网络,性价比最高的超算中心!
  16. C++校园疫情防控管理系统
  17. Java 8 字符串和时间相互转换
  18. table合并单元格 colspan(跨列)和rowspan(跨行)
  19. python人名最多统计_Python统计字符串中出现次数最多的人名
  20. MySQL面试:left join我要怎优化?

热门文章

  1. 支付项目:9、购物车模块
  2. 天下数据支招如何防范域名被劫持
  3. 随机取样已死,蓄水池抽样称王
  4. 输入三个字母,比较AISS码大小
  5. python3 beautifulsoup配合requests下载老司机图片
  6. 如何对市场进行深入了解,了解当前市场上的热销产品、消费者需求以及行业发展趋势?在哪里寻找专业报告、行业数据、市场分析文章等?
  7. 推荐一款超级好用的工具:uTools详解使用
  8. 互联网跟移动互联网_互联网如何变坏
  9. [databinding] ... ... conflicts with another tag that has the same ID
  10. 素数个数求解与素数的判定