Coupled cluster theory is a vital cornerstone of electronic structure theory and is being applied to ever-larger systems. Stochastic approaches to quantum chemistry have grown in importance and offer compelling advantages over traditional deterministic algorithms in terms of computational demands, theoretical flexibility, or lower scaling with system size. We present a highly parallelizable algorithm of the coupled cluster Monte Carlo method involving sampling of clusters of excitors over multiple time steps. The behavior of the algorithm is investigated on the uniform electron gas and the water dimer at coupled-cluster levels including up to quadruple excitations. We also describe two improvements to the original sampling algorithm, full non-composite, and multi-spawn sampling. A stochastic approach to coupled cluster results in an efficient and scalable implementation at arbitrary truncation levels in the coupled cluster expansion.

1.
R. J.
Bartlett
and
M.
Musiał
,
Rev. Mod. Phys.
79
,
291
(
2007
).
2.
K.
Kowalski
and
P.
Piecuch
,
J. Chem. Phys.
113
,
18
(
2000
).
3.
M. E.
Harding
,
J.
Vázquez
,
B.
Ruscic
,
A. K.
Wilson
,
J.
Gauss
, and
J. F.
Stanton
,
J. Chem. Phys.
128
,
114111
(
2008
).
4.
C.
Hättig
,
W.
Klopper
,
A.
Köhn
, and
D. P.
Tew
,
Chem. Rev.
112
,
4
(
2012
).
5.
A.
Karton
,
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
6
,
292
(
2016
).
6.
D. L.
Freeman
,
Phys. Rev. B
15
,
5512
(
1977
).
7.
R. F.
Bishop
and
K. H.
Lührmann
,
Phys. Rev. B
17
,
3757
(
1978
).
8.
R. F.
Bishop
and
K. H.
Lührmann
,
Phys. Rev. B
26
,
5523
(
1982
).
9.
J. J.
Shepherd
and
A.
Grüneis
,
Phys. Rev. Lett.
110
,
226401
(
2013
); e-print arXiv:1310.6059.
10.
A.
Roggero
,
A.
Mukherjee
, and
F.
Pederiva
,
Phys. Rev. B
88
,
115138
(
2013
); e-print arXiv:1304.1549.
11.
J. S.
Spencer
and
A. J. W.
Thom
,
J. Chem. Phys.
144
,
084108
(
2016
); e-print arXiv:1511.05752.
12.
J.
McClain
,
J.
Lischner
,
T.
Watson
,
D. A.
Matthews
,
E.
Ronca
,
S. G.
Louie
,
T. C.
Berkelbach
, and
G. K.-L.
Chan
,
Phys. Rev. B
93
,
235139
(
2016
); e-print arXiv:1512.04556.
13.
J. J.
Shepherd
,
J. Chem. Phys.
145
,
031104
(
2016
); e-print arXiv:1605.05699.
14.
V. A.
Neufeld
and
A. J. W.
Thom
,
J. Chem. Phys.
147
,
194105
(
2017
); e-print arXiv:1706.09923.
15.
S.
Hirata
,
I.
Grabowski
,
M.
Tobita
, and
R. J.
Bartlett
,
Chem. Phys. Lett.
345
,
475
(
2001
).
16.
S.
Hirata
,
R.
Podeszwa
,
M.
Tobita
, and
R. J.
Bartlett
,
J. Chem. Phys.
120
,
2581
(
2004
).
17.
F. R.
Manby
,
D.
Alfè
, and
M. J.
Gillan
,
Phys. Chem. Chem. Phys.
8
,
5178
(
2006
).
18.
S. J.
Nolan
,
M. J.
Gillan
,
D.
Alfè
,
N. L.
Allan
, and
F. R.
Manby
,
Phys. Rev. B
80
,
165109
(
2009
).
19.
A.
Grüneis
,
G. H.
Booth
,
M.
Marsman
,
J.
Spencer
,
A.
Alavi
, and
G.
Kresse
,
J. Chem. Theory Comput.
7
,
2780
(
2011
).
20.
G. H.
Booth
,
A.
Grüneis
,
G.
Kresse
, and
A.
Alavi
,
Nature
493
,
365
(
2013
).
21.
A.
Grüneis
,
J. Chem. Phys.
143
,
102817
(
2015
).
22.
K.
Liao
and
A.
Grüneis
,
J. Chem. Phys.
145
,
141102
(
2016
).
23.
P.
Schwerdtfeger
,
R.
Tonner
,
G. E.
Moyano
, and
E.
Pahl
,
Angew. Chem., Int. Ed.
55
,
12200
(
2016
).
24.
J.
McClain
,
Q.
Sun
,
G. K.-L.
Chan
, and
T. C.
Berkelbach
,
J. Chem. Theory Comput.
13
,
1209
(
2017
).
25.
T.
Gruber
,
K.
Liao
,
T.
Tsatsoulis
,
F.
Hummel
, and
A.
Grüneis
,
Phys. Rev. X
8
,
021043
(
2018
).
26.
K.
Raghavachari
,
G. W.
Trucks
,
J. A.
Pople
, and
M.
Head-Gordon
,
Chem. Phys. Lett.
157
,
479
(
1989
).
27.
T. J.
Lee
and
G. E.
Scuseria
,
Quantum Mechanical Electronic Structure Calculations with Chemical Accuracy
(
Springer
,
The Netherlands, Dordrecht
,
1995
), pp.
47
108
.
28.
MRCC, a quantum chemical program suite, written by
M.
Kllay
,
Z.
Rolik
,
J.
Csontos
,
P.
Nagy
,
G.
Samu
,
D.
Mester
,
J.
Cska
,
B.
Szab
,
I.
Ladjnszki
,
L.
Szegedy
,
B.
Ladczki
,
K.
Petrov
,
M.
Farkas
,
P. D.
Mezei
, and
B.
Hégel
, see http://www.mrcc.hu.
29.
E.
Agullo
,
J.
Demmel
,
J.
Dongarra
,
B.
Hadri
,
J.
Kurzak
,
J.
Langou
,
H.
Ltaief
,
P.
Luszczek
, and
S.
Tomov
,
J. Phys.: Conf. Ser.
180
,
012037
(
2009
).
30.
E.
Solomonik
,
D.
Matthews
,
J. R.
Hammond
,
J. F.
Stanton
, and
J.
Demmel
,
J. Parallel Distrib. Comput.
74
,
3176
(
2014
).
31.
N.
Flocke
and
R. J.
Bartlett
,
J. Chem. Phys.
121
,
10935
(
2004
).
32.
M.
Ziółkowski
,
B.
Jansík
,
T.
Kjærgaard
, and
P.
Jørgensen
,
J. Chem. Phys.
133
,
014107
(
2010
).
33.
C.
Riplinger
and
F.
Neese
,
J. Chem. Phys.
138
,
034106
(
2013
).
34.
C.
Riplinger
,
B.
Sandhoefer
,
A.
Hansen
, and
F.
Neese
,
J. Chem. Phys.
139
,
134101
(
2013
).
35.
W. M. C.
Foulkes
,
L.
Mitas
,
R. J.
Needs
, and
G.
Rajagopal
,
Rev. Mod. Phys.
73
,
33
(
2001
).
36.
A. J. W.
Thom
and
A.
Alavi
,
Phys. Rev. Lett.
99
,
143001
(
2007
).
37.
S. Y.
Willow
,
K. S.
Kim
, and
S.
Hirata
,
J. Chem. Phys.
137
,
204122
(
2012
).
38.
S. Y.
Willow
,
M. R.
Hermes
,
K. S.
Kim
, and
S.
Hirata
,
J. Chem. Theory Comput.
9
,
4396
(
2013
).
39.
R.
Baer
,
D.
Neuhauser
, and
E.
Rabani
,
Phys. Rev. Lett.
111
,
106402
(
2013
).
40.

