Pith. sign in

REVIEW 3 major objections 5 minor 15 references

A conditional Euclidean–Hamiltonian split can replace scalar sign reweighting by a monitored active-space approximation in finite fermion models.

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 · deepseek-v4-flash

2026-08-01 21:48 UTC pith:A5QC5KG2

load-bearing objection A clean, honestly scoped proposal for trading scalar sign reweighting for matrix harvesting, but the signful benchmark never exercises the harvesting loop — worth a serious referee. the 3 major comments →

arxiv 2607.15950 v1 pith:A5QC5KG2 submitted 2026-07-17 quant-ph cond-mat.str-elhep-lat

Conditional Euclidean-Hamiltonian reductions for sign-problem toy models

Tsogtgerel Gantumur This is my paper
classification quant-ph cond-mat.str-elhep-lat
keywords sign problemconditional Euclidean-Hamiltonian reductionoperator-valued kernelcorrelation-matrix GEVPprojected transfer matrixfinite-density fermionsHubbard ringeffective Hamiltonian
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.

Conditional Euclidean–Hamiltonian (CEH) reduction is a proposal for when to turn a residual quantum sector into a Monte Carlo weight: not before. Instead of tracing the active sector into a sign-changing determinant, CEH keeps that sector operator-valued, compresses it to a low-dimensional trial space, and harvests finite correlation or transfer matrices from a Monte-Carlo-tractable reference calculation. The paper shows in three finite benchmarks that this ordering produces usable effective Hamiltonians: a two-channel oscillator demonstrates conditional basis compression, a stochastic positive-measure calculation reproduces projected GEVP energies within errors, and a four-site Hubbard ring shows that deferred finite-density continuation works when the trial space contains the relevant number sectors. The authors are explicit that this is a controlled finite-model trade, not an end-to-end sign-problem solution and not a scaling claim.

Core claim

CEH studies the split H = G + A, where G is a reference Hamiltonian with a positive or otherwise controlled Euclidean representation and A is an operator-valued residual. For each sampled reference history, the active evolution is not traced into a scalar; it is compressed to an R-by-R matrix and averaged. Equivalently, the Euclidean input is the correlation-matrix family C_IJ(n) = <Phi_I, T^n Phi_J>, which defines a compressed transfer operator Teff = C(0)^(-1/2) C(1) C(0)^(-1/2); a generalized eigenvalue problem extracts projected energies. Finite density is applied afterward through Heff(mu) = Heff(0) - mu Neff, which exactly commutes with projection only for direct Galerkin compression,

What carries the argument

The central object is the operator split H = G + A with the residual sector left untraced. The workhorse identities are the projected correlation matrices C_IJ(n) = <Phi_I, T^n Phi_J>, the compressed transfer operator Teff = C(0)^(-1/2) C(1) C(0)^(-1/2), and the deferred finite-density prescription Heff(mu) = Heff(0) - mu Neff. Together they convert a Monte Carlo reference ensemble into a finite effective Hamiltonian without ever forming a scalar determinant weight for the active sector. The construction's validity is controlled by rank enlargement, metric conditioning, stochastic matrix-element errors, and number-sector diagnostics.

Load-bearing premise

The load-bearing premise is that the residual active dynamics are low-rank approximable and that the projected correlation or transfer matrices can be harvested from the reference calculation with usable overlap, variance, and conditioning; if the required rank grows toward the full Hilbert-space dimension or the harvested metric becomes too noisy and ill-conditioned, the CEH trade collapses.

What would settle it

Perform an end-to-end stochastic CEH calculation on a spatially extended Hubbard chain or ladder at finite mu, enlarging the trial space while tracking the smallest retained metric eigenvalue and the number of GEVP directions needed for stable energies. If the required rank approaches the Hilbert-space dimension, or if the metric condition number grows faster than the trial-space dimension and no stable projected spectrum is reached, the central claim is falsified.

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

