The metal–insulator transition (MIT) is a fundamental phenomenon in condensed matter physics and a hallmark of strong electronic correlations. Hydrogen-based systems offer a simple yet powerful model for investigating the MIT, as their insulating behavior arises purely from electron–electron interactions. In this work, we study finite hydrogen clusters with cubic geometries using Natural Orbital Functional Theory (NOFT), a method capable of accurately describing correlated systems beyond mean-field approaches. We focus on two key signatures of the MIT: the fundamental energy gap and the harmonic average of the atomic one-particle reduced density matrix. Our results show that NOFT captures the transition from insulating to metallic behavior as the interatomic distance decreases. By extrapolating the energy gap to the thermodynamic limit, we estimate a critical distance rc ~1.2 Å, in excellent agreement with quantum Monte Carlo benchmarks. These findings demonstrate the reliability of NOFT for describing strong correlation effects in large-scale models.
The metal–insulator transition (MIT) is a fundamental phenomenon in condensed matter physics and a hallmark of strong electronic correlations. Hydrogen-based systems offer a simple yet powerful model for investigating the MIT, as their insulating behavior arises purely from electron–electron interactions. In this work, we study finite hydrogen clusters with cubic geometries using Natural Orbital Functional Theory (NOFT), a method capable of accurately describing correlated systems beyond mean-field approaches. We focus on two key signatures of the MIT: the fundamental energy gap and the harmonic average of the atomic one-particle reduced density matrix. Our results show that NOFT captures the transition from insulating to metallic behavior as the interatomic distance decreases. By extrapolating the energy gap to the thermodynamic limit, we estimate a critical distance rc ~1.2 Å, in excellent agreement with quantum Monte Carlo benchmarks. These findings demonstrate the reliability of NOFT for describing strong correlation effects in large-scale models.
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.
Describing correlated electron systems near phase transitions has been a major challenge in computational condensed-matter physics. In this paper, we apply highly accurate fixed-node quantum Monte Carlo techniques, which directly work with many-body wave functions and simulate electron correlations, to investigate the metal-to-insulator transition of a correlated hydrogen lattice. By calculating spin and charge properties, and analyzing the low-energy Hilbert space, we identify the transition point and identify order parameters that can be used to detect the transition. Our results provide a benchmark for density functional theories seeking to treat correlated electron systems.
The clear distinction between metals and insulators is a well established and appartently simple concept which can be clarified by studying the behavior of the electrical conductivity. However, a theoretical description of the processes leading to one of the two extremes is no trivial task. The widely accepted band structure theory relies on an one-electron picture to characterize the metallic/insulating behavior of a material. However, in spite of the fact that it can be used to make successful predictions, this description does not capture one of the essential elements characterizing the nature of an electronic quantum system, i.e. the electron correlation. The key to the different behaviors is clearly contained in the correlated many- electron wave function, but the analysis of such a high-dimensional object in configuration space can be quite cumbersome. Therefore the development of powerful tools to simplify this task is necessary. The total position-spread (TPS) tensor, which was derived by Walter Kohn’s theory of the insulating state, plays a key role in the field of metal-insulator transitions (MITs). This tool provides information about the electron localization, which is strongly connected to the electrical conductivity. Therefore, in this work particular attention was paid to the calculation of the TPS tensor of small and extended systems. Moreover, part of this thesis focuses on a new formalism, which allows the study of spin mobility and spin entanglement and provides therefore deeper insight into the effects of electron correlation. In the framework of quantum information theory (QIT), the quantum entanglement can be quantified by means of measures such as the von Neumann entropy. Also, the analysis of the elements of the two-orbital reduced density matrix provides an overview of the complexity of the many-electron wave function. These insights were used in this work for the analysis of the metallic/insulating character of a system. However, in order to perform such investigastions, the problem of achieving an accurate description of the many- electron wave function has to be dealt with. Despite the tremedous progress in the efficiency of quantum-chemical methods, the treatement of extended systems using standard approaches is still a difficult (if not unfeasible) task. Therefore the development and testing of sophisticated, low-scaling methods is currently the subject of many investigations. In this work, two approaches were employed: 1) the ab-initio density matrix renormalization group (DMRG) approach which, based on a tensor product ansatz, provides very precise results and exploits QIT to improve its performance; 2) the method of increments (MoI) which, using a many-body expansion in terms of localized orbitals, yields accurate correlation energies. By using both methods on strongly correlated one-dimesional model systems, we worked on a new formalism of the method of increments which can be used to describe the whole dissociation curve without inconvenient discontinuities. Finally, aiming to calculate the total position-spread tensor of extended systems, we worked on the application of the MoI for quantities other than the correlation energy.
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.
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
For a long time strong electronic correlations in metals have mainly been associated with Mottness, the proximity to a Mott metal-insulator transition (MIT), where large Coulomb interactions induce the localization of charges. However, triggered by the discovery of the iron-based superconductors about ten years ago, it was realized that multi-orbital materials with only moderate Coulomb but sizeable Hund’s rule interactions – so-called Hund metals – allow for a distinct screening mechanism towards strong correlations: Hundness. Here, Hund’s rule constrains the spin rather than the charge dynamics. This discovery led to a vividly debated fundamental issue in the field of strongly correlated condensed matter systems, which is the main topic of the present thesis: what is the origin of strong correlations in the normal phase of Hund metals, Mottness or Hundness? And what are their decisive fingerprints? The goal of this dissertation is twofold. First, we present and advance our method: the numerical renormalization group (NRG) as viable real-frequency multi-band impurity solver for dynamical mean-field theory (DMFT), a common approach to tackle strongly correlated systems. Second, we apply DMFT+NRG to shed light on the Hund-metal problem raised above. In the first part of this thesis we present our state-of-the-art NRG solver, which offers direct access to data with unprecedented real-frequency spectral resolution at arbitrarily low energies and temperatures in contrast to commonly used Quantum Monte Carlo solvers. It is based on matrix product states and exploits non-abelian symmetries to reduce numerical costs. In the case of orbital symmetry, this allows us to treat multi-band models with more than two bands, and thus to tackle the Hund-metal problem for the first time with NRG. For multi-band models without orbital symmetry, an “interleaved” scheme of NRG (iNRG) was recently developed, dramatically increasing the numerical efficiency. Remarkably, the accuracy of iNRG is comparable to standard NRG, as we reveal in a detailed study. This finding establishes iNRG as a promising DMFT solver for material-specific model simulations. In the second part of this thesis we study a minimal toy model for Hund metals with DMFT+NRG, the orbital-symmetric three-band Hubbard-Hund model (3HHM) close to a lattice filling of 1/3. Our major insight is “spin-orbital separation” (SOS), a Hund’s-ruleinduced two-stage Kondo-type screening process, in which orbital screening occurs at much higher energies than spin screening. In Hund metals, i.e. far from a MIT phase boundary, SOS thus causes large electron masses by strongly reducing the coherence scale below which a Fermi liquid is formed. Further, it opens up a broad incoherent and strongly particle-hole asymmetric intermediate energy regime that reaches up to bare excitation scales. This SOS regime shows fractional power-law behavior and is characterized by resilient “Hund quasiparticles” with itinerant orbital degrees of freedom coupled non-trivially to quasi-free large spins. At zero temperature, the local density of states exhibits a two-tier quasiparticle peak on top of a broad incoherent background. In contrast, in Mott-correlated metals, i.e. close to the MIT phase boundary, the SOS regime becomes negligibly small and the Hubbard bands are well separated. These findings lead to distinct signatures of Hundness and Mottness in the temperature dependence of ARPES spectra, static local susceptibilities, resistivity, thermopower and entropy, many of which were also found in realistic simulations of the archetypal Hund- and Mott-correlated materials, Sr2RuO4 and V2O3. In summary, we provide evidence that and elucidate how Hundness evokes strong correlation effects in Hund metals. This might help to better interpret experimental results and guide superconducting theories.