24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

809011916

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 316.4
帖子: 55
在线: 42.3小时
虫号: 1501800
注册: 2011-11-20
性别: GG
专业: 化工热力学和基础数据

[求助] 【Lammps】LJ单位的交流已有1人参与

有关Lammps的LJ单位的在mannual中的解读如下：
mass = mass or m
distance = sigma, where x* = x / sigma
time = tau, where t* = t (epsilon / m / sigma^2)^1/2
energy = epsilon, where E* = E / epsilon
velocity = sigma/tau, where v* = v tau / sigma
force = epsilon/sigma, where f* = f sigma / epsilon
torque = epsilon, where t* = t / epsilon
temperature = reduced LJ temperature, where T* = T Kb / epsilon
pressure = reduced LJ pressure, where P* = P sigma^3 / epsilon
dynamic viscosity = reduced LJ viscosity, where eta* = eta sigma^3 / epsilon / tau
charge = reduced LJ charge, where q* = q / (4 pi perm0 sigma epsilon)^1/2
dipole = reduced LJ dipole, moment where *mu = mu / (4 pi perm0 sigma^3 epsilon)^1/2
electric field = force/charge, where E* = E (4 pi perm0 sigma epsilon)^1/2 sigma / epsilon
density = mass/volume, where rho* = rho sigma^dim

大体上知道了LJ相关单位是怎么回事，例如一个体系中有A,B,C,D四种原子，那么会有多个sigma和epsilon，利用混合规则可以得到原子间sigma和epsilon，此时如果规定了其中一个sigma为1，那么另外的sigma就会随之改变，由此定义了相关的pair参数。有两个问题：（1）以此理解类推，如果fix 1 all nvt 1.0 1.0 0.1，即是指温度保持不变做0.1T下的tdamp进行模拟，那此时的模拟是如何体现具体what温度下的平衡状态呢？这时T*=T*Kb / epsilon的意义是什么？难道要实现温度的意义就要定义real形式？（感觉不会是这样，应该还是有别的理解方法吧，个人觉得这个问题很基础但又很关键，希望路过的大神解释下）（2）之前根据大家的方法输出了data文件，此时如果原子定义成lj类型，那data文件中不同单位的坐标、参数岂不是都要随着进行改变？想知道大家如何解决这个问题的？

希望大家不吝赐教，谢谢！

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

选择最喜欢的，专注到无情。一个，就够了。

1楼 2015-09-07 04:10:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

alundilong

木虫 (正式写手)

应助: 129 (高中生)
金币: 2019.1
红花: 18
帖子: 309
在线: 56小时
虫号: 595468
注册: 2008-09-08
性别: GG
专业: 传热传质学

【答案】应助回帖

这里的温度应该是：T* = T Kb / epsilon。楼主应该多打了一个星号。

在转化为真实温度的时候，kb和epsilon都是知道的，对吧？

我想，楼主的问题集中在“规定了其中一个sigma为1”上面。

请楼主思考，为何平白无故的可以“规定了其中一个sigma为1”这个1怎么来的？
=======================================
现在看一个例子（都在lj作用势下）：
argon's epsilon = 0.2381 Kcal/mol
argon's sigma = 3.405 A

copper's epsilon = 9.566 Kcal/mol
copper's sigma = 2.2277 A

现在选择argon的sigma作为基准长度，argon的epsilon作为基准能量；
于是，argon的无量纲（units lj）下的长度是多少？argon's sigma/ argon's sigma = = 3.405/3.405 = 1 (到这里，我们就知道1是怎么来的了)
以此类推，copper的无量纲下的长度是 copper's sigma/argon's sigma = 2.2277/3.405

所有坐标也是以此类推
能量以此类推。。。

希望有解决楼主的困惑。