Pith. sign in

REVIEW 2 major objections 5 minor 22 references

Vacancy Effect on the Ground-State Energy of a Bose Gas Trapped by 1D Imperfect Artificial Crystal

T0 review · 2 major / 5 minor · reviewed 2026-07-13 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2607.08992 v1 pith:XDGSHQG7 submitted 2026-07-09 cond-mat.quant-gas cond-mat.mtrl-sci

E.I. Guerrero-Cruz , O.A. Rodríguez-López , M.A. Solís This is my paper · ORCID
classification cond-mat.quant-gascond-mat.mtrl-sci PACS 03.75.Hh05.30.Jp67.85.-d
keywords BosegasGross-PitaevskiiequationDiraccombvacanciesground-stateenergyone-dimensionalopticallatticemean-fieldtheorychemicalpotential
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

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.

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 (2)
  1. [Sec. IV.A, Fig. 6, Table IV] 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).
  2. [Secs. II–V] 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.
minor comments (5)
  1. [throughout] 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.
  2. [Fig. 6] 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.
  3. [Sec. IV] 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.
  4. [Sec. IV] The sampling protocol (“error percentage ho 2 %”) is mentioned only briefly; the number of independent vacancy realizations actually used for each data point should be reported.
  5. [References] 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.

Circularity Check

0 steps flagged · score 1.0 of 10

No load-bearing circularity: GSE values are obtained by direct numerical minimization of the GPE functional; exponential fits and 1/N extrapolations are post-hoc descriptions of those data, not inputs that force the claims.

full rationale

The derivation chain is self-contained numerical work. The stationary GPE (Eqs. 7–8, 19) is minimized by the standard Gradient-Flow with Discrete Normalization / imaginary-time method (Sec. III, citing Bao and external references), subject only to the external Dirac-comb potential and the interaction term. Ground-state energies, chemical potentials and densities are therefore outputs of that minimization (Tables I–V, Figs. 1–8). The free-gas benchmark is the external Lieb–Liniger formula (Sec. V); the perfect-crystal and free limits supply independent anchors. Exponential curves versus vacancy percentage and the regression versus 1/N_deltas (Fig. 6) are fitted after the fact to the already-computed energies; they do not enter the GPE solver or redefine the reported numbers. Self-citations to the group’s earlier ideal-gas vacancy theses appear only as background for the non-interacting gap and BEC motivation; they are not used to justify any interacting-gas result. No uniqueness theorem, ansatz, or definitional identity closes a logical loop. The unphysical overshoot of the one-vacancy infinite-system extrapolations for g>0 is a correctness/finite-size issue, not circularity.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The central claims rest on the standard mean-field GPE, the Dirac-comb model of an artificial crystal, the imaginary-time convergence theorem, and a handful of hand-chosen numerical parameters (P0, g grid, system sizes, 2 percent sampling cutoff). No new physical entities are postulated; free parameters appear only in post-hoc exponential fits and in the conventional choice of lattice strength. All load-bearing modeling choices are domain-standard rather than ad-hoc inventions.

free parameters (3)
  • delta strength P0 = 10
    Fixed by hand at P0 = 10 for all production runs; controls the depth of the artificial crystal and therefore the absolute energy scale of every reported gap and exponential.
  • exponential fit parameters (amplitude, decay constant, offset) = e.g. 0.1619 + 6.7307 exp(-VP/4.2414) for g=0
    Three-parameter exponential forms are fitted to GSE versus vacancy percentage for each g; the fitted numbers are used to summarize the central claim of exponential decrease.
  • sampling cutoff (2 percent energy error) = 2%
    Number of random vacancy realizations is increased until the relative error on GSE falls below 2 percent; this ad-hoc threshold determines which data points enter the reported averages and fits.
assumptions (4)
  • domain assumption Mean-field Gross-Pitaevskii equation accurately describes the ground state of a weakly interacting 1D Bose gas for the scanned range of g.
    Invoked from Sec. II onward; the entire numerical pipeline solves the GPE rather than the many-body Schrödinger equation.
  • standard math Imaginary-time gradient flow with discrete normalization converges to the true GPE ground state under periodic boundaries.
    Standard justification of the GFDN method (Sec. III); used without additional proof.
  • domain assumption A finite Dirac comb of strength P0 with randomly deleted deltas is a faithful model of an imperfect 1D artificial crystal.
    Defined in Sec. II and IV; vacancies are realized simply by zeroing selected delta coefficients.
  • ad hoc to paper Finite-size GSE data versus 1/N_deltas can be extrapolated by a simple exponential or linear fit to obtain the infinite-system limit.
    Performed in Sec. V and Fig. 6 for the one-vacancy case; the functional form of the extrapolation is chosen by the authors.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Vacancy Effect on the Ground-State Energy of a Bose Gas Trapped by 1D Imperfect Artificial Crystal." pith.science (2026). https://pith.science/paper/XDGSHQG7

