24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 865  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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~
已阅   回复此楼   关注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的回帖
查看全部 5 个回答

茹此无声

木虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
月只蓝: 金币+1, 鼓励交流 2014-01-19 09:42:56
修改fix computer 就行

» 本帖已获得的红花(最新10朵)

人类一思考,上帝就发笑。
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +10 junqihahaha 2026-03-26 11/550 2026-03-27 01:15 by yhmsz
[考研] 329求调剂 +6 钮恩雪 2026-03-25 6/300 2026-03-26 22:50 by fmesaito
[考研] 一志愿北化求调剂 +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv风雷
[考研] 327求调剂 +7 prayer13 2026-03-23 7/350 2026-03-26 20:48 by 不吃魚的貓
[考研] 281求调剂 +3 亚克西good 2026-03-26 5/250 2026-03-26 19:48 by 不吃魚的貓
[考研] 279求调剂 +6 红衣隐官 2026-03-21 6/300 2026-03-26 18:32 by 不吃魚的貓
[考研] 271求调剂 +6 生如夏花… 2026-03-22 6/300 2026-03-26 16:48 by 张凯十八号
[考研] 机械学硕310分,数一英一,一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
[考研] 资源与环境 调剂申请(333分) +9 holy J 2026-03-21 9/450 2026-03-26 15:47 by 161765490
[考研] 材料与化工328分调剂 +6 。,。,。,。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:03 by cmz0325
[考研] 285求调剂 +3 AZMK 2026-03-24 3/150 2026-03-25 12:23 by userper
[考研] 求调剂 +5 林之夕 2026-03-24 5/250 2026-03-24 17:16 by dick_runner
[考研] 300求调剂,材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
信息提示
请填处理意见