{"id":"e1826450-2b8d-4a69-a3b7-f07e7d0aaee8","arxiv_id":"2411.17336","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Analytic stability boundaries for pinned, split-pinned, sync-dot, and phase-locked states are derived for forced swarmalators in 2D and 3D periodic domains, extending previous 1D results.","lead":"This paper extends studies of periodically forced swarmalators from one-dimensional motion to two and three dimensions, deriving analytic stability boundaries for four static collective states. It also argues that simplified periodic-boundary models reproduce the qualitative dynamics of the original, harder-to-analyze model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The sync-dots stability boundary rests on an unproven reduction to Δy=0; if the stability of the full fixed-point manifold differs, Eqs. (19) and (32) are not exact.","rationale":"The reader's weakest-assumption analysis identifies exactly the step that determines whether the central analytic result is exact: the sync-dots stability boundary is derived from a special case of the fixed-point manifold and then asserted to hold for the full state on the basis of numerics. My independent reading of Section IV.A.3 confirms that the text explicitly says 'Numerics suggests' after setting Δy = 0, and the same procedure is repeated in 3D. This is not a minor gap: Eq. (19) and its 3D analogue Eq. (32) are among the paper's main analytic boundaries, and the p-dimensional extrapolation in Table II inherits the same assumption. The other static-state analyses (pinned, split pinned, phase locked) follow standard linear-stability arguments and do not present a comparable unproven reduction. The proposed test directly probes whether the Δy = 0 slice is lossless by comparing eigenvalues on the full fixed-point manifold. Until that check is performed, the correctness risk is real but not established, so the appropriate verdict remains CONDITIONAL, matching the reader's assessment; no verdict change is needed.","tokens_in":12938,"tokens_out":12639,"duration_ms":123401,"concrete_test":"For N = 2n with two equal clusters, linearize Eqs. (3) at a generic sync-dots fixed point (x1,y1,θ1) = (0,0,−π/4), (x2,y2,θ2) = (Δx,Δy,π/4), with cos Δx + cos Δy = −√2 F/K. For K = −3, choose F = 3.8, which lies inside the Eq. (19) interval, and compute the 3N×3N Jacobian eigenvalues at the Δy = 0 slice and at a generic point such as Δx = Δy with 2 cos Δx = −√2 F/K. Compare the maximal real part of the eigenvalues. If the sign changes across the manifold at fixed (K,F), Eq. (19) is not the exact boundary; if it does not, repeat the comparison at several (K,F) points on both sides of the claimed boundary to confirm the slice is representative.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section IV.A.3, after deriving the sync-dots fixed-point condition cos Δx + cos Δy = −√2 F/K (Eq. 17), the authors set one of the two spatial separations to zero, leaving a one-parameter slice of the fixed-point manifold, and then compute eigenvalues (18) for that slice. The text then states: 'Numerics suggests that the acquired stability condition due to the consideration of these fixed points provides the exact stability boundary for the sync dots state.' This is the load-bearing step for the advertised boundary Eq. (19), and the same unproven reduction is used in 3D in Eq. (32). The Jacobian at a generic sync-dots fixed point depends on Δx and Δy separately, not only on their cosine sum: the linearized position equations acquire coupling coefficients proportional to sin Δx and sin Δy, and the phase equation depends on cos Δx + cos Δy. No symmetry argument or coordinate transformation is supplied that would make Δy = 0, or Δy = Δz = 0 in 3D, a lossless reduction. Because the claimed stability region is central to the paper's advertised analytic contribution, and because the boundary is asserted rather than derived for the full fixed-point family, the central claim is not fully established. A numerical phase diagram can corroborate the boundary for selected parameter values but cannot establish that the slice is representative.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies forced swarmalator models in two and three spatial dimensions with periodic boundary conditions, extending a previously studied one-dimensional model. For the 2D and 3D models the authors derive analytic stability boundaries for four static collective states (pinned, split-pinned, sync dots, and phase locked) using linear stability analysis, and they support these boundaries with numerical phase diagrams. They also report that the simplified periodic-boundary models reproduce the collective states observed in a more realistic power-law swarmalator model, and they propose a p-dimensional extrapolation of the boundaries. The paper claims these are the first analytical results for forced swarmalators in two and three dimensions.","tokens_in":13205,"tokens_out":6258,"duration_ms":57771,"significance":"If the derived boundaries are correct, the paper makes a useful contribution by moving the analytic theory of forced swarmalators beyond one dimension and by identifying simple formulas that scale with the spatial dimension. The pinned and split-pinned analyses are clean, self-contained linear-stability derivations with no fitted parameters, and the numerical phase diagrams give concrete, testable predictions. The sync-dots and phase-locked boundaries, however, rest on an unproven reduction to a special submanifold of fixed points, so the central claim of exact analytical boundaries is not fully established as written. The paper is clearly organized and the numerical methodology is adequately documented, but the load-bearing analytical step for the sync-dots state needs either a rigorous justification or an explicit downgrade to a conjecture.","major_comments":[{"comment":"The stability boundary for the sync dots state is derived by setting one of the two spatial separations to zero, Delta y = 0, and the text asserts that numerics suggests this slice provides the exact stability boundary for the full sync-dots state. The Jacobian at a generic sync-dots fixed point, however, depends on sin Delta x and sin Delta y separately: for example, the phase equation linearizes to terms proportional to (K/2) sin Delta x and (K/2) sin Delta y, and the position-phase coupling contains analogous factors, so the eigenvalues are not determined solely by cos Delta x + cos Delta y. No symmetry argument or coordinate transformation is supplied that makes the Delta y = 0 slice lossless. Because Eq. (19) is the advertised stability boundary for sync dots and also provides one endpoint of the phase-locked boundary Eq. (24), this unproven reduction is load-bearing. Please either prove that the stability condition is independent of the slice (e.g., by an explicit eigenvalue computation for the full fixed-point family) or explicitly state Eq. (19) as a numerically supported conjecture rather than an exact analytical result.","section":"Section IV.A.3, Eqs. (17)-(19)"},{"comment":"The same unproven reduction is used in the three-dimensional analysis, where two of the three spatial separations are set to zero before computing the eigenvalues of the sync-dots Jacobian. The resulting boundary Eq. (32) is presented as exact, and the phase-locked boundary Eq. (35) inherits this issue. In 3D the Jacobian at a generic sync-dots fixed point will contain coupling terms proportional to sin Delta x, sin Delta y, and sin Delta z, so the reduction is not obviously lossless. The claim that Eq. (32) is the exact stability boundary therefore needs the same kind of justification as in the 2D case; otherwise it should be labeled as a conjecture supported by the numerics shown in Fig. 8.","section":"Section IV.B.3, Eqs. (30)-(32)"},{"comment":"The derivation of the sync-dots eigenvalues assumes the population splits evenly between the two clusters, with the statement 'without loss of generality' and a further assertion that numerics indicate stability is independent of the population split. This is another unproven assumption in the sync-dots analysis. If the exact stability boundary really is independent of the cluster population split, this should be demonstrated analytically within the reduced model or, failing that, explicitly listed as part of the numerical conjecture.","section":"Section IV.A.3, paragraph after Eq. (17)"}],"minor_comments":[{"comment":"In the Jacobian block M1 for the pinned state, the off-diagonal entry should be 2K/N rather than K/N, because at the pinned fixed point cos(x_j - x_i) + cos(y_j - y_i) = 2. As printed, the matrix is inconsistent with the stated eigenvalues -F and -F - 2K. The same factor-of-two issue appears in the off-diagonal entries of Eq. (14) for the split-pinned state. The eigenvalues listed in Eqs. (11) and (15) are consistent with the 2K/N off-diagonal, so I assume this is a typographical error, but it should be corrected.","section":"Eq. (10) and Eq. (14)"},{"comment":"The p-dimensional column of Table II is presented without derivation, and Section V says the framework 'allows us to predict analytical phase boundaries' in p dimensions. Since the formulas for p = 2 and p = 3 are themselves only conditionally established (see the sync-dots comments), the p-dimensional extrapolation should be explicitly labeled as conjectured, or a derivation should be supplied.","section":"Table II and Section V"},{"comment":"The claim that these are 'the first analytical results about swarmalator with forcing in two and three spatial dimensions' may overstate the scope, since the analysis applies to the periodic-boundary simplification rather than the original power-law model. I suggest wording such as 'first analytical stability results for the periodic-boundary forced swarmalator model in two and three dimensions'.","section":"Abstract and Introduction"},{"comment":"There are several typographical errors: 'Notie' in the Table I caption, 'Varaition' in the Fig. 7 caption, 'Mspilt pinned' in Section IV.A.2, 'The details analysis' in Section II, and 'C()' where C_x, C_y are meant in Section IV.A.3. These should be corrected in a revision.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of nlin.CD and the pinned/split-pinned analyses are solid, but the sync-dots stability boundary, which is central to the paper's advertised analytic contribution, rests on an unproven reduction to a lower-dimensional slice of the fixed-point manifold. This is fixable within the manuscript's scope: either the authors can prove the reduction (for example, by showing the relevant eigenvalues are independent of the slice) or they can reframe the sync-dots and phase-locked boundaries as numerically supported conjectures. The factor-of-two inconsistency in the printed Jacobians also suggests the manuscript needs a careful proofreading pass before it can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real advance here is a tractable generalization of the 1D forced swarmalator model to 2D and 3D periodic domains, with explicit stability curves for several static states. The pinned and split-pinned analyses are clean linear stability calculations, and the eigenvalues are internally consistent modulo a likely typo in the printed Jacobians (the off-diagonal entries in Eqs. 10 and 14 look like they should carry a factor 2). The 3D results are a mechanical extension of the 2D case, and the table of boundaries across dimensions is a nice summary, even if the p-dimensional extrapolation is a conjecture rather than a derivation.\n\nThe soft spot is the sync-dots boundary. The fixed-point condition leaves a one-parameter family, and the paper then sets one of the spatial separations to zero without justification. The Jacobian at a generic sync-dots fixed point depends on both sin Δx and sin Δy, so there is no obvious symmetry that makes the slice lossless. The text says 'numerics suggests' this gives the exact boundary, but that is an assertion, not a proof. This matters because the phase-locked lower boundary is derived from the sync-dots bifurcation, so two of the four advertised boundaries rest on this ansatz. I would not call the sync-dots boundary an analytic result as written; it is a numerically supported conjecture.\n\nAlso, the paper ships no code or data for the numerical phase diagrams, which makes independent verification slower. That is a minor issue for a theoretical paper, but it matters here because the numerical phase diagram is the only support for the sync-dots boundary.\n\nOn balance, the paper is worth engaging with. The pinned and split-pinned results are solid, the model is a reasonable simplification, and the authors are transparent about what they could not prove. But the sync-dots claim needs to be rephrased, or a proof of the slice reduction supplied, before it can be called analytic. I would send it to peer review with a request for revision focused on that point.","headline":"Useful extension of the 1D forced swarmalator analysis to 2D and 3D, but the advertised sync-dots stability boundary rests on an unproven slice ansatz and should be presented as numerical instead of analytic.","tokens_in":13746,"tokens_out":5135,"would_cite":true,"duration_ms":55105,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34C15","34D20"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives the first analytical stability boundaries for forced swarmalators moving in two and three dimensions with periodic boundary conditions, giving exact curves for pinned, split-pinned, sync-dots, and phase-locked states.","keywords":["swarmalators","external forcing","phase synchronization","collective states","stability boundaries","periodic boundary conditions","Kuramoto model","bifurcation analysis"],"falsifier":"Simulate Eq. (3) starting from two sync-dot clusters with both $\\Delta x$ and $\\Delta y$ nonzero and with unequal cluster sizes, slowly ramp $F$ at fixed $K<0$, and record the forcing values at which the two dots lose stability or the phase separation leaves $\\theta^* = \\pi/4$; if these values differ from Eq. (19) in 2D or Eq. (32) in 3D, the claimed exact boundary is falsified.","tokens_in":12701,"feed_emoji":"🧲","tokens_out":5445,"duration_ms":48782,"temperature":0.7,"pith_summary":"The paper asks what sinusoidal forcing does to swarmalators—oscillators that move and synchronize simultaneously—when the motion is not confined to a line but takes place in a plane or in three-dimensional space. It proposes tractable versions of the 2D and 3D forced swarmalator models with periodic boundary conditions and derives analytic conditions for the stability of four static collective states: pinned, split pinned, sync dots, and phase locked. If the derivation is right, it supplies the first analytical stability boundaries for forced swarmalators in two and three spatial dimensions, replacing purely numerical phase diagrams with exact curves in the $(K,F)$ parameter plane. The practical stake is control: systems such as magnetic colloids can be steered by external fields, and knowing exactly when a global pinned or locked state is stable tells you when that steering works.","feed_headline":"Forced swarmalators in 2D and 3D get exact stability lines","feed_subtitle":"Sinusoidal forcing pins, splits, and locks swarmalators; the paper maps each transition in the (K, F) plane.","key_machinery":"The technical workhorse is linear stability analysis at fixed points. For the pinned and split-pinned configurations, the Jacobian has a block structure that separates spatial and phase directions; its eigenvalues are computed exactly, and stability boundaries appear as zero-eigenvalue bifurcations, $F_c = -2K$ in 2D and $F_c = -3K$ in 3D. For sync dots and phase-locked states, the analysis reduces the swarmalator equations to fixed-point conditions such as $\\cos\\Delta x + \\cos\\Delta y = -\\sqrt{2}F/K$, and the sync-dots stability boundary is obtained by setting one of the spatial separations ($\\Delta y$ in 2D, two of the three in 3D) to zero, with the resulting condition asserted on numerical grounds to be the exact boundary for the full state.","core_discovery":"On its own terms, the paper's central claim is that a simplified forced swarmalator model on a periodic 2D torus, and its 3D analogue, reproduces the collective states of the full higher-dimensional model while remaining analytically tractable, and that the static states have computable stability boundaries. In 2D these are $F > -2K$ for the pinned state, $F < 2K$ for $K>0$ for split pinned, $-(2K+1)/\\sqrt{2} < F < -2K/\\sqrt{2}$ for sync dots, and $-2K/\\sqrt{2} < F < -2K$ for phase locked, with the sync-dots and phase-locked bands requiring $K<0$. In 3D the same boundaries hold with $K$ replaced by $3K$ where it multiplies the spatial dimension, and the sync-dots band shifted accordingly. The paper also arranges the 1D, 2D, and 3D results in Table II, exposing a $p$-dimensional pattern $F > -pK$ for pinning, and reports numerical evidence that the periodic-boundary model reproduces the behavior of the power-law-kernel forced swarmalator model.","pith_inferences":["The $p$-dimensional scaling visible in Table II is likely exact for all $p$, because the same block-Jacobian argument that gives $F=-2K$ in 2D and $F=-3K$ in 3D should give $F=-pK$ in general, so one could extend the analytic phase diagram to higher-dimensional tori by symmetry alone.","If the sync-dots boundary is truly independent of the spatial separation pattern, it may also control the stability of the coherent sync-dot clusters observed inside the unanalysed chimera state, which would give a first analytical handle on that state.","A direct numerical test of the special-case assumption is possible: integrate Eq. (3) with two sync-dot clusters separated in both $x$ and $y$, measure where the two dots lose stability, and compare with Eq. (19); a mismatch would show the claimed exactness is only approximate.","The reported equivalence with the power-law-kernel model is presently numerical; a stronger version would use the periodic-boundary phase diagram to predict, for each $(K,F)$, which of the six states the full model exhibits and how hard-shell repulsion shifts the state boundaries."],"forward_implications":["The analytical curves in the $(K,F)$ plane, given by Eqs. (12), (16), (19), and (24) in 2D and their 3D analogues, give future numerical studies of forced swarmalators a precise set of bifurcation points to test.","Table II's pattern suggests a general rule in $p$ spatial dimensions: pinning requires $F > -pK$, split pinning requires $F < pK$, and the sync-dots and phase-locked bands scale with $p$ as well.","Because the periodic-boundary model omits hard-shell repulsion yet still produces all the collective states found in the full 2D forced model, it can serve as a minimal reliable model for forced swarmalator dynamics.","The reported numerical agreement between the periodic-boundary model and the power-law-kernel model implies that analytically derived boundaries may transfer approximately to the more realistic model, giving predictions where only simulations existed before.","The stability analysis also fixes where each static state loses stability as forcing strength $F$ is ramped upward, which directly describes how an external field can be tuned to pin or unlock a swarmalator population."],"supporting_citations":[{"why":"It supplies the 1D forced-swarmalator analysis whose eigenvalue method and state classification are extended here to 2D and 3D.","marker":"[36]"},{"why":"It defines the unforced 2D swarmalator model with periodic boundary conditions and its exact phase diagram, the starting point for adding forcing.","marker":"[18]"},{"why":"It introduces the swarmalator model and the rainbow order parameters used throughout to characterize space-phase correlation.","marker":"[17]"},{"why":"It provides the numerical study of the full 2D forced swarmalator model that the simplified periodic-boundary models are intended to reproduce.","marker":"[21]"},{"why":"It motivates the 1D model as capturing the rotational piece of 2D swarmalator motion, justifying the periodic-boundary simplification.","marker":"[37]"}],"fun_headline_variants":["Swarmalators in 2D and 3D yield to exact stability maps","Exact stability lines for forced swarmalators in 2D and 3D","Higher-D forced swarmalators: exact transition boundaries","Map of swarmalator phases in 2D and 3D under forcing","Forced swarmalators: analytic stability in 2D and 3D"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the sync-dots stability boundary obtained by setting one spatial separation to zero, and then only checked numerically, is exactly the boundary of the full sync-dots state with both separations free.","fun_headline_variants_meta":{"raw":{"variants":["Swarmalators in 2D and 3D yield to exact stability maps","Exact stability lines for forced swarmalators in 2D and 3D","Higher-D forced swarmalators: exact transition boundaries","Map of swarmalator phases in 2D and 3D under forcing","Forced swarmalators: analytic stability in 2D and 3D"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001102,"raw_usage":{"total_tokens":4576,"prompt_tokens":906,"completion_tokens":3670,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":3568}},"tokens_in":522,"tokens_out":3670,"duration_ms":21014,"temperature":1.0,"reasoning_tokens":3568,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:13:35.214157+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate Eq. (3) starting from two sync-dot clusters with both $\\Delta x$ and $\\Delta y$ nonzero and with unequal cluster sizes, slowly ramp $F$ at fixed $K<0$, and record the forcing values at which the two dots lose stability or the phase separation leaves $\\theta^* = \\pi/4$; if these values differ from Eq. (19) in 2D or Eq. (32) in 3D, the claimed exact boundary is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the 1D forced-swarmalator analysis whose eigenvalue method and state classification are extended here to 2D and 3D."},{"cited_title":"Riedl, I","cited_arxiv_id":null,"evidence_quote":"It defines the unforced 2D swarmalator model with periodic boundary conditions and its exact phase diagram, the starting point for adding forcing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It introduces the swarmalator model and the rainbow order parameters used throughout to characterize space-phase correlation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It motivates the 1D model as capturing the rotational piece of 2D swarmalator motion, justifying the periodic-boundary simplification."}],"review_version":1}