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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps计算室温下Mo的能量有疑问
51/1返回列表
查看: 1108  |  回复: 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一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:49 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:25 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:12 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:57 by Aj1rhIDL5ixY
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:13 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:01 by Aj1rhIDL5ixY
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:00 by Aj1rhIDL5ixY
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
关闭
请填处理意见
关闭