{"id":"f8cd8cf8-4bdc-4821-b01d-882d974548b3","arxiv_id":"2507.16111","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A linearly coupled two-component defocusing system supports spontaneous symmetry breaking with exact asymmetric states, dark solitons at g=3, and shifted-core vortices.","lead":"This paper shows that two overlapping waves that repel each other and are weakly mixed can spontaneously become unequal, producing asymmetric continuous waves, dark solitons, and vortices. The analytical formulas and simulations cover optical fibers and two-state Bose-Einstein condensates, giving concrete states that experiments could look for.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Modulational stability of asymmetric CW is asserted via exchange-of-stability but not computed away from the SSB threshold; finite-q MI on the asymmetric branch could invalidate the stability of SSB, DS, and vortex states.","rationale":"The reader's weakest_assumption correctly identifies the inheritance of modulational stability as the key unproven step in the analytic part of the paper. I independently rechecked the exact DS solution (39) and confirmed that, with the tanh derivative computed correctly (the second derivative carries a factor tanh), the solution does satisfy the stationary equations at g=3, so the analytical soliton construction is sound. The energy comparison (24) and the existence of the asymmetric CW (17)-(18) are also correct. The remaining load-bearing gap is that the stability of the asymmetric CW is never computed for parameters away from the bifurcation point; only the symmetric CW at the threshold is shown to be stable (Eqs. 32-33). The exchange-of-stability principle covers the zero-wavenumber mode, but not finite-q modulational instabilities, which are a real concern for multicomponent defocusing systems with cross-repulsion. If such an instability exists, the claimed stable DSs and vortices built on the asymmetric background would be unstable, directly affecting the central claim. The proposed concrete test settles this by computing the BdG spectrum at representative parameters. Until such a test is performed, the CONDITIONAL verdict is appropriate.","tokens_in":15251,"tokens_out":26263,"duration_ms":230699,"concrete_test":"Numerically compute the Bogoliubov–de Gennes spectrum for the exact asymmetric CW (17)-(18) at g=2, κ=0.45, k=-1.1: linearize Eqs. (2)-(3) around the stationary solution with perturbations (δu,δv,δu*,δv*) e^{γt+iqx}, assemble the 4×4 BdG matrix, and check whether any eigenvalue has Im(γ)>0 for q in [0,q_max]. Also perform real-time simulations with random initial perturbations of amplitude 10^{-2} on a periodic domain with L=100, integrating to t_final=10^4 to see if the CW remains uniform. Repeat for a second point such as g=1.2, κ=0.06, k=-1.01 (Eq. 48). If no instability appears, the inheritance claim is supported; if it appears, the CONDITIONAL verdict should be upgraded to REJECT.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III B derives the MI dispersion for the symmetric CW only at the SSB-bifurcation point (Eqs. 32-33), then states that the asymmetric CW 'inherits the stability'. This is an exchange-of-stability argument for the zero-wavenumber SSB mode, but it does not rule out finite-q modulational instabilities specific to the asymmetric branch for -k strictly above 2κ/(g-1). The asymmetric branch has both components nonzero with different densities, so its Bogoliubov spectrum differs from the symmetric branch; no computation or simulation of that spectrum is reported. Because all claims of stable DSs on the asymmetric background (Section IV A) and stable shifted vortices on the asymmetric background (Section V) presuppose the modulational stability of this CW, an undetected finite-q MI would undermine the central physical predictions. This is not a disagreement with consensus, but a gap between what is proved and what is claimed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a two-component defocusing Gross-Pitaevskii / nonlinear-Schrödinger system with linear coupling κ and cross-repulsion g, in one and two dimensions. It reports spontaneous symmetry breaking (SSB) of continuous-wave (CW) states for g > 1, with an exact asymmetric CW solution (Eqs. 17–18), an energy argument showing this state is energetically favored (Eq. 24), and a modulational-stability analysis at the SSB bifurcation point (Eqs. 32–33). It also presents an exact dark-soliton (DS) solution at g = 3 (Eq. 39), numerically found stable DSs for g ≠ 3, including states with 'inner immiscibility' in the DS core on a miscible background, and stable 2D vortex states with inter-component core shifts. The central claim is that these symmetry-broken states are stable and observable in defocusing nonlinear fibers and Rabi-coupled BECs.","tokens_in":15467,"tokens_out":6454,"duration_ms":70391,"significance":"If correct, the paper establishes a new and simple setting for spontaneous symmetry breaking in a free-space defocusing two-component system, with exact analytical control over the CW and one special DS solution. The analytic derivations are self-contained and contain no fitted parameters: the asymmetric CW solution, the energy comparison, the MI dispersion at threshold, and the g = 3 DS obtained by symmetry from a known domain-wall solution are all explicit and checkable. The predicted outer and inner immiscibility thresholds (Eqs. 50 and 55) are falsifiable. The numerical results, if backed by quantitative evidence, would provide convincing support for the DS and vortex stability claims. The paper is likely to be of interest to the nonlinear-optics and BEC communities.","major_comments":[{"comment":"The statement that the asymmetric CW 'inherits the stability' from the symmetric state is an exchange-of-stability argument for the q = 0 SSB mode only. The dispersion relation (32) is derived for the symmetric CW at the bifurcation point, and it does not rule out finite-wavenumber modulational instabilities on the asymmetric branch for -k strictly above 2κ/(g-1). Because the asymmetric branch has unequal component densities, its Bogoliubov spectrum differs from that of the symmetric state. This gap is load-bearing, since the stability of the DS in Fig. 1(a) and the vortex in Fig. 4 presupposes a stable asymmetric CW background. The authors should either compute the Bogoliubov spectrum of the asymmetric CW (a constant-coefficient linearization, which should be tractable analytically or numerically) or provide quantitative numerical evidence of stability for the parameter values used in the simulations, including the perturbation amplitudes and evolution times.","section":"Section III.B, Eq. (32)-(33) and text after Eq. (33)"},{"comment":"The numerical stability claims are not supported by quantitative details. No grid spacing, time step, integration scheme, or final simulation time is reported, and no convergence tests for the imaginary-time propagation are described. The statement in Section IV.A that 'systematic simulations of perturbed evolution' verify stability is therefore difficult to assess. This is especially important because the paper claims stability for a range of g and κ values (e.g., g = 0.1 and 0.6 in Fig. 2, and the vortex states in Fig. 5) and also asserts instability of the PT-potential DS without showing details. The authors should report the numerical parameters, describe the perturbation procedure, and provide at least one convergence check against an exact solution (e.g., the g = 3 DS (39) or the asymmetric CW (17)-(18)).","section":"Section IV (numerical methodology, Figs. 1-5)"}],"minor_comments":[{"comment":"The phrase 'cf. Eq. (21)' after Eq. (21) is self-referential and likely intended to refer to Eq. (14) or Eq. (20); please correct the cross-reference.","section":"Eq. (21)"},{"comment":"The description of the solid and dashed lines in Fig. 3 is inconsistent: the text states that the dashed line corresponds to Eq. (50), while the caption states that the solid line with crosses is given by Eq. (50). Please reconcile the text and caption.","section":"Section IV.A and Fig. 3 caption"},{"comment":"The variational ansatz (52) assumes the symmetric DS profile as the background for the inner-immiscibility calculation. The paper should state more explicitly that this approximation applies only to the symmetric-CW-background case and cannot describe the outer-immiscibility regime.","section":"Section IV.B, Eq. (52)"},{"comment":"The abstract and title contain typographical errors ('dark so litons', 'a nd', 'e.g.'); these should be corrected in the final version.","section":"Abstract and title"},{"comment":"The numerical scheme used for imaginary-time and real-time propagation is not described. A brief statement of the discretization (e.g., split-step Fourier or finite difference) and boundary-condition treatment would improve reproducibility.","section":"Section IV.A and V"}],"recommendation":"major_revision","confidential_remarks":"This is a well-written paper from experienced authors with a clean analytical core. The main technical gap is the lack of a stability analysis of the asymmetric CW branch away from the SSB threshold; this is a fixable issue but it is load-bearing for the DS and vortex claims. The numerical stability evidence also needs to be quantified. I recommend major revision, not rejection, because the central analytic results are sound and the missing analysis is within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid SSB paper in the defocusing two-component GP class, and the main exact results hold up. The asymmetric CW solution (17) is new, the energy argument (24) is clean, and the g=3 dark soliton (39) follows legitimately from the earlier DW solution by a discrete symmetry. The most interesting numerical observation is the 'inner immiscibility' — DS and vortex cores split while the background stays mixed — and the plots make that credible. There is no fitting, no circular reasoning, and the derivations are self-contained; the self-citation to [36] is appropriate because that is the solution they build on.\n\nThe weak points are exactly the ones you noted. The MI calculation covers only the symmetric CW at the bifurcation point. The asymmetric branch 'inherits stability' in the zero-mode sense, but the paper does not compute Bogoliubov modes for the asymmetric CW away from threshold, and the authors explicitly say the analysis is too cumbersome. That is a real gap between what is proved and what is claimed. In practice, the numerical DS and vortex states on the asymmetric background suggest stability, but that is indirect evidence, not a proof. I would ask the authors for a numerical MI sweep or a linear-stability check of the asymmetric CW for representative parameters before accepting.\n\nThe numerics also lack standard details: no convergence checks, no perturbation amplitudes, no final times for the real-time evolutions, no code or data. For a paper whose new states are mostly numerical, that is a reproducibility problem, not a fatal one, but it needs to be fixed. The VA for the inner-immiscibility boundary is a crude approximation, as they admit, so I would not over-credit it.\n\nWho should read it: anyone working on coupled defocusing BECs, nonlinear fibers, or SSB in multimode systems. It deserves a serious referee. The revisions I would request are manageable: supply a finite-q stability calculation for the asymmetric CW, give the numerical details, and maybe comment on whether the exchange-of-stability argument can be made rigorous with a quick perturbative calculation. I would probably cite the exact CW and the inner-immiscibility observation in my own work.","headline":"Exact asymmetric CW and a clever g=3 dark soliton make this a solid SSB paper, though the asymmetric-branch stability proof and numerical reproducibility need work.","tokens_in":15954,"tokens_out":4186,"would_cite":true,"duration_ms":44034,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35B32","35C08","37K40"],"pacs":["03.75.Mn","42.65.Tg","05.45.Yv"],"model":"deepseek-v4-flash","headline":"In two repulsively coupled defocusing waves, the symmetric state gives way to an exact asymmetric one above a critical density, and the broken-symmetry background hosts dark solitons and vortices with shifted cores.","keywords":["spontaneous symmetry breaking","linearly coupled Gross-Pitaevskii equations","dark solitons","vortices","immiscibility","Bose-Einstein condensates","optical fibers","modulational stability"],"falsifier":"Directly linearize the system (2)-(3) around the asymmetric CW state (17)-(18) at a chemical potential $-k$ above the threshold (20) and check whether any perturbation wavenumber $q$ yields $\\mathrm{Re}\\,\\gamma>0$; the companion dynamical test is to evolve the symmetric state with small noise at $n>n_{\\rm thr}$ and verify that it relaxes to the stationary asymmetric state (17) rather than breaking into domain walls or oscillations.","tokens_in":15032,"feed_emoji":"🌀","tokens_out":18196,"duration_ms":156648,"temperature":0.7,"pith_summary":"This paper studies two wave components that repel each other more strongly than they repel themselves, are linearly coupled, and live in a defocusing (self-repulsive) medium: the setting of a two-component Bose-Einstein condensate with Rabi coupling, or of two circular polarizations of light in a birefringent fiber. It establishes that the symmetric equal-amplitude continuous-wave state stops being the energy minimum once the total density exceeds $n_{\\rm thr}=2\\kappa/(g-1)$ with $g>1$; an exact asymmetric state with $uv=\\kappa/(g-1)$ then has lower energy and is modulationally stable. The broken-symmetry background supports dark solitons, including an exact one at $g=3$, and in two dimensions it hosts vortices whose two components have mutually shifted cores that break isotropy. A distinct finding is 'inner' immiscibility: even when the uniform background remains fully symmetric and mixed, the core of a dark soliton or vortex can split between the two components. If the stability results are right, these states should be directly observable in Rabi-coupled condensates and in defocusing nonlinear fibers.","feed_headline":"Cross-repulsion tips two coupled waves into asymmetric states","feed_subtitle":"The lower-energy asymmetric state supports dark solitons and vortices with shifted cores.","key_machinery":"The argument runs on closed-form solutions of the stationary equations. For the uniform states, adding and subtracting the two cubic equations and factorizing produces the exact asymmetric amplitudes (17)-(18) with the fixed product $u_{\\rm as}v_{\\rm as}=\\kappa/(g-1)$; the energy-density comparison (22)-(24) then proves that the asymmetric branch is the lower-energy one above the threshold. Modulational stability is handled by a Madelung-form perturbation of the symmetric state at the bifurcation point, which yields the biquadratic dispersion relation (32) whose four roots (33) are purely imaginary; the asymmetric branch is thus taken to inherit stability by exchange of stability. For localized structures, the exact $g=3$ dark soliton (39) is a hyperbolic-tangent connection between two asymmetric backgrounds, obtained by applying the substitution $\\{u,v,\\kappa,x\\}\\to\\{u,-v,-\\kappa,-x\\}$ to the known domain-wall solution (37). The 'inner immiscibility' of soliton and vortex cores on symmetric backgrounds is explained by a variational approximation built on the tanh/$\\cosh$ ansatz (52), which yields the critical coupling $\\kappa_c=(7g-1)n/24$ of Eq. (55).","core_discovery":"The paper's central claim is that the linearly coupled Gross-Pitaevskii system (2)-(3) with self-defocusing and cross-repulsion of relative strength $g>1$ supports spontaneous symmetry breaking of its continuous-wave states. Adding and subtracting the stationary algebraic equations yields the exact asymmetric solution (17)-(18), with component product $u_{\\rm as}v_{\\rm as}=\\kappa/(g-1)$, which exists above the density threshold $n_{\\rm thr}=2\\kappa/(g-1)$ of Eq. (20). The energy-density difference (24) shows this state lies below the symmetric one for every $n>n_{\\rm thr}$, and the modulational-stability analysis at the bifurcation point, whose four dispersion roots (33) are purely imaginary, is invoked to conclude that the asymmetric branch inherits stability. At the special value $g=3$, an exact dark-soliton solution (39) connects two mirror-image asymmetric backgrounds and obeys the inversion symmetry $\\varphi(-x)=-\\psi(x)$. For general $g$, numerically stable dark solitons are found on asymmetric backgrounds, and the paper also finds that dark solitons and two-dimensional $S=1$ vortices can exhibit a shift between the components in their cores even when the supporting continuous-wave background is the fully mixed symmetric state.","pith_inferences":["Because the stability of the asymmetric branch away from threshold is inferred rather than computed directly, a targeted Bogoliubov analysis of the exact solution (17)-(18) for $n>n_{\\rm thr}$ would close the main gap in the analytic argument and would either confirm or overturn the exchange-of-stability assumption.","The inner-immiscibility mechanism is a general one: where the two components have opposite signs, the linear-mixing energy term changes sign, so localized cores can demix even when the homogeneous background cannot; this suggests looking for the same effect in discrete nonlinear lattices and in spin-orbit-coupled condensates.","A concrete experimental test in the fiber setting ($g=2$) is to launch the polarization-symmetric state just above threshold and measure the output polarization; the paper predicts a reproducible asymmetry, while the competing prediction of modulational breakup would produce a fluctuating or pulsing output."],"forward_implications":["In a Rabi-coupled binary condensate with repulsive interactions and $g>1$, a uniform mixed state should spontaneously develop unequal component densities once the total density exceeds $2\\kappa/(g-1)$.","In a self-defocusing birefringent fiber ($g=2$), the asymmetric polarization state is modulationally stable, so the CW asymmetry should persist over long propagation distances rather than breaking into bright-soliton chains.","Dark solitons exist on the asymmetric background for general $g$, with an exact stable solution at $g=3$; their hallmark is zero-crossing points in both components, which ordinary symmetric dark solitons do not have.","The core splitting ('inner immiscibility') can also occur in dark solitons and unit-charge vortices supported by a fully miscible symmetric background, and even for $g\\le 1$, with the variational prediction $\\kappa_c=(7g-1)n/24$."],"supporting_citations":[{"why":"Prior demonstration of symmetry breaking in a trapped two-component system with repulsive interactions and linear coupling, which the present work extends to free space and to dark solitons and vortices.","marker":"[22]"},{"why":"Supplies the immiscibility condition $g>1$ for the uncoupled limit, the key premise for the SSB threshold.","marker":"[30]"},{"why":"Original exact domain-wall solution at $g=3$ and $\\kappa=0$, from which the dark-soliton construction starts.","marker":"[35]"},{"why":"Exact domain-wall solutions for the coupled system with $g=3$ and $\\kappa\\ne 0$, whose substitution yields the new dark-soliton solution (39).","marker":"[36]"},{"why":"Defines the optical-fiber realization (two circular polarizations with $g=2$) and the modulational-instability context the present work avoids by using defocusing nonlinearity.","marker":"[9]"},{"why":"Provides the cross-phase-modulation instability mechanism that motivates the stability analysis of the CW states.","marker":"[31]"},{"why":"Background on dark solitons and vortices in defocusing media, including expectations of instability for higher-charge vortices.","marker":"[28]"},{"why":"The imaginary-time propagation method used to construct all numerical stationary states in 1D and 2D.","marker":"[39, 40]"},{"why":"Review of variational methods underpinning the variational approximation that explains inner immiscibility and yields Eq. (55).","marker":"[43]"}],"fun_headline_variants":["Cross-repulsion breaks symmetry, shifting soliton and vortex cores","Exact asymmetric waves support dark solitons and vortices","Symmetry-broken bimodal waves host shifted cores in solitons and vortices","Lower-energy asymmetric waves emerge for strong cross-repulsion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the asymmetric wave branch stays modulationally stable for every density above the threshold, argued only from the symmetric state's stability at the bifurcation point rather than from a direct stability check of the asymmetric branch.","fun_headline_variants_meta":{"raw":{"variants":["Cross-repulsion breaks symmetry, shifting soliton and vortex cores","Exact asymmetric waves support dark solitons and vortices","Symmetry-broken bimodal waves host shifted cores in solitons and vortices","Lower-energy asymmetric waves emerge for strong cross-repulsion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000911,"raw_usage":{"total_tokens":4001,"prompt_tokens":1117,"completion_tokens":2884,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":733,"completion_tokens_details":{"reasoning_tokens":2807}},"tokens_in":733,"tokens_out":2884,"duration_ms":25135,"temperature":1.0,"reasoning_tokens":2807,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:17:59.275390+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly linearize the system (2)-(3) around the asymmetric CW state (17)-(18) at a chemical potential $-k$ above the threshold (20) and check whether any perturbation wavenumber $q$ yields $\\mathrm{Re}\\,\\gamma>0$; the companion dynamical test is to evolve the symmetric state with small noise at $n>n_{\\rm thr}$ and verify that it relaxes to the stationary asymmetric state (17) rather than breaking into domain walls or oscillations.","supporting_citations":[{"cited_title":"Sakaguchi and B","cited_arxiv_id":null,"evidence_quote":"Prior demonstration of symmetry breaking in a trapped two-component system with repulsive interactions and linear coupling, which the present work extends to free space and to dark solitons and vortices."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the immiscibility condition $g>1$ for the uncoupled limit, the key premise for the SSB threshold."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Original exact domain-wall solution at $g=3$ and $\\kappa=0$, from which the dark-soliton construction starts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Exact domain-wall solutions for the coupled system with $g=3$ and $\\kappa\\ne 0$, whose substitution yields the new dark-soliton solution (39)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the optical-fiber realization (two circular polarizations with $g=2$) and the modulational-instability context the present work avoids by using defocusing nonlinearity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the cross-phase-modulation instability mechanism that motivates the stability analysis of the CW states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Background on dark solitons and vortices in defocusing media, including expectations of instability for higher-charge vortices."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Review of variational methods underpinning the variational approximation that explains inner immiscibility and yields Eq. (55)."}],"review_version":1}