Ground-State Properties of the Hydrogen Chain: Dimerization, Insulator-to-Metal Transition, and Magnetic Phases
Mário Motta, Claudio Genovese, Fengjie Ma, Zhi‐Hao Cui, Randy C. Sawaya, Garnet Kin‐Lic Chan, Natalia Chepiga, Phillip Helms, Carlos A. Jiménez-Hoyos, Andrew J. Millis, Ushnish Ray, Enrico Ronca, Hao Shi, Sandro Sorella, Edwin Miles Stoudenmire, Steven R. White, Shiwei Zhang
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 they require solving the grand-challenge problem of the many-electron Schrödinger 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 quantummechanical 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.
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 they require solving the grand-challenge problem of the many-electron Schrödinger 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 quantummechanical 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.
The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. In this work, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.
The linear chain of hydrogen atoms, a basic prototype for the transition from a metal to Mott insulator, is studied with a recent density functional theory model functional for nondynamic and strong correlation. The computed cohesive energy curve for the transition agrees well with accurate literature results. The variation of the electronic structure in this transition is characterized with a density functional descriptor that yields the atomic population of effectively localized electrons. These new methods are also applied to the study of the Peierls dimerization of the stretched even-spaced Mott insulator to a chain of H 2 molecules, a different insulator. The transitions among the two insulating states and the metallic state of the hydrogen chain system are depicted in a semiquantitative phase diagram. Overall, we demonstrate the capability of studying strongly correlated materials with a mean-field model at the fundamental level, in contrast to the general pessimistic view on such a feasibility.
Abstract: The metal–insulator transition (MIT) constitutes one of the most profound and persistent problems in condensed matter physics, sitting at the intersection of electronic structure theory, many-body physics, and emergent quantum phenomena. Conventional band theory, which successfully classifies materials as metals or insulators according to the filling of their electronic bands, fails dramatically when confronted with a class of compounds transition metal oxides chief among them that ought to be metallic by band-counting arguments yet are robustly insulating in the laboratory. The resolution of this paradox required abandoning the independent-electron picture entirely and recognizing that strong Coulomb repulsion between electrons can, by itself, drive localization. This insight, formalized by Mott and later distilled into the Hubbard model, established correlation-driven insulation as a phenomenon distinct from band insulation. The present paper offers a theoretical exploration of the MIT problem, tracing its conceptual development from the Mott criterion through the Hubbard Hamiltonian, the Brinkman–Rice scenario of quasiparticle mass divergence, and the watershed development of dynamical mean-field theory (DMFT), which rendered the Mott transition quantitatively tractable for the first time. We argue that the competition between the electronic bandwidth and the on-site interaction encoded in the ratio U/W provides the unifying organizing principle for the field. Building toward the modern frontier, we examine how the Mott transition connects to quantum phase transitions and quantum criticality, where zero-temperature instabilities governed by quantum rather than thermal fluctuations produce non-Fermi-liquid behavior and anomalous scaling. Case studies of vanadium oxides, rare-earth nickelates, and doped cuprates illustrate the diversity and subtlety of real correlated materials. We conclude that the MIT remains a central and unresolved problem whose continued study is increasingly entangled with the development of quantum materials and quantum technologies. Keywords: Metal–Insulator Transition; Mott Insulator; Strongly Correlated Electrons; Hubbard Model; DMFT; Quantum Criticality; Quantum Materials; Condensed Matter Physics
Quantum materials exhibit a rich dynamic of physical parameters, which, when combined, can lead to entirely different behaviors. These parameters constantly compete with each other, with the most influential parameters determining the state of the system. For example, in the case of metal–insulator transitions, electron–electron interactions compete with the kinetic energy of the electrons and disorder. Understanding these complex dynamics is crucial for both fundamental physics and the development of novel technological applications, particularly given the role of disorder in tuning critical temperatures, a property with significant potential benefit in the fabrication of new devices where temperature requirements are still the bottleneck. In this article, properties of the Mott metal–insulator transition within disordered electron systems are explored using the disordered Hubbard model, the simplest Hamiltonian for capturing the metal–insulator transition. The model solutions are obtained using the self-consistent statistical dynamical mean-field theory (statDMFT). statDMFT incorporates local electronic correlation effects while allowing for Anderson localization due to disorder.
Developing accurate and computationally efficient methods to calculate the electronic structure and total energy of correlated-electron materials has been a very challenging task in condensed matter physics and materials science. Recently, we have developed a correlation matrix renormalization (CMR) method which does not assume any empirical Coulomb interaction U parameters and does not have double counting problems in the ground-state total energy calculation. The CMR method has been demonstrated to be accurate in describing both the bonding and bond breaking behaviors of molecules. In this study, we extend the CMR method to the treatment of electron correlations in periodic solid systems. Using a linear hydrogen chain as a benchmark system, we show that the results from the CMR method compare very well with those obtained recently by accurate quantum Monte Carlo (QMC) calculations. We also study the equation of states of three-dimensional crystalline phases of atomic hydrogen. We show that the results from the CMR method agree much better with the available QMC data in comparison with those from density functional theory and Hartree-Fock calculations.
Correlated electron systems present a profound challenge in condensed matter physics due to the intricate interplay between charge, spin, orbital, and lattice degrees of freedom. Understanding their rich macroscopic phase diagrams, which include phenomena like high-temperature superconductivity, colossal magnetoresistance, and Mott insulation, necessitates a bridge between the microscopic quantum interactions and emergent collective behavior. This paper introduces and formalizes the concept of quantum fugacity as a critical thermodynamic parameter capable of establishing this connection. Extending the classical notion of fugacity—a measure of a real gas's effective partial pressure—to the quantum many-body realm, quantum fugacity accounts for both particle density and the profound influence of quantum statistics and strong electron correlations. We develop a theoretical framework to calculate quantum fugacity from microscopic Hamiltonians, demonstrating its utility in mapping quantum phase transitions and discerning the driving forces behind symmetry-breaking phases. Through this approach, we show how variations in quantum fugacity serve as indicators for the onset of novel phases, providing a unified perspective on the phase behavior of diverse correlated electron materials. This concept offers a powerful new tool for designing and predicting materials with desired quantum functionalities.