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

Physical review. B, Condensed matter

Static and dynamical properties of doped Hubbard clusters

Elbio R. Dagotto, Adriana Moreo, F. Ortolani, Didier Poilblanc, José A. Riera

Published May 1, 1992• 201 citations
Publisher

Abstract

We study the t-J and the Hubbard models at zero temperature using exact-diagonalization techniques on \ensuremath{\surd}10 \ifmmode\times\else\texttimes\fi{} \ensuremath{\surd}10 and 4\ifmmode\times\else\texttimes\fi{}4 sites clusters. Quantum Monte Carlo simulation results on larger lattices are also presented. All electronic fillings have been analyzed for the three models. We have measured equal-time correlation functions corresponding to various types of order (ranging from ``standard'' staggered spin order to more ``exotic'' possibilities like chiral order), as well as various dynamical properties of these models. Upper bounds for the critical hole doping (${\mathit{x}}_{\mathit{c}}$), where long-range antiferromagnetic order disappears, are presented. It was found that ${\mathit{x}}_{\mathit{c}}$ is very small in agreement with experiments for the high-${\mathit{T}}_{\mathit{c}}$ superconductors. For example, in the t-J model, ${\mathit{x}}_{\mathit{c}}$0.08 at J/t=0.4. However, short-distance spin correlations are important up to much higher dopings producing a sharp well-defined spin-wave-like peak in S(q=(\ensuremath{\pi},\ensuremath{\pi}),\ensuremath{\omega}). Regarding the possibility of phase separation in the Hubbard model, we have studied the behavior of the density of particles, 〈n〉, as a function of the chemical potential, using the Lanczos method on a 4\ifmmode\times\else\texttimes\fi{}4 Hubbard cluster, finding no indications of phase separation for any value of U/t. Then, we conclude that the t-J model at small J/t should not phase separate.In order to compare theoretical predictions with photoemission experiments, we evaluated the electronic density of states N(\ensuremath{\omega}), of the Hubbard and t-J models at several doping fractions. We found that upon doping the antiferromagnetic gap is filled for U\ensuremath{\sim}8t or smaller. The chemical potential moves across the insulating gap as one goes from electron to hole doping of the half-filled cluster, in agreement with x-ray absorption experiments but at variance with photoemission experiments. We have also calculated the optical conductivity, ${\mathrm{\ensuremath{\sigma}}}_{1}$(\ensuremath{\omega}), of the Hubbard and t-J models at all dopings on 4\ifmmode\times\else\texttimes\fi{}4 clusters. Results are compared with experiments and the weight of the Drude peak is presented as a function of couplings and dopings. Spectral weight found at small frequencies is associated with the mid-infrared band observed experimentally in ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$${\mathrm{Sr}}_{\mathit{x}}$${\mathrm{CuO}}_{4}$, and with the states filling the insulating gap in photoemission experiments. An overall good agreement with experiments in the normal state was found. Regarding the possibility of superconductivity n these models, we have studied s-, d-, and p-wave pairing correlations. Naively, the d-wave channel seems enhanced near half filling while the extended s-wave channel seems enhanced from half filling up to 40% doping.However, we found that the enhancement comes from short-distance effects and, thus, no numerical indications of superconductivity were found in these models. We emphasize the importance of analyzing the pairing correlations as a function of distance to distinguish between short- and long-distance effects in the susceptibilities. We also observed that spiral order is enhanced at small J/t and low doping. Uniform chiral rder is suppressed by dynamical holes while staggered chiral order may be enhanced, although with a small plaquette order parameter. We conclude with the observation that the simple one-band Hubbard model with intermediate values of U\ensuremath{\sim}8t may account for many of the ``anomalous'' properties of the normal state of the high-${\mathit{T}}_{\mathit{c}}$ superconductors.

DOI
10.1103/physrevb.45.10741

Related papers

QUANTUM MONTE CARLO SIMULATIONS FOR THE THREE-BAND HUBBARD MODEL

G. Dopf, Alejandro Muramatsu, Werner R. Hanke

We present in this article results of Quantum Monte Carlo simulations for the three-band Hubbard model obtained both with the grand canonical approach as well as with the T =0 formalism. A direct comparison of Quantum Monte Carlo results with experimental values for the charge-transfer gap leads to a set of parameters for the three-band Hubbard model that consistently describes salient features of the high- T c superconductors. We find long-range antiferromagnetic order in the undoped case as well as an incommensurate structure for a doping concentration of δ≈0.20. Moreover the magnetic susceptibility as a function of doping δ is in a qualitative agreement with experiments. Concerning dynamical properties, we have evaluated the spectral density [Formula: see text]. The dispersion of the state with the lowest excitation energy is in very good agreement with experimental photoemission and inverse photoemission experiments for a doping concentration of δ=0.25. In the superconducting regime, we find that the interaction vertex for the pairing correlation function in the extended s-channel shows a maximum for δ~0.20, resembling the dependence of T c on doping. This feature is not present for other symmetries of the order parameter or in other parameter regions.

