Pith. sign in

REVIEW 3 major objections 6 minor 18 references

Vortices in D-dimensional anisotropic Bose-Einstein condensates: dimensional perturbation theory with hypercylindrical symmetry

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

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

Pith's one-line read The paper derives semiclassical vortex formulas (density, chemical potential, energy, and critical rotation speed) for condensates in any number of dimensions from a single large-dimension limit of the Gross-Pitaevskii equation.

desk verdict A careful DPT extension to D-dimensional vortices that is honest about its limits, but the headline dimensional crossings rest on two uncontrolled expansions and may be artifacts. read the letter →

arxiv 2411.18355 v2 pith:NAANJ5VM submitted 2024-11-27 cond-mat.quant-gas math-phmath.MP

Maria Isabelle Fite , B. A. McKinney This is my paper · ORCID
classification cond-mat.quant-gasmath-phmath.MP PACS 03.75.Lm
keywords Bose-EinsteincondensatevortexGross-PitaevskiiequationdimensionalperturbationtheoryhypercylindricalsymmetryanisotropictrapsyntheticdimensionsThomas-Fermiapproximation
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 aims to show that the limit $\delta = 1/(D+2|m|-d) \to 0$ of the dimensionally scaled Gross-Pitaevskii equation — the mean-field equation for a dilute Bose-Einstein condensate — turns the vortex problem in a $D$-dimensional anisotropic trap into a solvable algebraic form. In that limit the condensate density is an inverted parabola with a hole at the vortex core, $$|\bar\psi(\bar r)|^2 = \frac{1}{\bar u_D}\left(2\bar\mu - \bar r_\$perp^{2}$ - \$lambda^{2}$\bar $z^{2}$ - \frac{\$kappa^{2}$}{4\bar r_\$perp^{2}$}\right),$$ and normalizing this density yields the paper's semiclassical approximations for the chemical potential, the energy per particle, and the critical rotation speed in arbitrary dimension $D$. The authors also observe that the chemical-potential curves for $D = 2,3,4,5$ cross one another as the interaction strength or the trap anisotropy $\lambda$ grows, so at strong parameters a higher-dimensional condensate can become chemically cheaper than a lower-dimensional one. The payoff would be an analytic footing for trapped vortices in arbitrary dimension, with applications the authors point to in synthetic dimensions — where internal states simulate extra spatial dimensions — and in holographic-style correspondences between a higher-dimensional system and a lower-dimensional boundary.

What carries the argument

The engine is the zeroth-order ($\delta \to 0$) limit of the dimensionally scaled hypercylindrical Gross-Pitaevskii equation (Eq. 5), with perturbation parameter $\delta = 1/\kappa$ and $\kappa = D + 2|m| - d$; here $D$ is the spatial dimension, $|m|$ the vortex quantum number, and $d$ a free 'reference dimension' that decides how much kinetic energy survives at zeroth order. A Jacobian transformation $\phi = r_\perp^{(D-2)/2}\psi$ removes the first-derivative terms from the Laplacian, and dimensionally scaled coordinates push the remaining kinetic terms into a prefactor $\delta^2$, so the $\delta\to0$ equation is algebraic and the density follows directly (Eq. 7). A second small parameter, $\bar\alpha = \kappa/(4\bar\mu)$ — the vortex-core radius scale relative to the condensate scale — brings the normalization, energy, and critical-velocity integrals to closed form in gamma functions (Eqs. 16, 24, 27), with the condition $\bar\mu \ge \kappa/2$ keeping the inner and outer radii real. The parameter $d$ connects the general formulas to known limits: $d = 3$ gives exact agreement with the $D = 3$ anisotropic results, and $d = D + 2|m|$ gives the Thomas-Fermi (zero-kinetic-energy) approximation.

What would settle it

Solve the full hypercylindrical Gross-Pitaevskii equation (Eq. 1) numerically for $D = 2,3,4,5$ with $m = 1$, $N = 1000$, and $\lambda = 1$, scanning the interaction strength $\bar u_D$, and compare the chemical-potential curves with Eq. (16) at $d = 3$: if the curves do not cross, or cross at noticeably different $\bar u_D$ values, the zeroth-order claim fails. A cheaper check targets the geometry: compare the numerically computed vortex-core radius as a function of $\bar u_D$ with $r'_{\perp,\min} = \sqrt{(1 - \sqrt{1 - 4\bar\alpha^2})/2}$; at $D = 2$, where $\delta = 1$, the core radius should disagree most strongly if the approximation breaks down exactly where it is applied.

Watch

Extended reading notes

Core claim

