Pericenter Distance Research Articles

Repeating Nuclear Transients from Repeating Partial Tidal Disruption Events: Reproducing ASASSN-14ko and AT2020vdq

Some electromagnetic outbursts from the nuclei of distant galaxies have been found to repeat on months-to-years timescales, and each of these sources can putatively arise from the accretion flares generated through the repeated tidal stripping of a star on a bound orbit about a supermassive black hole (SMBH), i.e., a repeating partial tidal disruption event (rpTDE). Here, we test the rpTDE model through analytical estimates and hydrodynamical simulations of the interaction between a range of stars, which differ from one another in mass and age, and an SMBH. We show that higher-mass (≳1M ⊙), evolved stars can survive many (≳10−100) encounters with an SMBH while simultaneously losing few × 0.01M ⊙, resulting in accretion flares that are approximately evenly spaced in time with nearly the same amplitude, quantitatively reproducing ASASSN-14ko. We also show that the energy imparted to the star via tides can lead to a change in its orbital period that is comparable to the observed decay in the recurrence time of ASASSN-14ko’s flares, Ṗ≃−0.0026 . Contrarily, lower-mass and less-evolved stars lose progressively more mass and produce brighter accretion flares on subsequent encounters for the same pericenter distances, leading to the rapid destruction of the star and cessation of flares. Such systems cannot reproduce ASASSN-14ko-like transients, but are promising candidates for recreating events such as AT2020vdq, which displayed a second and much brighter outburst compared to the first. Our results imply that the lightcurves of repeating transients are tightly coupled with stellar type.

Analytical Models for Secular Descents in Hierarchical Triple Systems

Three-body systems are prevalent in nature, from planetary to stellar to supermassive black hole scales. In a hierarchical triple system, oscillations of the inner orbit’s eccentricity and inclination can be induced on secular timescales. Over many cycles, the octupole-level terms in the secular equations of motion can drive the system to extremely high eccentricities via the eccentric Kozai–Lidov (EKL) mechanism. The overall decrease in the inner orbit’s pericenter distance has potentially dramatic effects for realistic systems, such as tidal disruption events. We present an analytical approximation in the test-particle limit to describe individual stepwise increases in eccentricity of the inner orbit. A second approximation, also in the test-particle limit, is obtained by integrating the equations of motion and calibrating to numerical simulations to estimate the overall octupole-level time evolution of the eccentricity. The latter approach is then extended beyond the test particle to the general case. The three novel analytical approximations are compared to numerical solutions to show that the models accurately describe the form and timescale of the secular descent from large distances to a close-encounter distance (e.g., the Roche limit). By circumventing the need for numerical simulations to obtain the long-term behavior, these approximations can be used to readily estimate properties of close encounters and descent timescales for populations of systems. We demonstrate this by calculating rates of EKL-driven migration for Hot Jupiters in stellar binaries.

Halo densities and pericenter distances of the bright Milky Way satellites as a test of dark matter physics

ABSTRACT We provide new constraints on the dark matter halo density profile of Milky Way (MW) dwarf spheroidal galaxies (dSphs) using the phase-space distribution function (DF) method. After assessing the systematics of the approach against mock data from the Gaia Challenge project, we apply the DF analysis to the entire kinematic sample of well-measured MW dwarf satellites for the first time. Contrary to previous findings for some of these objects, we find that the DF analysis yields results consistent with the standard Jeans analysis. In particular, in this study we rediscover (i) a large diversity in the inner halo densities of dSphs (bracketed by Draco and Fornax), and (ii) an anticorrelation between inner halo density and pericenter distance of the bright MW satellites. Regardless of the strength of the anticorrelation, we find that the distribution of these satellites in density versus pericenter space is inconsistent with the results of the high-resolution N-body simulations that include a disc potential. Our analysis motivates further studies on the role of internal feedback and dark matter microphysics in these dSphs.

Free-floating Planets, Survivor Planets, Captured Planets, and Binary Planets from Stellar Flybys

