Pith. sign in

REVIEW 2 major objections 5 minor 35 references

Iterative embedding and direct minimization of ghost-Gutzwiller, though formally equivalent, disagree on paramagnetic Mott insulators under a magnetic field.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-30 10:46 UTC pith:NPISKMQC

load-bearing objection Clear operational result: iterative ghost-GA embedding fails under Zeeman field in the paramagnetic Mott phase while direct minimization of the same functional does not. the 2 major comments →

arxiv 2607.27156 v1 pith:NPISKMQC submitted 2026-07-29 cond-mat.str-el

Direct minimization versus iterative embedding in the ghost-Gutzwiller method: a comparative study of magnetism in Mott insulators

Antonio Maria Tagliente , Ivan Pasqua , Michele Fabrizio This is my paper
classification cond-mat.str-el PACS 71.27.+a71.30.+h75.10.Lp
keywords ghost-GutzwillerMott insulatorquantum embeddingHubbard modelparamagnetismantiferromagnetismvariational methodsspin susceptibility
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Quantum embedding methods have long struggled to describe a Mott insulator that keeps full spin symmetry, because iterative solvers either break that symmetry or need hand-built fixes. This paper compares two ways to solve the ghost-Gutzwiller approximation for the single-band Hubbard model: a cheap iterative embedding loop, and direct minimization of the variational energy. In the paramagnetic Mott phase the iterative route needs ad-hoc recipes; in a Zeeman field those recipes fail, producing an energy jump and a spurious fully polarized insulator. Direct minimization avoids the artifacts, giving a smooth, only partly polarized paramagnetic solution with finite spin susceptibility. When antiferromagnetic order is allowed, the iterative scheme works and matches dynamical mean-field theory. The practical message is when the fast loop can be trusted and when the harder minimization is required.

Core claim

Across the Mott transition of the single-band Hubbard model, the iterative embedding scheme and direct minimization of the ghost-Gutzwiller functional behave very differently despite formal equivalence. The iterative scheme is efficient but fragile in the symmetry-invariant Mott phase: ad-hoc degeneracy recipes fail under a Zeeman field, yielding a discontinuous energy and a spurious fully polarized insulator. Direct minimization stabilizes a genuinely paramagnetic, only partially polarized Mott solution with continuous energy and finite zero-field spin susceptibility. When antiferromagnetism is allowed, the iterative scheme recovers the correct solution and aligns with dynamical mean-field

What carries the argument

The ghost-Gutzwiller variational energy functional, solved either by iterative embedding (self-consistent impurity ground state, analogous to dynamical mean-field theory) or by direct minimization of a constrained impurity wavefunction that keeps a decoupled Fermi-level bath and enforces half-filling and zero total Sz.

Load-bearing premise

The compact five-parameter impurity wavefunction built from a singly-occupied decoupled bath is assumed to span the true ghost-Gutzwiller minimum in the Mott phase, including under a magnetic field.

What would settle it

Unrestricted direct minimization of the full ghost-Gutzwiller functional in a finite Zeeman field: if that calculation recovers the iterative scheme's energy jump and fully polarized insulator, the claim that direct minimization is required for a physical paramagnetic Mott solution fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Iterative embedding can be trusted when symmetry breaking such as antiferromagnetism is allowed.
  • Direct minimization is required for symmetry-invariant Mott states and under perturbations that split atomic degeneracy.
  • Ghost-Gutzwiller with direct minimization yields finite zero-field spin susceptibility continuously across the Mott transition, unlike standard Gutzwiller.
  • The direct route lets one force fixed order parameters and explore metastable paramagnetic or spin-liquid-like trial states without artificial spin symmetrization.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same iterative fragility likely appears in other zero-temperature embedding methods whenever fields or multiplet structure split the atomic manifold.
  • Constrained impurity ansatzes that stay inside the representable density-matrix domain may become a practical fix whenever Lagrange-multiplier embedding hits a boundary minimum.
  • Direct variational control of the impurity wavefunction is a natural route to forced spin-charge fractionalization and multi-orbital spin-liquid trial states.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The manuscript compares two formally equivalent routes to the ghost-Gutzwiller (ghost-GA) saddle point for the half-filled single-band Hubbard model: an iterative quantum-embedding loop (analogous to DMFT) and direct minimization of the variational energy functional. In the paramagnetic sector the iterative scheme requires ad-hoc degeneracy resolution in the Mott phase and, under a small Zeeman field, collapses to a discontinuous energy and a spurious fully polarized insulator identical to ordinary GA. Direct minimization with a constrained five-parameter impurity ansatz yields a continuous energy, a only partially polarized Mott state, and a finite spin susceptibility across the transition. When antiferromagnetic order is allowed, the iterative scheme converges cleanly and agrees closely with DMFT (and with HF/GA at strong coupling). The authors conclude that iterative embedding is reliable once symmetry breaking is permitted, but that direct minimization is required for symmetry-invariant Mott states.

