Fast-moving stars around an intermediate-mass black hole in ω Centauri

Discussion of previous IMBH detections in ω Centauri

The debate about an IMBH in ω Cen dates back almost two decades but has remained controversial. Early dynamical modelling based on LOS integrated-light velocity dispersion measurements suggested an IMBH mass of \({M}_{{\bf{IMBH}}}=\left({4.0}_{-1.0}^{+0.75}\right)\times {10}^{4}{M}_{\odot }\) (ref. ³). These results were challenged with a precise redetermination of the centre of the cluster⁶ and dynamical modelling of proper motions⁷ measured from multi-epoch HST imaging observations that placed an upper limit of 1.2 × 10⁴M_☉ on the IMBH. Using additional integrated light observations and a centre based on the maximum LOS velocity dispersion, a best-fit IMBH mass of \({M}_{{\bf{IMBH}}}=\left(4.7\pm 1.0\right)\times {10}^{4}{M}_{\odot }\) was obtained in ref. ⁴. When assuming the AvdM10 centre, the IMBH mass was slightly lower, \({M}_{{\bf{IMBH}}}=\left(3.0\pm 0.4\right)\times {10}^{4}{M}_{\odot }\).

Subsequent comparisons of both proper motions and LOS velocities to N-body simulations continued to show evidence for an approximately 4.0 × 10⁴M_☉ IMBH^31,32. However, these observations were also shown to be fully consistent with a dark cluster of stellar mass black holes in the central region of ω Cen⁸. The lack of fast-moving stars in previous proper-motion catalogues supported this scenario over an IMBH⁹. Other works noted the influence of radial velocity anisotropy on dynamical mass estimates^33,34.

Most recently, the discovery of a counter-rotating core using VLT MUSE LOS velocity measurements of individual stars³⁵ highlighted once again the kinematic complexity of the centremost region of ω Cen. The centre of this counter-rotation coincides with the AvdM10 centre within about 5″ but is incompatible with the centres used in refs. ^3,4.

Previous accretion constraints in context

With the detection of a 10⁴⁻⁵M_☉ IMBH, the upper limits on any accretion signal at X-ray³⁶ and radio³⁷ wavelengths make this the most weakly accreting black hole known. Deep, 291 ksec Chandra observations place an upper limit of the 0.5–7 keV luminosity of around 10³⁰ ergs s⁻¹ (ref. ³⁶), roughly 12 orders of magnitude below the Eddington limit. The radio upper limit implies an even fainter source, with the 5 GHz upper limit of 1.3 × 10²⁷ ergs s⁻¹ (ref. ³⁷) corresponding to an implied X-ray luminosity using the fundamental plane of about 10²⁹ ergs s⁻¹ (ref. ³⁸). Assuming standard bolometric corrections of about 10 (ref. ³⁹), this X-ray luminosity upper limit suggests an Eddington ratio of \(\log ({L}_{{\rm{bol}}}/{L}_{{\rm{edd}}}) < -12\), far fainter than that for Sgr A*⁴⁰ or any other known black hole. This faint signal could be because of a combination of low surrounding gas density, a low accretion rate of that gas and/or a low radiative efficiency³⁷. Low-luminosity active galactic nuclei including Sgr A* are brightest at infra-red and sub-millimetre wavelengths most likely because of synchrotron emission from compact jets^41,42. Therefore, future observations with the James Webb Space Telescope or the Atacama Large Millimeter array would provide the highest sensitivity to any emission from the IMBH of ω Cen. Any detection would reveal the location of the IMBH as well as provide valuable constraints on the black hole accretion in this extremely faint source.

Proper-motion measurements and sample selection

Our proper-motion measurements are based on the reduction of archival HST data of the central region of ω Cen, taken over a time span of more than 20 years. We used the state-of-the-art photometry tool KS2 (ref. ³²) for the source detection and the astro-photometric measurements, and the established procedure described in refs. ^43,44,45,46 to measure proper motions relative to the bulk motion of the cluster. The result of this extensive study is a proper-motion catalogue with high-precision measurements for 1.4 million stars out to the half-light radius of ω Cen with a typical temporal baseline of more than 20 years. Owing to the large number of observations (in total we reduced over 500 images and some stars in the central region have up to 467 individual astrometric measurements) the catalogue reaches unprecedented depth and precision. The highest precision is achieved in the well-covered centre of the cluster, in which our proper motions have a median error of only about 6.6 µ as yr⁻¹ (0.17 km s⁻¹) per component for bright stars.

