Pericenter Passage Research Articles

The Dawes Review 9: The role of cold gas stripping on the star formation quenching of satellite galaxies

AbstractOne of the key open questions in extragalactic astronomy is what stops star formation in galaxies. While it is clear that the cold gas reservoir, which fuels the formation of new stars, must be affected first, how this happens and what are the dominant physical mechanisms involved is still a matter of debate. At least for satellite galaxies, it is generally accepted that internal processes alone cannot be responsible for fully quenching their star formation, but that environment should play an important, if not dominant, role. In nearby clusters, we see examples of cold gas being removed from the star-forming discs of galaxies moving through the intracluster medium, but whether active stripping is widespread and/or necessary to halt star formation in satellites, or quenching is just a consequence of the inability of these galaxies to replenish their cold gas reservoirs, remains unclear. In this work, we review the current status of environmental studies of cold gas in star-forming satellites in the local Universe from an observational perspective, focusing on the evidence for a physical link between cold gas stripping and quenching of the star formation. We find that stripping of cold gas is ubiquitous in satellite galaxies in both group and cluster environments. While hydrodynamical mechanisms such as ram pressure are important, the emerging picture across the full range of dark matter halos and stellar masses is a complex one, where different physical mechanisms may act simultaneously and cannot always be easily separated. Most importantly, we show that stripping does not always lead to full quenching, as only a fraction of the cold gas reservoir might be affected at the first pericentre passage. We argue that this is a key point to reconcile apparent tensions between statistical and detailed analyses of satellite galaxies, as well as disagreements between various estimates of quenching timescales. We conclude by highlighting several outstanding questions where we expect to see substantial progress in the coming decades, thanks to the advent of the Square Kilometre Array and its precursors, as well as the next-generation optical and millimeter facilities.

A homogeneous measurement of the delay between the onsets of gas stripping and star formation quenching in satellite galaxies of groups and clusters

ABSTRACT We combine orbital information from N-body simulations with an analytic model for star formation quenching and SDSS observations to infer the differential effect of the group/cluster environment on star formation in satellite galaxies. We also consider a model for gas stripping, using the same input supplemented with H i fluxes from the ALFALFA survey. The models are motivated by and tested on the Hydrangea cosmological hydrodynamical simulation suite. We recover the characteristic times when satellite galaxies are stripped and quenched. Stripping in massive ($M_{\rm vir}\sim 10^{14.5}\, {\rm M}_\odot$) clusters typically occurs at or just before the first pericentric passage. Lower mass ($\sim 10^{13.5}\, {\rm M}_\odot$) groups strip their satellites on a significantly longer (by $\sim 3\, {\rm Gyr}$) time-scale. Quenching occurs later: Balmer emission lines typically fade $\sim 3.5\, {\rm Gyr}$ ($5.5\, {\rm Gyr}$) after first pericentre in clusters (groups), followed a few hundred Myr later by reddenning in (g − r) colour. These ‘delay time-scales’ are remarkably constant across the entire satellite stellar mass range probed (∼109.5–$10^{11}\, {\rm M}_\odot$), a feature closely tied to our treatment of ‘group pre-processing’. The lowest mass groups in our sample ($\sim 10^{12.5}\, {\rm M}_\odot$) strip and quench their satellites extremely inefficiently: typical time-scales may approach the age of the Universe. Our measurements are qualitatively consistent with the ‘delayed-then-rapid’ quenching scenario advocated for by several other studies, but we find significantly longer delay times. Our combination of a homogeneous analysis and input catalogues yields new insight into the sequence of events leading to quenching across wide intervals in host and satellite mass.

Creating a galaxy lacking dark matter in a dark matter-dominated universe

ABSTRACT We use hydrodynamical cosmological simulations to show that it is possible to create, via tidal interactions, galaxies lacking dark matter (DM) in a DM-dominated universe. We select dwarf galaxies from the NIHAO project, obtained in the standard cold dark matter model and use them as initial conditions for simulations of satellite–central interactions. After just one pericentric passage on an orbit with a strong radial component, NIHAO dwarf galaxies can lose up to 80 per cent of their DM content, but, most interestingly, their central (≈8 kpc) DM-to-stellar mass ratio changes from a value of ∼25, as expected from numerical simulations and abundance matching techniques, to roughly unity as reported for NGC 1052-DF2 and NGC 1054-DF4. The stellar velocity dispersion drops from ∼30 $\, \rm km\, s^{-1}$ before infall to values as low as 6 ± 2 $\, \rm km\, s^{-1}$. These, and the half-light radius around 3 kpc, are in good agreement with observations from van Dokkum and collaborators. Our study shows that it is possible to create a galaxy without DM starting from typical dwarf galaxies formed in a DM-dominated universe, provided they live in a dense environment.

