The physics of strongly correlated materials poses one of the most challenging problems in condensed-matter sciences. Standard approximations applicable to wide classes of materials such as the local density approximation fail, due to the importance of the Coulomb repulsion between localized electrons. Instead, we resort to non-perturbative many-body methods. The calculations are, however, only feasible for rather small model systems. The full Hamiltonian of a real material is approximated by a model Hamiltonian comprising only the most important electronic degrees of freedom, while the effect of all other electrons is included in an average way by renormalizing the parameters. Realistic calculations of strongly correlated materials need to include sufficiently many of these electronic degrees of freedom. The new generation of massively parallel supercomputers allows for these realistic calculations. However, exploiting their computational power requires newly devised algorithms. As a solver we use the Lanczos method, which needs the full many-body state of the correlated system. It is thus limited by the available main memory. The foremost problem for a distributed-memory implementation is that the multiplication of the Hamiltonian to the many-body state leads to highly non-local memory access patterns. A simple yet important observation leads to an efficient solution: in the kinetic term of the Hamiltonian the electron-spin is conserved. Thus, writing the many-body vector as a matrix v(iup,idn), where the indices label spin-configurations, we find that the hopping term only connects vector elements that differ in one index. Hence, storing entire slices v(iup,:) on one node, the kinetic term for the spin-down electrons is local to that process. After transposing v, the same is true for the hopping of the spin-up electrons. With an efficient matrix transpose, implemented via MPI collectives, we thus obtain a version of the Lanczos method that efficiently well on JUBL, Jülich's Blue Gene/L. Using the code we treated two-orbital Hubbard Models which show interesting Mott-band insulator transitions. Kohn's criterium which determines the phase of a system by ground-state properties alone, namely the response of the ground-state energies to a vector potential, was employed. This vector potential can be equivalently regarded as a change in boundary conditions, requiring complex wave vectors. We found that half-filled systems become metallic in the region between a band and Mott insulating phase. To extrapolate to infinite system sizes and therefore getting access to an arbitrarily high resolution in reciprocal space, we use cluster perturbation theory (CPT). It works by diagonalizing a finite cluster exactly and then treating the hopping between identical clusters perturbatively in strong coupling perturbation theory. This leads to an effectively infinite chain. Due to our efficient Lanczos implementation it is possible to calculate relatively large clusters, minimizing finite-size effects. The combination of our Lanczos code and CPT allowed for a treatment of the the organic metal TTF-TCNQ. We were able to solve a long-standing problem in the interpretation of its experimental spectra. Being an interesting one-dimensional metal TTF-TCNQ has been studied thoroughly for more than 30 years. Theoretical estimates as well as experimental measurements suggest values for the parameters of the Hubbard model t and U. With these values Hubbard model calculations show qualitatively the same features as APRES experiments, namely the signatures of spin-charge separation at the Gamma-point. But in order to fit the experimental data to the numerical calculations t has to be doubled, which translates into a broadening of the spectra. This ad-hoc adjustment is not only unsatisfactory but also the temperature dependence of the spectral function is incorrectly reproduced. We resolved this problem by properly including nearest-neighbor Coulomb interactions. Recent DFT calculations show that indeed the commonly used values of t and U are correct. But Coulomb interactions between electrons on neighboring sites are also significant and must not be neglected. Thus, we studied the effect of the next-neighbor interaction parameter V being at about U/2. Including V in our Lanczos calculations for obtaining the angular-integrated spectral function shows that its effect is to broaden the spectra, comparable to increasing the value of t. Moreover, we found that doubling t indeed mimics the effect of Vapprox U/2. Surprisingly, the broadening of the spectrum as a result of V can be understood in first-order Rayleigh-Schrödinger perturbation theory.
The physics of strongly correlated materials poses one of the most challenging problems in condensed-matter sciences. Standard approximations applicable to wide classes of materials such as the local density approximation fail, due to the importance of the Coulomb repulsion between localized electrons. Instead, we resort to non-perturbative many-body methods. The calculations are, however, only feasible for rather small model systems. The full Hamiltonian of a real material is approximated by a model Hamiltonian comprising only the most important electronic degrees of freedom, while the effect of all other electrons is included in an average way by renormalizing the parameters. Realistic calculations of strongly correlated materials need to include sufficiently many of these electronic degrees of freedom. The new generation of massively parallel supercomputers allows for these realistic calculations. However, exploiting their computational power requires newly devised algorithms. As a solver we use the Lanczos method, which needs the full many-body state of the correlated system. It is thus limited by the available main memory. The foremost problem for a distributed-memory implementation is that the multiplication of the Hamiltonian to the many-body state leads to highly non-local memory access patterns. A simple yet important observation leads to an efficient solution: in the kinetic term of the Hamiltonian the electron-spin is conserved. Thus, writing the many-body vector as a matrix v(iup,idn), where the indices label spin-configurations, we find that the hopping term only connects vector elements that differ in one index. Hence, storing entire slices v(iup,:) on one node, the kinetic term for the spin-down electrons is local to that process. After transposing v, the same is true for the hopping of the spin-up electrons. With an efficient matrix transpose, implemented via MPI collectives, we thus obtain a version of the Lanczos method that efficiently well on JUBL, Jülich's Blue Gene/L. Using the code we treated two-orbital Hubbard Models which show interesting Mott-band insulator transitions. Kohn's criterium which determines the phase of a system by ground-state properties alone, namely the response of the ground-state energies to a vector potential, was employed. This vector potential can be equivalently regarded as a change in boundary conditions, requiring complex wave vectors. We found that half-filled systems become metallic in the region between a band and Mott insulating phase. To extrapolate to infinite system sizes and therefore getting access to an arbitrarily high resolution in reciprocal space, we use cluster perturbation theory (CPT). It works by diagonalizing a finite cluster exactly and then treating the hopping between identical clusters perturbatively in strong coupling perturbation theory. This leads to an effectively infinite chain. Due to our efficient Lanczos implementation it is possible to calculate relatively large clusters, minimizing finite-size effects. The combination of our Lanczos code and CPT allowed for a treatment of the the organic metal TTF-TCNQ. We were able to solve a long-standing problem in the interpretation of its experimental spectra. Being an interesting one-dimensional metal TTF-TCNQ has been studied thoroughly for more than 30 years. Theoretical estimates as well as experimental measurements suggest values for the parameters of the Hubbard model t and U. With these values Hubbard model calculations show qualitatively the same features as APRES experiments, namely the signatures of spin-charge separation at the Gamma-point. But in order to fit the experimental data to the numerical calculations t has to be doubled, which translates into a broadening of the spectra. This ad-hoc adjustment is not only unsatisfactory but also the temperature dependence of the spectral function is incorrectly reproduced. We resolved this problem by properly including nearest-neighbor Coulomb interactions. Recent DFT calculations show that indeed the commonly used values of t and U are correct. But Coulomb interactions between electrons on neighboring sites are also significant and must not be neglected. Thus, we studied the effect of the next-neighbor interaction parameter V being at about U/2. Including V in our Lanczos calculations for obtaining the angular-integrated spectral function shows that its effect is to broaden the spectra, comparable to increasing the value of t. Moreover, we found that doubling t indeed mimics the effect of Vapprox U/2. Surprisingly, the broadening of the spectrum as a result of V can be understood in first-order Rayleigh-Schrödinger perturbation theory.
Understanding the physics of strongly correlated materials is one of the grand-challenges in condensed-matter physics.Simple approximations such as the local density approximation fail, due to the importance of the Coulomb repulsion between localized electrons.Instead we have to resort to non-perturbative many-body techniques.Such calculations are, however, only feasible for quite small model systems.This means that the full Hamiltonian of a real material has to be approximated by a model Hamiltonian comprising only the most important electronic degrees of freedom, while the effect of all other electrons can merely be included in an average way.Realistic calculations of strongly correlated materials need to include as many of the electronic degrees of freedom as possible.We use the Lanczos method for calculating ground state and dynamical properties of such model Hamiltonians.In the method we have to handle the full many-body state of the correlated system.The method is thus limited by the available main memory.The principal problem for a distributed-memory implementation is that the central routine of the code, the application of the Hamiltonian to the many-body state, leads, due to the kinetic energy term, to very non-local memory access.We describe a solution for this problem and show that the new algorithm scales very efficiently on BlueGene/L.Moreover, our approach is not restricted to correlated electrons but can also be used to simulate quantum spin systems relevant for quantum computing.
In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
Since the advent of high-Tc cuprate superconductors in 1986, strongly correlated electron systems have attracted much attention. Since the cuprates are essentially two-dimensional, low-dimensional systems have moved into the focus of condensed-matter theory. From a theoretical point of view, one-dimensional systems are of particular interest because there are exact numerical and analytical methods which permit detailed studies and deep insights into the many-body problem. In the first part of this thesis, using the numerically exact methods Exact Diagonalization and the Density-Matrix Renormalization Group (DMRG), we investigate the properties of the Tomonaga–Luttinger liquid which is the generic metallic state of matter in one dimension. In particular, we concentrate on the investigation of the so-called Tomonaga–Luttinger liquid parameter which determines the critical exponent α for the density of states near the Fermi energy. Experimental results for some quasi one-dimensional materials report α≳1, which would imply Kρ≲0.17, a value which cannot be reconciled with the bare Hubbard model where KρH≥0.5, i.e., αH≲1/8. We develop new accurate numerical methods to obtain Kρ and investigate how to obtain such small values for Kρ for slightly doped charge-density-wave insulators. In the second part of this thesis, using the Thermodynamic Bethe Ansatz (TBA) as exact analytical method, we investigate the one-dimensional Hubbard model in the spindisordered regime, which is characterized by the temperature being much larger than the magnetic energy scale but small compared to the Mott–Hubbard gap. Our study is motivated by the controversy about the Mott–Hubbard insulator in infinite dimensions whose ground state is also spin-disordered. In this system the determination of the precise value of the critical interaction strength Uc where the Mott–Hubbard gap closes is still unsolved. Therefore, we provide an example of a Hubbard-type model with a disordered spin background which can be solved exactly. The thermodynamics of our model can be understood in terms of gapped charged excitations with an effective dispersion which are decoupled form the spin degrees of freedom; the latter contribute only entropically. An interpretation of this regime in terms of a putative interacting-electron system at zero temperature leads to a metal-insulator transition at a finite interaction strength above which the gap opens linearly. Our exact results indicate that the strong-coupling expansion of the ground-state energy cannot be used to locate Uc. However, the strong-coupling expansion of the gap permits a reliable extrapolation of the critical interaction strength.
Achieving materials-specific predictions for large, realistic molecules and solids with strong electron correlations remains a long-standing challenge in quantum chemistry. This dissertation develops scalable ab initio quantum many-body methods for crystalline materials, demonstrating that rigorous optimization of theoretical approximations and numerical algorithms enables predictive computation even in strongly correlated systems. The work tackles three interconnected challenges, namely describing material-specific properties of strongly correlated superconductors, achieving linear scaling with k-points for large-scale periodic calculations, and establishing reliable methods for electron-phonon interactions. The following chapters address each challenge in turn. Chapter 2 demonstrates that ab initio density matrix embedding theory, combined with symmetry-breaking strategies for superconducting order parameters, reproduces experimental trends in cuprate superconductivity, including pressure and layer effects on transition temperatures. By directly solving the electronic Schrödinger equation from material structures, these calculations capture multi-orbital covalency effects that simplified model Hamiltonians neglect, revealing their essential role in the spin fluctuations driving pairing. This predictive capability establishes a foundation for efficient computational screening of superconducting materials. Extending these calculations to larger k-point meshes and broader phase diagram explorations, however, demands more efficient computational algorithms. To overcome this computational bottleneck, Chapter 3 introduces interpolative separable density fitting framework. This development enables thermodynamic limit calculations with up to 1000 k-points, as validated on diamond, carbon dioxide, nickel monoxide, and cuprate systems. Integration with density matrix embedding theory and local natural orbital methods further yields converged ground-state energies from accurate correlated wavefunction methods such as CCSD(T) for crystalline materials in the thermodynamic limit. With computational efficiency established, the remaining challenge lies in extending the theoretical framework beyond purely electronic degrees of freedom. Chapter 4 addresses this need by establishing exponential ansatz wavefunctions as competitive approaches for electron-phonon systems. Unlike previous treatments based on perturbation theory, systematic benchmarking of coupled cluster theory and variational Lang-Firsov methods on the Holstein model demonstrates, for the first time, that these approaches accurately describe polaron formation across coupling regimes when combined with appropriate reference state optimization. These results provide a foundation for incorporating phonon effects into ab initio studies of superconducting materials. In addition to these methodological advances, GPU acceleration and MPI parallelization enhance computational throughput, enabling efficient large-scale calculations. Together, these developments constitute a computational framework for predictive ab initio quantum many-body calculations in crystalline materials, opening new avenues for understanding and designing strongly correlated quantum materials.
The main focus of this thesis is the detailed investigation of computational methods to tackle strongly correlated materials in which a rich variety of exotic phenomena are found. A many-body problem with sizable electronic correlations can no longer be explained by independent-particle approximations such as density functional theory (DFT) or tight-binding approaches. The influence of an electron to the others is too strong for each electron to be treated as an independent quasiparticle and consequently those standard band-structure methods fail even at a qualitative level. One of the most powerful approaches for strong correlations is the dynamical mean-field theory (DMFT), which has enlightened the understanding of the Mott transition based on the Hubbard model. For realistic applications, the dynamical mean-field theory is combined with various independent-particles approaches. The most widely used one is the DMFT combined with the DFT in the local density approximation (LDA), so-called LDA+DMFT. In this approach, the electrons in the weakly correlated orbitals are calculated by LDA while others in the strongly correlated orbitals are treated by DMFT. Recently, the method combining DMFT with Hedin's GW approximation was also developed, in which the momentum-dependent self-energy is also added. In this thesis, we discuss the application of those methodologies based on DMFT. First, we apply the dynamical mean-field theory to solve the 3-dimensional Hubbard model in Chap. 3. In this application, we model the interface between the thermodynamically coexisting metal and Mott insulator. We show how to model the required slab geometry and extract the electronic spectra. We construct an effective Landau free energy and compute the variation of its parameters across the phase diagram. Finally, using a linear mixture of the density and double-occupancy, we identify a natural Ising order parameter which unifies the treatment of the bandwidth and filling controlled Mott transitions. Secondly, we study the double-counting problem, a subtle issue that arises in LDA+DMFT. We propose a highly precise double-counting functional, in which the intersection of LDA and DMFT is calculated exactly, and implement a parameter-free version of the LDA+DMFT that is tested on one of the simplest strongly correlated systems, the H2 molecule. We show that the exact double-counting treatment along with a good DMFT projector leads to very accurate and total energy and excitation spectrum of H2 molecule. Finally, we implement various versions of GW+DMFT, in its fully self-consistent way, one shot GW approximation, and quasiparticle self-consistent scheme, and studied how well these combined methods perform on H2 molecule as compared to more established methods such as LDA+DMFT. We found that most flavors of GW+DMFT break down in strongly correlated regime due to causality violation. Among GW+DMFT methods, only the self-consistent quasiparticle GW+DMFT with static double-counting, and a new method with causal double-counting, correctly recover the atomic limit at large H-atom separation. While some flavors of GW+DMFT improve the single-electron spectra of LDA+DMFT, the total energy is best predicted by LDA+DMFT, for which the exact double-counting is known, and is static.
[ACCESS RESTRICTED TO THE UNIVERSITY OF MISSOURI AT AUTHOR'S REQUEST.] In this dissertation, we have presented our research on new developments in the Gutzwiller variational method as well as the applications of the method to two specific correlated electron systems of current interest: (i) the spinless fermion model with non-local Coulomb interaction in one dimension and (ii) the pressure-induced insulator-metal transition in the colossal magneto-resistive compound LaMnO3. Our results for the ground-state energy of the spinless fermion model are obtained by careful enumeration of the many-body configurations of the system and calculating the Gutzwiller reduction factor to the kinetic energy for general filling. We resolve the inconsistency in the existing literature for this problem as we find that our result for the half-filled case agrees with one of the three papers and attribute the discrepancy between them to the incorrect counting method and the assumption that the Gutzwiller approximation is equivalent to the mean-field slave-boson approach, which is although correct for the on-site Hubbard model, turns out not to be so in the present case. Compared with the exact diagonalization results, we show that the Gutzwiller method indeed offers a better solution than the slave-boson approach. As the second problem, we have studied the ground-state properties of paramagnetic LaMnO3 under hydrostatic pressure at room temperature. We apply the Gutzwiller method to an extended Hubbard model for spinless eg electrons on Mn3+ ion including all important ingredients such as cohesive energy as well as two main competing Coulomb and Jahn-Teller (JT) interactions. The plot of the energy of the paramagnetic LaMnO3 as a function of volume clearly confirms that 1) contrary to the existing accepted theoretical picture for vanishing lattice distortion as the metallic state emerges, the JT distortion survives even beyond the transition pressure into the metallic region and 2) the application of pressure leads to a structural phase separation and breaks the material into two blocks composed of domains of undistorted (metallic character) and distorted (insulating character) MnO6 octahedra; both of these findings have been already observed in experiments with no theoretical explanations. As the key unifying connection between both of these findings we suggest that the pressure-induced transition and the phase coexistence is indeed percolative. Finally, applying the idea of universality in percolation quantities near the threshold, we support our hypothesis by proposing a scaling law for the pressure-dependence of the resistance and use it as a best fit to recently-acquired experimental data. We find good agreement between our proposed model and the experimental results, confirming that the pressure-induced insulator-metal transition in LaMnO3 is indeed percolative.