Significance. The work cleanly isolates a practical failure mode of iterative embedding that is already known in zero-temperature DMFT and shows that the same pathology appears inside ghost-GA, while also demonstrating a concrete variational fix. The comparison is computational and falsifiable: energies, magnetization, Z and susceptibility are reported side-by-side for iterative embedding, direct minimization, standard GA, HF and DMFT (AF data from Amaricci). The result is useful for anyone using ghost-GA (or related embedding methods) to study paramagnetic Mott insulators, spin liquids or other symmetry-enforced states, and it delineates when the cheaper iterative loop can be trusted. The variational inequality under finite field—direct minimization already lower than the iterative solution—is a particularly strong point.

major comments (2)
  1. [Sec. III, Eqs. (18)–(20)] Sec. III, Eqs. (18)–(20): the direct-minimization results rest on a five-parameter constrained Lanczos ansatz built from a singly-occupied decoupled bath at half-filling and Sz=0. While the h=0 agreement with the embedding solution (Fig. 2) and the strict variational superiority under h=0.01 (Fig. 3, left) already prove that the iterative fully-polarized state is not the ghost-GA minimum, the manuscript should state more explicitly that the ansatz is a restricted submanifold and that any incompleteness can only lower the direct-min energy further. A short remark on whether an unconstrained optimization of the full impurity coefficients in (5) was attempted (or is feasible) would strengthen the claim that the partially polarized solution is the true ghost-GA minimum rather than an artifact of the parametrization.
  2. [Sec. II, Eq. (13)] Sec. II, Lagrangian (13) and the surrounding discussion of representability: the argument that the Mott minimum can lie on the boundary of the admissible region for R and Delta is central to why the embedding scheme fails. The paper would benefit from a sharper statement of which matrix elements of R vanish at Uc and how that maps onto the loss of interior critical points of L. A brief diagnostic (e.g., monitoring the smallest singular value of R across the transition) would make the boundary-minimum claim more concrete and transferable to multi-orbital or crystal-field cases mentioned in the text.
minor comments (5)
  1. [Abstract] Abstract, first sentence: typographical error “long-standing ing challenge” should be “long-standing challenge”.
  2. [Fig. 3] Fig. 3 right panel: the susceptibility is assembled from embedding (metal) and direct minimization (Mott). Please state explicitly in the caption how the metallic chi is extracted and whether the two branches join continuously at Uc within numerical resolution.
  3. [Sec. IV] Sec. IV: the AF direct-minimization ansatz (24) is less accurate than iterative embedding (Fig. 5). A one-sentence remark on why the constrained form is still useful (fixed-m scans, mutual information in Fig. 6) would help the reader.
  4. [Sec. III] Notation: N_ghosts is fixed to 2 throughout; a brief comment on whether the iterative failure under Zeeman field persists for larger bath sizes would be welcome, even if only as a qualitative expectation.
  5. [References] References: the unifying ghost-GA/DMFT work cited as arXiv:2603.20559 is central to the embedding formulation; ensure the citation is updated if a journal version appears before publication.

Circularity Check

0 steps flagged

No significant circularity: computational comparison of two solvers for the same variational functional, not a self-referential derivation.

full rationale

The paper’s load-bearing content is a side-by-side numerical comparison of two formally equivalent routes to the ghost-GA stationary point (iterative embedding of the Lagrange-multiplier functional L versus direct saddle-point search on the original energy functional F with a constrained impurity ansatz). Energies, magnetizations, Z, and susceptibility are obtained by solving those equations and are benchmarked against DMFT, HF, and standard GA. No parameter is fitted to data and re-labeled as a prediction; no uniqueness theorem is imported from the authors to forbid alternatives; the five-parameter Mott ansatz (Eqs. 18–20) is constructed explicitly in the text rather than smuggled in by citation. Self-citations point to the ghost-GA formalism and prior applications used as methods, not as the comparative result being established. The derivation chain is therefore self-contained and externally falsifiable by re-running the two solvers. Score 0 is appropriate.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The work sits entirely inside the existing ghost-GA / Gutzwiller-approximation framework on infinite-coordination lattices. Load-bearing inputs are the standard Hubbard model, the ghost-GA energy functional and constraints from prior literature, the choice N_ghosts=2, and a hand-built Mott-phase impurity ansatz. No new physical entities are postulated; free parameters are methodological (ghost count, field strength, ansatz angles), not fitted to external data to produce the claim.