Discovery of Extended Tidal Tails around the Globular Cluster Palomar 13

Abstract We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ∼5° on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar 13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars (RRLs) that are plausibly associated with the tails, in addition to four previously known in the cluster. From these RRLs, we find that the mean distance to the cluster and tails is 23.6 ± 0.2 kpc and estimate the total (initial) luminosity of the cluster to be , consistent with previous claims that its initial luminosity was higher than its current luminosity. Combined with previously determined proper motion and radial velocity measurements of the cluster, we find that Palomar 13 is on a highly eccentric orbit (e ∼ 0.8) with a pericenter of ∼9 kpc and an apocenter of ∼69 kpc, and a recent pericentric passage of the cluster ∼75 Myr ago. We note a prominent linear structure in the interstellar dust map that runs parallel to the candidate tidal features, but conclude that reddening due to dust is unlikely to account for the structure that we observe. If confirmed, the Palomar 13 stellar stream would be one of very few streams with a known progenitor system, making it uniquely powerful for studying the disruption of globular clusters, the formation of the stellar halo, and the distribution of matter within our Galaxy.

Joint gas and stellar dynamical models of WLM: an isolated dwarf galaxy within a cored, prolate DM halo

ABSTRACT We present multitracer dynamical models of the low-mass (M* ∼ 107), isolated dwarf irregular galaxy WLM in order to simultaneously constrain the inner slope of the dark matter (DM) halo density profile (γ) and flattening (qDM), and the stellar orbital anisotropy (βz, βr). For the first time, we show how jointly constraining the mass distribution from the H i gas rotation curve and solving the Jeans equations with discrete stellar kinematics lead to a factor of ∼2 reduction in the uncertainties on γ. The mass-anisotropy degeneracy is also partially broken, leading to reductions on uncertainty by ${\sim} 30{{\ \rm per\ cent}}$ on Mvir (and ${\sim} 70{{\ \rm per\ cent}}$ at the half-light radius) and ${\sim} 25{{\ \rm per\ cent}}$ on anisotropy. Our inferred value of γ = 0.3 ± 0.1 is robust to the halo geometry, and in excellent agreement with predictions of stellar feedback-driven DM core creation. The derived prolate geometry of the DM halo with qDM = 2 ± 1 is consistent with Lambda cold dark matter simulations of dwarf galaxy haloes. While self-interacting DM (SIDM) models with σ/mX ∼ 0.6 can reproduce this cored DM profile, the interaction events may sphericalize the halo. The simultaneously cored and prolate DM halo may therefore present a challenge for SIDM. Finally, we find that the radial profile of stellar anisotropy in WLM (βr) follows a nearly identical trend of increasing tangential anisotropy to the classical dwarf spheroidals, Fornax and Sculptor. Given WLM’s orbital history, this result may call into question whether such anisotropy is a consequence of tidal stripping in only one pericentric passage or if it instead is a feature of the largely self-similar formation and evolutionary pathways for some dwarf galaxies.

Evolution of inspiralling neutron star binaries: Effects of tidal interactions and orbital eccentricities

Neutron star (NS) binaries formed dynamically may have significant eccentricities while emitting gravitational waves (GWs) in the LIGO/VIRGO band. We study tidal effects in such eccentric inspiralling NS binaries using a set of hybrid equations. The NS is modelled as a compressible ellipsoid, which can deform nonlinearly due to tidal forces, while the orbit evolution is treated with the post-Newtonian (PN) theory up to 2.5-PN order. We find that in general, the tidal interaction can accelerate the inspiral, and cause orbital frequency and phase shifts. For circular inspirals, our calculations reproduce previous linear result at large binary separations, but incorporate the dynamical response of the NS at small separations. For eccentric inspirals, the frequency and phase shifts oscillate considerably near pericenter passages, and the oscillating amplitudes increase with eccentricities. As a result, the GW phase is also significantly influenced by the tidal effect. At merger, the cumulative GW phase shift can reach more than 10 radians (for typical NS mass $1.4M_\odot$ and radius 11.6 km), much larger than the circular inspiral case. Although the event rate of eccentric NS mergers is likely low, the detection of such mergers could provide a useful constraint on the NS equation of state.

