Pith. sign in

REVIEW 3 major objections 3 minor 33 references

Excitation spectrum of a double supersolid in a trapped dipolar Bose mixture

T0 review · 3 major / 3 minor · reviewed 2026-08-11 · deepseek-v4-flash

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

Pith's one-line read In the double-supersolid phase of a trapped dipolar Bose mixture, axial compression should split the superfluid breathing response into two single-component peaks, directly revealing two coexisting superfluids with different superfluid…

desk verdict A testable prediction—the compressional doublet—in a solid numerical study, but the generality claim rests on one parameter set and the manuscript needs cleanup. read the letter →

arxiv 2412.05215 v1 pith:YF3FR7B7 submitted 2024-12-06 cond-mat.quant-gas

Daniel Scheiermann , Albert Gallemí , Luis Santos This is my paper · ORCID
classification cond-mat.quant-gas PACS 03.75.Kk03.75.Mn67.85.-d
keywords dipolarBosemixturesdoublesupersolidBogoliubov-deGennesequationscollectiveexcitationsbreathingmodessuperfluidfractionincoherentdropletregime
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

The paper predicts a distinctive fingerprint for the double-supersolid phase of a trapped, miscible dipolar Bose mixture: the low-lying excitation spectrum rearranges into a three-mode structure in which the compressional (breathing) response splits into two superfluid modes, each carried almost entirely by one of the two components. Because a simple axial compression excites both modes, measuring that response would directly reveal that two interacting superfluids coexist and that their superfluid fractions differ markedly. The same spectra also show the roton, Higgs, and Goldstone modes acquiring a strong density-spin hybridization in an asymmetric mixture, and they let one watch a single component cross into the incoherent droplet regime while the other remains superfluid.

What carries the argument

The argument is carried by the extended Gross-Pitaevskii equations for the two components, which include the two-component Lee-Huang-Yang energy density from quantum fluctuations, together with their linearization into Bogoliubov-de Gennes equations. Each eigenmode is then classified by two observables: $Q$, the relative weight of total-density versus spin (relative-density) modulation, and $P$, the relative contribution of each component, with $P=1$ ($P=-1$) marking a mode of only component 1 (2). The phase-fluctuation strengths $\eta_\sigma$ and $\lambda_\sigma$ further detect the supersolid-to-droplet crossover, since a component entering the incoherent droplet regime develops phase variations confined between droplets, so $\lambda_\sigma/\eta_\sigma$ approaches zero. The doublet claim follows from applying an equal $x^2$ compression to both components and computing the normalized response $\bar S(\omega)$ separately for each component; the interpretation of the lower-energy peak as the lower superfluid fraction relies on a cited upper bound on superfluid fraction.

What would settle it

Measure the component-resolved axial compressional response of an asymmetric dipolar mixture in the double-supersolid regime: the claim requires two distinct superfluid breathing peaks with $P\simeq\pm1$ (each carried by one component) and the lower-frequency peak belonging to the more contrasted component. Observing comparable component weights in both peaks, a single unsplit peak, or the lower-frequency peak belonging to the less contrasted component would contradict the central prediction. The predicted location of the single-component incoherent transition ($a_{12}\simeq65$–$70\,a_0$ from $\lambda_1/\eta_1\to0$) can likewise be checked against the measured phase coherence between the two central droplets.

Watch

Extended reading notes

Core claim