@misc{pith2026260708992,
  author       = {Pith},
  title        = {Pith review of: Vacancy Effect on the Ground-State Energy of a Bose Gas Trapped by 1D Imperfect Artificial Crystal},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XDGSHQG7}},
  note         = {Machine review of arXiv:2607.08992}
}
abstract

For a weakly interacting Bose gas trapped by an imperfect one-dimensional artificial crystal, we study the effect of its punctual defects, i.e. vacancies, on the ground state properties of the system. In the framework of the mean field approximation, we numerically solve the corresponding Gross-Pitaevskii equation using the ``Gradient Flow with Discrete Normalization'' method, also known as the imaginary time method. The crystal is artificially produced by applying an external Dirac comb potential to the Bose gas where vacancies are created by randomly removing a predetermined number of deltas. We observe that as the number of randomly removed deltas increases, the ground state energy decreases exponentially from its value for the perfect crystal case until it reaches its value when the Bose gas is free. The ground state energy is reported for different magnitudes of the interaction between bosons and several system sizes which we extrapolate to infinity for the crystal with only one vacancy. Also, we observe the presence of an energy gap between the ground state energies of the perfect system and that with a vacancy, which is more noticeable for values of the particle interaction magnitude $ g \leq 0.1$, when the delta strength $P_0 = 10$. In addition, we report the boson distributions within the crystal, %inside a box with periodic boundary conditions, i.e. the probability density functions which show localization features around vacancies which disappear as $g$ increases. From the ground state energy, the chemical potential is obtained immediately.

Figures

Figures reproduced from arXiv: 2607.08992 by the authors.