free parameters (3)
  • N_ghosts = 2
    Fixed to 2 (impurity plus three baths) throughout; controls the quality of the ghost-GA approximation relative to DMFT and is not varied systematically here.
  • Zeeman field h = 0.01
    Working value h=0.01 used for all finite-field scans; small enough to probe linear response but chosen by hand.
  • Impurity ansatz angles/coefficients (φ, α_n)
    Variational parameters of the compact Mott-phase wavefunction (19); optimized, not fitted to external experiment, but the ansatz form itself is a restricted manifold chosen by the authors.
axioms (4)
  • domain assumption On infinite-coordination lattices the ghost-Gutzwiller expectation values and constraints (4)–(11) are exact evaluations of the variational wavefunction; at finite coordination the same formulas define the ghost-GA approximation.
    Stated in Sec. II; underpins every energy reported. Standard in GA/DMFT literature.
  • ad hoc to paper Promoting one-body density matrices Δ and R to free variables in the embedding Lagrangian (13) is valid only in the interior of the representability region; the Mott minimum can sit on the boundary where some R entries vanish.
    Sec. II paragraph after Eq. (13); this is the paper’s diagnosis of why iterative embedding fails, and is load-bearing for the interpretation (though the numerical failure itself is independent).
  • domain assumption Half-filled single-band Hubbard model on the infinite-coordination Bethe lattice with semicircular DOS (D=2, t=1) is the appropriate testbed for a hypothetical paramagnetic Mott insulator.
    Sec. III, Eqs. (14)–(15); standard benchmark, but excludes frustration needed for a true spin liquid.
  • ad hoc to paper The five-parameter constrained Lanczos ansatz (18)–(20) adequately represents the ghost-GA Mott insulator at Sz=0 with a decoupled Fermi-level bath.
    Sec. III; validated only by matching iterative m→0 energies above Uc at h=0, then used under finite h where iterative fails.

pith-pipeline@v1.2.0-daily-grok45 · 18148 in / 3169 out tokens · 65134 ms · 2026-07-30T10:46:47.924744+00:00 · methodology

0 comments
read the original abstract

Accurately describing a hypothetical symmetry-invariant Mott insulator presents a long-standing ing challenge in iterative quantum embedding methods. We address this issue within the ghost- Gutzwiller method, which can be solved either through an iterative embedding scheme, analogous to dynamical mean-field theory, or by directly minimizing its variational energy functional. Across the Mott transition of the single-band Hubbard model, these formally equivalent approaches behave very differently: the iterative scheme is computationally efficient but fragile, necessitating ad-hoc recipes in the Mott phase that fail in a Zeeman field, leading to a discontinuous energy and a spurious fully-polarized insulator. Direct minimization avoids these artifacts, stabilizing a genuinely paramagnetic solution. Conversely, when symmetry breaking is allowed, as in an antiferromagnetic phase, the iterative scheme yields the correct solution, closely aligning with dynamical mean-field theory. Our findings delineate the conditions under which the iterative embedding can be trusted and when direct minimization is instead required.

Figures

Figures reproduced from arXiv: 2607.27156 by Antonio Maria Tagliente, Ivan Pasqua, Michele Fabrizio.