Working in hypercylindrical coordinates — one Cartesian axis $z$ plus a $(D-1)$-dimensional hyperspherical subspace with hyperradius $r_\perp$ — the paper writes the Gross-Pitaevskii equation for a condensate in $D$ dimensions carrying a vortex of quantum number $|m|$ along $z$, and applies dimensional perturbation theory with $\delta = 1/\kappa$, where $\kappa = D + 2|m| - d$. At zeroth order ($\delta \to 0$) the derivative kinetic terms drop out and the density is read off directly: $|\bar\psi(\bar r)|^2 = (1/\bar u_D)(2\bar\mu - \bar r_\perp^2 - \lambda^2\bar z^2 - \kappa^2/(4\bar r_\perp^2))$ (Eq. 7), a truncated inverted parabola whose inner edge is the vortex core and whose outer edge is the condensate surface. Imposing the general-$D$ normalization condition and expanding in the second small parameter $\bar\alpha = \kappa/(4\bar\mu)$ yields the chemical potential (Eq. 16), the energy per particle (Eq. 24), and the critical vortex rotation speed (Eqs. 26–27) as explicit functions of $D$, $|m|$, the interaction $\bar u_D$, and the anisotropy $\lambda$. Setting $d = 3$ makes the formulas reduce exactly to the established $D = 3$ axially symmetric results, and setting $d = D + 2|m|$ recovers the general-$D$ Thomas-Fermi approximation (zero kinetic energy). With these formulas, the paper reports crossings of the $\bar\mu_D$-curves for $D = 2,3,4,5$: at weak interaction or anisotropy the chemical potential increases with $D$, as in a harmonic oscillator, while at strong parameters it decreases with $D$ through the Thomas-Fermi exponent $2/(D+2)$.

Load-bearing premise

The central assumption is that dropping the kinetic-energy terms because $\delta = 1/(D+2|m|-d)$ is small remains accurate at the parameters plotted, where for $D = 2,3,4,5$ with $m = 1$ and $d = 3$ this 'small' parameter is actually $1$, $1/2$, $1/3$, and $1/4$.

Editorial extensions

If this is right

  • For any $D \ge 2$, the vortex density, chemical potential, energy per particle, and critical rotation speed can be evaluated from the analytic formulas (Eqs. 7, 16, 24, 26) instead of solving the full Gross-Pitaevskii equation, with all dimensionality entering through $\kappa$ and $\bar\alpha$.
  • The free parameter $d$ interpolates between known limits: $d = 3$ reduces every general-$D$ formula to the established $D = 3$ axially symmetric results, and $d = D + 2|m|$ reduces the chemical potential to the general-$D$ Thomas-Fermi approximation.
  • Reality of the vortex core imposes a $D$- and $|m|$-dependent floor on the chemical potential, $\bar\mu \ge \kappa/2$, which grows with dimension and vorticity and acts like a dimensional zero-point energy.
  • In the zeroth-order density, the vortex-core radius decreases with interaction strength and the condensate density becomes more squeezed around the vortex axis as $D$ increases, giving quantitative predictions for the density profile (Figs. 1 and 2).
  • Within the zeroth-order approximation, raising the interaction or the anisotropy makes the chemical-potential curves for $D = 2,3,4,5$ cross, so that higher-dimensional condensates can have lower chemical potential than lower-dimensional ones; the paper notes that whether the crossings are physical or an artifact of the perturbation limit requires further study.

Reading between the lines

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

  • A testable consequence of the crossings, should they survive beyond zeroth order, is a non-monotonic ordering of vortex-nucleation thresholds: in a synthetic-dimension experiment with tunable interaction, the critical rotation speed for creating an $m=1$ vortex should no longer be monotonic in the effective dimension, and the ordering of thresholds across dimensions should invert as the interactio
  • Because $d$ acts as a variational dial that partially restores the kinetic energy discarded by $\delta \to 0$, an extension the paper leaves open is to optimize $d$ against direct numerical Gross-Pitaevskii solutions for each $(\lambda, \bar u_D, D)$; if the best-fit $d$ tracks $\lambda$ in the way Eqs. (35)–(37) suggest, the crossings in Figs. 3 and 4 may shift or vanish when computed with optimi
  • The pole structure of the $\bar\alpha$ expansion — no poles at second order, poles at odd $D$ appearing from fourth order upward — makes the second-order formulas the only fixed-order expansion valid at every odd dimension; for work at $D = 3$ one must stay at second order or use the exact hypergeometric normalization of Appendix A, since the fourth-order formula has a pole there.
  • The approximation is on its firmest ground when $\delta$ is genuinely small, i.e., for large $D$ or large $|m|$; checking the $D = 3$ formulas at $m = 2$ or $m = 3$, where $\delta$ is halved or thirded, would separate the large-$\kappa$ regime from the small-$D$ regime and indicate how far the 'arbitrary $D$' claim reaches.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 6 minor

