Pith. sign in

REVIEW 3 major objections 5 minor 41 references

A vortex state of ultralight dark matter suppresses dynamical friction for co-rotating globular clusters, naturally alleviating the Fornax timing problem.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-04 00:37 UTC pith:ZJYXMKEP

load-bearing objection A useful model-level study of vortex ULDM and Fornax globular clusters, but the central 'alleviated' claim outruns the calculation: local timescales are not infall times. the 3 major comments →

arxiv 2608.00258 v1 pith:ZJYXMKEP submitted 2026-07-31 astro-ph.GA

Vortex State of Ultralight Dark Matter and the Fornax Timing Problem

classification astro-ph.GA
keywords ultralight dark matterdynamical frictionglobular clustersFornax timing problemvortex solitonBose-Einstein condensate dark matterorbital decayGross-Pitaevskii-Poisson
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tries to establish that the old puzzle of Fornax's surviving globular clusters can be explained if its dark-matter core is a rotating vortex rather than a static ball. In a vortex state the dark matter flows around the center, so a cluster moving in the same direction sees a much smaller relative velocity; where the cluster speed matches the dark-matter flow, dynamical friction nearly vanishes and the orbital-decay time spikes. Computing these times for three observed Fornax clusters, the paper finds the heaviest one, GC3—the hardest to explain in the standard picture—is the one most rescued by the vortex. If this is right, the survival of these old clusters becomes a probe of the quantum state of dark matter rather than a contradiction needing modified gravity.

Core claim

On the paper's own terms: the dynamical friction acting on a globular cluster in ultralight dark matter is controlled by the relative velocity between the cluster and the local dark-matter flow, not by the cluster's speed alone. In the vortex soliton the dark matter rotates around the halo center, so a co-rotating cluster has |v_GC - u(r)| near zero at certain radii. There the gravitational wake is suppressed and the characteristic time T = v_GC M / F_fr develops sharp peaks, exceeding ten gigayears and reaching far higher values for GC3. The vortex's toroidal density profile also removes dark matter from the innermost region, further weakening the drag. The paper concludes that the Fornax t

What carries the argument

The central object is the s=1 vortex soliton of the Gross-Pitaevskii-Poisson system: a rotating toroidal condensate whose density vanishes on the axis and whose particle current defines a circulating velocity field u(r) = (hbar/m r_perp) e_phi. The load-bearing identity is the relative speed v = |v_GC - u(r)| that enters the dynamical-friction force; the force is evaluated pointwise at the cluster position through an angular-momentum expansion of the gravitational wake, controlled by the Mach number v/c_s and the local condensate density. The vortex alters both inputs to the drag—lowering the inner density and, for corotating orbits, reducing the relative velocity—which is what produces the

Load-bearing premise

The dynamical-friction formulas are derived for a homogeneous, stationary ultralight-dark-matter medium but are applied pointwise to the strongly varying vortex density and velocity fields; if gradients on the scale of the gravitational wake matter, the predicted peaks—and the resulting alleviation—are not guaranteed.

What would settle it

