Showing 1–50 of 104 results for author: Chan, G K

Quantum chemistry, classical heuristics, and quantum advantage

Authors: Garnet Kin-Lic Chan

Abstract: We describe the problems of quantum chemistry, the intuition behind classical heuristic methods used to solve them, a conjectured form of the classical complexity of quantum chemistry problems, and the subsequent opportunities for quantum advantage. This article is written for both quantum chemists and quantum information theorists. In particular, we attempt to summarize the domain of quantum chem… ▽ More We describe the problems of quantum chemistry, the intuition behind classical heuristic methods used to solve them, a conjectured form of the classical complexity of quantum chemistry problems, and the subsequent opportunities for quantum advantage. This article is written for both quantum chemists and quantum information theorists. In particular, we attempt to summarize the domain of quantum chemistry problems as well as the chemical intuition that is applied to solve them within concrete statements (such as a classical heuristic cost conjecture and a classification of different avenues for quantum advantage) in the hope that this may stimulate future analysis. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Introducing GPU-acceleration into the Python-based Simulations of Chemistry Framework

Authors: Rui Li, Qiming Sun, Xing Zhang, Garnet Kin-Lic Chan

Abstract: We introduce the first version of GPU4PySCF, a module that provides GPU acceleration of methods in PySCF. As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs effici… ▽ More We introduce the first version of GPU4PySCF, a module that provides GPU acceleration of methods in PySCF. As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs efficiently in the integral-direct Hartree-Fock Fock build and nuclear gradient construction. Benchmark calculations show a significant speedup of two orders of magnitude with respect to the multi-threaded CPU Hartree-Fock code of PySCF, and performance comparable to other GPU-accelerated quantum chemical packages including GAMESS and QUICK on a single NVIDIA A100 GPU. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Correlation Functions From Tensor Network Influence Functionals: The Case of the Spin-Boson Model

Authors: Haimi Nguyen, Nathan Ng, Lachlan P. Lindoy, Gunhee Park, Andrew J. Millis, Garnet Kin-Lic Chan, David R. Reichman

Abstract: We investigate the application of matrix product state (MPS) representations of the influence functionals (IF) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We… ▽ More We investigate the application of matrix product state (MPS) representations of the influence functionals (IF) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We examine three different IF approaches: one based on the Kadanoff-Baym contour targeting correlation functions at all times, one based on a complex contour targeting the correlation function at a single time, and a steady state formulation which avoids imaginary or complex times, while providing access to correlation functions at all times. We show that within the IF language, the steady state formulation provides a powerful approach to evaluate equilibrium correlation functions. △ Less

Submitted 22 June, 2024; originally announced June 2024.

Benchmarking the Exponential Ansatz for the Holstein model

Authors: Junjie Yang, Zhi-Hao Cui, Ankit Mahajan, Huanchen Zhai, David R. Reichman, Garnet Kin-Lic Chan

Abstract: Polarons are quasiparticles formed as a result of lattice distortions induced by charge carriers. The single-electron Holstein model captures the fundamentals of single polaron physics. We examine the power of the exponential ansatz for the polaron ground-state wavefunction in its coupled cluster, canonical transformation, and (canonically transformed) perturbative variants across the parameter sp… ▽ More Polarons are quasiparticles formed as a result of lattice distortions induced by charge carriers. The single-electron Holstein model captures the fundamentals of single polaron physics. We examine the power of the exponential ansatz for the polaron ground-state wavefunction in its coupled cluster, canonical transformation, and (canonically transformed) perturbative variants across the parameter space of the Holstein model. Our benchmark serves to guide future developments of polaron wavefunctions beyond the single-electron Holstein model. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 10 pages, 6 figures

Towards an exact electronic quantum many-body treatment of Kondo correlation in magnetic impurities

Authors: Tianyu Zhu, Linqing Peng, Huanchen Zhai, Zhi-Hao Cui, Garnet Kin-Lic Chan

Abstract: The Kondo effect is a prototypical quantum phenomenon arising from the interaction between localized electrons in a magnetic impurity and itinerant electrons in a metallic host. Although it has served as the testing ground for quantum many-body methods for decades, the precise description of Kondo physics with material specificity remains challenging. Here, we present a systematic ab initio approa… ▽ More The Kondo effect is a prototypical quantum phenomenon arising from the interaction between localized electrons in a magnetic impurity and itinerant electrons in a metallic host. Although it has served as the testing ground for quantum many-body methods for decades, the precise description of Kondo physics with material specificity remains challenging. Here, we present a systematic ab initio approach to converge towards an exact zero-temperature electronic treatment of Kondo correlations. Across a series of 3d transition metals, we extract Kondo temperatures matching the subtle experimental trends, with an accuracy far exceeding that of standard models. We further obtain microscopic insight into the origin of these trends. More broadly, we demonstrate the possibility to start from fully ab initio many-body simulations and push towards the realm of converged predictions. △ Less

Submitted 17 June, 2024; v1 submitted 28 May, 2024; originally announced May 2024.

Comments: 11 pages, 3 figures, with supplementary materials

Tensor Network Computations That Capture Strict Variationality, Volume Law Behavior, and the Efficient Representation of Neural Network States

Authors: Wen-Yuan Liu, Si-Jing Du, Ruojing Peng, Johnnie Gray, Garnet Kin-Lic Chan

Abstract: We introduce a change of perspective on tensor network states that is defined by the computational graph of the contraction of an amplitude. The resulting class of states, which we refer to as tensor network functions, inherit the conceptual advantages of tensor network states while removing computational restrictions arising from the need to converge approximate contractions. We use tensor networ… ▽ More We introduce a change of perspective on tensor network states that is defined by the computational graph of the contraction of an amplitude. The resulting class of states, which we refer to as tensor network functions, inherit the conceptual advantages of tensor network states while removing computational restrictions arising from the need to converge approximate contractions. We use tensor network functions to compute strict variational estimates of the energy on loopy graphs, analyze their expressive power for ground-states, show that we can capture aspects of volume law time evolution, and provide a mapping of general feed-forward neural nets onto efficient tensor network functions. Our work expands the realm of computable tensor networks to ones where accurate contraction methods are not available, and opens up new avenues to use tensor networks. △ Less

Submitted 21 May, 2024; v1 submitted 6 May, 2024; originally announced May 2024.

Comments: 4+9 pages, 13 figures

Convergence Analysis of the Stochastic Resolution of Identity: Comparing Hutchinson to Hutch++ for the Second-Order Green's Function

Authors: Leopoldo Mejía, Sandeep Sharma, Roi Baer, Garnet Kin-Lic Chan, Eran Rabani

Abstract: Stochastic orbital techniques offer reduced computational scaling and memory requirements to describe ground and excited states at the cost of introducing controlled statistical errors. Such techniques often rely on two basic operations, stochastic trace estimation and stochastic resolution of identity, both of which lead to statistical errors that scale with the number of stochastic realizations… ▽ More Stochastic orbital techniques offer reduced computational scaling and memory requirements to describe ground and excited states at the cost of introducing controlled statistical errors. Such techniques often rely on two basic operations, stochastic trace estimation and stochastic resolution of identity, both of which lead to statistical errors that scale with the number of stochastic realizations ($N_ξ$) as $\sqrt{N_ξ^{-1}}$. Reducing the statistical errors without significantly increasing $N_ξ$ has been challenging and is central to the development of efficient and accurate stochastic algorithms. In this work, we build upon recent progress made to improve stochastic trace estimation based on the ubiquitous Hutchinson's algorithm and propose a two-step approach for the stochastic resolution of identity, in the spirit of the Hutch++ method. Our approach is based on employing a randomized low-rank approximation followed by a residual calculation, resulting in statistical errors that scale much better than $\sqrt{N_ξ^{-1}}$. We implement the approach within the second-order Born approximation for the self-energy in the computation of neutral excitations and discuss three different low-rank approximations for the two-body Coulomb integrals. Tests on a series of hydrogen dimer chains with varying lengths demonstrate that the Hutch++-like approximations are computationally more efficient than both deterministic and purely stochastic (Hutchinson) approaches for low error thresholds and intermediate system sizes. Notably, for arbitrarily large systems, the Hutchinson-like approximation outperforms both deterministic and Hutch++-like methods. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Performant Automatic Differentiation of Local Coupled Cluster Theories: Response Properties and Ab Initio Molecular Dynamics

