pith. sign in

arxiv: 2511.04597 · v1 · pith:UYG4N37Anew · submitted 2025-11-06 · ⚛️ physics.comp-ph · cond-mat.stat-mech

Combining Harmonic Sampling with the Worm Algorithm to Improve the Efficiency of Path Integral Monte Carlo

Sourav Karmakar , Sutirtha Paul , Adrian Del Maestro , Barak Hirshberg This is my paper

Pith reviewed 2026-05-22 13:01 UTC · model grok-4.3

classification ⚛️ physics.comp-ph cond-mat.stat-mech
keywords path integral Monte Carloharmonic samplingworm algorithmquantum Monte Carloacceptance ratioautocorrelation timeanharmonic potential
0
0 comments X

The pith

Harmonic PIMC generates imaginary-time paths exactly from the harmonic part of the potential and accepts or rejects them using only the anharmonic remainder.

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

The paper develops H-PIMC by splitting the potential into harmonic and anharmonic contributions so that the harmonic paths can be sampled without rejection. This targets the low acceptance rates that slow standard path integral Monte Carlo in solids and confined liquids. For weakly to moderately anharmonic cases at βℏω=16 the approach raises acceptance ratios by factors of six to sixteen and shortens autocorrelation times by factors of seven to thirty. It also reaches convergence with two to four times fewer imaginary-time slices. The same partition extends to a mixed sampler for stronger anharmonicity and combines with the worm algorithm to retain the gains for indistinguishable particles.

Core claim

The central claim is that imaginary-time paths drawn exactly from the harmonic approximation to the potential can be accepted or rejected solely according to the anharmonic correction, producing an unbiased Markov chain whose statistical efficiency exceeds that of conventional PIMC proposals. Because the harmonic contribution no longer contributes to rejections, the dominant source of wasted moves disappears in nearly harmonic regimes. When the same harmonic sampling is restricted to the vicinity of local minima or combined with the worm algorithm, the efficiency improvement persists for both distinguishable and indistinguishable particles.

What carries the argument

The partition of the potential into a harmonic term whose imaginary-time paths are sampled exactly and an anharmonic remainder used only in the Metropolis acceptance step.

If this is right

  • Acceptance ratios rise by factors of 6-16 and autocorrelation times fall by factors of 7-30 for weakly to moderately anharmonic systems at βℏω=16.
  • Convergence is reached with two- to four-fold fewer imaginary-time slices.
  • M-PIMC restricts harmonic sampling to local minima and yields similar slice counts to standard PIMC while optimizing autocorrelation time for strongly anharmonic cases.
  • The harmonic-plus-worm combination produces comparable efficiency gains for systems of indistinguishable particles.

Where Pith is reading between the lines

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

  • The same potential partition could be tested in other path-integral based samplers to check whether comparable speedups appear outside the worm-algorithm setting.
  • Fewer required time slices may make simulations of larger particle numbers feasible on existing hardware.
  • The approach may extend naturally to time-dependent driving if a suitable instantaneous harmonic reference can be identified at each step.

Load-bearing premise

The potential can be partitioned into a harmonic term whose imaginary-time paths can be sampled exactly and an anharmonic remainder whose contribution enters only the acceptance step without introducing bias.

What would settle it

Apply both standard PIMC and H-PIMC to a chain of weakly anharmonic oscillators at βℏω=16, measure the bead-move acceptance ratio, and test whether the ratio increases by a factor between six and sixteen relative to the baseline method.

Figures

Figures reproduced from arXiv: 2511.04597 by Adrian Del Maestro, Barak Hirshberg, Sourav Karmakar, Sutirtha Paul.

