| 查看: 1054 | 回复: 4 | ||
[求助]
lammps计算室温下Mo的能量有疑问 已有1人参与
|
|
# LAMMPS-Mo units metal boundary p p p atom_style atomic lattice bcc 3.147 region Mo block 0 8 0 8 0 8 create_box 1 Mo create_atoms 1 region Mo group Mo region Mo mass 1 95.94 pair_style morse 6.0 pair_coeff 1 1 0.7702 1.4341 3.0122 neighbor 2.0 bin neigh_modify every 1 delay 0 check yes timestep 0.001 velocity Mo create 299 825577 dist gaussian rot yes thermo 100 thermo_style custom step temp enthalpy etotal pe ke press vol minimize 1e-6 1e-6 100000 100000 fix 1 Mo nve reset_timestep 0 run 20000 unfix 1 variable P equal pe variable T equal temp variable Press equal press variable V equal vol variable E equal etotal variable N equal atoms compute rdf Mo rdf 100 1 1 fix 2 Mo ave/time 100 10 1000 c_rdf[1] c_rdf[2] c_rdf[3] file Mo.rdf mode vector fix 1 Mo nvt 270 300 0.1 drag 0.2 run 50000 ---------------------------------------------------------------------------------------------------------------------------------------------------- ------------------------------------------------------------------------------------------------------------------------------------------------------ 因为我找到了室温也就是298K下Mo的比热容数据，所以想通过比热容公式 c=dE/(dT*m)来计算。经过各种单位换算得到c=(dE*4.18*10^3)/(dT*95.94*10^-3)，我想问一下，我上面的脚本需要怎么改才能实现我的目的？希望各位大神畅所欲言 @老虎大王 @月只蓝 |
回复此楼» 猜你喜欢
纯娱乐，不喜欢勿喷
已经有8人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
娱乐
已经有3人回复
-
面上提前没消息，有中的吗
已经有17人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】晶格常数大了以后 VASP的计算值准确吗
已经有5人回复
- 【讨论】lammps中NVE和NVT定温有什么区别
已经有7人回复
茹此无声
木虫 (小有名气)
- 应助: 12 (小学生)
- 金币: 1438.8
- 红花: 13
- 帖子: 244
- 在线: 65小时
- 虫号: 2065582
- 注册: 2012-10-16
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
ytudou2楼2014-01-19 00:19:50
3楼2014-01-19 09:50:25
4楼2014-02-21 17:08:29
茹此无声
木虫 (小有名气)
- 应助: 12 (小学生)
- 金币: 1438.8
- 红花: 13
- 帖子: 244
- 在线: 65小时
- 虫号: 2065582
- 注册: 2012-10-16
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2014-02-22 02:22:55
50