Summary. The paper applies dimensional perturbation theory (DPT) to the D-dimensional Gross-Pitaevskii equation in hypercylindrical coordinates with a vortex along the z-axis. With the perturbation parameter δ = 1/(D + 2|m| − d), the authors derive zeroth-order (δ→0) formulas for the condensate density (Eq. 7), chemical potential (Eq. 16), energy per particle (Eq. 24), and critical vortex velocity (Eq. 27). They study the dependence of these quantities on dimension D, interaction strength, and anisotropy λ, and report crossings of the chemical-potential curves for D = 2, 3, 4, 5 as interaction or anisotropy is increased. They also present a Thomas-Fermi limit, a D = 3 comparison to earlier work, and an exact hypergeometric normalization condition (Eq. A.1).

Significance. If the approximations are quantitatively reliable, the paper would supply the first general-D analytic description of vortices in an anisotropic trap, with potential applications to synthetic dimensions and holographic analogies. The derivation is algebraically careful, the D = 3 limit exactly recovers Ref. [9] when d = 3, and the Thomas-Fermi limit (Eq. 20) reduces to known results. The paper also includes a useful exact hypergeometric form of the normalization integral (Eq. A.1) that avoids the α-expansion used elsewhere. In addition, the authors candidly state that the dimensional crossings may be a mathematical artifact (Section 7), which is an honest limitation rather than an overclaim. However, the central quantitative results for D ≠ 3 rest on two uncontrolled approximations—the truncation at zeroth order in δ and at second order in α—so the significance is contingent on validation that the manuscript does not provide.

major comments (3)
  1. [Section 3, Eq. (16) and Figs. 3–4] The chemical potential curves and their dimensional crossings are obtained by expanding the normalization integrand to second order in α = κ/(4μ). For the parameters plotted (d = 3, m = 1, D = 2, 3, 4, 5), the reality condition μ ≥ κ/2 from Eqs. (14)–(15) implies α ≤ 1/2, and at weak interaction or weak anisotropy α is O(0.3–0.5), so the expansion parameter is not small anywhere in the figures. The crossings in Figs. 3 and 4 occur precisely in the weak-to-moderate regime where α is largest, so they may be artifacts of dropping O(α^3) terms. The paper contains an exact normalization condition (Eq. A.1) that avoids this expansion, but it is not used in any figure. Please either provide error estimates on the second-order truncation, implement Eq. A.1 to verify the μ-D curves, or benchmark against numerical GPE solutions for D = 2, 4, 5.
  2. [Section 3, Eq. (7) and definition of δ] The paper treats the δ→0 limit as the zeroth-order approximation and applies it to finite D with d = 3 and m = 1, where δ = 1/(D + 2|m| − d) takes values 1, 1/2, 1/3, 1/4 for D = 2, 3, 4, 5. For D = 2 and D = 3 the perturbation parameter is not small, so the validity of dropping the order-δ^2 kinetic derivative terms in Eq. (5) to obtain Eq. (7) is not established. The text frames the method as a large-D or large-|m| approximation, but the figures include parameter regimes where neither is large. Please state the expected domain of validity of the zeroth-order approximation, or include a comparison with finite-δ corrections for at least one of the plotted dimensions.
  3. [Section 6, Fig. 6] The free parameter d is shown to have a strong effect on the chemical potential: for λ = 1/100 the optimal d is between 1 and 2, while for λ = 100 the optimal d is between −68 and −90. Since d enters κ = D + 2|m| − d and hence the density and chemical potential, the dimensional crossings in Figs. 3 and 4, which are computed for d = 3, may not be robust to the choice of d. The paper should demonstrate that the qualitative crossings persist over a range of d, or should determine d from a variational or other principle rather than using a fixed value that is later tuned.
