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.
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.
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.
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.
Using the variational cluster approach, we study the transition from the antiferromagnetic to the superconducting phase of the two-dimensional Hubbard model at zero temperature. Our calculations are based on a method to evaluate the VCA grand potential which employs a modified Lanczos algorithm and avoids integrations over the real or imaginary frequency axis. Thereby, very accurate results are possible for cluster sizes not accessible to full diagonalization. This is important for an improved treatment of short-range correlations, including correlations between Cooper pairs in particular. We apply this improved method in order to investigate the cluster-size dependence of the phase-separation tendency that has been proposed recently on the basis of calculations for smaller clusters. While the energy barrier associated with phase separation rapidly decreases with increasing cluster size for both hole and electron doping, the extension of the phase-separation region behaves differently in the two cases. More specifically, our results suggest that phase separation remains persistent in the hole-doped case and disappears in the electron-doped case. We also study the evolution of the single-particle spectrum as a function of doping and point out the relevance of our results for experimental findings in electron and hole-doped materials.
Nematicity is a well-known property of liquid crystals and has been recently discussed in the context of strongly interacting electrons. An electronic nematic phase has been seen in many experiments in certain strongly correlated materials, in particular, in the pseudogap phase generic to many hole-doped cuprate superconductors. Recent measurements in high Tc superconductors have shown that even if the lattice is perfectly rotationally symmetric, the ground state can still have strongly nematic local properties. Our study of the two-dimensional one-band Hubbard model provides strong support for the recent experimental results on local rotational C4 symmetry breaking. The variational cluster approach is used here to show the possibility of an electronic nematic state and the proximity of the underlying symmetry-breaking ground state within the Hubbard model. We identify this nematic phase in the overdoped region and show that the local nematicity decreases with increasing electron filling. Our results also indicate that strong Coulomb interaction may drive the nematic phase into a phase similar to the stripe structure. The calculated spin (magnetic) correlation function in momentum space shows the effects resulting from real-space nematicity.
Many properties of the two-dimensional Hubbard model have been explored for the model parameters appropriate for strongly correlated electronic systems, especially cuprate superconductors.Most of the calculations are done using a well-established dynamical cluster quantum Monte Carlo method.Using this method, we investigate the effect of long-range hoppings on superconductivity with and without the presence of phonons on small clusters.The superconducting transition temperature, T c , is found to generally decrease with a negative next nearest neighbor hopping, t ′ .In the presence of the Holstein phonons, a finite t ′ enhances T c in the under-doped region for the hole-doped system, consistent with band structure calculations and experiment.The validity of the spin-susceptibility-mediated pairing in this model is studied and found to yield symmetries other than d-wave when a finite t ′ is considered.A new numerical algorithm for solving the embedded cluster problems is introduced, and used to calculate the thermodynamic properties of the model on larger clusters, especially those associated with the quantum critical behavior at finite doping.Our results suggest that the quantum critical point (QCP) which separates the Fermi liquid and pseudogap regions, is the second order terminus of the line of first order phase separation transition in the limit when a positive t ′ goes to zero.For small t ′ , the superconducting dome is centered at the QCP, suggesting that charge fluctuations might have a role in the pairing mechanism in this model.I am very thankful to University of Cincinnati for providing me with the opportunity to pursue my education and the support I received during my program.Part of this thesis is a result of my work as a researcher at Louisiana State University.I am also thankful to them for their support.I would like to thank my adviser, Mark Jarrell, who taught me how to be an effective researcher, provided me with the opportunities to interact with other professionals and fully supported me financially for most of the duration of the
We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half-filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalization of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions $U\ensuremath{\geqslant}3t$. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.