If this is right

  • In a two-channel oscillator, a channel-adapted conditional basis reduces ground-state error by orders of magnitude at small rank compared with an equal-dimensional undressed product basis.
  • The stochastic CEH loop—reference oscillator paths plus active-sector matrix observables—reproduces the deterministic projected GEVP energies within one standard error.
  • At finite budget the determinant-sign branch of the Hubbard benchmark shows severe sign cancellation, while the projected branch shows that thermal density requires weakly occupied neighboring number sectors.
  • If the CEH ordering holds, sign or phase reweighting is replaced by finite-rank matrix algebra plus diagnostics, changing the computational difficulty instead of eliminating it.
  • The paper claims no favorable scaling and no end-to-end stochastic CEH Hubbard solution; usefulness requires low-rank approximability and reference-harvestability.

Where Pith is reading between the lines

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

  • A natural next test, and the paper's stated decisive one, is a spatially extended model where the active basis is grown from reference Euclidean data; if the required rank grows gently with system size while conditioning stays under control, CEH could become a practical alternative in regimes where scalar sign reweighting fails.
  • CEH's better metric conditioning but worse energy error than a projected-density-matrix-style estimator suggests the two estimators could be combined—using CEH-correlated matrices for sector coverage and conditioning, and the other for energy extraction.
  • The deferred finite-density step implies a possible practical dividend: one reference calculation at mu = 0 could supply trial spaces for a range of chemical potentials, as long as the relevant number sectors are represented; the Hubbard sector-coverage result shows the risk when they are not.
  • In a mixed-scale application, CEH could be paired with a sign-free determinant sector: trace the sign-free part into a determinant while keeping only the signful part operator-valued, potentially reducing scalar sign fluctuations.

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

3 major / 5 minor

Summary. The paper proposes a conditional Euclidean–Hamiltonian (CEH) reduction for sign-problem toy models. The idea is to avoid the final scalarization step (trace or determinant) of the active sector in a Euclidean path integral by keeping that sector operator-valued, harvesting finite projected correlation or transfer matrices from a tractable reference calculation, and then applying finite-density dependence through a reduced effective Hamiltonian. The formalism is developed in Sections 2–3: operator-valued kernels, projected transfer matrices and GEVP (Prop. 3.1), conditional Galerkin spaces, and two finite-density continuation prescriptions. Three numerical benchmarks are presented: a two-channel oscillator for conditional basis compression, a positive-measure stochastic oscillator test for correlation-matrix harvesting and GEVP extraction with a PDMS-style comparison, and a four-site Hubbard ring with a determinant-sign stress test and an active-space sector-coverage study.

Significance. If validated end-to-end, the CEH approach could provide a practical way to convert a scalar sign problem into a monitored finite-rank matrix problem, with explicit diagnostics for rank, conditioning, and sector coverage. The paper is careful and honest about limitations: it explicitly disclaims a generic solution and identifies the prerequisites of low-rank approximability and reference-harvestability. The derivations are clean and the benchmarks are clearly described. However, the central utility claim—that CEH can replace scalar sign reweighting—is not fully demonstrated by the evidence presented, as the only signful benchmark does not exercise the stochastic CEH harvesting loop and the transfer-logarithm finite-density continuation is never tested. The paper is best viewed as a component-testing study with a promising but not-yet-complete proof of principle.

major comments (3)
  1. [§4.3, §5] The signful Hubbard benchmark never exercises the CEH harvesting loop. The determinant branch (§4.3.1) only demonstrates sign cancellation in a finite budget, while the projected branch (§4.3.2) builds trial spaces by exact diagonalization. The only stochastic CEH harvesting test (§4.2.2) is a sign-free two-channel oscillator. Thus the central claim that CEH replaces scalar sign reweighting by matrix harvesting is not demonstrated on any signful problem. The abstract's 'benchmarks support the proposed trade' is stronger than the evidence; the paper's own §5 correctly states this limitation, but the claim should be tempered or a signful harvesting test added.
  2. [§3.4, Eq. (3.19)] The deferred finite-density prescription H_eff(µ)=H_eff(0)-µN_eff is never tested. The Hubbard finite-density test in §4.3.2 uses direct Galerkin compression (Eq. 3.16), which commutes exactly with -µN on a fixed trial space. The correlation-matrix route (3.19) involves a transfer-logarithm approximation whose accuracy, conditioning, and number-sector leakage are not examined. Without a test of (3.19), the proposed finite-density CEH workflow remains purely formal.
  3. [§2.5, §5] The usefulness conditions of low-rank approximability and reference-harvestability are stated but not quantitatively assessed. The benchmarks provide partial evidence: the oscillator shows harvestability in a sign-free reference, and the Hubbard ring shows sector-coverage constraints. Neither demonstrates that a signful reference can supply the required projected matrices with usable overlap, variance, and conditioning. This is acknowledged as an open question, but it is the central premise of the method; the paper should either provide a proof-of-principle for a signful reference or explicitly frame the contribution as component tests only.