Authors: Xing Zhang, Chenghan Li, Hong-Zhou Ye, Timothy C. Berkelbach, Garnet Kin-Lic Chan

Abstract: In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNOCC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with… ▽ More In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNOCC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared (IR) spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer. △ Less

Submitted 2 June, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Image Super-resolution Inspired Electron Density Prediction

Authors: Chenghan Li, Or Sharir, Shunyue Yuan, Garnet K. Chan

Abstract: Drawing inspiration from the domain of image super-resolution, we view the electron density as a 3D grayscale image and use a convolutional residual network to transform a crude and trivially generated guess of the molecular density into an accurate ground-state quantum mechanical density. We find that this model outperforms all prior density prediction approaches. Because the input is itself a re… ▽ More Drawing inspiration from the domain of image super-resolution, we view the electron density as a 3D grayscale image and use a convolutional residual network to transform a crude and trivially generated guess of the molecular density into an accurate ground-state quantum mechanical density. We find that this model outperforms all prior density prediction approaches. Because the input is itself a real-space density, the predictions are equivariant to molecular symmetry transformations even though the model is not constructed to be. Due to its simplicity, the model is directly applicable to unseen molecular conformations and chemical elements. We show that fine-tuning on limited new data provides high accuracy even in challenging cases of exotic elements and charge states. Our work suggests new routes to learning real-space physical quantities drawing from the established ideas of image processing. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Tensor network influence functionals in the continuous-time limit: connections to quantum embedding, bath discretization, and higher-order time propagation

Authors: Gunhee Park, Nathan Ng, David R. Reichman, Garnet Kin-Lic Chan

Abstract: We describe two developments of tensor network influence functionals (in particular, influence functional matrix product states (IF-MPS)) for quantum impurity dynamics within the fermionic setting of the Anderson impurity model. The first provides the correct extension of the IF-MPS to continuous time by introducing a related mathematical object, the boundary influence functional MPS. The second c… ▽ More We describe two developments of tensor network influence functionals (in particular, influence functional matrix product states (IF-MPS)) for quantum impurity dynamics within the fermionic setting of the Anderson impurity model. The first provides the correct extension of the IF-MPS to continuous time by introducing a related mathematical object, the boundary influence functional MPS. The second connects the dynamics described by a compressed IF-MPS to that of a quantum embedding method with a time-dependent effective bath undergoing nonunitary dynamics. Using these concepts, we implement higher-order time propagators for the quench dynamics of the Anderson impurity model within the boundary IF-MPS formalism. The calculations illustrate the ability of the current formulation to efficiently remove the time step error in standard discrete-time IF-MPS implementations as well as to interface with state vector propagation techniques. They also show the advantages of IF-MPS dynamics, with its associated highly compact effective bath dynamics, over state vector propagation with a static bath discretization. △ Less

Submitted 2 July, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Journal ref: Physical Review B 110, 045104 (2024)

Fast emulation of fermionic circuits with matrix product states

Authors: Justin Provazza, Klaas Gunst, Huanchen Zhai, Garnet K. -L. Chan, Toru Shiozaki, Nicholas C. Rubin, Alec F. White

Abstract: We describe a matrix product state (MPS) extension for the Fermionic Quantum Emulator (FQE) software library. We discuss the theory behind symmetry adapted matrix product states for approximating many-body wavefunctions of spin-1/2 fermions, and we present an open-source, MPS-enabled implementation of the FQE interface (MPS-FQE). The software uses the open-source pyblock3 and block2 libraries for… ▽ More We describe a matrix product state (MPS) extension for the Fermionic Quantum Emulator (FQE) software library. We discuss the theory behind symmetry adapted matrix product states for approximating many-body wavefunctions of spin-1/2 fermions, and we present an open-source, MPS-enabled implementation of the FQE interface (MPS-FQE). The software uses the open-source pyblock3 and block2 libraries for most elementary tensor operations, and it can largely be used as a drop-in replacement for FQE that allows for more efficient, but approximate, emulation of larger fermionic circuits. Finally, we show several applications relevant to both near-term and fault-tolerant quantum algorithms where approximate emulation of larger systems is expected to be useful: characterization of state preparation strategies for quantum phase estimation, the testing of different variational quantum eigensolver Ansätze, the numerical evaluation of Trotter errors, and the simulation of general quantum dynamics problems. In all these examples, approximate emulation with MPS-FQE allows us to treat systems that are significantly larger than those accessible with a full statevector emulator. △ Less

Submitted 24 April, 2024; v1 submitted 29 December, 2023; originally announced December 2023.

Comments: 11 pages, 6 figures

AB-G$_0$W$_0$: A practical G$_0$W$_0$ method without frequency integration based on an auxiliary boson expansion

Authors: Johannes Tölle, Garnet Kin-Lic Chan

Abstract: Common G$_0$W$_0$ implementations rely on numerical or analytical frequency integration to determine the G$_0$W$_0$ self-energy, which results in a variety of practical complications. Recently, we demonstrated an exact connection between the G$_0$W$_0$ approximation and equation-of-motion (EOM) quantum chemistry approaches [J. Chem. Phys., 158, 124123 (2023)]. Based on this connection, we propose… ▽ More Common G$_0$W$_0$ implementations rely on numerical or analytical frequency integration to determine the G$_0$W$_0$ self-energy, which results in a variety of practical complications. Recently, we demonstrated an exact connection between the G$_0$W$_0$ approximation and equation-of-motion (EOM) quantum chemistry approaches [J. Chem. Phys., 158, 124123 (2023)]. Based on this connection, we propose a new method to determine G$_0$W$_0$ quasiparticle energies which completely avoids frequency integration and its associated problems. To achieve this, we make use of an auxiliary boson (AB) expansion. We name the new approach AB-G$_0$W$_0$ and demonstrate its practical applicability in a range of molecular problems. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Block2: a comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond

Authors: Huanchen Zhai, Henrik R. Larsson, Seunghoon Lee, Zhi-Hao Cui, Tianyu Zhu, Chong Sun, Linqing Peng, Ruojing Peng, Ke Liao, Johannes Tölle, Junjie Yang, Shuoxue Li, Garnet Kin-Lic Chan

Abstract: Block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matr… ▽ More Block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency, and to support integration with external numerical packages. Here we explain the design principles and currently supported features and present numerical examples in a range of applications. △ Less

Submitted 24 November, 2023; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 20 pages, 8 figures

Ab initio quantum many-body description of superconducting trends in the cuprates

Authors: Zhi-Hao Cui, Junjie Yang, Johannes Tölle, Hong-Zhou Ye, Huanchen Zhai, Raehyun Kim, Xing Zhang, Lin Lin, Timothy C. Berkelbach, Garnet Kin-Lic Chan

Abstract: Using a systematic ab initio quantum many-body approach that goes beyond low-energy models, we directly compute the superconducting pairing order of several doped cuprate materials and structures. We find that we can correctly capture two well-known trends: the pressure effect, where pairing order increases with intra-layer pressure, and the layer effect, where the pairing order varies with the nu… ▽ More Using a systematic ab initio quantum many-body approach that goes beyond low-energy models, we directly compute the superconducting pairing order of several doped cuprate materials and structures. We find that we can correctly capture two well-known trends: the pressure effect, where pairing order increases with intra-layer pressure, and the layer effect, where the pairing order varies with the number of copper-oxygen layers. From these calculations, we observe that the strength of superexchange and the covalency at optimal doping are the best descriptors of the maximal pairing order. Our microscopic analysis further identifies short-range copper spin fluctuations, together with multi-orbital charge fluctuations, as central to the pairing trends. Our work illustrates the possibility of a quantitative computational understanding of unconventional high-temperature superconducting materials. △ Less

Submitted 12 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 10 pages, 5 figures, with supplementary materials

An ab initio correction vector restricted active space approach to the L-edge XAS and 2p3d RIXS spectra of transition metal clusters

Authors: Seunghoon Lee, Huanchen Zhai, Garnet Kin-Lic Chan

