{"id":"d0c75948-3d11-416a-89a1-9682e98b130b","arxiv_id":"1908.07848","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A rotating two-component Bose-Einstein condensate can form non-overlapping parity-conjugate semicircular vortex lattices when the species have equal self-repulsion and strong mutual repulsion.","lead":"This paper uses computer simulations of a rotating binary Bose-Einstein condensate to show that, when the two atomic species push each other apart strongly, their vortex lattices settle into separate half-moon shaped regions that do not overlap. A generalist might read it to see how simple repulsive interactions can spontaneously break a perfectly round trap into two mirror-image halves.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Ground-state claim rests on a single imaginary-time descent with no energy comparison against symmetric, overlapping, or sheet competitors; the paper itself shows such states can be excited at g12=1100.","rationale":"The reader's weakest-assumption diagnosis is correct and is the most load-bearing issue in the paper. The central claim depends on the semicircular phase-separated vortex lattice being the global ground state of Eq. (9). Imaginary-time propagation can only certify a local minimum, and the paper provides no energy comparison against the natural competing states. The internal evidence strengthens this concern: at g12 = 1100 the paper explicitly states that the fully phase-separated semicircular states are excited states while the sheet structure is the ground state, demonstrating that the energy ordering between these topologies is nontrivial. At g12 = 2000, the same ordering is simply asserted. A concrete energy comparison would settle the issue. I do not see a more fundamental flaw: the equations are standard, the numerical method is established, the phase winding checks are appropriate, and the qualitative observation of non-overlapping vortex lattices is likely robust. The appropriate disposition remains CONDITIONAL, with the condition being an energy comparison against symmetric, overlapping, and sheet competitors. Since the reader already assigned CONDITIONAL, my stress-test does not change the verdict.","tokens_in":11315,"tokens_out":5245,"duration_ms":59766,"concrete_test":"At g1 = g2 = 1000, g12 = 2000, and fixed Ω (e.g., 0.6 and 0.95), run the same split-step imaginary-time solver to convergence from at least four initial states: (i) the published semicircular phase-separated state; (ii) a rotationally symmetric core-shell state with one component centered and the other annular; (iii) an overlapping state with both components filling the trap; (iv) the stripe/sheet state of Fig. 2(g)-(h) continued adiabatically to g12 = 2000. Evaluate the rotating-frame energy E in Eq. (9) for each converged solution. If the semicircular state is not strictly the lowest, the headline ground-state claim is false and the result must be reframed as a metastable spontaneously broken-symmetry state.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central assertion is that the semicircular phase-separated vortex lattices are spontaneously symmetry-broken ground states. The only evidence offered is imaginary-time propagation from a localized state with random phase, which is a local descent in the nonconvex energy landscape of Eq. (9). The paper states without supporting data: 'The fully phase-separated states are the ground states for all Ω' (Sec. 3, before Fig. 5). No energy of this state is compared with (i) a rotationally symmetric core-shell configuration, (ii) an overlapping binary vortex lattice, or (iii) the stripe/sheet state that the paper itself identifies as the ground state for g12 ≈ g1 = g2 (Fig. 2(g)-(h), and the sentence 'The fully phase-separated states in this case ... are excited states of higher energy'). At g12 = 2000 the same comparison is absent, yet the ground-state wording is retained. This matters because if the semicircular states are merely metastable, the phrase 'spontaneous symmetry breaking' overstates the finding; the robust, reproducible content would then be the existence and dynamical stability of non-overlapping vortex lattices, not their ground-state status. The phase checks in Fig. 7 establish unit circulation but do not address global energy minimality.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a rapidly rotating, harmonically trapped quasi-2D binary Bose-Einstein condensate with repulsive intra- and inter-species interactions, described by the coupled Gross-Pitaevskii equations (6) and (7). Using imaginary-time propagation from a localized random-phase initial state, the author reports that for sufficiently strong inter-species repulsion the two condensates phase-separate into non-overlapping regions, each carrying a vortex lattice. For equal intra-species interactions (g1 = g2 = 1000, g12 = 2000) the two domains are parity-conjugate semicircles with triangular vortex lattices; for nearly equal interactions (g1 = 1000, g2 = 900, g12 = 2000) the separation is complete but the shapes differ. The paper also computes a numerical critical inter-species coupling, compares it with the uniform-mixture threshold sqrt(g1 g2), and demonstrates dynamical stability of the symmetric case by real-time evolution.","tokens_in":11549,"tokens_out":8310,"duration_ms":89857,"significance":"If the ground-state claim is substantiated, the finding is significant: completely non-overlapping vortex lattices would directly address the experimental difficulty of observing vortices in overlapping binary condensates, and the predicted parity-conjugate semicircular shapes provide a clear signature of spontaneous rotational symmetry breaking. The paper has concrete strengths: the numerical method is explicitly specified, the vortex identification is checked through phase winding (Fig. 7), the symmetrical case is tested by long real-time propagation (Fig. 9), and the analytical threshold sqrt(g1 g2) is used as a benchmark rather than fitted. The central weakness is that the global-minimum status of the reported states is asserted rather than demonstrated, which matters because the abstract frames the result as spontaneous symmetry breaking of the ground state.","major_comments":[{"comment":"The statement 'The fully phase-separated states are the ground states for all Ω' is not supported by any energy comparison. The manuscript nowhere compares the energy of the semicircular phase-separated state with (i) a rotationally symmetric core-shell configuration, (ii) an overlapping binary vortex lattice, or (iii) the stripe/sheet states that the paper itself identifies as the ground state for g12 = 1100 (Fig. 2(g)-(h), and the sentence 'The fully phase-separated states in this case ... are excited states of higher energy'). Because the energy functional (9) is nonconvex, convergence of imaginary-time propagation from a single localized random-phase initial state is a local descent and does not by itself establish global minimality. This issue is load-bearing: the abstract's 'spontaneous symmetry breaking' and the ground-state phrasing depend on the semicircular states being true ground states, not merely stable excited states. Please provide energy differences, or at least total energies from Eq. (9), for the phase-separated state against these competitors at representative Ω and g12, and state explicitly for which parameter range the phase-separated state is the lowest-energy state.","section":"Sec. III, paragraph before Fig. 5"},{"comment":"The 'numerically computed gc12' plotted in Fig. 3 is never operationally defined. The text states that for g12 larger than the numerical gc12 a phase-separated vortex lattice is obtained, but the reader is not told what observable distinguishes a phase-separated vortex lattice from a stripe/sheet state, what tolerance is used, or how the critical value is extracted from the numerical runs. In addition, the claim that 'The numerical critical value gc12 is found to be independent of Ω (< 1)' is asserted without showing results for more than one Ω; the figure itself displays only Ω = 0.6. These details matter because the paper's parameter-domain recommendation (Sec. IV, g1g2/g12^2 ⪅ 0.75) and the phase-separation criterion rest on this numerical curve. Please specify the extraction criterion, the convergence criterion, and provide the Ω-dependence data or qualify the claim.","section":"Sec. III, Fig. 3"}],"minor_comments":[{"comment":"The caption entries 'g12−1200' and 'g12−1400' should presumably read 'g12 = 1200' and 'g12 = 1400'; the printed minus sign is confusing.","section":"Fig. 4 caption"},{"comment":"The phrase 'arranged usually in a Abrikosov triangular lattice' should read 'arranged usually in an Abrikosov triangular lattice'.","section":"Sec. I"},{"comment":"The statement that the asymmetric vortex lattice is 'only weakly stable' is not supported by any displayed data; either provide the corresponding real-time evolution or soften the claim.","section":"Sec. IV"},{"comment":"The overlaid contour plots of the two components would benefit from a legend, color definition, or explicit labeling of which grayscale/contour level corresponds to component 1 and component 2.","section":"Figs. 6 and 8"},{"comment":"The displayed condition appears typeset as 'g1g2/g2 12 < 1'; the intended condition is g12^2 > g1 g2. Please ensure the equation is unambiguous in the final version.","section":"Eq. (1)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a single-author numerical study using the author's own published OMP codes (Refs. [27,32,34]); this is legitimate, but independent verification or a clearer statement of numerical convergence would increase confidence. The main issue is the unsupported global-ground-state claim; if the author adds energy benchmarks against the relevant competitors and defines the numerical gc12 criterion, the paper could become acceptable. I recommend major revision rather than rejection because the existence, phase-winding, and dynamical-stability results are credible and the missing energy comparison is a targeted, feasible addition."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nHere is my take on Adhikari's binary BEC vortex paper. The useful core is real: for strong inter-species repulsion and equal intra-species strengths, the two condensates phase-separate into disjoint semicircular regions, each hosting a triangular vortex lattice, and the two lattices are parity conjugates. That geometry is genuinely new compared with the stripe/sheet/honeycomb configurations in Refs. [11,19,24]. The numerics are careful in the standard imaginary-time Crank-Nicolson style: phase plots confirm unit circulation, vortex counts scale with Ω, and a real-time run with a small frequency kick shows the symmetric lattice is dynamically stable. No code or data are shipped, but the method is well established and the qualitative pattern is probably robust.\n\nThe soft spot is exactly what the stress-test note flags. The paper claims these fully phase-separated states are ground states for all Ω, but the only evidence is a single imaginary-time descent from a random-phase localized state. There is no energy comparison against a rotationally symmetric core-shell state, an overlapping binary lattice, or the sheet/stripe states that the paper itself identifies as lower-energy for g12 close to g1=g2, see Fig. 2(g)-(h) and the sentence before Fig. 3. At g12=2000 the energy ordering might be correct, but it is not demonstrated. \"Spontaneous symmetry breaking\" is fair for the shapes themselves; \"ground state\" is stronger than the evidence. A related minor issue: no quantitative overlap measure; \"no overlap\" is judged by eye from contour plots.\n\nThe analytical threshold g12^c = sqrt(g1 g2) is imported from the uniform mixture and compared with the numerical boundary. That is a reasonable benchmark, not a fitted parameter, and the paper does note it is not rigorous in the trapped quasi-2D setting.\n\nWho gets value: anyone working on vortex lattices in binary or spinor condensates, especially experimental groups who want a geometry where the two components' vortices are spatially separated for imaging. It is an incremental but useful contribution, not a paradigm shift.\n\nRecommendation: it deserves a serious referee. The referee should ask for either (1) energy comparisons against the obvious competitor states, or (2) softening the ground-state wording to \"stationary states found by imaginary-time relaxation.\" Either way, the core observation is likely publishable.","headline":"Useful new geometry for phase-separated vortex lattices, but the ground-state claim outruns the evidence.","tokens_in":12120,"tokens_out":2252,"would_cite":true,"duration_ms":24549,"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":"This paper claims that rapidly rotating binary Bose-Einstein condensates with equal intra-species repulsion and strong inter-species repulsion settle into two completely separated, mirror-image semicircular vortex lattices, spontaneously…","keywords":["binary Bose-Einstein condensate","vortex lattice","Gross-Pitaevskii equation","spontaneous symmetry breaking","phase separation","rapid rotation","quasi-two-dimensional"],"falsifier":"Compute the rotating-frame energy (9) of the semicircular phase-separated state and of a rotationally symmetric annular phase-separated state (or an overlapping vortex lattice) at the same parameters, e.g., $g_1=g_2=1000$, $g_{12}=2000$, $\\Omega=0.6$; the central claim fails if the symmetric or overlapping configuration has lower energy. Experimentally, image both components after a rapid rotation ramp: for $g_1=g_2$ the two density profiles should be exact mirror images with no overlap, and any significant deviation would falsify the predicted parity-conjugate ground state.","tokens_in":11059,"feed_emoji":"🌀","tokens_out":16411,"duration_ms":132634,"temperature":0.7,"pith_summary":"The paper aims to show that a rapidly rotating two-component Bose-Einstein condensate, with equal (or nearly equal) intra-species repulsion and strong inter-species repulsion, spontaneously breaks the circular symmetry of its trap and settles into two completely separated vortex lattices. For exactly equal self-interaction strengths, the two lattices occupy identical semicircular regions that are parity conjugates of each other; for nearly equal strengths, the regions are still fully separated but have different shapes. The claim matters because non-overlapping vortex lattices in different spatial regions are far easier to observe experimentally than overlapping ones, and would allow direct study of how atomic interactions shape vortex-lattice order. The evidence is numerical solution of the mean-field Gross-Pitaevskii equation for a quasi-two-dimensional rotating binary condensate.","feed_headline":"Binary condensate splits into mirror-image vortex semicircles","feed_subtitle":"With equal self-repulsion and strong mutual repulsion, each component's vortices live on its own half of the trap.","key_machinery":"The load-bearing object is the coupled quasi-two-dimensional Gross-Pitaevskii system in the rotating frame, Eqs. (6)-(7), derived from the 3D rotating-frame equations by integrating out the tightly confined z-direction, together with its energy functional (9). The vortex lattices are stationary states of this functional, found by imaginary-time propagation with a random-phase-modulated initial state using a split-step Crank-Nicolson scheme. Two inequalities organize the physics: the uniform-mixture phase-separation condition $g_1g_2/g_{12}^2<1$ (Eq. 1), and the numerically refined critical value $g_{12}^c\\approx \\sqrt{g_1 g_2}$ that separates fully separated vortex-lattice states from stripe/sheet structures. The parity-conjugate semicircular geometry for $g_1=g_2$ is the signature that spontaneous symmetry breaking has occurred.","core_discovery":"The central discovery is that complete spatial separation of the two condensate components survives rapid rotation and organizes each component's vortices into its own lattice, with zero overlap between the vortex patterns. When $g_1=g_2$ and $g_{12}$ is sufficiently large, the ground state consists of two semicircular condensates, one the parity conjugate of the other, each carrying a triangular vortex lattice; when $g_1$ and $g_2$ are close but unequal, phase separation is still complete but the two lattices have different shapes and vortex counts. The paper further identifies a numerical critical inter-species strength $g_{12}^c$, independent of angular frequency $\\Omega$ and close to $\\sqrt{g_1 g_2}$, above which fully separated vortex lattices form; below it, stripe or sheet structures appear instead. The symmetric separated state is found to be dynamically stable under real-time evolution after a small change in $\\Omega$, while the asymmetric state is only weakly stable.","pith_inferences":["The mirror-image geometry suggests a concrete mechanism: with equal self-repulsion, swapping the two components is an exact symmetry of the energy, so a straight interface through the trap center costs the same energy for either assignment; the simulations realize this as two parity-conjugate semicircles.","The claimed $\\Omega$-independence of $g_{12}^c$ is a testable prediction: measuring the onset of complete phase separation in a rotating binary condensate across a range of rotation frequencies should produce a flat critical curve, while a strong $\\Omega$-dependence would mean rotation renormalizes the effective inter-species coupling beyond the mean-field picture.","The parity-conjugate relationship offers a built-in experimental calibration: at $g_1=g_2$, any imaging asymmetry between the two components can be separated from real physics because the ground state itself should be exactly mirror-symmetric.","The same non-overlap property should allow direct experimental counting of vortices in each component, turning vortex-lattice geometry into a sensitive probe of the ratio $g_{12}^2/(g_1 g_2)$."],"forward_implications":["For $g_1=g_2$ and $g_{12}$ above the numerical critical value, the two vortex lattices occupy disjoint semicircular domains that are parity conjugates, so vortices of one component never overlap the other component's density.","When $g_1$ and $g_2$ are close but unequal, phase separation remains complete but the two lattices have different shapes and, at higher rotation, different vortex counts (e.g., 21 = 11 + 10 at $\\Omega=0.7$ for $g_1=1000$, $g_2=900$).","The numerical critical inter-species strength for fully phase-separated vortex lattices is independent of rotation frequency $\\Omega$ and is well approximated by $\\sqrt{g_1g_2}$, refining the uniform-mixture criterion $g_1g_2/g_{12}^2<1$.","The symmetric phase-separated vortex lattice survives real-time evolution for 200 time units under a small change in $\\Omega$ (0.88 to 0.89), indicating dynamical stability; the asymmetric lattice is only weakly stable.","Below the critical $g_{12}$, rotating binary condensates form stripe or sheet structures rather than separated vortex lattices, so the fully separated regime is a distinct phase accessible by tuning scattering lengths via Feshbach resonance."],"supporting_citations":[{"why":"This supplies the phase-separation condition $g_1g_2/g_{12}^2<1$ (Eq. 1) that guides the parameter search and the analytical estimate $g_{12}^c=\\sqrt{g_1g_2}$.","marker":"[26]"},{"why":"This provides the rotating-frame formulation, the quantized-vortex picture, and the $\\Omega<\\omega$ mean-field validity limit that the model relies on.","marker":"[6]"},{"why":"This gives the coupled Gross-Pitaevskii equations (2)-(3) for a rotating binary condensate from which the quasi-2D model is derived.","marker":"[29]"},{"why":"This provides the quasi-2D reduction that integrates out the tight z-direction and yields the two-dimensional equations (6)-(7).","marker":"[30]"},{"why":"This supplies the split-step Crank-Nicolson imaginary-time method used to converge to stationary vortex-lattice states.","marker":"[27]"},{"why":"This supplies the OpenMP-adapted code used to simulate vortex lattices in the rapidly rotating condensate.","marker":"[34]"},{"why":"This describes the random-phase-modulated initial state that seeds vortex generation in the imaginary-time simulation, used here alongside [34].","marker":"[13]"},{"why":"This documents the stripe and sheet structures in rotating binary condensates that the paper must exclude; the fully separated vortex-lattice regime is defined against these.","marker":"[24]"}],"fun_headline_variants":["Mirror-image vortex halves from a binary BEC","Binary BEC splits into mirror-image vortex semicircles","Zero overlap: two vortex lattices in one rotating BEC","Phase-separated vortex lattices in a binary BEC, no overlap","Each condensate gets its own vortex lattice, no mixing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The picture rests on the assumption that the imaginary-time simulation, started from a random-phase wavefunction, converges to the true global ground state of the rotating-frame energy, so the mirrored semicircular states are genuine ground states rather than metastable local minima; the paper does not compare their energy against rotationally symmetric or overlapping configurations.","fun_headline_variants_meta":{"raw":{"variants":["Mirror-image vortex halves from a binary BEC","Binary BEC splits into mirror-image vortex semicircles","Zero overlap: two vortex lattices in one rotating BEC","Phase-separated vortex lattices in a binary BEC, no overlap","Each condensate gets its own vortex lattice, no mixing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000948,"raw_usage":{"total_tokens":4050,"prompt_tokens":950,"completion_tokens":3100,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":3015}},"tokens_in":566,"tokens_out":3100,"duration_ms":22072,"temperature":1.0,"reasoning_tokens":3015,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:32:26.396777+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the rotating-frame energy (9) of the semicircular phase-separated state and of a rotationally symmetric annular phase-separated state (or an overlapping vortex lattice) at the same parameters, e.g., $g_1=g_2=1000$, $g_{12}=2000$, $\\Omega=0.6$; the central claim fails if the symmetric or overlapping configuration has lower energy. Experimentally, image both components after a rapid rotation ramp: for $g_1=g_2$ the two density profiles should be exact mirror images with no overlap, and any significant deviation would falsify the predicted parity-conjugate ground state.","supporting_citations":[{"cited_title":"Mineev, Zh","cited_arxiv_id":null,"evidence_quote":"This supplies the phase-separation condition $g_1g_2/g_{12}^2<1$ (Eq. 1) that guides the parameter search and the analytical estimate $g_{12}^c=\\sqrt{g_1g_2}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This provides the rotating-frame formulation, the quantized-vortex picture, and the $\\Omega<\\omega$ mean-field validity limit that the model relies on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This gives the coupled Gross-Pitaevskii equations (2)-(3) for a rotating binary condensate from which the quasi-2D model is derived."},{"cited_title":"Salasnich, A","cited_arxiv_id":null,"evidence_quote":"This provides the quasi-2D reduction that integrates out the tight z-direction and yields the two-dimensional equations (6)-(7)."},{"cited_title":"Muruganandam, S","cited_arxiv_id":null,"evidence_quote":"This supplies the split-step Crank-Nicolson imaginary-time method used to converge to stationary vortex-lattice states."},{"cited_title":"Kishor Kumar, V","cited_arxiv_id":null,"evidence_quote":"This supplies the OpenMP-adapted code used to simulate vortex lattices in the rapidly rotating condensate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This describes the random-phase-modulated initial state that seeds vortex generation in the imaginary-time simulation, used here alongside [34]."},{"cited_title":"Kasamatsu, M","cited_arxiv_id":null,"evidence_quote":"This documents the stripe and sheet structures in rotating binary condensates that the paper must exclude; the fully separated vortex-lattice regime is defined against these."}],"review_version":1}