Entanglement asymmetry in the gapped XYZ spin-$\frac12$ chain

Technical Summary: Entanglement Asymmetry in the Gapped XYZ Spin-1/2 Chain

Problem Statement
The entanglement asymmetry, $\Delta S^{(n)}_A$, quantifies the degree to which a subsystem $A$ breaks a symmetry present in the global ground state. While analytic results for this quantity exist for free theories and, perturbatively, for the critical XXZ chain, there has been a lack of analytic characterization for interacting systems in gapped phases with explicitly broken continuous symmetries. Specifically, the behavior of the entanglement asymmetry in the gapped, $U(1)$-breaking phase of the interacting XYZ spin-1/2 chain remained unresolved. The challenge lies in determining the additive constant in the logarithmic growth of the asymmetry, which encodes the symmetry breaking, for an interacting model where standard free-fermion techniques fail.

Methodology
The authors address this problem by combining three distinct theoretical and numerical ingredients to compute the Rényi entanglement asymmetries for a large interval in the gapped phase of the XYZ chain:

1. Charged-Moment Identity and Susceptibility Relation: The authors prove a charged-moment identity for fermionic Gaussian states and injective matrix-product states (MPS) with gapped transfer matrices. This identity relates the entanglement asymmetry to the static susceptibility of the broken charge, $\chi_{zz}$. Specifically, they establish that the variance of the charge in the reduced density matrix, $V$, is directly proportional to the ordinary charge variance, $\text{Var}(Q_A)$, allowing the asymmetry to be expressed in terms of $\chi_{zz}$.
2. Sine-Gordon Sum Rules: The scaling limit of the $\gamma$-perturbed XXZ chain is identified with the sine-Gordon model. The authors derive a non-conservation sum rule to evaluate $\chi_{zz}$ using form factors of the operator $\sin(\beta\phi)$. This involves calculating contributions from the two-kink and one-breather ($B_1$) channels. These contributions provide a precise lower bound on the universal amplitude $c(\beta^2) = \chi_{zz}/M$ (where $M$ is the kink mass).
3. Bounding the Multi-Particle Tail: To account for contributions beyond the two-kink and one-breather sectors, the authors employ two moment sum rules (a static Kubo moment and an $f$-sum rule). These rules constrain the "tail" of multi-particle contributions (starting with the four-kink term) to a two-sided corridor on the attractive side of the free-fermion point and provide a lower bound on the repulsive side.
4. Exact Lattice Input and Numerical Validation: The kink mass and coupling constants are supplied by the exact Baxter–Johnson–Krinsky–McCoy (JKM) solution. These analytical inputs are combined with infinite-system density-matrix renormalization group (iDMRG) simulations. The simulations validate the master formula by measuring the susceptibility directly and by computing the asymmetry via charged moments of the MPS, confirming the theoretical predictions across various couplings and Rényi indices.

Key Contributions and Results

• Master Formula: The authors derive a master formula for the entanglement asymmetry in the gapped phase:
  $$\Delta S^{(n)}_A (\ell) = \frac{1}{2} \log (\pi \ell \chi_{zz}) + \frac{\log n}{2(n-1)} + O\left(\frac{1}{M\ell}\right)$$
  This formula links the asymptotic behavior of the asymmetry to the static susceptibility $\chi_{zz}$ and the kink mass $M$.
• Universal Amplitude Bounds: The paper provides a rigorous lower bound for the universal amplitude $c(\beta^2)$ composed of the two-kink and one-breather contributions ($c_{2k} + c_{B1}$).
• Tail Estimation: The deficit between the lower bound and the actual amplitude is attributed to multi-particle tails. The authors show that at the free-fermion point ($\beta^2 = 1/2$), the leading tail term (four-kink) vanishes quadratically, following a tangent law. Perturbative estimates and numerical bounds confirm that this term accounts for the small excess observed in simulations near the free-fermion point.
• Numerical Confirmation: iDMRG simulations, extrapolated to the $\gamma \to 0$ limit, reproduce the master formula. The measured susceptibility matches the values derived from the sum rules within systematic uncertainties. Furthermore, direct computation of the asymmetry from interval density matrices confirms the predicted subleading constant $c_1$.
• Behavior Near Free-Fermion Point: The analysis reveals that the amplitude follows a tangent law, $c(\beta^2) \approx 1 + (\frac{\pi}{2}-1)(\xi_{sG}-1)$, near the free-fermion point, with the slope determined by the cancellation between the two-kink threshold and the breather contribution.

Significance
The paper claims to provide the first analytic characterization of entanglement asymmetry for an interacting model with an explicitly broken $U(1)$ symmetry at equilibrium. By pinning the amplitude across the phase—from both sides of the free-fermion point for $\xi_{sG} < 1$ and bounding it for $\xi_{sG} > 1$—the work establishes a comprehensive understanding of how symmetry breaking manifests in entanglement measures for gapped interacting systems. The results highlight the role of the free-fermion point as a marginal case where higher-order particle channels switch off, and they offer a baseline for understanding post-quench dynamics and Mpemba-type behavior in such systems. The work bridges integrable field theory techniques (sine-Gordon form factors) with modern tensor network simulations (iDMRG/iMPS) to solve a non-trivial problem in quantum many-body physics.