In star clusters, close stellar encounters can strongly impact the architecture of a planetary system or even destroy it. We present a systematic study of the effects of stellar flybys on two-planet systems. When such a system experiences flybys, one or both planets can be ejected, forming free-floating planets (FFPs), captured planets (CPs) around the flyby star, and free-floating binary planets (BPs); the remaining single-surviving planets (SSPs) can have their orbital radii and eccentricities greatly changed. Through numerical experiments, we calculate the formation fractions (or branching ratios) of FFPs, SSPs, CPs, and BPs as a function of the pericenter distance of the flyby, and use them to derive analytical expressions for the formation rates of FFPs, SSPs, CPs and BPs in general cluster environments. We find that the production rates of FFPs and SSPs are similar (for the initial planet semimajor axis ratio a 1/a 2 = 0.6–0.8), while the rate for CPs is a few times smaller. The formation fraction of BPs depends strongly on a 1/a 2 and on the planet masses. For Jupiter-mass planets, the formation fraction of BPs is always less than 1% (for a 1/a 2 = 0.8) and typically much smaller (≲0.2% for a 1/a 2 ≲ 0.7). The fraction remains less than 1% when considering 4M J planets. Overall, when averaging over all flybys, the production rate of BPs is less than 0.1% of that of FFPs. We also derive the velocity distribution of FFPs produced by stellar flybys, and the orbital parameter distributions of SSPs, CPs, and BPs. These results can be used in future studies of exotic planets (including FFPs) and planetary systems.

Анализ регулярности/хаотичности динамики шаровых скоплений в центральной области Млечного Пути

We analyzed the regularity/chaoticity of the orbits of 45 globular clusters in the central region of the Galaxy with a radius of 3.5 kpc, which are subject to the greatest influence from an elongated rotating bar. Various analysis methods were used, namely, methods for calculating the maximum characteristic Lyapunov exponents, the method of Poincaré sections, the frequency method based on the calculation of fundamental frequencies, as well as the visual assessment method. Bimodality was discovered in the histogram of the distribution of positive Lyapunov exponents, calculated in the classical version, without renormalization of the shadow orbit, which makes it possible to implement a probabilistic method of GC classification. To construct the orbits of globular clusters, we used a gravitational potential model with a bar in the form of a triaxial ellipsoid, described in detail in the work of Bajkova et al., Izvestiya GAO in Pulkovo, 2023, 228, 1. The following bar parameters were adopted: mass 1010M⊙, semimajor axis length 5 kpc, bar viewing angle 25o, rotation speed 40 km s−1 kpc−1. To form the 6D-phase space required for orbit integration, the most accurate astrometric data to date from the Gaia satellite (EDR3) (Vasiliev, Baumgardt, 2021), as well as new refined average distances to globular clusters (Baumgardt, Vasiliev, 2021) were used. A classification is made of globular clusters with regular and chaotic dynamics. As the analysis showed, globular clusters with small pericentric distances and large eccentricities are most susceptible to the influence of the bar and demonstrate the greatest chaos. It is shown that the results of classifying globular clusters by the nature of their orbital dynamics, obtained using the various methods of analysis considered in the work, correlate well with each other.

The Orbital Structure and Selection Effects of the Galactic Center S-star Cluster

The orbital distribution of the S-star cluster surrounding the supermassive black hole in the center of the Milky Way is analyzed. A tight dependence of the pericenter distance r p on orbital eccentricity e ⋆ is found, log(rp)∼(1−e⋆) , which cannot be explained simply by a random distribution of semimajor axis and eccentricities. No stars are found in the region with high e ⋆ and large log(rp) or in the region with low e ⋆ and small log(rp) . Although the sample is still small, the G-clouds show a very similar distribution. The likelihood P(log(rp),(1−e⋆)) to determine the orbital parameters of S-stars is determined. P is very small for stars with large e ⋆ and large log(rp) . S-stars might exist in this region. To determine their orbital parameters, one however needs observations over a longer time period. On the other hand, if stars would exist in the region of low log(rp) and small e ⋆, their orbital parameters should by now have been determined. That this region is unpopulated therefore indicates that no S-stars exist with these orbital characteristics, providing constraints for their formation. We call this region, defined by log(rp/AU)<1.57+2.6(1−e⋆) , the zone of avoidance. Finally, it is shown that the observed frequency of eccentricities and pericenter distances is consistent with a random sampling of log(rp) and e ⋆ if one takes into account the fact that no stars exist in the zone of avoidance and that orbital parameters cannot yet be determined for stars with large r p and large e ⋆.

Modelling the orbital histories of satellites of Milky Way-mass galaxies: testing static host potentials against cosmological simulations