On its own terms, the paper establishes that for a miscible mixture of two dipolar condensates with unequal dipole moments (illustrated by $^{162}\mathrm{Dy}$ components with $\mu_1 = 10\,\mu_B$ and $\mu_2 = 9\,\mu_B$), reducing the interspecies scattering length $a_{12}$ drives the system into a two-droplet double supersolid. In that phase the two gapless Goldstone modes of the unmodulated mixture plus the softened roton reorganize into a triplet: an in-phase dipole mode at $\omega_x$, an out-of-phase dipole mode, and a supersolid Goldstone mode with strong spin-density hybridization and a dominant weight in the more dipolar component. The axial breathing response doubles: the hardening crystal mode is roughly equally shared, while the two softer superfluid modes become almost pure single-component modes, with the lower-energy mode dominated by the more contrasted component, which has the lower superfluid fraction. The paper further shows that the lowest mode's phase-fluctuation ratio $\lambda_1/\eta_1$ collapsing toward zero marks the transition of only component 1 into the incoherent droplet regime, leaving component 2 superfluid.

Load-bearing premise

The central prediction stands on the extended Gross-Pitaevskii equations with the two-component Lee-Huang-Yang quantum-fluctuation term being quantitatively accurate at the droplet densities of this strongly dipolar mixture, including the coherence of the droplet halo near the incoherent transition.

Editorial extensions

If this is right

  • In the double-supersolid regime, an axial compression applied equally to both components should produce two distinct superfluid breathing peaks in $\bar S(\omega)$, one dominated by component 1 and one by component 2.
  • Comparing the two peak energies gives a direct, component-resolved ordering of superfluid fractions: the more contrasted (more dipolar) component carries the lower-energy mode and the lower superfluid fraction.
  • Tracking the lowest mode's phase-fluctuation ratio $\lambda_1/\eta_1$ should reveal the onset of the incoherent droplet regime for only one component, with halo decoupling near $a_{12}\simeq70\,a_0$ and droplet decoherence near $a_{12}\simeq65\,a_0$.
  • A symmetric mixture hides this physics because the compressional perturbation couples only to density modes; the doublet is therefore a probe specific to asymmetric mixtures.
  • In the ID-supersolid regime the low-lying spectrum reduces to a two-mode structure, mirroring the two Goldstone modes expected when only one component remains superfluid.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, if the doublet splitting scales monotonically with the asymmetry between components (dipole moments, populations, or intraspecies scattering lengths), the same compressional measurement could serve as a quantitative two-fluid probe of superfluid-fraction differences well beyond the specific $10\mu_B/9\mu_B$ case studied here.
  • Beyond the paper, the near-single-component character of each superfluid breathing mode suggests a sum-rule picture in which the compressional response decomposes into two nearly independent one-component oscillator strengths; checking the measured peak weights against the relative particle numbers would test whether the single-component dominance is exact or only approximate.
  • Beyond the paper, the same $Q$, $P$, $\eta_\sigma$, $\lambda_\sigma$ classification could be exported to immiscible double supersolids or dipolar-non-dipolar mixtures, where spin-density hybridization is expected to be even stronger and where the three-mode structure may leave a different experimental signature.
  • Beyond the paper, since the halo decoupling appears as an abrupt hardening and disappearance of one breathing peak from $\bar S(\omega)$, time-resolved measurements after a quench of $a_{12}$ could reveal the dynamics of the incoherent transition rather than only its static signature.
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

3 major / 3 minor

Summary. The paper studies the low-lying excitation spectrum of a trapped miscible dipolar Bose mixture using the two-component extended Gross-Pitaevskii equations with a Lee-Huang-Yang correction and the associated Bogoliubov-de Gennes equations. It focuses on an asymmetric mixture of two 162Dy components with different dipole moments (µ1=10µB, µ2=9µB), equal intra-species scattering lengths, equal populations, and a fixed three-dimensional harmonic trap. The authors introduce diagnostics P, Q, ησ, and λσ to characterize the component weight, density/spin character, and phase-fluctuation strength of each mode. Their central results are that the asymmetric double supersolid displays a three-mode low-lying structure (in-phase dipole, out-of-phase spin dipole, and a supersolid Goldstone mode), and that axial compressional excitations produce a doublet of superfluid breathing modes, each almost fully dominated by one component. The paper also identifies spectral signatures of the transition in which one component becomes an incoherent droplet array while the other remains superfluid.

