The work sits off to the side among studies of black holes in dwarf galaxies: its few neighbors examine invisible intermediate-mass black holes, the shadows of giant black holes, and their connection to dark matter. Unlike most papers, which either look for signs of a black hole in radiation or model its influence on stars, this one tackles localization: where exactly to search. The Large Magellanic Cloud is the nearest massive dwarf galaxy, and if it hosts a supermassive black hole, that would be a unique chance to study how such objects form at the low-mass end. But no direct detection exists, and one reason is that the black hole may be offset from the photometric or kinematic center because of the recent interaction with the Small Magellanic Cloud and the tidal field of the Milky Way. Different methods of determining the center give conflicting results, which makes a targeted search uncertain.
The authors built a semi-analytic model in which the supermassive black hole is treated as an ensemble of test particles moving in a time-varying gravitational potential of the Large Magellanic Cloud, Small Magellanic Cloud, and Milky Way system. The potential of the Large Magellanic Cloud consists of a spherical dark matter halo, a stellar disk, and a rotating bar. The key element is a wobbling disk: its offset relative to the halo center mimics the response to the recent collision with the Small Magellanic Cloud, which excited the m=1 mode and created a displacement between the disk and the bar. Orbits were integrated with the galpy library in steps of about one million years over roughly 156 million years after the last close approach of the galaxies. Initial conditions for the black hole were drawn from probability distributions: a cylindrical radius from a lognormal distribution centered at 0.5 kpc, a vertical coordinate from a Gaussian with a dispersion of 100 pc, and velocities from an isotropic distribution with a dispersion of about 35 km/s. In total, 100,000 orbits were computed, taking dynamical friction into account via the Chandrasekhar formula. The positions were then projected onto the celestial sphere and a probability map was constructed.
The resulting probability map is asymmetric and elongated. The mean position of the predicted black hole has coordinates (α, δ) = (80.23°, −69.55°), and the 1σ error ellipse has semi-axes of 1.34° and 0.56° with a position angle of about 93.5° east of north. This mean lies roughly 6 arcminutes north of the accepted dynamical center, the bar center traced by red clumps. The probability peak (the mode) essentially coincides with the dynamical center and lies deep inside the 1σ region. All literature estimates of the center of the Large Magellanic Cloud fall within 2σ. So there is no statistically significant offset of the black hole from the center. An analytic estimate shows that the tidal field of the Milky Way currently exceeds that of the Small Magellanic Cloud by at least an order of magnitude, setting a preferred direction for vertical displacements, while disk asymmetry dominates in-plane displacements. Typical displacement amplitudes are tens of parsecs from external tides and hundreds of parsecs from internal asymmetry. The result is robust to changes in black hole mass over the range from 6×10⁴ to 6×10⁶ solar masses: the mean position shifts by no more than one arcminute.
The work differs from neighboring studies in that it does not look for signs of a black hole but predicts its likely position on the sky, turning an open-ended task into a targeted one. Unlike papers that study black holes through their shadows or their influence on dark matter, this is the first quantitative probability map for the Large Magellanic Cloud based on the dynamics of its interaction with the Small Magellanic Cloud and the tidal field of the Milky Way. That makes it possible to direct future observations to a specific region and shrink the search area. It also shows that even in a dynamically complex system, a testable prediction can be obtained by combining analytic estimates with numerical modeling.
The topic could advance through refined initial conditions: a broader distribution of the black hole's initial positions and velocities, including wandering black holes with offsets up to a kiloparsec, would enlarge the final probability map but would not change its preferred direction. A more detailed history of the interaction with the Small Magellanic Cloud and the non-stationarity of the Milky Way potential could also be included. On the observational side, future surveys such as SDSS-V Local Volume Mapper could test the prediction using emission lines, for example the coronal [Fe X] 6375 Å line, in the central region.
The result will affect the strategy for searching for supermassive black holes in dwarf galaxies and the planning of spectroscopic surveys of nearby systems.
The next step is to carry out targeted observations of the predicted region using integral-field spectroscopy and kinematic studies of candidates for nuclear star clusters.
The work is directly tied to the unsolved problem of how supermassive black holes form at the low-mass end: detecting a black hole in the Large Magellanic Cloud would provide key data on how such objects arise in dwarf galaxies. It also touches on the question of wandering black holes, namely how often they are offset from the centers of their galaxies by interactions and tidal fields.
The predicted region overlaps with the nuclear star cluster candidate NGC 1916, which makes this cluster a particularly interesting target in the search for an accreting black hole. The distance to the Large Magellanic Cloud is about 49.9 kiloparsecs, and the predicted mean position of the black hole turns out to be roughly 100 parsecs closer to us than the accepted distance to the galaxy.