minor comments (6)
  1. [Section 2] The text says the hypersphere reduces to a circle and line segment segment; the duplicated word 'segment' should be removed.
  2. [Figure 3 caption] The label '2D3D4D5D' appears to be a formatting artifact; it should be spaced or replaced by a proper legend.
  3. [Section 6] The word 'anistropy' in the caption of Fig. 6 is a typo and should be 'anisotropy.'
  4. [References] Reference [18] lists the journal as 'A VS Quantum Science'; the correct title is 'AVS Quantum Science.'
  5. [Section 3, Eq. (16)] The derivation of Eq. (16) from Eq. (13) by expanding to second order in α is central but not shown in detail; a brief appendix or inline derivation would make the truncation errors easier to assess.
  6. [Section 5, Eq. (33)] The D = 4 critical velocity expression in Eq. (33) is not obviously the D = 4 limit of Eq. (27); a consistency check or a note explaining the simplified form would improve confidence.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the DPT derivation is self-contained and benchmarked externally; self-citations are methodological, not load-bearing.

full rationale

The paper's derivation chain is largely self-contained: the zeroth-order density (Eq. 7) follows by taking the stated delta-to-zero limit in the dimensionally scaled GPE (Eq. 5), and the chemical potential (Eq. 16), energy (Eq. 24), and critical vortex speed (Eq. 27) are obtained by substituting that density into the normalization, energy, and critical-velocity integrals and expanding in alpha = kappa/(4 mu). The exact hypergeometric normalization condition (Eq. A.1) is also supplied, so the alpha expansion is an optional simplification rather than the definition of the result. No parameter is fitted to a quantity that is later relabeled as a prediction: the d values in Fig. 6 are explicitly treated as a variational comparison, and Figs. 1-4 use the fixed reference choice d = 3. The choice d = 3 makes the D = 3 limit coincide with Ref. [9] by construction, but this is an explicit parameter choice for benchmarking, and the D > 3 content is not forced by that choice. Self-citations to Refs. [7,8,15,16] supply the DPT technique and lower-dimensional coupling limits, but the current derivation of the scaled GPE and the zeroth-order equations is exhibited in the text, so these citations are methodological rather than load-bearing. The paper also benchmarks against the external variational results of Ref. [10] and against standard Thomas-Fermi limits. The acknowledged limitations -- delta = 1/(D + 2|m| - d) and alpha are not small over much of the plotted range, and the conclusion states 'More study is needed to determine whether the dimensional crossings are physical or a mathematical artifact, possibly due to the perturbation limit' -- are accuracy concerns, not circularity. No equation reduces to its own input by construction, so the circularity score is 0.

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

The central results rest on the GPE as a mean-field model, the hypercylindrical coordinate/eigenvalue structure, and the validity of a zeroth-order expansion in δ plus a truncated expansion in α. The only hand-tuned quantity is d, which changes the kinetic energy content and is adjusted in §6 to match variational data.

free parameters (1)
  • d = d=3 in Figs 1-4; d=1,2 for λ=1/100 and d=-68,-90 for λ=100 in Fig 6
    d is an arbitrary 'reference dimension' in κ=D+2|m|-d that controls how much centrifugal kinetic energy survives at zeroth order. Different d change μ and the μD-curves; in Fig 6 d is tuned to match the variational results of Ref [10]. The central results therefore depend on a hand-chosen parameter.
assumptions (4)
  • domain assumption The Gross-Pitaevskii mean-field equation accurately describes the dilute BEC in D dimensions with contact interaction u_D.
    Invoked as the starting point in Eq (1) without correction; standard but a modeling assumption.
  • standard math The eigenvalues of L^2_{D-2} are |m|(|m|+D-3), and the hypercylindrical coordinates have a (D-1)-hypersphere angular volume Ω_{D-1}.
    Used to write Eq (4) and the normalization condition Eq (8); relies on the convention in Refs [11,12]. For D=2 the hypersphere is degenerate, so this assumption is strained.
  • ad hoc to paper The zeroth-order δ→0 limit, dropping the δ^2 kinetic derivative terms, provides a valid approximation for the finite D and m values considered.
    δ=1/(D+2|m|-d) is not small for D=2,3 with m=1,d=3; the paper uses the limit anyway to produce Eqs (7), (16), and Figs 1-4.
  • ad hoc to paper The normalization integral can be truncated at second order in α=κ/(4μ), and higher-order poles are ignored.
    Eqs (16) and (24) use the second-order α expansion; higher orders have poles at odd D. No error bound is given.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Vortices in D-dimensional anisotropic Bose-Einstein condensates: dimensional perturbation theory with hypercylindrical symmetry." pith.science (2026). https://pith.science/paper/NAANJ5VM

@misc{pith2026241118355,
  author       = {Pith},
  title        = {Pith review of: Vortices in D-dimensional anisotropic Bose-Einstein condensates: dimensional perturbation theory with hypercylindrical symmetry},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NAANJ5VM}},
  note         = {Machine review of arXiv:2411.18355}
}
abstract

