Brownian Bridge for Coherent State Path Integral Monte Carlo

Riccardo Fantoni

REVIEW 3 major objections 4 minor 14 references

Brownian Bridge for Coherent State Path Integral Monte Carlo

T0 review · 3 major / 4 minor · reviewed 2026-07-13 · grok-4.5

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

Pith's one-line read A Brownian-bridge construction lets coherent-state path-integral Monte Carlo recover the thermodynamics of two-dimensional helium in agreement with the conventional plane-wave method.

desk verdict Solid technical fix for the missing Brownian-bridge move in CSPIMC, with clear tables, but continuum agreement still rests on two fitted parameters. read the letter →

arxiv 2607.08787 v1 pith:VODLBBUI submitted 2026-06-26 cond-mat.stat-mech cond-mat.mtrl-sciphysics.chem-phphysics.comp-phquant-ph

Riccardo Fantoni This is my paper · ORCID

classification cond-mat.stat-mechcond-mat.mtrl-sciphysics.chem-phphysics.comp-phquant-ph

keywords QuantumManyBodyCoherentStatesPathIntegralMonteCarloBrownianBridgeHeliumThermodynamics

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper supplies the missing sampling move that turns an earlier coherent-state path-integral Monte Carlo (CSPIMC) scheme into a practical algorithm for identical particles. The new multi-slice Brownian bridge exchanges both real particle coordinates and their associated ghost degrees of freedom, thereby generating the permutation sum required for Bose statistics. When the dimensionless parameter that links the coherent-state width to the imaginary-time step is held fixed near 0.7, the resulting kinetic and potential energies for sixteen helium atoms on a plane match those obtained from ordinary plane-wave path-integral Monte Carlo. The continuum limit is taken by simultaneously sending the time step to zero and the harmonic-oscillator stiffness that defines the coherent states to infinity, recovering the plane-wave theory. A sympathetic reader therefore obtains a numerically verified alternative representation of the thermal density matrix that can be used for the same low-temperature helium thermodynamics.

What carries the argument

The multi-slice Brownian bridge (Eqs. 4.1–4.2) that simultaneously reconstructs free-particle paths for two real particles and their four associated ghost coordinates, thereby sampling particle exchanges while leaving ghost momenta untouched.

What would settle it

Repeat the continuum-limit series of Tables VII–VIII without any fitted kinetic-energy prefactor; if the energies then fail to approach the plane-wave values, the claimed equivalence collapses.

Watch

Extended reading notes

Core claim

With the Brownian-bridge construction of Section IV and the coherent-state parameter held fixed at approximately 0.7, CSPIMC yields thermodynamic energies for N=16 two-dimensional 4He that agree with conventional plane-wave path-integral Monte Carlo; the continuum limit taken at fixed parameter recovers the plane-wave results.

Load-bearing premise

The kinetic-energy estimator of the coherent-state method must be multiplied by an adjustable factor of order 0.4 before it matches the plane-wave value, even though the bridge itself never moves the ghost momenta that enter that estimator.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

The new piece is the explicit Brownian-bridge construction (Sec. IV) that lets CSPIMC sample particle exchanges. The author already published the coherent-state algorithm; this note supplies the missing multislice move, the Metropolis ratio (4.6), and a set of tables for N=16 2-D 4He that show energies can be brought into statistical agreement with ordinary plane-wave PIMC once φ is held near 0.7.

What works: the bridge formulas are written cleanly enough to re-implement, the displace-versus-bridge consistency checks for Boltzmann statistics are sensible, and the tables (III–VIII) are internally consistent within the quoted errors. The continuum-limit discussion correctly insists that φ must stay fixed while τ o0 and ξ o∞. That is a genuine technical completion, not a restatement.

Soft spots are real but limited. The kinetic estimator (3.9) still carries an ad-hoc prefactor µ that is fitted run-by-run to force match with PWPIMC; its continuum value is still drifting (0.44 o0.30 as M goes from 250 to 1000). The bridge itself never moves the ghost momenta, so the second term in EK is untouched by the very move the paper advertises. φ is also chosen for best numerical agreement rather than derived. These are fitting steps, not contradictions; they simply mean the claimed recovery of plane-wave theory is empirical rather than parameter-free. No code or independent derivation of µ is supplied, which is a practical nuisance for anyone who wants to check the numbers.

