In this thesis, quantum cluster methods are used to calculate electronic properties of correlated-electron systems. A special focus lies in the determination of the ground state properties of a 3/4 filled triangular lattice within the one-band Hubbard model. At this filling, the electronic density of states exhibits a so-called van Hove singularity and the Fermi surface becomes perfectly nested, causing an instability towards a variety of spin-density-wave (SDW) and superconducting states. While chiral d+id-wave superconductivity has been proposed as the ground state in the weak coupling limit, the situation towards strong interactions is unclear. Additionally, quantum cluster methods are used here to investigate the interplay of Coulomb interactions and symmetry-breaking mechanisms within the nematic phase of iron-pnictide superconductors. The transition from a tetragonal to an orthorhombic phase is accompanied by a significant change in electronic properties, while long-range magnetic order is not established yet. The driving force of this transition may not only be phonons but also magnetic or orbital fluctuations. The signatures of these scenarios are studied with quantum cluster methods to identify the most important effects. Here, cluster perturbation theory (CPT) and its variational extention, the variational cluster approach (VCA) are used to treat the respective systems on a level beyond mean-field theory. Short-range correlations are incorporated numerically exactly by exact diagonalization (ED). In the VCA, long-range interactions are included by variational optimization of a fictitious symmetry-breaking field based on a self-energy functional approach. Due to limitations of ED, cluster sizes are limited to a small number of degrees of freedom. For the 3/4 filled triangular lattice, the VCA is performed for different cluster symmetries. A strong symmetry dependence and finite-size effects make a comparison of the results from different clusters difficult. The ground state in the weak-coupling limit is superconducting with chiral d+id-wave symmetry, in accordance to previous renormalization group approaches. In the regime of strong interactions SDW states are preferred over superconductivity and a collinaer SDW state with nonuniform spin moments on a quadrupled unit cell has the lowest grand potential. At strong coupling, inclusion of short-range quantum fluctuations turns out to favor this collinear state over the chiral phase predicted by mean-field theory. At intermediate interactions, no robust conclusion can be drawn from the results. Symmetry-breaking mechanisms within the nematic phase of the iron-pnictides are studied using a three-band model for the iron planes on a 4-site cluster. CPT allows a local breaking of the symmetry within the cluster without imposing long-range magnetic order. This is a crucial step beyond mean-field approaches to the magnetically ordered state, where such a nematic phase cannot easily be investigated. Three mechanisms are included to break the fourfold lattice symmetry down to a twofold symmetry. The effects of anisotropic magnetic couplings are compared to an orbital ordering field and anisotropic hoppings. All three mechanisms lead to similar features in the spectral density. Since the anisotropy of the hopping parameters has to be very large to obtain similar results as observed in ARPES, a phonon-driven transition is unlikely.
In this thesis, quantum cluster methods are used to calculate electronic properties of correlated-electron systems. A special focus lies in the determination of the ground state properties of a 3/4 filled triangular lattice within the one-band Hubbard model. At this filling, the electronic density of states exhibits a so-called van Hove singularity and the Fermi surface becomes perfectly nested, causing an instability towards a variety of spin-density-wave (SDW) and superconducting states. While chiral d+id-wave superconductivity has been proposed as the ground state in the weak coupling limit, the situation towards strong interactions is unclear. Additionally, quantum cluster methods are used here to investigate the interplay of Coulomb interactions and symmetry-breaking mechanisms within the nematic phase of iron-pnictide superconductors. The transition from a tetragonal to an orthorhombic phase is accompanied by a significant change in electronic properties, while long-range magnetic order is not established yet. The driving force of this transition may not only be phonons but also magnetic or orbital fluctuations. The signatures of these scenarios are studied with quantum cluster methods to identify the most important effects. Here, cluster perturbation theory (CPT) and its variational extention, the variational cluster approach (VCA) are used to treat the respective systems on a level beyond mean-field theory. Short-range correlations are incorporated numerically exactly by exact diagonalization (ED). In the VCA, long-range interactions are included by variational optimization of a fictitious symmetry-breaking field based on a self-energy functional approach. Due to limitations of ED, cluster sizes are limited to a small number of degrees of freedom. For the 3/4 filled triangular lattice, the VCA is performed for different cluster symmetries. A strong symmetry dependence and finite-size effects make a comparison of the results from different clusters difficult. The ground state in the weak-coupling limit is superconducting with chiral d+id-wave symmetry, in accordance to previous renormalization group approaches. In the regime of strong interactions SDW states are preferred over superconductivity and a collinaer SDW state with nonuniform spin moments on a quadrupled unit cell has the lowest grand potential. At strong coupling, inclusion of short-range quantum fluctuations turns out to favor this collinear state over the chiral phase predicted by mean-field theory. At intermediate interactions, no robust conclusion can be drawn from the results. Symmetry-breaking mechanisms within the nematic phase of the iron-pnictides are studied using a three-band model for the iron planes on a 4-site cluster. CPT allows a local breaking of the symmetry within the cluster without imposing long-range magnetic order. This is a crucial step beyond mean-field approaches to the magnetically ordered state, where such a nematic phase cannot easily be investigated. Three mechanisms are included to break the fourfold lattice symmetry down to a twofold symmetry. The effects of anisotropic magnetic couplings are compared to an orbital ordering field and anisotropic hoppings. All three mechanisms lead to similar features in the spectral density. Since the anisotropy of the hopping parameters has to be very large to obtain similar results as observed in ARPES, a phonon-driven transition is unlikely.
In order to understand various local phenomena related to electronic phase separation discovered recently in strongly correlated materials, quantum cluster methods are applied to a theoretical model, i.e., Hubbard model, which has successfully predicted basic properties of strongly correlated systems. Our calculations are performed at two levels: exact calculations based on the isolated small clusters via exact diagonalization and approximate calculations applied to large lattices using the variational cluster approximation (VCA).\nWe first study phase separation instabilities with the formation of a pairing gap in the two-dimensional Hubbard model for 8-site Betts clusters. The exact diagonalization method is applied to extract exact ground and excited states of the Hubbard model. The results show that the electronic states with one hole off half filling are unstable and they energetically prefer creating the spatial phase separation. The effect of next nearest hopping is also discussed in this part.\nWe studied the variation of the pairing when the out-of-plane correlation is present by adding an interacting apical site above the two-dimensional cluster plane. The calculations indicate that the out-of-plane correlation can be detrimental to the in-plane pairing effect. When it is not too strong to destroy electron pairs on the plane, the modulation of the apical site can drive the pairing gap to show a sinusoidal variation which is consistent with recent experimental discoveries.\nVariational cluster approximation is introduced to explore local properties of larger lattices. The local nematic state is studied in the square lattice using this method, which indicates that it is possible that the electronic states can locally break the C4 symmetry (90o$ rotational symmetry) to form a pattern where the electronic properties, such as charge and spin correlations, are different in x and y directions.\nThe spatial phase separation in the square and honeycomb lattices is also investigated by the variational cluster approximation. The study shows that Coulomb interaction and lattice geometry are both crucial for electronic phase separation. We propose that electronic phase separation is the result of an instability of the Fermi surface at the boundary of the first Brillouin zone of antiferromagnetic states.
Materials with strongly correlated electrons show a multitude of unconventional phenomena, such as superconductivity or magnetism, that cannot be explained in a single-particle picture and necessitate to take many-body effects into account. A route to explain the origin of these phenomena consists in constructing effective models for these materials that are then solved, often using numerical techniques. Paradigmatic models for two different classes of strongly correlated electron systems are investigated in this thesis: The Hubbard model is a paradigmatic model to describe the Mott metal-insulator transition (MIT) and the Kondo lattice model describes the low-energy physics in heavy fermion compounds. Both models are analyzed numerically using different variations of a quantum cluster technique, the variational cluster approximation (VCA). The first part focuses on the MIT, as it is investigated experimentally in quasi-two-dimensional charge-transfer salts. Its nature and associated universality class are still heavily discussed and both experiments and theoretical approaches come to different conclusions. In this thesis, the Mott transition is investigated on an anisotropic two-dimensional Hubbard model at half-filling. By using a control parameter, which induces this anisotropy, strong evidence for Mott quantum criticality is found in weakly coupled Hubbard chains. The results at zero temperature show that the second-order critical end point Tc of the interaction-driven metal-insulator transition can be tuned down to zero at strong anisotropy. Further results for the antiferromagnetic phase suggest a similar picture and motivate adding the anisotropy as a new axis in the phase diagram to account for a low-temperature critical end point of the Mott transition. The second part of this thesis focuses on unconventional phases that emerge in a paradigmatic model for heavy fermion systems, the Kondo lattice model. Using the VCA, the two-dimensional Kondo lattice model is investigated in the paramagnetic phase as well as in phases with broken symmetry. An antiferromagnetic phase at weak coupling and a phase with Kondo-singlet formation at strong coupling are analyzed. Within the antiferromagnetically ordered region, two phases with different Fermi-surface topology, which are separated by a discontinuous transition, are identified. The model is also tested for s-wave superconductivity as found recently using dynamical mean-field theory, but no indications for robust local pairing are seen within VCA. Instead, nodal d-wave superconductivity is found and analyzed in a large range of couplings and electron fillings, and its interplay with antiferromagnetic order at weak coupling is discussed. This motivates further studies on extended models that might allow for a close comparison with experiments on heavy fermion systems. Understanding the formation of unconventional phases in strongly correlated electron systems is not only of fundamental interest, but it is also essential for designing new functional materials. Although the results of approximate cluster techniques have to be handled with care, the new aspects presented in this thesis pave the way for further investigations in this direction.
In view of quite recent experimental activities on magnetic and superconducting properties of honeycomb and hexagonal lattice based materials, in this thesis we have used field-theoretic and many-body methods to investigate magnetic and superconducting properties of the large-U Hubbard model on the honeycomb lattice at half-filling and in the hole-doped regime. Within the framework of a functional-integral approach, we obtain the Lagrangian density associated with the charge (Grassmann fields) and spin [SU(2) gauge fields] degrees of freedom. The Hamiltonian related to the charge degrees of freedom is exactly diagonalized. In the strong-coupling regime, we derive a perturbative low-energy theory suitable to describe the (quantum) magnetic and superconducting phases at half-filling and in the hole-doped regime. At half-filling, we deal with the underlying spin degrees of freedom of the quantum antiferromagnetic (AF) Heisenberg model by employing a second-order spin-wave analysis, in which case we have calculated the ground-state energy and the staggered magnetization; the results are in very good agreement with previous studies. Further, in the continuum, we derive a nonlinear σ-model with a topological Hopf term that describes the AF-VBS (valence bond solid) competition. In the challenging hole-doped regime, our approach allows the derivation of a t-J Hamiltonian, and the analysis of the role played by charge and spin quantum fluctuations on the ground-state energy and, particularly, on the breakdown of the AF order at a critical hole doping; the results are benchmarked against recent Grassmann tensor product state simulations. In addition, we have performed an extensive study of the electronic structure of the doped system for each competing phase: AF, ferromagnetic (FM), and (spin-singlet pairing) s-, dx₂₋ʏ₂ – and idxʏ -wave superconducting (SC) state induced by purely electronic effects. In this context, an energetic analysis of the ground state of these phases reveal that the AF order prevails for low hole doping, while a dominantly chiral dx₂₋ʏ₂ + idxʏ superconducting state was found in the vicinity of the Van Hove singularity (high hole doping). We also stress that a thermodynamic analysis of the superconducting phase shows that the critical temperature is directly related to the exchange constant J = 4t²/U, in which t denotes the hopping amplitude and U the on-site Coulomb repulsion of the Hubbard model (purely electronic origin). Remarkably, the competition between the AF and dx₂₋ʏ₂ + idxʏ SC phases takes place by the occurrence of a first-order transition accompanied by a spatial phase separation of the referred phases.
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 following thesis is divided in two main parts. Chapters 2, 3 and 4 are devoted to the study of the so called pair-density-wave (PDW) superconducting state and some of its connections to electronic liquid crystal (ELC) phases, its topological aspects in a one dimensional model and its appearance in a quasi-one dimensional system. On the other hand, chapter 5 is focused on the investigation of the classical statistical mechanics properties of dimers, in particular, the dimer model on the Aztec diamond graph and its relation with the octahedron equation. In chapter 2 we present a theory of superconducting states where the Cooper pairs have a nonzero center-of-mass momentum, inhomogeneous superconducting states known as a pair-density-waves (PDWs) states. We show that in a system of spin-1/2 fermions in two dimensions in an electronic nematic spin-triplet phase where rotational symmetry is broken in both real and spin space PDW phases arise naturally in a theory that can be analysed using controlled approximations. We show that several superfluid phases that may arise in this phase can be treated within a controlled BCS mean field theory, with the strength of the spin-triplet nematic order parameter playing the role of the small parameter of this theory. We find that in a spin-triplet nematic phase, in addition to a triplet p-wave and spin-singlet d-wave (or s depending on the nematic phase) uniform superconducting states, it is also possible to have a d-wave (or s) PDW superconductor. The PDW phases found here can be either unidirectional, bidirectional, or tridirectional depending on the spin-triplet nematic phase and which superconducting channel is dominant. In addition, a triple-helix state is found in a particular channel. We show that these PDW phases are present in the weak-coupling limit, in contrast to the usual Fulde-Ferrell-Larkin-Ovchinnikov phases, which require strong coupling physics in addition to a large magnetic field (and often both). In chapter 3 we show that the pair-density-wave (PDW) superconducting state emergent in extended Heisenberg-Hubbard models in two-leg ladders is topological in the presence of an Ising spin symmetry and supports a Majorana zero mode (MZM) at an open boundary and at a junction with a uniform d-wave one-dimensional superconductor. Similarly to a conventional finite-momentum paired state, the order parameter of the PDW state is a charge-2e field with finite momentum. However, the order parameter here is a quartic electron operator and conventional mean-field theory cannot be applied to study this state. We use bosonization to show that the 1D PDW state has a MZM at a boundary. This superconducting state is an exotic topological phase supporting Majorana fermions with finite-momentum pairing fields and charge-4e superconductivity. In chapter 4 we provide a quasi-one-dimensional model which can support a PDW state. The model consists of an array of strongly-interacting one-dimensional systems, where the one-dimensional systems are coupled to each other by local interactions.Within the interchain mean-field theory (MFT), we find several SC states from the model, including a conventional uniform SC state, PDW SC state, and a coexisting phase of the uniform SC and PDW states. In this quasi-1D regime we can treat the strong correlation physics essentially exactly using bosonization methods and the crossover to the 2D system by means of interchain MFT. The resulting critical temperatures of the SC phases generically exhibit a power-law scaling with the coupling constants of the array, instead of the essential singularity found in weak-coupling BCS-type theories. Electronic excitations with an open Fermi surface, which emerge from the electronic Luttinger liquid systems below their crossover temperature to the Fermi liquid, are then coupled to the SC order parameters via the proximity effect. From the Fermi surface thus coupled to the SC order parameters, we calculate the quasiparticle spectrum in detail. We show that the quasiparticle spectrum can be fully gapped or nodal in the uniform SC phase and also in the coexisting phase of the uniform SC and PDW parameters. In the pure PDW state, the excitation spectrum has a reconstructed Fermi surface in the form of Fermi pockets of Bogoliubov quasiparticles. In chapter 5 we study the octahedron relation (also known as the A∞ T-system), obeyed in particular by the partition function for dimer coverings of the Aztec Diamond graph. For a suitable class of doubly periodic initial conditions, we find exact solutions with a particularly simple factorized form. For these, we show that the density function that measures the average dimer occupation of a face of the Aztec graph, obeys a system of linear recursion relations with periodic coefficients. This allows us to explore the thermodynamic limit of the corresponding dimer models and to derive exact “arctic” curves separating the various phases of the system.
This thesis investigates the interplay between electron-electron and electron-phonon interactions in low-dimensional lattice models, specifically within the Hubbard model incorporating Su-Schrieffer-Heeger (SSH) coupling with phonons. SSH interactions, which modulate the electron hopping amplitude through lattice displacements, have recently attracted renewed interest due to their potential relevance in strongly correlated materials. The main objective of this work is to determine whether the inclusion of SSH couplings in the Hubbard model can stabilize superconducting or spin-gapped metallic states, as well as other exotic phases of matter. While the half-filled SSH-Hubbard model has been extensively studied in one dimension and only partially examined in two dimensions, its phase diagram in the doped regime remains largely unexplored. This thesis aims to fill this gap, focusing on the emergence of correlated phases away from half-filling in both one and two spatial dimensions. We use the variational Monte Carlo method with a trial wave function defined in a hybrid electron-phonon Hilbert space. In addition to Jastrow factors that account for correlation effects, the wave function includes an antisymmetric component that captures backflow-inspired electron-phonon correlations. This approach enables us to accurately model both electron-electron and electron-phonon interactions. Moreover, since the variational method avoids the sign problem, we can simulate the system in the doped regime, an area typically inaccessible to exact quantum Monte Carlo methods (QMC), which are limited to half-filling. In one dimension, we observe the stabilization of a Luther-Emery liquid, a metallic phase with gapless charge excitations and a finite spin gap, upon doping a Peierls insulating phase. This finding resolves a long-standing debate regarding its realization in SSH chains. Benchmarking with density-matrix renormalization group results confirms the reliability of our approach. In two dimensions, we investigate the SSH-Hubbard model on the square lattice, obtaining results that are consistent with previous exact QMC studies at half-filling. Upon doping, we observe that SSH phonons can induce s-wave superconductivity in the weakly correlated regime. However, we suggest that the enhancement of d-wave superconductivity by SSH coupling in the strongly correlated regime may significantly depend on the specific implementation of the SSH interaction. In one particular case, we find that d-wave superconductivity is slightly suppressed by SSH coupling. These results are further confirmed by a Hartree-Fock analysis performed in the anti-adiabatic limit. Overall, this thesis highlights the crucial role of lattice effects in correlated electron systems, shedding light on the distinct contributions of SSH-type couplings in one and two dimensions. Our findings provide new insights into the microscopic mechanisms behind exotic metallic phases and set the stage for further studies of phonon-mediated phenomena in correlated systems.