Phase-dependent kink collisions and dual critical-velocity branches in the complex sine-Gordon model

Farnaz Eizadbaksh; Mohammad Mohammadi; Vahideh Bagheri

Complex solitons flip the collision rulebook under phase tuning

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 · glm-5.2

2026-07-10 01:46 UTC pith:SMJ36WY4

load-bearing objection Dual critical-velocity branches in CSG kink collisions are a genuine new finding, but the quantitative boundaries need energy conservation and convergence checks to be trustworthy. the 4 major comments →

arxiv 2607.08752 v2 pith:SMJ36WY4 submitted 2026-07-09 hep-th nlin.PS

Phase-dependent kink collisions and dual critical-velocity branches in the complex sine-Gordon model

Mohammad Mohammadi , Farnaz Eizadbaksh , Vahideh Bagheri This is my paper

classification hep-th nlin.PS PACS 11.10.Lm05.45.Yv11.10.Kk

keywords complexdynamicscollisioncollisionsmodelphasesine-gordonvelocity

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.

The complex sine-Gordon (CSG) model extends the classic sine-Gordon equation by letting the field be complex-valued, which gives each kink soliton an internal phase angle. The paper shows through numerical simulation that when two such kinks collide, the outcome depends jointly on how fast they hit each other and on the difference between their phases. The central discovery is that there is not one threshold velocity separating bounce-apart scattering from stick-together capture, but two distinct branches. In one phase range, speeding the kinks up makes them more likely to get captured, which is the opposite of the usual intuition and of real scalar-field models where faster means more likely to scatter. In another phase range, the conventional rule returns: faster means scattering. This dual structure arises because the relative phase acts as a genuine internal degree of freedom that controls how the real and imaginary components of each kink couple during the collision, opening channels for radiation emission, bion formation, and breather-like states that have no counterpart in the real sine-Gordon system.

Core claim

The paper identifies two phase-dependent critical-velocity branches in CSG kink-kink collisions: a red branch (roughly 0.16π < θ < 0.5π) where increasing velocity above the threshold drives capture, and a blue branch (roughly 0.55π < θ < π) where increasing velocity above the threshold restores scattering. This inverts the single-threshold picture known from real scalar-field models and establishes the relative phase as a dynamical parameter that restructures the collision-outcome landscape.

What carries the argument

The complex sine-Gordon Lagrangian replaces the real scalar field with a complex one, yielding a global U(1) symmetry and an associated conserved charge. Complex kink solutions carry a constant phase θ that splits into real and imaginary sub-kink components. The initial condition is a linear superposition of two boosted complex kinks with a constant shift ensuring the field modulus interpolates between adjacent vacua. The field equation is evolved via a second-order central finite-difference scheme. Collision outcomes are classified by tracking central energy density oscillations, emitted radiative-profile energy, and extreme values of energy density, kinetic energy, gradient energy, and the

Load-bearing premise

The initial condition uses a linear superposition of two complex kink solutions minus a constant shift, which is only an approximate solution to the nonlinear field equation. The paper does not quantify the error introduced by this approximation or test how sensitive the critical-velocity boundaries are to the initial separation distance. Since the collision outcomes depend on fine velocity-phase distinctions, small uncontrolled errors in the initial condition could shift the

What would settle it

If the two critical-velocity branches shift, merge, or disappear when the initial kink separation is systematically varied (while keeping velocity and phase fixed), the dual-branch structure would be an artifact of the approximate initial condition rather than a genuine feature of CSG dynamics.

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

If this is right

If the relative phase genuinely acts as an independent control parameter for collision outcomes, then multi-kink collisions in CSG could exhibit chaotic scattering patterns governed by a multi-dimensional phase-velocity parameter space, analogous to multi-bounce resonance windows in φ⁴ but with an added continuous internal coordinate.
The diagnostic power of extreme central quantities (energy density, field modulus) at critical transitions could be exported to other complex or multi-component field theories as a cheap numerical probe of internal-mode excitation during soliton collisions.
If a collective-coordinate or variational ansatz can be constructed for the phase degree of freedom, the red and blue critical-velocity branches might be derived analytically as level curves of an effective potential, providing a semi-analytical phase-velocity phase diagram.
The delayed energy transfer from long-lived bion structures into secondary radiative pulses suggests that post-collision radiation signatures could serve as an indirect probe of the internal phase configuration, even when the phase itself is not directly observable.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