Abstract: We describe an ab initio approach to simulate L-edge X-ray absorption (XAS) and 2p3d resonant inelastic X-ray scattering (RIXS) spectroscopies. We model the strongly correlated electronic structure within a restricted active space and employ a correction vector formulation instead of sum-over-states expressions for the spectra, thus eliminating the need to calculate a large number of intermediate… ▽ More We describe an ab initio approach to simulate L-edge X-ray absorption (XAS) and 2p3d resonant inelastic X-ray scattering (RIXS) spectroscopies. We model the strongly correlated electronic structure within a restricted active space and employ a correction vector formulation instead of sum-over-states expressions for the spectra, thus eliminating the need to calculate a large number of intermediate and final electronic states. We present benchmark simulations of the XAS and RIXS spectra of the iron complexes [FeCl4]^{-1/-2} and [Fe(SCH3)4]^{-1/-2} and interpret the spectra by deconvolving the correction vectors. Our approach represents a step towards simulating the X-ray spectroscopies of larger metal cluster systems that play a pivotal role in biology. △ Less

Submitted 31 May, 2023; v1 submitted 14 May, 2023; originally announced May 2023.

Multireference protonation energetics of a dimeric model of nitrogenase iron-sulfur clusters

Authors: Huanchen Zhai, Seunghoon Lee, Zhi-Hao Cui, Lili Cao, Ulf Ryde, Garnet Kin-Lic Chan

Abstract: Characterizing the electronic structure of the iron--sulfur clusters in nitrogenase is necessary to understand their role in the nitrogen fixation process. One challenging task is to determine the protonation state of the intermediates in the nitrogen fixing cycle. Here, we use a dimeric iron--sulfur model to study relative energies of protonation at C, S or Fe. Using a composite method based on c… ▽ More Characterizing the electronic structure of the iron--sulfur clusters in nitrogenase is necessary to understand their role in the nitrogen fixation process. One challenging task is to determine the protonation state of the intermediates in the nitrogen fixing cycle. Here, we use a dimeric iron--sulfur model to study relative energies of protonation at C, S or Fe. Using a composite method based on coupled cluster and density matrix renormalization group energetics, we converge the relative energies of four protonated configurations with respect to basis set and correlation level. We find that accurate relative energies require large basis sets, as well as a proper treatment of multireference and relativistic effects. We have also tested ten density functional approximations for these systems. Most of them give large errors in the relative energies. The best performing functional in this system is B3LYP, which gives mean absolute and maximum errors of only 10 and 13 kJ/mol with respect to our correlated wavefunction estimates, respectively. Our work provides benchmark results for the calibration of new approximate electronic structure methods and density functionals for these problems. △ Less

Submitted 24 November, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Comments: 13 pages, 8 figures

Hierarchical Clifford transformations to reduce entanglement in quantum chemistry wavefunctions

Authors: Ryan V. Mishmash, Tanvi P. Gujarati, Mario Motta, Huanchen Zhai, Garnet Kin-Lic Chan, Antonio Mezzacapo

Abstract: The performance of computational methods for many-body physics and chemistry is strongly dependent on the choice of basis used to cast the problem; hence, the search for better bases and similarity transformations is important for progress in the field. So far, tools from theoretical quantum information have been not thoroughly explored for this task. Here we take a step in this direction by prese… ▽ More The performance of computational methods for many-body physics and chemistry is strongly dependent on the choice of basis used to cast the problem; hence, the search for better bases and similarity transformations is important for progress in the field. So far, tools from theoretical quantum information have been not thoroughly explored for this task. Here we take a step in this direction by presenting efficiently computable Clifford similarity transformations for quantum chemistry Hamiltonians, which expose bases with reduced entanglement in the corresponding molecular ground states. These transformations are constructed via block diagonalization of a hierarchy of truncated molecular Hamiltonians, preserving the full spectrum of the original problem. We show that the bases introduced here allow for more efficient classical and quantum computation of ground state properties. First, we find a systematic reduction of bipartite entanglement in molecular ground states as compared to standard problem representations. This entanglement reduction has implications in classical numerical methods such as ones based on the density matrix renormalization group. Then, we develop variational quantum algorithms that exploit the structure in the new bases, showing again improved results when the hierarchical Clifford transformations are used. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 14 pages, 11 figures

Journal ref: J. Chem. Theory Comput. 19, 3194 (2023)

Exact relationships between the GW approximation and equation-of-motion coupled-cluster theories through the quasi-boson formalism

Authors: Johannes Tölle, Garnet Kin-Lic Chan

Abstract: We describe the relationship between the GW approximation and various equation-of-motion (EOM) coupled-cluster (CC) theories. We demonstrate the exact equivalence of the G$_0$W$_0$ approximation and the propagator theory for an electron-boson problem in a particular excitation basis. From there, we establish equivalence within the quasi-boson picture to the IP+EA-EOM unitary coupled-cluster propag… ▽ More We describe the relationship between the GW approximation and various equation-of-motion (EOM) coupled-cluster (CC) theories. We demonstrate the exact equivalence of the G$_0$W$_0$ approximation and the propagator theory for an electron-boson problem in a particular excitation basis. From there, we establish equivalence within the quasi-boson picture to the IP+EA-EOM unitary coupled-cluster propagator. We analyze the incomplete description of screening provided by the standard similarity-transformed IP+EA-EOM-CC and the recently introduced G$_0$W$_0$ Tamm-Dancoff approximation. We further consider the approximate decoupling of IP and EA sectors in EOM-CC treatments and devise the analogous particle-hole decoupling approach for the G$_0$W$_0$ approximation. Finally, we numerically demonstrate the exact relationships and magnitude of the approximations in calculations of a set of molecular ionization potentials and electron affinities. △ Less

Submitted 17 December, 2022; originally announced December 2022.

Interacting models for twisted bilayer graphene: a quantum chemistry approach

Authors: Fabian M. Faulstich, Kevin D. Stubbs, Qinyi Zhu, Tomohiro Soejima, Rohit Dilip, Huanchen Zhai, Raehyun Kim, Michael P. Zaletel, Garnet Kin-Lic Chan, Lin Lin

Abstract: The nature of correlated states in twisted bilayer graphene (TBG) at the magic angle has received intense attention in recent years. We present a numerical study of an interacting Bistritzer-MacDonald (IBM) model of TBG using a suite of methods in quantum chemistry, including Hartree-Fock, coupled cluster singles, doubles (CCSD), and perturbative triples (CCSD(T)), as well as a quantum chemistry f… ▽ More The nature of correlated states in twisted bilayer graphene (TBG) at the magic angle has received intense attention in recent years. We present a numerical study of an interacting Bistritzer-MacDonald (IBM) model of TBG using a suite of methods in quantum chemistry, including Hartree-Fock, coupled cluster singles, doubles (CCSD), and perturbative triples (CCSD(T)), as well as a quantum chemistry formulation of the density matrix renormalization group method (DMRG). Our treatment of TBG is agnostic to gauge choices, and hence we present a new gauge-invariant formulation to detect the spontaneous symmetry breaking in interacting models. To benchmark our approach, we focus on a simplified spinless, valleyless IBM model. At integer filling ($ν=0$), all numerical methods agree in terms of energy and $C_{2z} \mathcal{T}$ symmetry breaking. Additionally, as part of our benchmarking, we explore the impact of different schemes for removing ``double-counting'' in the IBM model. Our results at integer filling suggest that cross-validation of different IBM models may be needed for future studies of the TBG system. After benchmarking our approach at integer filling, we perform the first systematic study of the IBM model near integer filling (for $|ν|< 0.2$). In this regime, we find that the ground state can be in a metallic and $C_{2z} \mathcal{T}$ symmetry breaking phase. The ground state appears to have low entropy, and therefore can be relatively well approximated by a single Slater determinant. Furthermore, we observe many low entropy states with energies very close to the ground state energy in the near integer filling regime. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Multi-site Reaction Dynamics Through Multi-fragment Density Matrix Embedding

Authors: Chenghan Li, Junjie Yang, Xing Zhang, Garnet Kin-Lic Chan