Our implementation does not require canonical orbitals, so we have the full freedom to choose any transformation of the single particle orbitals, for example the study in Ref. 50.

41.
G. H.
Booth
,
A. J. W.
Thom
, and
A.
Alavi
,
J. Chem. Phys.
131
,
054106
(
2009
).
42.
D.
Cleland
,
G. H.
Booth
,
C.
Overy
, and
A.
Alavi
,
J. Chem. Theory Comput.
8
,
4138
(
2012
).
43.
J. J.
Shepherd
,
G.
Booth
,
A.
Grüneis
, and
A.
Alavi
,
Phys. Rev. B
85
,
081103
(
2012
).
44.
C.
Overy
,
G. H.
Booth
,
N. S.
Blunt
,
J. J.
Shepherd
,
D.
Cleland
, and
A.
Alavi
,
J. Chem. Phys.
141
,
244117
(
2014
).
45.
N. S.
Blunt
,
S. D.
Smart
,
G. H.
Booth
, and
A.
Alavi
,
J. Chem. Phys.
143
,
134117
(
2015
).
46.
N. S.
Blunt
,
G. H.
Booth
, and
A.
Alavi
,
J. Chem. Phys.
146
,
244105
(
2017
).
47.
A. J. W.
Thom
,
Phys. Rev. Lett.
105
,
263004
(
2010
).
48.
R. S. T.
Franklin
,
J. S.
Spencer
,
A.
Zoccante
, and
A. J. W.
Thom
,
J. Chem. Phys.
144
,
044111
(
2016
).
49.
C. J. C.
Scott
and
A. J. W.
Thom
,
J. Chem. Phys.
147
,
124105
(
2017
); e-print arXiv:1706.07017.
50.
V. A.
Neufeld
and
A. J. W.
Thom
, “
Exciting determinants in Quantum Monte Carlo: Loading the dice with fast, low memory weights
,”
J. Chem. Theory Comput.
(published online).
51.
J. E.
Deustua
,
J.
Shen
, and
P.
Piecuch
,
Phys. Rev. Lett.
119
,
223003
(
2017
).
52.
G. H.
Booth
,
S. D.
Smart
, and
A.
Alavi
,
Mol. Phys.
112
,
1855
(
2014
).
53.
J. S.
Spencer
,
N. S.
Blunt
, and
W. M.
Foulkes
,
J. Chem. Phys.
136
,
054110
(
2012
).
54.
F. R.
Petruzielo
,
A. A.
Holmes
,
H. J.
Changlani
,
M. P.
Nightingale
, and
C. J.
Umrigar
,
Phys. Rev. Lett.
109
,
230201
(
2012
).
55.