minor comments (5)
  1. [Title/Abstract] The title contains a typo: 'HAMIL TONIAN' should be 'HAMILTONIAN'.
  2. [§4.2.2, Eq. (4.10)] The factor e^{-Ωna/2} in the correlation-matrix estimator is stated without derivation. A brief explanation of how the Hermite-polynomial transition representation yields this normalization would improve reproducibility.
  3. [§4.3.1] The sign diagnostic qsign is clearly labeled as descriptive, but the text could note that the reported standard errors are over chains only and do not account for autocorrelation within chains, which is present by construction.
  4. [§4.3.2] The SVD threshold 10^{-10} max{D,R} σmax is chosen without sensitivity analysis. A sentence on the stability of the retained rank with respect to this threshold would be useful.
  5. [§5] The concluding sentence 'The benchmarks support the proposed trade' could be moved to match the more precise limitation statement in the same section: the benchmarks support the components, not the end-to-end trade in a signful setting.

Circularity Check

0 steps flagged

No significant circularity: the CEH construction is benchmarked against independent exact diagonalization, and the paper explicitly flags its own limitations.

full rationale

The paper does not fit parameters to reproduce target outputs, does not rename a fitted quantity as a prediction, and does not rely on a self-citation chain. The central construction is an explicit projection of H onto a trial space: Eq. (3.8) shows T_eff is the compression of e^{-aH}, and Eq. (3.16) shows direct Galerkin compression commutes exactly with adding -mu N. The finite-density continuation in the transfer-logarithm route is explicitly labeled as an effective model (Eq. 3.19: 'defines an effective projected model, not necessarily the direct Galerkin compression') and is never presented as exact. The benchmarks compare against independent exact diagonalization or deterministic projected references, so the claims are externally checked rather than circular. The paper also candidly identifies the gap between its central claim and its evidence, stating in §5 that the Hubbard benchmark is 'neither scalable nor an end-to-end stochastic CEH solution of the Hubbard sign problem' and in §2.5 that the approach 'is useful only when the residual dynamics are low-rank approximable and the required matrix data are reference-harvestable.' These are scope limitations, not circularity. The only resemblance to a circular step is that the stochastic CEH loop is tested on a sign-free oscillator while the signful Hubbard benchmark uses exact-diagonalization trial spaces, but the paper explicitly acknowledges this, and the absence of an end-to-end test does not make the derivation circular. No self-citations are load-bearing; all cited prior work is external and used for standard ingredients. Accordingly, the circularity score is 0.

Axiom & Free-Parameter Ledger

4 free parameters · 5 axioms · 0 invented entities

No parameter was fitted to reproduce the target energies; all constants are benchmark inputs or hand-set algorithmic settings. The model couplings (Ω, ε, t, g, U, t_hop) are Hamiltonian definitions, not fitted outputs. The chosen GEVP windows, regularization threshold, and Hubbard β/Δτ set the tests' operating points and affect the reported diagnostics, but the paper's qualitative conclusions do not rest on their precise values. No new physical entities are introduced; H_eff is a computational construct, not a postulated particle or force.

free parameters (4)
  • GEVP time window n0 and n = n0=2, n=8
    Hand-chosen reference and diagonalization times for the correlation-matrix GEVP; conditioning worsens with n0 (reported: cond grows from 4.76 at n0=1 to 514.25 at n0=4), so it is a user-set algorithmic parameter rather than a fitted constant.
  • Metric regularization threshold = 10^-10
    SVD cutoff for rank filtering in the multi-µ Hubbard active-space construction (§4.3.2). Chosen by hand; the paper reports retained ranks rR and says the construction is robust above the threshold.
  • Hubbard β and Δτ = β=8, Δτ=0.25
    Determinant-sign stress test settings (Nτ=32, one sweep attempts 128 flips). They define the sign diagnostic and are benchmark inputs, not fitted to output.
  • Euclidean step a in oscillator GEVP = 0.1
    Imaginary-time step used in the positive-measure harvesting benchmark (§4.2). A hand-set algorithmic choice that affects the symmetric-splitting error but not the qualitative conclusions.