Abstract: The practical description of disordered chemical reactions, where reactions involve multiple species at multiple sites, is presently challenge using correlated electronic structure methods. Here we describe the gradient theory of multi-fragment density matrix embedding theory, which potentially provides a minimal computational framework to model such processes at the correlated electron level. We… ▽ More The practical description of disordered chemical reactions, where reactions involve multiple species at multiple sites, is presently challenge using correlated electronic structure methods. Here we describe the gradient theory of multi-fragment density matrix embedding theory, which potentially provides a minimal computational framework to model such processes at the correlated electron level. We present the derivation and implementation of the gradient theory, its validation on model systems and chemical reactions using density matrix embedding, and its application to a molecular dynamics simulation of proton transport in a small water cluster, a simple example of multi-site reaction dynamics. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Is there evidence for exponential quantum advantage in quantum chemistry?

Authors: Seunghoon Lee, Joonho Lee, Huanchen Zhai, Yu Tong, Alexander M. Dalzell, Ashutosh Kumar, Phillip Helms, Johnnie Gray, Zhi-Hao Cui, Wenyuan Liu, Michael Kastoryano, Ryan Babbush, John Preskill, David R. Reichman, Earl T. Campbell, Edward F. Valeev, Lin Lin, Garnet Kin-Lic Chan

Abstract: The idea to use quantum mechanical devices to simulate other quantum systems is commonly ascribed to Feynman. Since the original suggestion, concrete proposals have appeared for simulating molecular and materials chemistry through quantum computation, as a potential ``killer application''. Indications of potential exponential quantum advantage in artificial tasks have increased interest in this ap… ▽ More The idea to use quantum mechanical devices to simulate other quantum systems is commonly ascribed to Feynman. Since the original suggestion, concrete proposals have appeared for simulating molecular and materials chemistry through quantum computation, as a potential ``killer application''. Indications of potential exponential quantum advantage in artificial tasks have increased interest in this application, thus, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation. We conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for quantum chemistry, it may be prudent to assume exponential speedups are not generically available for this problem. △ Less

Submitted 14 November, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Journal ref: Nat Commun 14, 1952 (2023)

Differentiable quantum chemistry with PySCF for molecules and materials at the mean-field level and beyond

Authors: Xing Zhang, Garnet Kin-Lic Chan

Abstract: We introduce an extension to the PySCF package which makes it automatically differentiable. The implementation strategy is discussed, and example applications are presented to demonstrate the automatic differentiation framework for quantum chemistry methodology development. These include orbital optimization, properties, excited-state energies, and derivative couplings, at the mean-field level and… ▽ More We introduce an extension to the PySCF package which makes it automatically differentiable. The implementation strategy is discussed, and example applications are presented to demonstrate the automatic differentiation framework for quantum chemistry methodology development. These include orbital optimization, properties, excited-state energies, and derivative couplings, at the mean-field level and beyond, in both molecules and solids. We also discuss some current limitations and directions for future work. △ Less

Submitted 4 August, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

A comparison between the one- and two-step spin-orbit coupling approaches based on the ab initio Density Matrix Renormalization Group

Authors: Huanchen Zhai, Garnet Kin-Lic Chan

Abstract: The efficient and reliable treatment of both spin-orbit coupling (SOC) and electron correlation is essential for understanding f-element chemistry. We analyze two approaches to the problem, the one-step approach where both effects are treated simultaneously, and the two-step state interaction approach. We report an implementation of the ab initio density matrix renormalization group (DMRG) with a… ▽ More The efficient and reliable treatment of both spin-orbit coupling (SOC) and electron correlation is essential for understanding f-element chemistry. We analyze two approaches to the problem, the one-step approach where both effects are treated simultaneously, and the two-step state interaction approach. We report an implementation of the ab initio density matrix renormalization group (DMRG) with a one-step treatment of the SOC effect which can be compared to prior two-step treatments on an equal footing. Using a dysprosium octahedral complex and bridged dimer as benchmark systems, we identify characteristics of problems where the one-step approach is beneficial for obtaining the low-energy spectrum. △ Less

Submitted 22 September, 2022; v1 submitted 6 July, 2022; originally announced July 2022.

Comments: 9 pages, 4 figures

The chromium dimer: closing a chapter of quantum chemistry

Authors: Henrik R. Larsson, Huanchen Zhai, C. J. Umrigar, Garnet Kin-Lic Chan

Abstract: The complex electronic structure and unusual potential energy curve of the chromium dimer have fascinated scientists for decades, with agreement between theory and experiment so far elusive. Here, we present a new ab initio simulation of the potential energy curve and vibrational spectrum that significantly improves on all earlier estimates. Our data support a shift in earlier experimental assignm… ▽ More The complex electronic structure and unusual potential energy curve of the chromium dimer have fascinated scientists for decades, with agreement between theory and experiment so far elusive. Here, we present a new ab initio simulation of the potential energy curve and vibrational spectrum that significantly improves on all earlier estimates. Our data support a shift in earlier experimental assignments of a cluster of vibrational frequencies by one quantum number. The new vibrational assignment yields an experimentally derived potential energy curve in quantitative agreement with theory across all bond lengths and across all measured frequencies. By solving this long-standing problem, our results raise the possibility of quantitative quantum chemical modeling of transition metal clusters with spectroscopic accuracy. △ Less

Submitted 24 August, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

Journal ref: J. Am. Chem. Soc., 2022, 144, 15932-15937

Using hyper-optimized tensor networks and first-principles electronic structure to simulate experimental properties of the giant {Mn84} torus

Authors: Dian-Teng Chen, Phillip Helms, Ashlyn R. Hale, Minseong Lee, Chenghan Li, Johnnie Gray, George Christou, Vivien S. Zapf, Garnet Kin-Lic Chan, Hai-Ping Cheng

Abstract: The single-molecule magnet {Mn84} is a challenge to theory due to its high nuclearity. Building on our prior work which characterized the structure of the spectrum of this magnet, we directly compute two experimentally accessible observables, the field-dependent magnetization up to 75 T and the temperature-dependent heat capacity, using parameter free theory. In particular, we use first principles… ▽ More The single-molecule magnet {Mn84} is a challenge to theory due to its high nuclearity. Building on our prior work which characterized the structure of the spectrum of this magnet, we directly compute two experimentally accessible observables, the field-dependent magnetization up to 75 T and the temperature-dependent heat capacity, using parameter free theory. In particular, we use first principles calculations to derive short- and long-range exchange interactions, while we compute the exact partition function of the resulting classical Potts and Ising spin models for all 84 Mn $S=2$ spins to obtain the observables. The latter computation is possible because of a simulation methodology that uses hyper-optimized tensor network contraction, borrowing from recent techniques developed to simulate quantum supremacy circuits. We also synthesize the magnet and measure its heat capacity and field-dependent magnetization. We observe good qualitative agreement between theory and experiment, identifying an unusual peak in the heat capacity in both, as well as a plateau in the magnetization. Our work also identifies some limitations of current theoretical modeling in large magnets, such as the sensitivity to small, long-range, exchange couplings. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: 7 pages, 8 figures. DTC and PH contributed equally

Entanglement in the quantum phases of an unfrustrated Rydberg atom array

Authors: Matthew J. O'Rourke, Garnet Kin-Lic Chan

Abstract: We report on the ground state phase diagram of interacting Rydberg atoms in the unfrustrated square lattice array. Using new tensor network algorithms, we scale to large systems in two dimensions while including all long-range interactions, revealing the phases in the bulk and their analogs in accessible finite arrays. We find a greatly altered phase diagram from earlier numerical and experimental… ▽ More We report on the ground state phase diagram of interacting Rydberg atoms in the unfrustrated square lattice array. Using new tensor network algorithms, we scale to large systems in two dimensions while including all long-range interactions, revealing the phases in the bulk and their analogs in accessible finite arrays. We find a greatly altered phase diagram from earlier numerical and experimental studies, and in particular, we uncover an emergent entangled quantum nematic phase that appears in the absence of frustration. Broadly our results yield a conceptual guide for future experiments, while our techniques provide a blueprint for converging numerical studies in other lattices. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: main text: 7 pages, 5 figures; supplementary info: 10 pages, 11 figures