Figure 1
Figure 1. Figure 1: FIG. 1. Ghost-GA energy [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Comparison between the energies per site obtained [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Ghost-GA results for the single-band Hubbard model in a magnetic field [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Total energy [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison between the energies and the magneti [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Comparison between the true ground state within the [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 3 linked inside Pith

  1. [1]

    Sun and G

    Q. Sun and G. K.-L. Chan, Quantum Embedding Theo- ries, Accounts of Chemical Research49, 2705 (2016)

    2016

  2. [2]

    Georges, G

    A. Georges, G. Kotliar, W. Krauth, and M. J. Rozen- berg, Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions, Rev. Mod. Phys.68, 13 (1996)

    1996

  3. [3]

    G. Kotliar, Landau theory of the Mott transition in the fully frustrated Hubbard model in infinite dimensions, The European Physical Journal B - Condensed Matter and Complex Systems11, 27 (1999)

    1999

  4. [4]

    Giuli, T.-H

    S. Giuli, T.-H. Lee, Y.-X. Yao, G. Kotliar, A. E. Ruck- enstein, O. Gingras, and N. Lanat` a, Unifying Vari- ational and Dynamical Quantum Embedding: From Ghost Gutzwiller Approximation to Dynamical Mean- Field Theory, arXiv:2603.20559 [cond-mat.str-el] (2026)

    arXiv 2026

  5. [5]

    Wen, S.-L

    J. Wen, S.-L. Yu, S. Li, W. Yu, and J.-X. Li, Experimen- tal identification of quantum spin liquids, npj Quantum Materials4, 12 (2019). 9

    2019

  6. [6]

    J. R. Chamorro, T. M. McQueen, and T. T. Tran, Chem- istry of Quantum Spin Liquids, Chemical Reviews121, 2898 (2021)

    2021

  7. [7]

    Matsuda, T

    Y. Matsuda, T. Shibauchi, and H.-Y. Kee, Kitaev quan- tum spin liquids, Rev. Mod. Phys.97, 045003 (2025)

    2025

  8. [8]

    Savary and L

    L. Savary and L. Balents, Quantum spin liquids: a re- view, Reports on Progress in Physics80, 016502 (2016)

    2016

  9. [9]

    Lanat` a, T.-H

    N. Lanat` a, T.-H. Lee, Y.-X. Yao, and V. Dobrosavljevi´ c, Emergent Bloch excitations in Mott matter, Phys. Rev. B96, 195126 (2017)

    2017

  10. [10]

    Lanat` a, Y

    N. Lanat` a, Y. Yao, C.-Z. Wang, K.-M. Ho, and G. Kotliar, Phase Diagram and Electronic Structure of Praseodymium and Plutonium, Phys. Rev. X5, 011008 (2015)

    2015

  11. [11]

    Guerci, M

    D. Guerci, M. Capone, and M. Fabrizio, Exciton Mott transition revisited, Phys. Rev. Mater.3, 054605 (2019)

    2019

  12. [12]

    M. S. Frank, T.-H. Lee, G. Bhattacharyya, P. K. H. Tsang, V. L. Quito, V. Dobrosavljevi´ c, O. Christiansen, and N. Lanat` a, Quantum embedding description of the Anderson lattice model with the ghost Gutzwiller ap- proximation, Phys. Rev. B104, L081103 (2021)

    2021

  13. [13]

    Mejuto-Zaera and M

    C. Mejuto-Zaera and M. Fabrizio, Efficient computa- tional screening of strongly correlated materials: Mul- tiorbital phenomenology within the ghost Gutzwiller ap- proximation, Phys. Rev. B107, 235150 (2023)

    2023

  14. [14]

    T.-H. Lee, N. Lanat` a, and G. Kotliar, Accuracy of ghost rotationally invariant slave-boson and dynamical mean field theory as a function of the impurity-model bath size, Phys. Rev. B107, L121104 (2023)

    2023

  15. [15]

    T.-H. Lee, C. Melnick, R. Adler, N. Lanat` a, and G. Kotliar, Accuracy of ghost-rotationally-invariant slave- boson theory for multiorbital Hubbard models and real- istic materials, Phys. Rev. B108, 245147 (2023)

    2023

  16. [16]

    Mejuto-Zaera, Quantum embedding for molecules us- ing auxiliary particles – the ghost Gutzwiller Ansatz, Faraday Discuss.254, 653 (2024)

    C. Mejuto-Zaera, Quantum embedding for molecules us- ing auxiliary particles – the ghost Gutzwiller Ansatz, Faraday Discuss.254, 653 (2024)

    2024

  17. [17]

    A. M. Tagliente, C. Mejuto-Zaera, and M. Fabrizio, Re- vealing spinons by proximity effect, Phys. Rev. B111, 125110 (2025)

    2025

  18. [18]

    S. Smit, M. Bluschke, P. Moen, N. Heinsdorf, E. Za- vatti, G. Bellomia, S. Giuli, S. K. Y. Dufresne, C. T. Suen, V. Zimmermann, C. Au-Yeung, S. Zhdanovich, J. I. Dadap, M. Zonno, S. Gorovikov, H. Lee, C.-T. Kuo, J.-S. Lee, D. Song, S. Ishida, H. Eisaki, B. Keimer, M. Michiardi, I. S. Elfimov, G. Levy, D. J. Jones, M. Capone, and A. Damascelli, Enhanced co...

    Pith/arXiv arXiv 2025

  19. [19]

    Giuli, C

    S. Giuli, C. Mejuto-Zaera, and M. Capone, Alter- magnetism from interaction-driven itinerant magnetism, Phys. Rev. B111, L020401 (2025)

    2025

  20. [20]

    Pasqua, A

    I. Pasqua, A. M. Tagliente, G. Bellomia, B. Monser- rat, M. Fabrizio, and C. Mejuto-Zaera, Band-structure picture for topology in strongly correlated systems with the ghost Gutzwiller ansatz, Phys. Rev. B113, 155133 (2026)

    2026

  21. [21]

    Bellomia, A

    G. Bellomia, A. Amaricci, and M. Capone, Local clas- sical correlations between physical electrons in Hubbard systems, Phys. Rev. B113, 155158 (2026)

    2026

  22. [22]

    Mejuto-Zaera and M

    C. Mejuto-Zaera and M. Fabrizio, Ghost Embed- ding Bridging Chemistry and One-Body Theories, arXiv:2602.17164 [cond-mat.str-el] (2026)

    arXiv 2026

  23. [23]

    Giuli, H

    S. Giuli, H. Hasan, B. Kloss, M. S. Frank, T.-H. Lee, O. Gingras, Y.-X. Yao, and N. Lanat` a, Linear foun- dation model for quantum embedding: Data-driven compression of the ghost Gutzwiller variational space, arXiv:2512.21666 [cond-mat.str-el] (2025)

    arXiv 2025

  24. [24]

    I.-C. Chen, A. Khindanov, C. M. Salazar, H. M. Barona, G. Harrabi, F. Zhang, C.-Z. Wang, T. Iadecola, N. Lanat` a, and Y.-X. Yao, Quantum-classical embedding via ghost Gutzwiller approximation for enhanced sim- ulations of correlated electron systems, npj Quantum In- formation12(2026)

    2026

  25. [25]

    Kazemi-Moridani, S

    A. Kazemi-Moridani, S. Giuli, T.-H. Lee, A. M. S. Trem- blay, M. Cˆ ot´ e, N. Lanat` a, and O. Gingras, Magnetic phases of the anisotropic triangular Hubbard model from the ghost-Gutzwiller approximation in the rotating spin- frame, arXiv:2604.24848 [cond-mat.str-el] (2026)

    Pith/arXiv arXiv 2026

  26. [26]

    Pasqua, C

    I. Pasqua, C. Mejuto-Zaera, and M. Fabrizio, Fraction- alized Fermi liquids with the ghost-Gutzwiller Ansatz, arXiv:2607.12410 [cond-mat.str-el] (2026)

    Pith/arXiv arXiv 2026

  27. [27]

    Fabrizio, Quantum fluctuations beyond the Gutzwiller approximation, Phys

    M. Fabrizio, Quantum fluctuations beyond the Gutzwiller approximation, Phys. Rev. B95, 075156 (2017)

    2017

  28. [28]

    A. A. Klyachko, Quantum marginal problem and N- representability, Journal of Physics: Conference Series 36, 72 (2006)

    2006

  29. [29]

    ˇZitko and M

    R. ˇZitko and M. Fabrizio,Z 2 gauge theory description of the Mott transition in infinite dimensions, Phys. Rev. B 91, 245130 (2015)

    2015

  30. [30]

    Bauer and A

    J. Bauer and A. C. Hewson, Field-dependent quasipar- ticles in the infinite-dimensional Hubbard model, Phys. Rev. B76, 035118 (2007)

    2007

  31. [31]

    The DMFT data have been kindly provided to us by Adriano Amaricci

  32. [32]

    Crippa, I

    L. Crippa, I. Krivenko, S. Giuli, G. Bellomia, A. Kowal- ski, F. Petocchi, A. Scazzola, M. Wallerberger, G. Mazza, L. de Medici, G. Sangiovanni, M. Capone, and A. Amar- icci, Next-generation EDIpack: A Lanczos-based package for quantum impurity models featuring general broken- symmetry phases, flexible bath topologies and multi- platform interoperability,...

    2025

  33. [33]

    Crippa, I

    L. Crippa, I. Krivenko, S. Giuli, G. Bellomia, A. Kowal- ski, F. Petocchi, A. Scazzola, M. Wallerberger, G. Mazza, L. de Medici, G. Sangiovanni, M. Capone, and A. Amar- icci, Codebase release 5.3 for EDIpack, SciPost Phys. Codebases 58 (2025)

    2025

  34. [34]

    Sangiovanni, A

    G. Sangiovanni, A. Toschi, E. Koch, K. Held, M. Capone, C. Castellani, O. Gunnarsson, S.-K. Mo, J. W. Allen, H.- D. Kim, A. Sekiyama, A. Yamasaki, S. Suga, and P. Met- calf, Static versus dynamical mean-field theory of Mott antiferromagnets, Phys. Rev. B73, 205121 (2006)

    2006

  35. [35]

    Taranto, G

    C. Taranto, G. Sangiovanni, K. Held, M. Capone, A. Georges, and A. Toschi, Signature of antiferromagnetic long-range order in the optical spectrum of strongly cor- related electron systems, Phys. Rev. B85, 085124 (2012)

    2012