A simple random walk model explains the disruption process of hierarchical, Eccentric three-body systems

ABSTRACT We study the disruption process of hierarchical three-body systems with bodies of comparable mass. Such systems have long survival times that vary by orders of magnitude depending on the initial conditions. By comparing with three-body numerical integrations, we show that the evolution and disruption of such systems can be statistically described as a simple random walk process in the outer orbit’s energy, where the energy exchange per pericenter passage (step size) is calculated from the initial conditions. In our derivation of the step size, we use previous analytic results for parabolic encounters, and average over the (Kozai–Lidov) oscillations in orbital parameters, which are faster then the energy diffusion time-scale. While similar random walk models were studied before, this work differs in two manners: (a) this is the first time that the Kozai–Lidov averaged step size is derived from first principles and demonstrated to reproduce the statistical evolution of numerical ensembles without fitting parameters, and (b) it provides a characteristic lifetime, instead of answering the binary question (stable/unstable), set by case-specific criteria.

Dwarfs in the Milky Way halo outer rim: first infall or backsplash satellites?

ABSTRACT Leo T is a gas-rich dwarf located at $414\, {\rm kpc}$ (1.4Rvir) distance from the Milky Way (MW) and it is currently assumed to be on its first approach. Here, we present an analysis of orbits calculated backwards in time for the dwarf with our new code delorean, exploring a range of systematic uncertainties, e.g. MW virial mass and accretion, M31 potential, and cosmic expansion. We discover that orbits with tangential velocities in the Galactic standard-of-rest frame lower than $| \vec{u}_{\rm t}^{\rm GSR}| \le 63^{+47}_{-39}\, {\rm km}\, {\rm s}^{\rm -1}$ result in backsplash solutions, i.e. orbits that entered and left the MW dark matter halo in the past, and that velocities above $| \vec{u}_{\rm t}^{\rm GSR}| \ge 21^{+33}_{-21}\, {\rm km}\, {\rm s}^{\rm -1}$ result in wide-orbit backsplash solutions with a minimum pericentre range of $D_{\rm min} \ge 38^{+26}_{-16}\, {\rm kpc}$, which would allow this satellite to survive gas stripping and tidal disruption. Moreover, new proper motion estimates overlap with our orbital solution regions. We applied our method to other distant MW satellites, finding a range of gas stripped backsplash solutions for the gasless Cetus and Eridanus II, providing a possible explanation for their lack of cold gas, while only first infall solutions are found for the H i-rich Phoenix I. We also find that the cosmic expansion can delay their first pericentre passage when compared to the non-expanding scenario. This study explores the provenance of these distant dwarfs and provides constraints on the environmental and internal processes that shaped their evolution and current properties.

Gigayear-timescale Destruction of High-eccentricity Asteroids by Spin and Why 2006 HY51 Has Been Spared

Abstract Asteroids and other small celestial bodies have markedly prolate shapes, and the perturbative triaxial torques that are applied during pericenter passages in highly eccentric orbits trigger and sustain a state of chaotic rotation. Because the prograde spin rate around the principal axis of inertia is not bounded from above, it can accidentally reach the threshold value corresponding to rotational breakup. Previous investigations of this process were limited to integrations of ∼103 orbits because of the stiff equation of motion. We present here a fast 1D simulation method to compute the evolution of this spin rate over ∼109 orbits. We apply the method to the most eccentric solar system asteroid known, 2006 HY51 (with e = 0.9684), and find that for any reasonably expected shape parameters it can never be accelerated to breakup speed. However, primordial solar system asteroids on more eccentric orbits may have already broken up from this type of rotational fission. The method also represents a promising opportunity to investigate the long-term evolution of extremely eccentric triaxial exo-asteroids (e > 0.99), which are thought to be common in white dwarf planetary systems.