This is for people already working with coherent-state or continuous-representation path integrals who need a working permutation sampler. It does not open a new physical question or reorganize the field, but it is formally grounded enough and the numerics are honest enough that a serious referee should see it. I would send it out for review; the referee can demand a first-principles expression for µ or a clearer statement of the residual freedom in φ. Worth a look if you care about the method; not urgent otherwise.

Referee Report

3 major / 4 minor

Summary. The manuscript supplies a concrete Brownian-bridge construction (Sec. IV, Eqs. 4.1–4.2 and Metropolis ratio 4.6) for the coherent-state path-integral Monte Carlo (CSPIMC) algorithm introduced in the author’s earlier work. The construction is used to sample particle exchanges for both Boltzmann and Bose statistics of N=16 two-dimensional 4He atoms interacting via a truncated Lennard-Jones potential. Thermodynamic energies obtained at fixed φ=ξ au/m≈0.7 are compared with conventional plane-wave PIMC (PWPIMC) at several temperatures and numbers of time slices (Tables III–VIII). The author concludes that the continuum limit τ o0 taken at fixed φ recovers the plane-wave theory and that the new bridge move correctly samples the permutation sum.

Significance. A working, permutation-capable coherent-state PIMC would enlarge the set of complete bases available for quantum many-body simulations and could open regimes inaccessible to plane-wave expansions. The paper’s explicit bridge formulas, Metropolis ratio, and side-by-side energy tables for a standard 2-D helium model constitute a concrete, falsifiable contribution that can be checked by independent implementations. The numerical agreement at fixed M=250 is useful even if the continuum recovery remains partly empirical.

major comments (3)

Eq. (3.9) and Tables V–VIII: the CSPIMC kinetic-energy estimator contains a free multiplicative factor µ that is adjusted run-by-run “in order to find agreement with the kinetic energy of the PWPIMC.” The continuum value of µ is still drifting (≈0.4 at M=250, ≈0.30 at M=1000). Without a first-principles derivation of µ the claimed continuum recovery of PWPIMC is an empirical fit rather than a controlled limit; either derive µ from the coherent-state measure or demonstrate that the unadjusted estimator converges to the same continuum value.
Sec. IV (after Eq. 4.2): the bridge move leaves the ghost momenta untouched. Because the second term of the kinetic estimator (3.9) depends only on those momenta, the permutation move that is the paper’s central technical contribution never re-samples the kinetic estimator. The manuscript must show that this choice does not bias the sampled distribution or the continuum limit, or else include the momenta in the bridge.
Tables VII–VIII: even at fixed φ=0.7 the potential energy continues to change with M and has not yet reached the PWPIMC continuum values quoted in the table captions. A quantitative extrapolation (or a clear statement that continuum agreement is only expected after µ is also adjusted) is needed before the claim “CSPIMC o PWPIMC as φ o1/√2” can be regarded as established.

minor comments (4)

