Pith. sign in

REVIEW 4 major objections 3 minor 1 cited by

A modified worm algorithm now computes the width-derivative of entanglement entropy in O(N) models at nonzero chemical potential.

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-02 20:33 UTC pith:K3UB7V6L

load-bearing objection Real algorithmic progress on boundary-deformed EE for O(N), but the sole quantitative check is an identity that would hold even for a biased simulation; needs a proper benchmark. the 4 major comments →

arxiv 2602.22881 v2 pith:K3UB7V6L submitted 2026-02-26 hep-lat hep-th

Lattice studies of entanglement entropy in O(N) models at finite densities

Aatu Rajala , Niko Jokela , Tobias Rindlisbacher This is my paper
classification hep-lat hep-th
keywords entanglement entropyO(N) modelfinite densityworm algorithmreplica trickboundary deformationlattice field theorydual variables
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.

Lattice Monte Carlo simulations of scalar field theories at finite density suffer from a sign problem; this paper bypasses it by rewriting the O(N) model in dual integer variables and sampling with a worm algorithm. The authors adapt the boundary-deformation method for entanglement entropy to this setting, introducing two new update moves — plaquette worms and defect-anti-defect worms — that keep local constraints intact when the entangling boundary is shifted. On the 3D nonlinear O(4) model, they estimate the derivative ∂ℓH2 ≈ ∂ℓS_EE and verify an exact identity linking this derivative to the charge density slope across the boundary. The reported agreement for several lattice sizes and chemical potentials leads them to conclude the algorithm is working as intended.

Core claim