The catalogue is larger than any other kinematic catalogue published for a globular cluster and significantly extends previous proper-motion catalogues for ω Cen^6,32,47. A detailed comparison with other proper-motion datasets is published along with the catalogue¹⁴. In a following section and in Extended Data Fig. 3, we compare the completeness of the different catalogues to show that it is plausible that the fast-moving stars have been missed in previous searches.

We use a high-quality subset of the proper-motion catalogue to search for real fast-moving stars and limit spurious astrometric measurements (for example, two sources that are falsely identified as one) that can have apparent high proper-motion measurements. Our criteria for this subset are based on the amount of available data for the measurements. Specifically, we used only sources that had at least 20 astrometric measurements covering a temporal baseline of at least 20 years and a fraction of rejected measurements (based on sigma clipping) of less than 15%. We also made cuts on the quality of the proper-motion fit requiring both a proper-motion error less than 0.194 mas yr⁻¹ ≈ 5 km s⁻¹ and a reduced χ² < 10 for the linear proper-motion fit for both the right ascension (RA) and declination (Dec) measurements. Apart from these quality selections, we also required the star to lie on the CMD sequence in an HST-based CMD (Fig. 2). These cuts help to drastically reduce the number of contaminants. These criteria are met uniformly out to a radius of about 90 arcsec; at larger radii they lead to selection effects due to reduced observational coverage. A total of 157,320 out of 241,133 (65.2%) entries of the proper-motion catalogue within r < 90″ match the combined criteria. Extended Data Table 1b shows the individual measured proper-motion components for the seven fast-moving stars. We note that for this analysis we have not applied the local a posteriori proper-motion corrections provided with the catalogue¹⁴, as we are studying the central region that is well dithered and observed with various rotation angles. We verified that applying these corrections would neither change our fast-moving star sample nor alter our conclusions.

Details on verification for fast-moving stars

The criteria detailed above should lead to a clean dataset with very few spurious proper-motion measurements. To ensure that the measurements for the fast-moving stars are reliable, we inspected each of them carefully.

As a first step, we tested the quality of the raw astrometric measurements by studying several goodness-of-fit parameters and photometric quality indicators for the point-spread-function fits used to measure stellar positions (Extended Data Fig. 3). We performed this analysis for the WFC3/UVIS F606W filter as it is the most used filter in the centre of ω Cen, and each star has at least 195 measurements in this filter. To verify the goodness of fit, we used the mean of the so-called quality-of-fit (QFIT) flag and the radial excess value, both of which take into account the residuals of the point-spread-function fit. Furthermore, we looked at the mean of the ratio of source flux with respect to the flux of neighbouring sources within its fit aperture. All five stars used for our analysis behave typically for well-measured stars of their magnitude and none of them show extreme values that would indicate problems with the photometry. It is noteworthy that the two stars excluded from our analysis based on their velocities being less than 3σ above the escape velocity show some deviations. Star B has a relatively low mean QFIT value and high radial excess. Star G is the only one in the sample for which the neighbour flux to star flux ratio is larger than 1.

As a second step, we looked at the stars in several stacked HST images taken at various epochs ranging from 2002 to 2023. The extensive multi-epoch imaging is demonstrated in Extended Data Fig. 2. The proper motion of a star at the escape velocity (62 km s⁻¹) is 2.41 mas yr⁻¹. Therefore, we expect to see a displacement of at least 50 mas over 21 years, which corresponds to 1.25 WFC3/UVIS pixels. This motion can be seen by eye for all seven stars in the multi-epoch images. Again the excluded stars B and G stick out, that is, that they are partially blended with the neighbouring stars, thus explaining their larger astrometric errors.

Finally, we tested the reliability of our proper-motion measurements by limiting the raw position measurements to different subsets and redoing both the linear and quadratic fits to the motion of the stars. The first test run included only high signal-to-noise measurements. The second test run included only measurements taken with the WFC3/UVIS F606W filter. By using only one filter, we are immune to colour-induced effects such as a partially resolved blend between two differently coloured stars. The proper motions of all fast-moving stars are consistent within the measurement uncertainties using both methods.

Comparison with other proper-motion datasets