Dual critical-velocity branches in CSG kink collisions are a genuine new finding, but the quantitative boundaries need energy conservation and convergence checks to be trustworthy.

read the letter

The main thing to know: this paper reports that complex sine-Gordon kink collisions have two distinct phase-dependent critical-velocity branches — a red branch where increasing velocity promotes capture, and a blue branch where it restores scattering. That dual structure is genuinely new and has no counterpart in real scalar-field models. The systematic two-parameter scan across velocity and phase, the bion/breather classification, and the radiative-profile energy analysis are all real contributions to the CSG literature. The qualitative phenomena are clearly demonstrated through energy-density spacetime plots and the extreme-value diagnostics at the collision center are a nice touch. Credit where it's due: the authors mapped out a genuinely two-dimensional parameter space and found structure nobody had seen before in this model. The stress-test concern about energy conservation and convergence testing is the right one to press. The red branch is physically counterintuitive — higher kinetic energy driving capture — and the paper offers no analytical mechanism for it. Simulations run to t~1600 with a single fixed discretization (h=0.02, k=0.019), and the distinction between bion annihilation and breather stabilization depends on long-time behavior of the central energy density. Yet there's no plot of total energy E(t) versus time, no grid-refinement check, and no quantification of how much energy leaks through the boundaries or accumulates from numerical drift. The reader's concern about initial-condition approximation errors is valid but secondary — that's a one-time perturbation, while energy drift is systematic and could shift the critical-velocity boundaries or even manufacture spurious transitions. The authors themselves acknowledge they can't determine critical velocities in two phase intervals and that high-velocity radiative-energy measurements carry significant uncertainty, which is honest but doesn't resolve the issue. The initial-condition superposition (Eq. 25) is standard practice for well-separated kinks and probably fine at b=-a=20, but a sensitivity check would be cheap and would strengthen the claim. None of this invalidates the qualitative findings — the dual branch structure is visible across multiple independent diagnostics (extreme values, radiative energy ratios, direct spacetime plots), which makes it unlikely to be a pure numerical artifact. But the quantitative boundaries in Fig. 2 need validation before they can be trusted as physical rather than methodological. This paper is for researchers working on soliton dynamics in field theories with internal degrees of freedom. A reader who wants to understand how phase moduli affect collision outcomes will get real value from it. It deserves a serious referee who asks for energy conservation plots and at least one convergence check before accepting the quantitative claims.

Referee Report

4 major / 8 minor

Summary. This paper presents a numerical study of complex kink-kink collisions in the complex sine-Gordon (CSG) model, focusing on how the collision outcome depends jointly on the initial velocity and the relative phase between the two kinks. The CSG model has a global U(1) symmetry, and its kink solutions carry an internal phase degree of freedom. The authors find that, unlike the integrable real sine-Gordon model with its elastic collisions, the CSG system exhibits scattering, capture, bion formation, breather-like states, and radiative profile emission. The central claim is the discovery of two distinct phase-dependent branches of critical velocity: a 'red' branch (0.16π < θ < 0.5π) where velocities above v_c lead to capture, and a 'blue' branch (0.55π < θ < π) where velocities above v_c lead to scattering. The authors also analyze the energy of emitted radiative profiles, the oscillation periods of bions, and extreme values of energy density and field modulus at the collision center as diagnostics of phase-controlled dynamics.

Significance. The identification of dual phase-dependent critical-velocity branches, particularly the counterintuitive 'red' branch where increasing velocity promotes capture, is a genuinely novel result that, if correct, substantially enriches the known phenomenology of soliton collisions in (1+1)-dimensional field theories. The systematic mapping of collision outcomes across the two-parameter velocity-phase space is a valuable contribution. The use of extreme-value diagnostics at the collision center (Sec. 6) provides a practical and sensitive tool for identifying critical transitions. The study is purely numerical and does not provide machine-checked proofs or reproducible code, but the phenomena described are clearly defined and falsifiable by independent simulation.

major comments (4)