We set psize(s)=1i=0s1psize(i) in order to ensure a normalised probability.

56.
C. J.
Umrigar
,
M. P.
Nightingale
, and
K. J.
Runge
,
J. Chem. Phys.
99
,
2865
(
1993
).
57.
J. M.
Turney
,
A. C.
Simmonett
,
R. M.
Parrish
,
E. G.
Hohenstein
,
F.
Evangelista
,
J. T.
Fermann
,
B. J.
Mintz
,
L. A.
Burns
,
J. J.
Wilke
,
M. L.
Abrams
,
N. J.
Russ
,
M. L.
Leininger
,
C. L.
Janssen
,
E. T.
Seidl
,
W. D.
Allen
,
H. F.
Schaefer
,
R. A.
King
,
E. F.
Valeev
,
C. D.
Sherrill
, and
T. D.
Crawford
,
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
2
,
556
(
2012
).
58.
Q.
Sun
,
T. C.
Berkelbach
,
N. S.
Blunt
,
G. H.
Booth
,
S.
Guo
,
Z.
Li
,
J.
Liu
,
J. D.
McClain
,
E. R.
Sayfutyarova
,
S.
Sharma
,
S.
Wouters
, and
G. K.-L.
Chan
,
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
8
,
e1340
(
2018
).
59.
J. M.
Foster
and
S. F.
Boys
,
Rev. Mod. Phys.
32
,
300
(
1960
).
60.
H.
Flyvbjerg
and
H. G.
Petersen
,
J. Chem. Phys.
91
,
461
(
1989
).
61.
W. A.
Vigor
,
J. S.
Spencer
,
M. J.
Bearpark
, and
A. J. W.
Thom
,
J. Chem. Phys.
142
,
104101
(
2015
); e-print arXiv:1407.1753.
62.
T. E.
Oliphant
,
Guide to NumPy
, 2nd ed. (
CreateSpace Independent Publishing Platform
,
USA
,
2015
).
63.
W.
McKinney
, in
Proceedings of the 9th Python in Science Conference
, edited by
S.
van der Walt
and
J.
Millman
(
SciPy
,
2010
), pp.
51
56
, available at http://conference.scipy.org/proceedings/scipy2010/mckinney.html.
65.
J. D.
Hunter
,
Comput. Sci. Eng.
9
,
90
(
2007
).
67.

