{"id":"6baf2b99-5329-4bb7-9d98-7e3bde7bcdc9","arxiv_id":"2411.18600","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For fractional sine-Gordon, Caputo time order 1<β<2 damps kink motion, while Riesz space order above 2 creates non-monotonic tails and an unstable kink-antikink saddle.","lead":"This paper studies what happens to sine-Gordon kinks when the usual derivatives are replaced by fractional ones. It finds that fractional time derivatives act like friction, while non-integer spatial derivatives can flip the force between a kink and an anti-kink from attraction to repulsion.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The α>2 saddle/repulsion may be an artifact of the periodic Riesz-domain implementation; the central infinite-line claim needs a domain-size check.","rationale":"The reader correctly spotted the erroneous Mittag-Leffler identity in the Appendix (E_β(λ^2 t) should be E_β(λ^2 t^β); E_2(z)=cosh(√z), not cosh z). That error, however, is a likely typo and does not directly bear on the main phenomenological claims, which are numerical. The more load-bearing issue for the central claim is the spatial periodic setting: all Riesz computations appear to use periodic boundary conditions, and the manuscript itself attributes the final equilibration at x=±L/4 to periodic images. The K–AK saddle for α>2 is the lynchpin of the 'long-range repulsion' conclusion, so if it is a finite-size branch rather than an infinite-line equilibrium, that conclusion is not established. This does not invalidate the paper; larger-domain tests and force-curve extraction would settle it. I therefore keep the reader's conditional verdict, with the added condition of a domain-size check.","tokens_in":13100,"tokens_out":13427,"duration_ms":152123,"concrete_test":"Recompute the stationary K–AK branch and the effective force F(d) for α=2.5 and α=2.01 on domains L=100, 200, and 400 with the same Fourier/Riesz discretization. If the saddle separation d*(α) and the zero of F(d) are L-independent (or extrapolate consistently to L→∞), the repulsion claim stands. If d* saturates at L/2 as α→2+ or shifts with L, the saddle is a periodic-image artifact and the infinite-line conclusion is unsupported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that α>2 induces long-range repulsion and a saddle K–AK equilibrium is supported mainly by stationary K–AK states computed with a spatially periodic implementation of the Riesz derivative. Section 5 explicitly notes that the Riesz derivative 'naturally stems from' periodic boundary conditions and that kinks settle at x=±L/4 because each member feels forces from periodic images. On a ring, even the standard α=2 sine-Gordon equation has a K–AK equilibrium at half the circumference, so a bound-state branch whose separation tends to L/2 as α→2+ is not by itself evidence of an intrinsic infinite-line saddle. The force zero-crossing inferred from the non-monotonic kink tail may therefore be a finite-size effect unless the separation and the force curve are shown to be independent of L for the values of α used.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the fractional sine-Gordon equation (1) in three regimes: Caputo time-fractional derivative with an ordinary Laplacian, Riesz space-fractional derivative with a second time derivative, and the combined case. Using numerical simulations with the fde_pi12_pc integrator for the Caputo part and a Fourier-spectral implementation for the Riesz part, the authors report that Caputo order 1<β<2 acts as dissipation, slowing and eventually stopping kinks and dampening breathers; that Riesz order α<2 produces monotonically attracting kinks, while α>2 produces non-monotonic kink tails and an unstable kink-antikink saddle equilibrium separating attraction from repulsion; and that the combined case inherits both effects. An appendix attempts to prove that the stability spectrum of static kinks is independent of β.","tokens_in":13261,"tokens_out":9069,"duration_ms":84358,"significance":"If the central claims hold, Eq. (1) offers a tunable dissipation mechanism through the Caputo order and tunable long-range attraction or repulsion through the Riesz order, which is of genuine interest for fractional nonlinear-wave applications. The paper is honest in labeling the 1/τ scaling law as a fit rather than a predictive derivation, and the saddle equilibrium is computed numerically rather than inferred from a fitted model. The manuscript contains extensive numerical explorations, spectral stability computations, and Floquet analysis for breathers. However, the validity of the α>2 saddle-and-repulsion picture is currently weakened by the use of a periodic-domain Riesz implementation without a demonstrated infinite-line limit, and the appendix contains a false Mittag-Leffler identity that undermines the stated β-independence proof.","major_comments":[{"comment":"The separation-of-variables step uses the identity ∂_t^β E_β(λ^2 t) = λ^2 E_β(λ^2 t), which is false for the Caputo derivative. The correct Caputo eigenfunction is E_β(λ^2 t^β), satisfying ∂_t^β E_β(λ^2 t^β) = λ^2 E_β(λ^2 t^β); the identity also fails for the argument written in Eq. (.2). In addition, the statement E_2(z) = cosh(z) is incorrect; the correct relation is E_2(z) = cosh(√z). Since Section 5 explicitly invokes this appendix to conclude that the stability spectrum of stationary kinks is independent of β, the proof as written is invalid. The correction is likely local and may preserve the formal conclusion, but the manuscript must be revised to state the correct Mittag-Leffler argument and to discuss the behavior of E_β(λ^2 t^β) for λ^2 < 0 and 1 < β < 2.","section":"Appendix, Eqs. (.2)-(.4)"},{"comment":"The central claim of an intrinsic infinite-line saddle equilibrium and long-range repulsion for α > 2 rests on stationary and dynamical kink-antikink computations performed with a spatially periodic implementation of the Riesz derivative. The paper itself notes in §5.2 that kinks equilibrate at x = ±L/4 because each member feels forces from periodic images, and Fig. 6 shows that the α → 2+ bound-state separation appears to tend to L/2, which is precisely the equilibrium separation of the periodic α = 2 sine-Gordon model on a ring. The force zero-crossing inferred from non-monotonic tails may therefore be a finite-size artifact. A domain-size independence study (for example, repeating the stationary K-AK computation and the force curve at L = 200 and L = 400, or using an infinite-line formulation of the Riesz derivative) is required to support the conclusion that the saddle and the repulsive branch survive in the infinite-line limit.","section":"§4.2 and §5.2, Figs. 6, 7, and 10"},{"comment":"The quantitative claims - the exponential decay of the kink velocity, the logarithmic scaling 1/τ = a log(b|β − 2| + 1), and the α-dependence of the saddle separation - are all numerical, but no grid-convergence or time-step refinement study is reported for either the fde_pi12_pc integrator or the spectral Riesz implementation. Because the long-range tails and the bifurcation near α = 2 are delicate, a resolution study (variation of spatial grid size and time step) should be added to rule out numerical artifacts and to give the reported fit parameters a defensible accuracy.","section":"§3.1 and §4.2"}],"minor_comments":[{"comment":"The initial condition in Eq. (4) appears to have an unbalanced parenthesis: the expression for (u0, u˙0) is missing a closing parenthesis after the second component.","section":"§3.2, Eq. (4)"},{"comment":"The caption states 'β = 2.15 (right)', but this section fixes β = 2 and varies α; the right panel is presumably α = 2.15 and the caption should be corrected.","section":"Fig. 8 caption"},{"comment":"There is a typo in 'ifv0 is above a certain threshold'; a space is missing after 'if'.","section":"§3.2, text before Fig. 4"},{"comment":"Reference [4] is incomplete: the journal name and volume are missing from the citation for the review on fractional calculus in biological modeling.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The main risk to publication is the periodic-domain artifact for the α > 2 saddle-and-repulsion claim; I would ask for explicit L-independence data before acceptance. The appendix identity error is embarrassing but readily fixable, and the rest of the phenomenology appears qualitatively robust. The paper is largely numerical and phenomenological, which is appropriate for the journal, but the authors should also state the well-posedness status of Eq. (1) at least briefly, since the manuscript currently does not address it."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives a systematic numerical picture of what happens to kinks, K–AK pairs, and breathers when you replace the time or space derivative in sine-Gordon with a fractional one. The Caputo result (kinks slow exponentially and stop, with a dissipation timescale that grows logarithmically as β→2) and the Riesz result (for α>2 kinks get non-monotonic tails and a K–AK saddle appears) are specific and not in the earlier work. The numerics look reasonably careful, and the 1/τ = a·log(b|β−2|+1) law is explicitly labeled a fit, not a predictive derivation. That is honest.\n\nThe soft spots are real but uneven. The appendix is wrong: Eq. (.4) is not the Caputo identity for E_β(λ^2 t); the standard identity applies to E_β(λ t^β), and E_2(z)=cosh(√z), not cosh(z). So the claimed proof that the static-kink stability spectrum is β-independent does not hold as written. That does not kill the numerical conclusion—the dissipative picture for β<2 is plausible—but the proof needs to be fixed.\n\nThe bigger question is the α>2 saddle. The Riesz derivative is implemented periodically, and the paper itself says kinks settle at ±L/4 because of periodic images. That makes a domain-size check essential: does the K–AK equilibrium separation and the force zero-crossing stay L-independent for the α values used? The paper says they varied domain lengths and got the same outcome, but no numbers are shown. The stress-test's claim that standard α=2 already has a K–AK equilibrium at L/2 does not hold up for me—in standard sG on a ring the kink–antikink force is attractive at all distances—but the general worry about finite-size contamination is fair and should be answered with a plot or table.\n\nAlso minor: no grid-convergence or well-posedness checks, and the exponential decay fits have no error bars. For a numerical survey that is not disqualifying, but worth asking for in revision.\n\nWho is this for? People working on fractional dispersive PDEs and coherent structures. It advances the program from [14] and [20,23] in a straightforward way. I would send it to a referee, with instructions to check the appendix, ask for domain-size dependence, and request convergence data. Conditional accept after revision.","headline":"A useful numerical map of fractional sine-Gordon dynamics whose main proof appendix is wrong and whose α>2 saddle needs a domain-size check.","tokens_in":13834,"tokens_out":5865,"would_cite":true,"duration_ms":52971,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q51","35R11"],"pacs":[],"model":"deepseek-v4-flash","headline":"For the fractional sine-Gordon equation, Caputo order $\\beta<2$ damps kinks to rest, while Riesz order $\\alpha>2$ turns their tails non-monotonic and makes kink-antikink forces repulsive at long range.","keywords":["sine-Gordon equation","kinks","breathers","fractional derivatives","Caputo derivative","Riesz derivative","kink-antikink interactions","fractional Klein-Gordon equation"],"falsifier":"Compute $\\partial_t^\\beta E_\\beta(\\lambda^2 t)$ numerically for $\\beta=1.5$ and $\\lambda^2<0$ and compare it with $\\lambda^2 E_\\beta(\\lambda^2 t)$; if the two differ, the spectral-independence claim is unsupported. Alternatively, solve the linearized fractional equation for a static kink at $\\beta=1.5$ and compare its eigenvalues with the $\\beta=2$ case: any $\\beta$-dependence of the spectrum would falsify the claim.","tokens_in":12921,"feed_emoji":"🌀","tokens_out":7394,"duration_ms":58266,"temperature":0.7,"pith_summary":"This paper studies the fractional sine-Gordon equation in which the usual second time derivative is replaced by a Caputo derivative of order $\\beta$ and the spatial Laplacian by a Riesz derivative of order $\\alpha$. It argues that the two fractional orders act as independent knobs: for $\\beta<2$ the Caputo derivative injects dissipation, so single kinks slow down exponentially and stop, and kink-antikink pairs either annihilate or separate depending on whether they clear a velocity threshold. For the spatial side, it claims that $\\alpha<2$ keeps kink tails monotonic and the kink-antikink force attractive, whereas $\\alpha>2$ makes the tails cross zero, creating a saddle equilibrium and long-range repulsion with short-range attraction. When both derivatives are active, the dissipative Caputo effect dominates the decay while the Riesz order sets the interaction landscape, producing breathers that decay and kink-antikink pairs that can settle into permanent separated states. The stakes are that a canonical soliton model gains tunable dissipation and tunable long-range interaction, directly relevant to experimental efforts to engineer dispersion in optical media.","feed_headline":"Time-fractional order damps kinks; space-fractional order repels them","feed_subtitle":"In one model, the two fractional orders separately control energy loss and long-range attraction or repulsion of solitons.","key_machinery":"The machinery is the pair of fractional operators in Eq. (1): the Caputo time derivative $\\partial_t^\\beta$, whose memory kernel produces effective damping for $1<\\beta<2$, and the Riesz space derivative $\\partial_x^\\alpha$, whose power-law kernel changes the decay and oscillation of kink tails with $\\alpha$. The paper also relies on a Mittag-Leffler separation-of-variables argument in the Appendix, intended to show that the static-kink stability spectrum is independent of $\\beta$, and on the tail-mediated-force picture, in which the interaction between two solitary waves is carried by their tails, to link the zero crossing of the kink tail to the sign change of the kink-antikink force.","core_discovery":"The central claim is that the fractional sine-Gordon equation (1) has a two-parameter phenomenology organized by the orders $\\beta$ and $\\alpha$. For $\\beta<2$, the Caputo time derivative acts like a damping term: kinks launched at finite velocity decelerate exponentially and come to rest; stationary kinks shed radiation that is progressively absorbed; and head-on kink-antikink collisions have a critical velocity below which the pair traps into a breather whose oscillations decay. For the spatial Riesz order, $\\alpha<2$ gives monotonic power-law tails and purely attractive kink-antikink forces, while $\\alpha>2$ produces non-monotonic tails with a single zero crossing on each side; the zero crossing converts the tail-mediated force to repulsion at large separation, so a saddle (unstable equilibrium) kink-antikink bound state exists for all $\\alpha>2$, bifurcating from infinite separation as $\\alpha\\to 2^+$. With both derivatives active, the Caputo effect dominates the overall decay, but the Riesz order sets the interaction landscape, yielding long-lived or permanent kink-antikink states that settle at fixed positions.","pith_inferences":["Beyond the paper's explicit claims, the generic dissipative effect of $\\beta<2$ suggests that the Caputo time operator could serve as a minimal phenomenological friction term for any soliton-bearing fractional PDE, without adding a separate damping term to the model.","The zero-crossing tail mechanism for $\\alpha>2$ points to a design principle: by tuning the spatial derivative order, one can switch solitary-wave interactions in power-law-tailed systems between purely attractive and attract-repel landscapes, which could be probed in engineered dispersive media.","The observed settling of kinks at $\\pm L/4$ under periodic boundary conditions is likely a finite-domain signature of the long-range interaction; using larger domains or absorbing boundaries should test whether the repulsion is intrinsic or enhanced by image forces.","If the claimed $\\beta$-independence of the static-kink spectrum holds, the Caputo-induced decay is not a spectral instability but a secular effect of the fractional time evolution, implying that the decay rate should be controlled by $\\beta$ alone rather than by the perturbation eigenfunction."],"forward_implications":["For $\\beta<2$, any kink launched with finite velocity eventually comes to rest, with exponential velocity decay whose characteristic time diverges as $\\beta\\to 2$.","Kink-antikink collisions in the Caputo case acquire a critical velocity depending on $\\beta$ and initial separation: below it the pair forms a damped breather and annihilates; above it they pass and subsequently decelerate.","For $\\alpha>2$, the kink-antikink interaction is attractive at short range and repulsive at long range, with an unstable stationary bound state at the separatrix; for $\\alpha<2$ the interaction is purely attractive.","For $\\alpha>2$, stationary kinks have non-monotonic tails with a single zero crossing on each side, and the band-edge resonance of the integrable case becomes an internal breathing mode.","For $1<\\beta<2$ and $\\alpha\\geq 2$, breather-initiated kink-antikink pairs separate and settle at fixed positions rather than annihilating, with the settling time growing as $\\alpha\\to 2^+$.","For $1<\\beta<2$ and $\\alpha\\geq 2$, breather-initiated kink-antikink pairs separate and settle at fixed positions rather than annihilating, with the settling time growing as $\\alpha\\to 2^+$."],"supporting_citations":[{"why":"Supplies the integrable $\\alpha=\\beta=2$ sine-Gordon baseline, the kink solution, Lorentz boost, and breather framework that the fractional study modifies.","marker":"[15]"},{"why":"Provides the spatially fractional $\\phi^4$ analogue whose kink-tail and spectral phenomenology the Riesz study mirrors and extends.","marker":"[14]"},{"why":"Introduces the dissipation-preserving numerical scheme for the Caputo-Riesz fractional wave equation and earlier observations of breather decay in the time-fractional sine-Gordon model.","marker":"[20]"},{"why":"Reports the dissipative effect of Caputo-time-fractional derivatives on nonlinear wave solutions, forming the basis for the $\\beta<2$ decay claim.","marker":"[23]"},{"why":"Supplies the numerical integrator fde12 used in the paper's Caputo-fractional simulations.","marker":"[22]"},{"why":"Gives the tail-mediated-force picture used to connect the kink-tail zero crossing to the sign change of the kink-antikink interaction.","marker":"[27]"},{"why":"Provides the Mittag-Leffler fractional-identity reference invoked in the Appendix's separation-of-variables stability argument.","marker":"[34]"}],"fun_headline_variants":["Fractional time damps kinks; fractional space flips their force","Caputo time adds drag, Riesz space can reverse pull","Sine-Gordon kinks: fractional order in time damps, in space repels","Two fractional derivatives tune kink damping and kink force sign","Fractional sine-Gordon: time order drags, space order can repel"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The Appendix's proof that kink stability spectra do not depend on $\\beta$ assumes that $\\partial_t^\\beta E_\\beta(\\lambda^2 t)=\\lambda^2 E_\\beta(\\lambda^2 t)$, but the standard Caputo identity for the Mittag-Leffler function applies to $E_\\beta(\\lambda t^\\beta)$ rather than $E_\\beta(\\lambda^2 t)$, so the separation-of-variables step is not a valid derivation as written.","fun_headline_variants_meta":{"raw":{"variants":["Fractional time damps kinks; fractional space flips their force","Caputo time adds drag, Riesz space can reverse pull","Sine-Gordon kinks: fractional order in time damps, in space repels","Two fractional derivatives tune kink damping and kink force sign","Fractional sine-Gordon: time order drags, space order can repel"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000685,"raw_usage":{"total_tokens":3105,"prompt_tokens":944,"completion_tokens":2161,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":2065}},"tokens_in":560,"tokens_out":2161,"duration_ms":15156,"temperature":1.0,"reasoning_tokens":2065,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:01:48.772982+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $\\partial_t^\\beta E_\\beta(\\lambda^2 t)$ numerically for $\\beta=1.5$ and $\\lambda^2<0$ and compare it with $\\lambda^2 E_\\beta(\\lambda^2 t)$; if the two differ, the spectral-independence claim is unsupported. Alternatively, solve the linearized fractional equation for a static kink at $\\beta=1.5$ and compare its eigenvalues with the $\\beta=2$ case: any $\\beta$-dependence of the spectrum would falsify the claim.","supporting_citations":[{"cited_title":"Cuevas-Maraver, P","cited_arxiv_id":null,"evidence_quote":"Supplies the integrable $\\alpha=\\beta=2$ sine-Gordon baseline, the kink solution, Lorentz boost, and breather framework that the fractional study modifies."},{"cited_title":"Fractional Solitons: A Homotopic Continuation from the Biharmonic to the Harmonic $\\phi^4$ Model","cited_arxiv_id":"2410.18426","evidence_quote":"Provides the spatially fractional $\\phi^4$ analogue whose kink-tail and spectral phenomenology the Riesz study mirrors and extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the dissipation-preserving numerical scheme for the Caputo-Riesz fractional wave equation and earlier observations of breather decay in the time-fractional sine-Gordon model."},{"cited_title":"Bountis, J","cited_arxiv_id":null,"evidence_quote":"Reports the dissipative effect of Caputo-time-fractional derivatives on nonlinear wave solutions, forming the basis for the $\\beta<2$ decay claim."},{"cited_title":"Garrappa, Numerical solution of fractional differential equations: A survey and a soft- ware tutorial, Mathematics 6 (2018) 16","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical integrator fde12 used in the paper's Caputo-fractional simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the tail-mediated-force picture used to connect the kink-tail zero crossing to the sign change of the kink-antikink interaction."},{"cited_title":"Mainardi, Fractional relaxation-oscillation and fractional diffusion-wave phenomena, Chaos, Solitons & Fractals 7 (9) (1996) 1461–1477","cited_arxiv_id":null,"evidence_quote":"Provides the Mittag-Leffler fractional-identity reference invoked in the Appendix's separation-of-variables stability argument."}],"review_version":1}