Figure 1
Figure 1. Figure 1: FIG. 1. A schematic of the H-PIMC (upper panel) and M [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Comparison between PIMC and H-PIMC for har [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. A comparison between PIMC and H-PIMC results is shown for a weakly anharmonic system (upper panel), a moderately [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. A comparison between PIMC and H-PIMC results for N = 2 indistinguishable particles is shown for a weakly anharmonic [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Comparison of wall-time to solution between the three [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. M-PIMC-PBC results for a sinusoidal potential [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. M-PIMC results for strong anharmonicity at [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
read the original abstract

We propose an improved Path Integral Monte Carlo (PIMC) algorithm called Harmonic PIMC (H-PIMC) and its generalization, Mixed PIMC (M-PIMC). PIMC is a powerful tool for studying quantum condensed phases. However, it often suffers from a low acceptance ratio for solids and dense confined liquids. We develop two sampling schemes especially suited for such problems by dividing the potential into its harmonic and anharmonic contributions. In H-PIMC, we generate the imaginary time paths for the harmonic part of the potential exactly and accept or reject it based on the anharmonic part. In M-PIMC, we restrict the harmonic sampling to the vicinity of local minimum and use standard PIMC otherwise, to optimize efficiency. We benchmark H-PIMC on systems with increasing anharmonicity, improving the acceptance ratio and lowering the auto-correlation time. For weakly to moderately anharmonic systems, at $\beta \hbar \omega=16$, H-PIMC improves the acceptance ratio by a factor of 6-16 and reduces the autocorrelation time by a factor of 7-30. We also find that the method requires a smaller number of imaginary time slices for convergence, which leads to another two- to four-fold acceleration. For strongly anharmonic systems, M-PIMC converges with a similar number of imaginary time slices as standard PIMC, but allows the optimization of the auto-correlation time. We extend M-PIMC to periodic systems and apply it to a sinusoidal potential. Finally, we combine H- and M-PIMC with the worm algorithm, allowing us to obtain similar efficiency gains for systems of indistinguishable particles.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 1 minor

Summary. The paper introduces Harmonic PIMC (H-PIMC) and Mixed PIMC (M-PIMC) by partitioning the potential into harmonic and anharmonic parts. In H-PIMC, imaginary-time paths are sampled exactly from the harmonic Boltzmann weight and accepted/rejected using only the anharmonic action difference via Metropolis. Benchmarks on systems with increasing anharmonicity at βℏω=16 report acceptance-ratio gains of 6-16 and autocorrelation-time reductions of 7-30, plus convergence with fewer slices (additional 2-4× speedup). M-PIMC restricts harmonic sampling near local minima for strongly anharmonic cases. The methods are extended to periodic systems and combined with the worm algorithm for indistinguishable particles.

Significance. If the efficiency claims hold, the approach offers a practical, parameter-light way to accelerate PIMC for solids and dense quantum liquids by removing harmonic Trotter error exactly and improving importance sampling. The worm-algorithm extension broadens applicability to bosonic/fermionic systems. These gains are load-bearing for the central claim of improved efficiency without bias.

major comments (3)
  1. [Methods / H-PIMC definition] The manuscript does not specify how the harmonic frequency ω is chosen or optimized for a general anharmonic potential (see abstract and § on H-PIMC construction). This choice directly controls the reported acceptance and autocorrelation improvements and must be made reproducible.
  2. [Results / benchmarks] Benchmark results (acceptance ratios 6-16, autocorrelation reductions 7-30, slice-count savings) are presented without statistical error bars or uncertainty quantification. This undermines quantitative assessment of the claimed factors.
  3. [Theory / detailed balance] While the detailed-balance argument via cancellation of harmonic factors is sketched, an explicit derivation (or reference to the exact harmonic sampler) confirming that the full distribution exp(−S_harmonic − S_anharmonic) is preserved without additional corrections should be added, e.g., as an appendix or in the theory section.
minor comments (1)
  1. [Abstract] Clarify in the abstract or introduction what physical systems correspond to the quoted βℏω=16 value and how it relates to the chosen harmonic frequency.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading, positive assessment of significance, and recommendation for minor revision. We address each major comment below and will incorporate the suggested clarifications and additions in the revised manuscript.

read point-by-point responses

  1. Referee: [Methods / H-PIMC definition] The manuscript does not specify how the harmonic frequency ω is chosen or optimized for a general anharmonic potential (see abstract and § on H-PIMC construction). This choice directly controls the reported acceptance and autocorrelation improvements and must be made reproducible.

    Authors: We agree that a clear specification of ω is needed for reproducibility. In the original manuscript the benchmarks are performed at the fixed value βℏω = 16 corresponding to the underlying harmonic oscillator frequency of the model potentials. For the revised version we will add a dedicated paragraph in the H-PIMC construction section stating that, for a general anharmonic potential, ω is obtained from the local harmonic approximation (second derivative of the potential at the minimum or a least-squares fit to the harmonic part). We will also briefly discuss how the acceptance ratio depends on the mismatch between this ω and the true local curvature. revision: yes

  2. Referee: [Results / benchmarks] Benchmark results (acceptance ratios 6-16, autocorrelation reductions 7-30, slice-count savings) are presented without statistical error bars or uncertainty quantification. This undermines quantitative assessment of the claimed factors.

    Authors: We concur that error bars would improve the quantitative presentation. In the revised manuscript we will recompute the reported acceptance ratios, autocorrelation times, and slice-count convergence metrics from multiple independent runs and include statistical uncertainties (standard errors of the mean) on all benchmark figures and tables. revision: yes

  3. Referee: [Theory / detailed balance] While the detailed-balance argument via cancellation of harmonic factors is sketched, an explicit derivation (or reference to the exact harmonic sampler) confirming that the full distribution exp(−S_harmonic − S_anharmonic) is preserved without additional corrections should be added, e.g., as an appendix or in the theory section.

    Authors: We thank the referee for this suggestion. The current sketch relies on the exact cancellation of the harmonic Boltzmann factors between the proposal and the acceptance step. In the revised manuscript we will add a short appendix that derives the detailed-balance condition explicitly for the H-PIMC update, showing that the harmonic contributions cancel identically and that the target distribution exp(−S_harmonic − S_anharmonic) is recovered without bias. We will also cite the standard exact harmonic path-integral sampler used to generate the proposals. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The H-PIMC construction generates proposals exactly from the known harmonic Boltzmann weight and applies a standard Metropolis step on the anharmonic action difference; the harmonic factors cancel by algebra in the acceptance ratio, preserving detailed balance for the full weight without any fitted parameters or redefinitions. This rests on the external fact of exact harmonic sampling in imaginary time plus the standard PIMC derivation, neither of which is supplied by the present paper. No self-citation is load-bearing for the central claim, no ansatz is smuggled, and no prediction reduces to a fit by construction. The efficiency numbers are reported benchmarks on test systems rather than derived results that collapse to inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the ability to partition any potential into harmonic plus anharmonic pieces and on the exact solvability of the harmonic oscillator in imaginary time; no new particles or forces are introduced and no parameters appear to be fitted to the target data.

axioms (2)
  • standard math The harmonic oscillator imaginary-time propagator is known exactly and can be sampled without approximation.
    Invoked when the paper states that paths for the harmonic part are generated exactly.
  • domain assumption The Metropolis acceptance step using only the anharmonic energy difference produces the correct equilibrium distribution for the full potential.
    This is the standard justification for splitting the potential; it is assumed rather than re-derived.

pith-pipeline@v0.9.0 · 5836 in / 1568 out tokens · 46270 ms · 2026-05-22T13:01:16.591296+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

64 extracted references · 64 canonical work pages · 1 internal anchor

  1. [1]

    Open works by breaking the link between two beads and sampling the position of a single new bead (called the head) which lies on the same time slice

    Open and Close These are complimentary updates which form a de- tailed balance pair. Open works by breaking the link between two beads and sampling the position of a single new bead (called the head) which lies on the same time slice. The old bead is referred to as the tail. (Histori- cally, they are known in the literature as Ira and Masha). The acceptan...

    work page

  2. [2]

    It consists of linking the head to a different worldline and creating a new head bead in the process which lies on a different world line

    Swap The Swap move is the key addition which allows us to sample permutations of identical particles efficiently. It consists of linking the head to a different worldline and creating a new head bead in the process which lies on a different world line. This allows us to sample con- figurations in different global winding sectors with only spatially local ...

    work page

  3. [3]

    The results show the energy convergence with the number of beads for both PIMC and H-PIMC atβℏω= 8

    For a harmonic potential Since H-PIMC samples the density matrix of the quan- tum harmonic oscillator exactly, it is an exact method for a harmonic trap without interactions: ˆV= 1 2 mω2ˆx2 (16) Figure 2 shows the H-PIMC results and its comparison with standard PIMC for a particle inside a harmonic trap at various temperatures (βℏω). The results show the ...

    work page

  4. [4]

    The an- harmonic potential is of the form ˆV= 1 2 mω2ˆx2 1 +c 3ˆx+c4ˆx2 ,(17) wherec 3,c 4 determine the level of anharmonicity

    For an anharmonic potential Next we apply the H-PIMC algorithm to an anhar- monic system with increasing anharmonicity. The an- harmonic potential is of the form ˆV= 1 2 mω2ˆx2 1 +c 3ˆx+c4ˆx2 ,(17) wherec 3,c 4 determine the level of anharmonicity. Here, we consider three different anharmonicity regimes: weekly, moderately and strongly anharmonic regime (...

    work page

  5. [5]

    ForN= 2indistinguishable particles H-PIMC provides similar improvements as in the sin- gle particle case also for indistinguishable particles. In Figure 4 we observe the same improvements in accep- tance ratio, autocorrelation time and convergence with the number of beads forN= 2 bosons trapped via weakly to moderately anharmonic one-dimensional po- tenti...

    work page

  6. [6]

    The speedup in acceptance ratio is defined as M-PIMC value PIMC value and the energy autocorrelation time ratio is defined as M-PIMC value PIMC value

    The M-PIMC calculations use 96 beads. The speedup in acceptance ratio is defined as M-PIMC value PIMC value and the energy autocorrelation time ratio is defined as M-PIMC value PIMC value . The shaded area indicates the optimal harmonic domain. results are expected since, for largerV 0, the worldlines withW= 0, which are the ones affected by M-PIMC- PBC, ...

    work page

  7. [7]

    Harmonic potential The harmonic potential is V(x) = 1 2 mω2x2 (A1) wherem= 1 atomic unit andℏω= 3 meV≈1.1025× 10−4 Hartree

    work page

  8. [8]

    The different anharmonic regimes considered are •Weak anharmonicity⇒c 3 = 0.0025 Bohr−1, •Moderate anharmonicity⇒c 3 = 0.0045 Bohr−1, •Strong anharmonicity⇒c 3 = 0.0055 Bohr−1

    Anharmonic potential The anharmonic potential of interest is of the form V(x) = 1 2 mω2x2 1 +c 3x+c 4x2 (A2) wherem= 1 atomic unit,ℏω= 3 meV≈1.1025×10 −4 Hartree,c 4 = 10 −5 Bohr−2. The different anharmonic regimes considered are •Weak anharmonicity⇒c 3 = 0.0025 Bohr−1, •Moderate anharmonicity⇒c 3 = 0.0045 Bohr−1, •Strong anharmonicity⇒c 3 = 0.0055 Bohr−1

    work page

  9. [9]

    1, T(x|y)/T ho(x|y) T(y|x)/T ho(y|x) e−τ Vanh(y) e−τ Vanh(x) # (B1) In case of the free particle Trotter splitting, it can be shown that A(y|x) = min

    Sinusoidal potential The sinusoidal potential is of the form V(x) =−V 0 cos 2πx L .(A3) The mass of the particle is,m= 4 amu≈7291.5 atomic unit. The length of the box is,L= 3.85 ˚A≈7.27 Bohr. The different barrier regimes are, •Low barrier⇒V 0 = 0.3 meV≈1.1025×10 −5 Hartree, •Intermediate barrier⇒V 0 = 0.6 meV≈2.205× 10−5 Hartree, •High barrier⇒V 0 = 1.0 ...

    work page

  10. [10]

    E. L. Pollock and D. M. Ceperley, Path-integral com- putation of superfluid densities, Phys. Rev. B36, 8343 (1987)

    work page 1987

  11. [11]

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

    work page 1995

  12. [12]

    Dmowski, S

    W. Dmowski, S. O. Diallo, K. Lokshin, G. Ehlers, G. Ferr´ e, J. Boronat, and T. Egami, Observation of dy- namic atom-atom correlation in liquid helium in real space, Nature Communications8, 10.1038/ncomms15294 (2017)

    work page doi:10.1038/ncomms15294 2017

  13. [13]

    D. M. Ceperley and E. Manousakis, Path integral monte carlo applications to quantum fluids in confined ge- ometries, The Journal of Chemical Physics115, 10111 (2001)

    work page 2001

  14. [14]

    Boninsegni, N

    M. Boninsegni, N. Prokof’ev, and B. Svistunov, Super- glass phase of he 4, Physical review letters96, 105301 (2006)

    work page 2006

  15. [15]

    Boninsegni, A

    M. Boninsegni, A. Kuklov, L. Pollet, N. Prokof’ev, B. Svistunov, and M. Troyer, Fate of vacancy-induced supersolidity in he 4, Physical review letters97, 080401 (2006). 14

    work page 2006

  16. [16]

    Boninsegni, A

    M. Boninsegni, A. B. Kuklov, L. Pollet, N. V. Prokof’ev, B. V. Svistunov, and M. Troyer, Luttinger liquid in the core of a screw dislocation in helium-4, Phys. Rev. Lett. 99, 035301 (2007)

    work page 2007

  17. [17]

    S¨ oyler, A

    G. S¨ oyler, A. B. Kuklov, L. Pollet, N. V. Prokof’ev, and B. V. Svistunov, Underlying mechanism for the giant iso- choric compressibility of solid 4He: Superclimb of dislo- cations, Phys. Rev. Lett.103, 175301 (2009)

    work page 2009

  18. [18]

    Rota and J

    R. Rota and J. Boronat, Path integral monte carlo cal- culation of momentum distribution in solid 4 he, Journal of low temperature physics162, 146 (2011)

    work page 2011

  19. [19]

    C. M. Herdman, A. Rommal, and A. Del Maestro, Quan- tum monte carlo measurement of the chemical potential of 4He, Phys. Rev. B89, 224502 (2014)

    work page 2014

  20. [20]

    E. J. Landinez Borda, W. Cai, and M. de Koning, Dis- location structure and mobility in hcp 4He, Phys. Rev. Lett.117, 045301 (2016)

    work page 2016

  21. [21]

    Del Maestro, M

    A. Del Maestro, M. Boninsegni, and I. Affleck, 4He lut- tinger liquid in nanopores, Phys. Rev. Lett.106, 105303 (2011)

    work page 2011

  22. [22]

    Kulchytskyy, G

    B. Kulchytskyy, G. Gervais, and A. Del Maestro, Local superfluidity at the nanoscale, Phys. Rev. B88, 064512 (2013)

    work page 2013

  23. [23]

    L. V. Marki´ c and H. R. Glyde, Superfluidity, BEC, and dimensions of liquid 4He in nanopores, Phys. Rev. B92, 064510 (2015)

    work page 2015

  24. [24]

    A. Nava, D. Giuliano, P. H. Nguyen, and M. Boninsegni, Quasi-one-dimensional 4He in nanopores, Phys. Rev. B 105, 085402 (2022)

    work page 2022

  25. [25]

    Happacher, P

    J. Happacher, P. Corboz, M. Boninsegni, and L. Pollet, Phase diagram of 4he on graphene, Phys. Rev. B87, 094514 (2013)

    work page 2013

  26. [26]

    S. W. Kim, M. M. Elsayed, N. S. Nichols, T. Lakoba, J. Vanegas, C. Wexler, V. N. Kotov, and A. Del Maestro, Atomically thin superfluid and solid phases for atoms on strained graphene, Phys. Rev. B109, 064512 (2024)

    work page 2024

  27. [27]

    Dornheim, S

    T. Dornheim, S. Groth, T. Sjostrom, F. D. Malone, W. Foulkes, and M. Bonitz, Ab initio quantum monte carlo simulation of the warm dense electron gas in the thermodynamic limit, Physical Review Letters117, 156403 (2016)

    work page 2016

  28. [28]

    Dornheim, S

    T. Dornheim, S. Groth, J. Vorberger, and M. Bonitz, Ab initio path integral monte carlo results for the dynamic structure factor of correlated electrons: From the electron liquid to warm dense matter, Physical review letters121, 255001 (2018)

    work page 2018

  29. [29]

    Dornheim, J

    T. Dornheim, J. Vorberger, and M. Bonitz, Nonlinear electronic density response in warm dense matter, Phys- ical Review Letters125, 085001 (2020)

    work page 2020

  30. [30]

    Pilati, K

    S. Pilati, K. Sakkos, J. Boronat, J. Casulleras, and S. Giorgini, Equation of state of an interacting bose gas at finite temperature: A path-integral monte carlo study, Phys. Rev. A74, 043621 (2006)

    work page 2006

  31. [31]

    Trotzky, L

    S. Trotzky, L. Pollet, F. Gerbier, U. Schnorrberger, I. Bloch, N. V. Prokof’ev, B. Svistunov, and M. Troyer, Suppression of the critical temperature for superfluid- ity near the mott transition, Nature Physics6, 998–1004 (2010)

    work page 2010

  32. [32]

    Cinti, T

    F. Cinti, T. Macr` ı, W. Lechner, G. Pupillo, and T. Pohl, Defect-induced supersolidity with soft-core bosons, Na- ture communications5, 3235 (2014)

    work page 2014

  33. [33]

    T. Dornheim, Path-integral monte carlo simulations of quantum dipole systems in traps: Superfluidity, quantum statistics, and structural properties, Physical Review A 102, 023307 (2020)

    work page 2020

  34. [34]

    Pascual and J

    G. Pascual and J. Boronat, Quasiparticle nature of the bose polaron at finite temperature, Physical review let- ters127, 205301 (2021)

    work page 2021

  35. [35]

    Ghosh, M

    R. Ghosh, M. Ciardi, R. Nath, and F. Cinti, Path integral monte carlo study of a doubly dipolar bose gas, Physical Review B110, 014513 (2024)

    work page 2024

  36. [36]

    S. Jang, S. Jang, and G. A. Voth, Applications of higher order composite factorization schemes in imag- inary time path integral simulations, The Journal of Chemical Physics115, 7832 (2001)

    work page 2001

  37. [37]

    Sprik, M

    M. Sprik, M. L. Klein, and D. Chandler, Staging: A sam- pling technique for the monte carlo evaluation of path integrals, Physical Review B31, 4234 (1985)

    work page 1985

  38. [38]

    G. E. P. Box and M. E. Muller, A note on the Generation of Random Normal Deviates, The Annals of Mathemat- ical Statistics29, 610 (1958)

    work page 1958

  39. [39]

    Boninsegni, N

    M. Boninsegni, N. V. Prokof’ev, and B. V. Svis- tunov, Worm algorithm for continuous-space path in- tegral Monte Carlo simulations, Phys. Rev. Lett.96, 070601 (2006)

    work page 2006

  40. [40]

    Boninsegni, N

    M. Boninsegni, N. V. Prokof’ev, and B. V. Svistunov, Worm algorithm and diagrammatic Monte Carlo: A new approach to continuous-space path integral Monte Carlo simulations, Phys. Rev. E74, 036701 (2006)

    work page 2006

  41. [41]

    R. E. Zillich, J. M. Mayrhofer, and S. A. Chin, Extrap- olated high-order propagators for path integral Monte Carlo simulations, The Journal of Chemical Physics132, 3297888 (2010)

    work page 2010

  42. [42]

    S. L. Mielke and D. G. Truhlar, A whole-path importance-sampling scheme for Feynman path integral calculations of absolute partition functions and free en- ergies, The Journal of Chemical Physics144, 4939869 (2016)

    work page 2016

  43. [43]

    Yilmaz, K

    A. Yilmaz, K. Hunger, T. Dornheim, S. Groth, and M. Bonitz, Restricted configuration path integral Monte Carlo, The Journal of Chemical Physics153, 0022800 (2020)

    work page 2020

  44. [44]

    Del Maestro, N

    A. Del Maestro, N. S. Nichols, T. R. Prisk, G. Warren, and P. E. Sokol, Experimental realization of one dimen- sional helium, Nat. Comm.13, 3168 (2022)

    work page 2022

  45. [45]

    Friedel oscillations in one-dimensional 4He

    B. Rosenow and A. Del Maestro, Friedel oscillations in one-dimensional 4He, arXiv:2411.13654 (2024)

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  46. [46]

    R. A. Friesner and R. M. Levy, An optimized harmonic reference system for the evaluation of discretized path in- tegrals, The Journal of Chemical Physics80, 4488 (1984)

    work page 1984

  47. [47]

    C. E. Chao and H. C. Andersen, Local parabolic refer- ence approximation of thermal Feynman path integrals in quantum Monte Carlo simulations, The Journal of Chem- ical Physics107, 10121 (1997)

    work page 1997

  48. [48]

    Robertson and S

    C. Robertson and S. Habershon, Harmonic-phase path- integral approximation of thermal quantum correlation functions, The Journal of Chemical Physics148, 5002189 (2017)

    work page 2017

  49. [49]

    S. G. Moustafa and A. J. Schultz, Harmonic Oscillator Staging Coordinates for Efficient Path Integral Simula- tions of Quantum Oscillators and Crystals, Journal of Chemical Theory and Computation , 4 (2024)

    work page 2024

  50. [50]

    R. P. Feynman and A. R. Hibbs,Quantum Mechanics and Path Integrals, emended edition, edited by daniel f. styer ed. (Dover Publications, Mineola, NY, 2010) originally published 1965 by McGraw–Hill

    work page 2010

  51. [51]

    Shao, Elementary derivation of the quantum propa- gator for the harmonic oscillator, American Journal of 15 Physics84, 770 (2016)

    J. Shao, Elementary derivation of the quantum propa- gator for the harmonic oscillator, American Journal of 15 Physics84, 770 (2016)

    work page 2016

  52. [52]

    L. F. Barrag´ an-Gil and R. Walser, Harmonic oscillator thermal density matrix: First-order differential equa- tions for the position representation, American Journal of Physics86, 22 (2018)

    work page 2018

  53. [53]

    Krauth,Statistical Mechanics: Algorithms and Com- putations, Vol

    W. Krauth,Statistical Mechanics: Algorithms and Com- putations, Vol. 13 (Oxford University Press, 2006)

    work page 2006

  54. [54]

    Spada, S

    G. Spada, S. Giorgini, and S. Pilati, Path-integral monte carlo worm algorithm for bose systems with periodic boundary conditions, Condensed Matter7, 30 (2022)

    work page 2022

  55. [55]

    Foreman-Mackey, D

    D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Good- man, emcee: The mcmc hammer, Publ. Astron. Soc. Pac. 125, 306 (2013)

    work page 2013

  56. [56]

    Feynman, A

    R. Feynman, A. Hibbs, and D. Styer,Quantum Mechan- ics and Path Integrals, Dover Books on Physics (Dover Publications, 2010)

    work page 2010

  57. [57]

    M. E. Tuckerman,Statistical mechanics: theory and molecular simulation(Oxford university press, 2023)

    work page 2023

  58. [58]

    N. S. Nichols, T. R. Prisk, G. Warren, P. Sokol, and A. Del Maestro, Dimensional reduction of helium-4 in- side argon-plated MCM-41 nanopores, Phys. Rev. B102, 144505 (2020)

    work page 2020

  59. [59]

    L. V. Marki´ c, K. Dˇ zelalija, and H. R. Glyde, Crossover from one to two dimensions in liquid 4He in a nanopore, Phys. Rev. B101, 104505 (2020)

    work page 2020

  60. [60]

    L. V. Marki´ c, H. Vrcan, Z. Zuhrianda, and H. R. Glyde, Superfluidity, Bose-Einstein condensation, and structure in one-dimensional Luttinger liquids, Phys. Rev. B97, 014513 (2018)

    work page 2018

  61. [61]

    Boninsegni, A

    M. Boninsegni, A. B. Kuklov, L. Pollet, N. V. Prokof’ev, B. V. Svistunov, and M. Troyer, Fate of vacancy-induced supersolidity in 4He, Phys. Rev. Lett.97, 080401 (2006)

    work page 2006

  62. [62]

    T. W. Whitfield and J. E. Straub, Enhanced sampling in numerical path integration: An approximation for the quantum statistical density matrix based on the nonex- tensive thermostatistics, Phys. Rev. E64, 066115 (2001)

    work page 2001

  63. [63]

    S. L. Mielke and D. G. Truhlar, Accelerating the Conver- gence and Reducing the Variance of Path Integral Calcu- lations of Quantum Mechanical Free Energies by Using Local Reference Potentials, Journal of Chemical Theory and Computation8, 1589 (2012)

    work page 2012

  64. [64]

    T. W. Whitfield and G. J. Martyna, Low variance en- ergy estimators for systems of quantum drude oscilla- tors: Treating harmonic path integrals with large sepa- rations of time scales, The Journal of chemical physics 126(2007)

    work page 2007