Variable-mass sine-Gordon with point defects: integrability, soliton transmission, and quasi-conservation

Technical Summary: Variable-mass sine-Gordon with point defects: integrability, soliton transmission, and quasi-conservation

Problem Statement
The sine-Gordon (SG) model is a paradigmatic integrable field theory supporting exact soliton dynamics and an infinite hierarchy of conserved charges. However, realistic physical media are often inhomogeneous, requiring models where coupling coefficients depend on space and time. The variable-mass sine-Gordon (vmSG) equation, characterized by light-cone dependent mass functions $\rho_0(\xi)$ and $\rho_1(\eta)$, admits accelerating and shape-changing solitons. While the bulk vmSG equation can be locally mapped to the constant-mass SG equation via coordinate transformation, this mapping does not trivialize the problem of point defects. A defect static in physical coordinates becomes a moving trajectory in the transformed coordinates, complicating the definition of sewing conditions and conserved charges.

The central problem addressed is the construction of integrable point defects within the vmSG framework. Specifically, the authors investigate whether the integrable structure (Lax pairs, Bäcklund transformations, and infinite conservation laws) can be preserved in the presence of inhomogeneous backgrounds and point defects. Furthermore, the paper explores the "quasi-integrability" of the system: whether deforming the integrable defect sewing conditions leads to a model where lower-order charges are asymptotically conserved (vanishing time-integrated anomalies) even if exact integrability is broken.

Methodology
The authors employ a systematic algebraic approach based on the zero-curvature representation and Bäcklund-gauge formulations:

1. Lax and Bäcklund-Gauge Formulations: Starting from the Lax pair of the vmSG equation, the authors derive type-I and type-II defect matrices ($K$-matrices). These matrices connect auxiliary linear problems on either side of the defect ($x<0$ and $x>0$).
2. Sewing Conditions: The defect matrices generate "sewing conditions" (matching conditions) at the defect location. For type-I defects, bulk fields interact directly; for type-II defects, an auxiliary localized degree of freedom is introduced.
3. Riccati Hierarchy: Using the auxiliary linear problems, the authors construct Riccati equations for auxiliary fields ($\Gamma_{12}, \Gamma_{21}$). Expanding these fields in powers of the spectral parameter $\lambda$ allows for the recursive generation of an infinite hierarchy of conserved charges, including bulk and defect contributions.
4. Deformation and Anomaly Analysis: The integrable sewing conditions are deformed by introducing parameters $\alpha$ and $\beta$. The time derivatives of the modified charges (momentum $P$ and energy $E$) are calculated to identify "anomalies" ($A_1, A_2$) that prevent exact conservation.
5. Analytical Soliton Solutions: One-soliton configurations are analyzed to determine transmission factors ($z$), spectral parameter preservation, and topological conversion (kink-to-kink vs. kink-to-antikink).
6. Numerical Simulations: A second-order finite-difference scheme is implemented on a two-half-line geometry. The bulk equations are evolved separately on each side, with the defect values updated via discretized sewing conditions. This allows for the observation of reflection, radiation, and the time-integration of anomalies in both integrable ($\alpha\beta=1$) and non-integrable ($\alpha\beta \neq 1$) regimes.

Key Contributions and Results

• Integrable Defect Construction: The authors successfully construct type-I and type-II defect matrices for the vmSG model.
  • Type-I: The defect matrix retains the algebraic form of the constant-mass case, but the Bäcklund equations are "dressed" by the background functions $\rho_0$ and $\rho_1$. Compatibility requires the background functions on both sides to satisfy $\bar{\rho}_j = \epsilon \rho_j$ (with $\epsilon = \pm 1$).
  • Type-II: The construction is more restrictive, requiring identical background profiles on both sides ($\bar{\rho}_0 = \rho_0, \bar{\rho}_1 = \rho_1$) to preserve the gauge structure and close the auxiliary field equations.
• Conserved Charges: An infinite hierarchy of charges is derived. The first non-trivial combinations define momentum ($P$) and energy ($E$) type quantities. In the homogeneous limit ($\rho_0=\rho_1=1$), these reduce to standard Noether charges. In inhomogeneous backgrounds, they are integrability-generated quantities that lack standard spacetime translation invariance.
• Soliton Transmission and Topological Conversion: An analytical transmission factor $z$ characterizes the interaction.
  • $z > 0$: Corresponds to kink-kink transmission (topological charge preserved).
  • $z < 0$: Corresponds to kink-antikink conversion (topological charge changes by 2, compensated by the defect).
  • Critical limits ($z \to 0$ or $z \to \infty$) correspond to absorption or emission configurations.
• Quasi-Conservation and Anomalies:
  • The deformation condition $\alpha\beta = 1$ is identified as the closure condition for exact integrability. When satisfied, anomalies vanish, and the deformation is absorbed into a single effective Bäcklund parameter.
  • When $\alpha\beta \neq 1$, the defect is genuinely non-integrable. However, for parity-centered soliton configurations, the first-order anomalies ($A_1, A_2$) are found to be odd functions of time (under specific background symmetries). Consequently, their time integrals vanish, leading to asymptotic conservation of the first-order charges $P$ and $E$.
  • Hierarchy Dependence: The paper demonstrates that this quasi-conservation is order-dependent. While first-order charges are asymptotically conserved, third-order anomalies generally possess non-vanishing time integrals. Thus, the system exhibits low-order quasi-conservation but does not form a full quasi-integrable hierarchy.
• Numerical Verification: Simulations confirm that in the non-integrable regime ($\alpha\beta \neq 1$), the dominant excitation remains transmitted, but weak reflected components and radiative tails are generated. The numerical data supports the analytical prediction that first-order integrated anomalies vanish while third-order ones do not.

Significance and Claims
The paper establishes a framework for analyzing soliton-defect interactions in inhomogeneous media using the variable-mass sine-Gordon model. Its primary significance lies in:

1. Defining Integrable Defects in Inhomogeneous Media: It provides the explicit algebraic machinery (Lax pairs, sewing conditions, charge hierarchies) for handling point defects in variable-mass backgrounds, showing that while the bulk equations can be mapped to constant mass, the defect theory remains non-trivial and requires specific background matching.
2. Hierarchy-Dependent Quasi-Integrability: The work elucidates how exact integrability deforms into charge-dependent quasi-conservation. It claims that while deforming defect conditions can preserve the asymptotic conservation of low-order charges (momentum and energy) for specific parity-symmetric processes, this property does not generically extend to higher-order charges.
3. Defect-Controlled Transport: The results suggest that defects can control not only the phase and position of solitons but also their topological sector (conversion) and absorption/emission dynamics.

The authors modestly conclude that while they have established a framework for defect-controlled transport and identified the conditions for low-order quasi-conservation, the construction of a fully quasi-integrable theory (where an infinite subsector of charges is asymptotically conserved) via correlated bulk-defect deformations remains an open question. The paper does not propose specific experimental setups but notes potential applications in non-uniform Josephson junctions, magnetic systems, nonlinear optics, and effective molecular/DNA models where such inhomogeneities and defects are relevant.