We introduce a novel simple algorithm for thermostatting path integral molecular dynamics (PIMD) with the Langevin equation. The staging transformation of path integral beads is employed for demonstration. The optimum friction coefficients for the staging modes in the free particle limit are used for all systems. In comparison to the path integral Langevin equation thermostat, the new algorithm exploits a different order of splitting for the phase space propagator associated to the Langevin equation. While the error analysis is made for both algorithms, they are also employed in the PIMD simulations of three realistic systems (the H2O molecule, liquid para-hydrogen, and liquid water) for comparison. It is shown that the new thermostat increases the time interval of PIMD by a factor of 4–6 or more for achieving the same accuracy. In addition, the supplementary material shows the error analysis made for the algorithms when the normal-mode transformation of path integral beads is used.

1.
R. P.
Feynman
,
Phys. Rev.
91
(
6
),
1291
1301
(
1953
).
2.
D.
Chandler
and
P. G.
Wolynes
,
J. Chem. Phys.
74
(
7
),
4078
4095
(
1981
).
3.
D. M.
Ceperley
,
Rev. Mod. Phys.
67
(
2
),
279
355
(
1995
).
4.
B. J.
Berne
and
D.
Thirumalai
,
Annu. Rev. Phys. Chem.
37
,
401
424
(
1986
).
5.
M. E.
Tuckerman
, in
Quantum Simulations of Complex Many-Body Systems: From Theory to Algorithms
, edited by
J.
Grotendorst
,
D.
Marx
, and
A.
Muramatsu
(
John von Neumann Institute for Computing
,
Jülich
,
2002
), Vol.
10
, pp.
269
298
.
6.
M.
Parrinello
and
A.
Rahman
,
J. Chem. Phys.
80
(
2
),
860
867
(
1984
).
7.
R. P.
Feynman
and
A. R.
Hibbs
,
Quantum Mechanics and Path Integrals
(
McGraw-Hill
,
New York
,
1965
).
8.
M. F.
Herman
,
E. J.
Bruskin
, and
B. J.
Berne
,
J. Chem. Phys.
76
(
10
),
5150
5155
(
1982
).
9.
J.
Cao
and
B. J.
Berne
,
J. Chem. Phys.
99
(
4
),
2902
2916
(
1993
).
10.
M. E.
Tuckerman
,
B. J.
Berne
,
G. J.
Martyna
, and
M. L.
Klein
,
J. Chem. Phys.
99
(
4
),
2796
2808
(
1993
).
11.
M. E.
Tuckerman
,
D.
Marx
,
M. L.
Klein
, and
M.
Parrinello
,
J. Chem. Phys.
104
(
14
),
5579
5588
(
1996
).
12.
R.
Hall
and
B. J.
Berne
,
J. Chem. Phys.
81
(
8
),
3641
3643
(
1984
).
13.
M. J.
Gillan
,
Phys. Rev. Lett.
58
(
6
),
563
566
(
1987
).
14.
K.
Singer
and
W.
Smith
,
Mol. Phys.
64
(
6
),
1215
1231
(
1988
).
15.
M. H.
Müser
,
Comput. Phys. Commun.
147
(
1-2
),
83
86
(
2002
).
16.
A. N.
Drozdov
and
P.
Talkner
,
J. Chem. Phys.
109
(
6
),
2080
2091
(
1998
).
17.
M.
Ceriotti
,
M.
Parrinello
,
T. E.
Markland
, and
D. E.
Manolopoulos
,
J. Chem. Phys.
133
(
12
),
124104
(
2010
).
18.
See supplementary material at http://dx.doi.org/10.1063/1.4954990 for the case of normal mode PIMD.
19.
E. L.
Pollock
and
D. M.
Ceperley
,
Phys. Rev. B
30
(
5
),
2555
2568
(
1984
).
20.
J.
Liu
,
J. Chem. Phys.
140
(
22
),
224107
(
2014
).
21.
J.
Liu
,
D.
Li
, and
X.
Liu
,
Sci. Sin. Chim.
46
(
1
),
27
37
(
2016
).
22.
M.
Tuckerman
,
B. J.
Berne
, and
G. J.
Martyna
,
J. Chem. Phys.
97
(
3
),
1990
2001
(
1992
).
23.
D.
Janežič
,
M.
Praprotnik
, and
F.
Merzel
,
J. Chem. Phys.
122
(
17
),
174101
(
2005
).
24.
A.
Ricci
and
G.
Ciccotti
,
Mol. Phys.
101
(
12
),
1927
1931
(
2003
).
25.
G.
Bussi
and
M.
Parrinello
,
Phys. Rev. E
75
(
5
),
056707
(
2007
).
26.
B.
Leimkuhler
and
C.
Matthews
,
J. Chem. Phys.
138
(
17
),
174102
(
2013
).
27.
B.
Leimkuhler
and
C.
Matthews
,
Appl. Math. Res. Express
2013
(
1
),
34
56
(
2013
).
28.
B.
Leimkuhler
,
C.
Matthews
, and
G.
Stoltz
,
IMA J. Numer. Anal.
36
(
1
),
13
79
(
2016
).
29.
J.
Liu
and
Z.
Zhang
,
J. Chem. Phys.
144
(
3
),
034307
(
2016
).
30.
H.
Partridge
and
D. W.
Schwenke
,
J. Chem. Phys.
106
(
11
),
4618
4639
(
1997
).
31.
I. F.
Silvera
and
V. V.
Goldman
,
J. Chem. Phys.
69
(
9
),
4209
4213
(
1978
).
32.
G. S.
Fanourgakis
and
S. S.
Xantheas
,
J. Chem. Phys.
128
(
7
),
074506
(
2008
).
33.
J.
Liu
, “
An efficient algorithm for path integral Liouville dynamics
” (unpublished).
34.
M.
Rossi
,
M.
Ceriotti
, and
D. E.
Manolopoulos
,
J. Chem. Phys.
140
(
23
),
234116
(
2014
).
35.
J.
Cao
and
G. A.
Voth
,
J. Chem. Phys.
99
(
12
),
10070
10073
(
1993
).
36.
M.
Ceriotti
,
D. E.
Manolopoulos
, and
M.
Parrinello
,
J. Chem. Phys.
134
(
8
),
084104
(
2011
).
37.
M.
Ceriotti
and
D. E.
Manolopoulos
,
Phys. Rev. Lett.
109
(
10
),
100604
(
2012
).
38.
F.
Brieuc
,
H.
Dammak
, and
M.
Hayoun
,
J. Chem. Theory Comput.
12
(
3
),
1351
1359
(
2016
).
39.
M.
Ceriotti
and
T. E.
Markland
,
J. Chem. Phys.
138
(
1
),
014112
(
2013
).
40.
J.
Shao
, private communication (
2016
).

Supplementary Material

You do not currently have access to this content.