A direct test: compute the drag on a point mass moving in the actual vortex density and velocity profile without the local-homogeneity approximation; if the characteristic-time peaks disappear, the alleviation does not hold. Observationally, measure GC3's orbit: the suppression requires the cluster to corotate with the halo flow near the peak radius, so a counter-rotating or strongly eccentric orbit would falsify the proposed mechanism.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Co-rotating globular clusters in a vortex dark-matter halo can survive more than 10 Gyr at radii where the same cluster in the ground-state soliton would spiral inward in a few Gyr.
  • The Fornax timing problem can be eased by the quantum state of the dark-matter core, with no need for baryonic feedback or modified gravity.
  • The orbital direction of a cluster becomes a measurable quantity: counter-rotating clusters feel enhanced drag while corotating ones are protected, so the spread in cluster orbits is itself a dark-matter diagnostic.
  • The vortex's depleted center slows orbital decay for all clusters at small radii, not just the corotating ones, because the local dark-matter density there is much lower than in the ground state.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Editorial extension: the same relative-velocity argument should apply to other massive tracers in rotating dark-matter cores, such as satellite galaxies and stellar bars; the paper points to bar torques as future work but does not quantify them.
  • Editorial extension: the protection is resonant—it occurs only where the cluster speed matches the local flow—so the net effect depends on how long the cluster lingers near that radius; a cluster passing through quickly would not be saved.
  • Editorial extension: a kinematic test follows directly: proper motions of the Fornax clusters should show the innermost cluster corotating with the inferred vortex flow, and counter-rotating clusters should sink faster.
  • Editorial extension: if realistic ultralight-dark-matter halos generically contain vortices, the well-known 'core stalling' of infalling satellites in cored halos may be partly a wave-mechanical effect rather than a purely classical cored-potential effect.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper studies dynamical friction acting on globular clusters embedded in an ultralight dark matter (ULDM) halo described by a rotating vortex soliton, and compares the resulting characteristic timescale for velocity change with that obtained for the non-rotating ground state. Using an analytic dynamical-friction formula from previous work, the authors compute the local timescale T = v_GC M / (F_fr,ULDM + F_fr,B) for three Fornax globular clusters (GC3, GC4, GC6) at their current projected radii. They find that in the corotating case the relative velocity between the cluster and the ULDM flow can vanish at certain radii, producing sharp peaks in T, and they conclude that the Fornax timing problem is 'naturally alleviated.' The paper is clearly written and builds on a substantial body of prior work by the same group and others, but the central claim rests on a local timescale rather than an actual orbital integration, which is acknowledged in the conclusions yet overstated in the abstract.

Significance. If the claimed suppression of dynamical friction in the vortex state were rigorously established, it would offer an interesting mechanism by which the quantum state of the dark matter halo could solve the Fornax timing problem without invoking baryonic feedback or modified gravity. The paper identifies a concrete and physically motivated effect—co-rotation reduces the relative velocity and thereby weakens the gravitational wake—and it provides an explicit comparison between ground-state and vortex-state predictions. The main analytical input (Eqs. 19–22) is drawn from the established literature on ULDM dynamical friction, and the paper is transparent about several limitations, including the need for future orbital-evolution studies. However, as presented, the quantitative claim of alleviation is not yet demonstrated; the leap from a local velocity-change timescale to an infall time is the key gap.

major comments (3)
  1. [Abstract & Sec. 4] The central claim that the Fornax timing problem is 'naturally alleviated' is not supported by the analysis. Eq. (18) defines T as a local velocity-change timescale evaluated at the cluster's current radius, but the survival of a globular cluster for >10 Gyr requires that its orbit does not decay on that timescale. The paper only computes T at present projected positions (Figs. 2–4); a cluster on an inspiraling orbit must first traverse regions where T is short before it can reach a corotation peak. The conclusions acknowledge this gap ('importance of further research on the evolution of globular clusters from a starting point to their current position'), yet the abstract and conclusions nonetheless claim the timing problem is alleviated. An orbital integration, or at least a conservative bound on the infall time that accounts for the radial profile of T, is needed to substantiate the he
  2. [Sec. 3, Eqs. (19)–(22)] The dynamical-friction formula (21)–(22) is derived for a homogeneous, stationary ULDM medium (Ref. [32]). It is applied pointwise to the vortex, with the medium entering only through the local density ρDM(r) and the relative velocity v=|v_GC−u(r)|. No justification is given that the density and velocity gradients of the vortex are small on the scale of the gravitational wake. Near the corotation peaks, where v→0, the wake can become large and the local approximation is especially questionable. The quantitative values of T—and hence the claimed alleviation—therefore depend on an unvalidated assumption. The authors should either justify the local approximation with an estimate of the wake size versus the gradient scale, or confront the result with a spatially resolved calculation.
  3. [Sec. 2, Eq. (4)] The vortex density profile ρ1(r)=ρc r⊥^2 e^{−r^2/R^2}/(r⊥^2+(2ξ)^2) is introduced as an ansatz, not derived as a stationary solution of the GPP system (1)–(2). The velocity field u(r) from Eq. (14) and all subsequent dynamical-friction predictions inherit this ansatz. While Ref. [18] is cited, the paper does not verify that the ansatz satisfies the GPP equations with the adopted baryonic potential, nor whether such a vortex is dynamically stable on the ~10 Gyr timescales of interest. Since the suppression peaks depend sensitively on the velocity field, this should be checked or at least explicitly discussed as a limitation.