Before our analysis, two other high-precision proper-motion catalogues based on HST data have been published^6,32 covering the centre of ω Cen. Both datasets were searched for central high proper-motion stars, but none of the stars in our sample have been reported before. To understand why this is the case, we compare the completeness of the different catalogues. Extended Data Fig. 3a shows histograms of the magnitudes of stars with measured proper motions. In the inner 20″, our new catalogue contains more than three times the number of stars of the literature catalogues and extends to significantly fainter magnitudes. The newly detected fast-moving stars all lie at faint magnitudes, in which the completeness of the older catalogues is significantly lower than in the new proper-motion catalogue. This is because of the larger amount of data and the updated source-finding algorithms in the new catalogue, which explains the previous non-detection of the fast-moving stars.

Discussion of photometric errors

Apart from the astrometric reliability, we also studied the quality of the photometric measurements used to locate the stars in the CMD. Although the innermost stars are faint, their statistical photometric errors are small because of the large number of individual photometric measurements combined to a weighted mean value. The statistical errors range from 0.004 mag to 0.037 mag and are given in Extended Data Table 1c. However, especially for faint stars, this statistical error is not able to capture systematic issues caused, for example, by the influence of brighter neighbouring stars. These issues can be identified only by verifying the quality of the point-spread-function fit used to determine the individual photometric measurements. We report the mean quality of fit, radial excess and neighbour flux to source flux ratio flags for both filters in Extended Data Table 1c and compare them with those of stars at similar (Δm < 0.5) magnitudes in Extended Data Fig. 3. All stars in the robustly measured sample show typical quality of fit for their respective magnitude. We note that stars B and G (which were excluded from the analysis) show comparatively poor QFIT. Stars E and G show a possible flux contribution from a neighbouring source (indicated by a high radial excess value and a high neighbour flux to source flux ratio). This can be confirmed by the stacked images shown in Extended Data Fig. 2, in which these stars show a close neighbour. Owing to the relatively bright magnitude of star E and the low astrometric scatter, we still consider its measurement valid.

Density of Milky Way contaminants versus fast-star background

To quantify our expected level of contamination from Milky Way foreground and background stars, we compared our results with those of a Besançon model¹⁸. Using the m1612 model, we simulate a 1 square degree patch centred on ω Cen retrieving Johnson colours and kinematics. We then transform the model Johnson V and I magnitudes into F606W and F814W magnitudes using linear relations fitted to Padova models⁴⁸ between V–I of 0 and 2. We use the same colour cuts and consider stars between F606W of 16–24. We then count the number of stars with a total proper motion above 2.41 mas yr⁻¹ (equivalent to our velocity cutoff at the escape velocity v_esc = 62 km s⁻¹). These would appear as contaminants in our fast-moving star sample. We find a density of 0.0039 stars per arcsec². This is somewhat higher than the 0.0026 ± 0.0003 stars per arcsec² found as the background level in our observations. This discrepancy is alleviated by considering only 65.2% of all stars within our catalogue meet the high-quality criteria used for our fast-moving star selection (see above). Correcting for this factor, we get an expected background of 0.0025 stars per arcsec², perfectly matching the observed background density (Extended Data Fig. 1). This suggests that our background level is consistent with being predominantly Milky Way contaminants. Relaxing our requirement that stars be more than 3σ above the escape velocity results in a higher observed background level of 0.0042 stars per arcsec², no longer consistent with the Milky Way background. This suggests that our stricter definition of a fast-moving star reduces contamination from poorly measured stars in ω Cen to a negligible level.

Other scenarios that could explain the fast-moving stars

A complete contamination of our sample by Milky Way foreground and background stars that are non-members of ω Cen can be ruled out statistically. We now explore and rule out alternative scenarios to the fast-moving stars being bound to an IMBH. One alternative explanation for stars with a high velocity is to have them bound in a close orbit with a stellar-mass black hole. This scenario can be ruled out for BHs less than 100M_☉, as the periods required to reach the observed velocities are less than 10 years, which is well within the 20-year span over which we have observed linear motions.