The origin of the black hole offset in M31

ABSTRACT Using state-of-the-art high-resolution fully GPU N-body simulations, we demonstrate for the first time that the infall of a dark matter-rich satellite naturally explains a present black hole offset by subparsecs in M31. Observational data of the tidal features provide stringent constraints on the initial conditions of our simulations. The heating of the central region of M31 by the satellite via dynamical friction entails a significant black hole offset after the first pericentric passage. After having reached its maximum offset, the massive black hole sinks towards the M31 centre due to dynamical friction and it is determined to be offset by subparsecs as derived by observations.

Simulations of the Milky Way’s Central Molecular Zone – II. Star formation

ABSTRACTThe Milky Way’s Central Molecular Zone (CMZ) has emerged in recent years as a unique laboratory for the study of star formation. Here we use the simulations presented in Tress et al. to investigate star formation in the CMZ. These simulations resolve the structure of the interstellar medium at sub-parsec resolution while also including the large-scale flow in which the CMZ is embedded. Our main findings are as follows. (1) While most of the star formation happens in the CMZ ring at $R\gtrsim 100 \, {\rm pc}$, a significant amount also occurs closer to Sgr A* at $R \lesssim 10\, {\rm pc}$. (2) Most of the star formation in the CMZ happens downstream of the apocentres, consistent with the ‘pearls-on-a-string’ scenario, and in contrast to the notion that an absolute evolutionary timeline of star formation is triggered by pericentre passage. (3) Within the time-scale of our simulations (∼100 Myr), the depletion time of the CMZ is constant within a factor of ∼2. This suggests that variations in the star formation rate are primarily driven by variations in the mass of the CMZ, caused, for example, by active galactic nuclei (AGN) feedback or externally induced changes in the bar-driven inflow rate, and not by variations in the depletion time. (4) We study the trajectories of newly born stars in our simulations. We find several examples that have age and 3D velocity compatible with those of the Arches and Quintuplet clusters. Our simulations suggest that these prominent clusters originated near the collision sites where the bar-driven inflow accretes on to the CMZ, at symmetrical locations with respect to the Galactic Centre, and that they have already decoupled from the gas in which they were born.

The edge of the Galaxy

ABSTRACT We use cosmological simulations of isolated Milky Way (MW)-mass galaxies, as well as Local Group (LG) analogues, to define the ‘edge’ – a caustic manifested in a drop in density or radial velocity – of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ∼1.4r200m, corresponding to the ‘splashback’ radius, and a second caustic located at ∼0.6r200m, which likely corresponds to the edge of the virialized material that has completed at least two pericentric passages. The splashback radius is ill defined in LG-type environments where the haloes of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and is a more useful definition for the boundary of the MW halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase–space distribution of dwarf galaxies in the LG. Using the current dwarf galaxy population, we predict the edge of the MW halo to be 292 ± 61 kpc.

The fate of disc galaxies in IllustrisTNG clusters

ABSTRACTWe study the stellar morphological evolution of disc galaxies within clusters in the TNG50 and TNG100 runs from the IllustrisTNG simulation suite. We select satellites of masses 109.7 ≤ M*, z = 0/M⊙ ≤ 1011.6 residing in clusters of masses 1014 ≲ M200c, z = 0/M⊙ ≤ 1014.6 at z = 0 and that were discs at accretion according to a kinematic morphology indicator (the circularity fraction). These are traced from the time of accretion to z = 0 and compared to a control sample of central galaxies mass-matched at accretion. Most cluster discs become non-discy by z = 0, in stark contrast with the control discs, of which a significant fraction remains discy over the same timescales. Cluster discs become non-discy accompanied by gas removal and star formation quenching, loss of dark matter, and little growth or a loss of stellar mass. In contrast, control discs transform while also losing gas mass and quenching, but growing significantly in dark matter and stellar mass. Most cluster satellites change morphologies on similar timescales regardless of stellar mass, in ∼0.5–4 Gyr after accretion. Cluster discs that experienced more numerous and closer pericentric passages show the largest change in morphology. Morphological change in all cases requires the presence of a gravitational perturbation to drive stellar orbits to non-discy configurations, along with gas removal/heating to prevent replenishment of the disc through continued star formation. For cluster discs, the perturbation is impulsive tidal shocking at pericentres and not tidal stripping of outer disc stellar material, whereas for control discs, a combination of mergers and feedback from active galactic nuclei appears to be the key driving force behind morphological transformations.

