Quantum criticality and incipient phase separation in the thermodynamic properties of the Hubbard model.

Transport measurements on the cuprates suggest the presence of a quantum critical point (QCP) hiding underneath the superconducting dome near optimal hole doping. We provide numerical evidence in support of this scenario via a dynamical cluster quantum Monte Carlo study of the extended two-dimensional Hubbard model. Single-particle quantities, such as the spectral function, the quasi-particle weight and the entropy, display a crossover between two distinct ground states: a Fermi liquid at low filling and a non-Fermi liquid with a pseudo-gap at high filling. Both states are found to cross over to a marginal Fermi-liquid state at higher temperatures. For finite next-nearest-neighbour hopping t', we find a classical critical point at temperature T(c). This classical critical point is found to be associated with a phase-separation transition between a compressible Mott gas and an incompressible Mott liquid corresponding to the Fermi liquid and the pseudo-gap state, respectively. Since the critical temperature T(c) extrapolates to zero as t' vanishes, we conclude that a QCP connects the Fermi liquid to the pseudo-gap region, and that the marginal Fermi-liquid behaviour in its vicinity is the analogue of the supercritical region in the liquid-gas transition.

Proximity of the superconducting dome and the quantum critical point in the two-dimensional Hubbard model.

We use the dynamical cluster approximation to understand the proximity of the superconducting dome to the quantum critical point in the two-dimensional Hubbard model. In a BCS formalism, T(c) may be enhanced through an increase in the d-wave pairing interaction (V(d)) or the bare pairing susceptibility (χ(0d)). At optimal doping, where V(d) is revealed to be featureless, we find a power-law behavior of χ(0d)(ω=0), replacing the BCS log, and strongly enhanced T(c). We suggest experiments to verify our predictions.

Cluster solver for dynamical mean-field theory with linear scaling in inverse temperature.

Dynamical mean-field theory and its cluster extensions provide a very useful approach for examining phase transitions in model Hamiltonians and, in combination with electronic structure theory, constitute powerful methods to treat strongly correlated materials. The key advantage to the technique is that, unlike competing real-space methods, the sign problem is well controlled in the Hirsch-Fye (HF) quantum Monte Carlo used as an exact cluster solver. However, an important computational bottleneck remains; the HF method scales as the cube of the inverse temperature, ß . This often makes simulations at low temperatures extremely challenging. We present here a method based on determinant quantum Monte Carlo which scales linearly in ß , with a quadratic term that comes in to play for the number of time slices larger than hundred, and demonstrate that the sign problem is identical to HF.