应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps计算室温下Mo的能量有疑问
51/1返回列表
查看: 1085  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ytudou

银虫 (小有名气)

[求助] 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)，我想问一下，我上面的脚本需要怎么改才能实现我的目的？希望各位大神畅所欲言

@老虎大王  @月只蓝
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
要给力啊，MS~LAMMPS~
1楼 2014-01-17 13:02:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茹此无声

木虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
月只蓝: 金币+1, 鼓励交流 2014-01-19 09:42:56
修改fix computer 就行
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

ytudou
人类一思考，上帝就发笑。
2楼2014-01-19 00:19:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytudou

银虫 (小有名气)
引用回帖:
2楼: Originally posted by 茹此无声 at 2014-01-19 00:19:50
修改fix computer 就行

那个，不好意思，能不能进一步问一下怎么改？有没有模板？或者近似的例子，因为本人编程能力实在太弱，所以可能只能做到修改……麻烦大神再指点一二~不胜感激~~
赞一下 回复此楼
要给力啊，MS~LAMMPS~
3楼2014-01-19 09:50:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytudou

银虫 (小有名气)
★
送红花一朵
月只蓝: 金币+1, 鼓励交流！ 2014-02-21 18:12:58
引用回帖:
2楼: Originally posted by 茹此无声 at 2014-01-19 00:19:50
修改fix computer 就行

关于lammps有不明白的地方，冒昧请教一下(^_^)，望不吝赐教～
dump 10 all xyz 500 dump.xyz
这句中10，500分别指的什么？
还有就是，这句输出的是谁的坐标？是全部原子一段时间内的坐标变换？还是什么？
真的困惑，望大神解救啊～

[ 发自小木虫客户端 ]
赞一下 回复此楼
要给力啊，MS~LAMMPS~
4楼2014-02-21 17:08:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茹此无声

木虫 (小有名气)

【答案】应助回帖

★
月只蓝: 金币+1, 感谢指导！ 2014-02-22 11:04:40
引用回帖:
4楼: Originally posted by ytudou at 2014-02-21 17:08:29
关于lammps有不明白的地方，冒昧请教一下(^_^)，望不吝赐教～
dump 10 all xyz 500 dump.xyz
这句中10，500分别指的什么？
还有就是，这句输出的是谁的坐标？是全部原子一段时间内的坐标变换？还是什么？
真的 ...

每隔500步输出所有原子的坐标，10是个id

[ 发自小木虫客户端 ]
赞一下 回复此楼
人类一思考，上帝就发笑。
5楼2014-02-22 02:22:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ytudou 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:12 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:59 by winsaint
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:53 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗？ +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
关闭
请填处理意见
关闭