Journal ref: Nat Commun 14, 5397 (2023)

Quantum harmonic free energies for biomolecules and nanomaterials

Authors: Alec F. White, Chenghan Li, Garnet Kin-Lic Chan

Abstract: Obtaining the free energy of large molecules from quantum mechanical energy functions is a longstanding challenge. We describe a method that allows us to estimate, at the quantum mechanical level, the harmonic contributions to the thermodynamics of molecular systems of unprecedented size, with modest cost. Using this approach, we compute the vibrational thermodynamics of a series of diamond nanocr… ▽ More Obtaining the free energy of large molecules from quantum mechanical energy functions is a longstanding challenge. We describe a method that allows us to estimate, at the quantum mechanical level, the harmonic contributions to the thermodynamics of molecular systems of unprecedented size, with modest cost. Using this approach, we compute the vibrational thermodynamics of a series of diamond nanocrystals, and show that the error per atom decreases with system size in the limit of large systems. We further show that we can obtain the vibrational contributions to the binding free energies of prototypical protein-ligand complexes where the exact computation is too expensive to be practical. Our work raises the possibility of routine quantum mechanical estimates of thermodynamic quantities in complex systems. △ Less

Submitted 23 November, 2021; originally announced November 2021.

On the pure state $v$-representability of density matrix embedding theory

Authors: Fabian M. Faulstich, Raehyun Kim, Zhi-Hao Cui, Zaiwen Wen, Garnet Kin-Lic Chan, Lin Lin

Abstract: Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition… ▽ More Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition and the vanishing low-level gap are related to the assumption that the high-level density matrix blocks are non-interacting pure-state $v$-representable (NI-PS-V), which assumes that the low-level density matrix is constructed following the Aufbau principle. In order to relax the NI-PS-V condition, we develop an augmented Lagrangian method to match the density matrix blocks without referring to the Aufbau principle. Numerical results for 2D Hubbard and hydrogen model systems indicate that in some challenging scenarios, the relaxation of the Aufbau principle directly leads to exact matching of the density matrix blocks, which also yields improved accuracy. △ Less

Submitted 13 January, 2022; v1 submitted 18 October, 2021; originally announced October 2021.

Matrix product states with large sites

Authors: Henrik R. Larsson, Huanchen Zhai, Klaas Gunst, Garnet Kin-Lic Chan

Abstract: We explore various ways to group orbitals into clusters in a matrix product state (MPS). We explain how a generic cluster MPS can often lead to an increase in computational cost and instead propose a special cluster structure, involving only the first and last orbitals/sites, with a wider scope for computational advantage. This structure is a natural formalism to describe correlated multireference… ▽ More We explore various ways to group orbitals into clusters in a matrix product state (MPS). We explain how a generic cluster MPS can often lead to an increase in computational cost and instead propose a special cluster structure, involving only the first and last orbitals/sites, with a wider scope for computational advantage. This structure is a natural formalism to describe correlated multireference (MR) theories. We demonstrate the flexibility and usefulness of this approach by implementing various uncontracted MR configuration interaction, perturbation and linearized coupled cluster theories using an MPS with large cluster sites. Applications to the nitrogen dimer, the chromium dimer, and benzene, including up to triple excitations in the external space, demonstrate the utility of an MPS with up to two large sites. We use our results to analyze the quality of different multireference approximations. △ Less

Submitted 21 January, 2022; v1 submitted 22 September, 2021; originally announced September 2021.

Journal ref: J. Chem. Theory Comput. 2022, 18, 2, 749-762

Conservation laws in coupled cluster dynamics at finite-temperature

Authors: Ruojing Peng, Alec F. White, Huanchen Zhai, Garnet Kin-Lic Chan

Abstract: We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [{\it J. Chem. Theory Comput.} \textbf{2019}, 15, 6137-6253] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest's theorem) for all one-particle properties, while remaining energy conserving for time-independent… ▽ More We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [{\it J. Chem. Theory Comput.} \textbf{2019}, 15, 6137-6253] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest's theorem) for all one-particle properties, while remaining energy conserving for time-independent Hamiltonians. We present the time-dependent orbital-optimized coupled cluster doubles method (Keldysh-OCCD) in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of non-equilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H$_2$, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model. △ Less

Submitted 4 June, 2021; originally announced June 2021.

Comments: 18 pages, 13 figures

The Fermionic Quantum Emulator

Authors: Nicholas C. Rubin, Klaas Gunst, Alec White, Leon Freitag, Kyle Throssell, Garnet Kin-Lic Chan, Ryan Babbush, Toru Shiozaki

Abstract: The fermionic quantum emulator (FQE) is a collection of protocols for emulating quantum dynamics of fermions efficiently taking advantage of common symmetries present in chemical, materials, and condensed-matter systems. The library is fully integrated with the OpenFermion software package and serves as the simulation backend. The FQE reduces memory footprint by exploiting number and spin symmetry… ▽ More The fermionic quantum emulator (FQE) is a collection of protocols for emulating quantum dynamics of fermions efficiently taking advantage of common symmetries present in chemical, materials, and condensed-matter systems. The library is fully integrated with the OpenFermion software package and serves as the simulation backend. The FQE reduces memory footprint by exploiting number and spin symmetry along with custom evolution routines for sparse and dense Hamiltonians, allowing us to study significantly larger quantum circuits at modest computational cost when compared against qubit state vector simulators. This release paper outlines the technical details of the simulation methods and key advantages. △ Less

Submitted 19 October, 2021; v1 submitted 28 April, 2021; originally announced April 2021.

Journal ref: Quantum 5, 568 (2021)

Low communication high performance ab initio density matrix renormalization group algorithms

Authors: Huanchen Zhai, Garnet Kin-Lic Chan

Abstract: There has been recent interest in the deployment of ab initio density matrix renormalization group computations on high performance computing platforms. Here, we introduce a reformulation of the conventional distributed memory ab initio DMRG algorithm that connects it to the conceptually simpler and advantageous sum of sub-Hamiltonians approach. Starting from this framework, we further explore a h… ▽ More There has been recent interest in the deployment of ab initio density matrix renormalization group computations on high performance computing platforms. Here, we introduce a reformulation of the conventional distributed memory ab initio DMRG algorithm that connects it to the conceptually simpler and advantageous sum of sub-Hamiltonians approach. Starting from this framework, we further explore a hierarchy of parallelism strategies, that includes (i) parallelism over the sum of sub-Hamiltonians, (ii) parallelism over sites, (iii) parallelism over normal and complementary operators, (iv) parallelism over symmetry sectors, and (v) parallelism over dense matrix multiplications. We describe how to reduce processor load imbalance and the communication cost of the algorithm to achieve higher efficiencies. We illustrate the performance of our new open-source implementation on a recent benchmark ground-state calculation of benzene in an orbital space of 108 orbitals and 30 electrons, with a bond dimension of up to 6000, and a model of the FeMo cofactor with 76 orbitals and 113 electrons. The observed parallel scaling from 448 to 2800 CPU cores is nearly ideal. △ Less

Submitted 22 June, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 14 pages, 6 figures; typos corrected

Journal ref: The Journal of Chemical Physics, 154, 224116 (2021)

Externally corrected CCSD with renormalized perturbative triples (R-ecCCSD(T)) and density matrix renormalization group and selected configuration interaction external sources

Authors: Seunghoon Lee, Huanchen Zhai, Sandeep Sharma, Cyrus J. Umrigar, Garnet Kin-Lic Chan

Abstract: We investigate the renormalized perturbative triples correction together with the externally corrected coupled-cluster singles and doubles (ecCCSD) method. We take the density matrix renormalization group (DMRG) and heatbath CI (HCI) as external sources for the ecCCSD equations. The accuracy is assessed for the potential energy surfaces of H2O, N2, and F2. We find that the triples correction signi… ▽ More We investigate the renormalized perturbative triples correction together with the externally corrected coupled-cluster singles and doubles (ecCCSD) method. We take the density matrix renormalization group (DMRG) and heatbath CI (HCI) as external sources for the ecCCSD equations. The accuracy is assessed for the potential energy surfaces of H2O, N2, and F2. We find that the triples correction significantly improves on ecCCSD and we do not see any instability of the renormalized triples with respect to dissociation. We explore how to balance the cost of computing the external source amplitudes with respect to the accuracy of the subsequent CC calculation. In this context, we find that very approximate wavefunctions (and their large amplitudes) serve as an efficient and accurate external source. Finally, we characterize the domain of correlation treatable using the externally corrected method and renormalized triples combination studied in this work via a well-known wavefunction diagnostic. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 20 pages, 4 figures