International Journal of Modern Physics C•Feb 1, 1992•2 citations

Anomalous low-doping phase of the Hubbard model

C. Gröber, R. Eder, Werner R. Hanke

We present results of a systematic quantum Monte Carlo study for the single-band Hubbard model. Thereby we evaluated single-particle spectra (PES and IPES), two-particle spectra (spin and density correlation functions), and the dynamical correlation function of suitably defined diagnostic operators, all as a function of temperature and hole doping. The results allow us to identify different physical regimes. Near half-filling we find an anomalous ``Hubbard-I phase,'' where the band structure is, up to some minor modifications, consistent with the Hubbard-I predictions. At lower temperatures, where the spin response becomes sharp, additional dispersionless ``bands'' emerge due to the dressing of electrons/holes with spin excitations. We present a simple phenomenological fit that reproduces the band structure of the insulator quantitatively. The Fermi surface volume in the low-doping phase, as derived from the single-particle spectral function, is not consistent with the Luttinger theorem, but qualitatively in agreement with the predictions of the Hubbard-I approximation. The anomalous phase extends up to a hole concentration of $\ensuremath{\approx}15%,$ i.e., the underdoped region in the phase diagram of high-${T}_{c}$ superconductors. We also investigate the nature of the magnetic ordering transition in the single-particle spectra. We show that the transition to a spin-density wave-like band structure is not accomplished by the formation of any resolvable ``precursor bands,'' but rather by a (spectroscopically invisible) band of spin-3/2 quasiparticles. We discuss implications for the ``remnant Fermi surface'' in insulating cuprate compounds and the shadow bands in the doped materials.

arXiv:cond-mat/0001366•Physical review. B, Condensed matter•Aug 15, 2000•137 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

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

Anomalous superconductivity and its competition with antiferromagnetism in doped Mott insulators

S. S. Kancharla, Bumsoo Kyung, David Sénéchal +4 more

