Search

We describe the ground- and excited-state electronic structure of bulk MnO and NiO, two prototypical correlated electron materials, using coupled cluster theory with single and double excitations (CCSD). As a corollary, this work also reports the first implementation of unrestricted periodic ab initio equation-of motion CCSD. Starting from a Hartree-Fock reference, we find fundamental gaps of 3.46 eV and 4.83 eV for MnO and NiO respectively for the 16 unit supercell, slightly overestimated compared to experiment, although finite-size scaling suggests that the gap is more severely overestimated in the thermodynamic limit. From the character of the correlated electronic bands we find both MnO and NiO to lie in the intermediate Mott/charge-transfer insulator regime, although NiO appears as a charge transfer insulator when only the fundamental gap is considered. While the lowest quasiparticle excitations are of metal 3d and O 2p character in most of the Brillouin zone, near the {\Gamma} point, the lowest conduction band quasiparticles are of s character. Our study supports the potential of coupled cluster theory to provide high level many-body insights into correlated solids.

A large collaboration carefully benchmarks 20 first principles many-body electronic structure methods on a test set of 7 transition metal atoms, and their ions and monoxides. Good agreement is attained between the 3 systematically converged methods, resulting in experiment-free reference values. These reference values are used to assess the accuracy of modern emerging and scalable approaches to the many-electron problem. The most accurate methods obtain energies indistinguishable from experimental results, with the agreement mainly limited by the experimental uncertainties. Comparison between methods enables a unique perspective on calculations of many-body systems of electrons., Comment: Simons collaboration on the many-electron problem

The evolution of the fireball in heavy ion collisions is an isentropic process, meaning that it follows a trajectory of constant entropy per baryon in the phase diagram of the strong interaction. Responsible for the collective acceleration of the fireball is the speed of sound of the system, while fluctuations of conserved charges are encoded in quark-number susceptibilities: together, they leave their imprint in final observables. Here, this isentropic evolution will be analysed within chiral effective models that account for both chiral and center symmetry breaking, two central aspects of QCD. Our discussion focusses on the impact on the isentropic trajectories of the treatment of high-momentum modes, of the meson contribution to thermodynamics and of the number of quark flavours., Comment: 6 pages, 5 figures, contribution to the proceedings of the FAIR next generation scientists - 6th Edition Workshop

We describe how to incorporate symmetries of the Hamiltonian into auxiliary-field quantum Monte Carlo calculations (AFQMC). Focusing on the case of Abelian symmetries, we show that the computational cost of most steps of an AFQMC calculation is reduced by $N_k^{-1}$, where $N_k$ is the number of irreducible representations of the symmetry group. We apply the formalism to a molecular system as well as to several crystalline solids. In the latter case, the lattice translational group provides increasing savings as the number of k points is increased, which is important in enabling calculations that approach the thermodynamic limit. The extension to non-Abelian symmetries is briefly discussed., Comment: 13 pages, 7 figures

Recent work has dramatically reduced the gate complexity required to quantum simulate chemistry by using linear combinations of unitaries based methods to exploit structure in the plane wave basis Coulomb operator. Here, we show that one can achieve similar scaling even for arbitrary basis sets (which can be hundreds of times more compact than plane waves) by using qubitized quantum walks in a fashion that takes advantage of structure in the Coulomb operator, either by directly exploiting sparseness, or via a low rank tensor factorization. We provide circuits for several variants of our algorithm (which all improve over the scaling of prior methods) including one with $\widetilde{\cal O}(N^{3/2} \lambda)$ T complexity, where $N$ is number of orbitals and $\lambda$ is the 1-norm of the chemistry Hamiltonian. We deploy our algorithms to simulate the FeMoco molecule (relevant to Nitrogen fixation) and obtain circuits requiring about seven hundred times less surface code spacetime volume than prior quantum algorithms for this system, despite us using a larger and more accurate active space., Comment: 44 pages, 17 figures, formatted for Quantum

The accurate computation of Hamiltonian ground, excited, and thermal states on quantum computers stands to impact many problems in the physical and computer sciences, from quantum simulation to machine learning. Given the challenges posed in constructing large-scale quantum computers, these tasks should be carried out in a resource-efficient way. In this regard, existing techniques based on phase estimation or variational algorithms display potential disadvantages; phase estimation requires deep circuits with ancillae, that are hard to execute reliably without error correction, while variational algorithms, while flexible with respect to circuit depth, entail additional high-dimensional classical optimization. Here, we introduce the quantum imaginary time evolution and quantum Lanczos algorithms, which are analogues of classical algorithms for finding ground and excited states. Compared to their classical counterparts, they require exponentially less space and time per iteration, and can be implemented without deep circuits and ancillae, or high-dimensional optimization. We furthermore discuss quantum imaginary time evolution as a subroutine to generate Gibbs averages through an analog of minimally entangled typical thermal states. Finally, we demonstrate the potential of these algorithms via an implementation using exact classical emulation as well as through prototype circuits on the Rigetti quantum virtual machine and Aspen-1 quantum processing unit., Comment: 18 pages, 7 figures; improved figures and text

We describe an algorithm to reduce the cost of auxiliary-field quantum Monte Carlo (AFQMC) calculations for the electronic structure problem. The technique uses a nested low-rank factorization of the electron repulsion integral (ERI). While the cost of conventional AFQMC calculations in Gaussian bases scales as $\mathcal{O}(N^4)$ where $N$ is the size of the basis, we show that ground-state energies can be computed through tensor decomposition with reduced memory requirements and sub-quartic scaling. The algorithm is applied to hydrogen chains and square grids, water clusters, and hexagonal BN. In all cases we observe significant memory savings and, for larger systems, reduced, sub-quartic simulation time., Comment: 14 pages, 13 figures, expanded dataset and text

