We study chaos and L\'evy flights in the general gravitational three-body problem. We introduce new metrics to characterize the time evolution and final lifetime distributions, namely Scramble Density $\mathcal{S}$ and the LF index $\mathcal{L}$, that are derived from the Agekyan-Anosova maps and homology radius $R_{\mathcal{H}}$. Based on these metrics, we develop detailed procedures to isolate the ergodic interactions and L\'evy flight interactions. This enables us to study the three-body lifetime distribution in more detail by decomposing it into the individual distributions from the different kinds of interactions. We observe that ergodic interactions follow an exponential decay distribution similar to that of radioactive decay. Meanwhile, L\'evy flight interactions follow a power-law distribution. L\'evy flights in fact dominate the tail of the general three-body lifetime distribution, providing conclusive evidence for the speculated connection between power-law tails and L\'evy flight interactions. We propose a new physically-motivated model for the lifetime distribution of three-body systems and discuss how it can be used to extract information about the underlying ergodic and L\'evy flight interactions. We discuss mass ejection probabilities in three-body systems in the ergodic limit and compare it to previous ergodic formalisms. We introduce a novel mechanism for a three-body relaxation process and discuss its relevance in general three-body systems.

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

We report the first discovery of a thick-disk planet, LHS 1815b (TOI-704b, TIC 260004324), detected in the TESS survey. LHS 1815b transits a bright (V = 12.19 mag, K = 7.99 mag) and quiet M dwarf located $ 29.87\pm0.02 pc$ away with a mass of $0.502\pm0.015 M_{\odot}$ and a radius of $0.501\pm0.030 R_{\odot}$. We validate the planet by combining space and ground-based photometry, spectroscopy, and imaging. The planet has a radius of $1.088\pm 0.064 R_{\oplus}$ with a $3 \sigma$ mass upper-limit of $8.7 M_{\oplus}$. We analyze the galactic kinematics and orbit of the host star LHS1815 and find that it has a large probability ($P_{thick}/P_{thin} = 6482$) to be in the thick disk with a much higher expected maximal height ($Z_{max} = 1.8 kpc$) above the Galactic plane compared with other TESS planet host stars. Future studies of the interior structure and atmospheric properties of planets in such systems using for example the upcoming James Webb Space Telescope (JWST), can investigate the differences in formation efficiency and evolution for planetary systems between different Galactic components (thick and thin disks, and halo).

We investigate the role of the Keplerian tidal field generated by a supermassive black hole (SMBH) on the three-body dynamics of stellar mass black holes. We consider two scenarios occurring close to the SMBH: the breakup of unstable triples and three-body encounters between a binary and a single. These two cases correspond to the hard and soft binary cases, respectively. The tidal field affects the breakup of triples by tidally limiting the system, so that the triples break earlier with lower breakup velocity, leaving behind slightly larger binaries (relative to the isolated case). The breakup direction becomes anisotropic and tends to follow the shape of the Hill region of the triple, favouring breakups in the radial direction. Furthermore, the tidal field can torque the system, leading to angular momentum exchanges between the triple and its orbit about the SMBH. This process changes the properties of the final binary, depending on the initial angular momentum of the triple. Finally, the tidal field also affects binary-single encounters: binaries tend to become both harder and more eccentric with respect to encounters that occur in isolation. Consequently, single-binary scattering in a deep Keplerian potential produces binaries with shorter gravitational wave merger timescales.

High-density cusps of compact remnants are expected to form around supermassive black holes (SMBHs) in galactic nuclei via dynamical friction and two-body relaxation. Due to the high density, binaries in orbit around the SMBH can frequently undergo close encounters with compact remnants from the cusp. This can affect the gravitational wave merger rate of compact binaries in galactic nuclei. We investigated this process by means of high accuracy few-body simulations, performed with a novel Monte Carlo approach. We find that, around a SgrA*-like SMBH, three-body encounters increase the number of mergers by a factor of 3. This occurs because close encounters can reorient binaries with respect to their orbital plane around the SMBH, increasing the number of Kozai-Lidov induced mergers. We obtain a binary black hole merger rate of 1.6x10^-6/yr per Milky Way-like nucleus.