axioms (5)
  • domain assumption Product Hilbert space decomposition H = L^2(X) ⊗ K and operator split H = G + A with A(x) acting on K.
    Used throughout §2 to define the reference/active split; the paper assumes H_x is self-adjoint and bounded below and that H_x + A(x̂) is a self-adjoint, bounded-below form sum (§2.1).
  • standard math Trotter product formula and operator-valued Feynman–Kac / heat-kernel representation.
    Eqs. (2.10)–(2.11) rely on the Trotter factorization; standard but central to the finite-slice construction of the operator-valued kernel.
  • standard math Rayleigh–Ritz / GEVP convergence: increasing trial-space dimension R converges isolated low-lying eigenvalues of T = e^{-aH}.
    Proposition 3.1 and §3.2 invoke the standard Ritz-value convergence result; the paper also relies on the GEVP stability at Euclidean times n0 and n.
  • domain assumption The reference sector admits a positive (or controllable-sign) Euclidean measure.
    Required for Monte Carlo harvesting in §4.2. The Hubbard determinant branch instead samples |det| and restores sign via reweighting (§4.3.1), so the reference measure assumption is only partially met there.
  • ad hoc to paper The deferred finite-density prescription H_eff(µ) = H_eff(0) − µ N_eff defines a reliable effective model.
    Eqs. (3.16)–(3.21): exact only for direct Galerkin compression with N-invariant trial space; in the transfer-log route it is only an effective generator, as the paper concedes in §5. This is the central approximation, not derived from a prior theorem.

pith-pipeline@v1.3.0-alltime-deepseek · 14712 in / 11605 out tokens · 125660 ms · 2026-08-01T21:48:34.798442+00:00 · methodology

0 comments
read the original abstract

Euclidean Monte Carlo methods are effective when the path-integral weight is real and nonnegative, but finite-density fermion systems often produce sign-changing or complex scalar weights after the fermionic sector is traced out. Hamiltonian formulations avoid this complex-weight sampling problem but face rapid Hilbert-space growth. This paper studies a conditional Euclidean-Hamiltonian (CEH) reduction that combines these two descriptions. The calculation is organized around a Monte Carlo-tractable reference problem and a residual active sector. Instead of tracing the active sector into a determinant or scalar weight, CEH keeps it operator-valued and uses the reference calculation to determine projected correlation or transfer matrices. These matrices define a finite effective Hamiltonian, with the remaining finite-density dependence introduced after projection. At finite rank, the result is an effective model whose accuracy must be tested through basis enlargement, metric conditioning, stochastic matrix-element errors, and number-sector diagnostics. The construction is examined in three finite benchmarks. A two-channel oscillator tests conditional basis compression; a positive-measure stochastic calculation tests correlation-matrix harvesting and GEVP extraction, including a comparison with a PDMS-style estimator; and a four-site Hubbard ring combines a finite-budget determinant-sign stress test with finite-density continuation in non-target active spaces. The benchmarks support the proposed trade from scalar sign reweighting to a monitored active-space approximation in structured finite models, but they do not provide an end-to-end stochastic CEH treatment of the Hubbard sign problem or establish favorable scaling. Usefulness requires both low-rank approximability and efficient extraction of the required projected matrix data.

Figures

Figures reproduced from arXiv: 2607.15950 by Tsogtgerel Gantumur.

