Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Extensions of Numerical Methods for Strongly Correlated Electron Systems

Karlis Mikelsons

Published Jan 1, 2009• 4 citations
OpenAlex

Abstract

This work presents extensions of the numerical methods for strongly correlated electron systems.The first part of the thesis discusses extensions and applications of the quantum cluster theories to the systems of classical spins.It is shown that such extensions can provide faster convergence through better estimation of the effects of fluctuations, yet they can also possess shortcomings which limit their application in the studies of the phase transitions.The second part of the thesis is dedicated to the numerical studies of the Hubbard model.Present Quantum Monte Carlo methods are reviewed and relationships among them are elucidated.The final part of the thesis contains the application of the developed numerical methods to investigate the phase diagram of the two-dimensional Hubbard model, especially the evidence of the Quantum Critical Point (QCP) at a finite doping.High accuracy results for thermodynamic quantities are presented in support of the existence of the QCP at a finite doping in two-dimensional Hubbard model.The relation of the QCP to the charge fluctuations is revealed and a mechanism that relates QCP to incipient phase separation is proposed.I am thankful to the University of Cincinnati for giving me the opportunity to pursue my studies, and for providing all the resources and support for my education and

Related papers

Numerical Methods for Two-particle Fluctuations and Real-time Dynamics of Strongly Correlated Electron Systems

Xinyang Dong

This thesis contains a series of numerical studies of strongly correlated electron systems. In these systems, interesting emergent properties are brought about by strong electron correlations, leading to unusual fluctuations and phase transitions. We start with a review of basic concepts in many-body physics from a field theory point of view, including second quantization and path integral formalisms. We then introduce extensions of the single-particle and two-particle Green’s function formalisms in the singlet superconducting state, which enable us to extract physical information of the systems in the symmetry-broken phase. Next, we proceed to the models and numerical methods we use to study the strongly correlated systems. We introduce a low energy effective model -- the Hubbard model, which contains a subset of the electron degrees of freedom and can be solved numerically using advanced many-body methods. With the dynamical cluster approximation and continuous-time quantum Monte Carlo impurity solver, we are able to study the competing fluctuations in the paramagnetic state and analyze the fluctuations behind superconductivity in the singlet superconducting state. We introduce self-consistent diagrammatic methods used in realistic material calculations in the last part of this thesis, with an outline of group theory concepts that can be used to optimize the simulations. Numerical representations and methods for effectively solving equations in realistic material calculations are discussed afterwards. We briefly review some of the developments in this field, and then introduce spectral methods that are based on mathematical properties of Legendre polynomials for solving both the imaginary- and real-time Dyson equations.

•Jan 1, 2022•0 citations

Criticality and Superconductivity in the Two-dimensional Hubbard Model of Strongly Correlated Electronic Systems

Ehsan Khatami

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

•Jan 1, 2009•0 citations

Local and Nonlocal Correlations in Interacting Electron Systems

Josef Kaufmann

Strong electronic correlations are the root of many fascinating phenomena of condensed matter physics. However, they also make calculations notoriously difficult. This thesis is concerned with the development and improvement of computational methods for strongly correlated electron systems. After a brief recapitulation of the Green’s function formalism,which is the theoretical basis of most such methods, we discuss the analytical continuation from imaginary to real frequencies. Especially the maximum entropy method is explained in detail and an implementation there of is presented.The remainder of this thesis deals with methods to compute Green’s functions on thei maginary frequency axis. First, the Anderson impurity model (AIM) is considered. It can be solved numerically – a task for which several methods exist. We focus on continuous-time quantum Monte Carlo with worm sampling, through which various kinds of correlation functions can be calculated. On this basis we derive and implement symmetric improved estimators that allow for the computation of the self-energy and vertex functions where the level of stochastic noise is largely reduced with respect to conventional estimators.As a next step the dynamical mean-field theory (DMFT) is introduced as a method to capture local correlations in the Hubbard model (HM) by self-consistently mapping it onto the AIM. We study the convergence behavior of this iterative procedure and the improvement of convergence due to mixing. Then, we compare electron scattering rates from DMFT in the weak-coupling regime of the square-lattice HM to results from Boltzmann scattering theory.While local correlations account for many interesting phenomena, such as the Mott-Hubbard metal-insulator transition, they are not sufficient to correctly describe, e.g., the physics of phase transitions in two- or three-dimensional systems. It is thus the final step in this thesis to also include nonlocal correlations by means of the dynamical vertex approximation (DΓA). For this we develop a new self-consistency scheme, which is then applied to several systems: We study the square-lattice HM, where we compare to benchmark methods in the weak-coupling regime, but also present results for intermediate and strong coupling. Furthermore, a two-band model with doping and a three-band Wannier model for strontium vanadate are studied. As an application of our new self-consistent DΓA, we study the effects of geometrical frustration in the kagome lattice, where we also compare our results to determinant quantum Monte Carlo.

