{"id":"53ebd7e2-fd23-4de8-b6b1-dc97e15a6fa4","arxiv_id":"2412.21170","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A spherical-harmonic reduction turns the spectral stability of one- and bi-frequency solitary waves in the 3D Soler model into radial ODE systems, with bi-frequency waves potentially more stable.","lead":"This paper reduces the spectral stability analysis of solitary waves in the 3D Soler model to a set of one-dimensional radial problems, organized by spherical harmonic degree. The same reduction is applied to bi-frequency solitary waves, showing that their stability analysis has the same structure and, under a stated condition, can be better.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bi-frequency spectral reduction is proved only for aligned ξ and η; the generic non-parallel case lacks the completeness lemma needed to rule out eigenvalues outside X_ℓ ⊕ Y_ℓ.","rationale":"The reduction to radial operators is the main contribution, and the aligned bi-frequency computation is a plausible extension of the one-frequency method. The reader's conditional verdict is appropriate. My stress-test focuses on the completeness of the invariant-subspace decomposition, which is the step that converts 'we found invariant subspaces' into 'there are no other unstable eigenvalues.' For one-frequency waves, Lemma 3.2 is asserted by inspection and would benefit from an explicit proof; for bi-frequency waves, the paper itself restricts the completeness proof to ξ,η parallel, while Theorem 4.1 is stated for arbitrary η. This mismatch means the genuinely new non-SU(1,1) bi-frequency regime is not fully covered. I do not see an internal contradiction or a reason to reject; the remedy is a completeness proof or an explicit counterexample. The improved-stability remark is explicitly conditional on a 2D pattern and does not by itself support a stronger conclusion. Hence the verdict should remain conditional pending the non-aligned completeness check.","tokens_in":26448,"tokens_out":26097,"duration_ms":295645,"concrete_test":"Fix ℓ and use normalized data ξ=e1, η=ε(cos θ e1 + sin θ e2) with 0<ε<1 and θ not 0 or π (i.e. η not parallel to ξ). Using the explicit formulas for σ_rΣ_Ω in (2.10) and (2.8), compute the rank of the finite-dimensional linear map that sends the coefficient values (A_{ℓ,m},B_{ℓ,m},P_{ℓ,m},Q_{ℓ,m},R_ℓ,S_ℓ) to the full 4(2ℓ+1)-dimensional space of four-spinor angular components of degree ℓ. If the rank is deficient for some ℓ, exhibit one perturbation in the cokernel, for example ζ(r)h_{ℓ,ℓ}(θ,φ)e2, and verify that it cannot be represented in X_ℓ ⊕ Y_ℓ; this would falsify the claimed reduction for non-parallel bi-frequency waves. If the map is surjective for all ℓ, write the inverse or triangular solve, which supplies the missing completeness proof for Lemma 5.1 in the generic case.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central reduction claim requires that the invariant subspaces cover all of L^2. For one-frequency waves this is Lemma 3.2, whose proof is only 'follows by inspection'; for bi-frequency waves the corresponding completeness statement is Lemma 5.1, and Section 5 explicitly proves it only when ξ and η are parallel. Theorem 4.1, however, is stated with ξ parallel to e1 and arbitrary η, and Section 4 announces a decomposition δ for general bi-frequency solitary waves without a proof. The non-parallel case is the generic regime in which the bi-frequency wave cannot be obtained from a one-frequency wave by the SU(1,1) symmetry, so this is exactly the novel regime in the abstract. If for such η there is an L^2 perturbation outside X_ℓ ⊕ Y_ℓ, the reduced systems (4.26)-(4.28) could miss an unstable eigenvalue. Remark 4.2 does not repair this: its improved-stability comparison is restricted to the aligned case and depends on a spectral pattern observed numerically in (2+1)D, not proved in 3D. The gap is load-bearing because completeness is what excludes all other spectral directions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a radial reduction of the linearized Soler model in (3+1) dimensions around one-frequency and bi-frequency solitary waves. For one-frequency waves, the authors exhibit invariant subspaces X_{l,m} and Y_l labeled by spherical harmonics, derive the reduced radial systems (3.11)-(3.15), and claim that spectral stability reduces to studying the radial operators A_{l,m}. For bi-frequency waves, they introduce analogous invariant subspaces X_l and Y_l, derive the reduced system (4.25) and its aligned special case (4.28), and argue via Remark 4.2 that bi-frequency waves may have stability comparable to or better than one-frequency waves. The paper is largely self-contained, with detailed algebraic derivations and appendices on spherical harmonics and the spin-orbit operator.","tokens_in":26635,"tokens_out":19647,"duration_ms":173416,"significance":"If the reduction is complete, this is a significant technical contribution: it is the first systematic reduction of the spectral stability problem for 3D Soler solitary waves to one-dimensional radial systems, enabling numerical computation of spectra for perturbations of arbitrary angular structure. The derivations are explicit and parameter-free, and the appendices provide useful identities for the spin-orbit operator. The one-frequency reduction and the aligned bi-frequency reduction are plausible and well-documented. However, the completeness of the invariant subspace decomposition is not established for the generic non-parallel bi-frequency case, which is exactly the novel regime not reachable by the SU(1,1) symmetry, so the central claim for bi-frequency waves is not fully supported as stated.","major_comments":[{"comment":"The completeness of the invariant subspaces for bi-frequency waves is proven only when ξ and η are parallel to e1. Theorem 4.1 is stated with ξ parallel to e1 and arbitrary η, and the paragraph preceding it announces a decomposition δ for general ξ, η, but no proof of surjectivity is given for the non-parallel case. The non-parallel, non-orthogonal pair (ξ,η) is precisely the generic case that cannot be obtained from a one-frequency solitary wave by an SU(1,1) transformation, and it is the case in which the claimed improvement of stability would be new. If for such η there is an L^2 perturbation outside X_l ⊕ Y_l, the reduced system (4.26)–(4.28) would miss the corresponding eigenmode. This gap is load-bearing for the central claim that spectral stability reduces to the radial operators, so the bi-frequency reduction is incomplete for the generic case.","section":"Section 5, Lemma 5.1 and Theorem 4.1"},{"comment":"The one-frequency reduction rests on the assertion that the union of the invariant subspaces X_{l,m} and Y_l spans L^2(R^3,C^4). The proof of Lemma 3.2 is given as 'The above lemma follows by inspection.' Given that this completeness is what excludes all other spectral directions, this is not a trivial statement; a detailed proof or a precise reference to a standard expansion in spinor spherical harmonics is needed. Without it, the reduction to the operators A_{l,m} is not fully justified.","section":"Lemma 3.2"},{"comment":"The second line of the system (4.25) reads i∂t P^h = −(A' + κ_l/r^2 B + (g−ω)P)h − 2g' u Re(...). However, the derivation from (4.14) and the analogous one-frequency system (3.11) give the coefficient −(g+ω)P, not −(g−ω)P. This sign error changes the reduced operator and would propagate to (4.28) and Remark 4.2, so it must be corrected.","section":"Eq. (4.25)"}],"minor_comments":[{"comment":"The abstract and introduction state that the technique is applied to bi-frequency solitary waves without qualification, whereas the completeness proof is restricted to the aligned case. Please add a caveat or adjust the claims accordingly.","section":"Abstract and Introduction"},{"comment":"The improved-stability statement is conditional on a spectral pattern observed numerically in (2+1)D, not proved in 3D. The remark should state explicitly that this is a heuristic transfer of a (2+1)D numerical observation, not a theorem in the present setting.","section":"Remark 4.2"},{"comment":"The radial functions A_{l,m}, B_{l,m}, etc. are said to belong to S(R^3,C^4), but they are functions of r only; this should be S(R_+) or similar.","section":"Lemma 3.2 and Lemma 5.1"},{"comment":"The system (4.25) uses the symbol (3−n)/r without defining n in Section 4; since the paper is set in n=3, the term vanishes and can be omitted or n should be specified.","section":"Eq. (4.25)"},{"comment":"In the proof of Lemma 5.1, the coefficients k_{l,m} are introduced without an explicit formula; please provide their definition or a reference.","section":"Lemma 5.1"},{"comment":"There is a parenthesis mismatch in the display: 'σrΣ¯hℓ,m)η' should be 'σrΣ(¯hℓ,mη)' or similar.","section":"Eq. (4.7)"}],"recommendation":"major_revision","confidential_remarks":"The paper has a solid algebraic core, and the one-frequency reduction plus the aligned bi-frequency case are likely publishable after revision. The main issue is the unproven completeness for generic non-parallel η, which is the regime that makes the bi-frequency analysis novel. The authors either need to prove completeness for arbitrary η or clearly restrict the main claims to the aligned case and adjust the title and abstract. The sign error in Eq. (4.25) should also be fixed. If the authors can supply the missing completeness argument, the paper would merit acceptance; otherwise, the claims should be scaled back accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the one-frequency radial reduction is the real meat here, and it looks solid. For the first time, spectral stability of one-frequency solitary waves in the 3D Soler model is reduced to radial operators A_l,m for every spherical harmonic degree, not just radial perturbations. The algebra in Section 3 is detailed and self-contained, and the Y_l subspace that is automatically spectrally neutral is a nice structural observation. This is a useful technical advance for anyone doing numerical stability computations in nonlinear Dirac models.\n\nThe bi-frequency part is more delicate. The aligned case (xi and eta parallel) is worked out cleanly, including the effective angular momentum m' = (1+2nu^2)m, and the claim that these waves can have stability similar to one-frequency waves is plausible. But the completeness of the invariant subspace decomposition is only proved for parallel xi and eta, in Section 5. Theorem 4.1 is stated for arbitrary eta while xi is parallel to e1, and the reduction for the generic non-parallel case is asserted without a completeness argument. That is load-bearing: without completeness, the reduced system could miss an unstable eigenvalue living outside X_l and Y_l. The authors do flag this themselves, calling the parallel case \"the most interesting one,\" so it is not hidden, but it does mean the abstract overstates the proven scope.\n\nTwo smaller things. Lemma 3.2, the one-frequency completeness, is dismissed with \"follows by inspection.\" It is probably true and can be proven with standard spinor spherical harmonic completeness, but a serious proof should be supplied. And Remark 4.2's improved-stability claim is explicitly conditional on a pattern observed numerically in (2+1)D, not proved in 3D; that is stated carefully enough, but readers should not walk away thinking 3D bi-frequency stability is established.\n\nIf I were refereeing, I would ask for (i) a real proof of Lemma 3.2, and (ii) either a completeness proof for the non-parallel bi-frequency case or an explicit narrowing of Theorem 4.1 to the proven aligned case. Neither request looks impossible. The paper deserves a serious referee and likely a conditional accept after revision; it is not ready in its current form.","headline":"The one-frequency radial reduction is a genuine advance; the bi-frequency completeness gap is real but openly acknowledged, so this deserves refereeing with a request for fixes.","tokens_in":27166,"tokens_out":1615,"would_cite":true,"duration_ms":17957,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q41","35B35","35C08"],"pacs":[],"model":"deepseek-v4-flash","headline":"Soler-model solitary wave stability reduces to a family of radial operators.","keywords":["Soler model","nonlinear Dirac equation","spectral stability","solitary waves","radial reduction","spherical harmonics","bi-frequency solitary waves","spin-orbit operator"],"falsifier":"Compute the full linearized spectrum in three dimensions for a single one-frequency solitary wave at a frequency near $\\omega_*\\approx 0.936m$; any eigenvalue with nonzero real part that does not show up in one of the radial systems $A_{\\ell,m}$ would disprove the claimed reduction. Equivalently, for a bi-frequency wave with $\\xi$ and $\\eta$ not parallel, solve the full linearization numerically and compare against the reduced system (4.25): a missing unstable mode would show that the non-parallel case is not covered.","tokens_in":26241,"feed_emoji":"🌀","tokens_out":6396,"duration_ms":59157,"temperature":0.7,"pith_summary":"The paper claims that spectral stability of solitary waves in the three-dimensional Soler model is exactly a radial problem. For a one-frequency solitary wave, the linearized Dirac operator splits the space of perturbations into invariant subspaces labeled by spherical-harmonic degree $\\ell$ and azimuthal order $m$; any eigenvalue with nonzero real part must appear in one of the radial operators $A_{\\ell,m}$ written in (3.15). The same decomposition is applied to bi-frequency solitary waves, yielding the reduced system (4.25)–(4.28); for polarizations parallel to $e_1$ the effective azimuthal order becomes $m'=(1+2\\nu^2)m$, so increasing the frequency mismatch $\\nu$ can lift the reduced system above the stability threshold of the one-frequency wave. The paper thus provides a route to numerical spectral computations in one radial variable for each angular sector.","feed_headline":"Stability of Soler solitons reduces to radial spectra","feed_subtitle":"A spherical-harmonic split turns 3D Dirac stability checks into one-dimensional radial problems.","key_machinery":"The load-bearing object is the angular decomposition of four-component spinors via spherical harmonics and the spin-orbit operator $S = r\\partial_r - r\\alpha_r \\alpha\\cdot\\nabla + \\frac{n-1}{2}$, whose commutation relations (Lemmas B.2–B.5) make the Dirac operator act within each angular sector. The one-frequency linearization $L = D_0 + g\\beta + 2g'\\beta\\varphi\\,\\mathrm{Re}(\\varphi^*\\beta\\,\\cdot) - \\omega$ leaves the subspaces $X_{\\ell,m}$ and $Y_\\ell$ invariant; projecting onto a sector yields the radial operator $L_0(\\omega)$ in (3.13) and the full reduced operator $A_{\\ell,m}$ in (3.15). The companion Lemma 2.1 computes the coefficients $C_{\\ell,m,k}$ that couple azimuthal orders, and in the parallel-polarization case these coefficients collapse to $-(|\\xi|^2+|\\eta|^2)m\\,\\delta_{m,k}$, producing the effective quantum number $m'=(1+2\\nu^2)m$. The $Y_\\ell$ sector is handled separately by a symmetric radial operator $A_{RS}$ whose spectrum is necessarily imaginary.","core_discovery":"On the paper's own terms, the central discovery is that the stability question for the Soler model in (3+1)D reduces to the spectra of the operators $A_{\\ell,m}$ defined in (3.15) for one-frequency solitary waves, and to the analogue defined by (4.25)–(4.28) for bi-frequency waves. The reduction is achieved by exhibiting invariant subspaces $X_{\\ell,m}$ and $Y_\\ell$ whose elements have definite angular structure: upper and lower spinor components are built from spherical harmonics $h_{\\ell,m}$ applied to fixed polarization vectors, with radial coefficients $A,B,P,Q$. The $Y_\\ell$ subspace is shown to contribute only purely imaginary spectrum, so unstable eigenvalues must lie in the $X$-sector. For bi-frequency waves with $\\xi$ and $\\eta$ parallel to $e_1$, Lemma 2.2 makes the angular coupling diagonal, and the reduced operator coincides with the one-frequency operator with azimuthal order $m'=(1+2\\nu^2)m$; this is the basis for the claim that bi-frequency waves can have stability properties the same as or better than one-frequency waves.","pith_inferences":["One testable extension: run the same radial reduction for a Dirac–Klein–Gordon model with Yukawa coupling, which shares the $SU(1,1)$ symmetry; if the same angular algebra holds, bi-frequency stability comparisons should carry over.","For non-parallel $\\xi,\\eta$, the coupling matrix $C_{m,k}$ mixes azimuthal orders instead of merely rescaling $m$; a finite-dimensional angular calculation would show whether the large-$\\nu$ stability gain survives away from the parallel case.","The reduction suggests a practical numerical protocol: discretize only the radial variable in each sector and compare against a full 3D spectral computation; agreement would close the completeness gap left by the 'by inspection' lemma.","The observed threshold $\\omega_*\\approx 0.936m$ from radial numerics becomes a prediction for every angular sector: if the reduced operators are computed for all $\\ell,m$, the threshold should be the same."],"forward_implications":["Spectral-stability checks for one-frequency Soler solitary waves reduce to solving the radial systems $A_{\\ell,m}$ for $\\ell\\in\\mathbb N_0$ and $0\\le m\\le\\ell$; the $Y_\\ell$ directions are automatically neutrally stable.","The eigenvalue $\\lambda=-2\\omega$ of geometric multiplicity two is forced by the $SU(1,1)$ symmetry and always sits in the spectrum; it does not by itself mean instability.","For bi-frequency waves with parallel polarizations, the effective azimuthal order is $(1+2\\nu^2)m$, so a wave with large $\\nu$ can have a stable reduced spectrum even if the corresponding one-frequency wave at the same $\\omega$ is unstable.","In the parallel-polarization case every Schwartz-class perturbation is decomposed into the invariant subspaces, so no unstable eigenmode can hide outside the reduced systems."],"supporting_citations":[{"why":"Establishes the SU(1,1) symmetry and the existence of bi-frequency solitary waves; the present work revisits and corrects its stability analysis.","marker":"[BC18]"},{"why":"Provides the numerical spectrum of the radial linearization in the one-frequency case and the observed stability threshold near omega*=0.936m that motivates the reduction.","marker":"[CMKS+16]"},{"why":"Supplies the existence theory for the solitary-wave profiles v(r,omega), u(r,omega) used in the ansatz.","marker":"[ES95]"},{"why":"Introduces the Soler model cubic self-interaction f(tau)=tau that the paper's stability analysis is built on.","marker":"[Sol70]"},{"why":"Derives the radial profile system for v and u whose solutions are the one-frequency solitary waves under study.","marker":"[BC17]"}],"fun_headline_variants":["Soler soliton stability reduces to radial spectra","Spherical harmonics cut Soler stability to 1D","Bi-frequency Soler waves match one-frequency stability","Radial reduction tames Soler spectral stability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole reduction rests on the completeness of the invariant-subspace decomposition, which the paper asserts by inspection for one-frequency waves and proves only for parallel polarizations in the bi-frequency case, leaving non-parallel bi-frequency waves an unproven assertion.","fun_headline_variants_meta":{"raw":{"variants":["Soler soliton stability reduces to radial spectra","Spherical harmonics cut Soler stability to 1D","Bi-frequency Soler waves match one-frequency stability","Radial reduction tames Soler spectral stability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00016,"raw_usage":{"total_tokens":1187,"prompt_tokens":854,"completion_tokens":333,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":470,"completion_tokens_details":{"reasoning_tokens":272}},"tokens_in":470,"tokens_out":333,"duration_ms":3666,"temperature":1.0,"reasoning_tokens":272,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:00:22.698351+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full linearized spectrum in three dimensions for a single one-frequency solitary wave at a frequency near $\\omega_*\\approx 0.936m$; any eigenvalue with nonzero real part that does not show up in one of the radial systems $A_{\\ell,m}$ would disprove the claimed reduction. Equivalently, for a bi-frequency wave with $\\xi$ and $\\eta$ not parallel, solve the full linearization numerically and compare against the reduced system (4.25): a missing unstable mode would show that the non-parallel case is not covered.","supporting_citations":[],"review_version":1}