-
Effective Hamiltonians for the study of real metals using quantum chemical theories
Authors:
Tina N. Mihm,
Tobias Schäfer,
Sai Kumar Ramadugu,
Andreas Grüneis,
James J. Shepherd
Abstract:
Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and…
▽ More
Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and doubles (CCSD) is one of the most important steps in moving towards chemical accuracy and, in recent years, its scope has successfully been expanded to the study of insulating surfaces and solids. Here, we show that CCSD theory can also be applied to real metals. In so doing, we overcome the limitation of needing extremely large supercells to capture long range electronic correlation effects. An effective Hamiltonian can be found using the transition structure factor--a map of electronic excitations from the Hartree--Fock wavefunction--which has fewer finite size effects than conventional periodic boundary conditions. This not only paves the way of applying coupled cluster methods to real metals but also reduces the computational cost by two orders of magnitude compared to previous methods. Our applications to phases of lithium and silicon show a resounding success in reaching the thermodynamic limit, taking the first step towards a truly universal quantum chemical treatment of solids.
△ Less
Submitted 10 March, 2021;
originally announced March 2021.
-
The sign problem in density matrix quantum Monte Carlo
Authors:
Hayley R. Petras,
William Z. Van Benschoten,
Sai Kumar Ramadugu,
James J. Shepherd
Abstract:
Density matrix quantum Monte Carlo (DMQMC) is a recently-developed method for stochastically sampling the $N$-particle thermal density matrix to obtain exact-on-average energies for model and \emph{ab initio} systems. We report a systematic numerical study of the sign problem in DMQMC based on simulations of atomic and molecular systems. In DMQMC, the density matrix is written in an outer product…
▽ More
Density matrix quantum Monte Carlo (DMQMC) is a recently-developed method for stochastically sampling the $N$-particle thermal density matrix to obtain exact-on-average energies for model and \emph{ab initio} systems. We report a systematic numerical study of the sign problem in DMQMC based on simulations of atomic and molecular systems. In DMQMC, the density matrix is written in an outer product basis of Slater determinants and has a size of space which is the square of the number of Slater determinants. In principle this means DMQMC needs to sample a space which scales in the system size, $N$, as $\mathcal{O}[(\exp(N))^2]$. In practice, there is a system-dependent critical walker population ($N_c$) which must be exceeded in order to remove the sign problem, and this imposes limitations by way of storage and computer time. We establish that $N_c$ for DMQMC is the square of $N_c$ for FCIQMC. By contrast, the minimum $N_c$ in the interaction picture modification of DMQMC (IP-DMQMC) only is directly proportionate to the $N_c$ for FCIQMC. We find that this comes from the asymmetric propagation of IP-DMQMC compared to the symmetric propagation of canonical DMQMC. An asymmetric mode of propagation is prohibitively expensive for DMQMC because it has a much greater stochastic error. Finally, we find that the equivalence between IP-DMQMC and FCIQMC seems to extend to the initiator approximation, which is often required to study larger basis sets and other systems. This suggests IP-DMQMC offers a way to ameliorate the cost of moving between a Slater determinant space and an outer product basis.
△ Less
Submitted 3 August, 2021; v1 submitted 29 January, 2021;
originally announced February 2021.
-
Using density matrix quantum Monte Carlo for calculating exact-on-average energies for ab-initio Hamiltonians in a finite basis set
Authors:
Hayley R. Petras,
Sai Kumar Ramadugu,
Fionn D. Malone,
James J. Shepherd
Abstract:
We here apply the recently developed initiator density matrix quantum Monte Carlo (i-DMQMC) to a wide range of chemical environments using atoms and molecules in vacuum. i-DMQMC samples the exact density matrix of a Hamiltonian at finite temperature and combines the accuracy of full configuration interaction quantum Monte Carlo (FCIQMC) - full configuration interaction (FCI) or exact energies in a…
▽ More
We here apply the recently developed initiator density matrix quantum Monte Carlo (i-DMQMC) to a wide range of chemical environments using atoms and molecules in vacuum. i-DMQMC samples the exact density matrix of a Hamiltonian at finite temperature and combines the accuracy of full configuration interaction quantum Monte Carlo (FCIQMC) - full configuration interaction (FCI) or exact energies in a finite basis set - with finite temperature. By way of exploring the applicability of i-DMQMC for molecular systems, we choose to study a recently developed test set by Rubenstein and coworkers: Be, H2O, and H10 at near-equilibrium and stretched geometries. We find that, for Be and H2O, i-DMQMC delivers energies which are sub-millihartree accuracy when compared with finite temperature FCI. For H2O and both geometries of H10 we examine the difference between FT-AFQMC and i-DMQMC which in turn is an estimate of the difference in canonical versus grand canonical energies. We close with a discussion of simulation parameters (initiator error and different basis sets) and by showing energy difference calculations in the form of specific heat capacity and ionization potential calculations.
△ Less
Submitted 2 December, 2019;
originally announced December 2019.
-
Fully quantum embedding with density functional theory for full configuration interaction quantum Monte Carlo
Authors:
Hayley R. Petras,
Daniel Graham,
Sai Kumar Ramadugu,
Jason D. Goodpaster,
James J. Shepherd
Abstract:
In common with many high-accuracy electronic structure methods, the initiator adaptation of full configuration interaction quantum Monte Carlo (i-FCIQMC) has difficulty treating realistic systems with large numbers of electrons. This barrier has prevented the application of i-FCIQMC to questions of catalysis that, even for the simplest of models, require high-accuracy modeling of several features…
▽ More
In common with many high-accuracy electronic structure methods, the initiator adaptation of full configuration interaction quantum Monte Carlo (i-FCIQMC) has difficulty treating realistic systems with large numbers of electrons. This barrier has prevented the application of i-FCIQMC to questions of catalysis that, even for the simplest of models, require high-accuracy modeling of several features of the electronic structure, such as strong and dynamic correlation, and localized vs. delocalized bonding. We here present a fully-quantum embedded version of i-FCIQMC , which we apply to calculate the bond dissociation energy of an ionic bond (LiH) and a covalent bond (HF) physisorbed to a benzene molecule. The embedding is performed using a recently-developed Huzinaga projection operator approach, which affords good synergy with i-FCIQMC by minimizing the number of orbitals in the calculation. We find that, without embedding, i-FCIQMC struggles to converge these calculations due to their substantial system sizes and a lack of error cancellation between reactants and products. With embedding, the i-FCIQMC calculation converges straightforwardly to CCSD(T) benchmarks. Our results suggest that embedded i-FCIQMC will be able treat system sizes well beyond our current reach (even though embedding introduces an error). We discuss how embedding might be improved (and thus the introduced error reduced) using i-FCIQMC energies as benchmarks.
△ Less
Submitted 14 October, 2019; v1 submitted 9 June, 2019;
originally announced June 2019.