{"id":"5b62a6b8-021c-4c24-8924-3131fdb6d5ba","arxiv_id":"2507.07851","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Two Kerr solitons in a cavity drift according to exponentially weak interactions that are either purely repulsive or oscillatory, with stable molecule positions only in the oscillatory-tail case.","lead":"The paper derives equations that describe, in the large-separation limit, how two optical solitons in a microresonator move toward or away from each other. The interaction is always exponentially weak, and it is either purely repulsive or oscillates between attraction and repulsion, which decides whether soliton molecules can form.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified. The asymptotic derivation is internally consistent; the main gap is reproducibility of the numerical verification, not a demonstrated correctness flaw.","rationale":"The reader's weakest assumption concerns the orderings of Φ and ∂Φ/∂t. I find that the gauge condition (11) actually forces the projection of ∂Φ/∂t onto the adjoint zero mode to be O(ε²), removing part of the concern. The remaining unproven part is the magnitude and localization of Φ, which is standard in the pulse-interaction literature and supported by the direct numerical simulations. The larger obstacle to full acceptance is reproducibility: no code, error bars, or coefficient tables are provided, so the theoretical curves in Figures 1-2 cannot be independently recomputed. I therefore recommend no change to the reader's CONDITIONAL verdict.","tokens_in":8988,"tokens_out":30689,"duration_ms":331383,"concrete_test":"For a numerically integrated two-soliton solution at several separations L=2ξ with the paper's parameters, compute the residual Φ by projecting the field onto the two adjoint zero modes via Eq. (11); verify that |Φ| decays like e^{-σr L} and that the projected contribution ⟨Ψ̄ξ1, ∂Φ/∂t⟩ in Eq. (14) is of order e^{-2σr L} relative to the leading drift term. If both scalings hold, the truncation in Section II.B is justified; if either fails, an O(ε) term was missed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"No load-bearing flaw found in the central argument. The two-soliton ansatz (10) with gauge conditions (11) is standard, and differentiating ⟨Ψ̄ξ1, Φ⟩=0 gives ⟨Ψ̄ξ1, ∂Φ/∂t⟩ = (dξ1/dt)⟨∂Ψ̄/∂x, Φ⟩ = O(ε²) once |Φ|=O(ε) and dξ1/dt=O(ε) are granted, so the neglected term in (14) is controlled. The expansion leading to Eq. (16) is exact: the only surviving cross terms contain at least one factor of each subtracted soliton, and all linear-in-Φ terms projected on the adjoint zero mode cancel or are O(ε²). The tail replacements in Eqs. (18)-(24) are justified because the integrals are dominated by x near the soliton center, where the neighbor is represented by its far tail. The soft spot is that the bound |Φ|=O(ε) with uniform localization is asserted in Section II.B rather than proved; this is a standard result for dissipative solitons with a spectral gap, and the numerical agreement in Figures 1-2 gives empirical support, but it is not exhibited rigorously. This is an omitted-proof gap, not a demonstrated error.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives effective equations of motion for the positions of two well-separated Kerr solitons in the Lugiato-Lefever equation. Using a two-soliton ansatz, a projection onto the adjoint translation zero mode, and the exponential tail asymptotics of the single soliton, it obtains explicit drift formulas: dξ/dt = b e^{-2σξ} for monotone tails and dξ/dt = B e^{-2σr ξ} sin(2σi ξ − β) for oscillatory tails, with b, B, β expressed as single-soliton integrals. The formulas are compared with drift velocities extracted from direct numerical solutions of Eq. (1) for several parameter sets.","tokens_in":9231,"tokens_out":16829,"duration_ms":179482,"significance":"If the result holds, the paper reduces the interaction of a Kerr soliton pair to a one-dimensional ODE with no fitted parameters, and it provides a concrete explanation of the repulsive/oscillatory interaction dichotomy and of soliton-molecule formation timescales. The main strengths are the absence of free parameters, the fact that all coefficients are computed from single-soliton data, the external validation against direct PDE simulation, and the falsifiable predictions about molecule formation. The method itself is not new, but its application to Kerr solitons with both monotone and oscillatory tails is a useful and timely contribution to the microresonator literature.","major_comments":[{"comment":"The derivation rests on the ordering assumptions Φ=O(ε), ∂Φ/∂t=O(ε^2), and dξ_i/dt=O(ε), but these are asserted rather than proved. In particular, the paper does not justify that the residual field Φ remains small and localized on the drift time scale, nor that the projection onto the adjoint zero mode captures all leading-order coupling. Since Eq. (15) and the cancellation of the linear-in-Φ terms in the class-2 expansion depend on these assumptions, the claim that the equation of motion is 'asymptotically exact' should either be accompanied by a precise set of hypotheses (with a reference to a rigorous pulse-interaction result) or be softened to 'formally asymptotically exact under the stated ordering assumptions.' This is a load-bearing gap in the presentation, though not a demonstrated error.","section":"Section II.B, Eqs. (10)-(15)"},{"comment":"The abstract states that the interaction is 'either purely repulsive, or alternates between attraction and repulsion, according to whether the decay of soliton tails is monotone or oscillatory.' However, Eq. (18) with b defined by Eq. (19) permits b<0, which would give pure attraction for monotone tails. The authors only report that b is positive 'in all the cases that we studied' (Section III.A) and later state that they 'have not been able to rule out this possibility' (Section IV). The abstract therefore overstates the proven classification; it should be qualified to the parameter regimes studied, or the paper should supply an argument that b>0 for all monotone-tail solitons. This does not invalidate the drift formulas, but it is a central qualitative claim and should be corrected.","section":"Abstract and Section III.A, Eqs. (18)-(19)"}],"minor_comments":[{"comment":"Equations (7) and (8) are identical; one of the two displays should be deleted.","section":"Section II.A"},{"comment":"The first line of Eq. (16) uses Ψs2 without a tilde while the second line uses ψ̃s2; since the ansatz is expressed in terms of subtracted solitons, the notation should be made consistent.","section":"Section II.B, Eq. (16)"},{"comment":"The sentence before Eq. (17) reads 'assumeξ1 = ξ = −ξ2'; with this literal choice the argument of ψ̃s in Eq. (17) should be x+ξ, not x+2ξ. The formulas are consistent if the origin is shifted to the first soliton, so that the second soliton is at −2ξ, and the text should state this shift explicitly.","section":"Section III, Eq. (17)"},{"comment":"The numerical validation does not state the discretization, the time-integration scheme, or the procedure used to extract soliton trajectories from the PDE solutions. These details should be provided, at least in an appendix or supplement, for reproducibility.","section":"Section III, Figs. 1-2"},{"comment":"Reference [16] has a formatting error: 'J. Carr and R. L. Pego, , Communications...' should include the title of the work.","section":"References"},{"comment":"The fixed-point stability statement after Eq. (24) relies on taking B as a positive modulus in the polar representation; the text should state this explicitly, since the assignment of stability to even/odd n would reverse if B were defined with the opposite sign.","section":"Section III.B, Eq. (24)"}],"recommendation":"major_revision","confidential_remarks":"This is a competent and useful application of the standard projection method. The derivation is internally consistent, and the numerical comparisons are convincing. The main issues are the unproved residual smallness assumption and an overstrong qualitative claim in the abstract; both are fixable within the scope of a revision, and I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing you should know: this is the first asymptotically exact derivation of two-soliton drift laws in the basic Lugiato-Lefever model, and I think the central argument holds. The explicit formulas (18) and (24), with coefficients defined by single-soliton integrals and no fitted parameters, are a real step beyond the variational treatments of [27,28]. The monotone-versus-oscillatory tail classification is clean, and the numerics in Figures 1-2 match the theory in the large-separation regime.\n\nThe paper does well to carry the residual field Phi explicitly and to state the orderings rather than hand-wave. The cancellation argument in Sec II.B is terse but correct in spirit: the O(1) terms cancel by N[Psi_s]=N[Psi_c]=0, the linear-in-Phi terms project to zero through the adjoint zero mode, and the surviving cross terms are those in Eq. (16). The stress-test note confirms that the O(epsilon^2) control on the time-derivative term follows once the stated orderings are granted.\n\nWhere are the soft spots? The main one is reproducibility: no code or data is shipped, and the numerical comparisons have no error bars. For a paper claiming asymptotic exactness, a reader cannot tell from the figures how large the separation must be for the theory to become quantitative. The orderings |Phi|=O(epsilon) and dPhi/dt=O(epsilon^2) are asserted, not proved; that is a standard gap in this literature, not a demonstrated flaw, but a rigorous referee could ask for a sharper statement. There is also a strong cancellation issue that the authors themselves note in the monotone-tail case: the integrand in (19) is odd-times-even, so the result is a small difference of large contributions, meaning the practical validity may start at quite large separations. It would help to quantify when the asymptotic law actually kicks in.\n\nThe citation pattern looks fair. Prior works [27,28] are variational and [29] adds high-order dispersion, so the novelty claim is accurate. Self-citations to [22,23] are for a related entrainment problem, not for the central result. There is a minor duplication of Eqs. (7)-(8) in the text, and the figures with multiplicative compensation factors are a bit hard to read, but these are cosmetic.\n\nWho is this for? Anyone modeling multisoliton states in Kerr microresonators, especially soliton molecules and crystals. The exponential slowing of binding is a useful caution for experiments. It deserves a serious referee. The math is internally consistent, the result is genuinely predictive, and the empirical check is there even if it needs hardening. I would send it to review, with a request for code/data and error bars, and I would probably cite it if I worked on soliton interactions.","headline":"A sound, parameter-free asymptotic derivation of Kerr soliton pair drift, with clean classification and credible numerics; the main gaps are reproducibility and asserted orderings, not correctness.","tokens_in":9757,"tokens_out":2277,"would_cite":true,"duration_ms":24457,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q51","35Q55","37K40","78A60"],"pacs":["42.65.Tg","05.45.Yv"],"model":"deepseek-v4-flash","headline":"This paper derives asymptotically exact equations for the drift velocity of two Kerr cavity solitons in the Lugiato-Lefever model, showing that tail shape alone decides whether the pair repels or forms a molecule.","keywords":["cavity solitons","Lugiato-Lefever equation","soliton interactions","drift velocity","soliton molecules","translational zero mode","tail asymptotics","Kerr microresonators"],"falsifier":"Run a direct numerical simulation of the Lugiato-Lefever equation for a parameter set with oscillatory tails, extract the two soliton positions over time for several large initial separations, and compare the measured drift velocity divided by $e^{-2\\sigma_r \\xi}$ with the predicted sinusoid $B\\sin(2\\sigma_i \\xi - \\beta)$; a systematic deviation that does not vanish as the separation grows would show that the leading-order projection misses a coupling.","tokens_in":8779,"feed_emoji":"🌊","tokens_out":8099,"duration_ms":81044,"temperature":0.7,"pith_summary":"This paper derives equations of motion for the positions of two bright solitons in a pumped Kerr cavity, described by the Lugiato-Lefever model, that become asymptotically exact as the soliton separation grows. The interaction arises because the exponential tails of the two solitons overlap, and this overlap couples to the global position variables associated with broken translation symmetry. For monotone tails the drift velocity is $d\\xi/dt = b e^{-2\\sigma \\xi}$ and the interaction is purely repulsive; for oscillatory tails it is $d\\xi/dt = B e^{-2\\sigma_r \\xi}\\sin(2\\sigma_i \\xi - \\beta)$, so attraction and repulsion alternate and stable fixed points correspond to soliton molecules. In both cases the coefficients are computed from single-soliton data alone. The paper checks these predictions against trajectories extracted from direct numerical solutions of the wave equation.","feed_headline":"Soliton pairs obey an exact drift law at large separation","feed_subtitle":"Tail shape decides whether two cavity solitons repel or bind, with coefficients from single-soliton data.","key_machinery":"The machinery is the tail-overlap projection method. A two-soliton waveform is written as two translated single solitons plus a constant pedestal plus a residual field $\\Phi$, with the decomposition fixed by requiring that $\\Phi$ be orthogonal to the two adjoint translational zero modes (gauge conditions (11)). The overlap parameter $\\varepsilon = e^{-\\sigma_r |\\xi_2-\\xi_1|}$ is the small parameter: the residual field is order $\\varepsilon$ and its time derivative is order $\\varepsilon^2$, so projecting the wave equation onto the adjoint zero modes leaves a single overlap integral (17). Replacing the far soliton by its tail asymptotics (3) or (4) converts that integral into the explicit drift laws (18) and (24), whose coefficients $b$, $B$, and $\\beta$ depend only on the single-soliton waveform.","core_discovery":"The central claim is that the slow drift of well-separated Kerr cavity solitons is governed by a closed, asymptotically exact effective equation for the soliton separation. Projecting the Lugiato-Lefever nonlinearity onto the adjoint zero mode of the translation symmetry, and using the known single-soliton tail asymptotics to simplify the resulting overlap integrals, the authors obtain Eq. (18) for monotone tails and Eq. (24) for oscillatory tails. The sign of the prefactor $b$ in the monotone case is positive in all studied examples, so the interaction there is purely repulsive; in the oscillatory case the prefactor $B$ multiplies a sine that changes sign periodically, producing two interlaced lattices of stable and unstable fixed points whose spacing is set by the imaginary part of the tail decay rate. Direct numerical simulation of the wave equation confirms the predicted velocity curves in the large-separation regime.","pith_inferences":["A direct extension the authors hint at but do not develop: in waveforms with more than two solitons, nearest-neighbor dominance means the two-body drift law can be summed to predict the spacing statistics of soliton crystals, giving a testable prediction for lattice uniformity.","Because the binding force decays exponentially with separation, the stable fixed points in the oscillatory-tail regime are shallow potential wells; an estimate of thermal or pump-noise-induced escape rates from these wells would tell whether the predicted high-order molecules are observable at all.","The formula implies a clean experimental signature: measured molecule separations should be spaced by approximately $\\pi/\\sigma_i$, so microresonator experiments with tunable detuning could verify the oscillatory-tail parameter directly.","One could extend the projection calculation to include weak high-order dispersion or pump modulation, following the same zero-mode projection, and check whether the drift law retains the same exponential-sinusoid form with renormalized coefficients."],"forward_implications":["Solitons with monotone tails always repel one another at large separation, so no bound two-soliton state can form in that parameter regime.","Solitons with oscillatory tails experience attraction and repulsion in alternating, equal-length windows of separation, with stable fixed points that correspond to soliton molecules.","Because the drift velocity falls off exponentially, the separation of a repelling pair grows only logarithmically in time.","Higher-order bound states are so weakly bound that in practice only tightly bound molecules form within experimentally accessible time scales.","For multi-soliton waveforms, the interaction is dominated by nearest-neighbor overlaps, so the two-body law governs the approach to uniformly spaced soliton crystals."],"supporting_citations":[{"why":"Supplies the Lugiato-Lefever equation that is the paper's model.","marker":"[26]"},{"why":"Supplies the classification of soliton tails as monotone or oscillatory and the asymptotic decay rates.","marker":"[30]"},{"why":"Provides the projection-on-adjoint-zero-mode treatment of pulse interactions treated as an integrability condition.","marker":"[17]"},{"why":"Origin of the perturbational approach for two-soliton systems that the paper generalizes.","marker":"[12]"},{"why":"Companion perturbational treatment of two-soliton systems that the paper extends to Kerr cavity solitons.","marker":"[13]"},{"why":"Generalizes the direct perturbation method to nonintegrable systems and kink and pulse interactions.","marker":"[14]"},{"why":"Supplies the normalization of the Lugiato-Lefever model used in the paper.","marker":"[23]"}],"fun_headline_variants":["Tail shape dictates soliton repulsion or binding","Exact law for slow drift of Kerr solitons","Monotone tails repel, oscillatory tails alternate","Soliton pair interaction from tail decay type","Cavity solitons: tail shape sets the force"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands on the assumption that the residual field left after subtracting two moving single solitons is as small as the tail overlap and changes slowly, so that its direct coupling to the soliton positions can be neglected at leading order.","fun_headline_variants_meta":{"raw":{"variants":["Tail shape dictates soliton repulsion or binding","Exact law for slow drift of Kerr solitons","Monotone tails repel, oscillatory tails alternate","Soliton pair interaction from tail decay type","Cavity solitons: tail shape sets the force"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000297,"raw_usage":{"total_tokens":1693,"prompt_tokens":889,"completion_tokens":804,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":727}},"tokens_in":505,"tokens_out":804,"duration_ms":8813,"temperature":1.0,"reasoning_tokens":727,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:30:38.673962+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a direct numerical simulation of the Lugiato-Lefever equation for a parameter set with oscillatory tails, extract the two soliton positions over time for several large initial separations, and compare the measured drift velocity divided by $e^{-2\\sigma_r \\xi}$ with the predicted sinusoid $B\\sin(2\\sigma_i \\xi - \\beta)$; a systematic deviation that does not vanish as the separation grows would show that the leading-order projection misses a coupling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Lugiato-Lefever equation that is the paper's model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classification of soliton tails as monotone or oscillatory and the asymptotic decay rates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the projection-on-adjoint-zero-mode treatment of pulse interactions treated as an integrability condition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Origin of the perturbational approach for two-soliton systems that the paper generalizes."},{"cited_title":"Lu, H.-J","cited_arxiv_id":null,"evidence_quote":"Companion perturbational treatment of two-soliton systems that the paper extends to Kerr cavity solitons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Generalizes the direct perturbation method to nonintegrable systems and kink and pulse interactions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the normalization of the Lugiato-Lefever model used in the paper."}],"review_version":1}