| 查看: 2040 | 回复: 2 | ||
809011916铜虫 (小有名气)
|
[求助]
【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文件中不同单位的坐标、参数岂不是都要随着进行改变？想知道大家如何解决这个问题的？ 希望大家不吝赐教，谢谢！ |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有10人回复
-
基金中了
已经有13人回复
-
关于filecode
已经有6人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
jiezhou5322
银虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2072.4
- 散金: 5
- 红花: 1
- 帖子: 390
- 在线: 332.6小时
- 虫号: 3093434
- 注册: 2014-03-27
- 专业: 岩土力学与岩土工程
2楼2016-03-08 21:30:14
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 所有坐标也是以此类推 能量以此类推。。。 希望有解决楼主的困惑。 |
3楼2016-03-09 00:23:03
100