{"id":"a6b149d2-6fd9-49f2-a869-76af0cdc2b76","arxiv_id":"2607.08992","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Ground-state energy of a weakly interacting 1D Bose gas in a Dirac-comb lattice falls exponentially with vacancy fraction, with a noticeable perfect-vs-one-vacancy gap at weak g and localization around defects that vanishes as g rises.","lead":"Numerical mean-field calculations show that vacancies in a 1D Dirac-comb lattice lower the ground-state energy of a weakly interacting Bose gas exponentially toward the free-gas value, with localization around defects that fades as interactions grow. The results quantify how point defects reshape energies and densities in artificial crystals relevant to optical-lattice experiments.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"One-vacancy infinite-system extrapolations for g>0 are unphysical: they overshoot the perfect-crystal E/N instead of approaching it from below as extensivity requires.","rationale":"The reader correctly flags mean-field reliability in 1D (already visible in the free-gas 17 % discrepancy with Lieb-Liniger) as a broad caveat that affects every quantitative claim. The more immediate, load-bearing flaw, however, is internal to the GPE calculations themselves: the thermodynamic-limit extrapolation for a single vacancy is performed and interpreted incorrectly once g>0. Because that extrapolation is explicitly part of the strongest claim and of the abstract, the finite-size exponential-with-vacancy-percentage results and the localization plots remain usable within mean-field, but the infinite-system numbers and the persistence of a gap must be discounted. This keeps the verdict CONDITIONAL; it does not push the work to REJECT. The two concerns are complementary rather than identical, hence partial agreement.","tokens_in":13331,"tokens_out":572,"duration_ms":80146,"concrete_test":"Recompute the GPE ground-state energy per particle for one vacancy at g=1, P0=10 and fixed density, using system sizes from ~10^3 up to at least 10^4–10^5 periods (with demonstrated grid convergence, e.g. ≥20 points per a). Plot E versus 1/N; if the data approach the perfect value 8.154 from below with a clean 1/N correction, the paper’s extrapolated entries (Table IV, Fig. 6) and the associated infinite-system gap claims are invalidated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"For any g>0 the ground-state energy per particle with a single vacancy must approach the perfect-crystal value in the thermodynamic limit at fixed density. A vacancy produces only a localized density enhancement whose integrated excess particle number remains O(1); the associated energy correction is therefore O(1) and vanishes as O(1/N) in E/N. The paper’s exponential fits versus 1/N_deltas (Fig. 6, Table IV) instead produce asymptotes lying above the perfect values (8.742 vs 8.154 for g=1; likewise for g=0.1, 0.5, 2). This shows that either the simulated sizes remain pre-asymptotic or the fitting ansatz is inappropriate, so the reported infinite one-vacancy energies and any residual gap are unreliable for interacting gases. (Only the g=0 case is consistent, because every particle occupies the single-particle bound state.)","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.5","summary":"The manuscript studies the zero-temperature ground-state energy (GSE), chemical potential, and probability density of a weakly interacting 1D Bose gas subject to a Dirac-comb external potential that models a finite artificial crystal. Vacancies are introduced by randomly deleting a prescribed fraction of the delta functions. Within the mean-field Gross–Pitaevskii equation the authors employ the gradient-flow-with-discrete-normalization (imaginary-time) method with a Fourier pseudo-spectral spatial discretization. They report that the GSE falls exponentially with vacancy percentage from the perfect-crystal value toward the free-gas value, that a single vacancy produces a noticeable energy lowering only for g ≤ 0.1 (at fixed delta strength P0 = 10), that the density localizes about vacancies for small g, and that finite-size data for one vacancy can be extrapolated to an infinite-system limit.","tokens_in":13582,"tokens_out":1134,"duration_ms":37278,"significance":"If the numerical trends survive scrutiny they supply concrete mean-field benchmarks for how point defects modify the GSE and density profile of a 1D Bose gas in an optical-lattice-like potential—information of interest for experiments with engineered defects and for the broader discussion of vacancy-assisted supersolidity. The free-gas comparison with the exact Lieb–Liniger solution, the systematic tables for several system sizes and interaction strengths, and the explicit localization plots are useful reference data. The work is a natural extension of the authors’ earlier ideal-gas studies of imperfect crystals.","major_comments":[{"comment":"Sec. IV.A, Fig. 6 and Table IV: for every g > 0 the extrapolated one-vacancy GSE lies above the perfect-crystal value (e.g. 8.742 versus 8.154 for g = 1). A single vacancy produces only a localized density excess of O(1) particles; the associated energy correction is therefore O(1) and must vanish as O(1/N) in the energy per particle. Consequently E/N must approach the perfect-crystal limit from below. The reported asymptotes, the statement that “the GSE tends to grow” with system size, and the infinite-system entries in Tables IV–V are unphysical. Either the exponential fitting ansatz is inappropriate or the simulated sizes remain pre-asymptotic. This directly undermines the abstract claim of an infinite-system extrapolation and any residual-gap discussion for interacting gases (the g = 0 case is consistent because every particle occupies the single-particle bound state).","section":"Sec. IV.A, Fig. 6, Table IV"},{"comment":"Secs. II–V and the free-gas benchmark of Sec. V: the Gross–Pitaevskii mean-field description is known to be only approximate in one dimension. The authors themselves record discrepancies of up to 17 % with the exact Lieb–Liniger free-gas energies. All quantitative statements—exponential decay constants, the size of the one-vacancy gap for g > 0, and the infinite-system extrapolations—must therefore be explicitly caveated as mean-field results. Without additional checks (e.g., DMRG or exact diagonalization on small systems) it remains unclear which features survive beyond mean field.","section":"Secs. II–V"}],"minor_comments":[{"comment":"Numerous typographical and formatting errors appear throughout (title page “and and M.A. Solís”, section headings “MA THEMA TICAL”, “V ACANCIES”, “tree cases”, “id est”, inconsistent spacing in equations). A careful proof-reading pass is required.","section":"throughout"},{"comment":"Fig. 6 caption quotes the fit 8.015 + 0.708 exp(−1105.995 ND) yet states the infinite-system limit as 8.724; the two numbers are inconsistent and should be reconciled.","section":"Fig. 6"},{"comment":"The precise relation between the dimensionless delta strength P0 = 10 and the physical parameter v0 (or the dimensionless V'ext) is never stated; a short clarifying sentence would help reproducibility.","section":"Sec. IV"},{"comment":"The sampling protocol (“error percentage \rho 2 %”) is mentioned only briefly; the number of independent vacancy realizations actually used for each data point should be reported.","section":"Sec. IV"},{"comment":"References [5] and [6] point to theses that are not readily accessible; either deposit them in a public repository or replace them with peer-reviewed citations where possible.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The unphysical thermodynamic-limit extrapolations for g > 0 are the single most serious technical flaw; once they are corrected or removed the remaining finite-size data and exponential trends look publishable after ordinary revision. The 17 % mean-field discrepancy with Lieb–Liniger is already acknowledged by the authors and need not be fatal if properly caveated. Scope is appropriate for a specialized quant-gas journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The core result is usable reference data: GPE ground-state energies for a weakly interacting 1D Bose gas in a Dirac-comb lattice drop exponentially with vacancy percentage, densities pile up at the missing deltas until g grows, and a perfect-vs-one-vacancy gap is visible mainly for g ≤ 0.1 at P0=10. They solve the stationary equation with standard imaginary-time GFDN plus DFT pseudo-spectral methods, give tables for several system sizes and g, recover the free-gas limit, and track Lieb-Liniger within the stated 17 percent. That is honest numerical progress inside the subfield and extends the perfect-comb work of Seaman et al. plus the group’s earlier ideal-gas vacancy papers.\n\nWhat they do well is the systematic scan and the density plots; the exponential fits to vacancy percentage are clean descriptive summaries, not circular. Chemical potentials follow immediately from the energies.\n\nThe soft spot that matters is the one-vacancy thermodynamic extrapolation. For any g>0 a single vacancy produces only an O(1) energy correction, so E/N must approach the perfect-crystal value from below as 1/N→0. Their fits (Fig. 6, Table IV) instead land above the perfect numbers (8.742 vs 8.154 for g=1, similarly for other g>0). Only the g=0 case is consistent. That means either the sizes are still pre-asymptotic or the exponential ansatz in 1/N is wrong; the reported infinite one-vacancy energies and residual gaps for interacting gases are therefore unreliable. Mean-field itself is already only approximate in 1D, which the free-gas discrepancy already flags. No code or raw data is supplied, but that is secondary.\n\nThis is for people who run optical-lattice GPE simulations and want concrete numbers for imperfect combs. It is not a conceptual breakthrough, yet the finite-size data and localization profiles are worth having once the extrapolation claim is fixed or dropped. I would send it to referees; the numerics are careful enough to deserve that scrutiny and a revision cycle on the thermodynamic limit.","headline":"Solid GPE numerics on random vacancies in a 1D Dirac comb, with clear exponential energy drop and localization, but the one-vacancy infinite extrapolations for g>0 are unphysical and overshoot the perfect crystal.","tokens_in":14192,"tokens_out":552,"would_cite":false,"duration_ms":24286,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.75.Hh","05.30.Jp","67.85.-d"],"model":"grok-4.5","headline":"Vacancies in a 1D Dirac-comb lattice lower the interacting Bose ground-state energy exponentially toward the free-gas value, open a clear gap at weak coupling, and pin density until interactions wash the pinning out.","keywords":["Bose gas","Gross-Pitaevskii equation","Dirac comb","vacancies","ground-state energy","one-dimensional optical lattice","mean-field theory","chemical potential"],"falsifier":"An exact or high-precision quantum Monte Carlo calculation of the ground-state energy for the same Dirac-comb-plus-vacancies problem at g ≤ 0.1 that fails to recover either the reported exponential decay with vacancy fraction or the finite gap between perfect and one-vacancy systems would falsify the central claim.","tokens_in":14217,"feed_emoji":"❄️","tokens_out":758,"duration_ms":6566,"temperature":0.7,"pith_summary":"This paper asks what a single missing scatterer, or a random fraction of missing scatterers, does to the ground-state energy and density of a weakly interacting one-dimensional Bose gas held by an artificial crystal. The crystal is built from a Dirac comb of delta-function barriers; vacancies are made simply by deleting some of those deltas. Solving the Gross-Pitaevskii equation by imaginary-time evolution shows that the ground-state energy falls exponentially from its perfect-crystal value to the free-gas value as the vacancy fraction grows. A single vacancy already opens a noticeable energy gap relative to the perfect lattice when the interaction strength is small (g ≤ 0.1 for barrier height P0 = 10). The same vacancy acts as an attractive centre that piles particles around the missing site; that localisation fades once interactions become stronger. The authors also extrapolate the one-vacancy energy to infinite system size and extract the chemical potential directly from the energy. The practical interest is that these defects restore the possibility of finite-temperature condensation even in one dimension, and they give a concrete, tunable handle on how disorder reshapes the lowest-energy state of a lattice Bose gas.","feed_headline":"Vacancies lower 1D Bose ground energy exponentially","feed_subtitle":"A single missing barrier opens a gap at weak coupling and pins density until interactions erase the pin.","key_machinery":"The stationary Gross-Pitaevskii equation for a Dirac-comb external potential, solved by the Gradient Flow with Discrete Normalization (imaginary-time) method under periodic boundary conditions; vacancies are introduced by randomly deleting a prescribed fraction of the delta barriers.","core_discovery":"As the fraction of randomly removed deltas increases, the ground-state energy of the weakly interacting Bose gas decreases exponentially from its perfect Dirac-comb value to the free-gas value. A single vacancy already produces a measurable energy gap that is largest for g ≤ 0.1 (at P0 = 10), and the probability density localises around vacancy sites until interactions wash the localisation away.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Vacancies cut 1D Bose ground energy exponentially to free-gas limit","Single vacancy opens weak-g energy gap in Dirac-comb Bose gas","Bose density pins at vacancies until interactions wash it out","More missing deltas lower ground energy from crystal to free value","One vacancy gaps Bose energy most for g≤0.1 at P0=10"],"cache_read_input_tokens":7552,"weakest_assumption_plain":"The mean-field Gross-Pitaevskii description remains quantitatively reliable for the one-dimensional interacting gas across the scanned range of interaction strengths, even though free-gas benchmarks already differ from the exact Lieb-Liniger energies by up to 17 percent.","fun_headline_variants_meta":{"raw":{"variants":["Vacancies cut 1D Bose ground energy exponentially to free-gas limit","Single vacancy opens weak-g energy gap in Dirac-comb Bose gas","Bose density pins at vacancies until interactions wash it out","More missing deltas lower ground energy from crystal to free value","One vacancy gaps Bose energy most for g≤0.1 at P0=10"]},"model":"grok-4.5","effort":"low","cost_usd":0.00345,"raw_usage":{"total_tokens":1193,"prompt_tokens":831,"num_sources_used":0,"completion_tokens":76,"cost_in_usd_ticks":34500000,"prompt_tokens_details":{"text_tokens":831,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":286,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":831,"tokens_out":76,"duration_ms":3993,"temperature":1.0,"reasoning_tokens":286,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T01:09:20.721833+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"An exact or high-precision quantum Monte Carlo calculation of the ground-state energy for the same Dirac-comb-plus-vacancies problem at g ≤ 0.1 that fails to recover either the reported exponential decay with vacancy fraction or the finite gap between perfect and one-vacancy systems would falsify the central claim.","supporting_citations":[],"review_version":1}