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

arXiv:2209.08092 · Physical review. B./Physical review. B

Fermi arcs from dynamical variational Monte Carlo

Peter Rosenberg, David Sénéchal, A.–M. S. Tremblay, M. Charlebois

Published Dec 19, 2022• 17 citations
arXivHTMLPDF

Abstract

Variational Monte Carlo is a many-body numerical method that scales well with system size. It has been extended to study the Green function only recently by Charlebois and Imada [Phys. Rev. X 10, 041023 (2020)]. Here we generalize the approach to systems with open boundary conditions in the absence of translational invariance. Removing these constraints permits the application of embedding techniques like cluster perturbation theory (CPT). Using this approach we provide robust evidence of the existence of Fermi arcs in the one-band Hubbard model, an enduring problem in the physics of the pseudogap in cuprate high-temperature superconductors. We study the behavior of the Fermi surface and of the density of states as a function of hole doping for clusters of up to 64 sites, well beyond the reach of modern exact diagonalization solvers. We observe that the technique reliably captures the transition from a Mott insulator at half filling to a pseudogap, evidenced by the formation of Fermi arcs, and finally to a metallic state at large doping. The ability to treat large clusters with quantum cluster methods helps to minimize potential finite-size effects and enables the study of systems with long-range orders, which will help extend the reach of these already powerful methods and provide important insights on the nature of various strongly correlated many-electron systems, including the high-${\mathrm{T}}_{c}$ cuprate superconductors.

DOI
10.1103/physrevb.106.245132

Related papers

Dynamical variational Monte Carlo as a quantum impurity solver: Application to cluster dynamical mean field theory

Peter Rosenberg, David Sénéchal, A.–M. S. Tremblay +1 more

Two of the primary sources of error in the Cluster dynamical mean field theory (CDMFT) technique arise from the use of finite size clusters and finite size baths, which makes the development of impurity solvers that can treat larger systems an essential goal. In this work we introduce an impurity solver based on the recently developed dynamical variational Monte Carlo (dVMC) method. Variational Monte Carlo possesses a favorable scaling as a function of system size, which enables the treatment of systems beyond the reach of current exact diagonalization solvers. To benchmark the technique, we perform a systematic set of CDMFT calculations on the single-band one-dimensional Hubbard model. We compare to results obtained with an exact diagonalization solver for small clusters, and against the exact solution in the thermodynamic limit obtained by Lieb and Wu for larger clusters. The development of improved impurity solvers will help extend the reach of quantum cluster methods, which can be applied to a wide range of strongly correlated electron systems, promising new insights into their emergent behavior.

arXiv:2307.15738•Physical review. B./Physical review. B•Dec 8, 2023•4 citations

Effect of hole doping on the electronic structure and the Fermi surface in the Hubbard model within norm-conserving cluster pertubation theory

Sergey V. Nikolaev, Sergei G. Ovchinnikov

The concentration dependences of the band structure, spectral weight, density of states, and Fermi surface in the paramagnetic state are studied in the Hubbard model within cluster pertubation theory with 2 × 2 clusters. Representation of the Hubbard X operators makes it possible to control conservation of the spectral weight in constructing cluster perturbation theory. The calculated value of the ground-state energy is in good agreement with the results obtained using nonperturbative methods such as the quantum Monte Carlo method, exact diagonalization of a 4 × 4 cluster, and the variational Monte Carlo method. It is shown that in the case of hole doping, the states in the band gap (in-gap states) lie near the top of the lower Hubbard band for large values of U and near the bottom of the upper band for small U . The concentration dependence of the Fermi surface strongly depends on hopping to second ( t ′) and third ( t ″) neighbors. For parameter values typical of HTSC cuprates, the existence of three concentration regions with different Fermi surfaces is demonstrated. It is shown that broadening of the spectral electron density with an energy resolution typical of contemporary ARPES leads to a pattern of arcs with a length depending on the concentration. Only an order-of-magnitude decrease in the linewidth makes it possible to obtain the true Fermi surface from the spectral density. The kinks associated with strong electron correlations are detected in the dispersion relation below the Fermi level.

Journal of Experimental and Theoretical Physics•Jan 1, 2012•12 citations

Fermi arcs and pseudogap in a lattice model of a doped orthogonal metal

Chuang Chen, Tian Yuan, Yang Qi +1 more

