We have studied the metal--insulator-like transition in pseudo-one-dimensional systems, i.e., lithium and beryllium rings, through the ab initio density matrix renormalization group (DMRG) method. Performing accurate calculations for different interatomic distances and using quantum information theory, we investigated the changes occurring in the wave function between a metallic-like state and an insulating state built from free atoms. We also discuss entanglement and relevant excitations among the molecular orbitals in the Li and Be rings and show that the transition bond length can be detected using orbital entropy functions. Also, the effect of different orbital bases on the effectiveness of the DMRG procedure is analyzed comparing the convergence behavior.
We have studied the metal--insulator-like transition in pseudo-one-dimensional systems, i.e., lithium and beryllium rings, through the ab initio density matrix renormalization group (DMRG) method. Performing accurate calculations for different interatomic distances and using quantum information theory, we investigated the changes occurring in the wave function between a metallic-like state and an insulating state built from free atoms. We also discuss entanglement and relevant excitations among the molecular orbitals in the Li and Be rings and show that the transition bond length can be detected using orbital entropy functions. Also, the effect of different orbital bases on the effectiveness of the DMRG procedure is analyzed comparing the convergence behavior.
We have studied transition metal clusters from a quantum information theory perspective using the density-matrix renormalization group (DMRG) method. We demonstrate the competition between entanglement and interaction localization and discuss the application of the configuration interaction-based dynamically extended active space procedure, which significantly reduces the effective system size and accelerates the speed of convergence for complicated molecular electronic structures. Our results indicate the importance of taking entanglement among molecular orbitals into account in order to devise an optimal DMRG orbital ordering and carry out efficient calculations on transition metal clusters. Apart from these algorithmic observations, which lead to a recipe for black-box DMRG calculations, our work provides physical understanding of electron correlation in molecular and cluster structures in terms of entropy measures of relevance also to recent work on tensor-network representations of electronic states. We also identify those molecular orbitals which are highly entangled and discuss the consequences for chemical bonding and for the structural transition from an dioxygen binding copper cluster to an bis-oxygen-bridged system with broken O--O bond.
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.
Abstract Context: In the realm of quantum chemistry, the accurate prediction of electronic structure and properties of nanostructures remains a formidable challenge. Density Functional Theory (DFT) and Density Matrix Renormalization Group (DMRG) have emerged as two powerful computational methods for addressing electronic correlation effects in diverse molecular systems. We compare ground-state energies (e0), density profiles (n) and average entanglement entropies (Sˉ) in metals, insulators and at the transition from metal to insulator, in homogeneous, superlattices and harmonically confined chains described by the fermionic one-dimensional Hubbard model. While for the homogeneous systems there is a clear hierarchy between the deviations, D, and all the deviations decrease with the chain size; for superlattices and harmonical confinement the relation among the deviations is less trivial and strongly dependent on the superlattice structure and the confinement strength considered. For the superlattices, in general increasing the number of impurities in the unit cell represents less precision on the DFT calculations. For the confined chains, DFT performs better for metallic phases, while the highest deviations appear for the Mott and band-insulator phases. This work provides a comprehensive comparative analysis of these methodologies, shedding light on their respective strengths, limitations, and applications. Methods: The DFT calculations were performed using the standard Kohn-Sham scheme within the BALDA approach. It integrated the numerical Bethe-Ansatz (BA) solution of the Hubbard model as the homogeneous density functional within a local-density approximation (LDA) for the exchange-correlation energy. The DMRG algorithms were implemented using the ITensor library, which is based on the Matrix Product States (MPS) ansatz. The calculations were performed until the energy reaches convergence of at least 10-8.
We investigate three different types of quantum phase transition occurring in quasi one-dimensional systems theoretically and numerically. First, we study the band-insulator to Mott-insulator transition occurring in charge-transfer complexes, for which the half-filled one-dimensional ionic Hubbard model is considered to be the prototype model. The study is carried out by first deriving an effective spin-one model, and then studying the model numerically using the density matrix renormalization group. We perform a careful finite-size scaling analysis of the mass gaps, order-parameters, and relative susceptibility. We confirm the existence of two quantum critical points. Analysis of the critical exponents confirms that the band-insulator-to-spontaneously-dimerized phase transition belong to the 2D Ising class. The spontaneously dimerized phase undergoes a phase transition to the Mott-insulator which is an infinite-order. Second, we investigate the Mott metal-insulator transition for the half-filled Hubbard model with both nearest-neighbor t and next-nearest-neighbor t′ hopping terms. We study the model using the bosonization approach and density matrix renormalization group simulations. An effective low-energy Hamiltonian that describes the insulator-metal transition is derived. We present results of density matrix renormalization group calculations of spin and charge distribution in various sectors of the phase diagram. The numerical results support the picture derived from the effective theory and give evidence for the complete separation of the transitions involving the spin and the charge degrees of freedom. Finally, we investigate quantum phase transitions phases in low-dimensional fermionic and spin models that go from uniform to spatially inhomogeneous, i.e., dimerized, trimerized, or incommensurate, phases. We propose a new approach based on studying the length dependence of the von Neumann entropy and its corresponding Fourier spectrum for finite segments in the ground state of finite chains. Peaks at a nonzero wave vector are indicators of oscillatory behavior in decaying correlation functions and also provide significant information about certain relevant features of the excitation spectrum; in particular, they can identify the wave vector of soft modes in critical models.
The accurate first-principles calculation of relative energies of transition metal complexes and clusters is still one of the great challenges for quantum chemistry. Dense lying electronic states and near degeneracies make accurate predictions difficult, and multireference methods with large active spaces are required. Often density functional theory calculations are employed for feasibility reasons, but their actual accuracy for a given system is usually difficult to assess (also because accurate ab initio reference data are lacking). In this work we study the performance of the density matrix renormalization group algorithm for the prediction of relative energies of transition metal complexes and clusters of different spin and molecular structure. In particular, the focus is on the relative energetical order of electronic states of different spin for mononuclear complexes and on the relative energy of different isomers of dinuclear oxo-bridged copper clusters.
This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr(2) molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mE(h) (the error from the extrapolated best value) when as many as 10,000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu(2)O(2)](2+) along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu(2)O(2)](2+) according to the one-electron density matrix analysis.