Figure 2
Figure 2. Ground state energy, in units of ℏ 2 /2ma2 , for a Bose gas within of 50, 100, 500 and 1000 total initial deltas where we randomly removed 10%, 20%, 30%, 40% and 50%, with P0 = 10 and g = 1.We did an exponential fit on the 1000 delta values, this is shown in the curve. Removed Initial number of deltas % 50 100 500 1000 0 8.1543 8.1543 8.1543 8.1543 10 6.4482 6.4598 6.4577 6.4637 20 5.1183 5.1084 5.0901 5.0623 30 4.0… view at source ↗
Figure 1
Figure 1. Ground state energy, in units of ℏ 2 /2ma2 , for an ideal Bose gas (g = 0) within of 50, 100, 500 and 1000 initial deltas where we randomly removed 10%, 20%, 30%, 40% y 50%, P0 = 10. The full curve is an exponential fit for the 1000 delta values. Removed Initial number of deltas % 50 100 500 1000 0 6.8942 6.8942 6.8942 6.8942 10 1.8338 1.2496 1.0210 0.7638 20 1.0511 0.5320 0.4836 0.3943 30 0.6741 0.4841 0.2827 0.259… view at source ↗
Figure 3
Figure 3. Ground state energy for systems with 50, 100, 500 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Probability density function of the ground state of [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 6
Figure 6. Figure 6: Extrapolation of the ground state energy to its [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 8
Figure 8. Figure 8: Enlargement of Fig [PITH_FULL_IMAGE:figures/full_fig_p007_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

22 extracted references

  1. [1]

    Bose–Einstein condensation for general dispersion rela- tions

    V. C. Aguilera-Navarro, M. de Llano and M. A. Solís, “Bose–Einstein condensation for general dispersion rela- tions”, Eur. J. Phys.201, 77 (1999)

    1999

  2. [2]

    Inves- tigation of Bose Condensation in Ideal Bose Gas Trapped under Generic Power Law Potential indDimension

    M. M. Faruk, M. S. Hossain and M. M. Rahman, “Inves- tigation of Bose Condensation in Ideal Bose Gas Trapped under Generic Power Law Potential indDimension", Commun. Theor. Phys.65(2016)

    2016

  3. [3]

    Bose- Einstein condensation in an external potential

    V. Bagnato, D. E. Pritchard and D. Kleppner, “Bose- Einstein condensation in an external potential”, Phys. Rev. A35, 10 (1986)

    1986

  4. [4]

    Bose gas with generalized dispersion relation plus an en- ergy gap

    J. G. Martínez-Herrera, J. García-Nila and M. A. Solís, “Bose gas with generalized dispersion relation plus an en- ergy gap”, Phys. Scr.94, 075002 (2019)

    2019

  5. [5]

    Bose Einstein condensation in crystals with vacancies

    J. García, “Bose Einstein condensation in crystals with vacancies", Master Thesis (In Spanish) (Universidad Na- cional Autónoma de México, Ciudad de México, 2019) Online PDF

    2019

  6. [6]

    Effect of vacancies on the critical BEC temperature of an ideal boson gas within imperfect crys- tals

    E. I. Guerrero, “Effect of vacancies on the critical BEC temperature of an ideal boson gas within imperfect crys- tals", Bachelor Thesis (In Spanish) (Universidad Na- cional Autónoma de México, Ciudad de México, 2020) Online PDF

    2020

  7. [7]

    Bose-Einstein Condensation in Multilayers

    P. Salas, M. Fortes, M. Llano, F. J. Sevilla and M. A. Solís, “Bose-Einstein Condensation in Multilayers" J. Low Temp. Phys.159, 540–548 (2010)

    2010

  8. [8]

    Pitaevskii and S

    L. Pitaevskii and S. Stringari,Bose-Einstein Conden- sation and Superfluidity(Oxford University Press, New York, 2016) 1st ed

    2016

Show all 22 references
  1. [9]

    Eighty years of superfluidity

    W. P. Halperin, “Eighty years of superfluidity", Nature 553, 413-414 (2018)

    2018

  2. [10]

    Quantum theory of crystal defects

    A. Andreev and I. Lifshitz, “Quantum theory of crystal defects", Sov. Phys. JETP29, 1107 (1969)

    1969

  3. [11]

    Speculations on Bose-Einstein conden- sation and quantum crystals

    G. V. Chester, “Speculations on Bose-Einstein conden- sation and quantum crystals" Phys. Rev. A2, 256-258 (1970)

    1970

  4. [12]

    Exact Analysis of an Inter- acting Bose Gas. The General Solution and the Ground State

    E. H. Lieb and W. Liniger,“Exact Analysis of an Inter- acting Bose Gas. The General Solution and the Ground State”, Phys. Rev.130, 1605 (1963)

    1963

  5. [13]

    Mobileimpuritiesinteract- ing with a few onedimensional lattice bosons

    V.R.YordanovandF.Isaule,“Mobileimpuritiesinteract- ing with a few onedimensional lattice bosons”, J. Phys. B: At. Mol. Opt. Phys.56045301 (2003)

    2003

  6. [14]

    Quantum States of ultra- cold bosons in optical lattices interacting via long-range interaction

    R. Panda and B. Chatterjee,“Quantum States of ultra- cold bosons in optical lattices interacting via long-range interaction ”, AIP Conf. Proc.3149, 130011 (2024)

    2024

  7. [15]

    Structural superfluid–Mott-insulator transition for a Bosegas in multirods

    O. A. Rodríguez-López, M. A. Solís and J. Boronat, “Structural superfluid–Mott-insulator transition for a Bosegas in multirods", Phys. Rev. A103, 013311 (2021)

    2021

  8. [16]

    One-dimensional Bose gas in opti- cal lattices of arbitrary strength,

    G. E. Astrakharchik, K. V. Krutitsky, M. Lewenstein, and F. Mazzanti, “One-dimensional Bose gas in opti- cal lattices of arbitrary strength," Phys. Rev. A93, 021605(R) (2016)

    2016

  9. [17]

    B. T. Seaman, L. D. Carr, and M. J. Holland, Phys. Rev. A 71, 033609 (2005)

    2005

  10. [18]

    Bose- Einstein condensates and the numerical solution of the Gross-Pitaevskii equation

    S. Succi, F. Toschi, M. P. Tosi and P. Vignolo, “Bose- Einstein condensates and the numerical solution of the Gross-Pitaevskii equation" in Computing in Science & Engineering,7, 6 (2005)

    2005

  11. [19]

    Quantum Mechanics of Electrons in Crystal Lattices

    R. de L. Kronig, W. G. Penney, “Quantum Mechanics of Electrons in Crystal Lattices”, Proy. Roy. Soc.130, 499 (1930)

    1930

  12. [20]

    Ground states and dynamics of multi- component Bose-Einstein condensates

    W. Bao, “Ground states and dynamics of multi- component Bose-Einstein condensates”, Multiscale Mod- eling and Simulation: a SIAM Interdisciplinary Journal, 2, 210 (2004)

    2004

  13. [21]

    Fortranprograms for the time-dependent Gross–Pitaevskii equation in a fully anisotropic trap

    S.K.AdhikariandP.Muruganandam, “Fortranprograms for the time-dependent Gross–Pitaevskii equation in a fully anisotropic trap", Computer Physics Communica- tions,180, 10 (2009)

    2009

  14. [22]

    Mathematical theory and numerical methods for Bose-Einstein Condensation

    W. Bao and Y. Cai, “Mathematical theory and numerical methods for Bose-Einstein Condensation”, Kinetic and Related Models6, 1 (2013)

    2013

Pith tools

Reviewed July 13, 2026 · model on record in the stance chip above.