ABSTRACT Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modelling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass haloes. When a typical surviving satellite fell in ($3.4\!-\!9.7\, \rm {Gyr}$ ago), the host halo mass and radius were typically 26–86 per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $r\lesssim 50\, \rm {kpc}$, opposite to dark matter only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at z = 0 and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $1.6\!-\!4.7\, \rm {Gyr}$ ago. Most orbital properties are minimally biased, ≲10 per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are ≈20 per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are ≈40–80 per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at z = 0.

Scattered disc dynamics: the mapping approach

ABSTRACT We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericentre of the parabolic orbit. The map shows excellent agreement with direct numerical integrations and elucidates how the dynamics depends on perturber mass and pericentre distance. We also use the map to explore the onset of chaos, statistical descriptions of chaotic transport, and sticking in mean-motion resonances. We discuss implications of our mapping model for the dynamical evolution of the Solar system’s scattered disc and other highly eccentric trans-Neptunian objects.

Tracking the evolution of satellite galaxies: mass stripping and dark-matter deficient galaxies

ABSTRACT Satellite galaxies undergo a variety of physical processes when they are accreted by groups and clusters, often resulting in the loss of baryonic and dark matter (DM) mass. In this work, we evaluate the predictions from the IllustrisTNG hydrodynamical simulation regarding the evolution of the matter content of satellites, focusing on a population that are accreted at z &amp;gt; 1 and retain their identity as satellites down to z = 0. At fixed host halo mass, the amount of DM and stellar mass stripped depends mostly on the pericentric distance, dperi, here normalized by host halo virial radius. The closest encounters result in significant loss of DM, with subhaloes retaining between 20 and a few per cent of their z = 1 mass. At fixed dperi, DM mass stripping seems more severe in lower mass haloes. Conversely, the average satellite in higher mass haloes has its stellar mass growth halted earlier, having lost a higher fraction of stellar mass by z = 0. We also show that mass stripping has a strong impact on the quenched fractions. The IllustrisTNG boxes are qualitatively consistent in these predictions, with quantitative differences mostly originating from the distinct subhalo mass ranges covered by the boxes. Finally, we have identified DM-deficient systems in all TNG boxes. These objects are preferentially found in massive clusters (Mhost ≳ 1013 M⊙), had very close encounters with their central galaxies ($d_{\rm peri}\simeq 0.05\, R_{\rm vir}$), and were accreted at high redshift (zinfall ≳ 1.4), reinforcing the notion that tidal stripping is responsible for their remarkable lack of DM.

Extreme Tidal Stripping May Explain the Overmassive Black Hole in Leo I: A Proof of Concept

A recent study found dynamical evidence of a supermassive black hole of ∼3 × 106 M ⊙ at the center of Leo I, the most distant dwarf spheroidal galaxy of the Milky Way. This black hole, comparable in mass to the Milky Way’s Sgr A*, places the system >2 orders of magnitude above the standard M •–M ⋆ relation. We investigate the possibility, from a dynamical standpoint, that Leo I’s stellar system was originally much more massive and, thus, closer to the relation. Extreme tidal disruption from one or two close passages within the Milky Way’s virial radius could have removed most of its stellar mass. A simple analytical model suggests that the progenitor of Leo I could have experienced a mass loss in the range 32%–57% from a single pericenter passage, depending on the stellar velocity dispersion estimate. This mass-loss percentage increases to the range 66%–78% if the pericenter occurs at the minimum distance allowed by current orbital reconstructions. Detailed N-body simulations show that the mass loss could reach ∼90% with up to two passages, again with pericenter distances compatible with the minimum value allowed by Gaia data. Despite very significant uncertainties in the properties of Leo I, we reproduce its current position and velocity dispersion, as well as the final stellar mass enclosed in 1 kpc (∼5 × 106 M ⊙) within a factor <2. The most recent tidal stream is directed along our line of sight toward Leo I, making it difficult to detect. Evidence from this extreme tidal disruption event could be present in current Gaia data in the form of extended tidal streams.

Dark Matter Halo Properties of the Galactic Dwarf Satellites: Implication for Chemo-dynamical Evolution of the Satellites and a Challenge to Lambda Cold Dark Matter

