24小时热门版块排行榜    

查看: 1230  |  回复: 0

x39988

新虫 (初入文坛)

[交流] 分子动力学的 ADMP 该如何设定温度

我目前尝试用Gaussian 中 ADMP进行 分子动力学的模拟

输入档如下

#b3lyp/6-31++g(d,p) admp=maxpoints=10 IOP(1/80=1000000,1/81=1,1/82=500,1/89=1)

title

1,1
C                  0.53346723    0.05849264    0.31010108
C                 -1.60293444    0.38889993    0.55333428
C                 -1.34189470   -0.54936854   -0.60971732
O                 -0.45541151    0.17114853    1.36443082
O                 -0.03254012   -0.16169929   -1.00646069
O                  1.77146835    0.14411290    0.51889036
H                 -1.63352428    1.40194629    0.21025800
H                 -2.51862395    0.17843365    1.06530760
H                 -2.05973139   -0.45187708   -1.39718535
H                 -1.34689706   -1.56487646   -0.27264368
Li                 2.98584505   -0.01568228   -0.92050351

因为我再ADMP的关键词中找不到关于温度的关键词
所以我参考了 ( http://muchong.com/bbs/viewthread.php?tid=5431433 )
这篇帖子中提供的设定


我有回去官网找IOP各项设定代表的意思
( http://www.gaussian.com/g_tech/g_iops/ov1.htm#1-80 )
例如
IOp(1/81)

Nuclear KE thermostat in ADMP -- temperature is checked and scaled every IOp(81) steps.

IOp(1/82)

L121: Temperature for nuclear KE thermostat.


其中IOP(1/82)=350的这个部分 按照我的逻辑应该是指350kelvin

但跑完之后 我从输出档内看到

Thermostat chosen for nuclear motion. Option =  1
   Thermostat temperature =  0.016626 hartrees
   Thermostat range =  0.016578 to  0.016673 hartrees

在里面写的温度 不是用kelvin标示
但是我把hartree经过换算成kelvin后 仍然不是350kelvin

想请问这究竟是怎么回事 ?

我已经尽量把所有信息都提供上来,还希望有人能替我解答解答

或者有人能告诉我ADMP该如何设定他的温度 ?

[ Last edited by x39988 on 2013-7-25 at 19:29 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 x39988 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +8 Kittylucky 2026-08-27 8/400 2026-08-27 18:33 by raiy123
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见