Since the discovery of the pseudogap and Fermi arc states in underdoped cuprates, the understanding of such non-Fermi-liquid states and the associated violation of Luttinger's theorem have been the central theme in correlated electron systems. However, still lacking is a well-accepted theoretical framework to unambiguously explain these metallic states that are clearly beyond Landau's Fermi liquid and Luttinger's theorem of a Fermi surface and electron filling. Here, we design a lattice model of orthogonal metals with fermion and Ising matter fields coupled to topological order and, by solving the model via unbiased quantum Monte Carlo simulation at generic electron fillings, find that the system gives birth to phenomena of the Fermi arc and pseudogap in the single-particle spectrum that go beyond the Luttinger sum rule with broken Fermi surface but no symmetry breaking. The pseudogap and Fermi arcs coexist with a background of a deconfined ${Z}_{2}$ gauge field, and we further find that the confinement transition of the gauge field triggers a superconductivity instability and that the hopping of the gauge-neutral fermions brings the ``large'' Fermi surface back from the Fermi arc state. Our unbiased numerical results provide a concrete model realization and theoretical framework for the coupling between gauge field and fermions and, in the process, generate the rich phenomena of the pseudogap, the Fermi arc, and superconductivity in generic correlated electron systems.

arXiv:2007.05543•Physical review. B./Physical review. B•Apr 22, 2021•13 citations

Nonperturbative approach to full Mott behavior

Tudor D. Stanescu, Philip Phillips

Though most fermionic Mott insulators order at low temperatures, ordering is ancillary to their insulating behavior. Our emphasis here is on disentangling ordering from the intrinsic strongly correlated physics of a doped half-filled band. To this end, we focus on the two-dimensional Hubbard model. Because the charge gap arises from on-site correlations, we have been refining the nonperturbative approach of Matsumoto and Mancini [Phys. Rev. B 55, 2095 (1997)] which incorporates local physics. Crucial to this method is a self-consistent two-site dynamical cluster expansion which builds in the nearest-neighbor energy scale J. At half-filling, we find that the spectral function possesses a gap of order U and is devoid of any coherent quasiparticle peaks although ordering or charge fractionalization are absent. At low temperatures, local antiferromagnetic correlations emerge. In the doped case, we find that the Fermi surface exceeds the Luttinger volume. The breakdown of Luttinger's theorem in the underdoped regime is traced both to the dynamically generated Mott gap as well to a nonvanishing of the imaginary part of the self-energy at the Fermi level. Spectral weight transfer across the Mott gap also emerges as a ubiquitous feature of a doped Mott insulator and suggests that high- and low-energy scales are inseparable. Additionally in the underdoped regime, we find that a pseudogap exists in the single-particle density of states as well as in the heat capacity. The pseudogap (which is set by the energy scale ${t}^{2}/U)$ is argued to be a ubiquitous feature of a lightly doped Mott state and simply represents the fact that hole transport involves double occupancy. In analogy with the Mott gap and antiferromagnetism, we propose that ordering may also accompany the formation of a pseudogap. We suggest a current pattern within a one-band model that preserves translational but breaks time-reversal symmetry along the canonical x and y axes but not along $x=\ifmmode\pm\else\textpm\fi{}y$ that is consistent with the experimental observations. Finally, we show that the Hall coefficient in a doped Mott insulator must change sign at a doping level $x<1/3.$ The sign change is tied to a termination of strong correlation physics in the doped Mott state.

arXiv:cond-mat/0301254•Physical Review B•Jun 15, 2004•38 citations

Correlated-Electron Systems and High-Temperature Superconductivity

Takashi Yanagisawa, Mitake Miyazaki, Yamaji K.

We present recent theoretical results on superconductivity in correlated-electron systems, especially in the two-dimensional Hubbard model and the three-band d-p model. The mechanism of superconductivity in high-temperature superconductors has been extensively studied on the basis of various electronic models and also electron-phonon models. In this study, we investigate the properties of superconductivity in correlated-electron systems by using numerical methods such as the variational Monte Carlo method and the quantum Monte Carlomethod. The Hubbard model is one of basic models for strongly correlated electron systems, and is regarded as the model of cuprate high temperature superconductors. The d-p model is more realistic model for cuprates. The superconducting condensation energy obtained by adopting the Gutzwiller ansatz is in reasonable agreement with the condensation energy estimated for YBa2Cu3O7. We show the phase diagram of the ground state using this method. We have further investigated the stability of striped and checkerboard states in the under-doped region. Holes doped in a half-filled square lattice lead to an incommensurate spin and charge density wave. The relationship of the hole density x and incommensurability δ, δ~x, is satisfied in the lower doping region, as indicated by the variationalMonte Carlocalculations for the two-dimensional Hubbard model. A checkerboard-like charge-density modulation with a roughly period has also been observed by scanning tunneling microscopy experiments in Bi2212 and Na-CCOC compounds. We have performed a variational Monte Carlo simulation on a two-dimensional t-t′-t″- U Hubbard model with a Bi-2212 type band structure and found that the period checkerboard spin modulation, that is characterized by multi Q vectors, is indeed stabilized. We have further performed an investigation by using a quantumMonte Carlomethod, which is a numerical method that can be used to simulate the behavior of correlated electron systems. We present a new algorithm of the quantum Monte Carlo diagonalization that is a method for the evaluation of expectation value without the negative sign problem. We compute pair correlation functions and show that pair correlation is indeed enhanced with hole doping.