•May 1, 2021•1 citation

Quantum Monte Carlo methods and strongly correlated electrons on honeycomb structures

Thomas C. Lang

In this thesis we apply recently developed, as well as sophisticated quantum Monte Carlo methods to numerically investigate models of strongly correlated electron systems on honeycomb structures. The latter are of particular interest owing to their unique properties when simulating electrons on them, like the relativistic dispersion, strong quantum fluctuations and their resistance against instabilities. This work covers several projects including the advancement of the weak-coupling continuous time quantum Monte Carlo and its application to zero temperature and phonons, quantum phase transitions of valence bond solids in spin-1/2 Heisenberg systems using projector quantum Monte Carlo in the valence bond basis, and the magnetic field induced transition to a canted antiferromagnet of the Hubbard model on the honeycomb lattice. The emphasis lies on two projects investigating the phase diagram of the SU(2) and the SU(N)-symmetric Hubbard model on the hexagonal lattice. At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. Previously elusive in experimentally relevant microscopic two-dimensional models, we show by means of large-scale quantum Monte Carlo simulations of the SU(2) Hubbard model on the honeycomb lattice, that a quantum spin-liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence bond liquid, akin to the one proposed for high temperature superconductors. Inspired by the rich phase diagrams of SU(N) models we study the SU(N)-symmetric Hubbard Heisenberg quantum antiferromagnet on the honeycomb lattice to investigate the reliability of 1/N corrections to large-N results by means of numerically exact QMC simulations. We study the melting of phases as correlations increase with decreasing N and determine whether the quantum spin liquid found in the SU(2) Hubbard model at intermediate coupling is a specific feature, or also exists in the unconstrained t-J model and higher symmetries.

•Jan 1, 2010•2 citations

Two-particle level diagrammatic approaches for strongly correlated systems

Herbert F. Fotso

This thesis presents the development of new numerical methods for the treatment of strongly correlated electron systems based on self-consistent approaches at both the one and the two-particle level such as the parquet formalism. The parquet formalism was solved for the first time on a two-dimensional cluster. When the fully irreducible vertex is approximated by the bare vertex, we obtain the parquet approximation. Its validity was investigated by comparing results that it produces to those of other conserving approximations such as the FLuctuation EXchange (FLEX) approximation or the Second Order Perturbation Theory (SOPT). We found that it provides a significant improvement of FLEX or SOPT and a satisfactory agreement with Quantum Monte Carlo results despite instabilities in the self-consistency at low temperatures and for strong Coulomb interaction. We use the parquet formalism to study the Quantum Critical Point at finite doping in the Hubbard model by decomposing the vertex into its contributions from different channels. We apply this decomposition to the pairing channel and we find that the dominant contribution to the vertex originates in the spin channel even at the quantum critical doping. Furthermore, we explore the divergence of the two parts of the pairing matrix at optimal doping and observe that the irreducible vertex decreases monotonically as the doping is increased while the bare susceptibility exhibits an algebraic divergence at the quantum critical doping supporting the Quantum Critical BCS scenario proposed by She and Zaanen. To circumvent the instabilities in the iteration of the parquet formalism, we explored the dual fermion approach introduced by Rubtsov et al. Here, we extended the formalism to the Dynamical Cluster Approximation, in the process introducing a small parameter in the dual fermions perturbation theory. We demonstrate the quality of the resulting Dual Fermion DCA through a systematic study of the cluster size dependence and of the different perturbative approximations. These efforts represent the initial steps in the development of the Multi-Scale Many Body approach that appropriately treats correlations at different length scales.