The origin of the molecular clouds orbiting SgrA* is one of the most debated questions about our Galactic Centre. Here, we present the hydrodynamic simulation of a molecular cloud infalling towards SgrA*, performed with the adaptive-mesh-refinement code RAMSES. Through such simulation, we propose that the circumnuclear disc and the +20 km/s cloud originated from the same tidal disruption episode, occurred less than 1 Myr ago. We also show that recent star formation is to be expected in the +20 km/s cloud, as also suggested by recent observations.

In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs.

Studying the formation and evolution of black hole binaries (BHBs) is essential for the interpretation of current and forthcoming gravitational wave (GW) detections. We investigate the statistics of BHBs that form from isolated binaries, by means of a new version of the SEVN population-synthesis code. SEVN integrates stellar evolution by interpolation over a grid of stellar evolution tracks. We upgraded SEVN to include binary stellar evolution processes and we used it to evolve a sample of $1.5\times{}10^8$ binary systems, with metallicity in the range $\left[10^{-4};4\times 10^{-2}\right]$. From our simulations, we find that the mass distribution of black holes (BHs) in double compact-object binaries is remarkably similar to the one obtained considering only single stellar evolution. The maximum BH mass we obtain is $\sim 30$, $45$ and $55\, \mathrm{M}_\odot$ at metallicity $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. A few massive single BHs may also form ($\lesssim 0.1\%$ of the total number of BHs), with mass up to $\sim 65$, $90$ and $145\, \mathrm{M}_\odot$ at $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. These BHs fall in the mass gap predicted from pair-instability supernovae. We also show that the most massive BHBs are unlikely to merge within a Hubble time. In our simulations, merging BHs like GW151226 and GW170608, form at all metallicities, the high-mass systems (like GW150914, GW170814 and GW170104) originate from metal poor ($Z\lesssim{}6\times 10^{-3}$) progenitors, whereas GW170729-like systems are hard to form, even at $Z = 10^{-4}$. The BHB merger rate in the local Universe obtained from our simulations is $\sim 90 \mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$, consistent with the rate inferred from LIGO-Virgo data.

Recent spectroscopic analysis has set an upper limit to the age of the S-stars, the ~30 B-type stars in highly eccentric orbits around the supermassive black hole (SMBH) in the Galactic center. The inferred age (<15 Myr) is in tension with the binary break-up scenario proposed to explain their origin. However, the new estimate is compatible with the age of the disk of O-type stars that lies at a farther distance from the SMBH. Here we investigate a new formation scenario, assuming that both S-stars and the O-type stars were born in the same disk around SgrA*. We simulate encounters between binaries of the stellar disk and stellar black holes from a dark cusp around SgrA*. We find that B-type binaries can be easily broken up by the encounters and their binary components are kicked into highly eccentric orbits around the SMBH. In contrast, O-type binaries are less frequently disrupted and their members remain in low eccentricity orbits. This mechanism can reproduce 12 S-stars just by assuming that the binaries initially lie within the stellar disk as observed nowadays. To reproduce all the S-stars, the original disk must have been extended down to 0.006 pc. However in this case many B- and O-type stars remain in low eccentricity orbits below 0.03 pc, in contrast with the observations. Therefore, some other mechanism is necessary to disrupt the disk below 0.03 pc. This scenario can also explain the high eccentricity of the G-objects, if they have a stellar origin.

Planet host stars with well-constrained ages provide a rare window to the time domain of planet formation and evolution. The NASA K2 mission has enabled the discovery of the vast majority of known planets transiting stars in clusters, providing a valuable sample of planets with known ages and radii. We present the discovery of two planets transiting K2-264, an M2 dwarf in the intermediate age (600-800 Myr) Praesepe open cluster (also known as the Beehive Cluster, M44, or NGC 2632), which was observed by K2 during Campaign 16. The planets have orbital periods of 5.8 and 19.7 days, and radii of $2.2 \pm 0.2 $ and $2.7 \pm 0.2$ $R_\oplus$, respectively, and their equilibrium temperatures are $496 \pm 10$ and $331 \pm 7$ $K$, making this a system of two warm sub-Neptunes. When placed in the context of known planets orbiting field stars of similar mass to K2-264, these planets do not appear to have significantly inflated radii, as has previously been noted for some cluster planets. As the second known system of multiple planets transiting a star in a cluster, K2-264 should be valuable for testing theories of photoevaporation in systems of multiple planets. Follow-up observations with current near-infrared (NIR) spectrographs could yield planet mass measurements, which would provide information about the mean densities and compositions of small planets soon after photoevaporation is expected to have finished. Follow-up NIR transit observations using Spitzer or large ground-based telescopes could yield improved radius estimates, further enhancing the characterization of these interesting planets.