Booth et al. used a custom hash function based upon the list of occupied orbitals. We find hashing the bit string representation of the determinant simpler and computationally more efficient, whilst giving at least as good distribution over processors when a hash function of sufficient quality is used, and therefore use the MurmurHash2 function.

68.

A node may consist of a single processor or multiple processors. Within the MPI paradigm, we distribute over MPI ranks, where each rank contains one or more threads.

69.
R. M.
Martin
,
Electronic Structure
(
Cambridge University Press
,
Cambridge
,
2004
), pp.
100
118
.
70.
G.
Giuliani
and
G.
Vignale
,
Quantum Theory of the Electron Liquid
(
Cambridge University Press
,
Cambridge
,
2005
), pp.
1
68
.
71.
P.-F.
Loos
and
P. M. W.
Gill
,
Wiley Interdiscip. Rev.: Comput. Mol. Sci.
6
,
410
(
2016
); e-print arXiv:1601.03544.
72.
L. K.
Wagner
and
D. M.
Ceperley
,
Rep. Prog. Phys.
79
,
094501
(
2016
); e-print arXiv:1602.01344.
73.

We note that available computational resources rarely grow polynomially with system size!.

74.
J. R.
Lane
,
J. Chem. Theory Comput.
9
,
316
(
2013
).
75.
E.
Papajak
and
D. G.
Truhlar
,
J. Chem. Theory Comput.
7
,
10
(
2011
).
76.
Z.
Rolik
,
L.
Szegedy
,
I.
Ladjnszki
,
B.
Ladczki
, and
M.
Kllay
,
J. Chem. Phys.
139
,
094105
(
2013
).
77.
A. A.
Holmes
,
H. J.
Changlani
, and
C. J.
Umrigar
,
J. Chem. Theory Comput.
12
,
1561
(
2016
).
78.
T. H.
Dunning
,
J. Chem. Phys.
90
,
1007
(
1989
).
79.

Using MRCC,28 a single iteration was not completed within a week on a 32-core node on the Cambridge Service for Data Driven Discovery (CSD3) operated by the University of Cambridge Research Computing Service (http://www.csd3.cam.ac.uk/). This provides a lower estimate of 2 node-months, or 46000 core hours.

80.

Terminology from Ref. 50 has been used. Single excitations were sampled uniformly, double excitations with the heat bath excitation generator. Idea by and mentioned in Holmes et al.77 

81.

These approaches were conceived and implemented some time before the work in Ref. 49, and briefly referred to in the same, and are reported in fuller detail now for completeness.

82.
H.
Jagode
,
A.
Danalis
, and
J.
Dongarra
,
Int. High Perform. Comput. Appl.
32
,
540
(
2018
).
You do not currently have access to this content.