Significance. If the numerical results are correct, the paper offers an experimentally accessible signature of the double-supersolid phase and of component-dependent superfluid fractions in dipolar mixtures: the compressional-response doublet shown in Fig. 6. The study uses physical input parameters (scattering lengths, dipole moments, trap frequencies, and particle number) with no parameters fitted to the output, and the predicted three-mode structure is falsifiable with existing experimental techniques. The mode-character diagnostics P, Q, ησ, and λσ are a useful extension of tools developed for single-component dipolar supersolids. The main limitations are the reliance on the approximate two-component LHY functional and the fact that the central prediction is demonstrated for a single two-droplet parameter set; neither limitation is internally inconsistent, but together they bound the strength of the paper's generality claims.

major comments (3)
  1. [Sec. VI, Figs. 3 and 6; Sec. VII] The central claim that the compressional doublet 'would generally be a clear proof' of the two-fluid character of the double supersolid is supported by exactly one parameter set: µ1=10µB, µ2=9µB, a11=a22=100a0, N1=N2=N/2, and one fixed trap geometry with a two-droplet ground state. The text in Sec. VI states that the results 'are to a large extend representative of other asymmetric mixtures,' and Sec. VII repeats the generality claim while acknowledging the restriction to two-droplet supersolids 'due to numerical complexity.' This generality is not demonstrated. Section V shows that in the symmetric limit the x2 perturbation is exactly orthogonal to the spin superfluid mode, so the appearance of two compressional superfluid modes is not symmetry-enforced; it depends on the degree of spin-density hybridization, which is parameter dependent, as the strong variation of P with a12 in Fig. 3 indicates. Please either add calculations for at least one additional asymmetry, population imbalance, droplet number, or trap geometry, or explicitly restrict the central prediction to the demonstrated regime.
  2. [Secs. IV and VI] The quantitative predictions, including the breathing-doublet frequencies, the mode-softening locations, and the value a12≈65–70a0 for the onset of the incoherent-droplet transition, are presented without any numerical methods or convergence checks. The manuscript does not specify the computational grid or box size, the imaginary-time or relaxation scheme used to obtain the ground states, the basis or eigensolver used for the BdG matrix in Eqs. (6)–(8), the number of retained modes, or tests showing that the frequencies and the diagnostics P, Q, ησ, and λσ are converged with respect to these choices. No code or data are provided. Without this information the numbers in Figs. 3, 6, and 7 cannot be independently assessed, and the claimed agreement at the transition cannot be verified.
  3. [Sec. III, Fig. 1; Sec. VI.D] The phase diagram is constructed using ad hoc density-contrast thresholds, C1=0.99 and C2=0.1 and 0.99, and no sensitivity analysis is reported. The later assignment of the double-supersolid to ID-supersolid transition at a12≈65–70a0 (Sec. VI.D, Fig. 7) is described as being in 'very good agreement' with the contrast-based phase boundary, but that comparison is only meaningful if the boundary does not shift appreciably with the chosen thresholds. Please quantify the threshold dependence or provide an independent phase criterion for the spectral identification.
minor comments (3)
  1. [Throughout, especially Secs. V–VI and Figs. 2–5] The submitted manuscript text contains extensive extraneous material on 'Topological hole localization in binary Bose mixtures in spin-dependent ladders' interleaved with the main text, together with repeated garbled equations and figure captions. This must be removed and all figures reproduced cleanly; as submitted, it prevents a reliable check of the displayed results.
  2. [Eq. (7)] The notation ̃μσ is introduced for the chemical potential obtained from the ground-state calculation, but the tilde is not used consistently elsewhere in the paper; please unify the notation.
  3. [Sec. VI.C] In the discussion of Fig. 6, it would be helpful to state explicitly how the separate response spectra for the two components are normalized and whether the relative heights in the two panels can be compared quantitatively, given the description that each S̄(ω) is normalized to the maximum value in either component.