The central quantitative claim — the dual critical-velocity branches in Fig. 2 — rests entirely on simulations with a single fixed discretization (h=0.02, k=0.019, Eq. 26). No convergence test or grid-sensitivity analysis is provided. Given that the red branch is physically counterintuitive (higher velocity promoting capture) and that the distinction between bion annihilation and breather stabilization depends on long-time evolution (t~1600 in Fig. 4), it is essential to demonstrate that the critical-velocity boundaries v_c(θ) are stable under spatial and temporal refinement (e.g., halving h and k). Without this, the quantitative boundaries could be numerical artifacts. The authors should provide at least a spot-check convergence study for representative phase values in both the red and blue branches.
The paper reports no monitoring of total energy conservation. For this Hamiltonian system, E(t) = ∫T^00 dx should be constant. Energy drift from the discretization scheme could accumulate systematically over the long simulation times used (up to t~1600) and could shift critical-velocity boundaries or create spurious capture/scattering transitions. The authors should report the relative energy drift |E(t)-E(0)|/E(0) over the simulation window for representative cases, particularly for the long-lived bion simulations in Sec. 5.2.
The initial condition (Eq. 25) uses a linear superposition of two complex kink solutions minus a constant shift, which is only an approximate solution to the nonlinear field equation. The paper states this is accurate 'if the initial separation is sufficiently large' but does not quantify the error or test sensitivity to the initial separation b=-a=20. Since the collision outcomes depend on fine velocity-phase distinctions, uncontrolled initial-condition errors could shift the critical-velocity boundaries. The authors should either estimate the initial energy imbalance or demonstrate that results are insensitive to varying the initial separation.
The paper offers no physical or analytical mechanism for the red branch phenomenon (higher velocity promoting capture). While a purely numerical discovery is acceptable, some discussion of why this occurs — for instance, whether higher-velocity collisions excite internal modes that enhance energy trapping, or whether the radiative emission channel opens differently at different velocities — would strengthen the claim and help rule out numerical artifacts. The authors acknowledge in Sec. 7 that a collective-coordinate analysis might be possible; even a qualitative discussion of the mechanism would improve the paper.

minor comments (8)

Sec. 2, Eq. (10): The Noether current j^μ includes a normalization constant η, but its value is never specified. This should be stated.
Sec. 5.1: The phase ranges 0 < θ < 0.16π and 0.5π < θ < 0.55π are excluded due to 'computational constraints' with v_c tending to 0 or 1. It would help to clarify what specifically prevents numerical determination in these ranges (e.g., v_c too close to 0 or c to resolve within the simulation domain).
Fig. 2: The data points are sparse in some regions. The caption should clarify how v_c was determined at each phase (e.g., bisection in velocity with what tolerance).
Sec. 5.3: The radiative profile detection criterion |R-2π| < 0.001 is stated, but the sensitivity of the results in Figs. 7-9 to this threshold is only briefly mentioned. A brief quantitative statement (e.g., 'varying the threshold by a factor of 10 changes E_r by less than X%') would be reassuring.
Fig. 7 caption: The top and bottom panels are described as measurements at 100 and 220 time units post-collision, but the text says t=(20/v)+100 and t=(20/v)+220. The caption should match the text.
Sec. 5.2, Fig. 4: The green line representing the 'nonzero mean' is mentioned in the text but not clearly labeled in the figure. Please label it.
The term 'breather-like state' is used throughout but never precisely defined. A brief characterization (e.g., localized, periodic, non-decaying within the simulation window) would help distinguish it from a bion.
References: The self-citation [19] provides the model foundation. Consider citing any prior work on complex kink interactions or phase-dependent soliton scattering in other models for broader context.

Circularity Check

0 steps flagged

No circularity detected: numerical study with independently verifiable model and no fitted-parameter predictions

full rationale

