| 查看: 2077 | 回复: 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文件中不同单位的坐标、参数岂不是都要随着进行改变？想知道大家如何解决这个问题的？ 希望大家不吝赐教，谢谢！ |
回复此楼» 猜你喜欢
2026年叶企孙基金
已经有7人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
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
30