Circularity Check

0 steps flagged · score 1.0 of 10

No circular derivation: the excitation spectra are computed from stated physical inputs with no fitted parameters; self-citations are contextual only.

full rationale

The central results (the three-mode low-lying structure, the Higgs/Goldstone character, and the breathing doublet) are obtained by solving the extended Gross-Pitaevskii equations with the stated two-component LHY functional, linearizing them into Bogoliubov-de Gennes equations (Eqs. (5)-(8)), and diagonalizing for a fixed set of physical parameters (masses, dipole moments, scattering lengths, trap frequencies, atom number). No parameter is fitted to the computed spectra; the superfluid-mode doublet is a falsifiable output of the model rather than an input. The symmetric case (Sec. V) shows by exact decoupling that the x^2 compressional perturbation is orthogonal to the pure spin mode, which makes the appearance of the doublet in the asymmetric case a nontrivial consequence of spin-density hybridization, not a definitional artifact. Citations to prior work, including Refs. [23], [30], and [31], supply the model and interpretative vocabulary (e.g., the catalyzation effect and Leggett bounds), but the quantitative spectra in Figs. 2-6 are computed in this paper from those stated inputs. No self-citation is used to force the predicted mode structure, and no 'prediction' reduces by construction to a fitted quantity. The main caveats are the quantitative validity of the LHY approximation and the single-parameter-set demonstration of generality, which are correctness and scope risks rather than circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities. It uses one modeling assumption, the extended Gross-Pitaevskii approximation with LHY corrections, which is standard in this field. The diagnostic observables Q, P, eta, and lambda are computed from the BdG modes rather than fitted. The main assumptions are the validity of the eGPE model and the representativeness of the one explored parameter set.

assumptions (5)
  • domain assumption Extended Gross-Pitaevskii equations with the two-component Lee-Huang-Yang term describe the dipolar mixture quantitatively (Eqs. (2)-(4)).
    The entire ground-state and spectral analysis rests on this model, which is a one-loop local-density approximation for quantum fluctuations; Ref. [23] is cited as its source.
  • domain assumption Both components have equal mass m (Sec. II).
    The authors state the mass is assumed equal for simplicity; real isotope mixtures may have unequal masses, which could alter the mode hybridization.
  • domain assumption The mixture remains miscible in all calculations performed.
    The spectrum analysis is restricted to the miscible regime; the immiscible double supersolid is left to Ref. [29].
  • ad hoc to paper Phase boundaries can be determined from density contrast thresholds C1 = 0.99 and C2 = 0.1 and 0.99 (Sec. III).
    The iso-contrast lines are used to define the approximate phase transitions; these thresholds are chosen by hand rather than derived.
  • domain assumption The two-droplet supersolid scenario is representative of the general double-supersolid regime.
    Stated in Sec. VI and in the Conclusions as an expectation, but the simulations cover a single parameter set with unequal dipole moments.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Excitation spectrum of a double supersolid in a trapped dipolar Bose mixture." pith.science (2026). https://pith.science/paper/YF3FR7B7

@misc{pith2026241205215,
  author       = {Pith},
  title        = {Pith review of: Excitation spectrum of a double supersolid in a trapped dipolar Bose mixture},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YF3FR7B7}},
  note         = {Machine review of arXiv:2412.05215}
}
read the original abstract

Dipolar Bose-Einstein condensates are excellent platforms for studying supersolidity, characterized by coexisting density modulation and superfluidity. The realization of dipolar mixtures opens intriguing new scenarios, most remarkably the possibility of realizing a double supersolid, composed by two interacting superfluids. We analyze the complex excitation spectrum of a miscible trapped dipolar Bose mixture, showing that it provides key insights about the double supersolid regime. We show that this regime may be readily probed experimentally by monitoring the appearance of a doublet of superfluid compressional modes, linked to the different superfluid character of each component. Additionally, the dipolar supersolid mixture exhibits a non-trivial spin nature of the dipolar rotons, the Higgs excitation, and the low-lying Goldstone modes. Interestingly, the analysis of the lowest-lying modes allows for monitoring the transition of just one of the components into the incoherent droplet regime, whereas the other remains coherent, highlighting their disparate superfluid properties.

