Pith. sign in

REVIEW 3 major objections 5 minor 32 references

Shift of the Bose-Einstein condensation transition in the presence of a second atomic species

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

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

Pith's one-line read The critical temperature for Bose-Einstein condensation is shifted by the presence of a second atomic species, and this paper derives analytic expressions for that shift, predicting a measurable effect in a sodium-potassium mixture.

desk verdict Useful extension of Giorgini's critical-temperature shift to two-species mixtures, but the printed coordinate transformation in Eq. (14) doesn't match Eq. (15) — factor and α inversion — so the derivation needs a fix. read the letter →

arxiv 2602.12880 v2 pith:5CEFH6BM submitted 2026-02-13 cond-mat.quant-gas physics.atom-ph

Pedro M. Gaspar , Vanderlei S. Bagnato , Patricia C. M. Castilho This is my paper · ORCID
classification cond-mat.quant-gasphysics.atom-ph PACS 03.75.Hh05.30.Jp67.85.-d
keywords Bose-Einsteincondensationcriticaltemperatureshiftbosonicmixturesinterspeciesinteractionssodium-potassiummixturethermalcloudThomas-Fermicondensatephasetransitioncontrol
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

Atomic interactions between two bosonic species affect ultracold mixtures, but the specific effect on the critical temperature of the Bose-Einstein transition had only been estimated in symmetric, fully thermal cases. This paper derives general analytical expressions for the critical-temperature shift of a primary species caused by a secondary bosonic species, treating the two regimes where the secondary species is thermal or already condensed. The formulas are built from a mean-field expansion around the ideal-gas critical point and reduce to double sums that depend on the trap-frequency ratio and the secondary species' chemical potential. Applied to a sodium-potassium mixture with realistic parameters, the shift is predicted to be comparable in magnitude to the known intraspecies interaction shift, and thus observable in current experiments. If the predictions hold, they turn the secondary species into a quantitative control knob for the phase transition, enabling transitions driven by changing atom numbers or interaction strengths.

What carries the argument

The machinery is a first-order perturbation of the coupled Gross-Pitaevskii equations: the interacting density n1(r) is expanded in terms of the ideal-gas densities of both species (Eq. 10), with the interspecies term carrying half the exchange factor of the intraspecies term. The key step is a coordinate rescaling (Eq. 14) that makes the two trapping potentials proportional when the trap-frequency ratios match, turning the spatial integral in Eq. (11) into a sum over polylogarithm-like terms. For the condensed secondary species, the Thomas-Fermi inverted parabola is inserted as the condensed density, giving the additional integral in Eq. (18).

What would settle it

Compare the predicted shift, Eqs. (15) and (18), against a self-consistent Hartree-Fock calculation that allows species-2's density to be modified by species-1's mean-field, or against a direct measurement of the sodium-potassium critical temperature as a function of potassium atom number at fixed sodium number and temperature; a deviation in the shift curve larger than the finite-size correction would indicate the density ansatz or the frequency-proportionality assumption is violated.

Watch

Extended reading notes

Core claim

The paper's central claim is that the first-order correction to the critical temperature of species-1 due to species-2 is given by Eq. (11), and that when species-2 is thermal this reduces to Eq. (15), a double sum over Bose factors; when species-2 is condensed, the Thomas-Fermi part adds the second term in Eq. (18). The authors show that for a potassium-sodium mixture the resulting shift can be made equal in size to the intraspecies shift by choosing the atom-number ratio, and that crossing the secondary species' critical point introduces a distinct feature in the shift curve.

Load-bearing premise

The load-bearing premise is that the secondary species keeps the ideal-gas (or ideal-gas-plus-Thomas-Fermi) density profile of an isolated cloud, unaffected by the primary species' mean field, and that the two trapping potentials are strictly proportional; if the secondary cloud is visibly deformed by the primary species, or if the trap frequencies are not proportional, the quantitative shift predictions will not hold.

Editorial extensions