The central claim is that entanglement entropy is now a computable observable in O(N) models at finite density, with the replica trick implemented through a boundary-deformation worm algorithm. The key technical step is the construction of two defect-free boundary updates: plaquette worms, which alter the flux on a temporal plaquette to make the two sides of the boundary swap compatible with charge-conservation and evenness constraints, and defect-anti-defect worms, which move a constraint violation around until it meets and annihilates its partner. Both are designed to respect detailed balance. The authors validate their implementation using the identity ∂²H2/(∂μ∂ℓ) = -4 N_t N_s^{d-1} ∂ℓ n(

What carries the argument

The central object is the ratio of replica partition functions Z(ℓ,2)/Z(ℓ+1,2), whose negative logarithm gives the lattice derivative ∂ℓ H2. Because directly sampling this ratio would require overlapping ensembles, the method connects the two partition functions by a sequence of local boundary deformations, with each step changing the temporal boundary condition at one spatial site. The two new worm moves — plaquette worms (worm heads restricted to a temporal plaquette that adjust flux variables to avoid constraint violations) and defect-anti-defect worms (moving a defect until it annihilates with its anti-defect) — make these deformations reversible. The identity (12) acts as the internal c

Load-bearing premise

The two new worm moves must preserve detailed balance and be able to reach every allowed configuration; the paper gives a plausibility argument but no formal proof, and the only numerical check is the internal identity (12).

What would settle it

Compute the ratio Z(ℓ,2)/Z(ℓ+1,2) exactly on a small lattice (for example N_t = 2, N_s = 4, O(2) or O(4) with enumerated dual configurations) and compare the exact value of ∂ℓH2 with the histogram estimate from this algorithm; a mismatch beyond statistical error would show the moves bias the sampling.

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

If this is right

  • Entanglement entropy becomes a usable observable in dual-variable simulations of O(N) models with a chemical potential, where sign problems previously blocked direct study.
  • The density of states across the entangling boundary, ∂ℓ n(ℓ,2), is measured along the way, giving a second physical quantity from the same histograms.
  • The observed turnover of ∂ℓ H2 as μ crosses its critical value suggests the entropy derivative tracks the finite-density phase transition; the authors intend to use this to extract critical exponents.
  • The identity (12) provides a built-in self-test that future boundary-deformation simulations can monitor to catch sampling biases.
  • In the large-ℓ limit, ∂ℓ H2 may be compared with thermal entropy density, extending the connection proposed for gauge theories.

Where Pith is reading between the lines

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

  • Because the only validation is an internal relation computed from the same histograms, an independent check on small lattices where the ratio Z(ℓ,2)/Z(ℓ+1,2) can be computed exactly would decisively confirm unbiasedness.
  • The plaquette-worm and defect-worm ideas are not O(N)-specific: any dual model with hard local constraints (gauge-Higgs, CP^{N-1}) faces the same defect problem when boundary conditions are swapped, so the recipe may transfer directly.
  • Monitoring the identity (12) as a function of lattice size and chemical potential could serve as a practical convergence diagnostic, flagging parameter regions where the boundary-deformation sampling becomes unreliable.
  • The histogram-based approach may suffer from slow mixing over the deformation space; measuring autocorrelations would quantify how efficiently the algorithm traverses boundary configurations.

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

4 major / 3 minor

Summary. This proceedings paper adapts the boundary-deformation method for computing derivatives of Rényi entanglement entropy from SU(N) gauge theories to O(N) models at finite chemical potential, using the dual-variable/worm formulation. The replica trick gives ∂ℓ S_EE ≈ ∂ℓ H2 = −∂ℓ log Z(ℓ,2); since changing ℓ by one has little configuration overlap, the authors deform the boundary locally and collect histograms. Two new defect-avoiding worm moves are introduced: plaquette worms that equalize temporal-link flux configurations before a boundary flip, and defect-anti-defect worms that move defects to their antiparticle and annihilate them. Results are presented for the 3D nonlinear O(4) model, with ∂ℓ H2 as a function of ℓ and μ, and an identity (12) is used as a numerical cross-check. The paper concludes that agreement with (12) 'strongly indicates that our algorithm is working as intended.'

Significance. If the proposed algorithm is unbiased, it would provide a practical route to entanglement-entropy derivatives in finite-density O(N) models, where sign-problem-free dual formulations exist, and could give new probes of phase transitions. The derivation of Eq. (12) is transparent, and the idea of treating boundary deformations as a sequence of local histogrammed updates is potentially efficient. However, the only quantitative validation reported in the paper is an internal consistency check that cannot detect a biased update. The significance of the paper therefore hinges on additional external validation or on a substantially more cautious interpretation of the reported results.

major comments (4)
  1. [§4, Eq. (12) and Fig. 6] Eq. (12) is an exact algebraic identity for any differentiable function log Z(ℓ,2): it states that ∂μ∂ℓ log Z(ℓ,2) can be computed in two ways. Both sides are measured from the same boundary-deformation histograms and the same n(ℓ,2) samples. Consequently, if the new worm moves sample a biased stationary ensemble, both sides will still agree (up to statistical noise) because the bias cancels in this internal comparison. Fig. 6 therefore does not test the correctness of the plaquette or defect-anti-defect moves. To support the central claim, provide at least one external benchmark: exact enumeration on a small lattice, comparison with an independent Metropolis/local update on a small system, or an analytic limit such as large-N or free-field behavior.
  2. [§3, Fig. 3 and Fig. 4] Detailed balance and ergodicity of the two new defect-avoidance moves are asserted only by verbal reversibility arguments. For a boundary update that consists of a sequence of several plaquette worms followed by the reverse sequence, the proposal probability of the forward sequence must equal that of the reverse sequence; this is not demonstrated when the number or order of plaquettes depends on the configuration. Likewise, the defect-anti-defect worm, which is restricted to avoid the two temporal links whose endpoints are swapped, must be shown to be reversible and irreducible on the relevant configuration space. The manuscript does not provide a proof, an exact small-system test, or a code release that would allow independent verification. This is the load-bearing step for the unbiasedness of Z(ℓ,2)/Z(ℓ+1,2).
  3. [§4, Figs. 5–6] No statistical uncertainties are shown in any figure, although the text states that jack-knife resampling is used. Without error bars, the claimed Nt and μ dependencies and the 'agreement' in Fig. 6 cannot be quantitatively assessed. In addition, Eq. (9) approximates the derivative ∂ℓ H2 by a one-lattice-spacing finite difference, and the systematic O(a) error from this approximation is not estimated or discussed. The reported plateau values and the statement that ∂ℓH2 decreases after the critical μ are therefore preliminary. Please add error bars and an estimate of the finite-difference systematic error.
  4. [§4, final paragraph, and §5] The conclusion that 'our algorithm is working as intended' is stronger than the evidence presented. The only check is an internal identity, and the visual agreement in Fig. 6 has no error bars. Given the absence of an external benchmark, the claim should either be supported by an additional validation or softened to state that the algorithm is consistent with the exact identity but that unbiasedness remains to be established. As written, the abstract and conclusion present the algorithm's correctness as established, which is not supported by the reported data.
minor comments (3)
  1. [Eq. (2)] The notation φ_x·φ_x appears without definition; presumably it denotes the O(N)-invariant norm squared, φ_x·φ_x. Please state this explicitly and check the sign of the source term (j·φ_x) against the dual formulation used in Eq. (4).
  2. [Reference [28]] The title of Ref. [28] as printed in the bibliography is duplicated: 'Worm algorithm for the Worm algorithm for the CP^{N−1} model model'. This should be corrected.
  3. [Figs. 5 and 6] Axis labels and legends are incomplete. The right panel of Fig. 5 shows 'j=0, j3=0.2' but the text specifies only j3=0.2; the curves are not labeled with the corresponding Nt values. Fig. 6 labels are similarly unclear. Please make the figures self-contained so the reader can reproduce the parameter settings.

Circularity Check

1 steps flagged

Validation via Eq. (12) is a definitional identity, so the claimed algorithm check cannot detect a biased update.

specific steps
  1. self definitional [Section 4, Eqs. (10)-(12) and Fig. 6]
    "∂²H2/(∂μ∂ℓ) = −4NtNs^{d−1} ∂n(ℓ,2)/∂ℓ.(12) ... The plots match quite well and (12) seems to hold. ... This strongly indicates that our algorithm is working as intended."

    Equation (12) is not an independent test: using (10) ∂ℓH2 = −∂ℓ logZ(ℓ,2) and (11) n = (4NtNs^{d−1})^{-1} ∂μ logZ(ℓ,2), both sides of (12) equal −∂μ∂ℓ logZ(ℓ,2) by commutation of mixed partial derivatives. Both plotted quantities are estimated from the same boundary-deformation histograms, so the agreement is automatic for whatever distribution the simulation samples, including a biased one. Therefore the agreement cannot validate the detailed-balance/ergodicity of the new plaquette and defect-anti-defect worms; it is a self-consistency check that passes by construction.

full rationale

The paper's derivation of ∂ℓH2 from histogram ratios is logically self-contained: the boundary-deformation method and the worm updates are described, and the plotted results follow from measured histograms. The main issue is the central validation step. The paper claims that agreement with Eq. (12) 'strongly indicates that our algorithm is working as intended.' But Eq. (12) is an exact identity following from Eqs. (10) and (11) for any differentiable log Z(ℓ,2), and both sides are computed from the same boundary-deformation histograms. A violation of detailed balance in the new plaquette or defect-anti-defect worms would lead to a biased estimate of log Z(ℓ,2), yet both sides of (12) would still agree because they are derived from the same biased log Z. Thus the check is tautological with respect to the algorithmic correctness claim. The paper provides a plausible verbal reversibility argument for the new moves, but no formal proof, no exact small-system benchmark, no comparison with published results, and no code release. These are correctness risks, but the circularity specifically lies in using a definitional identity as evidence for unbiasedness. The self-citations to prior boundary-deformation work [6,7] are not load-bearing in a circular way: they introduce a method that is then adapted, and the adaptation's correctness is the issue. I therefore score 6: a central validation step reduces by construction, while the raw measurements and algorithmic framework retain independent content.

Axiom & Free-Parameter Ledger

0 free parameters · 6 axioms · 0 invented entities

The central claim involves no fitted free parameters; κ=1.2 and j3=0.2 are simulation choices rather than fitted constants. No new physical entities are introduced; 'defects' are algorithmic bookkeeping for violations of local constraints.

axioms (6)
  • standard math Replica identity tr(ρ_A^r)=Z(ℓ,r)/Z^r
    Used in Eq. (7) to define Rényi entropy; standard QFT/lattice result.
  • domain assumption Finite-difference Eq. (9) approximates ∂ℓ S_EE by −log Z(ℓ+1,2)+log Z(ℓ,2)
    The UV-divergent part cancels; this lattice approximation is inherited from earlier studies.
  • domain assumption Dual-variable worm representation (4) samples the complex-action O(N) theory correctly
    Taken from refs [19-21]; the sign problem is avoided via the flux/monomer representation.
  • domain assumption Boundary-deformation sequence with histograms gives unbiased ratios of partition functions
    Inherited from [6,7] and required for Eq. (9).
  • ad hoc to paper Plaquette-worm and defect-anti-defect-worm moves preserve detailed balance and ergodicity
    Core new algorithmic claim; argued verbally but only tested via identity (12), not independently established.
  • standard math Mixed-derivative identity (12) is exact by commuting derivatives
    Follows from Eqs. (10)-(11); used for validation.

pith-pipeline@v1.3.0-alltime-deepseek · 7696 in / 13589 out tokens · 130158 ms · 2026-08-02T20:33:01.716453+00:00 · methodology

0 comments
read the original abstract

As a characteristic property of all quantum systems, entanglement participates in many important quantum phenomena. In this proceeding, we employ it in the study of quantum field theories at finite density. We incorporate evaluations of entanglement entropy using the replica trick into MC simulations of O(N) models at finite density with the worm algorithm and present some initial results for the nonlinear O(4) model in 3 dimensions.

Figures

Figures reproduced from arXiv: 2602.22881 by Aatu Rajala, Niko Jokela, Tobias Rindlisbacher.

Figure 1
Figure 1. Figure 1: An illustration of the worm update in the 𝑘-sector. The first step is the insertion of the external source-sink pair 𝜙 +𝜙 − . Then the head 𝜙 + is moved to a neighboring site 𝑥 ′ = 𝑥 + b𝜈 and 𝑘 𝑥,𝜈 is shifted by ±1 (depending on whether 𝜈 is a positive or negative direction) to compensate for the charge displacement. This is continued until the head and tail are again on the same site when they can be remo… view at source ↗
Figure 2
Figure 2. Figure 2: Illustration of the topologies of 𝑍 and 𝑍(𝑙, 𝑟 = 2). The fields are 1/𝑇 periodic in the temporal direction for 𝑍. For 𝑍(ℓ, 2), in region 𝐴 of width ℓ, the fields are 2/𝑇 periodic. In region 𝐵, they are 1/𝑇 periodic. We have chosen to set 𝑟 = 2 here to provide a more clear example of 𝑍(ℓ, 𝑟) and because it is the partition function used to approximate 𝜕ℓ𝑆EE. Then 𝜕ℓ𝑆EE becomes 𝜕𝑆EE(ℓ ′ ) 𝜕ℓ′ [PITH_FULL_IMA… view at source ↗
Figure 3
Figure 3. Figure 3: A visualization of the plaquette boundary update. Plaquette worms are worm updates that are constrained to move along a temporal plaquette. They change the values of the 𝑘𝑖 = 𝑘 𝑥𝑖−b𝑡,𝑡 and 𝜒𝑖 = 𝜒𝑥𝑖−b𝑡,𝑡 such that Δ𝑘 = 𝑘1 − 𝑘0 = 0 and Δ𝜒 = 𝜒1 − 𝜒0 = 0 (mod 2), which will allow the boundary conditions to be changed without producing defects. After the change, the same plaquettes are done in reverse to restor… view at source ↗
Figure 4
Figure 4. Figure 4: A pair of defects introduced by the change of temporal boundary conditions consists of a defect (purple disk) on sites 𝑥0 and an anti-defect (purple circle) on site 𝑥1 (middle panel). A worm (red) can move the defect from site 𝑥0 around till it eventually reaches the site 𝑥1 (right panel), after which the defects can be removed. 4. Results In the following, we report on some results from our simulations of… view at source ↗
Figure 5
Figure 5. Figure 5: Left: 𝜕ℓ𝐻2 ≈ 𝜕ℓ𝑆𝐸𝐸 as a function of ℓ for 𝑁𝑡 = {5, 8} and 𝜇 = {0.15; 0.30}. Right: 𝜕ℓ𝐻2 as a function of 𝜇 for 𝑁𝑡 = {5, 6, 7, 8, 9, 10} and ℓ = 17.5. From (8), it is apparent that 𝜕ℓ𝐻2 = −𝜕ℓ log 𝑍(ℓ, 2) . (10) Therefore, if 𝜕ℓ𝐻2 is evaluated at different 𝜇 and a numerical derivative with respect to 𝜇 is taken, estimates of the quantity −𝜕𝜇𝜕ℓ log 𝑍(ℓ, 2) are obtained. Another important quantity obtained via… view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of results for −𝜕𝜇𝜕ℓ log 𝑍(ℓ, 2) as computed from 𝜕𝜇𝜕ℓ𝐻2 and −4𝑁𝑡𝑁 𝑑−1 𝑠 𝜕ℓ𝑛 for 𝑁𝑡 = {5, 7} and ℓ = 17.5. Both cases show a clear agreement between the two evaluations. 5. Conclusion and outlook In this proceeding, we have introduced a modified boundary deformation method for measuring 𝜕ℓ𝑆EE in lattice 𝑂(𝑁) models at finite 𝜇 and presented results for the nonlinear 𝑂(4) model in 𝑑 = 3. Our resu… view at source ↗

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Thermal and chemical response from entanglement entropy

    hep-th 2026-03 conditional novelty 5.0

    The derivative of entanglement entropy with respect to region size equals the thermal entropy density, and a generalized Maxwell relation connects it to charge density — tested nonperturbatively in the 3D O(4) model.

Reference graph

Works this paper leans on

32 extracted references · 1 canonical work pages · cited by 1 Pith paper

  1. [1]

    P. V. Buividovich and M. I. Polikarpov,Numerical study of entanglement entropy in SU(2) lattice gauge theory,Nucl. Phys. B802(2008) 458 [0802.4247]

    Pith/arXiv arXiv 2008

  2. [2]

    Nakagawa, A

    Y. Nakagawa, A. Nakamura, S. Motoki and V. I. Zakharov,Entanglement entropy of SU(3) Yang-Mills theory,PoSLAT2009(2009) 188 [0911.2596]. 8 Lattice studies of entanglement entropy in𝑂(𝑁)models at finite densitiesAatu Rajala

    Pith/arXiv arXiv 2009

  3. [3]

    Nakagawa, A

    Y. Nakagawa, A. Nakamura, S. Motoki and V. I. Zakharov,Quantum entanglement in SU(3) lattice Yang-Mills theory at zero and finite temperatures,PoSLATTICE2010(2010) 281 [1104.1011]

    Pith/arXiv arXiv 2010

  4. [4]

    E. Itou, K. Nagata, Y. Nakagawa, A. Nakamura and V. I. Zakharov,Entanglement in Four-Dimensional SU(3) Gauge Theory,PTEP2016(2016) 061B01 [1512.01334]

    Pith/arXiv arXiv 2016

  5. [5]

    Rabenstein, N

    A. Rabenstein, N. Bodendorfer, P. Buividovich and A. Schäfer,Lattice study of Rényi entanglement entropy in𝑆𝑈(𝑁𝑐)lattice Yang-Mills theory with𝑁𝑐=2,3,4,Phys. Rev. D 100(2019) 034504 [1812.04279]

    Pith/arXiv arXiv 2019

  6. [6]

    Rindlisbacher, N

    T. Rindlisbacher, N. Jokela, A. Pönni, K. Rummukainen and A. Salami,Improved lattice method for determining entanglement measures in SU(N) gauge theories,PoS LATTICE2022(2022) 031 [2211.00425]

    Pith/arXiv arXiv 2022

  7. [7]

    Jokela, K

    N. Jokela, K. Rummukainen, A. Salami, A. Pönni and T. Rindlisbacher,Progress in the lattice evaluation of entanglement entropy of three-dimensional Yang-Mills theories and holographic bulk reconstruction,JHEP12(2023) 137 [2304.08949]

    Pith/arXiv arXiv 2023

  8. [8]

    Alba,Out-of-equilibrium protocol for Rényi entropies via the Jarzynski equality,Phys

    V. Alba,Out-of-equilibrium protocol for Rényi entropies via the Jarzynski equality,Phys. Rev. E95(2017) 062132 [1609.02157]

    Pith/arXiv arXiv 2017

  9. [9]

    Bulgarelli and M

    A. Bulgarelli and M. Panero,Entanglement entropy from non-equilibrium Monte Carlo simulations,JHEP06(2023) 030 [2304.03311]

    Pith/arXiv arXiv 2023

  10. [10]

    Bulgarelli and M

    A. Bulgarelli and M. Panero,Duality transformations and the entanglement entropy of gauge theories,JHEP06(2024) 041 [2404.01987]

    Pith/arXiv arXiv 2024

  11. [11]

    Bulgarelli, E

    A. Bulgarelli, E. Cellini, K. Jansen, S. Kühn, A. Nada, S. Nakajima, K. A. Nicoli et al., Flow-Based Sampling for Entanglement Entropy and the Machine Learning of Defects,Phys. Rev. Lett.134(2025) 151601 [2410.14466]

    Pith/arXiv arXiv 2025

  12. [12]

    Bulgarelli, E

    A. Bulgarelli, E. Cellini, K. Jansen, S. Kühn, A. Nada, S. Nakajima, K. A. Nicoli et al., Computing quantum entanglement with machine learning, in42th International Symposium on Lattice Field Theory, 12, 2025,2512.11389

    arXiv 2025

  13. [13]

    Amorosso, S

    R. Amorosso, S. Syritsyn and R. Venugopalan,Entanglement entropy of a color flux tube in (2+1)D Yang-Mills theory,JHEP12(2024) 177 [2410.00112]

    Pith/arXiv arXiv 2024

  14. [14]

    Amorosso, S

    R. Amorosso, S. Syritsyn and R. Venugopalan,Entanglement entropy of a color flux tube in (1+1)D Yang–Mills theory,Phys. Lett. B868(2025) 139806 [2411.12818]

    Pith/arXiv arXiv 2025

  15. [15]

    Amorosso, S

    R. Amorosso, S. Syritsyn and R. Venugopalan,Entanglement enabled tomography of flux tubes in (2+1)d yang-mills theory, 2026

    2026

  16. [16]

    Calabrese and J

    P. Calabrese and J. L. Cardy,Entanglement entropy and quantum field theory,J. Stat. Mech. 0406(2004) P06002 [hep-th/0405152]. 9 Lattice studies of entanglement entropy in𝑂(𝑁)models at finite densitiesAatu Rajala

    Pith/arXiv arXiv 2004

  17. [17]

    Calabrese and J

    P. Calabrese and J. Cardy,Entanglement entropy and conformal field theory,J. Phys. A42 (2009) 504005 [0905.4013]

    Pith/arXiv arXiv 2009

  18. [18]

    Prokof’ev and B

    N. Prokof’ev and B. Svistunov,Worm algorithms for classical statistical models,Physical review letters87(2001) 160601

    2001

  19. [19]

    Rindlisbacher and P

    T. Rindlisbacher and P. de Forcrand,Lattice simulation of the SU(2) chiral model at zero and non-zero pion density,PoSLATTICE2015(2016) 171 [1512.05684]

    Pith/arXiv arXiv 2016

  20. [20]

    Rindlisbacher, O

    T. Rindlisbacher, O. Åkerlund and P. de Forcrand,Sampling of General Correlators in Worm Algorithm-based Simulations,Nucl. Phys. B909(2016) 542 [1602.09017]

    Pith/arXiv arXiv 2016

  21. [21]

    Rindlisbacher,Lattice Field Theory at Finite Density, Ph.D

    T. Rindlisbacher,Lattice Field Theory at Finite Density, Ph.D. thesis, Zurich, ETH, 2017. 10.3929/ethz-b-000167730

    work page doi:10.3929/ethz-b-000167730 2017

  22. [22]

    Gattringer and T

    C. Gattringer and T. Kloiber,Lattice study of the Silver Blaze phenomenon for a charged scalar𝜙 4 field,Nucl. Phys. B869(2013) 56 [1206.2954]

    Pith/arXiv arXiv 2013

  23. [23]

    Gattringer and T

    C. Gattringer and T. Kloiber,Spectroscopy in finite density lattice field theory: An exploratory study in the relativistic Bose gas,Phys. Lett. B720(2013) 210 [1212.3770]

    Pith/arXiv arXiv 2013

  24. [24]

    M. G. Endres,Method for simulating O(N) lattice models at finite density,Phys. Rev. D75 (2007) 065012 [hep-lat/0610029]

    Pith/arXiv arXiv 2007

  25. [25]

    Bruckmann, C

    F. Bruckmann, C. Gattringer, T. Kloiber and T. Sulejmanpasic,Dual lattice representations for O(N) and CP(N−1) models with a chemical potential,Phys. Lett. B749(2015) 495 [1507.04253]

    Pith/arXiv arXiv 2015

  26. [26]

    Schmidt, Y

    A. Schmidt, Y. Delgado Mercado and C. Gattringer,Monte Carlo simulation of abelian gauge-Higgs lattice models using dual representation,PoSLATTICE2012(2012) 098 [1211.1573]

    Pith/arXiv arXiv 2012

  27. [27]

    Delgado Mercado, C

    Y. Delgado Mercado, C. Gattringer and A. Schmidt,Surface worm algorithm for abelian Gauge-Higgssystemsonthelattice,Comput.Phys.Commun.184(2013)1535[1211.3436]

    Pith/arXiv arXiv 2013

  28. [28]

    Rindlisbacher and P

    T. Rindlisbacher and P. de Forcrand,Worm algorithm for the Worm algorithm for the𝐶𝑃𝑁−1 model model,Nucl. Phys. B918(2017) 178 [1610.01435]

    Pith/arXiv arXiv 2017

  29. [29]

    R. H. Swendsen and J.-S. Wang,Nonuniversal critical dynamics in monte carlo simulations, Phys. Rev. Lett.58(1987) 86

    1987

  30. [30]

    Pelissetto and E

    A. Pelissetto and E. Vicari,Critical phenomena and renormalization group theory,Phys. Rept.368(2002) 549 [cond-mat/0012164]

    Pith/arXiv arXiv 2002

  31. [31]

    M. A. Metlitski, C. A. Fuertes and S. Sachdev,Entanglement Entropy in the O(N) model, Phys. Rev. B80(2009) 115122 [0904.4477]

    Pith/arXiv arXiv 2009

  32. [32]

    Goldstein and E

    M. Goldstein and E. Sela,Symmetry-resolved entanglement in many-body systems,Phys. Rev. Lett.120(2018) 200602 [1711.09418]. 10

    Pith/arXiv arXiv 2018