minor comments (5)
  1. [Eq. (5)] The notation is inconsistent: 'as' appears in Eq. (5) and in the text, while the subscripted form 'a_s' is used elsewhere. Use one convention throughout. Also the dimensionless argument of the fourth root is not immediately clear and should be parenthesized.
  2. [Figs. 3 and 4 captions] Figure 3 caption says 'v = u vGC', which is garbled; it should be v = v_GC + u for counter-rotation. Figure 4 caption says 'v = u - vGC', but the text uses v = |v_GC − u|; the absolute value should be indicated.
  3. [Sec. 3, after Eq. (18)] The circular velocity v_GC used in Eq. (18) is not defined explicitly. Presumably it is the circular velocity v_1B from Eq. (12), but this should be stated for reproducibility.
  4. [Introduction] Typo: 'a ultralight dark matter' should be 'an ultralight dark matter'.
  5. [References] Reference [23] contains a corrupted author string ('E. L. /suppress Lokas'); this should be corrected.

Circularity Check

0 steps flagged

No significant circularity: the Fornax result is an application of externally sourced dynamical-friction formulas, not a fit or self-referential derivation; minor self-citations are non-load-bearing.

full rationale

Walking the derivation chain: the ULDM ground and vortex states are taken from external variational/BEC literature (Refs. [18,20]), and the Fornax parameters (R, rho_c, M_B, b) come from observational constraints (Refs. [22,23,26]). The dynamical-friction force is explicitly attributed to external linear-response calculations (Refs. [29,30,32]), with the relative-velocity dependence cited to Ref. [31]. The paper's own Refs. [10] and [36] provide prior context and numerical/analytical technique, but the governing formulas are not sourced exclusively to the authors, so the self-citation is not definitionally load-bearing. The corotation peaks in Fig. 4 follow from the stated dependence of the friction force on the relative velocity v=|v_GC-u| and the vortex velocity field; this is a physical consequence of the model, not a fitted parameter renamed as a prediction. No parameter is fitted to the survival of the Fornax clusters. The main gap is that the paper computes the local instantaneous timescale T, not a full orbital-decay integration, and the manuscript itself acknowledges in Sec. 4 that further research on cluster evolution from a starting point to the current position is needed. That is an overclaim/scope issue, not circularity. Overall, the central derivation is self-contained and the Fornax-specific statement is an application rather than an identity with the inputs.

Axiom & Free-Parameter Ledger

2 free parameters · 5 axioms · 0 invented entities

The central claim rests on a specific vortex ansatz and a homogeneous-medium DF formula applied locally. The only scanned physical parameter is the ULDM particle mass m; g is derived from the mass-radius relation. No new particles or forces are invented.

free parameters (2)
  • ULDM particle mass m = Scanned over 3.09×10^-22 to 10^-21 eV
    The paper presents results for five values; for each, the self-interaction g is fixed via Eq. (17). The characteristic times and peak locations depend on m.
  • Self-interaction coupling g=4πℏ²a_s/m = 6.52×10^-50 to 2.1×10^-46 eV^-2 across the m scan
    Determined from the mass-radius relation once m is chosen; not independently measured. Included because the model's dynamical friction depends on it.