If this is right

  • The critical temperature of a species becomes a controllable function of the secondary species' atom number, so a transition can be induced by adding or removing atoms of the second species at fixed temperature.
  • The analytic formulas are written for arbitrary conservative traps and any bosonic mixture, so they can be reused for new species pairs once scattering lengths and trap frequencies are known.
  • Because the shift depends on g12, tuning the interspecies scattering length (e.g., via a Feshbach resonance) changes the critical temperature continuously, offering a second control parameter beyond atom number.
  • When the secondary species is condensed, the Thomas-Fermi contribution produces a qualitatively new term, so the shift curve has a signature at the secondary species' critical point that can be looked for experimentally.

Reading between the lines

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

  • If the predicted shift is measured, the temperature at which the primary species condenses could serve as a high-precision readout of the interspecies mean-field potential, effectively turning the BEC transition into a probe of interspecies interactions.
  • The single-iteration density ansatz likely underestimates the shift in strongly interacting or near-miscibility regimes; a self-consistent version that updates species-2's density in the field of species-1 would test the robustness of the quantitative claims while preserving the predicted qualitative control.
  • The same first-order derivation could be translated to Bose-Fermi mixtures by replacing the Bose-Einstein functions with Fermi functions, yielding analogous closed-form shifts for fermionic impurities near the BEC transition.
  • The missing citation for the proportional-trap-frequency assumption suggests this geometric condition may fail in some experimental trap geometries; in those cases the spherical reduction in Eq. (14) would need to be replaced by a full anisotropic integration, which could change the numeric prefactor but not the overall structure of the shift.
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 / 5 minor

Summary. The paper derives first-order mean-field expressions for the shift of the BEC critical temperature of one bosonic species due to interspecies interactions with a second bosonic species, treating the second species either as a thermal gas (Eq. 15) or as a partially condensed cloud (Eq. 18). The calculation follows the established single-species method of Giorgini et al. and is applied to a 23Na-39K mixture, with the claim that the interspecies shift can be comparable to the intraspecies one and hence measurable. A consistency check for identical species is presented in Fig. 1, where the interspecies shift reaches half the intraspecies value at N2=N1.

Significance. If the derivation holds, this is a useful analytic extension of a standard result to Bose-Bose mixtures, with no fitted parameters and a clean limiting check. The formulas are simple enough to be used by experimental groups and the paper explicitly identifies a realistic Na-K system. However, the result is a first-order perturbative mean-field estimate that assumes species-2 retains its isolated-cloud density profile; the quantitative predictions in Fig. 2 should be read with that caveat. Overall, the paper makes a modest but real contribution, provided the technical inconsistencies below are corrected.

major comments (3)
  1. [Eq. (17), §III.B] The Thomas-Fermi radius is defined as R_{i,2}=sqrt(μ2/(m2 ω_{i,2}^2)). For the standard TF density n=(μ-V)/g, the correct relation is R_i^2=2μ/(mω_i^2). The printed definition is a factor √2 too small. This is not merely cosmetic: when Eq. (17) is used in Eq. (11), the printed R_i gives an exponent n α μ2/(2 k_B T) in the condensed contribution, whereas Eq. (18) has n α μ2/(k_B T), which corresponds to the standard 2μ convention. Please correct Eq. (17) (R_i=sqrt(2μ2/(m2 ω_{i,2}^2))) and confirm that Eq. (18) was obtained with this corrected radius.
  2. [Eq. (14) and text following it, §III.A] The transformation in Eq. (14) as printed gives V1/(k_B T)=r'^2/2 and V2/(k_B T)=r'^2/(2α), not V1=α k_B T r'^2 and V2=k_B T r'^2 as stated. The stated simplification requires a different scaling. I have independently checked that Eq. (15) is nevertheless the correct result of applying Eq. (14): the 2^{3/2} from the angular integral combines with (λ_T2)^{-3}, and the resulting α^{-3/2} matches the explicit α in Eq. (15). Thus the concrete prefactor concern raised in the stress-test note does not land; the actual defect is the misleading sentence about the simplified potentials. Please correct the wording so a reader can follow the derivation.
  3. [§III.A, condition after Eq. (14)] The reduction to spherical coordinates requires ω_{x,1}/ω_{x,2}=ω_{y,1}/ω_{y,2}=ω_{z,1}/ω_{z,2}. This condition is stated with a missing citation ('[REF]') and is not automatically satisfied by arbitrary optical or magnetic traps. If it fails, V2 becomes anisotropic in the scaled coordinates and Eq. (15) is no longer the correct integral reduction. Since the abstract claims extension to 'arbitrary conservative traps,' the authors should either provide the generalized expression (possibly with a sum over the three axes) or explicitly restrict the claimed applicability. As written, the scope is narrower than advertised.