Another scenario could be that the stars are actually unbound from the cluster and have recently been accelerated by three-body interactions, either with stellar-mass black hole binaries or with an IMBH. Ejection by an IMBH in the centre of the cluster can be ruled out by the very high rate of ejections necessary to sustain the observed number of fast-moving stars within the centre and the absence of observed fast-moving stars at larger radii. To sustain a density of 0.18 fast-moving stars per arcsec² in the inner 3 arcsec (equivalent to our conservative sample of five stars) moving with at least 2.4 mas yr⁻¹ would require ejections with a rate of 0.004 stars per year. This would lead to about 117 additional fast-moving stars at larger radii (20″ < r < 90″), apart from around 60 foreground stars expected from the (completeness corrected) Besançon Milky Way model. In our dataset, we find 61 fast-moving stars between 20″ < r < 90″, consistent with the expected Milky Way background but not consistent with a substantial number of additional ejected stars. Moreover, a high hypothetical ejection rate of 0.004 stars per year would deplete all of the about 10 million stars of ω Cen in just 2.5 Gyr. If no IMBH is present, accelerations of stars above the escape velocity are still possible by three- or four-body interactions between stellar or compact object binaries⁴⁹. However, these interactions would not be limited to the innermost few arcseconds of the cluster because of the slowly varying stellar density in the core of ω Cen. Furthermore, the expected rate of these ejection events is of the order of less than one ejection per 1 million years, around 1,000 times lower than needed to explain the observed number of fast-moving stars in the centre of ω Cen^49,50.

Search for the fast-moving stars in recent LOS velocity data

Line-of-sight (LOS) velocities of the fast-moving stars could help to exclude contaminants and provide further constraints on the orbits of the stars and the mass and position of the IMBH. The deepest and most extensive spectroscopic catalogue of stars in ω Cen is part I of our recently published oMEGACat¹³. This catalogue was created using a large mosaic of observations with the VLT MUSE integral field spectrograph and contains both LOS velocity measurements and metallicities for more than 300,000 stars within the half-light radius of ω Cen. Although we could successfully cross-match five of the seven fast-moving stars, their signal-to-noise ratio is typically too low (about 2) for reliable velocity measurements, in particular for the four fastest, innermost stars.

We could, however, obtain an LOS velocity value for star E (v_LOS = 261.7 ± 2.7 km s⁻¹) and star F (v_LOS = 232.5 ± 4.0 km s⁻¹). These velocities are very close to the systemic LOS velocity of ω Cen (232.99 ± 0.06 km s⁻¹; ref. ¹³), confirming their membership in the cluster, as the Milky Way foreground is centred at v_LOS ≈ 0 with a dispersion of 70 km s⁻¹ (ref. ¹⁸). However, as the relative LOS velocity with respect to the cluster is low and we have only those two velocities for these outer stars, the LOS velocities do not add stronger constraints on the IMBH. For this reason, we did not include them in the rest of our analysis.

Testing the robustness of the assumed escape velocity

Varying the parameters of the N-body models

Because we use the escape velocity of ω Cen (assuming no IMBH is present) as the threshold for determining whether a star is considered fast or not, it is important to verify the robustness of the escape velocity value. We adopt an escape velocity v_esc = 62 km s⁻¹ (ref. ¹⁶); we have verified this value based on fitting similar N-body models to several state-of-the-art datasets, including MUSE LOS velocity dispersion measurements⁵¹ and HST proper-motion-based dispersion measurements⁵² for the central kinematics and Gaia DR3⁵³ measurements at larger radii using an assumed distance of 5.43 kpc. We varied both the assumed initial stellar mass function (using either the canonical Kroupa IMF⁵⁴ or the bottom-light IMF derived in ref. ⁵⁵) and the black hole retention fraction (assuming values of 10%, 30%, 50% or 100%). Despite changes to the central mass-to-light (M/L) ratio between the models, the central escape velocity changes only minimally, with a range of values from 61.1 km s⁻¹ to 64.8 km s⁻¹. Adopting any of these values leaves our sample of seven central stars above the escape velocity unchanged.

An independent test using surface-brightness profiles

As a second test, independent of the N-body models, we calculated an escape velocity profile based on a surface brightness profile using various surface brightness profiles and dynamical models from the literature. We started by parameterizing the surface brightness profile using multi-Gaussian expansion (MGE)⁵⁶ models. We then converted the surface brightness to a mass density using several mass-to-light ratios and distances used in the literature. From the mass density, we can derive the gravitational potential (Φ(r)). The escape velocity profile is then given by

