Geminal-based high-accuracy quantum mechanics for few-body systems

愛理 川﨑, Airi Kawasaki

Sandro Sorella, Nicolas Devaux, Mario Dagrada +2 more

We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.

James K. Freericks, Jeffrey Cohn, Khadijeh Najafi +1 more

Andrey M. Tokmachev

Geminals are counterparts of two‐electron chemical bonds and lone pairs in the realm of wave functions. Antisymmetrized products of geminals provide a solid framework for studies of electron pairing in molecular systems. Natural advantages of geminal wave functions such as the correct description of bond breaking and formation make them a powerful model for electronic structure calculations. The concerted efforts to develop geminal‐based methods concentrate on increasing their accuracy and universality. The downside is the high computational cost limiting applications to small and toy systems. In contrast, the inherent local structure of optimal geminals opens pathways to development of methods for large systems. In particular, perspectives for linear scaling and hybrid quantum/classical approaches using geminals are discussed. The same principles are proposed for development of computational schemes for covalently bound and molecular crystals. © 2015 Wiley Periodicals, Inc.

D. J. Thouless, H. S. W. Massey, R. W. Hellwarth