Portably parallel construction of a configuration‐interaction wave function from a matrix–product state using the Charm++ framework

The construction of configuration‐interaction (CI) expansions from a matrix product state (MPS) involves numerous matrix operations and the skillful sampling of important configurations in a large Hilbert space. In this work, we present an efficient procedure for constructing CI expansions from MPS employing the parallel object‐oriented Charm++ programming framework, upon which automatic load‐balancing and object migrating facilities can be employed. This procedure was employed in the MPS‐to‐CI utility (Moritz et al., J. Chem. Phys. 2007, 126, 224109), the sampling‐reconstructed complete active‐space algorithm (SR‐CAS, Boguslawski et al., J. Chem. Phys. 2011, 134, 224101), and the entanglement‐driven genetic algorithm (EDGA, Luo et al., J. Chem. Theory Comput. 2017, 13, 4699). It enhances productivity and allows the sampling programs to evolve to their population‐expansion versions, for example, EDGA with population expansion (PE‐EDGA). Further, examples of 1,2‐dioxetanone and firefly dioxetanone anion (FDO−) molecules demonstrated the following: (a) parallel efficiencies can be persistently improved by simply by increasing the proportions of the asynchronous executions and (b) a sampled CAS‐type CI wave function of a bi‐radical‐state FDO− molecule utilizing the full valence (30e,26o) active space can be constructed within a few hours with using thousands of cores.

[1]  Wenjian Liu,et al.  iCI: Iterative CI toward full CI. , 2016, Journal of chemical theory and computation.

[2]  Markus Reiher,et al.  Exact decoupling of the Dirac Hamiltonian. II. The generalized Douglas-Kroll-Hess transformation up to arbitrary order. , 2004, The Journal of chemical physics.

[3]  Garnet Kin-Lic Chan,et al.  Spin-adapted density matrix renormalization group algorithms for quantum chemistry. , 2012, The Journal of chemical physics.

[4]  Markus Reiher,et al.  Spin-adapted matrix product states and operators. , 2016, The Journal of chemical physics.

[5]  Laxmikant V. Kale,et al.  Charm++ and AMPI: Adaptive Runtime Strategies via Migratable Objects , 2009 .

[6]  P. Zimmerman,et al.  Evaluation of full valence correlation energies and gradients. , 2019, The Journal of chemical physics.

[7]  Jörg Kussmann,et al.  Employing OpenCL to Accelerate Ab Initio Calculations on Graphics Processing Units. , 2017, Journal of chemical theory and computation.

[8]  Daniel S. Levine,et al.  Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method. , 2018, Journal of chemical theory and computation.

[9]  G. Chan,et al.  Entangled quantum electronic wavefunctions of the Mn₄CaO₅ cluster in photosystem II. , 2013, Nature chemistry.

[10]  Takeshi Yanai,et al.  Multireference configuration interaction theory using cumulant reconstruction with internal contraction of density matrix renormalization group wave function. , 2013, The Journal of chemical physics.

[11]  Laxmikant V. Kalé,et al.  NAMD: a Parallel, Object-Oriented Molecular Dynamics Program , 1996, Int. J. High Perform. Comput. Appl..

[12]  Markus Reiher,et al.  Construction of CASCI-type wave functions for very large active spaces. , 2011, The Journal of chemical physics.

[13]  Sebastian Wouters,et al.  Longitudinal static optical properties of hydrogen chains: finite field extrapolations of matrix product state calculations. , 2012, The Journal of chemical physics.

[14]  Dingxing Wang,et al.  Analysis of factors affecting execution performance of openMP programs , 2005 .

[15]  Sebastian Mai,et al.  OpenMolcas: From source code to insight. , 2019, Journal of chemical theory and computation.

[16]  Ali Alavi,et al.  Communications: Survival of the fittest: accelerating convergence in full configuration-interaction quantum Monte Carlo. , 2010, The Journal of chemical physics.

[17]  J. Olsen,et al.  Pushing configuration-interaction to the limit: Towards massively parallel MCSCF calculations. , 2017, The Journal of chemical physics.

[18]  Yingjin Ma Elucidating the multi-configurational character of the firefly dioxetanone anion and its prototypes in the biradical region using full valence active spaces. , 2020, Physical chemistry chemical physics : PCCP.

[19]  Celestino Angeli,et al.  Multireference Perturbation Theory with Cholesky Decomposition for the Density Matrix Renormalization Group , 2016, Journal of chemical theory and computation.

[20]  Laxmikant V. Kalé,et al.  OpenAtom: Scalable Ab-Initio Molecular Dynamics with Diverse Capabilities , 2016, ISC.

[21]  Takeshi Yanai,et al.  Multistate Complete-Active-Space Second-Order Perturbation Theory Based on Density Matrix Renormalization Group Reference States. , 2017, Journal of chemical theory and computation.

[22]  Basile F. E. Curchod,et al.  GPU-Accelerated State-Averaged Complete Active Space Self-Consistent Field Interfaced with Ab Initio Multiple Spawning Unravels the Photodynamics of Provitamin D3. , 2016, The journal of physical chemistry letters.