This paper is a numerical study of complex kink collisions in the CSG model. The Lagrangian (Eq. 6), field equation (Eq. 7), and kink solutions (Eq. 14) are stated explicitly and can be independently verified. The central claims — dual critical-velocity branches, bion/breather formation, radiative profile emission — are all extracted directly from simulation outputs, not derived from an analytical formula that could reduce to its inputs. The critical velocities in Fig. 2 are read off from simulation data across the velocity-phase plane; no parameter is fitted to a subset of data and then 'predicted' on another subset. The radiative-profile energy ratio diagnostic (E_r/K in Fig. 9) is defined from first principles (energy density Eq. 9 integrated over detected radiative regions) and used as a classification tool, not as a prediction. The sole self-citation (Ref. [19], Mohammadi & Riazi 2014) provides the model definition and static solutions, which are independently reproducible from the Lagrangian; it is not invoked as a uniqueness theorem or to smuggle in an ansatz. The initial condition (Eq. 25) is a standard approximate superposition with a stated justification (large separation), not a circular definition. No step in the derivation chain reduces to its own inputs by construction. The skeptic's concerns about energy conservation monitoring and convergence testing are correctness/robustness issues, not circularity.

Axiom & Free-Parameter Ledger

6 free parameters · 4 axioms · 0 invented entities

The paper introduces no new physical entities, particles, forces, or dimensions. All solutions (complex kinks, radiative profiles, Q-balls) are derived from the standard CSG Lagrangian (Eq. 6). The 'red' and 'blue' critical-velocity branches are phenomenological classifications of simulation outcomes, not new theoretical constructs.

free parameters (6)

h (spatial step) = 0.02
Numerical discretization parameter chosen by the authors; no convergence justification provided.
k (temporal step) = 0.019
Numerical discretization parameter chosen by the authors; no stability or convergence analysis provided.
b=-a (initial kink separation) = 20
Initial separation chosen to justify the approximate superposition initial condition; sensitivity to this choice is not tested.
C (topological charge normalization) = 1/2π
Stated as a 'positive normalization constant' chosen by convention; does not affect dynamics.
η (Noether current normalization) = unspecified
Real normalization constant in the Noether current (Eq. 10-11); value not specified but does not affect collision dynamics.
Radiative profile detection threshold = |R-2π|<0.001
Ad hoc numerical criterion for identifying radiative profiles; authors state moderate variations do not qualitatively change results but do not show this.

axioms (4)

domain assumption Linear superposition of two complex kink solutions minus a constant shift provides an accurate approximate initial condition when kinks are sufficiently separated.
Stated in Sec. 4 (Eq. 25): 'if the initial separation is sufficiently large, such a linear combination, augmented by a constant shift, provides an accurate approximate initial condition.' No error bound is given.
domain assumption Second-order central finite-difference scheme with h=0.02, k=0.019 is sufficiently accurate to resolve critical-velocity boundaries.
Used throughout Sec. 4-6. No convergence test or comparison with alternative schemes is provided.
domain assumption Collision outcomes depend only on the phase difference θ=θ₂-θ₁, not on absolute phases.
Stated in Sec. 5: 'as expected from the global U(1) symmetry...the collision outcomes depend on the phase difference.' This is physically reasonable and confirmed numerically.
domain assumption The spatial domain is wide enough that boundary reflections do not affect results during the simulated time.
Stated in Sec. 4: 'the spatial domain is chosen wide enough that the boundaries remain beyond the furthest distance traveled.' Domain size is not specified.

pith-pipeline@v1.1.0-glm · 15935 in / 2812 out tokens · 210750 ms · 2026-07-10T01:46:25.815985+00:00 · methodology

0 comments

read the original abstract

The complex sine-Gordon (CSG) model contains an internal phase degree of freedom that strongly modifies the dynamics of its solitary-wave solutions. We present a numerical study of complex kink--kink collisions and determine how the final state depends jointly on the initial velocity and relative phase. In contrast with the elastic collisions of the real sine-Gordon model, the CSG system exhibits scattering, capture, long-lived bion formation, breather-like states, and emission of radiative profiles. The simulations reveal two distinct phase-dependent branches of critical velocity. In one branch, increasing the initial velocity promotes capture, whereas in the other it restores scattering. This dual structure highlights the rich velocity--phase dependence of the collision dynamics. We also compute the energy carried by radiative profiles and examine extreme values of the energy density, kinetic and gradient contributions, and field modulus at the collision center. These quantities show sharp transitions at critical points and provide sensitive diagnostics of phase-controlled dynamics. These results suggest that the relative phase behaves as an effective internal degree of freedom that plays an important role in the collision dynamics of complex solitons.