Elucidating dark matter density profiles in Galactic dwarf satellites is essential to understanding not only the quintessence of dark matter, but also the evolution of the satellites themselves. In this work, we present the current constraints on dark matter densities in Galactic ultrafaint dwarf (UFD) and diffuse galaxies. Applying our constructed nonspherical mass models to the currently available kinematic data of the 25 UFDs and two diffuse satellites, we find that whereas most of the galaxies have huge uncertainties on the inferred dark matter density profiles, Eridanus II, Segue I, and Willman 1 favor cuspy central profiles even when considering effects of a prior bias. We compare our results with the simulated subhalos on the plane between the dark matter density at 150 pc and the pericenter distance. We find that the most observed satellites and the simulated subhalos are similarly distributed on this plane, except for Antlia 2, Crater 2, and Tucana 3, which are less than one-tenth of the density. Despite considerable tidal effects, the subhalos detected by commonly used subhalo finders have difficulty explaining such a huge deviation. We also estimate the dynamical mass-to-light ratios of the satellites and confirm the ratio is linked to stellar mass and metallicity. Tucana 3 deviates largely from these relations, while it follows the mass–metallicity relation. This indicates that Tucana 3 has a cored dark matter halo, despite a significant uncertainty in its ratios.

Backsplash galaxies and their impact on galaxy evolution: a three-stage, four-type perspective

ABSTRACT We study the population of backsplash galaxies at z = 0 in the outskirts of massive, isolated clusters of galaxies taken from the mdpl2-sag semi-analytical catalogue. We consider four types of backsplash galaxies according to whether they are forming stars or passive at three stages in their lifetimes: before entering the cluster, during their first incursion through the cluster, and after they exit the cluster. We analyse several geometric, dynamic, and astrophysical aspects of the four types at the three stages. Galaxies that form stars at all stages account for the majority of the backsplash population (58 per cent) and have stellar masses typically below $M_\star \sim 3\times 10^{10} \ h^{-1}\, {\rm {M_\odot }}$ that avoid the innermost cluster’s regions and are only mildly affected by it. In a similar mass range, galaxies that become passive after exiting the cluster (26 per cent) follow orbits characterized by small pericentric distance and a strong deflection by the cluster potential well while suffering a strong loss of both dark matter and gas content. Only a small fraction of our sample (4 per cent) becomes passive while orbiting inside the cluster. These galaxies have experienced heavy pre-processing and the cluster’s tidal stripping and ram pressure provide the final blow to their star formation. Finally, galaxies that are passive before entering the cluster for the first time (12 per cent) are typically massive and are not affected significantly by the cluster. Using the bulge/total mass ratio as a proxy for morphology, we find that a single incursion through a cluster does not result in significant morphological changes in all four types.

Шаровые скопления в центральной области Млечного Пути I. Влияние бара на параметры орбит по данным Gaia EDR3

The work is devoted to the analysis of the influence of the galactic bar on the orbital motion of globular clusters in the central region of the Galaxy. For this task, 45 globular clusters were selected, 34 of which belong to the bulge and 11 to the disk. The most accurate astrometric data from the Gaia satellite (Vasiliev and Baumgardt, 2021), as well as new refined average distances (Baumgardt and Vasiliev, 2021), were used to form the 6D phase space required for orbit integration. The orbits of globular clusters are obtained both in an axisymmetric potential and in a potential including a bar. In this case, the mass, rotation speed, shape and scale length of the bar were varied. A comparison is made of such orbital parameters as apocentric and pericentric distances, eccentricity and maximum distance from the galactic plane. It is shown that the mass of the bar exerts the greatest influence on the orbital motion, which is expressed mainly in an increase in both the apocentric and pericentric distances in the vast majority of globular clusters. The eccentricities of the orbits in the overwhelming majority also change significantly, and there is a change both upward and downward, especially in the range of values from 0.2 to 0.8. The greatest changes in parameters are observed in globular clusters with high radial velocities and small pericentric distances. The change in orbital parameters depending on the bar rotation speed is less pronounced. The influence of the geometric parameters of the bar is insignificant in the accepted range of their changes. Several examples show that globular clusters in the bulge are more affected by the bar than those belonging to the disk.

Merging galaxy clusters in IllustrisTNG