minor comments (5)
  1. [After Eq. (14)] The missing '[REF]' citation should be supplied; the claim that the proportional-frequency condition 'is usually the case' needs support.
  2. [§III.B, before Eq. (18)] The text says 'when T^0_{c,2} < T^0_{c,1}' but the condensed case under discussion is T^0_{c,2} > T^0_{c,1}. This is presumably a typo and should be corrected.
  3. [Abstract] Typo: 'cases in with the second species' should be 'cases in which the second species'.
  4. [Fig. 2 caption] Typo: 'in therms of' should be 'in terms of'.
  5. [After Eq. (16)] The phrase 'negative and much larger than the thermal energy' is ambiguous; it should read 'negative with magnitude much larger than k_B T' or equivalent.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the shift formulas are first-order perturbation results with the secondary-species density as an independent input.

full rationale

The core derivation is self-contained. The interspecies shift (δTc,1/T0c,1)12 in Eq. (11) is a first-order perturbative expression in the interspecies coupling g12, and the secondary species enters only through its unperturbed density profile: the ideal trapped-gas profile (12) in the thermal case or the ideal thermal-plus-Thomas-Fermi profile (17) in the condensed case. These profiles are inputs, not quantities fitted to the shift being predicted; no parameter of Eq. (15) or Eq. (18) is calibrated against the output. The check against the identical-species limit uses the external single-species result Eq. (3) as a benchmark, and the half-factor relation is a consequence of the stated factor-of-two difference between intra- and interspecies mean-field terms, not an imposed fit. The standard results cited for the ideal-gas critical temperature, the single-species interaction shift, and the Thomas-Fermi profile are externally established and not unique to this paper; the self-citation [16] is used only for the textbook ideal-gas formula and is not load-bearing. The coordinate-transformation issue in Eq. (14) — where the printed rescaling does not produce the claimed simplified potentials — is a technical/correctness defect, not a circular-reasoning defect: the derived quantity is not equivalent to an input by construction. No fitted-input-called-prediction, self-definitional reduction, or load-bearing self-citation chain is present.

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

The central claim rests on the standard mean-field/HF treatment of trapped Bose gases, on the assumption that the secondary species keeps its isolated-cloud density profile, and on a trap-proportionality condition. No new particles, forces, or conserved quantities are introduced, and no parameters are fitted to the target result.

assumptions (5)
  • domain assumption Hartree-Fock / local-density-approximation density correction n1 = n0_1 - 2 g1 n0_1 ∂n0_1/∂µ1 - g12 n0_2 ∂n0_1/∂µ1 (Eq. 10).
    The entire derivation inherits the single-species HF/LDA treatment from ref. [17]; if this approximation fails near the critical point, the predicted shift is wrong.
  • domain assumption Species-2 density is the unperturbed isolated-cloud density: ideal gas for the thermal case (Eq. 12) and ideal thermal plus Thomas-Fermi for the condensed case (Eq. 17).
    Back-action of species-1 on the species-2 profile and condensate mean-field depletion of the thermal fraction are neglected. This directly sets the overlap integral in Eq. (11).
  • domain assumption Trapping potentials of the two species are proportional, so that after rescaling V1/(kBT)=αr'^2 and V2/(kBT)=r'^2 (Eq. 14).
    This condition is load-bearing for the spherical-coordinate reduction to Eqs. (15) and (18); the paper states it is 'usually the case' but gives a missing [REF].
  • domain assumption Ideal-gas relations for the critical atom number N_c = ζ(3)(kBT/(ℏω))^3 and condensate fraction N_BEC/N = 1 - (T/T_c)^3.
    Used in Eqs. (16) and (17) to fix µ2 and the condensed fraction; interaction corrections to these relations are ignored.
  • domain assumption First-order expansion around T_c = T_c0 + δT_c and µ = 0.
    Truncation at first order in the interaction and in the temperature shift requires the shift to be small; the paper's own numbers suggest a few percent, which is acceptable.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Shift of the Bose-Einstein condensation transition in the presence of a second atomic species." pith.science (2026). https://pith.science/paper/5CEFH6BM

