REVIEW 2 major objections 4 minor 48 references
Suppression of capillary instability in a confined quantum liquid filament
T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Transverse harmonic confinement fully suppresses the Rayleigh–Plateau instability of a quantum liquid filament once the trap frequency exceeds the inverse capillary time.
desk verdict A solid numerical study showing confinement suppresses capillary instability in a quantum filament, with a quantitative criterion that needs one more benchmark. 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 generalized capillary dispersion relation of Eq. (26), adapted from a classical hydrodynamic result for a liquid cylinder subjected to a radial force. It turns the infinite-filament Rayleigh–Plateau criterion $kR<1$ into $kR < \sqrt{1-(\Omega\tau_c)^2}$, making the capillary growth rate a sensitive function of the confinement frequency. Supporting that relation are the Bogoliubov–de Gennes linear-response equations (14)–(17) built from an effective single-component functional with a Lee–Huang–Yang quantum-fluctuation term, which the paper benchmarks against coupled two-component Gross–Pitaevskii simulations.
What would settle it
A direct experimental test would be to measure the growth rate of axial density modulations in a 41K–87Rb droplet filament inside a waveguide as a function of trap frequency; observing any imaginary-frequency mode (breakup) at trap frequencies above $\Omega_c = \tau_c^{-1}$, or a growth-rate spectrum inconsistent with Eq. (26) at fixed fitted radius and capillary time, would refute the central claim.
Extended reading notes
Core claim
The paper establishes that the unstable spectrum of a quantum liquid filament in free space matches the classical Rayleigh–Plateau dispersion relation once the fitted radius and capillary time are used, and that under transverse harmonic confinement the growth rate obeys a generalized relation in which the trap frequency enters as $\omega^2_{RP} = \tau_c^{-2} [I_1(kR)/I_0(kR)]\, kR\, [1-(kR)^2-(\Omega\tau_c)^2]$. As a result the unstable band shrinks to $kR < \sqrt{1-(\Omega\tau_c)^2}$ and vanishes when $\Omega = \tau_c^{-1}$. The paper interprets this as the static-harmonic-confinement analogue of the classical dynamic stabilization of a liquid cylinder under an oscillatory radial force, and it validates the BdG spectra against imaginary-time plus real-time Gross–Pitaevskii simulations of the two-component mixture.
Load-bearing premise
The paper's key premise is that the classical dispersion relation for a liquid cylinder under an oscillatory radial force remains valid when the oscillatory force is replaced by a static harmonic trap and applied to a compressible, inhomogeneous quantum filament.
Editorial extensions
If this is right
- The free-space filament spectra collapse onto the classical universal Rayleigh–Plateau curve once rescaled by the fitted radius and capillary time, confirming the capillary interpretation even at low linear densities where no flat-top bulk exists.
- Increasing the transverse trap frequency narrows the unstable wavevector band, shifting the most-unstable mode toward smaller wavenumbers and smaller growth rates.
- Above the critical frequency $\Omega_c=\tau_c^{-1}$ the filament has no imaginary-frequency BdG modes, so it is fully stable against capillary breakup.
- The critical frequency decreases with increasing linear density, since $\tau_c$ grows with $N/L$.
- The fitted radius and capillary time from the confined spectra agree with estimates from the density profiles, supporting the hydrodynamic analogy even in the presence of quantum corrections and confinement.
Reading between the lines
- If the predicted threshold $\Omega_c=\tau_c^{-1}$ holds quantitatively, a waveguide experiment could map the full transition by varying the trap frequency and measuring the breakup time, giving a direct measurement of the capillary time from the stabilization point.
- The stabilizing radial-force mechanism suggests that time-periodic or anharmonic radial potentials, or even rotation of an anisotropic trap, could also shift the capillary spectrum; the paper leaves these geometries open.
- The analogy with classical dynamic stabilization hints that a sufficiently strong radial force might stabilize not only the axisymmetric mode but also higher azimuthal modes, which the current axial-polarization analysis does not probe.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies an infinite cylindrical self-bound Bose-Bose mixture filament (41K-87Rb) described by an effective single-component GP functional with the LHY quantum correction, with and without transverse harmonic confinement. The authors solve the BdG equations for axial modes, compare the free-space unstable spectra with full two-component GP simulations, and show that the spectra collapse onto the classical Rayleigh-Plateau dispersion when rescaled by fitted radius R and capillary time tau_c. For a radially trapped filament, they find that increasing trap frequency Omega narrows the unstable band and eventually removes it, which they interpret via a generalized dispersion relation (Eq. 26) adapted from a classical inviscid-liquid result (Ref. [26]) and predict complete stabilization for Omega > Omega_c = tau_c^{-1} (Eq. 28). The fitted R and tau_c in the confined case are compared with density-profile estimates (Fig. 8).
Significance. The free-space part of the paper is convincing and valuable: the BdG spectra are explicitly benchmarked against full two-component GP simulations (Fig. 3), and the fitted R and tau_c agree with independent density-profile estimates (Fig. 5), giving a quantitative connection between quantum liquid filaments and classical capillary instability. If the central confined-filament claim holds, the paper provides a simple, experimentally testable criterion (Omega_c = tau_c^{-1}) for completely suppressing Rayleigh-Plateau breakup in quantum liquid filaments by transverse waveguide confinement. However, the confined-case analysis is not on the same footing: the key dispersion relation Eq. (26) is taken from a classical hydrodynamic model for an incompressible uniform-density liquid and adapted to a static harmonic trap without derivation, and the trapped BdG spectra are fitted using that same expression, so the predicted stabilization threshold is not independently validated.
major comments (2)
- [Section V, Eq. (26)] The generalized dispersion relation is asserted by adapting Eq. (2.3) of Ref. [26] from oscillatory forcing to static harmonic confinement, but no derivation is provided for the present compressible, inhomogeneous quantum filament. Since Eq. (26) is the direct source of the central prediction Omega_c = tau_c^{-1} (Eq. (28)), and since the BdG spectra in Fig. 6 are fitted using Eq. (26) with R and tau_c as free parameters, the agreement shown in Fig. 6 cannot independently validate the functional form. The comparison in Fig. 8 checks only that the best-fit R and tau_c are close to density-based estimates, which would also be consistent with a different Omega-dependence if the fit compensates through R and tau_c. The authors should derive Eq. (26) within the effective single-component framework or benchmark the confined BdG results against full two-component GP simulations near the predicted Omega_c, as was done for the free-space case in Sec. IV B.
- [Section V, Eq. (27) and Appendix A] The full-stabilization claim requires resolving the k -> 0 limit. According to Eq. (27), as Omega approaches Omega_c the unstable band shrinks to kR -> 0, i.e., the instability is restricted to arbitrarily long wavelengths. The numerical BdG calculations use a periodic supercell of finite axial length L, so the resolved wavevectors are quantized (k_n = 2 pi n / L) unless an explicit continuum-k procedure is used. The manuscript should demonstrate that the vanishing of unstable modes at large Omega is not a finite-box artifact, by showing convergence of the BdG spectra in L (or in k resolution) for frequencies near Omega_c. Ideally, a dynamical GP simulation in a sufficiently long box at Omega just above the predicted Omega_c would directly confirm that the filament does not break up.
minor comments (4)
- [Fig. 6 caption] The caption reads 'The solid lines rare fits using Eq. (26)'; this should be 'are fits'.
- [Section V heading] The heading 'RADIALL Y CONFINED FILAMENT' contains a spacing typo and should read 'RADIALLY CONFINED FILAMENT'.
- [Section V, final paragraph] The phrase 'This scenario clear contrasts' should be 'This scenario clearly contrasts'.
- [Fig. 7 caption] The caption states that the curves are 'obtained from Eq. 26'; it would be clearer to state explicitly that these are the classical predictions, not the BdG numerical results, so the reader does not mistake them for independent evidence.
Circularity Check
No construction-level circularity: the BdG spectra, not the fits, show that confinement suppresses the instability; Eq. 26 is a transparent external-model adaptation, and the fitted Rayleigh-Plateau parameters are benchmarked against independent density-based estimates.
full rationale
The central result, namely the progressive narrowing and eventual disappearance of the unstable band under transverse confinement, is an output of the BdG eigenvalue problem (Eqs. 14-17) solved on the stationary filament obtained from Eq. (12); no fitting parameter enters the spectra displayed in Fig. 6. The Rayleigh-Plateau formula (Eq. 22) and its confined variant (Eq. 26) are introduced afterwards as interpretive fits: as the paper states, 'the spectra of the unstable modes of the trapped filament can be fitted using Eq. (26) to extract the parameters R and tau_c, which can then be compared with the corresponding estimates obtained from the calculated density profiles.' Because R and tau_c are free parameters in the fit, the collapse onto the classical curve is a shape and consistency check, not a construction of the spectrum; the independent estimates from the density profile (Eqs. 24-25) provide external validation. Equation 26 is adapted from Ref. [26] (Patankar, Basak, and Dasgupta, JFM 2022), an external classical-fluid work, and the adaptation is explicit ('By adapting the results of Ref. [26]...'), including the caveat that it 'was originally derived for an incompressible liquid filament with uniform density.' Thus the critical frequency Omega_c = tau_c^{-1} is an analytic consequence of a transparent model assumption, not a self-referential definition. The self-citations used for the single-component functional and surface-tension estimates are not the sole support of the central claim, because the single-component BdG approach is benchmarked against full two-component GP simulations (Fig. 1(b) and Fig. 3(a)). The main scientific risk is that the Omega-dependence of Eq. 26 is not independently validated for a compressible, inhomogeneous quantum filament, but that is a correctness and robustness question, not circularity.
Assumptions & free parameters
free parameters (3)
- Filament radius R (fit parameter) =
Values vary with linear density; for example, roughly 2 to 7 micrometers for the highest density and below 1…
- Capillary time tau_c (fit parameter) =
Values vary with linear density and confinement; roughly 1 to 8 ms across the studied range.
- Surface tension prefactor alpha =
2
assumptions (5)
- domain assumption The density ratio between components stays at the bulk equilibrium value eta = sqrt(g22/g11) throughout the filament, including the surface region.
- standard math The LHY energy functional can be represented by the fitting function f(m2/m1,1,xi) of Eq. (4) with coefficients from Ref [37].
- standard math BdG linearization around the stationary filament gives the full excitation spectrum, and the periodic-supercell method with a finite axial length L represents the infinite filament.
- domain assumption The classical confined Rayleigh-Plateau spectrum of Ref [26] Eq. (2.3), derived for oscillatory forcing, applies to static harmonic confinement with the form of Eq. (26).
- domain assumption Surface tension can be estimated from the quantum kinetic energy of the density variation via Eq. (25) with alpha=2.
Cite this review
Pith. "Pith review of Suppression of capillary instability in a confined quantum liquid filament." pith.science (2026). https://pith.science/paper/CFNN2YSC
@misc{pith2026250711223,
author = {Pith},
title = {Pith review of: Suppression of capillary instability in a confined quantum liquid filament},
year = {2026},
howpublished = {\url{https://pith.science/paper/CFNN2YSC}},
note = {Machine review of arXiv:2507.11223}
}
read the original abstract
Quantum Bose-Bose mixtures with strong attraction can form self-bound, liquid-like droplets stabilized by quantum fluctuations. Despite equilibrium densities much lower than those of classical liquids, these droplets exhibit finite surface tension and liquid-like behaviors. Recent experiments have demonstrated Rayleigh-Plateau instability in elongated droplets confined in an optical waveguide. Here we consider the case of an infinite filament and extend the theoretical description to include transverse harmonic confinement. By solving the Bogoliubov-deGennes equations within a single-component framework, benchmarked against full Gross-Pitaevskii simulations, we show that increasing confinement progressively suppresses the instability, leading to complete stabilization beyond a critical trap frequency.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[31]
F. Ancilotto, M. Barranco, and M. Pi, Breakup of quan- tum liquid filaments into droplets, Phys. Rev. A 107, 063312 (2023)
work page 2023
-
[26]
S. Patankar, S. Basak, and R. Dasgupta, Dynamic stabil- isation of Rayleigh-Plateau modes on a liquid cylinder, J. Fluid Mech. 946, A2 (2022)
work page 2022
-
[1]
D. S. Petrov, Quantum Mechanical Stabilization of a Col- lapsing Bose-Bose Mixture, Phys. Rev. Lett. 115, 155302 (2015). 10
work page 2015
- [2]
-
[3]
Ferrier-Barbut, H
I. Ferrier-Barbut, H. Kadau, M. Schmitt, M. Wenzel, and T. Pfau, Observation of Quantum Droplets in a Strongly Dipolar Bose Gas, Phys. Rev. Lett. 116, 215301 (2016)
2016
-
[4]
Schmitt, M
M. Schmitt, M. Wenzel, F. B¨ ottcher, I. Ferrier-Barbut, and T. Pfau, Self-bound droplets of a dilute magnetic quantum liquid, Nature 539, 259 (2016)
2016
- [5]
-
[6]
I. Ferrier-Barbut, M. Schmitt, M. Wenzel, H. Kadau, and T. Pfau, Liquid quantum droplets of ultracold magnetic atoms, J. Phys. B: At. Mol. Opt. Phys.49, 214004 (2016)
work page 2016
Show all 48 references
-
[7]
Wenzel, F
M. Wenzel, F. B¨ ottcher, T. Langen, I. Ferrier-Barbut, and T. Pfau, Striped states in a many-body system of tilted dipoles, Phys. Rev. A 96, 053630 (2017)
2017
-
[8]
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 (2018)
2018
-
[9]
Semeghini, G
G. Semeghini, G. Ferioli, L. Masi, C. Mazzinghi, L. Wol- swijk, F. Minardi, M. Modugno, G. Modugno, M. In- guscio, and M. Fattori, Self-Bound Quantum Droplets of Atomic Mixtures in Free Space, Phys. Rev. Lett. 120, 235301 (2018)
2018
-
[10]
D’Errico, A
C. D’Errico, A. Burchianti, M. Prevedelli, L. Salasnich, F. Ancilotto, M. Modugno, F. Minardi, and C. Fort, Ob- servation of quantum droplets in a heteronuclear bosonic mixture, Phys. Rev. Res. 1, 033155 (2019)
2019
-
[11]
Burchianti, C
A. Burchianti, C. D’Errico, M. Prevedelli, L. Salasnich, F. Ancilotto, M. Modugno, F. Minardi, and C. Fort, A Dual-Species Bose-Einstein Condensate with Attractive Interspecies Interactions, Condens. Matter 5, 21 (2020)
2020
-
[12]
Z. Guo, F. Jia, L. Li, Y. Ma, J. M. Hutson, X. Cui, and D. Wang, Lee-Huang-Yang effects in the ultracold mixture of 23Na and 87Rb with attractive interspecies interactions, Phys. Rev. Res. 3, 033247 (2021)
2021
-
[13]
Cavicchioli, C
L. Cavicchioli, C. Fort, F. Ancilotto, M. Modugno, F. Mi- nardi, and A. Burchianti, Dynamical Formation of Mul- tiple Quantum Droplets in a Bose-Bose Mixture, Phys. Rev. Lett. 134, 093401 (2025)
2025
-
[14]
Plateau, Experimental and theoretical researches on the figures of equilibrium of a liquid mass withdrawn from the action of gravity, Philos
J. Plateau, Experimental and theoretical researches on the figures of equilibrium of a liquid mass withdrawn from the action of gravity, Philos. Mag. 14, 431 (1857)
-
[15]
Rayleigh, On The Instability Of Jets, Proc
L. Rayleigh, On The Instability Of Jets, Proc. Lond. Math. Soc. s1-10, 4 (1878)
-
[16]
Rayleigh, On the Stability, or Instability, of certain Fluid Motions, Proc
L. Rayleigh, On the Stability, or Instability, of certain Fluid Motions, Proc. Lond. Math. Soc. s1-11, 57 (1879)
-
[17]
Rayleigh, On the instability of a cylinder of viscous liquid under capillary force, Philos
L. Rayleigh, On the instability of a cylinder of viscous liquid under capillary force, Philos. Mag. 34, 145 (1892)
-
[18]
A. A. Castrej´ on-Pita, J. R. Castrej´ on-Pita, and I. M. Hutchings, Breakup of Liquid Filaments, Phy. Rev. Lett. 108, 074506 (2012)
2012
-
[19]
Tomotika, On the instability of a cylindrical thread of a viscous liquid surrounded by another viscous fluid, Proc
S. Tomotika, On the instability of a cylindrical thread of a viscous liquid surrounded by another viscous fluid, Proc. R. Soc. Lond. 150, 322 (1935)
1935
-
[20]
S. L. Goren, The instability of an annular thread of fluid, J. Fluid Mech. 12, 309–319 (1962)
1962
-
[21]
P. S. Hammond, Nonlinear adjustment of a thin annu- lar film of viscous fluid surrounding a thread of another within a circular cylindrical pipe, J. Fluid Mech. 137, 363–384 (1983)
1983
-
[22]
C. Zhao, Y. Zhang, and T. Si, Slip-enhanced Rayleigh- Plateau instability of a liquid film on a fibre, J. Fluid Mech. 954, A46 (2023)
2023
-
[23]
Brosa, S
U. Brosa, S. Grossmann, and A. M¨ uller, Nuclear scission, Phys. Rep. 197, 167 (1990)
1990
-
[24]
Shabahang, J
S. Shabahang, J. J. Kaufman, D. S. Deng, and A. F. Abouraddy, Observation of the Plateau-Rayleigh cap- illary instability in multi-material optical fibers, Appl. Phys. Lett. 99, 161909 (2011)
2011
-
[25]
Eggers and E
J. Eggers and E. Villermaux, Physics of liquid jets, Rep. Prog. Phys. 71, 036601 (2008)
2008
-
[27]
Patrascu and C
C. Patrascu and C. Balan, The stabilizing effect of con- finement on a liquid jet in a viscous outer fluid, U.P.B. Sci. Bull., Series A 81, 85 (2019)
2019
-
[28]
N. B. Speirs, K. R. Langley, P. Taborek, and S. T. Thoroddsen, Jet breakup in superfluid and normal liq- uid 4He, Phys. Rev. Fluids 5, 044001 (2020)
2020
-
[29]
Ancilotto, M
F. Ancilotto, M. Barranco, and M. Pi, Nanoscopic jets and filaments of superfluid 4He at zero temperature: A DFT study, J. Chem. Phys. 158, 144306 (2023)
2023
-
[30]
Sasaki, N
K. Sasaki, N. Suzuki, and H. Saito, Capillary instabil- ity in a two-component Bose-Einstein condensate, Phys. Rev. A 83, 053606 (2011)
2011
-
[32]
Hernandez-Rajkov, N
D. Hernandez-Rajkov, N. Grani, F. Scazza, G. Del Pace, W. J. Kwon, M. Inguscio, K. Xhani, C. Fort, M. Mod- ugno, F. Marino, and G. Roati, Connecting shear flow and vortex array instabilities in annular atomic superflu- ids, Nat. Phys. 20, 939 (2024)
2024
-
[33]
S. Huh, W. Yun, G. Yun, S. Hwang, K. Kwon, J. Hur, S. Lee, H. Takeuchi, S. K. Kim, and J. Choi, Beyond skyrmion spin texture from quantum Kelvin-Helmholtz instability, arXiv:2408.11217 (2024)
2024 arXiv
-
[34]
Y. Geng, J. Tao, M. Zhao, S. Mukherjee, S. Eckel, G. K. Campbell, and I. B. Spielman, The Rayleigh-Taylor in- stability in a binary quantum fluid, arXiv:2411.19807 (2024)
2024 arXiv
-
[35]
The true ground state would, instead, correspond to a spher- ical droplet in free space or to an elongated droplet in a radial waveguide
The infinitely extented filament is a stationary state of the system and is therefore subject to instability. The true ground state would, instead, correspond to a spher- ical droplet in free space or to an elongated droplet in a radial waveguide
-
[36]
Ancilotto, M
F. Ancilotto, M. Barranco, M. Guilleumas, and M. Pi, Self-bound ultradilute Bose mixtures within local density approximation, Phys. Rev. A 98, 053623 (2018)
2018
-
[37]
Minardi, F
F. Minardi, F. Ancilotto, A. Burchianti, C. D’Errico, C. Fort, and M. Modugno, Effective expression of the Lee-Huang-Yang energy functional for heteronuclear mixtures, Phys. Rev. A 100, 063636 (2019)
2019
-
[38]
Dalfovo, S
F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Theory of Bose-Einstein condensation in trapped gases, Rev. Mod. Phys. 71, 463 (1999)
1999
-
[39]
N. B. Jørgensen, G. M. Bruun, and J. J. Arlt, Dilute fluid governed by quantum fluctuations, Phys. Rev. Lett.121, 173403 (2018)
2018
-
[40]
Cikojevi´ c, E
V. Cikojevi´ c, E. Poli, F. Ancilotto, L. Vranjeˇ s-Marki´ c, 11 and J. Boronat, Dilute quantum liquid in a K-Rb Bose mixture, Phys. Rev. A 104, 033319 (2021)
2021
-
[41]
A cylindrical filament–rather than a spherical droplet–can be obtained by choosing suitable initial conditions
Stationary configurations of the energy functional (3) are obtained using a steepest descent algorithm [47], which iteratively propagates the two condensate components in imaginary time. A cylindrical filament–rather than a spherical droplet–can be obtained by choosing suitabl...
-
[42]
The time-dependent GP equations are solved using a split-step Fourier method; see, e.g., Ref. [48]
-
[43]
Using a different definition of the radial size does not affect the results
-
[44]
Ao and S
P. Ao and S. T. Chui, Binary Bose-Einstein condensate mixtures in weakly and strongly segregated phases, Phys. Rev. A 58, 4836 (1998)
1998
-
[45]
Jimbo and H
K. Jimbo and H. Saito, Surfactant behavior in three- component Bose-Einstein condensates, Physical Review A 103, 063323 (2021), 2104.00313
2021 arXiv
-
[46]
It is worth noting that the variations of the surface ten- sion with the filament’s linear density, as well as those arising from the harmonic radial confinement (see Sec.V), are found to be negligible within the investigated param- eter range
-
[47]
W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes: The Art of Scientific Com- puting, 3rd ed. (Cambridge University Press, USA, 2007)
2007
-
[48]
Jackson, J
B. Jackson, J. F. McCann, and C. S. Adams, Output coupling and flow of a dilute Bose-Einstein condensate, J. Phys. B: At. Mol. Opt. Phys. 31, 4489 (1998)
1998
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.