24小时热门版块排行榜

返回列表

x39988

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 61.5
帖子: 2
在线: 12.5小时
虫号: 2545853
注册: 2013-07-15

[交流] 分子动力学的 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 ]

回复此楼