Constructing Tensor Network Influence Functionals for General Quantum Dynamics

Authors: Erika Ye, Garnet Kin-Lic Chan

Abstract: We describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of the bond dimensions and time-like entanglement in the tensor network description. We s… ▽ More We describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of the bond dimensions and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation. △ Less

Submitted 25 July, 2021; v1 submitted 14 January, 2021; originally announced January 2021.

Comments: 12 pages, 10 figures

Journal ref: J. Chem. Phys. 155, 044104 (2021)

A coupled cluster framework for electrons and phonons

Authors: Alec F. White, Yang Gao, Austin J. Minnich, Garnet Kin-Lic Chan

Abstract: We describe a coupled cluster framework for coupled systems of electrons and phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard-Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations which generally perform well for weak to moderate coupling. Finally, we report progres… ▽ More We describe a coupled cluster framework for coupled systems of electrons and phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard-Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations which generally perform well for weak to moderate coupling. Finally, we report progress towards an implementation for {\it ab initio} calculations on solids, and present some preliminary results on finite-size models of diamond. We also report the implementation of electron-phonon coupling matrix elements from crystalline Gaussian type orbitals (cGTO) within the PySCF program package. △ Less

Submitted 13 October, 2020; v1 submitted 28 September, 2020; originally announced September 2020.

The Ground State Electronic Energy of Benzene

Authors: Janus J. Eriksen, Tyler A. Anderson, J. Emiliano Deustua, Khaldoon Ghanem, Diptarka Hait, Mark R. Hoffmann, Seunghoon Lee, Daniel S. Levine, Ilias Magoulas, Jun Shen, Norman M. Tubman, K. Birgitta Whaley, Enhua Xu, Yuan Yao, Ning Zhang, Ali Alavi, Garnet Kin-Lic Chan, Martin Head-Gordon, Wenjian Liu, Piotr Piecuch, Sandeep Sharma, Seiichiro L. Ten-no, C. J. Umrigar, Jürgen Gauss

Abstract: We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$ζ$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire of… ▽ More We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$ζ$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around $-863$ m$E_{\text{H}}$. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 m$E_{\text{H}}$), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward. △ Less

Submitted 7 October, 2020; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 29 pages, 1 figure, 2 tables. SI as an ancillary file

Journal ref: J. Phys. Chem. Lett., 11, 8922 (2020)

Automatic transformation of irreducible representations for efficient contraction of tensors with cyclic group symmetry

Authors: Yang Gao, Phillip Helms, Garnet Kin-Lic Chan, Edgar Solomonik

Abstract: Tensor contractions are ubiquitous in computational chemistry and physics, where tensors generally represent states or operators and contractions express the algebra of these quantities. In this context, the states and operators often preserve physical conservation laws, which are manifested as group symmetries in the tensors. These group symmetries imply that each tensor has block sparsity and ca… ▽ More Tensor contractions are ubiquitous in computational chemistry and physics, where tensors generally represent states or operators and contractions express the algebra of these quantities. In this context, the states and operators often preserve physical conservation laws, which are manifested as group symmetries in the tensors. These group symmetries imply that each tensor has block sparsity and can be stored in a reduced form. For nontrivial contractions, the memory footprint and cost are lowered, respectively, by a linear and a quadratic factor in the number of symmetry sectors. State-of-the-art tensor contraction software libraries exploit this opportunity by iterating over blocks or using general block-sparse tensor representations. Both approaches entail overhead in performance and code complexity. With intuition aided by tensor diagrams, we present a technique, irreducible representation alignment, which enables efficient handling of Abelian group symmetries via only dense tensors, by using contraction-specific reduced forms. This technique yields a general algorithm for arbitrary group symmetric contractions, which we implement in Python and apply to a variety of representative contractions from quantum chemistry and tensor network methods. As a consequence of relying on only dense tensor contractions, we can easily make use of efficient batched matrix multiplication via Intel's MKL and distributed tensor contraction via the Cyclops library, achieving good efficiency and parallel scalability on up to 4096 Knights Landing cores of a supercomputer. △ Less

Submitted 25 September, 2022; v1 submitted 15 July, 2020; originally announced July 2020.

All-electron Gaussian-based $G_0W_0$ for Valence and Core Excitation Energies of Periodic Systems

Authors: Tianyu Zhu, Garnet Kin-Lic Chan

Abstract: We describe an all-electron $G_0W_0$ implementation for periodic systems with $k$-point sampling implemented in a crystalline Gaussian basis. Our full-frequency $G_0W_0$ method relies on efficient Gaussian density fitting integrals and includes both analytic continuation and contour deformation schemes. Due to the compactness of Gaussian bases, no virtual state truncation is required as is seen in… ▽ More We describe an all-electron $G_0W_0$ implementation for periodic systems with $k$-point sampling implemented in a crystalline Gaussian basis. Our full-frequency $G_0W_0$ method relies on efficient Gaussian density fitting integrals and includes both analytic continuation and contour deformation schemes. Due to the compactness of Gaussian bases, no virtual state truncation is required as is seen in many plane-wave formulations. Finite size corrections are included by taking the $q \to 0$ limit of the Coulomb divergence. Using our implementation, we study quasiparticle energies and band structures across a range of systems including molecules, semiconductors, rare gas solids, and metals. We find that the $G_0W_0$ band gaps of traditional semiconductors converge rapidly with respect to the basis size, even for the conventionally challenging case of ZnO. Using correlation-consistent bases of polarized triple-zeta quality, we find the mean absolute relative error of the extrapolated $G_0W_0$@PBE band gaps to be only 5.2% when compared to experimental values. For core excitation binding energies (CEBEs), we find that $G_0W_0$ predictions improve significantly over those from DFT if the $G_0W_0$ calculations are started from hybrid functionals with a high percentage of exact exchange. △ Less

Submitted 22 December, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 14 pages, 5 figures, 6 tables

Journal ref: J. Chem. Theory Comput. 17, 727-741 (2021)

Minimal matrix product states and generalizations of mean-field and geminal wavefunctions

Authors: Henrik R. Larsson, Carlos A. Jiménez-Hoyos, Garnet Kin-Lic Chan

Abstract: Simple wavefunctions of low computational cost but which can achieve qualitative accuracy across the whole potential energy surface (PES) are of relevance to many areas of electronic structure as well as to applications to dynamics. Here, we explore a class of simple wavefunctions, the minimal matrix product state (MMPS), that generalizes many simple wavefunctions in common use, such as projected… ▽ More Simple wavefunctions of low computational cost but which can achieve qualitative accuracy across the whole potential energy surface (PES) are of relevance to many areas of electronic structure as well as to applications to dynamics. Here, we explore a class of simple wavefunctions, the minimal matrix product state (MMPS), that generalizes many simple wavefunctions in common use, such as projected mean-field wavefunctions, geminal wavefunctions, and generalized valence bond states. By examining the performance of MMPSs for PESs of some prototypical systems, we find that they yield good qualitative behavior across the whole PES, often significantly improving on the aforementioned ansätze. △ Less

Submitted 26 June, 2020; v1 submitted 7 May, 2020; originally announced May 2020.

Journal ref: J. Chem. Theory Comput. 2020, 16, 5057-5066

Novel algorithms and high-performance cloud computing enable efficient fully quantum mechanical protein-ligand scoring

Authors: Narbe Mardirossian, Yuhang Wang, David A. Pearlman, Garnet Kin-Lic Chan, Toru Shiozaki

