Validating the d-density-wave Hartree–Fock account of the cuprate phase diagram: direct DDW detection, glassy domain-wall ground state, gauge-invariant response, and parameter provenance

A 2014 Hartree–Fock computation by Laughlin (Phys. Rev. B 89, 035134) derives the high-$T_c$ cuprate phase diagram from a three-parameter Fermi-liquid Hamiltonian with three coexisting orders — spin density wave (SDW), d-wave superconductivity, and d-density wave (DDW) — identifying the pseudogap as DDW order and reproducing seven experimental observables. The paper itself acknowledges a set of unresolved issues that jointly constitute its unfinished validation program: the missing direct detection of DDW magnetic Bragg peaks (possibly masked by glassy domain formation), the unknown true disordered ground state beyond the doubled-unit-cell constraint, the omitted vertex corrections needed for a gauge-invariant electromagnetic response, the sensitivity of certain numerical predictions to $U$ and the half-filling SDW moment, and the absence of a first-principles derivation of the effective parameters. The load-bearing core — whether the pseudogap is DDW order of the kind the theory posits — remains open: detection is contested rather than settled, the glassy domain-wall ground state is uncharacterized, and the gauge-invariant response has not been completed; only the parameter-provenance item has been largely addressed in general form by adjacent ab-initio work.

Background

High-temperature (high-$T_c$) cuprate superconductors are layered copper-oxide materials that become superconducting when carriers are doped into an antiferromagnetic insulating state (a Mott insulator, in which electron localization arises from Coulomb repulsion rather than band filling) at half filling. The phase diagram versus hole doping $p$ contains, besides the superconducting dome, a pseudogap regime at underdoping in which the low-energy electronic density of states is partially suppressed far above the superconducting transition. What the pseudogap is — a genuine broken-symmetry order, a crossover, or a consequence of competing orders — remains one of the central open questions of cuprate physics.

Laughlin (Phys. Rev. B 89, 035134 (2014); companion paper Phys. Rev. Lett. 112, 017004 (2014)) constructed a quantitative mean-field account of this phase diagram. The starting point is a Fermi-liquid Hamiltonian that adiabatically deforms the density-functional band structure and retains three effective interaction parameters: a copper-site repulsion $U$, a superexchange-like bond interaction $J$, and an intersite Cooper-pair repulsion $V_c$ (two further symmetry-allowed parameters are set to zero by experiment). The equations are solved at Hartree–Fock level — each electron moving in the self-consistent average field of the others — allowing three coexisting order parameters: (i) a spin density wave (SDW), the antiferromagnetic order that doubles the unit cell; (ii) d-wave superconductivity (DWS); and (iii) a d-density wave (DDW), also called orbital antiferromagnetism, in which circulating charge currents alternate from bond to bond, breaking translational symmetry with a doubled unit cell and breaking time-reversal symmetry — the 'staggered flux' state earlier proposed as cuprate hidden order by Chakravarty, Laughlin, Morr, and Nayak. The three parameters are fixed by three experimental inputs — the half-filling antiferromagnetic moment, the superconducting gap at optimal doping, and the maximum pseudogap, which the theory identifies as DDW — and the computation then reproduces the observed quantum phase transitions near $p\approx0.05$ and $p\approx0.16$, the insulating behavior below $p\approx0.05$, the d-wave pseudogap quasiparticle spectrum, the doping dependences of the pseudogap and superconducting gaps, the superconducting transition temperature $T_c(p)$, the London penetration depth $\lambda(p)$ (the depth to which a magnetic field penetrates a superconductor, fixed by the superfluid density), and the spin-wave velocity.

The paper itself, however, identifies a set of issues that its validation program explicitly leaves open. First, the required DDW order has not been directly detected: neutron searches for its characteristic magnetic Bragg peaks (elastic scattering at the doubled-cell wave vectors, produced by the tiny orbital moments) have not produced unambiguous signals, and the paper argues that pseudogap 'glassiness' — spatially disordered, sample-dependent domain formation seen in scanning tunneling microscopy — may mask them, a claim that itself needs proof. Second, the computations constrain SDW and DDW order to be periodic in the doubled unit cell; relaxing this constraint allows domain walls (interfaces between symmetry-related domains) that trap carriers and drive insulation, a glassy phenomenon the translationally periodic Hartree–Fock treatment cannot capture, so the true ground state in the presence of disorder and the energetics and dynamics of domain walls are unknown. Third, the DDW order parameter is gauge-variant (it carries charge current and changes under the phase rotation defining electromagnetic gauge invariance) while the Hamiltonian is not, and the paper's penetration-depth evaluation omits the vertex corrections (correction terms in the linear-response kernel that enforce gauge invariance, i.e. charge conservation), so the microscopic justification of its response formula is incomplete. Fourth, several numerical predictions — the chemical-potential jump between hole- and electron-type doping and the magnon spectrum near the Brillouin-zone edge — are sensitive to the chosen $U$ and the half-filling SDW moment. Fifth, the parameters $(U,J,V_c)$ are rationalized only by a bond-level perturbative argument — superexchange mediated by the bonding oxygen atom of the Cu–O plane — and not derived from the real material chemistry by first principles.

