第一原理トランスコリレーテッド(transcorrelated)法の開発および凝縮体へのその応用

Benjamin C. Corbett, Akimasa Miyake

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. We examine systems of size up to $12\ifmmode\times\else\texttimes\fi{}12$ lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by $2.4\ifmmode\times\else\texttimes\fi{}$--$14\ifmmode\times\else\texttimes\fi{}$, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Thimo Grotz, Libby Heaney, T Strunz Walter

啓史 松本

Alexey Uvarov, Artur F. Izmaylov

In electronic structure calculations, the transcorrelated method consists in transforming the Hamiltonian so as to remove the Coulomb cusp in its eigenfunctions. As a result, the wavefunction can be described more accurately without increasing the size of the basis set. However, the transcorrelated Hamiltonian is non-Hermitian and non-normal, which makes many common quantum algorithms inapplicable. Recently, a quantum eigenvalue estimation algorithm (QEVE) was proposed for non-Hermitian Hamiltonians with real spectra [FOCS 65, 1051 (2024)]. Although the asymptotic scaling of this algorithm with the desired accuracy is shown to be optimal, the constant factor in its complexity scaling has not been analyzed. Here we investigate the cost of QEVE applied to transcorrelated electronic Hamiltonians of second-row atoms and compare it to the cost of applying standard qubitization to non-transcorrelated Hamiltonians. We find that, with the xTC approximation, the T gate count of QEVE in the minimal STO-6G basis is between those of standard qubitization in the cc-pVTZ and cc-pVQZ bases. The accuracy of the transcorrelated energy differs between systems: for Li and Be, it is more accurate than the cc-pVQZ energy, while for larger atoms, the error gradually increases, exceeding the cc-pVDZ level for O, F, and Ne.

Alberto Baiardi, Michał Lesiuk, Markus Reiher

In this work, we present the first implementation of the transcorrelated electronic Hamiltonian in an optimization procedure for matrix product states by the density matrix renormalization group (DMRG) algorithm. In the transcorrelation ansatz, the electronic Hamiltonian is similarity-transformed with a Jastrow factor to describe the cusp in the wave function at electron–electron coalescence. As a result, the wave function is easier to approximate accurately with the conventional expansion in terms of one-particle basis functions and Slater determinants. The transcorrelated Hamiltonian in first quantization comprises up to three-body interactions, which we deal with in the standard way by applying robust density fitting to two- and three-body integrals entering the second-quantized representation of this Hamiltonian. The lack of hermiticity of the transcorrelated Hamiltonian is taken care of along the lines of the first work on transcorrelated DMRG [ J. Chem. Phys. 2020, 153, 164115] by encoding it as a matrix product operator and optimizing the corresponding ground state wave function with imaginary-time time-dependent DMRG. We demonstrate our quantum chemical transcorrelated DMRG approach at the example of several atoms and first-row diatomic molecules. We show that transcorrelation improves the convergence rate to the complete basis set limit in comparison to conventional DMRG. Moreover, we study extensions of our approach that aim at reducing the cost of handling the matrix product operator representation of the transcorrelated Hamiltonian.

Alberto Baiardi, Markus Reiher

We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG with the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.