Abstract: Ranking the binding of small molecules to protein receptors through physics-based computation remains challenging. Though inroads have been made using free energy methods, these fail when the underlying classical mechanical force fields are insufficient. In principle, a more accurate approach is provided by quantum mechanical density functional theory (DFT) scoring, but even with approximations, t… ▽ More Ranking the binding of small molecules to protein receptors through physics-based computation remains challenging. Though inroads have been made using free energy methods, these fail when the underlying classical mechanical force fields are insufficient. In principle, a more accurate approach is provided by quantum mechanical density functional theory (DFT) scoring, but even with approximations, this has yet to become practical on drug discovery-relevant timescales and resources. Here, we describe how to overcome this barrier using algorithms for DFT calculations that scale on widely available cloud architectures, enabling full density functional theory, without approximations, to be applied to protein-ligand complexes with approximately 2500 atoms in tens of minutes. Applying this to a realistic example of 22 ligands binding to MCL1 reveals that density functional scoring outperforms classical free energy perturbation theory for this system. This raises the possibility of broadly applying fully quantum mechanical scoring to real-world drug discovery pipelines. △ Less

Submitted 18 April, 2020; originally announced April 2020.

Comments: 15 pages, 5 figures, 1 table

Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems

Authors: Alec F. White, Garnet Kin-Lic Chan

Abstract: We discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for… ▽ More We discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for applying the method in practice. The method is then applied to the 1D Hubbard model, the uniform electron gas at warm, dense conditions, and some simple materials. The performance on model systems at high temperatures is encouraging: for the 1-dimensional Hubbard model FT-CCSD provides a qualitatively accurate description of finite-temperature correlation effects even at $U = 8$, and it allows for the computation of systematically improvable exchange-correlation energies of the warm, dense UEG over a wide range of conditions. We highlight the obstacles that remain in using the method for realistic {\it ab initio} calculations on materials. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Ab Initio Full Cell GW+DMFT for Correlated Materials

Authors: Tianyu Zhu, Garnet Kin-Lic Chan

Abstract: Quantitative prediction of electronic properties in correlated materials requires simulations without empirical truncations and parameters. We present a method to achieve this goal through a new ab initio formulation of dynamical mean-field theory (DMFT). Instead of using small impurities defined in a low-energy subspace, which require complicated downfolded interactions which are often approximat… ▽ More Quantitative prediction of electronic properties in correlated materials requires simulations without empirical truncations and parameters. We present a method to achieve this goal through a new ab initio formulation of dynamical mean-field theory (DMFT). Instead of using small impurities defined in a low-energy subspace, which require complicated downfolded interactions which are often approximated, we describe a full cell $GW$+DMFT approach, where the impurities comprise all atoms in a unit cell or supercell of the crystal. Our formulation results in large impurity problems, which we treat here with efficient quantum chemistry impurity solvers that work on the real-frequency axis, combined with a one-shot $G_0W_0$ treatment of long-range interactions. We apply our full cell approach to bulk Si, two antiferromagnetic correlated insulators NiO and $α$-Fe$_2$O$_3$, and the paramagnetic correlated metal SrMoO$_3$, with impurities containing up to 10 atoms and 124 orbitals. We find that spectral properties, magnetic moments, and two-particle spin correlation functions are obtained in good agreement with experiments. In addition, in the metal oxide insulators, the balanced treatment of correlations involving all orbitals in the cell leads to new insights into the orbital character around the insulating gap. △ Less

Submitted 15 November, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 10 pages, 7 figures

Journal ref: Phys. Rev. X 11, 021006 (2021)

Recent developments in the PySCF program package

Authors: Qiming Sun, Xing Zhang, Samragni Banerjee, Peng Bao, Marc Barbry, Nick S. Blunt, Nikolay A. Bogdanov, George H. Booth, Jia Chen, Zhi-Hao Cui, Janus Juul Eriksen, Yang Gao, Sheng Guo, Jan Hermann, Matthew R. Hermes, Kevin Koh, Peter Koval, Susi Lehtola, Zhendong Li, Junzi Liu, Narbe Mardirossian, James D. McClain, Mario Motta, Bastien Mussard, Hung Q. Pham , et al. (24 additional authors not shown)

Abstract: PYSCF is a Python-based general-purpose electronic structure platform that both supports first-principles simulations of molecules and solids, as well as accelerates the development of new methodology and complex computational workflows. The present paper explains the design and philosophy behind PYSCF that enables it to meet these twin objectives. With several case studies, we show how users can… ▽ More PYSCF is a Python-based general-purpose electronic structure platform that both supports first-principles simulations of molecules and solids, as well as accelerates the development of new methodology and complex computational workflows. The present paper explains the design and philosophy behind PYSCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PYSCF as a development environment. We then summarize the capabilities of PYSCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PYSCF across the domains of quantum chemistry, materials science, machine learning and quantum information science. △ Less

Submitted 10 July, 2020; v1 submitted 27 February, 2020; originally announced February 2020.

Journal ref: J. Chem. Phys. 153, 024109 (2020)

Ground-state phase diagram of the three-band Hubbard model from density matrix embedding theory

Authors: Zhi-Hao Cui, Chong Sun, Ushnish Ray, Bo-Xiao Zheng, Qiming Sun, Garnet Kin-Lic Chan

Abstract: We determine the ground-state phase diagram of the three-band Hubbard model across a range of model parameters using density matrix embedding theory. We study the atomic-scale nature of the antiferromagnetic (AFM) and superconducting (SC) orders, explicitly including the oxygen degrees of freedom. All parametrizations of the model display AFM and SC phases, but the decay of AFM order with doping i… ▽ More We determine the ground-state phase diagram of the three-band Hubbard model across a range of model parameters using density matrix embedding theory. We study the atomic-scale nature of the antiferromagnetic (AFM) and superconducting (SC) orders, explicitly including the oxygen degrees of freedom. All parametrizations of the model display AFM and SC phases, but the decay of AFM order with doping is too slow compared to the experimental phase diagram, and further, coexistence of AFM and SC orders occurs in all parameter sets. The local magnetic moment localizes entirely at the copper sites. The magnetic phase diagram is particularly sensitive to $Δ_{pd}$ and $t_{pp}$, and existing estimates of the charge transfer gap $Δ_{pd}$ appear too large in so-called minimal model parametrizations. The electron-doped side of the phase diagram is qualitatively distinct from hole-doped side and we find an unusual two-peak structure in the SC in the full model parametrization. Examining the SC order at the atomic scale, within the larger scale $d_{x^2 - y^2}$-wave SC pairing order between Cu-Cu and O-O, we also observe a local $p_{x (y)}$ [or $d_{xz (yz)}$]-symmetry modulation of the pair density on the Cu-O bonds. Our work highlights some of the features that arise in a three-band versus one-band picture, the role of the oxygen degrees of freedom in new kinds of atomic-scale SC orders, and the necessity of re-evaluating current parametrizations of the three-band Hubbard model. △ Less

Submitted 30 November, 2020; v1 submitted 14 January, 2020; originally announced January 2020.

Comments: 14 pages, 11 figures

Journal ref: Phys. Rev. Research 2, 043259 (2020)

Quantum algorithms for quantum chemistry and quantum materials science

Authors: Bela Bauer, Sergey Bravyi, Mario Motta, Garnet Kin-Lic Chan

Abstract: As we begin to reach the limits of classical computing, quantum computing has emerged as a technology that has captured the imagination of the scientific world. While for many years, the ability to execute quantum algorithms was only a theoretical possibility, recent advances in hardware mean that quantum computing devices now exist that can carry out quantum computation on a limited scale. Thus i… ▽ More As we begin to reach the limits of classical computing, quantum computing has emerged as a technology that has captured the imagination of the scientific world. While for many years, the ability to execute quantum algorithms was only a theoretical possibility, recent advances in hardware mean that quantum computing devices now exist that can carry out quantum computation on a limited scale. Thus it is now a real possibility, and of central importance at this time, to assess the potential impact of quantum computers on real problems of interest. One of the earliest and most compelling applications for quantum computers is Feynman's idea of simulating quantum systems with many degrees of freedom. Such systems are found across chemistry, physics, and materials science. The particular way in which quantum computing extends classical computing means that one cannot expect arbitrary simulations to be sped up by a quantum computer, thus one must carefully identify areas where quantum advantage may be achieved. In this review, we briefly describe central problems in chemistry and materials science, in areas of electronic structure, quantum statistical mechanics, and quantum dynamics, that are of potential interest for solution on a quantum computer. We then take a detailed snapshot of current progress in quantum algorithms for ground-state, dynamics, and thermal state simulation, and analyze their strengths and weaknesses for future developments. △ Less