Civil War Book Review•Jul 1, 2011•1 citation

Numerical and analytical approaches to strongly correlated electron systems

C. Knecht

In this thesis we consider three different models for strongly correlated electrons, namely a multi-band Hubbard model as well as the spinless Falicov-Kimball model, both with a semi-elliptical density of states in the limit of infinite dimensions d, and the attractive Hubbard model on a square lattice in d=2.\n\nIn the first part, we study a two-band Hubbard model with unequal bandwidths and anisotropic Hund's rule coupling (J_z-model) in the limit of infinite dimensions within the dynamical mean-field theory (DMFT). Here, the DMFT impurity problem is solved with the use of quantum Monte Carlo (QMC) simulations. Our main result is that the J_z-model describes the occurrence of an orbital-selective Mott transition (OSMT), in contrast to earlier findings. We investigate the model with a high-precision DMFT \nalgorithm, which was developed as part of this thesis and which supplements QMC with a high-frequency expansion of the self-energy.\nThe main advantage of this scheme is the extraordinary accuracy of the numerical solutions, which can be obtained already with moderate computational effort, so that studies of multi-orbital systems within the DMFT+QMC are strongly improved. We also found that a suitably defined\nFalicov-Kimball (FK) model exhibits an OSMT, revealing the close connection of the Falicov-Kimball physics to the J_z-model in the OSM phase.\n\nIn the second part of this thesis we study the attractive Hubbard model in two spatial dimensions within second-order self-consistent perturbation theory.\nThis model is considered on a square lattice at finite doping and at low temperatures. Our main result is that the predictions of first-order perturbation theory (Hartree-Fock approximation) are renormalized by a factor of the order of unity even at arbitrarily weak interaction (U->0). The \nrenormalization factor q can be evaluated as a function of the filling n for 0<n<1. In the limit n->0, the q-factor vanishes, signaling the divergence of self-consistent perturbation theory in this limit. Thus we present the first asymptotically exact results at weak-coupling for the negative-U Hubbard model in d=2 at finite doping.