@misc{pith2026260212880,
  author       = {Pith},
  title        = {Pith review of: Shift of the Bose-Einstein condensation transition in the presence of a second atomic species},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5CEFH6BM}},
  note         = {Machine review of arXiv:2602.12880}
}
abstract

Atomic interactions play an important role in the properties of ultracold atomic gases. In single component bosonic systems, its effect is already present at the critical point for the Bose-Einstein condensate phase transition by shifting it to lower temperatures as a consequence of effective repulsion between the atoms. When considering atomic bosonic mixtures, interesting effects arise from the competition between intra- and interspecies interactions such as the miscible-immiscible phase transition and the particular case of self-bounded quantum droplets. In such a scenario, it is natural to expect that these interactions will also affect the critical point of each species composing the mixture. In this paper, we obtain analytical expressions for the critical temperature shift of the phase transition to a Bose-Einstein condensate in the presence of a second species. We treat differently the cases in with the second species is above or below its own critical temperature and apply the obtained relations to the case of a $^{23}$Na-$^{39}$K bosonic mixture which can be realized in current running experimental setups. Our findings can be easily extended to other atomic mixtures trapped by arbitrary conservative traps.

Figures

Figures reproduced from arXiv: 2602.12880 by the authors.

Figure 1
Figure 1. FIG. 1. Comparing the shift due to single species and due [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Shift in the critical temperature of species-1 in therms [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 1 linked inside Pith

  1. [1]

    Observation of bose- einstein condensation in a dilute atomic vapor

    Mike H Anderson, Jason R Ensher, Michael R Matthews, Carl E Wieman, and Eric A Cornell. Observation of bose- einstein condensation in a dilute atomic vapor. science, 269(5221):198–201, 1995

    1995

  2. [2]

    Bose-einstein condensation in a gas of sodium atoms

    Kendall B Davis, M-O Mewes, Michael R Andrews, Nico- laas J van Druten, Dallin S Durfee, DM Kurn, and Wolf- gang Ketterle. Bose-einstein condensation in a gas of sodium atoms. Physical review letters, 75(22):3969, 1995

    1995

  3. [3]

    the shift could be controlled by the strength and sign of the interspecies interaction IV

    The dashed line represents the BEC transition point for species-2. the shift could be controlled by the strength and sign of the interspecies interaction IV. CONCLUSIONS In ultracold gases, understanding how the critical tem- perature varies with the interaction strength, whether in single-species systems or mixtures, is extremely use- ful for designing e...

    arXiv 2013

  4. [4]

    Observation of quantum hawking radi- ation and its entanglement in an analogue black hole

    Jeff Steinhauer. Observation of quantum hawking radi- ation and its entanglement in an analogue black hole. Nature Physics, 12(10):959–965, 2016

    2016

  5. [5]

    Evidence of bose-einstein condensation in an atomic gas with attractive interactions

    Cl C Bradley, CA Sackett, JJ Tollett, and Randall G Hulet. Evidence of bose-einstein condensation in an atomic gas with attractive interactions. Physical review letters, 75(9):1687, 1995

    1995

  6. [6]

    Bose-Einstein condensation and superfluidity, volume 164

    Lev Pitaevskii and Sandro Stringari. Bose-Einstein condensation and superfluidity, volume 164. Oxford Uni- versity Press, 2016

    2016

  7. [7]

    Quantum simu- lations with ultracold atoms in optical lattices

    Christian Gross and Immanuel Bloch. Quantum simu- lations with ultracold atoms in optical lattices. Science, 357(6355):995–1001, 2017

    2017

  8. [8]

    Dy- namics of single neutral impurity atoms immersed in an ultracold gas

    Nicolas Spethmann, Farina Kindermann, Shincy John, Claudia Weber, Dieter Meschede, and Artur Widera. Dy- namics of single neutral impurity atoms immersed in an ultracold gas. Physical review letters, 109(23):235301, 2012

    2012

Show all 32 references
  1. [9]

    Self-trapping of im- purities in bose-einstein condensates: Strong attractive and repulsive coupling

    M Bruderer, W Bao, and D Jaksch. Self-trapping of im- purities in bose-einstein condensates: Strong attractive and repulsive coupling. Europhysics Letters, 82(3):30004, 2008

    2008

  2. [10]

    Binary mixtures of bose condensates of alkali atoms

    Tin-Lun Ho and VB Shenoy. Binary mixtures of bose condensates of alkali atoms. Physical review letters, 77(16):3276, 1996

    1996

  3. [11]

    Ultracold dipolar gas of fermionic na 23 k 40 molecules in their absolute ground state

    Jee Woo Park, Sebastian A Will, and Martin W Zwier- lein. Ultracold dipolar gas of fermionic na 23 k 40 molecules in their absolute ground state. Physical review letters, 114(20):205302, 2015

    2015

  4. [12]

    Properties of two-species bose condensates

    H Pu and NP Bigelow. Properties of two-species bose condensates. Physical review letters, 80(6):1130, 1998

    1998

  5. [13]

    Hartree-fock theory for double conden- sates

    BD Esry, Chris H Greene, James P Burke Jr, and John L Bohn. Hartree-fock theory for double conden- sates. Physical Review Letters, 78(19):3594, 1997

    1997

  6. [14]

    (18) The previous expression can also be solved numerically considering a specific set of parameters defining the Bose- Bose mixture

    Then, the shift in the critical temperature due to the interaction with the second species, whenT 0 c,2 < T0 c,1 is δTc,1 T 0 c,1 12 =− g12(λT 0 c,1,1)−3 3N1kBT 0 c,1 kBT 0 c,1 ℏω2 3 ×   ∞X n,n′=1 1 n1/2n′3/2 1 (αn)3/2 − 1 (αn+n ′)3/2   − 5g12(λT 0 c,1,1)−3 2N1kBT 0 c,1 NB...

  7. [15]

    Topology of the ground state of two interacting bose-einstein conden- sates

    Francesco Riboli and Michele Modugno. Topology of the ground state of two interacting bose-einstein conden- sates. Physical Review A, 65(6):063614, 2002

    2002

  8. [16]

    Phase diagram of two-component bosons on an optical lattice

    Ehud Altman, Walter Hofstetter, Eugene Demler, and Mikhail D Lukin. Phase diagram of two-component bosons on an optical lattice. New Journal of Physics, 5(1):113, 2003

    2003

  9. [17]

    Superfluid-insulator transitions of two-species bosons in an optical lattice

    A Isacsson, Min-Chul Cha, K Sengupta, and SM Girvin. Superfluid-insulator transitions of two-species bosons in an optical lattice. Physical Review B—Condensed Matter and Materials Physics, 72(18):184507, 2005

    2005

  10. [18]

    Bose-einstein condensation in an external po- tential

    Vanderlei Bagnato, David E Pritchard, and Daniel Kleppner. Bose-einstein condensation in an external po- tential. Physical Review A, 35(10):4354, 1987

    1987

  11. [19]

    Con- densate fraction and critical temperature of a trapped interacting bose gas

    Stefano Giorgini, Lev P Pitaevskii, and S Stringari. Con- densate fraction and critical temperature of a trapped interacting bose gas. Physical Review A, 54(6):R4633, 1996

    1996

  12. [20]

    Critical temperature of bose– einstein condensation in trapped atomic bose–fermi mix- tures

    Alexander P Albus, Stefano Giorgini, Fabrizio Illumi- nati, and Luciano Viverit. Critical temperature of bose– einstein condensation in trapped atomic bose–fermi mix- tures. Journal of Physics B: Atomic, Molecular and Optical Physics, 35(23):L511, 2002

    2002

  13. [21]

    Large- nexpansion for condensation and stability of bose- bose mixtures at finite temperatures

    Orest Hryhorchak and Volodymyr Pastukhov. Large- nexpansion for condensation and stability of bose- bose mixtures at finite temperatures. arXiv preprint arXiv:1904.12351, 2019

    1904 arXiv

  14. [22]

    A compact experimental machine for study- ing tunable bose–bose superfluid mixtures

    Patr ´ ıcia Christina Marques Castilho, E Pedrozo-Pe˜ nafiel, EM Gutierrez, PL Mazo, G Roati, KM Farias, and VS Bagnato. A compact experimental machine for study- ing tunable bose–bose superfluid mixtures. Laser Physics Letters, 16(3):035501, 2019

    2019

  15. [23]

    Miscibility regimes in a 23na–39k quantum mixture

    Emmanuel Mercado Gutierrez, Gustavo Alves de Oliveira, Kilvia Mayre Farias, Vanderlei Salvador Bagnato, and Patricia Christina Marques Castilho. Miscibility regimes in a 23na–39k quantum mixture. Applied Sciences, 11(19):9099, 2021

    2021

  16. [24]

    Feshbach spectroscopy and dual-species bose-einstein condensation of na 23- k 39 mixtures

    Torben A Schulze, Torsten Hartmann, Kai K Voges, Matthias W Gempel, Eberhard Tiemann, Alessandro Zenesini, and Silke Ospelkaus. Feshbach spectroscopy and dual-species bose-einstein condensation of na 23- k 39 mixtures. Physical Review A, 97(2):023623, 2018

    2018

  17. [25]

    Bose–Einstein condensation in dilute gases

    Christopher J Pethick and Henrik Smith. Bose–Einstein condensation in dilute gases. Cambridge university press, 2008

    2008

  18. [26]

    Although the chemical potential is zero for the thermal fraction, it is nonzero and positive for the condensed one

  19. [27]

    New experimental system to study coupled vortices in a two-species Bose-Einstein condensate 23Na-41K with tunable interactions

    Patricia Christina Marques Castilho. New experimental system to study coupled vortices in a two-species Bose-Einstein condensate 23Na-41K with tunable interactions. PhD thesis, Universidade de S˜ ao Paulo, 2017

    2017

  20. [28]

    Binary alloy phase diagrams

    Hiroaki Okamoto, TB Massalski, et al. Binary alloy phase diagrams. ASM International, Materials Park, OH, USA, 12:3528–3531, 1990

    1990

  21. [29]

    Quantum phase transition from a superfluid to a mott insulator in a gas of ultracold atoms

    Markus Greiner, Olaf Mandel, Tilman Esslinger, Theodor W H¨ ansch, and Immanuel Bloch. Quantum phase transition from a superfluid to a mott insulator in a gas of ultracold atoms. nature, 415(6867):39–44, 2002

    2002

  22. [30]

    Ob- servation of quantum droplets in a heteronuclear bosonic mixture

    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. Physical Review Research, 1(3):033155, 2019

    2019

  23. [31]

    Quantum liquid droplets in a mixture of bose-einstein condensates

    CR Cabrera, L Tanzi, J Sanz, B Naylor, P Thomas, P Cheiney, and Leticia Tarruell. Quantum liquid droplets in a mixture of bose-einstein condensates. Science, 6 359(6373):301–304, 2018

    2018

  24. [32]

    Self-bound quantum droplets of atomic mixtures in free space

    G Semeghini, G Ferioli, Leonardo Masi, C Mazzinghi, Louise Wolswijk, F Minardi, M Modugno, G Modugno, M Inguscio, and M Fattori. Self-bound quantum droplets of atomic mixtures in free space. Physical review letters, 120(23):235301, 2018

    2018

Pith tools

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