Problem Statement

Complete — or refute — the validation program of the d-density-wave Hartree–Fock account of the cuprate phase diagram (Laughlin, Phys. Rev. B 89, 035134 (2014)), whose common scientific objective is to establish whether the cuprate pseudogap is d-density-wave order of the kind the theory posits. Three issues are jointly load-bearing, each testing a different consequence of the same hypothesis. (1) Detection: either detect the DDW order directly — magnetic Bragg peaks, or equivalent orbital-moment signatures, at the doubled-cell wave vectors with the predicted orbital-current form factor and moment magnitude — or establish a quantitative upper bound on such orbital moments that excludes the DDW account of the pseudogap, adjudicating in the process whether glassy domain formation masks the signal as the source paper conjectures. (2) Ground state: determine the true ground state of the effective Hamiltonian once the doubled-unit-cell periodicity constraint is relaxed, including whether it contains SDW and DDW domain walls; if such walls occur, characterize their pattern, energetics, and carrier-trapping dynamics and the role of quenched disorder, and if the ground state is instead uniform, containing no domain walls, that is likewise a complete outcome of (2). In either case, determine whether the ground state is consistent with the insulating, glassy pseudogap phenomenology observed experimentally. (3) Response: construct the gauge-invariant electromagnetic response of the coexisting SDW/DDW/DWS state, executing the vertex corrections the source paper defers, and determine whether its penetration-depth and oscillator-strength predictions used for comparison with experiment survive. Supporting the same computations, though not entering the validation criterion below, an answer may also tighten (4) the effective $U$ and the half-filling SDW moment enough to remove the stated sensitivity of the chemical-potential jump and zone-edge magnon predictions, and (5) derive the bond-level parameters $(U,J,V_c)$ — in particular the oxygen-mediated superexchange of the bond argument — from the real material chemistry by first-principles computation, and test the fitted values ($U=0.76t$, $J=0.75t$, $V_c=0.87t$, $t=0.19$ eV) against them. The account is validated when the DDW order is detected (or its non-detection is shown to follow from glassiness established under (2)), the ground state of (2) is characterized and found consistent with the observed insulating, glassy pseudogap phenomenology, and the response theory of (3) is gauge-invariant — satisfying the Ward identities after the vertex correction — with penetration-depth and oscillator-strength predictions that agree with the measured data; it is refuted when the orbital moments are bounded below every value compatible with the pseudogap energy scale, when the ground state or response structure is shown inconsistent with the theory's own premises, or when the ground state or the penetration-depth and oscillator-strength predictions are shown inconsistent with the observed pseudogap phenomenology.

References

