请问用分子动力学方法计算氩晶体的行为的时候,咋计算氩晶体压强
http://lammps.sandia.gov/doc/compute_pressure.html看页面内公式及文献
■网友
那个虽然已经帮不上什么忙了,随便写写吧。刚好前一段在写分子动力学的算法,当时我做第一件事就是算LJ 势能下的氩气。问题不外乎出现在reduced unit上,里面有一部分单独有讲。强推 Frenkel Molecualr dynamics这本书里面有很详细的讲解,不过里面的code是用Fortran写的。但问题是 这本书NVE 的算法给的很详细,假如想写NVT 或者NPT的体系的话需要另找paper。。
■网友
占坑。回去给你发一个paper。里面有比较详细的思路。你们的Loneard—Jones势能是怎么算的?用的理论推导还是编程?本人恰好写过类似的python代码。不过其实你谷歌,uiuc phy 485的官网,六个作业中的第二个,就是这套问题的解决方案。可以给我看下思路么?有这个势能以后求解很多物理常量应该非常容易。目测你们Loneard—Jones势能的计算或者边界条件(这个里头边界的处理比较tricky)搞错了。哥们顺便说一下工作单位?目测是清北科浙复交或者国外top。以及如果是本科生请接受我的膜拜。
■网友
这个不难,我博士阶段模拟过这个。
建议你用这个公式。其中
是体系位力的时间平均。
这些内容在很多分子动力学的教材上都有的。
■网友
没给算压强的方式那你是怎么算的=。=
假设你压强计算方式的是正确的,那你可以考虑一下你的参数是不是有问题(温度和LJ参数,通常LJ势能只能模拟理想气体)
■网友
我其实没搞懂问题。。。压强大就是体积小的呗。。如果用的nve的话。。
【请问用分子动力学方法计算氩晶体的行为的时候,咋计算氩晶体压强】 所以你是怎么判断能量稳定的? 你就先用npt等energy converge以后再nve就好啦。再搞不定就问你们上课的老师呗。。
推荐阅读
- 银行高管的薪酬结构图表最新的,请问在哪里可以收集
- 在上海陆家嘴附近工作,请问去哪边租房会比较好
- 坐标合肥,请问在哪里能捕捉到程序员这种生物他们大都出现在哪里呢
- 孕妇到了预产期,小孩还没有入盆,请问是该等下去还是直接剖腹产
- 想做手机网游代充,请问要去哪里找货源
- 为啥无听过电子钱咭被不法分子破解?
- 请问有哪些机器人大赛是比较权威的
- 请问国内最好的数据恢复公司是哪个朋友结婚视频不小心被误删,急求帮忙
- 请问杨毅微博中的这两人是谁
- 通信专业大二在读,以后想从事软件方向。请问现在应该准备些啥