•Jan 1, 2006•1 citation
This work presents extensions of the numerical methods for strongly correlated electron systems.The first part of the thesis discusses extensions and applications of the quantum cluster theories to the systems of classical spins.It is shown that such extensions can provide faster convergence through better estimation of the effects of fluctuations, yet they can also possess shortcomings which limit their application in the studies of the phase transitions.The second part of the thesis is dedicated to the numerical studies of the Hubbard model.Present Quantum Monte Carlo methods are reviewed and relationships among them are elucidated.The final part of the thesis contains the application of the developed numerical methods to investigate the phase diagram of the two-dimensional Hubbard model, especially the evidence of the Quantum Critical Point (QCP) at a finite doping.High accuracy results for thermodynamic quantities are presented in support of the existence of the QCP at a finite doping in two-dimensional Hubbard model.The relation of the QCP to the charge fluctuations is revealed and a mechanism that relates QCP to incipient phase separation is proposed.I am thankful to the University of Cincinnati for giving me the opportunity to pursue my studies, and for providing all the resources and support for my education and
This thesis contains a series of numerical studies of strongly correlated electron systems. In these systems, interesting emergent properties are brought about by strong electron correlations, leading to unusual fluctuations and phase transitions. We start with a review of basic concepts in many-body physics from a field theory point of view, including second quantization and path integral formalisms. We then introduce extensions of the single-particle and two-particle Green’s function formalisms in the singlet superconducting state, which enable us to extract physical information of the systems in the symmetry-broken phase. Next, we proceed to the models and numerical methods we use to study the strongly correlated systems. We introduce a low energy effective model -- the Hubbard model, which contains a subset of the electron degrees of freedom and can be solved numerically using advanced many-body methods. With the dynamical cluster approximation and continuous-time quantum Monte Carlo impurity solver, we are able to study the competing fluctuations in the paramagnetic state and analyze the fluctuations behind superconductivity in the singlet superconducting state. We introduce self-consistent diagrammatic methods used in realistic material calculations in the last part of this thesis, with an outline of group theory concepts that can be used to optimize the simulations. Numerical representations and methods for effectively solving equations in realistic material calculations are discussed afterwards. We briefly review some of the developments in this field, and then introduce spectral methods that are based on mathematical properties of Legendre polynomials for solving both the imaginary- and real-time Dyson equations.
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
Strong electronic correlations are the root of many fascinating phenomena of condensed matter physics. However, they also make calculations notoriously difficult. This thesis is concerned with the development and improvement of computational methods for strongly correlated electron systems. After a brief recapitulation of the Green’s function formalism,which is the theoretical basis of most such methods, we discuss the analytical continuation from imaginary to real frequencies. Especially the maximum entropy method is explained in detail and an implementation there of is presented.The remainder of this thesis deals with methods to compute Green’s functions on thei maginary frequency axis. First, the Anderson impurity model (AIM) is considered. It can be solved numerically – a task for which several methods exist. We focus on continuous-time quantum Monte Carlo with worm sampling, through which various kinds of correlation functions can be calculated. On this basis we derive and implement symmetric improved estimators that allow for the computation of the self-energy and vertex functions where the level of stochastic noise is largely reduced with respect to conventional estimators.As a next step the dynamical mean-field theory (DMFT) is introduced as a method to capture local correlations in the Hubbard model (HM) by self-consistently mapping it onto the AIM. We study the convergence behavior of this iterative procedure and the improvement of convergence due to mixing. Then, we compare electron scattering rates from DMFT in the weak-coupling regime of the square-lattice HM to results from Boltzmann scattering theory.While local correlations account for many interesting phenomena, such as the Mott-Hubbard metal-insulator transition, they are not sufficient to correctly describe, e.g., the physics of phase transitions in two- or three-dimensional systems. It is thus the final step in this thesis to also include nonlocal correlations by means of the dynamical vertex approximation (DΓA). For this we develop a new self-consistency scheme, which is then applied to several systems: We study the square-lattice HM, where we compare to benchmark methods in the weak-coupling regime, but also present results for intermediate and strong coupling. Furthermore, a two-band model with doping and a three-band Wannier model for strontium vanadate are studied. As an application of our new self-consistent DΓA, we study the effects of geometrical frustration in the kagome lattice, where we also compare our results to determinant quantum Monte Carlo.
In this thesis we apply recently developed, as well as sophisticated quantum Monte Carlo methods to numerically investigate models of strongly correlated electron systems on honeycomb structures. The latter are of particular interest owing to their unique properties when simulating electrons on them, like the relativistic dispersion, strong quantum fluctuations and their resistance against instabilities. This work covers several projects including the advancement of the weak-coupling continuous time quantum Monte Carlo and its application to zero temperature and phonons, quantum phase transitions of valence bond solids in spin-1/2 Heisenberg systems using projector quantum Monte Carlo in the valence bond basis, and the magnetic field induced transition to a canted antiferromagnet of the Hubbard model on the honeycomb lattice. The emphasis lies on two projects investigating the phase diagram of the SU(2) and the SU(N)-symmetric Hubbard model on the hexagonal lattice. At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. Previously elusive in experimentally relevant microscopic two-dimensional models, we show by means of large-scale quantum Monte Carlo simulations of the SU(2) Hubbard model on the honeycomb lattice, that a quantum spin-liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence bond liquid, akin to the one proposed for high temperature superconductors. Inspired by the rich phase diagrams of SU(N) models we study the SU(N)-symmetric Hubbard Heisenberg quantum antiferromagnet on the honeycomb lattice to investigate the reliability of 1/N corrections to large-N results by means of numerically exact QMC simulations. We study the melting of phases as correlations increase with decreasing N and determine whether the quantum spin liquid found in the SU(2) Hubbard model at intermediate coupling is a specific feature, or also exists in the unconstrained t-J model and higher symmetries.
This thesis presents the development of new numerical methods for the treatment of strongly correlated electron systems based on self-consistent approaches at both the one and the two-particle level such as the parquet formalism. The parquet formalism was solved for the first time on a two-dimensional cluster. When the fully irreducible vertex is approximated by the bare vertex, we obtain the parquet approximation. Its validity was investigated by comparing results that it produces to those of other conserving approximations such as the FLuctuation EXchange (FLEX) approximation or the Second Order Perturbation Theory (SOPT). We found that it provides a significant improvement of FLEX or SOPT and a satisfactory agreement with Quantum Monte Carlo results despite instabilities in the self-consistency at low temperatures and for strong Coulomb interaction. We use the parquet formalism to study the Quantum Critical Point at finite doping in the Hubbard model by decomposing the vertex into its contributions from different channels. We apply this decomposition to the pairing channel and we find that the dominant contribution to the vertex originates in the spin channel even at the quantum critical doping. Furthermore, we explore the divergence of the two parts of the pairing matrix at optimal doping and observe that the irreducible vertex decreases monotonically as the doping is increased while the bare susceptibility exhibits an algebraic divergence at the quantum critical doping supporting the Quantum Critical BCS scenario proposed by She and Zaanen. To circumvent the instabilities in the iteration of the parquet formalism, we explored the dual fermion approach introduced by Rubtsov et al. Here, we extended the formalism to the Dynamical Cluster Approximation, in the process introducing a small parameter in the dual fermions perturbation theory. We demonstrate the quality of the resulting Dual Fermion DCA through a systematic study of the cluster size dependence and of the different perturbative approximations. These efforts represent the initial steps in the development of the Multi-Scale Many Body approach that appropriately treats correlations at different length scales.
In this thesis we consider three different models for strongly correlated electrons, namely a multi-band Hubbard model as well as the spinless Falicov-Kimball model, both with a semi-elliptical density of states in the limit of infinite dimensions d, and the attractive Hubbard model on a square lattice in d=2.\n\nIn the first part, we study a two-band Hubbard model with unequal bandwidths and anisotropic Hund's rule coupling (J_z-model) in the limit of infinite dimensions within the dynamical mean-field theory (DMFT). Here, the DMFT impurity problem is solved with the use of quantum Monte Carlo (QMC) simulations. Our main result is that the J_z-model describes the occurrence of an orbital-selective Mott transition (OSMT), in contrast to earlier findings. We investigate the model with a high-precision DMFT \nalgorithm, which was developed as part of this thesis and which supplements QMC with a high-frequency expansion of the self-energy.\nThe main advantage of this scheme is the extraordinary accuracy of the numerical solutions, which can be obtained already with moderate computational effort, so that studies of multi-orbital systems within the DMFT+QMC are strongly improved. We also found that a suitably defined\nFalicov-Kimball (FK) model exhibits an OSMT, revealing the close connection of the Falicov-Kimball physics to the J_z-model in the OSM phase.\n\nIn the second part of this thesis we study the attractive Hubbard model in two spatial dimensions within second-order self-consistent perturbation theory.\nThis model is considered on a square lattice at finite doping and at low temperatures. Our main result is that the predictions of first-order perturbation theory (Hartree-Fock approximation) are renormalized by a factor of the order of unity even at arbitrarily weak interaction (U->0). The \nrenormalization factor q can be evaluated as a function of the filling n for 0<n<1. In the limit n->0, the q-factor vanishes, signaling the divergence of self-consistent perturbation theory in this limit. Thus we present the first asymptotically exact results at weak-coupling for the negative-U Hubbard model in d=2 at finite doping.