Mergers between galaxy clusters are an important stage in the formation of the large-scale structure of the Universe. Some of the mergers show a spectacular bow shock that formed as a result of recent passage of a smaller cluster through a bigger one, the classic example of this being the so-called bullet cluster. In this paper, I describe ten examples of interacting clusters identified among 200 of the most massive objects, with total masses above 1.4 × 1014 M⊙, from the IllustrisTNG300 simulation by searching for prominent bow shocks in their temperature maps. Despite different mass ratios of the two merging clusters, the events are remarkably similar in many respects. In all cases, the companion cluster passed close to the main one only once, between 0.9 and 0.3 Gyr ago, with the pericenter distance of 100–530 kpc and a velocity of up to 3400 km s−1. The subcluster, typically an order of magnitude smaller in mass than the main cluster before the interaction, loses most of its dark matter and gas in the process. The displacement between the collisionless part of the remnant and the bow shock is such that the remnant typically lags behind the shock or coincides with it, with a single exception of the merger occurring with the largest velocity. Usually about 1% of the gas cells in the merging clusters are shocked, and the median Mach numbers of these gas cells are around two. Due to the relatively small size of the simulation box, no close analog of the bullet cluster was found, but I identified one case that is similar in terms of mass, velocity, and displacement. The presented cases bear more resemblance to less extreme observed interacting clusters such as A520 and Coma.

Deviation in stellar trajectory induced by asymmetry in partial tidal disruption

ABSTRACT We study partial tidal disruption and present a quantitative analysis of the orbital dynamics of the remnant self-bound core. We perform smoothed particle hydrodynamical simulations to show that partial disruption of a star due to the tidal field of a black hole leads to a jump in the specific orbital energy and angular momentum of the core. It directly leads to deviation in the core’s trajectory apart from getting a boost in its velocity. Our analysis shows that the variations in the specific orbital energy and angular momentum are higher when the pericenter distance is lower. We conclude that higher mass asymmetry of the two tidal tails increases the magnitude of the trajectory deviations. We show that observable deviations are only possible when mass ratio $q \lesssim \mathcal {O}(10^3)$ , which indicates the range of intermediate-mass black holes.

On the anticorrelation between pericentric distance and inner dark matter density of Milky Way’s dwarf spheroidal galaxies

ABSTRACT An anticorrelation between the central density of the dark matter (DM) halo (ρ150, DM) and the pericentric distances (rp) of the Milky Way’s (MW’s) dwarf spheroidal galaxies (dSphs) has been reported in the literature. The existence and origin of such anticorrelation are, however, controversial, one possibility being that only the densest dSphs can survive the tidal field towards the centre of our Galaxy. In this work, we place particular emphasis on quantifying the statistical significance of such anticorrelation, by using available literature data in order to explore its robustness under different assumptions on the MW gravitational potential, and for various derivations of ρ150 and rp. We consider models in which the MW is isolated and has low ($8.8\times 10^{11}\, {\rm M}_{\odot }$ ) and high ($1.6\times 10^{12}\, {\rm M}_{\odot }$ ) halo masses, respectively, as well as configurations in which the MW’s potential is perturbed by a Large Magellanic Cloud (LMC) infall. We find that, while data generally support models in which the dSphs’ central DM density decreases as a function of their pericentric radius, this anticorrelation is statistically significant at 3σ level only in ${\sim} 12~{{ \rm per\ cent}}$ of the combinations of ρ150 and rp explored. Moreover, including the impact of the LMC’s infall on to the MW weakens or even washes away this anticorrelation, with respect to models in which the MW is isolated. Our results suggest that the strength and existence of such anticorrelation are still debatable: exploring it with high-resolution simulations including baryonic physics and different DM flavours will help us to understand its emergence.

The Mass Fallback Rate of the Debris in Relativistic Stellar Tidal Disruption Events

Highly energetic stellar tidal disruption events (TDEs) provide a way to study black hole characteristics and their environment. We simulate TDEs with the Phantom code in a general relativistic and Newtonian description of a supermassive black hole’s gravity. Stars, which are placed in parabolic orbits with different β parameters, are constructed with the stellar evolution code MESA and therefore have realistic stellar density profiles. We study the mass fallback rate of the debris and its dependence on β, stellar mass, and age as well as the supermassive black hole’s spin and the choice of the gravity description. We calculate the peak value , time to peak t peak, duration of the super-Eddington phase t Edd, time during which , early rise-time τ rise, and late-time slope n ∞. We recover the trends of , t peak, τ rise, and n ∞ with β, stellar mass, and age which were obtained in previous studies. We find that t Edd, at a fixed β, scales primarily with the stellar mass, while scales with the compactness of stars. The effect of the SMBH’s rotation depends on the orientation of its rotational axis relative to the direction of the stellar motion in the initial orbit. Encounters in prograde orbits result in narrower curves with higher , while the opposite occurs for retrograde orbits. We find that disruptions, at the same pericenter distance, are stronger in a relativistic tidal field than in a Newtonian one. Therefore, relativistic curves have higher , and shorter t peak and t Edd.

