路径积分分子动力学(PIMD)方法提供了一种方便的方法来计算凝聚相系统的量子力学结构和热力学性质,但代价是在系统的物理振动之上引入了一组额外的高频简正模式。对如此宽的频率范围进行有效采样是一个相当大的恒温挑战。在这里,我们介绍了一个简单的随机路径积分朗之万方程(PILE)恒温器,它利用了自由路径积分法向模频率的分析知识。我们还应用了最近开发的基于广义朗之万方程(GLE)的有色噪声恒温器,该恒温器可自动实现类似的频率优化采样。将这些恒温器的采样效率与更传统的Nosé–Hoover链(NHC)恒温器进行比较,以了解液态水和氢-钯系统的一些物理相关特性。几乎在所有情况下,新的PILE恒温器都能与NHC恒温器一样工作,同时可以实现更高效的计算。GLE恒温器也被证明非常稳健,在所有考虑的情况下都能提供接近最佳的采样效率。我们怀疑,这些简单的随机恒温器将因此在许多未来的PIMD模拟中找到有用的应用。

1
钢筋混凝土。
费曼
A.R.公司。
希布斯
,
量子力学与路径积分
(
麦格劳-希尔
,
纽约
,
1965
).
2.
钢筋混凝土。
费曼
,
统计力学
(
出版商
,
纽约
,
1972
).
三。
D。
钱德勒
P.G.公司。
Wolynes公司
,
化学杂志。物理学。
74
,
4078
(
1981
).
4
J.A.公司。
巴克
,
化学杂志。物理学。
70
,
2914
(
1979
).
5
M。
帕里内洛
答:。
拉赫曼
,英寸
量子问题中的蒙特卡罗方法
,编辑人
M.H.先生。
卡罗斯
(
雷德尔
,
多德雷赫特
,
1984
),第页。
105
.
6
M。
帕里内洛
答:。
拉赫曼
,
化学杂志。物理学。
80
,
860
(
1984
).
7
相对湿度。
霍尔
B.J.公司。
伯尔尼
,
化学杂志。物理学。
81
,
3641
(
1984
).
8
高压断路器。
安徒生
,
化学杂志。物理学。
72
,
2384
(
1980
).
9
B.J.公司。
伯尔尼
D。
Thirulamai公司
,
每年。物理版。化学。
37
,
401
(
1986
).
10
美国。
没有
,
化学杂志。物理学。
81
,
511
(
1984
).
11
W.G.公司。
胡佛
,
物理学。版本A
31
,
1695
(
1985
).
12
G·J。
马蒂娜
,
M·L。
克莱因
、和
机械工程师。
塔克曼
,
化学杂志。物理学。
97
,
2635
(
1992
).
13
机械工程师。
塔克曼
,
B.J.公司。
伯尔尼
,
G·J。
马蒂纳
、和
M.L.公司。
克莱因
,
化学杂志。物理学。
99
,
2796
(
1993
).
14
G.公司。
布西
M。
帕里内洛
,
物理学。版本E
75
,
056707
(
2007
).
15
G.公司。
布西
,
D。
多纳迪奥
、和
M。
帕里内洛
,
化学杂志。物理学。
126
,
014101
(
2007
).
16
M。
塞里奥蒂
,
G.公司。
布西
、和
M。
帕里内洛
,
物理学。修订稿。
102
,
020601
(
2009
).
17
M。
塞里奥蒂
,
G.公司。
布西
、和
M。
帕里内洛
,
化学杂志。理论计算。
6
,
1170
(
2010
).
18
H。
Trotter公司
,
程序。美国数学。Soc公司。
10
,
545
(
1959
).
19
洛杉矶。
舒尔曼
,
路径集成技术及其应用
(
多佛
,
纽约
,
2005
).
20
M.F.公司。
赫尔曼
,
E.J.公司。
布鲁斯金
、和
B.J.公司。
伯尔尼
,
化学杂志。物理学。
76
,
5150
(
1982
).
21
D.M.博士。
塞佩利
,
修订版Mod。物理学。
67
,
279
(
1995
).
22
爱尔兰共和国。
克雷格
D.E.博士。
马诺洛波洛斯
,
化学杂志。物理学。
121
,
3368
(
2004
).
23
B.J.公司。
布拉姆斯
D.E.博士。
马诺洛波洛斯
,
化学杂志。物理学。
125
,
124105
(
2006
).
24.
J·M·。
桑兹·塞尔纳
M.P.公司。
卡尔沃
,
数值哈密顿问题
(
查普曼和霍尔
,
伦敦
,
1994
).
25
B。
莱姆库勒
美国。
帝国
,
哈密顿动力学模拟
(
剑桥大学出版社
,
剑桥
,
2004
).
26
L。
Verlet公司
,
物理学。版次。
159
,
98
(
1967
).
27
在没有恒温器的情况下,通过将动量保持在正常模式表示中,并将力转换为该表示,以便在方程式中进行更新,每个时间步长的傅里叶变换次数可以从每自由度四次减少到三次。(21)–(25); 这也消除了等式中PILE恒温器中傅里叶变换的需要。(27)和(29).
28
答:。
里奇
G.公司。
奇科蒂
,
摩尔物理学。
101
,
1927
(
2003
).
29
R。
茨万齐格
,
非平衡统计力学
(
牛津大学出版社
,
牛津
,
2001
).
30
C.W.公司。
加德纳
,
随机方法手册
,第3版。(
施普林格
,
柏林
,
2003
).
31
M。
塞里奥蒂
,
G.公司。
布西
、和
M。
帕里内洛
,
物理学。修订稿。
103
,
030603
(
2009
).
32
G·J。
马蒂娜
,
机械工程师。
塔克曼
,
D.J.博士。
托拜厄斯
、和
M.L.公司。
克莱因
,
摩尔物理学。
87
,
1117
(
1996
).
33
美国。
总会计师。
沃斯人
,
化学杂志。物理学。
107
,
9514
(
1997
).
34
G.公司。
布西
M。
帕里内洛
,
计算。物理学。Commun公司。
179
,
26
(
2008
).
35
注册会计师。
库哈斯基
P.J.公司。
罗斯基
,
化学杂志。物理学。
82
,
5164
(
1985
).
36
答:。
瓦尔奎斯特
B.J.公司。
伯尔尼
,
化学。物理学。莱特。
117
,
214
(
1985
).
37
美国。
哈伯逊
,
T.E.公司。
马克兰
、和
D.E.博士。
马诺洛波洛斯
,
化学杂志。物理学。
131
,
024501
(
2009
).
38
T.E.公司。
马克兰
D.E.博士。
马诺洛波洛斯
,
化学杂志。物理学。
129
,
024105
(
2008
).
39
T.E.公司。
马克兰
D.E.博士。
马诺洛波洛斯
,
化学。物理学。莱特。
464
,
256
(
2008
).
40.
S.L.公司。
阿德勒
,
物理学。版次D
23
,
2901
(
1981
).
41
X·W。
,
J.A.公司。
齐默尔曼
,
B.M.公司。
Wong(王)
、和
J·J。
霍伊特
,
J.马特。物件。
23
,
704
(
2008
).
42
合适的选择τ0可以在计算之前为GLE恒温器自动确定,如下所示。自由环聚合物的最大内部频率为ω最大值=2/(βn个)。由于我们使用的GLE适合于4个数量级的频率范围(见图。1),这意味着我们应该ω0在等式中(41)将成为ω0=102ω最大值=0.02/(βn个),给予τ0=1/(2ω0)=25βn个.对于298 K的液态水n个=32环状聚合物珠τ0=20英尺对于钯中的氢原子,在350 K下,有十个珠子,它给出τ0=55英尺; 检查图。2–5表明这些值τ0确实在每种情况下都能为GLE提供近乎最佳的采样。
43
在我们对液态水的模拟中,发现使用PILE恒温器比NVE模拟引入了5%的开销。GLE恒温器的开销(需要生成更多随机数并执行矩阵-向量乘法)约为30%。NHC方程的数值积分需要计算每个自由度和链变量的几个指数。然后,为了避免守恒量的任何漂移,需要进行多时间步长积分,这进一步放大了这一努力,我们发现这会导致我们模拟中250%的总开销。然而,我们应该强调的是,只有当作用在环状聚合物每个珠上的物理力能够以极低的成本进行评估时,才能获得如此大的恒温开销,就像使用环状聚合物收缩技术的经验水模型一样。从头算PIMD模拟NHC开销通常只占计算总成本的一小部分。
您当前无权访问此内容。