with the tidal radius r_tidal = 48.6′ ≈ 74.6 pc from refs. ^57,58.

These tests showed that the central escape velocity does not depend strongly on the stellar mass distribution in the centremost region, instead it is dominated by the global M/L ratio and the assumed distance. The early dynamical models in the IMBH debate assumed both a distance of 4.8 kpc (ref. ⁵⁹) and an M/L of 2.6 (vdMA10; ref. ⁷) or 2.7 (N08; refs. ^3,4). With these values our tests give a central escape velocity of 55.4 km s⁻¹ (vdMA10) and 56.9 km s⁻¹ (N08). If we would also use the 4.8 kpc distance to scale the proper motions, this gives a cutoff of 2.43 mas yr⁻¹ (vdMA10) and 2.48 mas yr⁻¹ (N08), close to the adopted cutoff at 2.41 mas yr⁻¹ and not changing the sample of the seven fast-moving stars detected. Owing to the parallax measurements of the Gaia satellite and updated kinematic distance measurements, the distance to ω Cen was robustly redetermined and larger values have been found (5.24 ± 0.11 kpc; ref. ⁶⁰; 5.43 ± 0.05; ref. ¹²). A dynamical model using the same surface brightness profile as ref. ³ but a larger distance of \({5.14}_{-0.24}^{+0.25}\,{\rm{kpc}}\) was presented in ref. ⁸; this study found an M/L of \({2.55}_{-0.28}^{+0.35}\). Using these values, the central escape velocity derived from the surface brightness profile is 61.1 km s⁻¹ (equivalent to a proper motion of 2.51 mas yr⁻¹); again not changing our fast-moving stars sample. Finally, varying the distance of any model by 0.2 kpc while holding the M/L constant results in variation in the escape velocity of the order of ±3 km s⁻¹. These results show that our fast-moving star limit (v_esc = 62 km s⁻¹ at a distance of 5.43 kpc) is consistent with the escape velocity values directly derived from surface-brightness profiles and several dynamically estimated M/L values. To visualize the escape velocity, we calculated the escape velocity using the surface-brightness profile in ref. ³, an M/L of 2.4, and a distance of 5.43 kpc as found from the N-body models with a cluster of stellar mass black holes. The resulting profile is shown in Extended Data Fig. 5f. The predicted escape velocity is flat out to approximately 50″, a property shared by all of the calculated escape velocity profiles. This makes the detection of the fast-moving stars only in the central few arcseconds more compelling.

An empirical confirmation of the central escape velocity

We make one final, and relatively model-independent, empirical confirmation of the central escape velocity based on the distribution of 2D velocities in the innermost region of ω Cen (Extended Data Fig. 4). As we have seen in the analysis above, the escape velocity varies only slightly within the inner about 50″ of the core of ω Cen. Furthermore, the velocity dispersion profile is relatively flat in the innermost 10″, with a value of about 20 km s⁻¹ (refs. ^6,35,52). Therefore, we would expect rather similar distributions of stellar velocities in both the centre (0″ < r < 3″) and an outer ring at (3″ < r < 10″). Although we observe a clear excess of fast-moving stars in the inner 3 arcsec, there is a sharp cutoff very close to the adopted escape velocity in the (3″ < r < 10″) bin. Although there is a total of 2,090 stars, there is only one star with a velocity significantly faster than the escape velocity (instead of 17 stars expected from a 2D Maxwell–Boltzmann distribution with σ_1D = 20 km s⁻¹). This suggests that the stars with these velocities have escaped the central region. This outer fast-moving star has a 2D velocity of 75.8 km s⁻¹ and is at a radius of r = 9.5″. From the density of Milky Way contaminants with apparent velocities above the escape velocity, we would expect about 0.7 foreground stars in the (3″ < r < 10″) region; therefore, this fast-moving star is consistent with being a Milky Way foreground star.

The escape velocity provides a minimum black hole mass

The escape velocity for an isolated black hole is given by

In ω Cen, we have to take into account the potential of the globular cluster as well. If we assume this to be constant over the very small region in which we found the fast-moving stars (an assumption that agrees with the published surface brightness profiles, see Extended Data Fig. 6f), we obtain

If a star at the distance of r_3D with a velocity v_3D is bound to the black hole, we can calculate the following lower limit on the black hole mass:

