24小时热门版块排行榜    

查看: 7106  |  回复: 24

chinaboy798

银虫 (小有名气)


[交流] 学习lammps半年的部分经验

开头:想到哪里就说到哪里,新手借鉴,高手批评指正。

1.lammps的动力学计算的原理是牛顿运动方程,所以在初学lammps的时候,建议阅读一下《理论力学》,另外,《热力学和统计力学》,“分子动力学的原理相关书籍”也需要看一下。
2.安装lammps,最好在linux系统下使用,因为随后走大规模计算的时候一般都是linux的平台,何不以linux开始?所以要掌握一些linux系统的基本使用常识,网上相关知识也很多。
3.学会编译lammps,对于如何编译网上帖子很多,关键在于要预安装fftw和mpich2两个软件,而且要选择这两个软件合适的版本才能编译成功。另外详细的我编译成功的经验我以后再发出来分享。
4.window版本的lammps可以用,但是用来跑跑小模型,测试等。
5.数据处理软件,一般会用oringin,matlab,或者用其他计算机语言编程处理数据,可视化上一般会使用atomeye,VMD,OVITO等。具体怎么分析数据,可以查找自己问题的相关文献去测试。所以前面提到的一些也得学会使用。
6.安装好程序后,首先要看的,就是手册了,lammps的手册虽然看起来有1200多页,但是实际需要重点关心的也就400页左右,第1,2,3,6章强烈建议多看几遍,其他的部分等熟练掌握了再看不迟。命令部分我写过另外一篇帖子,http://muchong.com/bbs/viewthread.php?tid=8149677,里面带星号的命令多看几遍,都是非常常用的命令。另外compute,fix这两个命令需要额外关注,连manual自己都把这两个命令单列出来了,可想而知他俩的重要性。
7.conpute命令是用来定义计算方式的,比如怎么算温度,压力等;fix命令式用来加载其他限制条件的,比如系综,力等。而且fix命令的加载是全时间的,代表的是只要你不说停,这个命令的效果就一直存在。
8.基于牛顿运动方程,我们得到的最原始的信息都是每一个原子的位置,速度,受力(加速度)。需要我们去调控的,也是原子的位置,速度,受力。所有的命令的出发点以及核心,我们力所能及的,都是去定义这三个中某一个或多个的初值或者控制他们的变化方式。
9.lammps的学习的来源主要有,基本理论书籍,lammps的使用手册(英文版的),lammps的官网,论坛(有一个分子模拟论坛,以及这里小木虫),以及lammps官网上附加的其他相关网站,网上公开的入门级的input文件,lammps的QQ交流群,还有一个自己课题的相关的用lammps的文献。这几部分应该说在学习的时候都应该要关注到。
10.重点提一下官网,特别是初学的时候,里面有大量的介绍性的入门资料,当然,manual也在里面。
11.使用lammps的最终目的,就是学会写出一个合理的脚本in文件,所以,学习的时候就应该落到这里,看再多,还是要试试才知道。最好自己亲手写一个小的能跑的in文件,慢慢根据其他的资料摸索,尝试其中命令对结果的影响,毕竟手册的英文的,有时候理解难免有偏颇之处,去试,才是最让自己信服的结果。
12.更进一步,就是知了lammps的核心原理,懂一些c++,自己到lammps的源程序中看看它到底是怎么工作的,如果还能提出自己的一些更好的原理建议,甚至自己改编,就最好不过了。
13.上述的只是lammps部分的需要关注的问题,另外还有两个问题是初始模型的构建,以及势函数。对于复杂的问题,这两块儿需要另外一些讨论,复杂的建模一般使用MS(需要购买版权),或者自己编程。如果需要处理的问题势函数没有人提到过,或者别人的势函数不适用你的问题,势函数的拟合问题就略微有点复杂,这里也不再具体讨论。
14.因此,在计算lammps的时候,至少需要一个in文件,lammps的程序,如果势函数需要另外的指定势函数,那么还需要一个势函数文件;如果模型比较复杂,lammps自带的命令不能构建,就还需要一个你构建的模型的额外的data文件。
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学 MD分子动力学 分子模拟 好东西啊
分子动力学模拟过程 lammps LAMMPS怕不怕,不怕! 分子动力学-初级问题资源
lammps

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

992648352

银虫 (小有名气)


你好,请问一下lammps能不能固定分子的键长和 键角,并且我想特指某一对原子间的非键相互作用参数呢?
12楼2014-12-03 14:10:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chinaboy798

银虫 (小有名气)


引用回帖:
12楼: Originally posted by 992648352 at 2014-12-03 14:10:26
你好,请问一下lammps能不能固定分子的键长和 键角,并且我想特指某一对原子间的非键相互作用参数呢?

高分子我还没关注过,非常抱歉!
13楼2014-12-04 13:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ktwyw

木虫 (小有名气)


引用回帖:
12楼: Originally posted by 992648352 at 2014-12-03 14:10:26
你好,请问一下lammps能不能固定分子的键长和 键角,并且我想特指某一对原子间的非键相互作用参数呢?

可以。

[ 发自手机版 http://muchong.com/3g ]
15楼2014-12-17 05:57:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liwl926

金虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
非常感谢你的分享,我是lammps新手,请教一下lammps可以用来优化晶胞参数么?minimize命令还需要什么样的参数或命令才能使优化过程中晶胞参数也向更合理的方向发生变化?多谢!
16楼2014-12-17 09:57:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hust_fd

新虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
您好,请问您是如何不用lammps自带的建模命令建立data文件的? 我使用MS然后转化的,但是中间出现了很多问题,能否给我一些建议,我对这个伤透了脑筋!
18楼2014-12-19 15:09:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖仅楼主可见
20楼2015-08-25 22:08:27
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

douniuniu

新虫 (初入文坛)


21楼2015-09-15 15:47:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

214337458

捐助贵宾 (初入文坛)


非常感谢!写得好!
22楼2015-09-16 21:06:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Awcse

新虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
您好,我做的是第一性原理的MD,最近用VMD处理结果的时候遇到一个问题想计算RMSD,却总是提示我未选定原子,不知道您是否有这方面的经验能告知一下?
24楼2015-12-24 09:26:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr. J

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,请问如何添加一个新的力场,求解答

发自小木虫IOS客户端
25楼2016-05-09 23:18:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
dmbb2楼
2014-12-03 12:57   回复  
chinaboy798(金币+1): 谢谢参与
nono20093楼
2014-12-03 12:59   回复  
chinaboy798(金币+1): 谢谢参与
2014-12-03 13:07   回复  
chinaboy798(金币+1): 谢谢参与
xiejf5楼
2014-12-03 13:26   回复  
chinaboy798(金币+1): 谢谢参与
1
2014-12-03 13:29   回复  
chinaboy798(金币+1): 谢谢参与
gorgan7楼
2014-12-03 13:31   回复  
chinaboy798(金币+1): 谢谢参与
tzynew8楼
2014-12-03 13:31   回复  
chinaboy798(金币+1): 谢谢参与
neu2349楼
2014-12-03 13:42   回复  
chinaboy798(金币+1): 谢谢参与
time8810楼
2014-12-03 13:52   回复  
chinaboy798(金币+1): 谢谢参与
假大空11楼
2014-12-03 14:09   回复  
chinaboy798(金币+1): 谢谢参与
fenqin14楼
2014-12-16 21:37   回复  
三丁丁17楼
2014-12-18 00:42   回复  
2015-08-25 14:03   回复  
2015-09-21 15:21   回复  
相关版块跳转 我要订阅楼主 chinaboy798 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 72/3600 2026-08-16 13:15 by 医学老男孩
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见