学习lammps模拟,不能只看教程和手册,就算把手册背下来,遇到实际问题也不一定能够解决,必须要有实践过程。
最好的方法就是找到已经发表的论文,按照论文中介绍的参数和模拟步骤,亲自写一个lammps代码进行模拟,并与论文的模拟结果进行对比。
对于新手来说,这个过程肯定会比较痛苦,可能会遇到各种各样意想不到的问题,当逐渐解决这些问题之后,你会发现你的水平已经有了很大的提高,回过头来再看看曾经遇到的那些难题已能轻松应对。
本号除了陆续推出lammps基础教程外,每周尽量推出1-2篇实例教程,教程以已经发表的论文为例,复现论文的模拟过程,共大家学习和参考。
建议先下载论文,仔细阅读论文,从文中找出模拟参数和模拟步骤,例如模型尺寸、温度大小,采用npt系综还是nvt系综,等等。
自己尝试编写代码,最后再对比我给出案例代码,即时反馈才是最有效的学习方式。
如果在这个过程中有不明白的问题,可以加微信公众号(lammps加油站)一起讨论,互相进步。
这篇推文以《Atomic simulations of the formation of twist grain boundary and mechanical properties of graphene/aluminum nanolaminated composites》为例,尝试复现文章铝/石墨烯复合物结晶过程。
模拟的结晶过程如下图:
模拟势能变化曲线对比如下:
源代码下载请关注微信公众号获取:lammps加油站