[23]  S. White,et al.  Measuring orbital interaction using quantum information theory , 2005, cond-mat/0508524.

[24]  Haibo Ma,et al.  Externally-Contracted Multireference Configuration Interaction Method Using a DMRG Reference Wave Function. , 2018, Journal of chemical theory and computation.

[25]  Sheng Guo,et al.  A projected approximation to strongly contracted N-electron valence perturbation theory for DMRG wavefunctions. , 2016, The Journal of chemical physics.

[26]  C J Umrigar,et al.  Heat-Bath Configuration Interaction: An Efficient Selected Configuration Interaction Algorithm Inspired by Heat-Bath Sampling. , 2016, Journal of chemical theory and computation.

[28]  Ali Alavi,et al.  Semistochastic Heat-Bath Configuration Interaction Method: Selected Configuration Interaction with Semistochastic Perturbation Theory. , 2016, Journal of chemical theory and computation.

[29]  M. Hanrath Wavefunction quality and error estimation of single- and multi-reference coupled-cluster and CI methods: The H4 model system , 2008 .

[30]  Takeshi Yanai,et al.  Computational Evidence of Inversion of (1)La and (1)Lb-Derived Excited States in Naphthalene Excimer Formation from ab Initio Multireference Theory with Large Active Space: DMRG-CASPT2 Study. , 2016, Journal of chemical theory and computation.

[31]  Alejandro Duran,et al.  The Design of OpenMP Tasks , 2009, IEEE Transactions on Parallel and Distributed Systems.

[32]  Sandeep Sharma,et al.  Full-valence density matrix renormalisation group calculations on meta-benzyne based on unrestricted natural orbitals. Revisit of seamless continuation from broken-symmetry to symmetry-adapted models for diradicals , 2017 .

[33]  Markus Reiher,et al.  New Approaches for ab initio Calculations of Molecules with Strong Electron Correlation. , 2015, Chimia.

[34]  J P Coe,et al.  Machine Learning Configuration Interaction. , 2018, Journal of chemical theory and computation.

[35]  David E. Bernholdt,et al.  Parallel computational chemistry made easier: The development of NWChem , 1995 .

[36]  Laxmikant V. Kalé,et al.  Automatic MPI to AMPI Program Transformation Using Photran , 2010, Euro-Par Workshops.

[37]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 3. Analytical Energy Gradients, Geometry Optimization, and First Principles Molecular Dynamics. , 2009, Journal of chemical theory and computation.

[38]  Qiming Sun,et al.  N-Electron Valence State Perturbation Theory Based on a Density Matrix Renormalization Group Reference Function, with Applications to the Chromium Dimer and a Trimer Model of Poly(p-Phenylenevinylene). , 2015, Journal of chemical theory and computation.

[39]  Reinhold Schneider,et al.  Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer , 2018, Journal of chemical theory and computation.

[40]  K. Morokuma,et al.  Multireference Ab Initio Density Matrix Renormalization Group (DMRG)-CASSCF and DMRG-CASPT2 Study on the Photochromic Ring Opening of Spiropyran. , 2013, Journal of chemical theory and computation.

[41]  Markus Reiher,et al.  Automated Selection of Active Orbital Spaces. , 2016, Journal of chemical theory and computation.

[42]  Markus Reiher,et al.  The generalized Douglas–Kroll transformation , 2002 .

[43]  Ali Alavi,et al.  Multireference linearized coupled cluster theory for strongly correlated systems using matrix product states. , 2015, The Journal of chemical physics.

[44]  G. Chan,et al.  Communication: An efficient stochastic algorithm for the perturbative density matrix renormalization group in large active spaces. , 2018, The Journal of chemical physics.

[45]  EigenmannRudolf,et al.  OpenMP to GPGPU , 2009 .

[46]  Matthias Troyer,et al.  An efficient matrix product operator representation of the quantum chemical Hamiltonian. , 2015, The Journal of chemical physics.

[47]  G. Chan,et al.  A Perturbative Density Matrix Renormalization Group Algorithm for Large Active Spaces. , 2018, Journal of chemical theory and computation.

[48]  Hai-bo Ma,et al.  Efficient Reconstruction of CAS-CI-Type Wave Functions for a DMRG State Using Quantum Information Theory and a Genetic Algorithm. , 2017, Journal of chemical theory and computation.

[49]  R. Cimiraglia,et al.  n-electron valence state perturbation theory: A spinless formulation and an efficient implementation of the strongly contracted and of the partially contracted variants , 2002 .

[50]  J. Ivanic Direct configuration interaction and multiconfigurational self-consistent-field method for multiple active spaces with variable occupations. I. Method , 2003 .

[51]  Quan Manh Phung,et al.  Cumulant Approximated Second-Order Perturbation Theory Based on the Density Matrix Renormalization Group for Transition Metal Complexes: A Benchmark Study. , 2016, Journal of chemical theory and computation.