Figure 1
Figure 1. Figure 1: Low-energy errors versus R per channel. For the conditional ground state, the residual ∥HΨR − ERΨR∥ decreases from approximately 2.00 × 10−1 at R = 2 to 1.81 × 10−2 at R = 16. Thus the experiment demonstrates a low-rank basis-compression effect within ordinary Rayleigh–Ritz approximation. The gain is largest when the active labels select substantially different reference-sector geometries. 4.2. Positive-me… view at source ↗
Figure 2
Figure 2. Figure 2: Stochastic GEVP energies versus sample count. 4.2.3. Comparison with a PDMS-style estimator. We also compare the CEH estimator with a PDMS-style construction on the same two-channel Hamiltonian. The PDMS-style branch estimates projected thermal density and energy matrices, Zac(β) = ⟨ψa, e −βHψc⟩, Eac(β) = ⟨ψa, He −βHψc⟩, (4.11) and extracts projected energies from the corresponding generalized matrix probl… view at source ↗
Figure 3
Figure 3. Figure 3 [PITH_FULL_IMAGE:figures/full_fig_p019_3.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

15 extracted references · 2 linked inside Pith

  1. [1]

    Nagata, Finite-density lattice QCD and sign problem: current status and open problems,Prog

    K. Nagata, Finite-density lattice QCD and sign problem: current status and open problems,Prog. Part. Nucl. Phys.127:103991, 2022

    2022

  2. [2]

    Alexandru, G

    A. Alexandru, G. Ba¸ sar, P. F. Bedaque, and N. C. Warrington, Complex paths around the sign problem, Rev. Mod. Phys.94:015006, 2022

    2022

  3. [3]

    Aarts and D

    G. Aarts and D. Sexty, Lattice field theories with a sign problem, arXiv:2604.24290 [hep-lat], 2026

    Pith/arXiv arXiv 2026

  4. [4]

    C. W. Bauer et al., Quantum simulation for high-energy physics,PRX Quantum4:027001, 2023

    2023

  5. [5]

    Hersh and G

    R. Hersh and G. Papanicolaou, Non-commuting random evolutions, and an operator-valued Feynman-Kac formula,Comm. Pure Appl. Math.25:337–367, 1972

    1972

  6. [6]

    Brydges and P

    D. Brydges and P. Federbush, A semi-Euclidean approach to boson-fermion model theories,J. Math. Phys. 15:730–732, 1974

    1974

  7. [7]

    L¨ uscher and U

    M. L¨ uscher and U. Wolff, How to calculate the elastic scattering matrix in two-dimensional quantum field theories by numerical simulation,Nucl. Phys. B339:222–252, 1990

    1990

  8. [8]

    Blossier, M

    ALPHA collaboration, B. Blossier, M. Della Morte, G. von Hippel, T. Mendes, and R. Sommer, On the generalized eigenvalue method for energies and matrix elements in lattice field theory,JHEP04:094, 2009

    2009

  9. [9]

    Jirari, H

    H. Jirari, H. Kr¨ oger, X. Q. Luo, and K. J. M. Moriarty, Monte Carlo Hamiltonian,Phys. Lett. A258:6–14, 1999

    1999

  10. [10]

    Karna, H

    A. Karna, H. S. Wu, S. Chandrasekharan, and R. K. Kaul, Projected density matrix sampling for lattice Hamiltonians,Phys. Rev. D113:074511, 2026

    2026

  11. [11]

    F. L. Temmen, M. Gisti, D. J. Luitz, T. Luu, and J. Ostmeyer, Hamiltonian Monte Carlo enhanced by exact diagonalization, arXiv:2603.17788 [cond-mat.str-el], 2026

    arXiv 2026

  12. [12]

    Bollmark, S

    G. Bollmark, S. Mardazad, J. S. Hofmann, and A. Kantian, Fusing matrix-product states with quantum Monte Carlo: reducing entanglement and sign problem at the same time, arXiv:2411.00480 [cond-mat.str-el], 2024

    Pith/arXiv arXiv 2024

  13. [13]

    Z. Wang, Z. Wang, B.-B. Mao, and Z. Yan, Generalized reduced-density-matrix quantum Monte Carlo gives access to more, arXiv:2603.10948 [cond-mat.str-el], 2026

    arXiv 2026

  14. [14]

    Hicks and D

    C. Hicks and D. Lee, Trimmed sampling algorithm for the noisy generalized eigenvalue problem,Phys. Rev. Res.5:L022001, 2023

    2023

  15. [15]

    F. F. Assaad and H. G. Evertz, World-line and determinantal quantum Monte Carlo methods for spins, phonons and electrons, inComputational Many-Particle Physics, Lecture Notes in Physics739, pp. 277–356, Springer, 2008. McGill University, Montr ´eal, Canada National University of Mongolia, Ulaanbaatar, Mongolia Mongolian Academy of Sciences, Institute of M...

    2008