Tidal evolution of galaxies in the most massive cluster of IllustrisTNG-100

We study the tidal evolution of galaxies in the most massive cluster of the IllustrisTNG-100 simulation. For the purpose of this work, we selected 112 galaxies with the largest stellar masses at present and followed their properties over time. Using their orbital history, we divided the sample into unevolved (infalling), weakly evolved (with one pericenter passage), and strongly evolved (with multiple pericenters). The samples are clearly separated by the value of the integrated tidal force from the cluster the galaxies experienced during their entire evolution and their properties depend strongly on this quantity. As a result of tidal stripping, the galaxies of the weakly evolved sample lost between 10 and 80% of their dark mass and less than 10% of stars, while those in the strongly evolved one lost more than 70% of dark mass and between 10 and 55% of stellar mass, and are significantly less, or not at all dark-matter dominated. While 33% of the infalling galaxies do not contain any gas, this fraction increases to 67% for the weakly evolved sample, and to 100% for the strongly evolved sample. The strongly evolved galaxies lose their gas earlier and faster (within 2–6 Gyr), but the process can take up to 4 Gyr from the first pericenter passage. These galaxies are redder and more metal rich, and at redshift z = 0.5, the population of galaxies in the cluster becomes predominantly red. As a result of tidal stirring, the morphology of the galaxies evolves from oblate to prolate and their rotation is diminished, thus the morphology-density relation is reproduced in the simulated cluster. The strongly evolved sample contains at least six convincing examples of tidally induced bars and six more galaxies that had their bars enhanced by their interaction with the cluster.

Dynamical self-friction: how mass loss slows you down

ABSTRACT We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some points in a subhalo’s orbit the torque of the stripped material can boost the orbital angular momentum of the remnant, the net effect over the long term is orbital decay regardless of the initial orbital parameters or subhalo mass. In order to quantify the strength of self-friction, we run a suite of simulations spanning typical host-to-subhalo mass ratios and orbital parameters. We find that the time-scale for self-friction, defined as the exponential decay time of the subhalo’s orbital angular momentum, scales with mass ratio and orbital circularity similar to standard dynamical friction. The decay time due to self-friction is roughly an order of magnitude longer, suggesting that self-friction only contributes at the 10 per cent level. However, along more radial orbits, self-friction can occasionally dominate over dynamical friction close to pericentric passage, where mass stripping is intense. This is also the epoch at which the self-friction torque undergoes large and rapid changes in both magnitude and direction, indicating that self-friction is an important process to consider when modelling pericentric passages of subhaloes and their associated satellite galaxies.

The recurrent impact of the Sagittarius dwarf on the star formation history of the Milky Way

Satellites orbiting disk galaxies can induce phase space features such as spirality, vertical heating and phase-mixing in their disks. Such features have also been observed in our own Galaxy, but the complexity of the Milky Way disk has only recently been fully mapped by Gaia Data Release 2 (DR2) data. This complex behaviour is mainly ascribed to repeated perturbations induced by the Sagittarius dwarf galaxy (Sgr) along its orbit, pointing to this satellite as the main dynamical architect of the Milky Way disk. Here, we model Gaia DR2-observed colour–magnitude diagrams to obtain a detailed star formation history of the ~2 kpc bubble around the Sun. It reveals three conspicuous and narrow episodes of enhanced star formation that we can precisely date as having occurred 5.7, 1.9 and 1.0 Gyr ago. The timing of these episodes coincides with proposed Sgr pericentre passages according to (1) orbit simulations, (2) phase space features in the Galactic disk and (3) Sgr stellar content. These findings most probably suggest that Sgr has also been an important actor in the build-up of the stellar mass of the Milky Way disk, with the perturbations from Sgr repeatedly triggering major episodes of star formation. Galactic close encounters can induce gravitational effects in the participants. Here Ruiz-Lara et al. have reconstructed the star formation history of the region of our Galaxy close to the Sun, finding that three recent visits of the neighbouring Sagittarius dwarf galaxy have resulted in well-defined episodes of star formation.