Orbital evolution of satellite galaxies in self-interacting dark matter models

Dark matter self-interactions can leave distinctive signatures on the properties of satellite galaxies around Milky Way--like hosts through their impact on tidal stripping, ram pressure, and gravothermal collapse. We delineate the regions of self-interacting dark matter parameter space---specified by interaction cross section and a velocity scale---where each of these effects dominates and show how the relative mass loss depends on the satellite's initial mass, density profile, and orbit. We obtain novel, conservative constraints in this parameter space using Milky Way satellite galaxies with notably high central densities and small pericenter distances. Our results for self-interacting dark matter models, in combination with constraints from clusters of galaxies, favor either velocity-dependent cross sections that lead to gravothermal core collapse in the densest satellites or small cross sections which more closely resemble cold and collisionless dark matter.

Stellar ages, masses, extinctions, and orbital parameters based on spectroscopic parameters ofGaiaDR3

Context.Gaia’s third data release provides radial velocities for 33 million stars and spectroscopically derived atmospheric parameters for more than 5 million targets. When combined with the astrometric data, these allow us to derive orbital and stellar parameters that are key to understanding the stellar populations of the Milky Way and to perform Galactic archaeology.Aims.We used the calibrated atmospheric parameters, 2MASS andGaia-EDR3 photometry, and parallax-based distances to compute the ages, initial stellar masses, and reddenings for the stars with spectroscopic parameters. We also derived the orbits for all of the stars with measured radial velocities and astrometry, adopting two sets of line-of-sight distances from the literature.Methods.Four different sets of ages, masses, and absolute magnitudes in different photometric bands are obtained through an isochrone fitting method, considering different combinations of input parameters. The reddenings are obtained by comparing the observed colours with those obtained from the isochrone projection. Finally, the orbits are computed adopting an axisymmetric potential of the Galaxy.Results.Comparisons with reference catalogues of field and cluster stars suggest that reliable ages are obtained for stars younger than 9–10 Gyr when the estimated relative age uncertainty is &lt; 50 per cent. For older stars, ages tend to be underestimated. The most reliable stellar type for age determination are turn-off stars, even when the input atmospheric parameters have large uncertainties. Ages for giants and main-sequence stars are retrieved with uncertainties of the order of 2 Gyr when extinction towards the star’s line of sight is less thanAV ≲ 2.5 mag.Conclusions.The catalogue of ages, initial stellar masses, reddenings, galactocentric positions and velocities, as well as the stellar orbital actions, eccentricities, apocentre, pericentre and maximum distance from the Galactic plane reached during their orbits, is made publicly available to be downloaded. With this catalogue, the full chemo-dynamical properties of the extended solar neighbourhood unfold and allow us to better identify the properties of the spiral arms, to parametrise the dynamical heating of the disc, and to thoroughly study the chemical enrichment of the Milky Way.

Planes of Satellites around Simulated Disk Galaxies. II. Time-persistent Planes of Kinematically Coherent Satellites in ΛCDM

We use two zoom-in ΛCDM hydrodynamical simulations of massive disk galaxies to study the possible existence of fixed satellite groups showing a kinematically coherent behavior across evolution (angular momentum conservation and clustering). We identify three such groups in the two simulations, defining kinematically coherent persistent planes (KPPs) that last at least from virialization to z = 0 (more than 7 Gyr). This proves that orbital pole clustering is not necessarily set in at low redshift, representing a long-lived property of galaxy systems. KPPs are thin and oblate, represent ∼25%–40% of the total number of satellites in the system, and are roughly perpendicular to their corresponding central disk galaxies during certain periods, consistently with Milky Way z = 0 data. KPP satellite members are statistically distinguishable from satellites outside KPPs: they show higher specific orbital angular momenta, orbit more perpendicularly to the central disk galaxy, and have larger pericentric distances than the latter. We numerically prove, for the first time, that KPPs and the best-quality positional planes share the same space configuration across time, such that KPPs act as “skeletons” preventing the latter from being washed out in short timescales. In one of the satellite−host systems, we witness the late capture of a massive dwarf galaxy endowed with its own satellite system, also organized into a KPP configuration prior to its capture. We briefly explore the consequences this event has on the host’s KPP and on the possible enhancement of the asymmetry in the number of satellites rotating in one sense or the opposite within the KPP.