REVIEW 3 major objections 6 minor 28 references
Universal scaling law for quantum droplet formation
T0 review · 3 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read Quantum droplet size in a many-droplet condensate follows the universal law $r_d \propto v^{-1/3}$, where $v$ is the rate at which parameters change during droplet formation.
desk verdict A genuinely new and testable scaling law for quantum droplet size, backed by decent numerics, but the derivation as written does not justify the 1/3 exponent; the simulations carry the paper. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the nonlinear Bogoliubov mode equation in the small-$k$ limit, $m\delta\ddot\psi - g n k^2 v t \,\delta\psi + O(\delta\psi^2)=0$, obtained from the Lee-Huang-Yang-corrected Gross-Pitaevskii equation. Rescaling time by $\tilde t = v^{1/3}t$ and rescaling $\delta\psi$ removes all $v$ dependence from the mode dynamics, so each mode becomes nonlinear at a fixed rescaled time; this fixes the freeze-out time $t_f \propto v^{-1/3}$. The length scale is then read off the dispersion relation $\omega^2 = (k^2/2m)(\hbar^2 k^2/2m + 2\bar n V''_{\rm LHY})$, whose most unstable mode has $k_{\max} \propto |V''_{\rm LHY}|^{1/2} \propto (v t)^{1/2}$; evaluated at $t=t_f$, this gives droplet separation and radius proportional to $v^{-1/3}$.
What would settle it
A many-droplet experiment with potassium-39 ramping the magnetic field through the Feshbach resonance near 56.85 G at rates from 0.004 to 0.02 ms$^{-1}$ should show droplet radii shrinking by about a factor of $5^{1/3}\approx 1.7$ as the rate increases; a fitted exponent clearly outside $0.327$--$0.375$, a non-power-law, or radii set by total atom number rather than ramp rate would refute the claimed scaling.
Extended reading notes
Core claim
The paper's central claim is that, at zero temperature, a multiple-droplet system has a preferred scale proportional to $v^{-1/3}$, where $v$ is the rate of change of parameters at droplet formation, and that this scaling is independent of the detailed form of the quantum corrections and conjectured to be independent of dimension. The author derives the exponent by identifying the droplet separation with the wavelength of the fastest-growing Bogoliubov mode evaluated at freeze-out: the freeze-out time follows from rescaling time as $\tilde t = v^{1/3}t$ in the nonlinear mode equation, and the most unstable mode has $k_{\max} \propto (v t)^{1/2}$, giving $r_s \propto (v t_f)^{-1/2} \propto v^{-1/3}$. Numerical simulations of a two-dimensional Gross-Pitaevskii system with beyond-mean-field corrections support the law, with fitted exponents $d\in(0.327,0.375)$ for droplet radii and separations over a range of ramp rates.
Load-bearing premise
The argument assumes that the droplet separation is set by the wavelength of the fastest-growing Bogoliubov mode at freeze-out, with freeze-out defined by the time-rescaled nonlinear equation; if mergers, thermal noise, or nonlinear saturation set the dominant scale instead, the $v^{-1/3}$ exponent would not be robust.
Editorial extensions
If this is right
- The same $v^{-1/3}$ should appear in both droplet radius and droplet separation, so either observable can test the prediction in a many-droplet run.
- The scaling holds for ramp rates within a window where many droplets form; at very slow ramps the droplet number becomes small and at very fast ramps the freeze-out picture breaks down, so the law's range is limited by these constraints.
- Because the exponent is derived without relying on the detailed form of the quantum correction, the law should transfer to other droplet-forming atomic mixtures and, by conjecture, to three dimensions.
- At zero temperature the droplets form from quantum fluctuations rather than thermal noise, making the scaling a quantum analogue of the underdamped Kibble-Zurek mechanism for ordering dynamics.
Reading between the lines
- Thermal fluctuations would introduce a temperature-dependent freeze-out time, so an experiment near the few-nanokelvin limits quoted in the proposal should show a crossover away from $v^{-1/3}$; measuring that crossover would test how purely quantum the mechanism is.
- If three-body losses act faster than the freeze-out scale in three dimensions, the observed droplet radii could shrink after formation, shifting the fitted exponent below $1/3$ even while the formation law itself stays $1/3$.
- The freeze-out argument implies that more than the mean radius inherits the $v^{-1/3}$ scale, so the full droplet-size distribution and its higher moments should also scale with $v^{-1/3}$ if the law is universal.
- Balanced mixtures reduce the dynamics to one field; for imbalanced mixtures the single-mode reduction fails, and whether a $v^{-1/3}$ law survives would require a separate multi-component analysis.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes that when a binary Bose mixture is ramped through the droplet-instability threshold at a rate v, the characteristic droplet separation and radius scale as v^{-1/3} at zero temperature. The argument combines a freeze-out time obtained from a long-wavelength, time-dependent Bogoliubov mode equation with a length scale obtained from the fastest-growing mode of the full dispersion. The claim is tested by two-dimensional truncated-Wigner simulations of a potassium-39 mixture with parameters ramped across a Feshbach resonance; the fits yield exponents d in (0.327, 0.375), which the author reports as consistent with d = 1/3. An experimental realization is outlined.
Significance. If the 1/3 exponent were established, it would be a clean, falsifiable universal scaling law linking quantum droplet formation to the Kibble-Zurek mechanism, with consequences for droplet-size distributions and for analogies with oscillons. The numerical data are openly available, the simulations check lattice-size convergence, and the concrete experimental parameters are a useful strength. The main weakness is that the theoretical derivation contains an internal inconsistency in the freeze-out/length-scale matching; until that is repaired, the theoretical claim is not supported, and the numerical fits alone do not discriminate 1/3 from neighboring exponents.
major comments (3)
- [Droplet formation, Eqs. (3)-(4) and footnote 2] The central derivation is internally inconsistent. In Eq. (3), after the substitution t = v^{-1/3} t̃, the linear term becomes g n k^2 v^{2/3} t̃ δψ, so v is not removed for a fixed k; the statement that 'the solutions for any given k mode are independent of v' does not follow from the rescaling as written. If one instead uses the mode-dependent time t̃ = (g n k^2 v/m)^{1/3} t, then the freeze-out time for each mode is t_f ∝ v^{-1/3} k^{-2/3} times a dimensionless constant. Imposing the self-consistency condition k_max ∝ (v t_f)^{1/2} together with t_f ∝ v^{-1/3} k_max^{-2/3} gives t_f ∝ v^{-1/2} and k_max ∝ v^{1/4}, hence r_s ∝ v^{-1/4}, not v^{-1/3}. Footnote 2 acknowledges the mixing of the long-wavelength freeze-out with the full-dispersion length scale, but this mixing is exactly where the exponent changes.
- [Simulation, Figure 2] The claimed agreement between the numerics and the prediction d = 1/3 is not quantified: the four quoted exponents 0.327, 0.345, 0.375, and 0.343 have no error bars, and the spread around 1/3 is as large as +0.042. The authors should report confidence intervals for d, goodness-of-fit statistics, and a sensitivity study with respect to the droplet identification threshold nd/2; without these, the data support a power law but do not establish consistency with d = 1/3.
- [Droplet formation, paragraph after Eq. (4)] The step 'We may interpret the corresponding length scale as the droplet separation' is asserted rather than demonstrated. The simulations extract droplet separations at the final time, after mergers and nonlinear saturation have occurred, so the measured scale need not correspond to the inverse fastest-growing mode at freeze-out. A direct diagnostic, such as the time-resolved structure-factor peak position during the ramp, is needed to connect the measured separation to the freeze-out length scale; otherwise the exponent could be set by coarsening or by the initial fluctuation spectrum.
minor comments (6)
- [Abstract] The phrase 'has is a preferred scale' is ungrammatical and should read 'has a preferred scale'.
- [Various] There are several typographical errors: 'refered' should be 'referred' in the Conclusions, 'transtion' should be 'transition', 'timesstep' should be 'time-step', and 'resuts' should be 'results'.
- [Droplet formation, text before Eq. (3)] The name 'Sazuki' should be 'Suzuki'.
- [Eq. (1)] The name 'Gross-Pitaevsky' should be 'Gross-Pitaevskii'.
- [Simulation and Table I] The symbol n_d is used both for the density at the instability threshold and for the final droplet density in Table I and the simulation section; this overloaded notation is confusing and should be clarified.
- [Eq. (13)] The range dB/dt ∈ (0.3, 0.12) mG ms^{-1} is written in descending order, which is unconventional; the authors should also briefly explain how Eq. (13) is derived from the scattering-length dependence.
Circularity Check
No significant circularity: the d = 1/3 prediction is derived from the linearized Bogoliubov scaling and then tested against independent simulations; the only self-citation is an analogy and is not load-bearing.
full rationale
The paper's central claim is that multiple quantum droplets have a preferred scale r_d ∝ v^{-1/3}. The derivation in the 'Droplet formation' section is self-contained: Eq. (3) removes v by rescaling time, giving t_f ∝ v^{-1/3}; Eq. (2) supplies k_max ∝ (vt)^{1/2}; combining these gives r_s ∝ v^{-1/3}. The numerical simulations independently measure droplet radii and fit the exponent d, obtaining d ∈ (0.327, 0.375), which is then compared with the predicted 1/3. This is not a fitted parameter renamed as a prediction: the prediction is made before the fits and does not use the fitted d as input. The only self-citation, Ref. [7], is used as an analogy to dark-matter oscillons and plays no role in the derivation. Footnote 2 acknowledges that the freeze-out time is computed in a long-wavelength approximation while the length scale uses k_max from the full dispersion; that is a consistency or correctness concern, not a circularity, because neither step is defined in terms of the simulated droplet radii. No equation reduces to another by construction, and no load-bearing premise is imported from the author's own prior work. The derivation may be semi-quantitative and open to the skeptic's power-counting objection, but circularity is not present.
Assumptions & free parameters
free parameters (2)
- power-law exponent d of droplet radius scaling =
d = 0.327, 0.345, 0.375, 0.343 in the four fits of Fig. 2, i.e. d in (0.327, 0.375)
- droplet identification density threshold =
nd/2
assumptions (5)
- domain assumption Near the instability, the curvature V''_LHY can be replaced by -g v t.
- domain assumption The droplet separation is set by the inverse of the fastest-growing Bogoliubov mode at freeze-out, r_s proportional to 1/k_max.
- domain assumption The mean density n-bar and droplet density n_d are independent of v, so r_d is proportional to r_sep.
- domain assumption The effective 2D Gross-Pitaevskii equation with logarithmic quantum correction describes balanced Bose-Bose mixtures near instability.
- domain assumption Truncated Wigner Gaussian initial conditions faithfully represent quantum vacuum fluctuations.
Cite this review
Pith. "Pith review of Universal scaling law for quantum droplet formation." pith.science (2026). https://pith.science/paper/DX3DONAD
@misc{pith2026250421641,
author = {Pith},
title = {Pith review of: Universal scaling law for quantum droplet formation},
year = {2026},
howpublished = {\url{https://pith.science/paper/DX3DONAD}},
note = {Machine review of arXiv:2504.21641}
}
abstract
Given the right set of circumstances, ultracold quantum gases are able to change character and condense into a liquid state of quantum droplets. The size distribution of the droplets is determined dynamically in the condensation process. A semi-quantitative argument is presented which suggests that, at zero temperature, a multiple droplet system has is a preferred scale $\propto v^{-1/3}$, where $v$ is the rate of change of parameters at the time of droplet formation. Numerical simulations of two dimensional systems strongly support a power law $v^{-d}$, with an exponent $d\in(0.327,0.375)$.
Figures
Reference graph
Works this paper leans on
-
[1]
Quantum mechanical stabilization of a col- lapsing bose-bose mixture,
D. S. Petrov, “Quantum mechanical stabilization of a col- lapsing bose-bose mixture,” Phys. Rev. Lett.115, 155302 (2015)
2015
-
[2]
Ultradilute low-dimensional liquids,
D. S. Petrov and G. E. Astrakharchik, “Ultradilute low-dimensional liquids,” Phys. Rev. Lett. 117, 100401 (2016)
work page 2016
-
[3]
Self-bound quantum droplets of atomic mixtures in free space,
G. Semeghini, G. Ferioli, L. Masi, C. Mazzinghi, L. Wol- swijk, F. Minardi, M. Modugno, G. Modugno, M. Ingus- cio, and M. Fattori, “Self-bound quantum droplets of atomic mixtures in free space,” Physical Review Letters 120 (2018), 10.1103/physrevlett.120.235301
-
[5]
Quantum liquid droplets in a mixture of bose-einstein condensates,
C. R. Cabrera, L. Tanzi, J. Sanz, B. Naylor, P. Thomas, P. Cheiney, and L. Tarruell, “Quantum liquid droplets in a mixture of bose-einstein condensates,” Science 359, 301–304 (2018)
work page 2018
-
[6]
Colli- sions of self-bound quantum droplets,
Giovanni Ferioli, Giulia Semeghini, Leonardo Masi, Gio- vanni Giusti, Giovanni Modugno, Massimo Inguscio, Al- bert Gallem´ ı, Alessio Recati, and Marco Fattori, “Colli- sions of self-bound quantum droplets,” Phys. Rev. Lett. 122, 090401 (2019)
work page 2019
-
[7]
Ian G. Moss, “Dark matter droplets,” (2024), arXiv:2407.13243 [astro-ph.CO]
arXiv 2024
-
[8]
Novel pulsons (or stability from instability),
V. G. Makhankov, G. Kummer, and A. B. Shvachka, “Novel pulsons (or stability from instability),” Physica D: Nonlinear Phenomena 3, 344–349 (1979)
work page 1979
-
[9]
Phase transitions with subcritical bubbles,
Marcelo Gleiser, Edward W. Kolb, and Richard Watkins, “Phase transitions with subcritical bubbles,” Nucl. Phys. B 364, 411–450 (1991)
work page 1991
Show all 28 references
-
[10]
Pseudostable bubbles,
Marcelo Gleiser, “Pseudostable bubbles,” Phys. Rev. D 49, 2978–2981 (1994), arXiv:hep-ph/9308279
1994 arXiv
-
[11]
Topology of cosmic domains and strings,
T W B Kibble, “Topology of cosmic domains and strings,” Journal of Physics A: Mathematical and Gen- eral 9, 1387 (1976)
1976
-
[12]
Cosmological Experiments in Superfluid Helium?
W. H. Zurek, “Cosmological Experiments in Superfluid Helium?” Nature 317, 505–508 (1985)
1985
-
[13]
Universality of phase transition dynamics: Topological defects from symmetry breaking,
Adolfo del Campo and Wojciech H. Zurek, “Universality of phase transition dynamics: Topological defects from symmetry breaking,” International Journal of Modern Physics A 29, 1430018 (2014)
2014
-
[14]
Dynamics of the order parameter in symmetry breaking phase transi- tions,
Fumika Suzuki and Wojciech H. Zurek, “Dynamics of the order parameter in symmetry breaking phase transi- tions,” (2024), arXiv:2412.15568 [cond-mat.stat-mech]
2024
-
[15]
Off- equilibrium scaling behaviors across first-order transi- tions,
Haralambos Panagopoulos and Ettore Vicari, “Off- equilibrium scaling behaviors across first-order transi- tions,” Phys. Rev. E92, 062107 (2015), arXiv:1508.02503 [cond-mat.stat-mech]
2015 arXiv
-
[16]
Emergence and scaling of spin turbulence in quenched antiferromagnetic spinor bose-einstein con- densates,
Seji Kang, Sang Won Seo, Joon Hyun Kim, and Y. Shin, “Emergence and scaling of spin turbulence in quenched antiferromagnetic spinor bose-einstein con- densates,” Physical Review A 95 (2017), 10.1103/phys- reva.95.053638
2017 doi
-
[17]
Renormalization-group the- ory for cooling first-order phase transitions in potts mod- els,
Ning Liang and Fan Zhong, “Renormalization-group the- ory for cooling first-order phase transitions in potts mod- els,” Phys. Rev. E 95, 032124 (2017)
2017
-
[18]
Dynamics of the quantum search and quench- induced first-order phase transitions,
Ivan B. Coulamy, Andreia Saguia, and Marcelo S. Sarandy, “Dynamics of the quantum search and quench- induced first-order phase transitions,” Phys. Rev. E 95, 022127 (2017)
2017
-
[19]
Observation of gener- alized Kibble-Zurek mechanism across a first-order quan- tum phase transition in a spinor condensate,
L. Y. Qiu, H. Y. Liang, Y. B. Yang, H. X. Yang, T. Tian, Y. Xu, and L. M. Duan, “Observation of gener- alized Kibble-Zurek mechanism across a first-order quan- tum phase transition in a spinor condensate,” Science Advances 6, eaba7292 (2020), arXiv:2001.10210 [cond- mat.quant-gas]
2020 arXiv
-
[20]
Eigenvalues and eigenfunctions of a bose system of hard spheres and its low-temperature properties,
T. D. Lee, Kerson Huang, and C. N. Yang, “Eigenvalues and eigenfunctions of a bose system of hard spheres and its low-temperature properties,” Phys. Rev. 106, 1135– 1145 (1957)
1957
-
[21]
Three-body losses in trapped bose-einstein-condensed gases,
Yeong E. Kim and Alexander L. Zubarev, “Three-body losses in trapped bose-einstein-condensed gases,” Phys. Rev. A 69, 023602 (2004). 5
2004
-
[22]
Dimensional crossover for the beyond-mean-field correction in bose gases,
Tobias Ilg, Jan Kumlin, Luis Santos, Dmitry S. Petrov, and Hans Peter B¨ uchler, “Dimensional crossover for the beyond-mean-field correction in bose gases,” Phys. Rev. A 98, 051604 (2018)
2018
-
[23]
Phases and dynamics of quantum droplets in the crossover to two- dimensions,
Jose Carlos Pelayo, George Bougas, Thom´ as Fogarty, Thomas Busch, and Simeon I. Mistakidis, “Phases and dynamics of quantum droplets in the crossover to two- dimensions,” (2025), arXiv:2407.16383 [cond-mat.quant- gas]
2025 arXiv
-
[24]
Feshbach resonances in ultracold39k,
Chiara D’Errico, Matteo Zaccanti, Marco Fattori, Gia- como Roati, Massimo Inguscio, Giovanni Modugno, and Andrea Simoni, “Feshbach resonances in ultracold39k,” New Journal of Physics 9, 223–223 (2007)
2007
-
[25]
Test of the universality of the three-body efimov parameter at narrow feshbach resonances,
Sanjukta Roy, Manuele Landini, Andreas Trenkwalder, Giulia Semeghini, Giacomo Spagnolli, Andrea Simoni, Marco Fattori, Massimo Inguscio, and Giovanni Mod- ugno, “Test of the universality of the three-body efimov parameter at narrow feshbach resonances,” Phys. Rev. Lett. 111, 0...
2013
-
[26]
Feshbach resonances in potassium bose-bose mixtures,
L. Tanzi, C. R. Cabrera, J. Sanz, P. Cheiney, M. Tomza, and L. Tarruell, “Feshbach resonances in potassium bose-bose mixtures,” Physical Review A 98 (2018), 10.1103/physreva.98.062712
2018 doi
-
[27]
Trapped imbalanced quantum droplets,
Thomas A. Flynn, Nick A. Keepfer, Nick G. Parker, and Thomas P. Billam, “Trapped imbalanced quantum droplets,” (2023), arXiv:2309.04300 [cond-mat.quant- gas]
2023 arXiv
-
[28]
Data supporting this publication,
“Data supporting this publication,” 10.25405/data.ncl.28956770. SUPPLEMENT AR Y MA TERIAL Renormalisation It is instructive to analyse the renormal- isation of the quantum vacuum energy in two dimensions which has been used in the main text. We start with sep- arate number den...
-
[29]
Figure 3
Numerical simulations were also run with white noise initial conditions⟨ ¯ψ(k)ψ(k′)⟩ = 1 2δkk′ with no significant change to the results. Figure 3. The histogram shows the probability distribution of atom numbers in each lump. (Compiled from five separate runs with v = 0.01ms−...
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.