We present a drag force model for evolving chaotic few-body interactions with the inclusion of orbital energy losses, such as tidal dissipation and gravitational wave (GW) emission. The main effect from such losses is the formation of two-body captures, that for compact objects result in GW mergers, and for stars lead to either compact binaries, mergers or disruptions. Studying the inclusion of energy loss terms in few-body interactions is therefore likely to be important for modeling and understanding the variety of transients that soon will be observed by current and upcoming surveys. However, including especially tides in few-body codes has been shown to be technically difficult and computationally heavy, which has lead to very few systematic tidal studies. In this paper we derive a drag force term that can be used to model the effects from tidal, as well as other, energy losses in few-body interactions, if the two-body orbit averaged energy loss is known a priori. This drag force model is very fast to evolve, and gives results in agreement with other approaches, including the impulsive and affine tide approximations.

We investigate the formation of circumnuclear gas structures from the tidal disruption of molecular clouds in galactic nuclei, by means of smoothed particle hydrodynamics simulations. We model galactic nuclei as composed of a supermassive black hole (SMBH) and a nuclear star cluster (NSC) and consider different mass ratios between the two components. We find that the relative masses of the SMBH and the NSC have a deep impact on the morphology of the circumnuclear gas. Extended disks form only inside the sphere of influence of the SMBH. In contrast, compact rings naturally form outside the SMBH's sphere of influence, where the gravity is dominated by the NSC. This result is in agreement with the properties of the Milky Way's circumnuclear ring, which orbits outside the SMBH sphere of influence. Our results indicate that compact circumnuclear rings can naturally form outside the SMBH sphere of influence.

Recent observations suggest ongoing planet formation in the innermost parsec of the Galactic center (GC). The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. Photoevaporation by the ultraviolet field of young stars, combined with ongoing tidal disruption, could enhance the near-infrared luminosity of such starless planets, making their detection possible even with current facilities. In this paper, we investigate the chance of planet tidal captures by means of high-accuracy N-body simulations exploiting Mikkola's algorithmic regularization. We consider both planets lying in the clockwise (CW) disk and planets initially bound to the S-stars. We show that tidally captured planets remain on orbits close to those of their parent star. Moreover, the semi-major axis of the planet orbit can be predicted by simple analytic assumptions in the case of prograde orbits. We find that starless planets that were initially bound to CW disk stars have mild eccentricities and tend to remain in the CW disk. However, we speculate that angular momentum diffusion and scattering with other young stars in the CW disk might bring starless planets on low-angular momentum orbits. In contrast, planets initially bound to S-stars are captured by the SMBH on highly eccentric orbits, matching the orbital properties of the G1 and G2 clouds. Our predictions apply not only to planets but also to low-mass stars initially bound to the S-stars and tidally captured by the SMBH.

Recent observations suggest ongoing planet formation in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. We investigate the chance of planet tidal captures by running three-body encounters of SMBH-star-planet systems with a high-accuracy regularized code. We show that tidally captured planets have orbits close to those of their parent star. We conclude that the final periapsis distance of the captured planet from the SMBH will be much larger than 200 AU, unless its parent star was already on a highly eccentric orbit.