Proximity to a Mott insulating phase is likely to be an important physical ingredient of a theory that aims to describe high-temperature superconductivity in the cuprates. Quantum cluster methods are well suited to describe the Mott phase. Hence, as a step toward a quantitative theory of the competition between antiferromagnetism and $d$-wave superconductivity in the cuprates, we use cellular dynamical mean-field theory to compute zero-temperature properties of the two-dimensional square lattice Hubbard model. The $d$-wave order parameter is found to scale like the superexchange coupling $J$ for on-site interaction $U$ comparable to or larger than the bandwidth. The order parameter also assumes a dome shape as a function of doping, while, by contrast, the gap in the single-particle density of states decreases monotonically with increasing doping. In the presence of a finite second neighbor hopping ${t}^{\ensuremath{'}}$, the zero-temperature phase diagram displays the electron-hole asymmetric competition between antiferromagnetism and superconductivity that is observed experimentally in the cuprates. Adding realistic third neighbor hopping ${t}^{\ensuremath{''}}$ improves the overall agreement with the experimental phase diagram. Since band parameters can vary depending on the specific cuprate considered, the sensitivity of the theoretical phase diagram to band parameters challenges the commonly held assumption that the doping vs ${T}_{c}∕{T}_{c}^{\mathit{max}}$ phase diagram of the cuprates is universal. The calculated angle-resolved photoemission spectrum displays the observed electron-hole asymmetry. The tendency to homogeneous coexistence of the superconducting and antiferromagnetic order parameters is stronger than observed in most experiments but consistent with many theoretical results and with experiments in some layered high-temperature superconductors. Clearly, our calculations reproduce important features of $d$-wave superconductivity in the cuprates that would otherwise be considered anomalous from the point of view of the standard Bardeen--Cooper--Schrieffer approach. At strong coupling, $d$-wave superconductivity and antiferromagnetism naturally appear as two equally important competing instabilities of the normal phase of the same underlying Hamiltonian.

arXiv:cond-mat/0508205•Physical Review B•May 23, 2008•184 citations

Density of states of doped Hubbard clusters

Elbio R. Dagotto, Adriana Moreo, F. Ortolani +2 more

Numerical results for the density of states to remove or add electrons to a finite Hubbard cluster are presented for both positive and negative on-site interactions. The behavior of the filling versus the chemical potential is discussed. At half filling there is a gap in the density of states. When the repulsive-U Hubbard model is doped spectral weight is transferred into the gap, and the chemical potential is found to move across the gap as one goes from electron to hole doping. We discuss the relationship of these results to experiment. For the negative-U Hubbard model, the chemical potential shifts continuously and lies in the middle of the superconducting gap.

Physical Review Letters•Sep 30, 1991•105 citations
We study the t-J and the Hubbard models at zero temperature using exact-diagonalization techniques on \ensuremath{\surd}10 \ifmmode\times\else\texttimes\fi{} \ensuremath{\surd}10 and 4\ifmmode\times\else\texttimes\fi{}4 sites clusters. Quantum Monte Carlo simulation results on larger lattices are also presented. All electronic fillings have been analyzed for the three models. We have measured equal-time correlation functions corresponding to various types of order (ranging from standard'' staggered spin order to more exotic'' possibilities like chiral order), as well as various dynamical properties of these models. Upper bounds for the critical hole doping (xc{\mathit{x}}_{\mathit{c}}xc​), where long-range antiferromagnetic order disappears, are presented. It was found that xc{\mathit{x}}_{\mathit{c}}xc​ is very small in agreement with experiments for the high-Tc{\mathit{T}}_{\mathit{c}}Tc​ superconductors. For example, in the t-J model, xc{\mathit{x}}_{\mathit{c}}xc​0.08 at J/t=0.4. However, short-distance spin correlations are important up to much higher dopings producing a sharp well-defined spin-wave-like peak in S(q=(\ensuremath{\pi},\ensuremath{\pi}),\ensuremath{\omega}). Regarding the possibility of phase separation in the Hubbard model, we have studied the behavior of the density of particles, 〈n〉, as a function of the chemical potential, using the Lanczos method on a 4\ifmmode\times\else\texttimes\fi{}4 Hubbard cluster, finding no indications of phase separation for any value of U/t. Then, we conclude that the t-J model at small J/t should not phase separate.In order to compare theoretical predictions with photoemission experiments, we evaluated the electronic density of states N(\ensuremath{\omega}), of the Hubbard and t-J models at several doping fractions. We found that upon doping the antiferromagnetic gap is filled for U\ensuremath{\sim}8t or smaller. The chemical potential moves across the insulating gap as one goes from electron to hole doping of the half-filled cluster, in agreement with x-ray absorption experiments but at variance with photoemission experiments. We have also calculated the optical conductivity, \ensuremathσ1{\mathrm{\ensuremath{\sigma}}}_{1}\ensuremathσ1​(\ensuremath{\omega}), of the Hubbard and t-J models at all dopings on 4\ifmmode\times\else\texttimes\fi{}4 clusters. Results are compared with experiments and the weight of the Drude peak is presented as a function of couplings and dopings. Spectral weight found at small frequencies is associated with the mid-infrared band observed experimentally in {\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$${\mathrm{Sr}}_{\mathit{x}}$${\mathrm{CuO}}_{4}, and with the states filling the insulating gap in photoemission experiments. An overall good agreement with experiments in the normal state was found. Regarding the possibility of superconductivity n these models, we have studied s-, d-, and p-wave pairing correlations. Naively, the d-wave channel seems enhanced near half filling while the extended s-wave channel seems enhanced from half filling up to 40% doping.However, we found that the enhancement comes from short-distance effects and, thus, no numerical indications of superconductivity were found in these models. We emphasize the importance of analyzing the pairing correlations as a function of distance to distinguish between short- and long-distance effects in the susceptibilities. We also observed that spiral order is enhanced at small J/t and low doping. Uniform chiral rder is suppressed by dynamical holes while staggered chiral order may be enhanced, although with a small plaquette order parameter. We conclude with the observation that the simple one-band Hubbard model with intermediate values of U\ensuremath{\sim}8t may account for many of the ``anomalous'' properties of the normal state of the high-Tc{\mathit{T}}_{\mathit{c}}Tc​ superconductors.
We present in this article results of Quantum Monte Carlo simulations for the three-band Hubbard model obtained both with the grand canonical approach as well as with the T =0 formalism. A direct comparison of Quantum Monte Carlo results with experimental values for the charge-transfer gap leads to a set of parameters for the three-band Hubbard model that consistently describes salient features of the high- T c superconductors. We find long-range antiferromagnetic order in the undoped case as well as an incommensurate structure for a doping concentration of δ≈0.20. Moreover the magnetic susceptibility as a function of doping δ is in a qualitative agreement with experiments. Concerning dynamical properties, we have evaluated the spectral density [Formula: see text]. The dispersion of the state with the lowest excitation energy is in very good agreement with experimental photoemission and inverse photoemission experiments for a doping concentration of δ=0.25. In the superconducting regime, we find that the interaction vertex for the pairing correlation function in the extended s-channel shows a maximum for δ~0.20, resembling the dependence of T c on doping. This feature is not present for other symmetries of the order parameter or in other parameter regions.
We present results of a systematic quantum Monte Carlo study for the single-band Hubbard model. Thereby we evaluated single-particle spectra (PES and IPES), two-particle spectra (spin and density correlation functions), and the dynamical correlation function of suitably defined diagnostic operators, all as a function of temperature and hole doping. The results allow us to identify different physical regimes. Near half-filling we find an anomalous ``Hubbard-I phase,'' where the band structure is, up to some minor modifications, consistent with the Hubbard-I predictions. At lower temperatures, where the spin response becomes sharp, additional dispersionless ``bands'' emerge due to the dressing of electrons/holes with spin excitations. We present a simple phenomenological fit that reproduces the band structure of the insulator quantitatively. The Fermi surface volume in the low-doping phase, as derived from the single-particle spectral function, is not consistent with the Luttinger theorem, but qualitatively in agreement with the predictions of the Hubbard-I approximation. The anomalous phase extends up to a hole concentration of $\ensuremath{\approx}15%,$ i.e., the underdoped region in the phase diagram of high-Tc{T}_{c}Tc​ superconductors. We also investigate the nature of the magnetic ordering transition in the single-particle spectra. We show that the transition to a spin-density wave-like band structure is not accomplished by the formation of any resolvable ``precursor bands,'' but rather by a (spectroscopically invisible) band of spin-3/2 quasiparticles. We discuss implications for the ``remnant Fermi surface'' in insulating cuprate compounds and the shadow bands in the doped materials.
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.
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
Proximity to a Mott insulating phase is likely to be an important physical ingredient of a theory that aims to describe high-temperature superconductivity in the cuprates. Quantum cluster methods are well suited to describe the Mott phase. Hence, as a step toward a quantitative theory of the competition between antiferromagnetism and ddd-wave superconductivity in the cuprates, we use cellular dynamical mean-field theory to compute zero-temperature properties of the two-dimensional square lattice Hubbard model. The ddd-wave order parameter is found to scale like the superexchange coupling JJJ for on-site interaction UUU comparable to or larger than the bandwidth. The order parameter also assumes a dome shape as a function of doping, while, by contrast, the gap in the single-particle density of states decreases monotonically with increasing doping. In the presence of a finite second neighbor hopping ${t}^{\ensuremath{'}}$, the zero-temperature phase diagram displays the electron-hole asymmetric competition between antiferromagnetism and superconductivity that is observed experimentally in the cuprates. Adding realistic third neighbor hopping ${t}^{\ensuremath{''}}$ improves the overall agreement with the experimental phase diagram. Since band parameters can vary depending on the specific cuprate considered, the sensitivity of the theoretical phase diagram to band parameters challenges the commonly held assumption that the doping vs Tc∕Tcmax{T}_{c}∕{T}_{c}^{\mathit{max}}Tc​∕Tcmax​ phase diagram of the cuprates is universal. The calculated angle-resolved photoemission spectrum displays the observed electron-hole asymmetry. The tendency to homogeneous coexistence of the superconducting and antiferromagnetic order parameters is stronger than observed in most experiments but consistent with many theoretical results and with experiments in some layered high-temperature superconductors. Clearly, our calculations reproduce important features of ddd-wave superconductivity in the cuprates that would otherwise be considered anomalous from the point of view of the standard Bardeen--Cooper--Schrieffer approach. At strong coupling, ddd-wave superconductivity and antiferromagnetism naturally appear as two equally important competing instabilities of the normal phase of the same underlying Hamiltonian.
Numerical results for the density of states to remove or add electrons to a finite Hubbard cluster are presented for both positive and negative on-site interactions. The behavior of the filling versus the chemical potential is discussed. At half filling there is a gap in the density of states. When the repulsive-U Hubbard model is doped spectral weight is transferred into the gap, and the chemical potential is found to move across the gap as one goes from electron to hole doping. We discuss the relationship of these results to experiment. For the negative-U Hubbard model, the chemical potential shifts continuously and lies in the middle of the superconducting gap.