24小时热门版块排行榜    

查看: 1659  |  回复: 7

skplin

铜虫 (小有名气)

[求助] 分子动力学我设置的恒温为什么降温了 已有1人参与

ISTART=0
ICHARG=2
NPAR=4

PREC=Accurate
ENCUT=300
LREAL=Auto
ISMEAR=-1
SIGMA=0.05
ISPIN=2

IBRION=-1
ISIF=2
NSW=0
POTIM=0.5

APACO=16
NPACO=256

TEBEG=1000
TEEND=1000
SMASS=1

LWAVE=.FALSE.
LCHARG=.TRUE.
哪里错了?????

分子动力学我设置的恒温为什么降温了
$D{68I7D_7G8302HUJ4HKH6.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
skplin: 金币+10 2016-05-18 19:52:14
如果体系是优化过的,温度大概会降到你设置的初始值的一半然后慢慢回升到你的heatbath的温度.
看看virial theorem就明白了.
2楼2016-05-10 18:09:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

【答案】应助回帖

还有IBRION是不是-1吧,你仔细看看manual
3楼2016-05-10 18:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skplin

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by luzihen at 2016-05-10 18:09:50
如果体系是优化过的,温度大概会降到你设置的初始值的一半然后慢慢回升到你的heatbath的温度.
看看virial theorem就明白了.

该怎么办
smass如果等于2会有好转吗?
4楼2016-05-18 19:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

引用回帖:
4楼: Originally posted by skplin at 2016-05-18 19:52:51
该怎么办
smass如果等于2会有好转吗?...

读读手册再来问,IBRION先设置成0试试再说。你有没有看过轨迹,原子核动了没有。

6.21.1 IBRION=-1
No update; ions are not moved, but NSW outer loops are performed. In each outer loop the electronic degrees of freedom are
re-optimized (for NSW>0 this obviously does not make much sense, except for test purposes). If no ionic update is required
use NSW=0 instead.
5楼2016-05-20 15:48:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skplin

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by luzihen at 2016-05-20 15:48:58
读读手册再来问,IBRION先设置成0试试再说。你有没有看过轨迹,原子核动了没有。

6.21.1 IBRION=-1
No update; ions are not moved, but NSW outer loops are performed. In each outer loop the electronic d ...

分子动力学怎么看轨迹,,,手册上说分子动力学应该选IBRION=0   然后原子动了
6楼2016-05-21 10:18:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skplin

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by luzihen at 2016-05-20 15:48:58
读读手册再来问,IBRION先设置成0试试再说。你有没有看过轨迹,原子核动了没有。

6.21.1 IBRION=-1
No update; ions are not moved, but NSW outer loops are performed. In each outer loop the electronic d ...

IBRION=0   nsw= 500   smass= 1   温度都是1000   但是温度变化范围还是很大
7楼2016-05-21 10:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

引用回帖:
7楼: Originally posted by skplin at 2016-05-21 10:20:44
IBRION=0   nsw= 500   smass= 1   温度都是1000   但是温度变化范围还是很大...

什么叫大?
有波动很正常,没有fluctuation就不是NVT系综了
8楼2016-05-21 14:12:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 skplin 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[基金申请] fileCode有新解读? +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见