{"id":"cff31c27-d5e6-475f-8e59-499290808c14","arxiv_id":"2507.16277","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In 3D two-component Rydberg-dressed BECs, tuning the inter-component soft-core range relative to the intra-component one yields FCC, SC, segregated tubular/planar, and immiscible ground states.","lead":"A theoretical study maps the crystal-like and layered states of a mixture of two ultracold gases whose atoms are dressed to interact over long distances. It predicts simple cubic, tubular, and planar density patterns that can form spontaneously from a uniform initial state.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Ground-state phase diagram rests on a variational search restricted to cubic lattices; absence of convergence tests or non-cubic competitors leaves the global-minimum claim unproven.","rationale":"The paper's central claim is that the predicted FCC/SC/ST/SP (and FCC/BCC/PW at equal radii) are ground states, and that these states are dynamically accessible. The most load-bearing requirement for this claim is that the reported structures are global energy minima of the coupled GPE. The available evidence is a Gaussian variational ansatz restricted to cubic Bravais lattices and a small number of GPE simulations with undisclosed numerical parameters and no convergence analysis. This leaves open the possibility that a different lattice (e.g., a triangular tube lattice for the ST phase) or a periodic box of different shape would produce a lower-energy state, which would change the phase diagram and the ground-state language. The reader's weakest assumption identifies exactly this gap, and I find no other concern that is more load-bearing. The concern is about missing evidence, not a demonstrated error, so the reader's CONDITIONAL verdict is appropriate and does not need adjustment.","tokens_in":21333,"tokens_out":12871,"duration_ms":130508,"concrete_test":"At a representative point in the ST region (e.g., R_c^(12)/R_c = 0.95, as/Rc = 2.3e-3), extend the variational energy calculation to a 2D triangular lattice of Gaussian tubes (and, in the SP region, to planes with a finite in-plane modulation), and compare the minimized energy against the SC/ST/SP energies from the current ansatz. Alternatively, run imaginary-time GPE in a rectangular box with an aspect ratio commensurate with a triangular lattice and with random initial conditions; if a triangular-tube structure emerges with lower energy than the square-tube ST state, the reported ground-state sequence is incomplete.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The phase diagram in Fig. 2(a) is built from imaginary-time GPE and a variational ansatz that only considers Gaussian profiles on SC, FCC, and BCC lattices, plus the ST/SP limits with one or two divergent widths (Supplemental Sec. II). The accessible text gives no box-size, grid-resolution, or initialization details for the GPE runs, and no comparison against non-cubic Bravais lattices (HCP, diamond, hex-tetragonal) or aperiodic patterns. Because a cubic periodic box cannot represent a hexagonal or other non-cubic lattice without distortion, the numerics may bias the system toward the reported cubic arrangements. This directly undermines the central claim that the states labeled FCC, SC, ST, SP are the true ground states: if a triangular arrangement of tubes (for ST) or a different stacking of planes has lower energy at the same parameters, the phase boundaries and even the existence of the reported sequence are not established. The variational method is, by the authors' own description, only a 'rough estimation' (Supplemental Sec. II), so it cannot certify global optimality on its own.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a three-dimensional two-component Bose-Einstein condensate with Rydberg-dressed soft-core interactions. The authors derive an effective mean-field model from a coupled light-atom Hamiltonian in the Supplemental Material, including explicit formulas for the soft-core potential parameters and a table of feasible atomic configurations. Within this model, they combine imaginary-time GPE evolution and a Gaussian variational ansatz restricted to SC, FCC, BCC, ST, and SP density profiles. The central claim is a ground-state phase diagram in which, for imbalanced inter-component versus intra-component blockade radii, the sequence FCC → SC → SP and FCC/SC → ST → SP occurs as the s-wave scattering length increases, while for equal radii the single-component-like FCC → BCC → plane-wave sequence is recovered. The paper also presents real-time GPE simulations showing spontaneous formation of SC, ST, and SP structures.","tokens_in":21587,"tokens_out":6598,"duration_ms":72606,"significance":"If the phase diagram is correct, the paper identifies a new family of three-dimensional supersolid and smectic-like phases in an experimentally accessible platform and strengthens the analogy to pasta phases in neutron-star matter. Strengths include the microscopic derivation of the effective soft-core potentials in Supplement I, the concrete Cs/Rb implementation schemes in Table I, the variational energy landscapes that make the ST/SP mechanisms transparent, and the independent roton-wavelength estimate of the SP periodicity. The dynamical-formation simulations are a useful addition. However, because the ground-state claim rests on a restricted variational basis and on numerical searches whose convergence and completeness are not documented, the significance is currently conditional on strengthening the global-minimum evidence.","major_comments":[{"comment":"The global-minimum claim is not established by the evidence presented. The Gaussian variational ansatz is explicitly described as \"a relatively rough estimation\" (Supplemental Sec. II) and is restricted to SC, FCC, BCC, ST, and SP trial profiles; the imaginary-time GPE is run in a periodic box that cannot represent non-cubic Bravais lattices or aperiodic patterns without distortion. No box-size or aspect-ratio convergence tests and no comparison against alternative structures (e.g., HCP, diamond, or larger supercells allowing incommensurate modulations) are reported. Since the abstract and conclusion state that these are the ground states and the solid phase boundaries in Fig. 2(a) are variational, the present evidence only establishes optimality within the considered ansatz. I request either (i) an expanded numerical search over a wider class of candidate unit cells with convergence checks, or (ii) a revised, more cautious claim of stable phases within the considered variational and numerical space, with corresponding changes to the abstract and conclusion.","section":"Ground-state phase diagram (Fig. 2(a)) and Supplemental Material Sec. II"},{"comment":"The quantitative relationship between the numerical GPE results and the variational phase boundaries is not documented. The text states that the solid lines are in \"reasonable agreement\" with the accurate full GPE simulation results, but the figure shows numerical points only along one path, and the manuscript does not report the numerical energies per particle, the numerical critical scattering lengths, or the grid resolution used. Without this information it is not possible to judge the accuracy of the variational phase boundaries, which carry the phase diagram, or the location of the claimed multicritical point at as/Rc ≈ 1.5 × 10^-3. Please provide representative numerical energy comparisons and critical values along cuts such as Rc^(12)/Rc = 0.95, and state explicitly what the dashed line in Fig. 2(a) denotes.","section":"Fig. 2(a) and the numerical 'Ground-state phase diagram'"},{"comment":"The roton-wavelength comparison in Fig. 2(b) mixes parameters: the caption fixes as/Rc = 0.5 × 10^-3 for the FCC state and as/Rc = 2.3 × 10^-3 for the SC, ST, and SP states. The comparison of arot with the SP lattice constant is therefore not made at a single set of interaction parameters across the phases, and it is unclear whether arot is computed at the same value of as for all Rc^(12)/Rc. To make this consistency check meaningful, the authors should plot arot and the lattice constants along the same parameter path used for the phase boundaries, and specify explicitly which value of as each curve corresponds to.","section":"Fig. 2(b) and the roton-wavelength comparison"}],"minor_comments":[{"comment":"The adiabatic elimination of the third-order correlators in Eqs. (6)–(8) neglects several terms without stating the ordering assumption; a sentence identifying the small parameters (e.g., ratios of Rabi frequencies to detunings and the dilute-gas mean-field conditions) would help readers assess the validity of the derivation.","section":"Supplemental Material Sec. I"},{"comment":"The phrase \"second-order phase transition\" for the SC→ST and ST→SP transitions appears without an order parameter or quantitative numerical evidence; please either define the order parameter and show its behavior across the transitions or describe these as continuous-looking transitions.","section":"Ground-state phase diagram"},{"comment":"Please label which species is shown in each density panel and define the color scale; currently the figure is difficult to read in the accessible version.","section":"Fig. 4"},{"comment":"The axis labels and the meaning of the dashed line and of the label \"cf. Fig. 3\" are not explained in the caption; please define the vertical and horizontal axes explicitly, including the scaling of as by 10^-3.","section":"Fig. 2(a)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of the journal and the microscopic derivation in the Supplemental Material is a genuine strength. The main risk is precisely the one identified in the stress-test note: the global-minimum claim is supported only by a restricted variational ansatz and by numerical simulations whose box-size convergence and candidate-lattice completeness are not documented. This is fixable either by adding a broader numerical search or by softening the ground-state language to 'stable phases found within the considered structures.' Because the latter is a substantive but tractable revision, I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid mean-field prediction paper with a genuinely new result—the 3D binary RdBEC phase diagram with SC, ST, SP, and immiscible phases—but the ground-state claim is stronger than the evidence. The variational ansatz only explores cubic lattices, and the numerical GPE runs are shown along a limited parameter path, so the global-minimum statement needs a qualification.\n\nWhat is actually new: the step from 2D binary or 3D single-component systems to 3D binary with unequal inter- and intra-component blockade radii is not incremental. The SC state (two interpenetrating FCC lattices), the segregated tubular and planar states, and the BCC at equal radii are concrete additions. The mechanism—roton instability of the out-of-phase branch—is physically sensible, and the roton-wavelength match to the SP lattice constant in Fig. 2(b) is a good independent check. The Supplemental derivation of the effective soft-core potentials, with concrete Rb and Cs level schemes and ARC-calculated C6 coefficients, is useful and makes this more than a toy model. The paper is also honest about dissipation and cites recent experiments on long-lived RdBECs.\n\nWhere it is soft: First, the phase diagram in Fig. 2(a) is almost entirely variational. The ansatz restricts the search to Gaussian profiles on SC/FCC/BCC plus the ST/SP limits obtained by letting widths diverge on those lattices. That excludes triangular tube arrangements, hexagonal close packing, or any aperiodic structure, and a cubic periodic box in the GPE runs would bias against non-cubic order. The Supplement itself calls the method 'relatively rough,' so the boundaries in Fig. 2(a) should be read as variational estimates, not certified transitions. This is fixable with box-size convergence checks, a couple of non-cubic competitor runs, and a softer statement about global optimality. Second, the GPE numerics appear only along one dashed line; grid resolution, box size, and initialization details are not given. Third, the 'new BCC' claim at equal radii deserves a careful check against existing single-component results—[23] may already contain a BCC transition. These are moderate issues, not fatal ones.\n\nWho it's for: anyone working on Rydberg-dressed BECs, 3D supersolids, or the neutron-star pasta analogy. It will get cited. It deserves a serious referee rather than a desk reject; the referee should ask for the missing numerical documentation and a less strong ground-state label.\n\nRecommendation: send it to review. With convergence tests and a non-cubic lattice check, this could be a solid reference for binary RdBEC phases.","headline":"A solid extension of RdBEC phase diagrams to 3D binary mixtures with real new phases, but the ground-state claim outruns the variational search.","tokens_in":22107,"tokens_out":4146,"would_cite":true,"duration_ms":45717,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two-component Rydberg-dressed Bose-Einstein condensates with imbalanced inter- and intra-component blockade radii have a ground-state sequence that goes from face-centered-cubic crystals through segregated tubes to segregated planes as…","keywords":["Rydberg-dressed Bose-Einstein condensate","binary Bose mixture","soft-core interaction","supersolid","ground-state phase diagram","segregated planar and tubular states","face-centered cubic lattice"],"falsifier":"At parameters where the paper predicts an SC ground state (for example $R_c^{(12)}/R_c = 0.8$ and $a_s/R_c = 1.0 \\times 10^{-3}$), run imaginary-time Gross-Pitaevskii evolution in boxes of at least two different sizes and test trial densities for hexagonal and tetragonal lattices; if any alternative structure has lower energy, or the SC state is not reproduced, the phase diagram needs revision.","tokens_in":21130,"feed_emoji":"⚛️","tokens_out":12622,"duration_ms":114810,"temperature":0.7,"pith_summary":"Two-component Rydberg-dressed Bose-Einstein condensates in three dimensions, where the inter-species and intra-species soft-core repulsions have different ranges, are shown to order into a sequence of supersolid phases as the s-wave scattering length grows. At weak contact interactions the two species crystallize into interpenetrating face-centered-cubic lattices whose combined density is simple cubic. Stronger contact repulsion lowers the dimensionality of the density modulation, first to segregated tubular patterns and then to alternating planar sheets. When the inter-component blockade radius exceeds the intra-component one, the mixture phase-separates instead. The ordered states are not only energy minima; real-time evolution from a uniform noisy condensate reaches them within tens of milliseconds, which makes the phases experimentally accessible.","feed_headline":"Imbalanced interactions drive Rydberg BECs from FCC crystals to planes","feed_subtitle":"Varying the scattering length lowers the symmetry-breaking dimension from 3D lattices to tubes to sheets.","key_machinery":"The central object is the coupled Gross-Pitaevskii energy functional for two components interacting through contact terms and soft-core potentials $V_{\\alpha\\beta}(r) = \\tilde{C}_6 / (1 + (r/R_c^{(\\alpha\\beta)})^6)$, with the ratio $R_c^{(12)}/R_c$ as the key control parameter. The argument is carried by a Gaussian variational ansatz in which each lattice site is an anisotropic Gaussian with widths $\\sigma_x$, $\\sigma_y$, $\\sigma_z$; the divergence of one or two widths signals the transition from three-dimensional crystals to segregated tubes or planes. The roton minimum of the lower Bogoliubov branch $\\omega_-(k)$ provides the length scale of the density modulation, and imaginary-time evolution in a periodic box confirms the variational phase boundaries.","core_discovery":"The paper claims that a 3D binary Rydberg-dressed BEC with symmetric intra-component interactions but an imbalanced inter-component blockade radius, $R_c^{(12)} \\neq R_c$, has a ground-state phase diagram controlled by the s-wave scattering length $a_s$. For $R_c^{(12)} < R_c$ and small $a_s$, each component forms a face-centered-cubic lattice displaced by half a lattice constant from the other, so the two-species density is a simple cubic crystal. Increasing $a_s$ makes large local densities energetically costly, and the system responds not by melting into a uniform state but by reducing the dimensionality of its symmetry breaking: it first forms segregated tubes, where the density is modulated in a plane and uniform along the tube axis, and then segregated planes, modulated along one direction only. The modulation period of these lower-dimensional states is set by the roton minimum of the out-of-phase Bogoliubov branch of the miscible mixture. When $R_c^{(12)} > R_c$ the two species phase-separate, and at $R_c^{(12)} = R_c$ the model reduces to a single-component RdBEC with a first-order FCC-to-BCC transition followed by a plane-wave state. The paper also shows that all these states can be reached dynamically from a homogeneous condensate seeded with noise.","pith_inferences":["If the variational ansatz were extended to other Bravais lattices (for example hexagonal close-packed or tetragonal) at the same parameters, some of the phase boundaries could move, but the qualitative dimension-reduction mechanism (3D crystal to tubes to planes under increasing contact repulsion) would likely survive.","The same imbalance mechanism should operate in other long-range interacting binary condensates, such as dipolar mixtures, whenever the inter- and intra-species interaction ranges differ, suggesting that planar and tubular supersolid phases are generic rather than specific to Rydberg dressing.","A testable prediction from the paper's Bogoliubov analysis is that the planar and tubular periodicity should equal $2\\pi/k_{\\rm rot}$ of the out-of-phase branch; a measurement of the density correlation function during the dynamic evolution could verify this directly.","The claim that the superfluid fraction drops at the SC-ST and ST-SP transitions could be tested by time-of-flight or Bragg spectroscopy, giving a clear experimental signature of the lower-dimensional phases."],"forward_implications":["Binary Rydberg-dressed BECs become a tunable 3D platform for supersolidity, with crystal structures (FCC, SC, BCC) and lower-dimensional density-modulated states all exhibiting a finite superfluid fraction.","The SC-to-ST-to-SP sequence gives a cold-atom analogue of the 'pasta' phases proposed for neutron-star crusts, with the scattering length playing the role of pressure.","Because the states form dynamically from a uniform condensate with noise, experimental realization only requires preparing a miscible binary mixture and then ramping the scattering length, rather than engineering the target lattice.","At equal blockade radii the model reproduces the single-component FCC-to-BCC transition, so the binary system provides a new route to study BCC supersolidity in 3D.","The multicritical point at $R_c^{(12)}/R_c = 1$ and $a_s/R_c \\approx 1.51 \\times 10^{-3}$ ties together all the crystalline, tubular, planar, and phase-separated phases, offering a single parameter setting to test the whole diagram."],"supporting_citations":[{"why":"Supplies the single-component 3D Rydberg-dressed BEC FCC ground state that the binary model extends.","marker":"[22]"},{"why":"Gives the soft-core single-component FCC result that the equal-radius limit reproduces.","marker":"[23]"},{"why":"Introduces the two-component Rydberg-dressed BEC model with equal blockade radii that this work generalizes to imbalanced radii.","marker":"[25]"},{"why":"Reports prior two-dimensional binary Rydberg-dressed BEC phases that motivate the three-dimensional search.","marker":"[39]"},{"why":"Provides the Gaussian variational trial-density method used to map energy landscapes and locate the SC, ST, SP, FCC, and BCC phases.","marker":"[57–60]"}],"fun_headline_variants":["Rydberg BEC mixtures crystallize then flatten","From 3D crystals to planes: tuning Rydberg BEC mixtures","Binary Rydberg BECs: FCC lattice to sheets","Rydberg BECs: from FCC crystals to tubes and planes","Scattering length drives Rydberg BEC dimensionality reduction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The phase diagram is only as reliable as the assumption that the five trial crystal families plus the periodic-box imaginary-time search exhaust all lower-energy density patterns, so an omitted lattice could change the boundaries.","fun_headline_variants_meta":{"raw":{"variants":["Rydberg BEC mixtures crystallize then flatten","From 3D crystals to planes: tuning Rydberg BEC mixtures","Binary Rydberg BECs: FCC lattice to sheets","Rydberg BECs: from FCC crystals to tubes and planes","Scattering length drives Rydberg BEC dimensionality reduction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000336,"raw_usage":{"total_tokens":1843,"prompt_tokens":909,"completion_tokens":934,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":848}},"tokens_in":525,"tokens_out":934,"duration_ms":8635,"temperature":1.0,"reasoning_tokens":848,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:14:01.627578+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At parameters where the paper predicts an SC ground state (for example $R_c^{(12)}/R_c = 0.8$ and $a_s/R_c = 1.0 \\times 10^{-3}$), run imaginary-time Gross-Pitaevskii evolution in boxes of at least two different sizes and test trial densities for hexagonal and tetragonal lattices; if any alternative structure has lower energy, or the SC state is not reproduced, the phase diagram needs revision.","supporting_citations":[{"cited_title":"Hsueh, Y.-C","cited_arxiv_id":null,"evidence_quote":"Introduces the two-component Rydberg-dressed BEC model with equal blockade radii that this work generalizes to imbalanced radii."},{"cited_title":"Han, X.-F","cited_arxiv_id":null,"evidence_quote":"Reports prior two-dimensional binary Rydberg-dressed BEC phases that motivate the three-dimensional search."}],"review_version":1}