arXiv:1307.2456•Journal of Modern Physics•Jan 1, 2013•23 citations

Charge ordering in extended Hubbard models: Variational cluster approach

Markus Aichhorn, Hans Gerd Evertz, Wolfgang von der Linden +1 more

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.

arXiv:cond-mat/0402580•Physical Review B•Dec 3, 2004•99 citations
Variational Monte Carlo is a many-body numerical method that scales well with system size. It has been extended to study the Green function only recently by Charlebois and Imada [Phys. Rev. X 10, 041023 (2020)]. Here we generalize the approach to systems with open boundary conditions in the absence of translational invariance. Removing these constraints permits the application of embedding techniques like cluster perturbation theory (CPT). Using this approach we provide robust evidence of the existence of Fermi arcs in the one-band Hubbard model, an enduring problem in the physics of the pseudogap in cuprate high-temperature superconductors. We study the behavior of the Fermi surface and of the density of states as a function of hole doping for clusters of up to 64 sites, well beyond the reach of modern exact diagonalization solvers. We observe that the technique reliably captures the transition from a Mott insulator at half filling to a pseudogap, evidenced by the formation of Fermi arcs, and finally to a metallic state at large doping. The ability to treat large clusters with quantum cluster methods helps to minimize potential finite-size effects and enables the study of systems with long-range orders, which will help extend the reach of these already powerful methods and provide important insights on the nature of various strongly correlated many-electron systems, including the high-Tc{\mathrm{T}}_{c}Tc​ cuprate superconductors.
Two of the primary sources of error in the Cluster dynamical mean field theory (CDMFT) technique arise from the use of finite size clusters and finite size baths, which makes the development of impurity solvers that can treat larger systems an essential goal. In this work we introduce an impurity solver based on the recently developed dynamical variational Monte Carlo (dVMC) method. Variational Monte Carlo possesses a favorable scaling as a function of system size, which enables the treatment of systems beyond the reach of current exact diagonalization solvers. To benchmark the technique, we perform a systematic set of CDMFT calculations on the single-band one-dimensional Hubbard model. We compare to results obtained with an exact diagonalization solver for small clusters, and against the exact solution in the thermodynamic limit obtained by Lieb and Wu for larger clusters. The development of improved impurity solvers will help extend the reach of quantum cluster methods, which can be applied to a wide range of strongly correlated electron systems, promising new insights into their emergent behavior.
The concentration dependences of the band structure, spectral weight, density of states, and Fermi surface in the paramagnetic state are studied in the Hubbard model within cluster pertubation theory with 2 × 2 clusters. Representation of the Hubbard X operators makes it possible to control conservation of the spectral weight in constructing cluster perturbation theory. The calculated value of the ground-state energy is in good agreement with the results obtained using nonperturbative methods such as the quantum Monte Carlo method, exact diagonalization of a 4 × 4 cluster, and the variational Monte Carlo method. It is shown that in the case of hole doping, the states in the band gap (in-gap states) lie near the top of the lower Hubbard band for large values of U and near the bottom of the upper band for small U . The concentration dependence of the Fermi surface strongly depends on hopping to second ( t ′) and third ( t ″) neighbors. For parameter values typical of HTSC cuprates, the existence of three concentration regions with different Fermi surfaces is demonstrated. It is shown that broadening of the spectral electron density with an energy resolution typical of contemporary ARPES leads to a pattern of arcs with a length depending on the concentration. Only an order-of-magnitude decrease in the linewidth makes it possible to obtain the true Fermi surface from the spectral density. The kinks associated with strong electron correlations are detected in the dispersion relation below the Fermi level.
Since the discovery of the pseudogap and Fermi arc states in underdoped cuprates, the understanding of such non-Fermi-liquid states and the associated violation of Luttinger's theorem have been the central theme in correlated electron systems. However, still lacking is a well-accepted theoretical framework to unambiguously explain these metallic states that are clearly beyond Landau's Fermi liquid and Luttinger's theorem of a Fermi surface and electron filling. Here, we design a lattice model of orthogonal metals with fermion and Ising matter fields coupled to topological order and, by solving the model via unbiased quantum Monte Carlo simulation at generic electron fillings, find that the system gives birth to phenomena of the Fermi arc and pseudogap in the single-particle spectrum that go beyond the Luttinger sum rule with broken Fermi surface but no symmetry breaking. The pseudogap and Fermi arcs coexist with a background of a deconfined Z2{Z}_{2}Z2​ gauge field, and we further find that the confinement transition of the gauge field triggers a superconductivity instability and that the hopping of the gauge-neutral fermions brings the ``large'' Fermi surface back from the Fermi arc state. Our unbiased numerical results provide a concrete model realization and theoretical framework for the coupling between gauge field and fermions and, in the process, generate the rich phenomena of the pseudogap, the Fermi arc, and superconductivity in generic correlated electron systems.
Though most fermionic Mott insulators order at low temperatures, ordering is ancillary to their insulating behavior. Our emphasis here is on disentangling ordering from the intrinsic strongly correlated physics of a doped half-filled band. To this end, we focus on the two-dimensional Hubbard model. Because the charge gap arises from on-site correlations, we have been refining the nonperturbative approach of Matsumoto and Mancini [Phys. Rev. B 55, 2095 (1997)] which incorporates local physics. Crucial to this method is a self-consistent two-site dynamical cluster expansion which builds in the nearest-neighbor energy scale J. At half-filling, we find that the spectral function possesses a gap of order U and is devoid of any coherent quasiparticle peaks although ordering or charge fractionalization are absent. At low temperatures, local antiferromagnetic correlations emerge. In the doped case, we find that the Fermi surface exceeds the Luttinger volume. The breakdown of Luttinger's theorem in the underdoped regime is traced both to the dynamically generated Mott gap as well to a nonvanishing of the imaginary part of the self-energy at the Fermi level. Spectral weight transfer across the Mott gap also emerges as a ubiquitous feature of a doped Mott insulator and suggests that high- and low-energy scales are inseparable. Additionally in the underdoped regime, we find that a pseudogap exists in the single-particle density of states as well as in the heat capacity. The pseudogap (which is set by the energy scale t2/U){t}^{2}/U)t2/U) is argued to be a ubiquitous feature of a lightly doped Mott state and simply represents the fact that hole transport involves double occupancy. In analogy with the Mott gap and antiferromagnetism, we propose that ordering may also accompany the formation of a pseudogap. We suggest a current pattern within a one-band model that preserves translational but breaks time-reversal symmetry along the canonical x and y axes but not along $x=\ifmmode\pm\else\textpm\fi{}y$ that is consistent with the experimental observations. Finally, we show that the Hall coefficient in a doped Mott insulator must change sign at a doping level x<1/3.x<1/3.x<1/3. The sign change is tied to a termination of strong correlation physics in the doped Mott state.
We present recent theoretical results on superconductivity in correlated-electron systems, especially in the two-dimensional Hubbard model and the three-band d-p model. The mechanism of superconductivity in high-temperature superconductors has been extensively studied on the basis of various electronic models and also electron-phonon models. In this study, we investigate the properties of superconductivity in correlated-electron systems by using numerical methods such as the variational Monte Carlo method and the quantum Monte Carlomethod. The Hubbard model is one of basic models for strongly correlated electron systems, and is regarded as the model of cuprate high temperature superconductors. The d-p model is more realistic model for cuprates. The superconducting condensation energy obtained by adopting the Gutzwiller ansatz is in reasonable agreement with the condensation energy estimated for YBa2Cu3O7. We show the phase diagram of the ground state using this method. We have further investigated the stability of striped and checkerboard states in the under-doped region. Holes doped in a half-filled square lattice lead to an incommensurate spin and charge density wave. The relationship of the hole density x and incommensurability δ, δ~x, is satisfied in the lower doping region, as indicated by the variationalMonte Carlocalculations for the two-dimensional Hubbard model. A checkerboard-like charge-density modulation with a roughly period has also been observed by scanning tunneling microscopy experiments in Bi2212 and Na-CCOC compounds. We have performed a variational Monte Carlo simulation on a two-dimensional t-t′-t″- U Hubbard model with a Bi-2212 type band structure and found that the period checkerboard spin modulation, that is characterized by multi Q vectors, is indeed stabilized. We have further performed an investigation by using a quantumMonte Carlomethod, which is a numerical method that can be used to simulate the behavior of correlated electron systems. We present a new algorithm of the quantum Monte Carlo diagonalization that is a method for the evaluation of expectation value without the negative sign problem. We compute pair correlation functions and show that pair correlation is indeed enhanced with hole doping.
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.