Simulations of gas sloshing induced by a newly discovered gas poor substructure in galaxy cluster Abell 1644

ABSTRACT Collision events lead to peculiar morphologies in the intracluster gas of galaxies clusters. That seems to be the case of Abell 1644, a nearby galaxy cluster, composed of three main structures: the southern cluster that exhibits a spiral-like morphology, A1644S; the northern cluster seen in X-ray observations, A1644N1; and the recently discovered substructure, A1644N2. By means of N-body hydrodynamical simulations, we attempt to reconstruct the dynamical history of this system. These simulations resulted in two specific scenarios: (i) The collision between A1644S and A1644N2. Our best model has an inclination between the merger plane and the plane of the sky of 30°, and reaches the best morphology 1.6 Gyr after the pericentric passage. At this instant A1644N2 is gas poor, becoming nearly undetectable in X-ray emission. This model shows a good agreement with observations; (ii) The collision between A1644S and A1644N1. This approach did not give rise to results as satisfactory as the first scenario, due to great disturbances in density and mismatching temperature maps. As a complementary study, we perform a three-cluster simulation using as base the best-fitting model to reproduce the current state of A1644 with the three main structures. This scenario presented a good agreement to the global morphology of the observations. Thus, we find that the more likely scenario is a collision between A1644S and the newly discovered A1644N2, where A1644N1 may be present as long as it does not greatly interfere in the formation of the spiral feature.

A Link between Ram Pressure Stripping and Active Galactic Nuclei

Abstract The dense environment of a galaxy cluster can radically transform the content of in-falling galaxies. Recent observations have found a significant population of active galactic nuclei (AGN) within “jellyfish galaxies,” galaxies with trailing tails of gas and stars that indicate significant ram pressure stripping. The relationship between AGN and ram pressure stripping is not well understood. In this Letter, we investigate the connection between AGN activity and ram pressure in a fully cosmological setting for the first time using the RomulusC simulation, one of the highest resolution simulations of a galaxy cluster to date. For lower mass galaxies (with stellar masses M * ≲ 109.5 M ⊙) both star formation and black hole accretion are suppressed by ram pressure before they reach pericenter, whereas for more massive galaxies accretion onto the black hole is enhanced during pericentric passage. Our analysis also indicates that as long as the galaxy retains gas, AGN with higher Eddington ratios are more likely to be found in galaxies experiencing higher ram pressure. We conclude that prior to quenching star formation, ram pressure triggers enhanced accretion onto the black hole, which in our model then produces heating and outflows due to AGN feedback. These simulations suggest that AGN feedback may in turn serve to aid in the quenching of star formation in tandem with ram pressure.

Refracting into ultra-diffuse galaxy NGC 1052-DF2 by passing near the centre of NGC 1052

ABSTRACT The recent observation of the ultra-diffuse galaxy NGC 1052-DF2 shows a galaxy may lack dark matter, which becomes a challenge to the standard structure formation theory. Here, we show that such phenomena can be a natural consequence if the NGC 1052-DF2 had experienced a single passage within a few kpc to the centre of the galaxy NGC 1052. The tidal effect of NGC 1052 in the encounter will heat the NGC 1052-DF2, and stretch the previous dwarf galaxy significantly into its current size. The relative lack of dark matter in the observed region is a natural consequence of the dark matter limited total amount in the corresponding small central region before the encounter, together with a systematic underestimation of the trace mass estimator method during revirialization after the encounter. The observed flat distribution of the ultra-diffuse galaxy can be reproduced with a cored initial star profile, which is a major improvement compared with the previous work. Our results show no need for introducing any new physical mechanism, as well as an alternative origin of an ultra-diffuse galaxy without repeated pericentre passage.

Detection of the Schwarzschild precession in the orbit of the star S2 near the Galactic centre massive black hole

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star’s radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this Letter we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2’s orbit. Owing to its highly elliptical orbit (e = 0.88), S2’s SP is mainly a kink between the pre-and post-pericentre directions of motion ≈±1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017−2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2’s position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states (“flares”) of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter fSP for the SP (fSP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, δϕ ≈ 12′ per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find fSP = 1.10 ± 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2’s orbit cannot exceed ≈0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M⊙.