{"id":"0b3c466d-1c06-4fc6-a821-91847fc092b6","arxiv_id":"2607.29292","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In a 1D two-component quantum droplet, a central repulsive barrier drives a sharp atom-number-dependent transition between one-sided (polarized) and symmetric (unpolarized) ground states, with a discontinuous jump and near-zero modes in the excitation spectrum.","lead":"This paper uses computer simulations to study quantum droplets in a 1D gas of two atomic species when a laser-made dimple or bump sits in the middle, mapping how droplet shape, low-energy vibrations, and left-right symmetry change with atom number. It finds a sharp 'polarization transition' for a central bump and very soft, near-zero vibration modes for a dimple - information that could guide experiments on ultracold droplets.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Polarized-branch ground-state claim rests on unverified global minimization; no energy comparison or hysteresis scan supports the N_cr=126 transition.","rationale":"The central claim is that the repulsive defect drives a ground-state polarization transition at N_cr and that this is reflected in a BdG spectral discontinuity. The reader's weakest assumption was that imaginary-time propagation finds the global energy minimum for each N. I agree this is the most load-bearing unverified premise: without an energy comparison between the polarized and symmetric branches, or a hysteresis scan, the reported N_cr=126 could be a numerical artifact of initial-condition-dependent convergence rather than a true ground-state crossing. The spectral discontinuity in Fig. 5(c) is consistent with this but does not discriminate between a genuine transition and a metastable branch switch. This concern directly affects the truth of the central claim. The abstract/conclusion contradiction about the g12-dependence of N_cr is also serious and must be corrected, but it is a presentation error that does not by itself invalidate the numerics; the global-minimum issue is more fundamental. The paper uses standard methods (eGPE imaginary-time propagation, BdG diagonalization) and cites relevant prior work, which lends some credibility, but the missing energy comparison and hysteresis check are explicit gaps that can be closed with a modest computational test. I therefore recommend no change to the CONDITIONAL verdict: the paper should be published only if the authors provide the energy comparison/hysteresis data and fix the contradictory statements in the abstract and conclusion.","tokens_in":14327,"tokens_out":4117,"duration_ms":40656,"concrete_test":"For fixed g12=-0.9, V0=0.3, sigma=1, lambda=0, compute the total energy per particle from Eq. (1) for the fully polarized and the symmetric two-well states at N = 80,100,120,126,140,160,200. Obtain the polarized state by imaginary-time propagation seeded with a single-well localized guess and the unpolarized state by a symmetric seed; verify each is a stationary point (imaginary-time gradient norm < threshold). Plot E_pol(N) - E_unpol(N). The transition is genuine only if this difference changes sign at the same N_cr where the solver switches. Additionally, run an upward sweep N=50->200 and a downward sweep 200->50, using the converged field at each N as the initial condition for the next N; if the transition N differs by more than, say, 5 atoms, the reported N_cr is not a well-defined ground-state threshold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Imaginary-time propagation from 'a random initial state and a superposition of symmetric and asymmetric functions' (Sec. III, opening) is used to obtain the ground state, but no energy comparison between the polarized (N<N_cr) and symmetric/unpolarized (N>=N_cr) candidates is shown. For a first-order-like transition, both branches can be local minima of the eGPE energy functional; the reported N_cr=126 may be the point at which the solver's basin of attraction switches rather than the thermodynamic crossing. The discontinuity in the BdG spectrum (Fig. 5(c)) is consistent with either a genuine ground-state transition or a metastable branch switching. No hysteresis scan (sweeping N up/down with the previous solution as initial guess) is provided, so bistability cannot be excluded. If the polarization transition is a metastability artifact, the central phase diagram (Fig. 2) and the associated spectral softening are not ground-state properties.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a one-dimensional two-component Bose-Bose quantum droplet in a harmonic trap with a localized Gaussian defect, using the extended Gross-Pitaevskii equation with the exact LHY correction and Bogoliubov–de Gennes theory. It reports a polarization transition as a function of atom number in the presence of a repulsive central barrier: for N < N_cr the ground state is fully polarized (the droplet lies in one well), while for N >= N_cr it becomes symmetric. The transition is claimed to appear as a discontinuity and softening of the low-lying collective spectrum, with up to six near-zero modes in the unpolarized regime. The paper also studies the dependence of N_cr on the intercomponent attraction, the role of harmonic confinement, and the real-time dynamics after quenches of the interaction and of the defect strength.","tokens_in":14492,"tokens_out":3425,"duration_ms":37894,"significance":"If the reported transition is a genuine ground-state property, the paper would add a useful example of how a localized defect controls droplet self-binding, polarization, and collective excitations. The theoretical framework is standard and imported from established references (Petrov 2015; Ilg et al. 2018; Mistakidis et al. 2023), and the paper computes outputs rather than fitting parameters. The main claims are falsifiable by direct numerical energy comparison and by experiments with tunable barriers. However, the manuscript currently contains a direct internal contradiction about the direction of the N_cr(g12) trend, and the central ground-state transition is not verified against metastability or basis-convergence artifacts. These issues are load-bearing for the main phase-diagram claim.","major_comments":[{"comment":"The abstract states that \"the critical number of atoms for the transition decreases as the attractive intercomponent interaction increases\" (i.e., stronger attraction -> smaller N_cr), and the Conclusions state that \"the unpolarized state becomes energetically favourable at a smaller number of atoms for larger intercomponent interactions.\" In direct opposition, Sec. III.A (text near Fig. 2) states \"with less attractive g12, N_cr decreases,\" and the discussion says \"with less attractive g12, N_cr decreases... attributed to the attraction that inhibits the transition.\" The body and Fig. 2 therefore claim weaker attraction -> smaller N_cr. Both cannot be correct. This is a central claim of the paper, not a minor typo, because the abstract and conclusions advertise the opposite trend from the main figure.","section":"Abstract; Sec. III.A; Sec. IV (Conclusions)"},{"comment":"The ground state is found by imaginary-time propagation starting from \"a random initial state and a superposition of symmetric and asymmetric functions.\" For a first-order-like transition, both the polarized (N < N_cr) and symmetric (N >= N_cr) solutions can be local minima of the eGPE energy functional. The paper does not compare the energies of these two candidate states, nor does it perform a hysteresis scan (e.g., sweeping N upward and downward while seeding with the previous solution). Thus N_cr = 126 at g12 = -0.9 is not established as a genuine ground-state transition; it may reflect the point where the solver's basin of attraction changes. This concern is load-bearing for the phase diagram and the interpretation of the BdG discontinuity in Fig. 5(c). A direct E_polarized(N) vs. E_symmetric(N) comparison and a hysteresis test should be reported.","section":"Sec. III, opening; Sec. III.A; Fig. 2"},{"comment":"The Bogoliubov–de Gennes matrix is diagonalized in a basis of only 200 harmonic-oscillator states, with no convergence test as a function of basis size. The central spectral claims include up to six modes clustered near zero energy and a spectral discontinuity at the transition. If the basis is insufficient for the broad, flat-top droplets at large N or for the strongly localized fragments in a double well, the near-zero modes could be numerical artifacts. The authors should report a convergence check for at least a few representative points, e.g., N = 150 and N = 200 at V0 = 0.3 and V0 = -0.5, with larger bases (300-500 states) and possibly a spatial-grid check.","section":"Sec. III, opening; Sec. III.B; Figs. 5-6"}],"minor_comments":[{"comment":"The caption labels the V0 = -0.5 case as an \"attractive barrier.\" A negative V0 is a potential well/dimple, not a barrier. Please use consistent terminology.","section":"Fig. 1(a) caption and text"},{"comment":"The inset axis labels are unclear: the horizontal axis is presumably N and the vertical axis is |I|, but no axis titles are shown. Also the jump from |I| = 1 to 0 at N around 126 would be easier to assess with a marker or a vertical dashed line.","section":"Fig. 2 inset"},{"comment":"The text says that for small N the dipole mode \"possesses a small finite energy, which overlaps with zero-energy modes.\" It would help to state explicitly in the figure or text that the green star at small N is not exactly zero, since the legend otherwise appears to show a zero-energy branch.","section":"Sec. III.B, Fig. 5(c) discussion"},{"comment":"The text says the unpolarized state has six modes clustered around zero, but the list in Sec. III.B mentions two U(1) modes and four low-lying modes (i-iv). That totals six, yet the subsequent sentence says \"consists of six modes clustered around zero\" and then \"in addition, two nearly degenerate breathing modes emerge.\" Please clarify whether the six include the breathing pair or are only the symmetry-related modes.","section":"Sec. III.B, mode counting"}],"recommendation":"major_revision","confidential_remarks":"The contradiction between the abstract/conclusion and the body regarding N_cr(g12) is striking and must be resolved; it is probably a wording error but it directly affects the paper's headline claim. The more substantive issue is that the polarization transition is not validated against metastability: an energy comparison and hysteresis scan are necessary before N_cr can be called a ground-state transition. The BdG basis-size check is also needed. These are fixable within the manuscript's scope, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bob,\n\nQuick take: the paper is a competent, standard-tool numerical study that is undermined by a self-contradiction in its central quantitative claim. The body and Fig. 2 show the critical atom number for the polarization transition decreases as g12 becomes less attractive; the abstract and conclusion state the opposite. That is not a cosmetic mismatch—it is the headline result.\n\nWhat is genuinely new: for a 1D two-component droplet with a central Gaussian barrier, they map out N_cr(g12), identify up to six near-zero BdG modes in the unpolarized phase, show the dipole softens as V0 or λ changes, and give quench dynamics that align with recent fragmentation experiments. The eGPE+BdG machinery is standard and imported from solid references; the numerics look internally consistent, and the six-mode structure is physically plausible for two weakly coupled droplets.\n\nThe soft spots, in order:\n1. The abstract/conclusion vs. body contradiction. This has to be fixed before anything else.\n2. No energy comparison between the polarized and symmetric candidates, and no hysteresis scan. The stress-test note is right: the transition at N_cr=126 could be a basin-of-attraction switch in imaginary-time propagation rather than a true ground-state crossing. I don't think that's likely, but the paper needs to show it.\n3. The 200-state BdG basis has no convergence test. Minor but easy to add.\n4. No code or data deposition. That makes independent verification harder.\n5. The claim that the LHY energy is valid 'for any value' is an overstatement; the 1D correction has a known range of validity (p > -0.618). It doesn't affect the results presented, but it should be qualified.\n\nWho is this for? People working on droplets in optical potentials, especially 39K experiments. It is not a breakthrough, but it is a useful mapping of a regime that has not been systematically explored.\n\nMy recommendation: send it to peer review. The internal contradiction is severe but fixable, and the missing energy comparison is exactly what a referee should request. The core physics is plausible enough to warrant that effort.","headline":"A conventional but useful droplet-defect study whose headlined result is stated backwards in the abstract and conclusion, and whose ground-state transition lacks an energy comparison.","tokens_in":15094,"tokens_out":3523,"would_cite":false,"duration_ms":36364,"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 central repulsive barrier drives a one-dimensional binary quantum droplet through a polarization transition at a critical atom number, and the excitation spectrum exposes the switch as a discontinuity and six near-zero modes.","keywords":["quantum droplets","two-component Bose mixtures","Lee-Huang-Yang correction","polarization transition","Bogoliubov-de Gennes spectrum","double-well potential","population imbalance","quench dynamics"],"falsifier":"Directly compare the energies of the polarized and symmetric stationary states for N around 126 at g12 = -0.9, V0 = 0.3, sigma = 1, lambda = 0; if the symmetric state has lower energy below 126, the claimed transition is a solver artifact. A second check is to scan N upward and downward for hysteresis: a genuine ground-state transition should reproduce Ncr from both directions, while a metastable branch would jump at different N depending on the direction.","tokens_in":104,"feed_emoji":"💧","tokens_out":6813,"duration_ms":133137,"temperature":0.7,"pith_summary":"This paper studies a one-dimensional two-component ultradilute quantum droplet — a self-bound state stabilized by the balance between attractive mean-field interactions and repulsive Lee-Huang-Yang quantum fluctuations — placed in a harmonic trap with a central Gaussian defect. Its central claim is that a repulsive barrier makes the ground state polarize spontaneously: for fewer than a critical number of atoms (Ncr = 126 for g12 = -0.9, V0 = 0.3, sigma = 1, lambda = 0) the droplet sits entirely in one well of the double-well potential, while at and above Ncr it becomes symmetric and occupies both wells. That polarization transition appears in the Bogoliubov excitation spectrum as a discontinuity and a strong softening of low-lying modes, with up to six modes clustering near zero energy in the unpolarized phase. The authors also examine how Ncr shifts with intercomponent attraction and how sudden or adiabatic quenches of the interaction and defect produce breathing, localized filaments, and fragmentation into multiple droplets. If the picture is right, the droplet in a double well is a clean, controllable example of an atom-number-driven symmetry-breaking transition with a measurable spectral fingerprint.","feed_headline":"Droplet in a double well flips from one side to both at 126 atoms","feed_subtitle":"A repulsive barrier forces a 1D quantum droplet to choose one well, then spreads it symmetrically — a switch visible in its excitation spect","key_machinery":"The key machinery is the extended Gross-Pitaevskii equation (eGPE) with the exact one-dimensional Lee-Huang-Yang correction, solved by imaginary-time propagation, and linearized through Bogoliubov-de Gennes (BdG) equations on a basis of 200 harmonic-oscillator states. The central diagnostic is the population imbalance I = (NL - NR)/(NL + NR), which acts as the order parameter: it jumps from 1 (polarized) to 0 (unpolarized) at Ncr. The spectral signatures — softening of the dipole mode, a discontinuity in mode frequencies, and the clustering of six modes near zero — connect the ground-state transition to observable collective dynamics.","core_discovery":"On the paper's own terms, the discovery is that the competition between barrier-induced localization and interaction-driven delocalization produces a sharp ground-state transition in a one-dimensional two-component quantum droplet. Below a critical atom number Ncr (126 for g12 = -0.9, V0 = 0.3, sigma = 1, lambda = 0), the LHY-stabilized droplet is fully polarized in one well of the double-well potential, with population imbalance |I| = 1. At Ncr it discontinuously rearranges into an unpolarized state with equal occupation of both wells (I = 0). The transition is not only static: the Bogoliubov-de Gennes quasiparticle spectrum shows a discontinuity at Ncr, and the unpolarized phase supports u","pith_inferences":["Editorial: The stated direction of the Ncr(g12) trend is internally inconsistent: the abstract and conclusion say stronger attraction lowers Ncr, while Sec. III and the Fig. 2 caption say weaker attraction lowers Ncr. The figure as drawn appears to follow the body's version; the discrepancy is unresolved in the text.","Editorial: Because no energy comparison between the polarized and symmetric stationary states is shown, and no hysteresis scan is reported, the transition boundary Ncr(g12) should be read as a property of the solver's converged state rather than a proven ground-state phase boundary.","Editorial: The six near-zero modes depend on a BdG basis of 200 oscillator states with no reported convergence test; checking against larger bases or against real-time evolution spectra would confirm they are physical rather than numerical.","Editorial: A natural extension is to map Ncr as a function of barrier height V0 and width sigma, and to compare the predicted spectral discontinuity with sum-rule or many-body estimates for the 1D Bose mixture."],"forward_implications":["If the transition is real, a 1D two-component droplet in a double-well trap is a clean model system for an atom-number-driven symmetry-breaking transition, with Ncr as a tunable control parameter.","The six near-zero modes predicted in the unpolarized phase imply slow inter-well dynamics, including Josephson-type oscillations and coupled dipole excitations, that could be seen in real-time imaging of droplet mixtures.","The softening of the dipole mode with increasing N, and as the defect is tuned from attractive to zero, provides a measurable precursor of the transition and of the restoration of translational invariance.","The quench protocols — breathing for weak quenches, filaments and fragmentation for strong quenches — offer experimental routes to controllably create localized fragments or multi-droplet arrays.","The discontinuity in the quasiparticle spectrum at Ncr gives an observable, frequency-domain marker of the ground-state switch, complementary to density-imaging measurements of population imbalance."],"fun_headline_variants":["Quantum droplet flips wells at 126 atoms","Barrier splits 1D droplet into both wells above threshold","Critical atom number sets droplet polarization switch","Two-well droplet transitions from one side to both at N=126","Attractive interactions lower the droplet's critical atom number"],"cache_read_input_tokens":16256,"weakest_assumption_plain":"The load-bearing premise is that imaginary-time propagation from random and mixed-symmetry initial states converges to the global energy minimum at every atom number, so the polarized branch below Ncr is the true ground state rather than a metastable artifact of the solver.","fun_headline_variants_meta":{"raw":{"variants":["Quantum droplet flips wells at 126 atoms","Barrier splits 1D droplet into both wells above threshold","Critical atom number sets droplet polarization switch","Two-well droplet transitions from one side to both at N=126","Attractive interactions lower the droplet's critical atom number"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000238,"raw_usage":{"total_tokens":1322,"prompt_tokens":696,"completion_tokens":626,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":440,"completion_tokens_details":{"reasoning_tokens":549}},"tokens_in":440,"tokens_out":626,"duration_ms":7127,"temperature":1.0,"reasoning_tokens":549,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T09:51:22.480687+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly compare the energies of the polarized and symmetric stationary states for N around 126 at g12 = -0.9, V0 = 0.3, sigma = 1, lambda = 0; if the symmetric state has lower energy below 126, the claimed transition is a solver artifact. A second check is to scan N upward and downward for hysteresis: a genuine ground-state transition should reproduce Ncr from both directions, while a metastable branch would jump at different N depending on the direction.","supporting_citations":[],"review_version":1}