{"id":"95af8266-86ee-4743-a548-77389c0d2598","arxiv_id":"2608.07939","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Planar spin-orbit coupling stabilizes three-dimensional semi-vortex and mixed-mode solitons in atomic-molecular condensates with quadratic nonlinearity, with vortex components carrying over half the norm.","lead":"This paper reports numerical evidence that stable three-dimensional solitons, including vortex-carrying states, can exist in a spin-orbit-coupled atomic-molecular Bose-Einstein condensate with quadratic interactions. The result matters because stable 3D vortex solitons in free space have been a long-standing goal in nonlinear optics and ultracold matter.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability of the 3D vortex solitons is asserted from finite-time real-time runs with no stated numerical parameters, no full-3D perturbation protocol, and no Bogoliubov-de Gennes spectrum; in quadratic media the azimuthal splitting instability is exactly the mechanism that must be ruled out.","rationale":"The paper proposes a plausible mechanism, and the equations and qualitative cascading argument are coherent. The strongest claim, however, is a stability claim, and the evidence presented is finite-time real-time evolution without numerical details or a linear stability analysis. The reader's weakest assumption identifies exactly this gap. I sharpen it: the ansatz itself builds in the azimuthal phase winding, while the known failure mode for vortex solitons in quadratic media is azimuthal fragmentation. If the numerical work intentionally or accidentally preserved axisymmetry, the reported stability would be an artifact of the symmetry of the initial condition rather than a property of the full 3D system. The Vakhitov-Kolokolov criterion is not a substitute for probing symmetry-breaking modes. A Bogoliubov-de Gennes calculation, or at minimum a documented full-3D run with explicit symmetry-breaking perturbation and stated parameters, would settle the issue. This does not require rejection; it requires conditional acceptance pending that verification. The reader's CONDITIONAL verdict is therefore appropriate and unchanged.","tokens_in":12101,"tokens_out":9730,"duration_ms":118900,"concrete_test":"Run a full three-dimensional Bogoliubov-de Gennes stability analysis for the stationary SV soliton at (N,alpha)=(200,0) and at a point near the claimed boundary, e.g. (200,0.9), linearizing Eqs. (1)-(3) about the converged imaginary-time solution and diagonalizing over all azimuthal indices m=0, +/- 1, +/- 2. If any eigenfrequency has positive imaginary part above numerical tolerance, the stability claim fails; if the spectrum is purely real, the finite-time runs are corroborated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that planar SOC alone stabilizes 3D vortex solitons in the quadratic system—depends on the stability classifications in Fig. 3, which the text supports only by 'long-term real-time propagation' in the captions of Figs. 1 and 2. The paper does not report the integration time, grid resolution, boundary conditions, or perturbation amplitude, nor does it state whether the evolution was performed in full 3D or with imposed axisymmetry. This matters specifically because the SV ansatz (10) enforces the azimuthal structure phi2, psi ~ e^{i theta}: if the code or the initial noise preserved rotational symmetry, the azimuthal splitting instability known to plague vortex solitons in quadratic media [12,13,24,31,32] could not appear. The Vakhitov-Kolokolov criterion d mu/dN < 0 cited in the discussion of Fig. 4 is necessary but not sufficient for vortex stability; it does not probe m != 0 perturbations. Without a Bogoliubov-de Gennes spectrum or a documented symmetry-breaking real-time test, the stability regions in Fig. 3 are not established in exactly the regime where the claim is nontrivial.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a mean-field model of an atomic-molecular Bose-Einstein condensate with planar Rashba spin-orbit coupling applied to the atomic components and quadratic three-wave coupling to a molecular component. Using imaginary-time propagation to construct stationary states, it reports semi-vortex (SV) and mixed-mode (MM) solitons in free 3D space, studies their dependence on the atomic-molecular mismatch alpha, explains the existence cutoff via a cascading approximation, and claims that the SV solitons carry a dominant vortex fraction of the norm. The central physical claim is that planar SOC alone stabilizes 3D vortex-carrying solitons in a quadratic medium, a regime previously thought to be dominated by azimuthal splitting instability. The existence and qualitative structure of the soliton families are plausible, but the stability claim is supported only by finite-time real-time propagation with no numerical parameters, no linear stability analysis, and no documented symmetry-breaking perturbation protocol.","tokens_in":12313,"tokens_out":14707,"duration_ms":173022,"significance":"If the stability claim is correct, the result is significant: it would identify a mechanism for stable vortex solitons in free space with quadratic nonlinearity, going beyond the usual reliance on competing cubic nonlinearities or external potentials. The model and cascading reduction are standard, and the existence families in Figs. 3-5 would be useful additions to the soliton literature. However, the manuscript currently does not supply the evidence needed to distinguish a genuinely stable 3D vortex soliton from a metastable state whose azimuthal instability has not yet developed over the numerical horizon. The central claim therefore needs to be either proven or carefully downgraded before publication.","major_comments":[{"comment":"Stability is the load-bearing claim, but the evidence is not quantitative. The captions of Figs. 1 and 2 state that the solitons are \"confirm[ed]\" to be stable after long-term real-time propagation, and Fig. 3 draws stable and unstable regions from that criterion, yet the paper reports no integration time, grid resolution, boundary conditions, perturbation amplitude, or whether the evolution was carried out in full 3D without imposed rotational symmetry. The azimuthal splitting instability of vortex solitons in quadratic media (Refs. [12,13,24,31,32]) is precisely a non-axisymmetric process; if the solver or the initial state preserves continuous rotational symmetry, that instability cannot develop. Please add a Bogoliubov-de Gennes linear stability analysis with the azimuthal-mode spectrum, or at least full-3D real-time runs initialized with explicit non-axisymmetric perturbations, for representative points in each region of Fig. 3, and state the quantitative stability criterion and the numerical protocol.","section":"Figs. 1-3 and related captions"},{"comment":"The Vakhitov-Kolokolov condition dμ/dN<0 mentioned in the discussion of Fig. 4 is necessary but not sufficient for stability of vortex solitons, because it only constrains rotationally symmetric perturbations and cannot detect the m≠0 azimuthal modes that are dangerous in quadratic media. Since the text itself calls it an indicator of \"plausible stability,\" it cannot justify the stable/unstable boundary in Fig. 3. The stability classification should be based on the linear spectrum or on explicitly documented symmetry-breaking dynamics.","section":"Fig. 4 and discussion"}],"minor_comments":[{"comment":"\"This is scheme for realizing stable vortex solitons\" should read \"This provides a scheme for realizing stable vortex solitons.\"","section":"Abstract"},{"comment":"The angular-momentum integral is written as dxdy, but the fields depend on three coordinates; it should be dxdydz, or the z-integration convention should be explained.","section":"Eq. (8)"},{"comment":"The vortex structure of the molecular component ψ in the SV ansatz is imposed by the phase-matching condition ψ(0)=φ1φ2 and the e^{iθ} factor in ansatz (10); the numerical observation in Fig. 1 is therefore a consistency check of the stationary solution, and the conclusion should not imply that this vortex structure is an emergent prediction.","section":"Eqs. (9)-(10) and Fig. 1"},{"comment":"There are several typographical errors: \"eliminaring\" near Eq. (13), \"repsectively\" in the Fig. 1 caption, \"solitons do not exists\" in the Fig. 4 caption, and \"region is includes\" in the Fig. 3 caption.","section":"Throughout"},{"comment":"Reference [5] is cited as \"Mod. Rev. Phys.\"; the journal is Reviews of Modern Physics.","section":"Reference [5]"},{"comment":"Please state the branch-continuation procedure used to produce the families in Fig. 4, including how N and α are swept and how the imaginary-time convergence criterion is defined.","section":"Fig. 4"},{"comment":"Equation (14) defines Fv using N2+N3, and N3 carries the factor 2 in Eq. (5); please state explicitly that this factor is included and discuss whether the conclusion that the vortex share exceeds 50% depends on that normalization convention.","section":"Eq. (14) and Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The deficit in this manuscript is concentrated in the evidence for the central stability claim rather than in the model formulation. The authors are clearly capable of supplying the missing numerical details and a Bogoliubov-de Gennes analysis; if those are provided and the stability regions survive, the paper would be a solid contribution. If the linear stability analysis reveals azimuthal instabilities, the conclusions will need to be substantially softened. The manuscript should not be accepted without this additional evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper is proposing a genuinely new combination—planar spin-orbit coupling plus quadratic atomic-molecular nonlinearity in 3D—and the idea is plausible. The cascading approximation explains the existence cutoff cleanly, the 39K calibration makes it experimentally addressable, and the claim that >50% of the norm sits in vortex components is a concrete, falsifiable prediction. That part is worth taking seriously.\n\nThe soft spot is the stability claim. The reader's concern is correct: the paper infers stability from 'long-term real-time propagation' without reporting the integration time, grid resolution, boundary conditions, perturbation amplitude, or whether the runs were full 3D or restricted by symmetry. In quadratic media the azimuthal splitting instability is the known killer of vortex solitons, and if the numerics preserved rotational symmetry, that mechanism cannot appear. The Vakhitov–Kolokolov criterion cited for Fig. 4 is necessary but not sufficient; it doesn't see m≠0 perturbations. So the stability regions in Fig. 3, which are the central result, are not established as presented. The vortex structure of the molecular component is also partly built into the ansatz via ψ=φ1φ2, so that observation is less of an independent finding.\n\nThat said, the math is not wrong as far as it goes. The GPE system, conserved quantities, and the SV/MM ansätze are standard, and the cascading-limit argument for α_c is coherent. The missing pieces are reproducibility and a proper stability check. I don't see a load-bearing flaw in the idea itself, just in the evidence for it.\n\nWho is this for: specialists in multidimensional solitons and SOC BECs. A serious referee should be engaged, but the paper needs major revision: a Bogoliubov–de Gennes spectrum or at least documented full-3D real-time runs with randomized symmetry-breaking perturbations, plus all numerical parameters. The typo in the abstract ('This is scheme') is easy to fix.\n\nRecommendation: send to peer review, but expect a heavy revision. I would not cite it in its current form.","headline":"Plausible new mechanism for stable 3D vortex solitons, but the stability claim is not yet backed by the numerics as presented.","tokens_in":12915,"tokens_out":2048,"would_cite":false,"duration_ms":22490,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A flat spin-orbit coupling stabilizes 3D vortex solitons","keywords":["atomic-molecular condensates","parametric interconversion","spin-orbit coupling","semi-vortices","mixed-mode solitons","cascading approximation","soliton stability","quadratic nonlinearity"],"falsifier":"Take the stationary SV and MM solutions shown in Figs. 1 and 2, add a random phase perturbation with amplitude $\\epsilon$ between $10^{-4}$ and $10^{-2}$ of the local density, and evolve Eqs. (1)-(3) for at least an order of magnitude longer than the reported run; if the vortex ring fragments into several filaments or the angular momentum $M_z$ changes by more than numerical error, the stability claim is disproved. A cheaper quantitative check is the Bogoliubov spectrum of the linearized system around these stationary states: a positive imaginary part of any eigenfrequency would falsify asymptotic stability.","tokens_in":2040,"feed_emoji":"🌀","tokens_out":2408,"duration_ms":105333,"temperature":0.7,"pith_summary":"This paper claims that a planar (two-dimensional) spin-orbit coupling applied only to the atomic components of an atomic-molecular Bose-Einstein condensate is enough to stabilize fully three-dimensional solitons in free space, despite the quadratic (three-wave) nonlinearity that would otherwise make vortex states split apart. The solitons come in two types: semi-vortices, in which one atomic component is vortex-free and the other carries unit circulation, with the molecular field inheriting the vortex; and mixed modes, superpositions of zero- and unit-vorticity components with zero total angular momentum. Both families exist only up to a critical detuning $\\alpha_c$, and in the semi-vortex case the vortex-carrying components hold more than half of the total norm. If correct, this addresses an open problem: stable vortex solitons in uniform quadratic media, previously thought to require competing nonlinearities or external structure.","feed_headline":"A flat spin-orbit coupling stabilizes 3D vortex solitons","feed_subtitle":"Atomic-molecular condensates can now hold free-space vortex solitons, with over half their norm in the vortex components.","key_machinery":"The load-bearing ingredient is the planar Rashba coupling $\\lambda(\\partial_x-i\\partial_y)\\varphi_2$ and $\\lambda(\\partial_x+i\\partial_y)\\varphi_1$ acting on the atomic components, combined with the phase-matched quadratic coupling $-\\psi\\varphi_1\\varphi_2$. This pairing forces the molecular field to inherit the phase winding of the vortex atomic component, so the vortex topology is shared among all three fields. The conserved $z$-component of angular momentum $M_z$ ties the spin imbalance $|\\varphi_1|^2-|\\varphi_2|^2$ to the orbital angular momentum, providing the topological bookkeeping that makes the vortex-holding molecular component a natural partner. Stationary solutions are produced by imaginary-time propagation from SV and MM inputs, and the cascading approximation, which adiabatically eliminates $\\psi$ at large negative detuning, is used to map the system onto effective two-component cubic equations, explaining both the formation mechanism and the nonexistence above $\\alpha_c$.","core_discovery":"The central discovery is that the mean-field system of two atomic components $\\varphi_1,\\varphi_2$ and a molecular field $\\psi$, coupled by the quadratic three-wave interaction, supports stable free-space 3D solitons when a planar Rashba spin-orbit coupling acts on the atomic components alone. The solitons come in two symmetry classes: semi-vortices (SV), where $\\varphi_1$ has zero vorticity and $\\varphi_2$ and $\\psi$ carry unit circulation, and mixed modes (MM), where each component superposes vorticities $0$ and $\\pm1$ so that the total angular momentum vanishes. The paper reports existence up to a critical detuning $\\alpha_c>0$, explains that boundary with the cascading approximation, and verifies stability by imaginary-time preparation followed by real-time propagation, supported by the monotonic decrease of chemical potential with norm. The vortex-carrying components in the SV states account for more than half of the total norm, a distinguishing feature relative to earlier SOC-supported solitons.","pith_inferences":["A Bogoliubov–de Gennes linear-stability calculation around the stationary SV and MM solutions would sharpen the regime boundaries: if unstable modes with small growth rates exist, the blue stability regions in Fig. 3 would shrink to metastability on the numerical horizon.","The cascading reduction suggests that for large negative detuning the atomic subsystem alone obeys a two-component cubic SOC model with attraction proportional to $1/|\\alpha|$; comparing full soliton profiles with this reduced model's predictions would provide a quantitative test of the stabilization mechanism.","The same geometry may transfer to photonic systems with synthetic spin-orbit coupling and three-wave mixing, so that the optical-bullet counterpart would inherit the stabilization; the paper names this as a future direction but does not demonstrate it.","The conserved angular momentum expression contains a spin-imbalance term $\\frac12(|\\varphi_1|^2-|\\varphi_2|^2)$, so a targeted perturbation that changes the relative populations of the two atomic components should have a measurable effect on the vortex stability, which could be tested in simulations today."],"forward_implications":["Stable free-space 3D vortex solitons become available in quadratic nonlinear media without a competing cubic nonlinearity, removing the azimuthal-instability obstruction that previously forced vortex fragmentation.","The same planar SOC mechanism should be realizable in a $^{39}$K condensate with Raman-laser-induced SOC and photoassociation, with the quoted scales giving about 6500 atoms for $N=100$ and a physical evolution time near one second for $t=1000$.","The negative slope $d\\mu/dN$ persists across the studied parameter domains, so the predicted families should be identifiable as attractors in quench experiments that suddenly turn on the SOC and molecular coupling.","Both soliton families disappear above a critical mismatch: for the SV case the reported cutoff is $\\alpha>1.06$ at the studied parameters, while for the MM case the critical detuning grows with atom number.","The large vortex-norm share, exceeding 50%, means the angular-momentum-carrying part dominates the wavefunction, which is favorable for applications that require a strong vortex response."],"supporting_citations":[{"why":"Defines the SV and MM soliton types and the ansatz for spin-orbit-coupled binary BECs that this paper transplants into three dimensions with quadratic nonlinearity.","marker":"[63]"},{"why":"Shows that in the cubic SOC system 3D solitons are only metastable, providing the contrast that makes planar SOC in the quadratic system a new result.","marker":"[64]"},{"why":"Supplies the cascading approximation used to adiabatically eliminate the molecular field and to explain the existence boundary in detuning.","marker":"[9]"},{"why":"Establishes that the quadratic nonlinearity can support stable 3D fundamental solitons, against which this paper adds vortex states.","marker":"[10]"},{"why":"Documents the azimuthal instability of vortex solitons in quadratic media, the problem the paper claims to overcome.","marker":"[13]"},{"why":"Demonstrates stabilization of vortex solitons by competing quadratic and cubic nonlinearities, the baseline that the present paper removes the need for.","marker":"[37]"},{"why":"Provides the experimental Raman-laser scheme for planar SOC in Bose-Einstein condensates used for the proposed realization.","marker":"[70]"},{"why":"Provides the photoassociation mechanism for atomic-molecular conversion used to realize the quadratic coupling in the proposed setup.","marker":"[71]"}],"fun_headline_variants":["Stable 3D vortex solitons from spin-orbit coupling","Spin-orbit coupling tames 3D solitons in free space","Atomic-molecular condensates host stable 3D vortices","Half-vortex solitons stabilized by flat SOC","Planar SOC enables stable 3D vortex solitons"],"cache_read_input_tokens":14976,"weakest_assumption_plain":"The paper's stability verdict rests on watching imaginary-time-prepared solitons survive \"long-term real-time propagation\" without reporting the integration duration, perturbation size, or grid resolution, and without computing a linear-stability spectrum; if the numerical horizon is shorter than the growth time of an azimuthal instability, the blue stability regions would not represent true asymptotic stability.","fun_headline_variants_meta":{"raw":{"variants":["Stable 3D vortex solitons from spin-orbit coupling","Spin-orbit coupling tames 3D solitons in free space","Atomic-molecular condensates host stable 3D vortices","Half-vortex solitons stabilized by flat SOC","Planar SOC enables stable 3D vortex solitons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000591,"raw_usage":{"total_tokens":2764,"prompt_tokens":927,"completion_tokens":1837,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":1762}},"tokens_in":543,"tokens_out":1837,"duration_ms":14073,"temperature":1.0,"reasoning_tokens":1762,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:38:32.752715+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the stationary SV and MM solutions shown in Figs. 1 and 2, add a random phase perturbation with amplitude $\\epsilon$ between $10^{-4}$ and $10^{-2}$ of the local density, and evolve Eqs. (1)-(3) for at least an order of magnitude longer than the reported run; if the vortex ring fragments into several filaments or the angular momentum $M_z$ changes by more than numerical error, the stability claim is disproved. A cheaper quantitative check is the Bogoliubov spectrum of the linearized system around these stationary states: a positive imaginary part of any eigenfrequency would falsify asymptotic stability.","supporting_citations":[{"cited_title":"optical bullets","cited_arxiv_id":null,"evidence_quote":"Defines the SV and MM soliton types and the ansatz for spin-orbit-coupled binary BECs that this paper transplants into three dimensions with quadratic nonlinearity."},{"cited_title":"Sakaguchi and B","cited_arxiv_id":null,"evidence_quote":"Shows that in the cubic SOC system 3D solitons are only metastable, providing the contrast that makes planar SOC in the quadratic system a new result."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the cascading approximation used to adiabatically eliminate the molecular field and to explain the existence boundary in detuning."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that the quadratic nonlinearity can support stable 3D fundamental solitons, against which this paper adds vortex states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the azimuthal instability of vortex solitons in quadratic media, the problem the paper claims to overcome."},{"cited_title":"Mihalache, D","cited_arxiv_id":null,"evidence_quote":"Demonstrates stabilization of vortex solitons by competing quadratic and cubic nonlinearities, the baseline that the present paper removes the need for."},{"cited_title":"Dalfovo, S","cited_arxiv_id":null,"evidence_quote":"Provides the experimental Raman-laser scheme for planar SOC in Bose-Einstein condensates used for the proposed realization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the photoassociation mechanism for atomic-molecular conversion used to realize the quadratic coupling in the proposed setup."}],"review_version":1}