axioms (5)
  • domain assumption The Gross-Pitaevskii-Poisson system (Eqs. 1-2) with weak repulsive self-interaction is the correct description of ULDM halos.
    Standard framework for ULDM, justified by large occupation numbers; cited from [7,17].
  • domain assumption The vortex density profile ρ1(r)=ρc r⊥² e^{-r²/R²}/(r⊥²+(2ξ)²) (Eq. 4) represents a stable stationary ULDM halo core.
    Adopted from prior BEC vortex literature [18]; not derived from GPP in this paper, and its applicability to Fornax is assumed.
  • ad hoc to paper The dynamical friction force formulas (Eqs. 19-22) derived for a homogeneous, stationary medium can be applied locally to the inhomogeneous vortex, with the flow entering only through the relative velocity in the Mach number.
    The paper uses local density ρDM(r) and relative velocity |vGC-u(r)| in a homogeneous-medium formula; the validity of this 'local density approximation' is not assessed.
  • domain assumption Globular clusters move on circular orbits in the equatorial plane z=0 with velocity equal to the local circular velocity of the combined potential.
    Stated in Sec. 3; the analysis evaluates T at fixed r, so non-circular or inclined orbits are ignored.
  • domain assumption The baryonic matter follows a Plummer sphere with MB=3×10^7 M⊙ and b=668 pc, and the DF from baryons is Chandrasekhar with σ=10 km/s.
    Standard adopted values from [21-23,33].

pith-pipeline@v1.3.0-alltime-deepseek · 11593 in / 14011 out tokens · 135841 ms · 2026-08-04T00:37:09.611128+00:00 · methodology

0 comments
read the original abstract

We investigate the impact of the vortex state of the ultralight dark matter (ULDM) on the dynamical friction acting on moving globular clusters. Comparing this force with that for the solitonic ground state, it is shown that the internal structure and rotation of the ULDM core strongly affect the orbital decay of globular clusters. In particular, co-directional rotation in a vortex state can lead to significant suppression of dynamic friction at certain distances where globular clusters and ULDM velocities match. Applying these findings to the Fornax dwarf galaxy, it is found that the Fornax timing problem is naturally alleviated.

Figures

Figures reproduced from arXiv: 2608.00258 by A. Zaporozhchenko, E. Gorbar, K. Korshynska, O. Barabash, O. Teslyk, T. Gorkavenko, V. Gorkavenko.