The abstract and introduction state “numerically exact” results; given residual M-dependence and the fitted µ this phrasing should be softened.
Notation for the multi-index ghost variables (Q^l_{α,k}, P^l_{α,k}) is dense; a short clarifying sentence or diagram would help the reader follow the path-integral measure (3.8).
Typographical inconsistencies appear (e.g., “Universit` a”, “timeslices”, “adimensional”); a careful copy-edit is needed.
The value φ≈0.7 is presented as ≈2^{-1/2}; a brief remark on whether this choice is unique or merely convenient would be useful.

Circularity Check

2 steps flagged · score 5.0 of 10

µ in the CSPIMC kinetic estimator is fitted run-by-run to force EK agreement with PWPIMC; continuum recovery of plane-wave results therefore depends on the tuned values of µ and φ rather than a parameter-free limit.

fitted input called prediction [Eq. (3.9) and Tables V–VIII (Sec. V)]
"The parameter µ, in the CSPIMC case, is necessary in order to find agreement with the kinetic energy of the PWPIMC, as explained in Ref. [1]. … In the Table EK is taken from Table III and used to determine µ. EP is the total potential energy."

µ is solved for so that the CSPIMC kinetic-energy expression exactly reproduces the PWPIMC EK already listed in the preceding tables. Kinetic agreement is therefore forced by construction; the subsequent claim of “very good agreement” for the thermodynamics (and the continuum recovery of PWPIMC) rests on this fitted parameter rather than an independent prediction.
fitted input called prediction [Sec. V and Conclusions (Tables VII–VIII, φ = 0.7)]
"We found favorable match between PWPIMC and CSPIMC results for φ = 0.7 ≈ 2^{-1/2}. … CSPIMC φ o1/√2 −−−−−−−−−−−−−−−−→ PWPIMC"

The value φ ≈ 0.7 is selected because it produces numerical agreement with PWPIMC; the continuum-limit arrow is then written with that same fitted value. The recovery of plane-wave theory is therefore an empirical statement about the chosen fitting point, not a parameter-free derivation.

full rationale

The paper’s central numerical claim is that the new Brownian-bridge construction plus fixed φ yields thermodynamic energies that match conventional PWPIMC for 2-D 4He and that the continuum limit at fixed φ recovers the plane-wave theory. The kinetic estimator used for that comparison (Eq. 3.9) contains an extra free parameter µ that is introduced expressly “in order to find agreement with the kinetic energy of the PWPIMC” and is then determined, table by table, by setting the CSPIMC EK equal to the already-computed PWPIMC EK. Because the bridge move leaves the ghost momenta untouched, the second term of the estimator is never re-sampled by the paper’s main technical contribution. Consequently the reported kinetic agreement is true by construction once µ is chosen, and the residual EP agreement (and the statement “CSPIMC → PWPIMC as φ → 1/√2”) is obtained only after two parameters have been tuned. The bridge construction itself is independent and non-circular, so the circularity is partial rather than total; the score of 5 reflects that the quantitative continuum claim rests on fitted inputs while the algorithmic novelty does not.

Assumptions & free parameters 2 free parameters · 4 assumptions · 1 invented entities

The central numerical claim rests on two free parameters (φ and µ) that are adjusted so CSPIMC matches PWPIMC, on the standard Trotter factorization and Metropolis algorithm, and on the coherent-state representation introduced in the author’s prior work. No new physical entities are postulated; the “ghosts” are auxiliary variables of the coherent-state expansion.

free parameters (2)

φ = ξ au/m = ≈0.7 ≈ 1/√2
Dimensionless stiffness of the coherent-state oscillator; held fixed at ≈0.7 to obtain numerical agreement with plane-wave PIMC. Continuum limit requires φ fixed while au o0.
µ (kinetic-energy prefactor) = ≈0.4 (continuum)
Multiplicative factor inserted into the CSPIMC kinetic-energy estimator (Eq. 3.9) and adjusted so that EK matches the PWPIMC value; continuum value drifts toward ≈0.4.

assumptions (4)

standard math Trotter product formula for the short-time density matrix
Used to factor the thermal density matrix into M links (Sec. III).
standard math Metropolis algorithm with detailed balance samples the correct path measure
Invoked for both the displace and bridge moves (Sec. IV).
domain assumption Coherent-state resolution of the identity and the associated short-time Green function of Ref. [1]
The entire CSPIMC weight (Eqs. 3.1–3.7) is taken from the author’s previous paper without re-derivation.
ad hoc to paper Leaving ghost momenta unmoved by the bridge does not bias the kinetic-energy estimator
Stated in Sec. IV without proof; required for the claim that the bridge alone samples Bose statistics correctly.

invented entities (1)

Ghost coordinates and momenta (Q,P) attached to each real particle
purpose: Auxiliary variables that realize the coherent-state expansion of the density matrix
Introduced already in Ref. [1]; they are integration variables, not new physical degrees of freedom. No independent experimental handle is claimed.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Brownian Bridge for Coherent State Path Integral Monte Carlo." pith.science (2026). https://pith.science/paper/VODLBBUI

@misc{pith2026260708787,
  author       = {Pith},
  title        = {Pith review of: Brownian Bridge for Coherent State Path Integral Monte Carlo},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VODLBBUI}},
  note         = {Machine review of arXiv:2607.08787}
}

read the original abstract

We propose a new Brownian bridge construction for our newly devised Coherent States Path Integral Monte Carlo algorithm. We apply it to the numerically exact calculation of the thermodynamic properties of the Helium fluid on a plane at low non zero temperature. We find very good agreement with the conventional plane waves path integral Monte Carlo results.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 1 linked inside Pith

[1]

a priori

and our newly devised Coherent States PIMC (CSPIMC) algorithm. The two PIMC differ for the expression of the hot kinetic density matrix at an imaginary timestepτ=β/MwithMa large number of timeslices. The short imaginary time density matrix CSPIMC expression is reviewed in the next section. The PWPIMC requires a multidimensional integral overdN Mcoordinate...

arXiv 2026
[2]

Fantoni, Coherent State Path Integral Monte Carlo, Eur

R. Fantoni, Coherent State Path Integral Monte Carlo, Eur. Phys. J. D79, 146 (2025)

2025
[3]

J. R. Klauder, Continuous-Representation Theory. I. Postulates of Continuous-Representation Theory, J. Math. Phys.4, 1055 (1963)

1963
[4]

J. R. Klauder, Continuous-Representation Theory. II. Generalized Relation between Quantum and Classical Dynamics, J. Math. Phys.4, 1058 (1963)

1963
[5]

J. R. Klauder, ContinuousRepresentation Theory. III. On Functional Quantization of Classical Systems, J. Math. Phys.5, 177 (1964)

1964
[6]

McKenna and J

J. McKenna and J. R. Klauder, ContinuousRepresentation Theory. IV. Structure of a Class of Function Spaces Arising from Quantum Mechanics, J. Math. Phys.5, 878 (1964)

1964
[7]

J. R. Klauder and J. McKenna, ContinuousRepresentation Theory. V. Construction of a Class of Scalar Boson Field Continuous Representations, J. Math. Phys.6, 68 (1965)

1965
[8]

J. R. Klauder and B. Skagerstam,Coherent States(World Scientific, Singapore, 1985)

1985

Show all 14 references

[9]

D. M. Ceperley, Path integrals in the theory of condensed helium, Rev. Mod. Phys.67, 279 (1995)

1995
[10]

M. H. Kalos and P. A. Whitlock,Monte Carlo Methods(John Wiley & Sons Inc., New York, 1986)

1986
[11]

M. C. Gordillo and D. M. Ceperley, Path-integral calculation of the two-dimensional 4He phase diagram, Phys. Rev. B58, 6447 (1998)

1998
[12]

W. L. McMillan, Ground State of Liquid He 4, Phys. Rev.138, 442 (1964)

1964
[13]

H. F. Trotter, On the Product of Semi-Groups of Operators, Proc. Am. Math. Soc.10, 545 (1959)

1959
[14]

Metropolis, A

N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. M. Teller, and E. Teller, Equation of state calculations by fast computing machines, J. Chem. Phys.1087, 21 (1953)

1953

Pith tools

Reviewed July 13, 2026 · model on record in the stance chip above.

[1] [1]

a priori

and our newly devised Coherent States PIMC (CSPIMC) algorithm. The two PIMC differ for the expression of the hot kinetic density matrix at an imaginary timestepτ=β/MwithMa large number of timeslices. The short imaginary time density matrix CSPIMC expression is reviewed in the next section. The PWPIMC requires a multidimensional integral overdN Mcoordinate...

arXiv 2026

[2] [2]

Fantoni, Coherent State Path Integral Monte Carlo, Eur

R. Fantoni, Coherent State Path Integral Monte Carlo, Eur. Phys. J. D79, 146 (2025)

2025

[3] [3]

J. R. Klauder, Continuous-Representation Theory. I. Postulates of Continuous-Representation Theory, J. Math. Phys.4, 1055 (1963)

1963

[4] [4]

J. R. Klauder, Continuous-Representation Theory. II. Generalized Relation between Quantum and Classical Dynamics, J. Math. Phys.4, 1058 (1963)

1963

[5] [5]

J. R. Klauder, ContinuousRepresentation Theory. III. On Functional Quantization of Classical Systems, J. Math. Phys.5, 177 (1964)

1964

[6] [6]

McKenna and J

J. McKenna and J. R. Klauder, ContinuousRepresentation Theory. IV. Structure of a Class of Function Spaces Arising from Quantum Mechanics, J. Math. Phys.5, 878 (1964)

1964

[7] [7]

J. R. Klauder and J. McKenna, ContinuousRepresentation Theory. V. Construction of a Class of Scalar Boson Field Continuous Representations, J. Math. Phys.6, 68 (1965)

1965

[8] [8]

J. R. Klauder and B. Skagerstam,Coherent States(World Scientific, Singapore, 1985)

1985