Figures

Figures reproduced from arXiv: 2412.05215 by the authors.

Figure 2
Figure 2. FIG. 2. Modes of the symmetric mixture as a function of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 1
Figure 1. FIG. 1: Snake-like enumeration of the sites of the ladder (top) and e FIG. 1: Sna [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 7
Figure 7. FIG. 7. Ratio [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figures from the paper (2 more)
Figure 6
Figure 6. Figure 6: FIG. 6. Strength [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 1
Figure 1. Figure 1: FIG. 1: Snake-like enumeration of the sites of the ladder (top) and e [PITH_FULL_IMAGE:figures/full_fig_p007_1.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

33 extracted references · 14 canonical work pages

  1. [1]

    Boninsegni and N

    M. Boninsegni and N. V. Prokof’ev, Rev. Mod. Phys. 84, 759 (2012)

    2012

  2. [2]

    Kim and M

    E. Kim and M. H. W. Chan, Nature 427, 225 (2004)

    work page 2004

  3. [3]

    Balibar, Nature 464, 176 (2010)

    S. Balibar, Nature 464, 176 (2010)

    work page 2010

  4. [4]

    D. Y. Kim and M. H. W. Chan, Phys. Rev. Lett. 109, 155301 (2012)

    work page 2012

  5. [5]

    J.-R. Li, J. Lee, W. Huang, S. Burchesky, B. Shteynas, F. C ¸ . Top, A. O. Jamison, and W. Ketterle, Nature543, 91 (2017)

    2017

  6. [6]

    K. T. Geier, G. I. Martone, P. Hauke, and S. Stringari, Phys. Rev. Lett. 127, 115301 (2021)

    work page 2021

  7. [7]

    L´ eonard, A

    J. L´ eonard, A. Morales, P. Zupancic, T. Esslinger, and T. Donner, Nature 543, 87 (2017)

    2017

  8. [8]

    Chomaz, I

    L. Chomaz, I. Ferrier-Barbut, F. Ferlaino, B. Laburthe- Tolra, B. L. Lev, and T. Pfau, Reports on Progress in Physics 86, 026401 (2022)

    2022

Show all 33 references
  1. [9]

    Recati and S

    A. Recati and S. Stringari, Nature Reviews Physics 5, 735 (2023)

    2023

  2. [10]

    Tanzi, E

    L. Tanzi, E. Lucioni, F. Fam` a, J. Catani, A. Fioretti, C. Gabbanini, R. N. Bisset, L. Santos, and G. Modugno, Phys. Rev. Lett. 122, 130405 (2019)

    2019

  3. [11]

    B¨ ottcher, J.-N

    F. B¨ ottcher, J.-N. Schmidt, M. Wenzel, J. Hertkorn, M. Guo, T. Langen, and T. Pfau, Phys. Rev. X 9, 011051 (2019)

    2019

  4. [12]

    Chomaz, D

    L. Chomaz, D. Petter, P. Ilzh¨ ofer, G. Natale, A. Traut- mann, C. Politi, G. Durastante, R. M. W. van Bijnen, A. Patscheider, M. Sohmen, M. J. Mark, and F. Fer- laino, Phys. Rev. X 9, 021012 (2019)

    2019

  5. [13]

    D. S. Petrov, Phys. Rev. Lett. 115, 155302 (2015)

    2015

  6. [14]

    S. M. Roccuzzo and F. Ancilotto, Phys. Rev. A 99, 041601 (2019)

    2019

  7. [15]

    M. Guo, F. B¨ ottcher, J. Hertkorn, J.-N. Schmidt, M. Wenzel, H. P. B¨ uchler, T. Langen, and T. Pfau, Nature 574, 386 (2019)

    2019

  8. [16]

    Tanzi, S

    L. Tanzi, S. M. Roccuzzo, E. Lucioni, F. Fam` a, A. Fioretti, C. Gabbanini, G. Modugno, A. Recati, and S. Stringari, Nature 574, 382 (2019)

    2019

  9. [17]

    Natale, R

    G. Natale, R. M. W. van Bijnen, A. Patscheider, D. Pet- ter, M. J. Mark, L. Chomaz, and F. Ferlaino, Phys. Rev. Lett. 123, 050402 (2019)

    2019

  10. [18]

    Hertkorn, F

    J. Hertkorn, F. B¨ ottcher, M. Guo, J. N. Schmidt, T. Lan- gen, H. P. B¨ uchler, and T. Pfau, Phys. Rev. Lett. 123, 193002 (2019)

    2019

  11. [19]

    Watanabe and T

    H. Watanabe and T. c. v. Brauner, Phys. Rev. D 85, 085010 (2012)

    2012

  12. [20]

    Trautmann, P

    A. Trautmann, P. Ilzh¨ ofer, G. Durastante, C. Politi, M. Sohmen, M. J. Mark, and F. Ferlaino, Phys. Rev. Lett. 121, 213601 (2018)

    2018

  13. [21]

    Durastante, C

    G. Durastante, C. Politi, M. Sohmen, P. Ilzh¨ ofer, M. J. Mark, M. A. Norcia, and F. Ferlaino, Phys. Rev. A 102, 033330 (2020)

    2020

  14. [22]

    Politi, A

    C. Politi, A. Trautmann, P. Ilzh¨ ofer, G. Durastante, M. J. Mark, M. Modugno, and F. Ferlaino, Phys. Rev. A 105, 023304 (2022)

    2022

  15. [23]

    R. N. Bisset, L. A. Pe˜ na Ardila, and L. Santos, Phys. Rev. Lett. 126, 025301 (2021)

    2021

  16. [24]

    Bland, E

    T. Bland, E. Poli, L. A. Pe˜ na Ardila, L. Santos, F. Fer- laino, and R. N. Bisset, Phys. Rev. A106, 053322 (2022)

    2022

  17. [25]

    S. Li, U. N. Le, and H. Saito, Phys. Rev. A 105, L061302 (2022)

    2022

  18. [26]

    Arazo, A

    M. Arazo, A. Gallem ´ ı, M. Guilleumas, R. Mayol, and L. Santos, Phys. Rev. Res. 5, 043038 (2023)

    2023

  19. [27]

    Kirkby, T

    W. Kirkby, T. Bland, F. Ferlaino, and R. N. Bisset, SciPost Phys. Core 6, 084 (2023)

    2023

  20. [28]

    Halder, S

    S. Halder, S. Das, and S. Majumder, Phys. Rev. A 107, 063303 (2023)

    2023

  21. [29]

    Kirkby, A.-C

    W. Kirkby, A.-C. Lee, D. Baillie, T. Bland, F. Ferlaino, P. B. Blakie, and R. N. Bisset, Phys. Rev. Lett. 133, 103401 (2024)

    2024

  22. [30]

    Scheiermann, L

    D. Scheiermann, L. A. Pe˜ na Ardila, T. Bland, R. N. Bis- set, and L. Santos, Phys. Rev. A 107, L021302 (2023)

    2023

  23. [31]

    J. C. Smith, D. Baillie, and P. B. Blakie, Phys. Rev. Lett. 126, 025302 (2021)

    2021

  24. [32]

    Ticknor, Phys

    C. Ticknor, Phys. Rev. A 88, 013623 (2013)

    2013

  25. [33]

    A. J. Leggett, Phys. Rev. Lett. 25, 1543 (1970)

    1970

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.