A lower limit can also be calculated if the LOS velocity and distance are not known, as v_3D ≥ v_2D and r_3D  ≥ r_2D. As we do not know the exact 2D position of the black hole relative to the fast-moving stars, we calculated this lower limit for all stars and a grid of assumed 2D locations around the AvdM10 centre⁶. Each individual star alone would allow for a very low mass, as the location of the black hole could coincide with the star (Extended Data Fig. 5a–d). However, combining these limits for all stars gives a higher minimum black hole mass (Extended Data Fig. 5e). If we assume that all five robustly detected stars are bound to the black hole, the lower limit is about 8,200M_☉ and the minimum mass location is only 0.3″ away from the AvdM10 (ref. ⁶) centre at the location RA = 201.6967370° and Dec = −47.4795066°. If we assume that the two most constraining stars are just random foreground contaminants, which is ruled out at the 3σ level (P = 0.0026), this limit drops to about 4,100M_☉, still well within the IMBH range.

Acceleration measurements

The astrometric analysis in the catalogue¹⁴ considered only linear motions of the stars. If there is a massive black hole present near the centre, we might also be able to measure accelerated motion of the closest stars, allowing for a direct mass measurement of the black hole. With an IMBH mass of 40,000M_☉ and at a radius of 0.026 pc (1″ on the sky), the acceleration of a star would be 0.25 km s⁻¹ yr⁻¹ (or 0.01 mas yr⁻²). This is at the limit of the precision of our current dataset. With a 20-year baseline, we expect only a deviation of 0.05 pixel from a linear motion. For bright stars, the astrometric uncertainty can be as low as 0.01 pixel; however, for the faint fast-moving stars we have detected, the errors are significantly larger.

To constrain these possible accelerations, we repeated the fit of the motion of each star enabling the addition of a quadratic component. The results for this fit are shown in Extended Data Table 1. The accelerations of the stars are consistent with zero within 3σ, but two stars have more than 2σ acceleration measurements. The errors on our acceleration measurements lie between 0.004 mas yr⁻² and 0.03 mas yr⁻² and are, therefore, of a magnitude similar to the expected acceleration signal. The strongest acceleration is shown by star B, which has been excluded from the robust subset of fast-moving stars because its proper motion is not 3σ above the escape velocity. Owing to the proximity of a bright neighbour star, we do not deem this acceleration measurement to be reliable.

As the LOS distances to the fast-moving stars are unknown, it is not possible to place direct constraints on the IMBH mass using the upper limits on accelerations. If no acceleration is detected, as is the case for the centremost star A, this could mean that either the black hole is not very massive or the LOS distance of star A to the black hole is large. Combining the measurements for the ensemble of fast-moving stars and making some assumptions on their spatial distribution still enables us to use the acceleration limits to place further constraints on the black hole mass and its location. This is described in the next section.

Markov chain Monte Carlo fitting of the acceleration data

Assuming the fast-moving stars are bound to the IMBH, we can model the stars as being on Keplerian orbits around the IMBH. We used Bayesian analysis to sample the posterior distribution for the unknown mass and position of the black hole. In this analysis, there were eight free parameters: black hole mass; its on-sky x– and y-positions; and five LOS distances between the black hole and each fast-moving star. This analysis makes use of the available astrometric observations but stops short of modelling individual stellar orbits that would introduce more free parameters.

We use a likelihood function with these eight free parameters and give the likelihood based on the observed on-sky x– and y-acceleration, proper motion and position of the five robustly measured fast-moving stars. For each star, we calculated a first likelihood term based on the modelled acceleration a_modelled using a Gaussian distribution with mean a_observed and width equal to the acceleration uncertainty. The second term in the likelihood accounts for the escape velocity constraints and is kept constant if the observed 2D velocity of the star is below the modelled escape velocity. For stars with 2D velocities above the modelled escape velocity, the likelihood is a Gaussian distribution with mean v_2D − v_esc total and width equal to the uncertainty in observed proper motion.

We make these assumptions about the model: (1) The black hole mass is between 1M_☉ and 100,000M_☉, because a black hole mass beyond this upper limit is ruled out by our N-body models. (2) The black hole is located within the distribution of the fast-moving stars. We use a Gaussian prior in the black hole x– and y-positions with a mean equal to the mean position of the fast-moving stars and width equal to their one-dimensional positional standard deviation, σ_stars = 0.0221 pc; we also use a cutoff at ±0.16 pc. (3) The stellar positions are isotropically distributed around the black hole. We model the LOS positions of the stars relative to the black hole using a Gaussian distribution with mean 0 and width σ_stars.