Submitted 10 July, 2020; v1 submitted 10 January, 2020; originally announced January 2020.

Comments: 33 pages, 21 figures

Journal ref: Chemical Reviews, 2020

Theoretical Prediction of Magnetic Exchange Coupling Constants from Broken-Symmetry Coupled Cluster Calculations

Authors: Henry F. Schurkus, Dian-Teng Chen, Hai-Ping Cheng, Garnet Kin-Lic Chan, John F. Stanton

Abstract: Exchange coupling constants ($J$) are fundamental to the understanding of spin spectra of magnetic systems. Here we investigate the broken-symmetry (BS) approaches of Noodleman and Yamaguchi in conjunction with coupled cluster (CC) methods to obtain exchange couplings. $J$ values calculated from CC in this fashion converge smoothly towards the FCI result with increasing level of CC excitation. We… ▽ More Exchange coupling constants ($J$) are fundamental to the understanding of spin spectra of magnetic systems. Here we investigate the broken-symmetry (BS) approaches of Noodleman and Yamaguchi in conjunction with coupled cluster (CC) methods to obtain exchange couplings. $J$ values calculated from CC in this fashion converge smoothly towards the FCI result with increasing level of CC excitation. We compare this BS-CC scheme to the complementary EOM-CC approach on a selection of bridged molecular cases and give results from a few other methodologies for context. △ Less

Submitted 9 January, 2020; originally announced January 2020.

Exploring the magnetic properties of the largest single molecule magnets

Authors: Henry F. Schurkus, Dianteng Chen, Matthew J. O'Rourke, Hai-Ping Cheng, Garnet K. -L. Chan

Abstract: The giant $\{ \mathrm{Mn}_{70} \}$ and $\{ \mathrm{Mn}_{84} \}$ wheels are the largest nuclearity single-molecule magnets synthesized to date and understanding their magnetic properties poses a challenge to theory. Starting from first principles calculations, we explore the magnetic properties and excitations in these wheels using effective spin Hamiltonians. We find that the unusual geometry of t… ▽ More The giant $\{ \mathrm{Mn}_{70} \}$ and $\{ \mathrm{Mn}_{84} \}$ wheels are the largest nuclearity single-molecule magnets synthesized to date and understanding their magnetic properties poses a challenge to theory. Starting from first principles calculations, we explore the magnetic properties and excitations in these wheels using effective spin Hamiltonians. We find that the unusual geometry of the superexchange pathways leads to weakly coupled $\{ \mathrm{Mn}_{7} \}$ subunits carrying an effective $S=2$ spin. The spectrum exhibits a hierarchy of energy scales and massive degeneracies, with the lowest energy excitations arising from Heisenberg-ring-like excitations of the $\{ \mathrm{Mn}_{7} \}$ subunits around the wheel, at energies consistent with the observed temperature dependence of the magnetic susceptibility. We further suggest an important role for weak longer-range couplings in selecting the precise spin ground-state of the $\mathrm{Mn}$ wheels out of the nearly degenerate ground-state band. △ Less

Submitted 3 January, 2020; v1 submitted 1 January, 2020; originally announced January 2020.

Comments: 7 pages, 8 figures. HFS and DTC contributed equally

Ground-state properties of the hydrogen chain: insulator-to-metal transition, dimerization, and magnetic phases

Authors: Mario Motta, Claudio Genovese, Fengjie Ma, Zhi-Hao Cui, Randy Sawaya, Garnet Kin-Lic Chan, Natalia Chepiga, Phillip Helms, Carlos Jimenez-Hoyos, Andrew J. Millis, Ushnish Ray, Enrico Ronca, Hao Shi, Sandro Sorella, Edwin M. Stoudenmire, Steven R. White, Shiwei Zhang

Abstract: Accurate and predictive computations of the quantum-mechanical behavior of many interacting electrons in realistic atomic environments are critical for the theoretical design of materials with desired properties, and require solving the grand-challenge problem of the many-electron Schrodinger equation. An infinite chain of equispaced hydrogen atoms is perhaps the simplest realistic model for a bul… ▽ More Accurate and predictive computations of the quantum-mechanical behavior of many interacting electrons in realistic atomic environments are critical for the theoretical design of materials with desired properties, and require solving the grand-challenge problem of the many-electron Schrodinger equation. An infinite chain of equispaced hydrogen atoms is perhaps the simplest realistic model for a bulk material, embodying several central themes of modern condensed matter physics and chemistry, while retaining a connection to the paradigmatic Hubbard model. Here we report a combined application of cutting-edge computational methods to determine the properties of the hydrogen chain in its quantum-mechanical ground state. Varying the separation between the nuclei leads to a rich phase diagram, including a Mott phase with quasi long-range antiferromagnetic order, electron density dimerization with power-law correlations, an insulator-to-metal transition and an intricate set of intertwined magnetic orders. △ Less

Submitted 13 July, 2020; v1 submitted 4 November, 2019; originally announced November 2019.

Comments: 9 pages, 5 figures, supplemental information included as ancillary file

Journal ref: Phys. Rev. X 10, 031058 (2020)

Exact Parameterization of Fermionic Wave Functions via Unitary Coupled Cluster Theory

Authors: Francesco A. Evangelista, Garnet Kin-Lic Chan, Gustavo E. Scuseria

Abstract: A formal analysis is conducted on the exactness of various forms of unitary coupled cluster (UCC) theory based on particle-hole excitation and de-excitation operators. Both the conventional single exponential UCC parameterization and a disentangled (factorized) version are considered. We formulate a differential cluster analysis to determine the UCC amplitudes corresponding to a general quantum st… ▽ More A formal analysis is conducted on the exactness of various forms of unitary coupled cluster (UCC) theory based on particle-hole excitation and de-excitation operators. Both the conventional single exponential UCC parameterization and a disentangled (factorized) version are considered. We formulate a differential cluster analysis to determine the UCC amplitudes corresponding to a general quantum state. The exactness of conventional UCC (ability to represent any state) is explored numerically and it is formally shown to be determined by the structure of the critical points of the UCC exponential mapping. A family of disentangled UCC wave functions are shown to exactly parameterize any state, thus showing how to construct Trotter-error-free parameterizations of UCC for applications in quantum computing. From these results, we derive an exact disentangled UCC parameterization that employs an infinite sequence of particle-hole or general one- and two-body substitution operators. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Constructing Auxiliary Dynamics for Nonequilibrium Stationary States by Variance Minimization

Authors: Ushnish Ray, Garnet Kin-Lic Chan

Abstract: We present a strategy to construct guiding distribution functions (GDFs) based on variance minimization. Auxiliary dynamics via GDFs mitigates the exponential growth of variance as a function of bias in Monte Carlo estimators of large deviation functions. The variance minimization technique exploits the exact properties of eigenstates of the tilted operator that defines the biased dynamics in the… ▽ More We present a strategy to construct guiding distribution functions (GDFs) based on variance minimization. Auxiliary dynamics via GDFs mitigates the exponential growth of variance as a function of bias in Monte Carlo estimators of large deviation functions. The variance minimization technique exploits the exact properties of eigenstates of the tilted operator that defines the biased dynamics in the nonequilibrium system. We demonstrate our techniques in two classes of problems. In the continuum, we show that GDFs can be optimized to study interacting driven diffusive systems where the efficiency is systematically improved by incorporating higher correlations into the GDF. On the lattice, we use a correlator product state ansatz to study the 1D WASEP. We show that with modest resources we can capture the features of the susceptibility in large systems that marks the phase transition from uniform transport to a traveling wave state. Our work extends the repertoire of tools available to study nonequilibrium properties in realistic systems. △ Less

Submitted 19 December, 2019; v1 submitted 25 September, 2019; originally announced September 2019.

Comments: 7 pages, 5 figures and 1 table (including appendix)

Journal ref: J. Chem. Phys. 152, 104107 (2020);