Figure 1
Figure 1. Figure 1: The gravitational potential with and without bary [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The characteristic time when GC changes notably it [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The characteristic time when GC changes notably it [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The characteristic time when GC changes notably it [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 1 canonical work pages

  1. [1]

    A. B. Pace, M. G. Walker, S. E. Koposov, N. Caldwell, M. Mat eo, E. W. Olszewski et al., Spectroscopic confirmation of the sixth globular cluster in the fornax dwarf spheroidal galaxy, Astrophys. J. 923 (2021) 77

  2. [2]

    K. Oh, D. Lin and H. Richer, Globular clusters in the Fornax dwarf spheroidal galaxy , The Astrophysical Journal 531 (2000) 727

  3. [3]

    Meadows, J

    N. Meadows, J. Navarro, I. Santos-Santos, A. Benitez-Ll ambay and C. Frenk, Cusp or core? Revisiting the globular cluster timing problem in For nax, Mon. Not. Roy. Astron. Soc. 491 (2020) 3336

  4. [4]

    N. Bar, D. Blas, K. Blum and H. Kim, Assessing the Fornax globular cluster timing problem in different models of dark matter , Phys. Rev. D 104 (2021) 043021 [2102.11522]

  5. [5]

    D. J. E. Marsh, Axion Cosmology, Phys. Rept. 643 (2016) 1 [1510.07633]

  6. [6]

    Lee, Brief History of Ultra-light Scalar Dark Matter Models , EPJ Web Conf

    J.-W. Lee, Brief History of Ultra-light Scalar Dark Matter Models , EPJ Web Conf. 168 (2018) 06005 [1704.05057]

  7. [7]

    E. G. M. Ferreira, Ultra-light dark matter , Astron. Astrophys. Rev. 29 (2021) 7 [2005.03254]

  8. [8]

    Matos, L

    T. Matos, L. A. Ure˜ na L´ opez and J.-W. Lee,Short review of the main achievements of the scalar field, fuzzy, ultralight, wave, BEC dark matter model , Front. Astron. Space Sci. 11 (2024) 1347518 [2312.00254]

  9. [9]

    Schive, Fuzzy dark matter simulations , Liv

    H.-Y. Schive, Fuzzy dark matter simulations , Liv. Rev. Comput. Astrophys. 12 (2026) 1 [2509.23231]

  10. [10]

    E. V. Gorbar, O. V. Barabash, V. M. Gorkavenko, K. Korshy nska, A. I. Momot and A. O. Zaporozhchenko, Damping of dynamical friction force in self-interacting ul tralight dark matter and fornax timing problem , Class. Quant. Grav. 43 (2026) 155003 [2511.06123]

  11. [11]

    Chandrasekhar, Dynamical Friction

    S. Chandrasekhar, Dynamical Friction. I. General Considerations: The Coeffici ent of Dynamical Friction, The Astrophysical Journal 97 (1943) 255

  12. [12]

    Bondi and F

    H. Bondi and F. Hoyle, On the Mechanism of Accretion by Stars , Mon. Not. Roy. Astron. Soc. 104 (1944) 273

  13. [13]

    Dokuchaev, Emission of Magnetoacoustic Waves in the Motion of Stars in C osmic Space, Soviet Astronomy 8 (1964) 23

    V. Dokuchaev, Emission of Magnetoacoustic Waves in the Motion of Stars in C osmic Space, Soviet Astronomy 8 (1964) 23

  14. [14]

    M. A. Ruderman and E. A. Spiegel, Galactic Wakes, Astrophys. J. 165 (1971) 1

  15. [15]

    Rephaeli and E

    Y. Rephaeli and E. E. Salpeter, Flow Past a Massive Object and the Gravitational Drag , Astrophys. J. 240 (1980) 20

  16. [16]

    E. C. Ostriker, Dynamical Friction in a Gaseous Medium , The Astrophysical Journal 513 (1999) 252

  17. [17]

    Calmet, ed., Quantum aspects of black holes

    X. Calmet, ed., Quantum aspects of black holes . Springer, 2015, 10.1007/978-3-319-10852-0

  18. [18]

    Y. O. Nikolaieva, A. O. Olashyn, Y. I. Kuriatnikov, S. I. Vilchynskii and A. I. Yakimenko, Stable vortex in Bose-Einstein condensate dark matter , Low Temp. Phys. 47 (2021) 684 [2103.07856]. 12

  19. [19]

    Korshynska, Y

    K. Korshynska, Y. M. Bidasyuk, E. V. Gorbar, J. Jia and A. I. Yakimenko, Dynamical galactic effects induced by solitonic vortex structure in bo sonic dark matter , Eur. Phys. J. C 83 (2023) 451 [2301.13110]

  20. [20]

    Chavanis, Mass-radius relation of Newtonian self-gravitating Bose- Einstein condensates with short-range interactions: I

    P.-H. Chavanis, Mass-radius relation of Newtonian self-gravitating Bose- Einstein condensates with short-range interactions: I. Analytical results, Phys. Rev. D 84 (2011) 043531 [1103.2050]

  21. [21]

    H. C. Plummer, On the problem of distribution in globular star clusters , Monthly Notices of the Royal Astronomical Society 71 (1911) 460

  22. [22]

    Cappellari, R

    M. Cappellari, R. Bacon, M. Bureau, M. C. Damen, R. L. Dav ies, P. T. De Zeeuw et al., The SAURON Project–IV. The Mass-to-Light Ratio, the Virial Mass Estimator and the Fundamental Plane of Elliptical and Lenticular Galaxies , Monthly Notices of the Royal Astronomical Society 366 (2006) 1126

  23. [23]

    Kowalczyk, A

    K. Kowalczyk, A. del Pino, E. L. /suppress Lokas and M. Valluri,Schwarzschild dynamical model of the fornax dwarf spheroidal galaxy , Monthly Notices of the Royal Astronomical Society 482 (2019) 5241

  24. [24]

    Chavanis, Mass-radius relation of self-gravitating Bose-Einstein c ondensates with a central black hole , The European Physical Journal Plus 134 (2019) 352

    P.-H. Chavanis, Mass-radius relation of self-gravitating Bose-Einstein c ondensates with a central black hole , The European Physical Journal Plus 134 (2019) 352

  25. [25]

    A. S. Dmitriev, D. G. Levkov, A. G. Panin, E. K. Pushnaya a nd I. I. Tkachev, Instability of rotating Bose stars , Phys. Rev. D 104 (2021) 023504 [2104.00962]

  26. [26]

    J. R. Jardel and K. Gebhardt, The dark matter density profile of the Fornax dwarf , The Astrophysical Journal 746 (2012) 89

  27. [27]

    V. M. Gorkavenko, A. I. Yakimenko, A. O. Zaporozhchenko and E. V. Gorbar, Dynamical friction in ultralight dark matter: Plummer sphere perspec tive, Phys. Scripta 100 (2025) 075039 [2412.15428]

  28. [28]

    O. V. Barabash, T. V. Gorkavenko, V. M. Gorkavenko, O. M. Teslyk, N. S. Yakovenko, A. O. Zaporozhchenko et al., Analytic Calculation of Dynamical Friction for Plummer Sphere in Ultralight Dark Matter , Ukr. J. Phys. 70 (2025) 759 [2504.06448]

  29. [29]

    Desjacques, A

    V. Desjacques, A. Nusser and R. Buehler, Analytic solution to the dynamical friction acting on circularly moving perturbers , The Astrophysical Journal 928 (2022) 64

  30. [30]

    Buehler and V

    R. Buehler and V. Desjacques, Dynamical friction in fuzzy dark matter: Circular orbits , Phys. Rev. D 107 (2023) 023516 [2207.13740]

  31. [31]

    L. Hui, J. P. Ostriker, S. Tremaine and E. Witten, Ultralight scalars as cosmological dark matter, Phys. Rev. D 95 (2017) 043541 [1610.08297]

  32. [32]

    Berezhiani, G

    L. Berezhiani, G. Cintia, V. De Luca and J. Khoury, Dynamical friction in dark matter superfluids: The evolution of black hole binaries , JCAP 06 (2024) 024 [2311.07672]

  33. [33]

    M. G. Walker, M. Mateo, E. W. Olszewski, O. Y. Gnedin, X. W ang, B. Sen et al., Velocity Dispersion Profiles of Seven Dwarf Spheroidal Galaxies , Astrophys. J. Lett. 667 (2007) L53 [0708.0010]

  34. [34]

    N. Bar, S. Danieli and K. Blum, Dynamical Friction in Globular Cluster-rich Ultra-diffuse Galaxies: The Case of NGC5846-UDG1 , Astrophys. J. Lett. 932 (2022) L10 [2202.10179]. 13

  35. [35]

    Binney and S

    J. Binney and S. Tremaine, Galactic Dynamics . Princeton University Press, Princeton, NJ, 2 ed., 2008

  36. [36]

    V. M. Gorkavenko, O. V. Barabash, T. V. Gorkavenko, O. M. Teslyk, A. O. Zaporozhchenko, J. Jia et al., Dynamical friction in rotating ultralight dark matter galactic cores, Class. Quant. Grav. 41 (2024) 235013 [2408.00104]

  37. [37]

    J. A. Petts, A. Gualandris and J. I. Read, A semi-analytic dynamical friction model that reproduces core stalling, Monthly Notices of the Royal Astronomical Society 454 (2015) 3778

  38. [38]

    Kaur and N

    K. Kaur and N. C. Stone, Density wakes driving dynamical friction in cored potentia ls, Monthly Notices of the Royal Astronomical Society 515 (2022) 407

  39. [39]

    Banik and F

    U. Banik and F. C. van den Bosch, A self-consistent, time-dependent treatment of dynamical friction: New insights regarding core stalling a nd dynamical buoyancy , The Astrophysical Journal 912 (2021) 43

  40. [40]

    Di Cintio and B

    P. Di Cintio and B. Marcos, A disturbance in the force: How force fluctuations hinder dynamical friction and induce core stalling , Astronomy & Astrophysics 700 (2025) A230

  41. [41]

    M. D. Weinberg, Evolution of barred galaxies by dynamical friction , Monthly Notices of the Royal Astronomical Society 213 (1985) 451. 14