The posterior was sampled using a Markov chain Monte Carlo (MCMC) ensemble sampler implemented using the package emcee⁶¹ using recommended burn-in and autocorrelation corrections. We show the posterior distribution for the black hole mass in Extended Data Fig. 6a. The 99% confidence lower limit (21,100M_☉) is significantly higher than that derived from escape velocity constraints alone, whereas the upper limit on the mass is not well constrained. We also find a position for the black hole east of the AvdM10 centre, with \(\varDelta x={-0.017}_{-0.031}^{+0.017}\,{\rm{pc}}\) and \(\varDelta y={0.011}_{-0.025\,}^{+0.011}{\rm{p}}{\rm{c}}\) (Extended Data Fig. 6b). The coordinates of the MCMC based centre estimate are RA = 201.6970128° and Dec = −47.4794533°. We note that the black hole location estimate is dominated by the marginal 2σ acceleration signal of star D, which is the faintest star in the sample; follow-up studies are required to obtain more precise acceleration measurements.

N-body models

Apart from the analysis of stars with velocities above the escape velocity, we also used a set of existing N-body models with and without central IMBHs to get further constraints on the IMBH mass. We compared the simulations to the full velocity dispersion and surface density profile of ω Cen to determine the best-fitting model and the mass of a central IMBH. The set of models and the details of the fitting procedure are described in detail in refs. ^5,9. We note that these models have been presented already in the literature, but the fits to these models have been updated to incorporate the most recent Gaia DR3 data.

In short, the models started from King profiles⁶² with central concentrations between c = 0.2 and c = 2.5 and initial half-mass radii between r_h = 2 pc and r_h = 35 pc. In the models with an IMBH, we varied the mass of the IMBH so that it contains either 0.5%, 1%, 2% or 5% of the cluster mass at T = 12 Gyr when the simulations were stopped. The models with an IMBH assumed a retention fraction of stellar-mass black holes of 10%, whereas in the models without an IMBH we varied the assumed retention fraction of stellar-mass black holes between 10% and 100%. At the end of the simulations, we calculated surface density and velocity dispersion profiles for each N-body model and then determined the best-fitting model by interpolation in our grid of models and using χ² minimization against the observed velocity and surface density profile of ω Cen.

The velocity distributions from observations and the models are shown in Extended Data Fig. 7. We compare the distribution of measured 2D stellar velocities in the inner 10″ of ω Cen with the various models using a Kolmogorov–Smirnov test. Moreover, we compare the fraction of fast-moving stars in the innermost 3 arcsec (Extended Data Table 2). Models without an IMBH and with a 20,000M_☉ IMBH are both strongly excluded by both the overall velocity distribution and the complete lack of fast-moving stars. The overall velocity distribution is in best agreement with the 47,000M_☉ distribution, whereas the fraction of fast-moving stars is best matched by the 39,000M_☉ simulation. This tension might be alleviated in future models that contain both an IMBH and a cluster of stellar mass black holes. We caution that these simulations have smaller numbers of stars than observed and that there can be substantial variations in the distribution of central stars due to strong encounters with remnants and binaries. Nonetheless, the simulations suggest that black holes with masses of M ≲ 50,000M_☉ are consistent with the observed distribution of central velocities and fast-moving stars, whereas the no-IMBH case and significantly more massive black holes are disfavoured because of an overprediction of fast-moving stars. Updated N-body models fit to the oMEGACat kinematic data and dynamical modelling of these same datasets with Jeans models are currently underway.

Comparison with S-stars in the Galactic centre

The black hole indicated by our fast-moving star detection is only the second after Sgr A*, for which we can study the motion of multiple individual bound stellar companions. Therefore, the extensively measured stars around Sgr A* provide a unique comparison point to our fast-moving star sample. We compare the motions of the stars in the S-star catalogue in ref. ¹⁷ with our ω Cen fast-moving star sample in Extended Data Fig. 8. When taking into account the different distances and the approximate black hole mass ratio of 100, the motions show similar amplitudes. However, the density of tracers in ω Cen is significantly lower, despite the greater depth of the observations.