We investigate D-dimensional atomic Bose-Einstein condensates in a hypercylindrical trap with a vortex core along the z-axis and quantized circulation $\hbar m$. We analytically approximate the hypercylindrical Gross-Pitaevskii equation using dimensional perturbation theory with perturbation parameter $\delta=1/(D+2|m|-d)$, \textcolor{black}{where $d$ controls the contribution of kinetic energy at zeroth order}. We derive the zeroth-order ($\delta \to 0$) semiclassical approximations for the condensate energy, density, chemical potential, and critical vortex rotation speed in arbitrary dimensions. We investigate the effect of trap anisotropy on lower effective dimensionality and compute properties of vortices in higher dimensions motivated by the study of synthetic dimensions and holographic duality, where a higher-dimensional gravitational model corresponds to a lower-dimensional quantum model. In the zeroth-order approximation, we observe crossings between energy levels for different dimensions as a function of interaction strength and anisotropy parameters.

Figures

Figures reproduced from arXiv: 2411.18355 by the authors.

Figure 1
Figure 1. FIG. 1: Side view cross section of zeroth-order density [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Inner (vortex core) radius [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Chemical potential ¯µ [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Chemical potential ¯µ [PITH_FULL_IMAGE:figures/full_fig_p013_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: Contour plots of the zeroth-order density [PITH_FULL_IMAGE:figures/full_fig_p017_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p018_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references · 16 canonical work pages

  1. [9]

    Sinha, Phys

    S. Sinha, Phys. Rev. A 55, 4325 (1997)

    work page 1997

  2. [1]

    A. L. Fetter and A. A. Svidzinsky, Journal of Physics: Condensed Matter 13, R135 (2001)

    work page 2001

  3. [2]

    G¨ orlitz, J

    A. G¨ orlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, T. L. Gustavson, J. R. Abo-Shaeer, A. P. Chikkatur, S. Gupta, S. Inouye, T. Rosenband, and W. Ketterle, Phys. Rev. Lett. 87, 130402 (2001)

    work page 2001

  4. [3]

    McCanna and H

    B. McCanna and H. M. Price, Phys. Rev. Res. 3, 023105 (2021)

    work page 2021

  5. [4]

    Boada, A

    O. Boada, A. Celi, J. I. Latorre, and M. Lewenstein, Phys. Rev. Lett. 108, 133001 (2012)

    2012

  6. [5]

    Sugawa, F

    S. Sugawa, F. Salces-Carcoba, A. R. Perry, Y. Yue, and I. B. Spielman, Science 360, 1429 (2018)

    work page 2018

  7. [6]

    Emergence of Large-Scale Structures in Holographic Superfluid Turbulence,

    W.-C. Yang, C.-Y. Xia, Y. Tian, M. Tsubota, and H.-B. Zeng, “Emergence of Large-Scale Structures in Holographic Superfluid Turbulence,” (2024), arXiv:2402.17980 [hep-th]

    arXiv 2024

  8. [7]

    B. A. McKinney, M. Dunn, and D. K. Watson, Phys. Rev. A 69, 053611 (2004)

    work page 2004

Show all 18 references
  1. [8]

    B. A. McKinney and D. K. Watson, Phys. Rev. A 65, 033604 (2002)

    2002

  2. [10]

    K. K. Das, Phys. Rev. A 66, 053612 (2002)

    2002

  3. [11]

    in Dimensional Scaling in Chemical Physics, edited by D. R. Herschbach, J. Avery, and O. Goscinski (Kluwer Academic Publishers, Dordrecht, 1992) pp. 61–80

    1992

  4. [12]

    D. D. Frantz and D. R. Herschbach, The Journal of Chemical Physics 92, 6668 (1990)

    1990

  5. [13]

    D. R. Herschbach, Journal of Chemical Physics 84, 838 (1986)

    1986

  6. [14]

    Baym and C

    G. Baym and C. J. Pethick, Phys. Rev. Lett. 76, 6 (1996)

    1996

  7. [15]

    L. G. McKinney and B. A. McKinney, Physica Scripta 98, 015404 (2023)

    2023

  8. [16]

    T. T. Le, Z. Osman, D. Watson, M. Dunn, and B. A. McKinney, Physica Scripta 94, 065203 24 (2019)

    2019

  9. [17]

    Mumford, D

    J. Mumford, D. Kamp, and D. H. J. O’Dell, Physical Review A 110, 043310 (2024)

    2024

  10. [18]

    Tononi, L

    A. Tononi, L. Salasnich, and A. Yakimenko, A VS Quantum Science 6, 030502 (2024). 25

    2024

Pith tools

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