[Paper Review] Interpretation of topologically restricted measurements in lattice sigma-models

Irais Bautista, Wolfgang Bietenholz+4 more|arXiv (Cornell University)|Feb 11, 2014

Theoretical and Computational Physics4 citations

TL;DR

This paper presents and tests two analytical methods to extract full physical expectation values and topological susceptibility from Monte Carlo simulations trapped in single topological sectors—common in lattice field theories with slow topological transitions. Using correlation functions of topological charge density and sector-weighted fitting, the approach accurately estimates observables and χₜ even with restricted data, validated numerically in 2d O(3) σ-models with high precision.

ABSTRACT

We consider models with topological sectors, and difficulties with their Monte Carlo simulation. In particular we are concerned with the situation where a simulation has an extremely long auto-correlation time with respect to the topological charge. Then reliable numerical measurements are possible only within single topological sectors. The challenge is to assemble such restricted measurements to obtain an approximation for the full-fledged result, which corresponds to the correct sampling over the entire set of configurations. Under certain conditions this is possible, and it provides in addition an estimate for the topological susceptibility chi_t. Moreover, the evaluation of chi_t might be feasible even from data in just one topological sector, based on the correlation of the topological charge density. Here we present numerical test results for these techniques in the framework of non-linear sigma-models.

Motivation & Objective

• To address the challenge of Monte Carlo simulations being trapped in single topological sectors due to long autocorrelation times.
• To develop a method for extracting full physical expectation values ⟨O⟩ from measurements restricted to individual topological sectors.
• To enable estimation of the topological susceptibility χₜ even when full sampling over sectors is infeasible.
• To validate the accuracy of these approximations in non-linear σ-models with known exact or numerical results.

Proposed method

• Uses the topological charge density correlation function to estimate χₜ within a single topological sector, based on a formula from Ref. [6].
• Applies a sector-weighted fitting procedure (Eq. 5.1 from Ref. [14]) to combine restricted measurements across multiple volumes and sectors to reconstruct ⟨O⟩.
• Employs the geometric definition of topological charge on the lattice to ensure integer Q values for each configuration.
• Performs numerical tests in 2d O(3) and O(2) σ-models using standard actions and cluster algorithms for full sampling comparison.
• Uses least-squares fitting to estimate ⟨S⟩/V and χₜ from restricted data in volumes L = 16 to 32.
• Validates results against directly measured values from full Monte Carlo histories.

Experimental results

Research questions

• RQ1 Can physical expectation values be reliably reconstructed from Monte Carlo data restricted to a single topological sector?
• RQ2 Is the topological susceptibility χₜ estimable from topological charge density correlations within a single sector?
• RQ3 How accurate are the approximations for ⟨O⟩ and χₜ when only limited sector data is available?
• RQ4 What range of topological charges and ⟨Q²⟩ values allows reliable estimation using the proposed methods?
• RQ5 Can these methods be extended to full QCD simulations with dynamical fermions?

Key findings

• The method based on topological charge density correlations successfully estimates χₜ within a single sector, with results agreeing well with full-sampling values.
• The sector-weighted fitting procedure (Eq. 5.1) reconstructs ⟨S⟩/V with high accuracy, matching the directly measured value of 1.24008(5) at β=1 and L=32.
• Estimates of χₜ from restricted data (e.g., 0.0164(5) for L=16–32) are consistent with the full-sampling result of 0.01721(4).
• The approximations work best when ⟨Q²⟩ ≳ 1.5 and only sectors with |Q| ≤ 2 are included in the analysis.
• The approach remains robust even with modest system sizes and significant splitting between ξ₀, ξ₁, and ξ₂, indicating stability under finite-size effects.
• The results demonstrate feasibility for applications in full QCD, especially with finer lattices where topological freezing becomes more severe.