The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermost parsec. We investigate the perturbations exerted by a gas ring on a nearly-Keplerian stellar disk orbiting a SMBH by means of combined direct N-body and smoothed particle hydrodynamics simulations. We simulate the formation of gas rings through the infall and disruption of a molecular gas cloud, adopting different inclinations between the infalling gas cloud and the stellar disk. We find that a CNR-like ring is not efficient in affecting the stellar disk on a timescale of 3 Myr. In contrast, a gas ring in the innermost 0.5 pc induces precession of the longitude of the ascending node Omega, significantly affecting the stellar disk inclination. Furthermore, the combined effect of two-body relaxation and Omega-precession drives the stellar disk dismembering, displacing the stars from the disk. The impact of precession on the star orbits is stronger when the stellar disk and the inner gas ring are nearly coplanar. We speculate that the warm gas in the inner cavity might have played a major role in the evolution of the clockwise disk.

Several thousand solar masses of molecular, atomic and ionized gas lie in the innermost ~10 pc of our Galaxy. The most relevant structure of molecular gas is the circumnuclear ring (CNR), a dense and clumpy ring surrounding the supermassive black hole (SMBH), with a radius of ~2 pc. We propose that the CNR formed through the tidal disruption of a molecular cloud, and we investigate this scenario by means of N-body smoothed-particle hydrodynamics simulations. We ran a grid of simulations with different cloud mass (4X10^4, 1.3X10^5 solar masses), different initial orbital velocity (v_in=0.2-0.5 v_esc, where v_esc is the escape velocity from the SMBH), and different impact parameter (b=8, 26 pc). The disruption of the molecular cloud leads to the formation of very dense and clumpy gas rings, containing most of the initial cloud mass. If the initial orbital velocity of the cloud is sufficiently low (v_in<0.4 v_esc, for b=26 pc) or the impact parameter is sufficiently small (b<10 pc, for v_in>0.5 v_esc), at least two rings form around the SMBH: an inner ring (with radius ~0.4 pc) and an outer ring (with radius ~2-4 pc). The inner ring forms from low-angular momentum material that engulfs the SMBH during the first periapsis passage, while the outer ring forms later, during the subsequent periapsis passages of the disrupted cloud. The inner and outer rings are misaligned by ~24 degrees, because they form from different gas streamers, which are affected by the SMBH gravitational focusing in different ways. The outer ring matches several properties (mass, rotation velocity, temperature, clumpiness) of the CNR in our Galactic centre. We speculate that the inner ring might account for the neutral gas observed in the central cavity.

We have run direct N-body simulations to investigate the impact of stellar evolution and dynamics on the structural properties of young massive (3x10^4 solar masses) star clusters (SCs) with different metallicities (Z=1, 0.1, 0.01 solar metallicity). Metallicity drives the mass loss by stellar winds and supernovae (SNe), with SCs losing more mass at high metallicity. We have simulated three sets of initial conditions, with different initial relaxation timescale. We find that the evolution of the half-mass radius of SCs depends on how fast two-body relaxation is with respect to the lifetime of massive stars. If core collapse is slow in comparison with stellar evolution, then mass loss by stellar winds and SNe is the dominant mechanism driving SC evolution, and metal-rich SCs expand more than metal-poor ones. In contrast, if core collapse occurs on a comparable timescale with respect to the lifetime of massive stars, then SC evolution depends on the interplay between mass loss and three-body encounters: dynamical heating by three-body encounters (mass loss by stellar winds and SNe) is the dominant process driving the expansion of the core in metal-poor (metal-rich) SCs. As a consequence, the half-mass radius of metal-poor SCs expands more than that of metal-rich ones. We also find core radius oscillations, which grow in number and amplitude as metallicity decreases.

The early evolution of a dense young star cluster (YSC) depends on the intricate connection between stellar evolution and dynamical processes. Thus, N-body simulations of YSCs must account for both aspects. We discuss N-body simulations of YSCs with three different metallicities (Z=0.01, 0.1 and 1 Zsun), including metallicity-dependent stellar evolution recipes and metallicity-dependent prescriptions for stellar winds and remnant formation. We show that mass-loss by stellar winds influences the reversal of core collapse. In particular, the post-collapse expansion of the core is faster in metal-rich YSCs than in metal-poor YSCs, because the former lose more mass (through stellar winds) than the latter. As a consequence, the half-mass radius expands more in metal-poor YSCs. We also discuss how these findings depend on the total mass and on the virial radius of the YSC. These results give us a clue to understand the early evolution of YSCs with different metallicity.