[52]  Garnet Kin-Lic Chan,et al.  The ab-initio density matrix renormalization group in practice. , 2015, The Journal of chemical physics.

[53]  Frank Neese,et al.  Coupled Cluster Method with Single and Double Excitations Tailored by Matrix Product State Wave Functions. , 2016, The journal of physical chemistry letters.

[54]  David A Mazziotti,et al.  Active-space two-electron reduced-density-matrix method: complete active-space calculations without diagonalization of the N-electron Hamiltonian. , 2008, The Journal of chemical physics.

[55]  M. Reiher,et al.  Decomposition of density matrix renormalization group states into a Slater determinant basis. , 2007, The Journal of chemical physics.

[56]  Sandeep Sharma,et al.  The density matrix renormalization group in quantum chemistry. , 2011, Annual review of physical chemistry.

[57]  Mark S. Gordon,et al.  New Multithreaded Hybrid CPU/GPU Approach to Hartree-Fock. , 2012, Journal of chemical theory and computation.

[58]  Martin Head-Gordon,et al.  A deterministic alternative to the full configuration interaction quantum Monte Carlo method. , 2016, The Journal of chemical physics.

[59]  U. Schollwoeck The density-matrix renormalization group in the age of matrix product states , 2010, 1008.3477.

[60]  Garnet Kin-Lic Chan,et al.  Strongly contracted canonical transformation theory. , 2010, The Journal of chemical physics.

[61]  Garnet Kin-Lic Chan,et al.  Multireference quantum chemistry through a joint density matrix renormalization group and canonical transformation theory. , 2010, The Journal of chemical physics.

[62]  Francesco A Evangelista,et al.  Adaptive multiconfigurational wave functions. , 2014, The Journal of chemical physics.

[63]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation. , 2008, Journal of chemical theory and computation.

[64]  Markus Reiher,et al.  Exact decoupling of the Dirac Hamiltonian. I. General theory. , 2004, The Journal of chemical physics.

[65]  Christine M. Isborn,et al.  Excited-State Electronic Structure with Configuration Interaction Singles and Tamm–Dancoff Time-Dependent Density Functional Theory on Graphical Processing Units , 2011, Journal of chemical theory and computation.

[66]  Qiming Sun,et al.  Stochastic Multiconfigurational Self-Consistent Field Theory. , 2015, Journal of chemical theory and computation.

[67]  Giovanni Li Manni,et al.  The generalized active space concept in multiconfigurational self-consistent field methods. , 2011, The Journal of chemical physics.

[68]  Mark F. Adams,et al.  Charm++ and MPI: Combining the Best of Both Worlds , 2015, 2015 IEEE International Parallel and Distributed Processing Symposium.

[69]  Per-Olof Widmark,et al.  Density matrix averaged atomic natural orbital (ANO) basis sets for correlated molecular wave functions , 1995 .

[70]  David E. Bernholdt,et al.  High performance computational chemistry: An overview of NWChem a distributed parallel application , 2000 .

[71]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[72]  Wataru Mizukami,et al.  Density matrix renormalization group for ab initio calculations and associated dynamic correlation methods: A review of theory and applications , 2015 .

[73]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

[74]  Kenneth A. Hawick,et al.  Parallel graph component labelling with GPUs and CUDA , 2010, Parallel Comput..

[75]  Mickaël G. Delcey,et al.  Orbital entanglement and CASSCF analysis of the Ru–NO bond in a Ruthenium nitrosyl complex , 2015, Physical chemistry chemical physics : PCCP.

[76]  S. Knecht,et al.  Large-scale parallel configuration interaction. I. Nonrelativistic and scalar-relativistic general active space implementation with application to (Rb-Ba)+. , 2008, The Journal of chemical physics.

[77]  C J Umrigar,et al.  Semistochastic projector Monte Carlo method. , 2012, Physical review letters.

[78]  M. Reiher,et al.  Entanglement Measures for Single- and Multireference Correlation Effects. , 2012, The journal of physical chemistry letters.

[79]  Alistair P. Rendell,et al.  The restricted active space self-consistent-field method, implemented with a split graph unitary group approach , 1990 .

[80]  Mark Bull,et al.  Development of mixed mode MPI / OpenMP applications , 2001, Sci. Program..

[81]  Francesco Aquilante,et al.  SplitGAS Method for Strong Correlation and the Challenging Case of Cr2. , 2013, Journal of chemical theory and computation.

[82]  Takeshi Yanai,et al.  Fully Internally Contracted Multireference Configuration Interaction Theory Using Density Matrix Renormalization Group: A Reduced-Scaling Implementation Derived by Computer-Aided Tensor Factorization. , 2015, Journal of chemical theory and computation.

[83]  O. Legeza,et al.  Quantum data compression, quantum information generation, and the density-matrix renormalization group method , 2004, cond-mat/0401136.

[84]  M. Hanrath,et al.  Multi-reference coupled-cluster study of the potential energy surface of the hydrogen fluoride dissociation including excited states , 2009 .