The quantum simulation of quantum chemistry is a promising application of quantum computers. However, for N molecular orbitals, the $\mathcal{O}(N^4)$ gate complexity of performing Hamiltonian and unitary Coupled Cluster Trotter steps makes simulation based on such primitives challenging. We substantially reduce the gate complexity of such primitives through a two-step low-rank factorization of the Hamiltonian and cluster operator, accompanied by truncation of small terms. Using truncations that incur errors below chemical accuracy, we are able to perform Trotter steps of the arbitrary basis electronic structure Hamiltonian with $\mathcal{O}(N^3)$ gate complexity in small simulations, which reduces to $\mathcal{O}(N^2 \log N)$ gate complexity in the asymptotic regime, while our unitary Coupled Cluster Trotter step has $\mathcal{O}(N^3)$ gate complexity as a function of increasing basis size for a given molecule. In the case of the Hamiltonian Trotter step, these circuits have $\mathcal{O}(N^2)$ depth on a linearly connected array, an improvement over the $\mathcal{O}(N^3)$ scaling assuming no truncation. As a practical example, we show that a chemically accurate Hamiltonian Trotter step for a 50 qubit molecular simulation can be carried out in the molecular orbital basis with as few as 4,000 layers of parallel nearest-neighbor two-qubit gates, consisting of fewer than 100,000 non-Clifford rotations. We also apply our algorithm to iron-sulfur clusters relevant for elucidating the mode of action of metalloenzymes., Comment: 8 pages, 4 figures

In this paper we provide a consistent framework to address the notorious decomposition of the single-photon total angular momentum (TAM) into a spin (SAM) and an orbital (OAM) component. In particular, we find that the canonical SAM and OAM, generators of internal and spatial rotations respectively in the space of physical states, are mutually compatible but unsharp quantum observables, therefore POVM (Positive Operator-Valued Measures) describe their joint measurements. On the other hand, a non-canonical SAM and OAM can be constructed that is mutually incompatible but represent sharp quantum observables, thus PVM (Projector-Valued Measurements), reflecting their consistency with the transversality condition characterizing singlephoton wavefunctions. We discuss the implementations on joint measurements for both decompositions and provide an explicit calculation of all these quantities for circularly polarized Gaussian single-photon states., Comment: 10 + 3 pages, 1 figure

We propose an algorithm for accurate, systematic and scalable computation of interatomic forces within the auxiliary-field Quantum Monte Carlo (AFQMC) method. The algorithm relies on the Hellman-Fenyman theorem, and incorporates Pulay corrections in the presence of atomic orbital basis sets. We benchmark the method for small molecules by comparing the computed forces with the derivatives of the AFQMC potential energy surface, and by direct comparison with other quantum chemistry methods. We then perform geometry optimizations using the steepest descent algorithm in larger molecules. With realistic basis sets, we obtain equilibrium geometries in agreement, within statistical error bars, with experimental values. The increase in computational cost for computing forces in this approach is only a small prefactor over that of calculating the total energy. This paves the way for a general and efficient approach for geometry optimization and molecular dynamics within AFQMC., Comment: 5 pages, 4 figures

The auxiliary-field quantum Monte Carlo (AFQMC) method provides a computational framework for solving the time-independent Schroedinger equation in atoms, molecules, solids, and a variety of model systems. AFQMC has recently witnessed remarkable growth, especially as a tool for electronic structure computations in real materials. The method has demonstrated excellent accuracy across a variety of correlated electron systems. Taking the form of stochastic evolution in a manifold of non-orthogonal Slater determinants, the method resembles an ensemble of density-functional theory (DFT) calculations in the presence of fluctuating external potentials. Its computational cost scales as a low-power of system size, similar to the corresponding independent-electron calculations. Highly efficient and intrinsically parallel, AFQMC is able to take full advantage of contemporary high-performance computing platforms and numerical libraries. In this review, we provide a self-contained introduction to the exact and constrained variants of AFQMC, with emphasis on its applications to the electronic structure in molecular systems. Representative results are presented, and theoretical foundations and implementation details of the method are discussed., Comment: 22 pages, 11 figures

We address the computation of ground-state properties of chemical systems and realistic materials within the auxiliary-field quantum Monte Carlo method. The phase constraint to control the fermion phase problem requires the random walks in Slater determinant space to be open-ended with branching. This in turn makes it necessary to use back-propagation (BP) to compute averages and correlation functions of operators that do not commute with the Hamiltonian. Several BP schemes are investigated and their optimization with respect to the phaseless constraint is considered. We propose a modified BP method for the computation of observables in electronic systems, discuss its numerical stability and computational complexity, and assess its performance by computing ground-state properties for several substances, including constituents of the primordial terrestrial atmosphere and small organic molecules., Comment: 13 pages, 9 figures

We present numerical results for the equation of state of an infinite chain of hydrogen atoms. A variety of modern many-body methods are employed, with exhaustive cross-checks and validation. Approaches for reaching the continuous space limit and the thermodynamic limit are investigated, proposed, and tested. The detailed comparisons provide a benchmark for assessing the current state of the art in many-body computation, and for the development of new methods. The ground-state energy per atom in the linear chain is accurately determined versus bondlength, with a confidence bound given on all uncertainties., Comment: 29 pages, 16 figures -- revised and updated version