1. R. B. Laughlin, 'Hartree-Fock computation of the high-Tc cuprate phase diagram', Physical Review B 89, 035134 (2014). DOI: 10.1103/PhysRevB.89.035134, https://doi.org/10.1103/PhysRevB.89.035134 (preprint arXiv:1306.5359, https://arxiv.org/abs/1306.5359).
2. R. B. Laughlin, 'Fermi-Liquid Computation of the Phase Diagram of High-Tc Cuprate Superconductors with an Orbital Antiferromagnetic Pseudogap', Physical Review Letters 112, 017004 (2014). DOI: 10.1103/PhysRevLett.112.017004, https://doi.org/10.1103/PhysRevLett.112.017004
3. Sudip Chakravarty, R. B. Laughlin, Dirk K. Morr, Chetan Nayak, 'Hidden order in the cuprates', Physical Review B 63, 094503 (2001). DOI: 10.1103/PhysRevB.63.094503, https://doi.org/10.1103/PhysRevB.63.094503
4. S. Strässle, B. Graneli, M. Mali, J. Roos, H. Keller, 'Absence of Orbital Currents in Superconducting YBa2Cu4O8 Using a Zeeman-Perturbed Nuclear-Quadrupole-Resonance Technique', Physical Review Letters 106, 097003 (2011). DOI: 10.1103/PhysRevLett.106.097003, https://doi.org/10.1103/PhysRevLett.106.097003
5. Kazuhiro Fujita, Mohammad H. Hamidian, Stephen D. Edkins, Chung Koo Kim, Yuhki Kohsaka, Masaki Azuma, Mikio Takano, Hidenori Takagi, Hiroshi Eisaki, Shin-ichi Uchida, Andrea Allais, Michael J. Lawler, Eun-Ah Kim, Subir Sachdev, J. C. Séamus Davis, 'Direct phase-sensitive identification of a d-form factor density wave in underdoped cuprates', Proceedings of the National Academy of Sciences 111, E3026–E3032 (2014). DOI: 10.1073/pnas.1406297111, https://doi.org/10.1073/pnas.1406297111
6. A. M. Tsvelik, A. V. Chubukov, 'Composite charge order in the pseudogap region of the cuprates', Physical Review B 89, 184515 (2014). DOI: 10.1103/PhysRevB.89.184515, https://doi.org/10.1103/PhysRevB.89.184515
7. M. H. Hamidian, S. D. Edkins, Chung Koo Kim, J. C. Davis, A. P. Mackenzie, H. Eisaki, S. Uchida, M. J. Lawler, E.-A. Kim, S. Sachdev, K. Fujita, 'Atomic-scale electronic structure of the cuprate d-symmetry form factor density wave state', Nature Physics 12, 150–156 (2016). DOI: 10.1038/nphys3519, https://doi.org/10.1038/nphys3519
8. T. P. Croft, E. Blackburn, J. Kulda, Ruixing Liang, D. A. Bonn, W. N. Hardy, S. M. Hayden, 'No evidence for orbital loop currents in charge-ordered YBa2Cu3O6+x from polarized neutron diffraction', Physical Review B 96, 214504 (2017). DOI: 10.1103/PhysRevB.96.214504, https://doi.org/10.1103/PhysRevB.96.214504
9. P. Bourges, Y. Sidis, L. Mangin-Thro, 'Comment on "No evidence for orbital loop currents in charge-ordered YBa2Cu3O6+x from polarized neutron diffraction"', Physical Review B 98, 016501 (2018). DOI: 10.1103/PhysRevB.98.016501, https://doi.org/10.1103/PhysRevB.98.016501
10. T. P. Croft, E. Blackburn, J. Kulda, Ruixing Liang, D. A. Bonn, W. N. Hardy, S. M. Hayden, 'Reply to "Comment on 'No evidence for orbital loop currents in charge-ordered YBa2Cu3O6+x from polarized neutron diffraction'"', Physical Review B 98, 016502 (2018). DOI: 10.1103/PhysRevB.98.016502, https://doi.org/10.1103/PhysRevB.98.016502
11. Geremia Massarelli, Bryce Wu, Arun Paramekanti, 'Orbital Edelstein effect from density-wave order', Physical Review B 100, 075136 (2019). DOI: 10.1103/PhysRevB.100.075136, https://doi.org/10.1103/PhysRevB.100.075136
12. Motoaki Hirayama, Youhei Yamaji, Takahiro Misawa, Masatoshi Imada, 'Ab initio effective Hamiltonians for cuprate superconductors', Physical Review B 98, 134501 (2018). DOI: 10.1103/PhysRevB.98.134501, https://doi.org/10.1103/PhysRevB.98.134501
13. Imam Makhfudz, 'd-Density Wave Scenario Description of the New Hidden Charge Order in Cuprates', Journal of the Physical Society of Japan 85, 064701 (2016). DOI: 10.7566/JPSJ.85.064701, https://doi.org/10.7566/JPSJ.85.064701
14. S. W. Lovesey, D. D. Khalyavin, 'Neutron scattering by Dirac multipoles', Journal of Physics: Condensed Matter 29, 215603 (2017). DOI: 10.1088/1361-648X/aa5ad8, https://doi.org/10.1088/1361-648X/aa5ad8
15. Shengtao Jiang, Douglas J. Scalapino, Steven R. White, 'Density matrix renormalization group based downfolding of the three-band Hubbard model: Importance of density-assisted hopping', Physical Review B 108, L161111 (2023). DOI: 10.1103/PhysRevB.108.L161111, https://doi.org/10.1103/PhysRevB.108.L161111
16. R. S. Markiewicz, A. Bansil, 'Theory of cuprate pseudogap as antiferromagnetic order with charged domain walls', Physical Review B 109, 045116 (2024). DOI: 10.1103/PhysRevB.109.045116, https://doi.org/10.1103/PhysRevB.109.045116
17. J. M. Tranquada, P. M. Lozano, Juntao Yao, G. D. Gu, Qiang Li, 'From nonmetal to strange metal at the stripe-percolation transition in La2−xSrxCuO4', Physical Review B 109, 184510 (2024). DOI: 10.1103/PhysRevB.109.184510, https://doi.org/10.1103/PhysRevB.109.184510