Our Galactic center contains young stars, including the few million year old clockwise disk between 0.05 and 0.5 pc from the Galactic center, and the S-star cluster of B-type stars at a galactocentric distance of ~0.01 pc. Recent observations suggest the S-stars are remnants of tidally disrupted binaries from the clockwise disk. In particular, Koposov et al. 2020 discovered a hypervelocity star that was ejected from the Galactic center 5 Myr ago, with a velocity vector consistent with the disk. We perform a detailed study of this binary disruption scenario. First, we quantify the plausible range of binary semimajor axes in the disk. Dynamical evaporation of such binaries is dominated by other disk stars rather than by the isotropic stellar population. For the expected range of semimajor axes in the disk, binary tidal disruptions can reproduce the observed S-star semimajor axis distribution. Reproducing the observed thermal eccentricity distribution of the S-stars requires an additional relaxation process. The flight time of the Koposov star suggests that this process must be effective within 10 Myr. We consider three possibilities: (i) scalar resonant relaxation from the observed isotropic star cluster, (ii) torques from the clockwise disk, and (iii) an intermediate-mass black hole. We conclude that the first and third mechanisms are fast enough to reproduce the observed S-star eccentricity distribution. Finally, we show that the primary star from an unequal-mass binary would be deposited at larger semimajor axes than the secondary, possibly explaining the dearth of O stars among the S-stars.

Axisymmetric disks of high eccentricity, low mass bodies on near-Keplerian orbits are unstable to an out-of-plane buckling. This "inclination instability" exponentially grows the orbital inclinations, raises perihelia distances and clusters in argument of perihelion. Here we examine the instability in a massive primordial scattered disk including the orbit-averaged gravitational influence of the giant planets. We show that differential apsidal precession induced by the giant planets will suppress the inclination instability unless the primordial mass is $\gtrsim 20$ Earth masses. We also show that the instability should produce a "perihelion gap" at semi-major axes of hundreds of AU, as the orbits of the remnant population are more likely to have extremely large perihelion distances ($\mathcal{O}(100~\rm{AU})$) than intermediate values.

Disks of low-mass bodies on high-eccentricity orbits in near-Keplerian potentials can be dynamically unstable to buckling out of the plane. In this letter, we present $N$-body simulations of the long-term behavior of such a system, finding apsidal clustering of the orbits in the disk plane. The timescale over which the clustering is maintained increases with number of particles, suggesting that lopsided configurations are stable at large $N$. This discovery may explain the observed apsidal ($\varpi$) clustering of extreme trans-Neptunian Objects in the outer solar system.

Eccentric nuclear disks (ENDs) are a type of star cluster in which the stars lie on eccentric, apsidally-aligned orbits in a disk around a central supermassive black hole (SMBH). These disks can produce a high rate of tidal disruption events (TDEs) via secular gravitational torques. Previous studies of ENDs have included stars with only one mass. Here, we present the first study of an eccentric nuclear disk with two stellar species. We show that ENDs show radial mass segregation consistent with previous results from other cluster types. Additionally, ENDs show vertical mass segregation by which the heavy stars sink to lower inclinations than light stars. These two effects cause heavy stars to be more susceptible to tidal disruption, which can be seen in the higher fraction of heavy stars that are disrupted compared to light stars.

Massive bodies undergo orbital eccentricity oscillations when embedded in an axisymmetric disk of smaller mass orbits. These eccentricity oscillations are driven by secular torques that seek to equalize the apsidal precession rates of all the orbits in the disk. We investigate this mechanism within the context of detached objects in the outer Solar System, but we find that the oscillation timescale is too long for it to be dynamically important. It could however be interesting for phenomenon a bit farther from home; namely, feeding supermassive blackholes and polluting the surfaces of white dwarf stars.

An eccentric nuclear disk consists of stars moving on apsidally-aligned orbits around a central black hole. The secular gravitational torques that dynamically stabilize these disks can also produce tidal disruption events (TDEs) at very high rates in Newtonian gravity. General relativity, however, is known to quench secular torques via rapid apsidal precession. Here we show that for a disk to black hole mass ratio $\gtrsim 10^{-3}$, the system is in the full loss cone regime. The magnitude of the torque per orbital period acting on a stellar orbit means that general relativistic precession does not have a major effect on the dynamics. Thus we find that TDE rates from eccentric nuclear disks are not affected by general relativistic precession. Furthermore, we show that orbital elements between successive TDEs from eccentric nuclear disks are correlated, potentially resulting in unique observational signatures.

Axisymmetric disks of eccentric Kepler orbits are vulnerable to an instability which causes orbits to exponentially grow in inclination, decrease in eccentricity, and cluster in their angle of pericenter. Geometrically, the disk expands to a cone shape which is asymmetric about the mid-plane. In this paper, we describe how secular gravitational torques between individual orbits drive this "inclination instability". We derive growth timescales for a simple two-orbit model using a Gauss $N$-ring code, and generalize our result to larger $N$ systems with $N$-body simulations. We find that two-body relaxation slows the growth of the instability at low $N$ and that angular phase coverage of orbits in the disk is important at higher $N$. As $N \to \infty$, the e-folding timescale converges to that expected from secular theory.

In some galaxies, the stars orbiting the supermassive black hole take the form of an eccentric nuclear disk, in which every star is on a coherent, apsidally-aligned orbit. The most famous example of an eccentric nuclear disk is the double nucleus of Andromeda, and there is strong evidence for many more in the local universe. Despite their apparent ubiquity however, a dynamical explanation for their longevity has remained a mystery: differential precession should wipe out large-scale apsidal-alignment on a short timescale. Here we identify a new dynamical mechanism which stabilizes eccentric nuclear disks, and explain for first time the negative eccentricity gradient seen in the Andromeda nucleus. The stabilizing mechanism drives oscillations of the eccentricity vectors of individual orbits, both in direction (about the mean body of the disk) and in magnitude. Combined with the negative eccentricity gradient, the eccentricity oscillations push some stars near the inner edge of the disk extremely close to the black hole, potentially leading to tidal disruption events. Order of magnitude calculations predict extremely high rates in recently-formed eccentric nuclear disks ($\sim0.1 - 1$ ${\rm yr}^{-1} {\rm gal}^{-1}$). Unless the stellar disks are replenished, these rates should decrease with time as the disk depletes in mass. If eccentric nuclear disks form during gas-rich major mergers, this may explain the preferential occurrence of tidal disruption events in recently-merged and post-merger (E+A/K+A) galaxies.

We find that clouds of optically-thin, pressure-confined gas are prone to fragmentation as they cool below $\sim10^6$ K. This fragmentation follows the lengthscale $\sim{c}_{\text{s}}\,t_{\text{cool}}$, ultimately reaching very small scales ($\sim{0.1} \text{pc}/n$) as they reach the temperature $\sim10^4$ K at which hydrogen recombines. While this lengthscale depends on the ambient pressure confining the clouds, we find that the column density through an individual fragment $N_{\text{cloudlet}}\sim10^{17} \text{cm}^{-3}$ is essentially independent of environment; this column density represents a characteristic scale for atomic gas at $10^4$ K. We therefore suggest that "clouds" of cold, atomic gas may in fact have the structure of a mist or a fog, composed of tiny fragments dispersed throughout the ambient medium. We show that this scale emerges in hydrodynamic simulations, and that the corresponding increase in the surface area may imply rapid entrainment of cold gas. We also apply it to a number of observational puzzles, including the large covering fraction of diffuse gas in galaxy halos, the broad line widths seen in quasar and AGN spectra, and the entrainment of cold gas in galactic winds. While our simulations make a number of assumptions and thus have associated uncertainties, we show that this characteristic scale is consistent with a number of observations, across a wide range of astrophysical environments. We discuss future steps for testing, improving, and extending our model.

We study the dynamical evolution of the putative gas clouds G1 and G2 recently discovered in the Galactic center. Following earlier studies suggesting that these two clouds are part of a larger gas streamer, we combine their orbits into a single trajectory. Since the gas clouds experience a drag force from background gas, this trajectory is not exactly Keplerian. By assuming the G1 and G2 clouds trace this trajectory, we fit for the drag force they experience and thus extract information about the accretion flow at a distance of thousands of Schwarzschild radii from the black hole. This range of radii is important for theories of black hole accretion, but is currently unconstrained by observations. In this paper we extend our previous work by accounting for radial forces due to possible inflow or outflow of the background gas. Such radial forces drive precession in the orbital plane, allowing a slightly better fit to the G1 and G2 data. This precession delays the pericenter passage of G2 by 4-5 months relative to estimates derived from a Keplerian orbital fit; if it proves possible to identify the pericenter time observationally, this enables an immediate test of whether G1 and G2 are gas clouds part of a larger gas streamer. If G2 is indeed a gas cloud, its closest approach likely occurred in late summer 2014, after many of the observing campaigns monitoring G2's anticipated pericenter passage ended. We discuss how this affects interpretation of the G2 observations.

Disks of bodies orbiting a much more massive central object are extremely common in astrophysics. When the orbits comprising such disks are eccentric, we show they are susceptible to a new dynamical instability. Gravitational forces between bodies in the disk drive exponential growth of their orbital inclinations and clustering in their angles of pericenter, expanding an initially thin disk into a conical shape by giving each orbit an identical 'tilt' with respect to the disk plane. This new instability dynamically produces the unusual distribution of orbits observed for minor planets beyond Neptune, suggesting that the instability has shaped the outer Solar System. It also implies a large initial disk mass (1-10 Earth masses) of scattered bodies at hundreds of AU; we predict increasing numbers of detections of minor planets clustered in their angles of pericenter with high inclinations.

The massive black hole in our galactic center, Sgr A*, accretes only a small fraction of the gas available at its Bondi radius. The physical processes determining this accretion rate remain unknown, partly due to a lack of observational constraints on the gas at distances between ~10 and ~10$^5$ Schwarzschild radii (Rs) from the black hole. Recent infrared observations identify low-mass gas clouds, G1 and G2, moving on highly eccentric, nearly co-planar orbits through the accretion flow around Sgr A*. Although it is not yet clear whether these objects contain embedded stars, their extended gaseous envelopes evolve independently as gas clouds. In this paper we attempt to use these gas clouds to constrain the properties of the accretion flow at ~10$^3$ Rs. Assuming that G1 and G2 follow the same trajectory, we model the small differences in their orbital parameters as evolution resulting from interaction with the background flow. We find evolution consistent with the G-clouds originating in the clockwise disk. Our analysis enables the first unique determination of the rotation axis of the accretion flow: we localize the rotation axis to within 20 degrees, finding an orientation consistent with the parsec-scale jet identified in x-ray observations and with the circumnuclear disk, a massive torus of molecular gas ~1.5 pc from Sgr A*. This suggests that the gas in the accretion flow comes predominantly from the circumnuclear disk, rather than the winds of stars in the young clockwise disk. This result will be tested by the Event Horizon Telescope within the next year. Our model also makes testable predictions for the orbital evolution of G1 and G2, falsifiable on a 5-10 year timescale.

We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to the one of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (i) G1 was originally on the G2 orbit preceding G2 by 13 years and (ii) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.

We present three-dimensional magnetohydrodynamic simulations of magnetized gas clouds accelerated by hot winds. We initialize gas clouds with tangled internal magnetic fields and show that this field suppresses the disruption of the cloud: rather than mixing into the hot wind as found in hydrodynamic simulations, cloud fragments end up co-moving and in pressure equilibrium with their surroundings. We also show that a magnetic field in the hot wind enhances the drag force on the cloud by a factor ~(1+v_A^2/v_wind^2)$, where v_A is the Alfven speed in the wind and v_wind measures the relative speed between the cloud and the wind. We apply this result to gas clouds in several astrophysical contexts, including galaxy clusters, galactic winds, the Galactic center, and the outskirts of the Galactic halo. Our results can explain the prevalence of cool gas in galactic winds and galactic halos and how such cool gas survives in spite of its interaction with hot wind/halo gas. We also predict that drag forces can lead to a deviation from Keplerian orbits for the G2 cloud in the galactic center.

The center of our galaxy is home to a massive black hole, SgrA*, and a nuclear star cluster containing stellar populations of various ages. While the late type stars may be too old to have retained memory of their initial orbital configuration, and hence formation mechanism, the kinematics of the early type stars should reflect their original distribution. In this contribution we present a new statistic which uses directly-observable kinematical stellar data to infer orbital parameters for stellar populations, and is capable of distinguishing between different origin scenarios. We use it on a population of B-stars in the Galactic center that extends out to large radii (0.5 pc) from the massive black hole. We find that the high K-magnitude population form an eccentric distribution, suggestive of a Hills binary-disruption origin.

We present a new directly-observable statistic which uses sky position and proper motion of stars near the Galactic center massive black hole to identify populations with high orbital eccentricities. It is most useful for stars with large orbital periods for which dynamical accelerations are difficult to determine. We apply this statistic to a data set of B-stars with projected radii 0."1 < p < 25" (~0.004 - 1 pc) from the massive black hole in the Galactic center. We compare the results with those from N-body simulations to distinguish between scenarios for their formation. We find that the scenarios favored by the data correlate strongly with particular K-magnitude intervals, corresponding to different zero-age main-sequence (MS) masses and lifetimes. Stars with 14 < mK < 15 (15 - 20 solar masses, t_{MS} = 8-13 Myr) match well to a disk formation origin, while those with mK > 15 (<15 solar masses, t_{MS} >13 Myr), if isotropically distributed, form a population that is more eccentric than thermal, which suggests a Hills binary-disruption origin.

We identify a gravitational-dynamical process in near-Keplerian potentials of galactic nuclei that occurs when an intermediate-mass black hole (IMBH) is migrating on an eccentric orbit through the stellar cluster towards the central supermassive black hole (SMBH). We find that, apart from conventional dynamical friction, the IMBH experiences an often much stronger systematic torque due to the secular (i.e., orbit-averaged) interactions with the cluster's stars. The force which results in this torque is applied, counterintuitively, in the same direction as the IMBH's precession and we refer to its action as "secular-dynamical anti-friction" (SDAF). We argue that SDAF, and not the gravitational ejection of stars, is responsible for the IMBH's eccentricity increase seen in the initial stages of previous N-body simulations. Our numerical experiments, supported by qualitative arguments, demonstrate that (1) when the IMBH's precession direction is artificially reversed, the torque changes sign as well, which decreases the orbital eccentricity, (2) the rate of eccentricity growth is sensitive to the IMBH migration rate, with zero systematic eccentricity growth for an IMBH whose orbit is artificially prevented from inward migration, and (3) SDAF is the strongest when the central star cluster is rapidly rotating. This leads to eccentricity growth/decrease for the clusters rotating in the opposite/same direction relative to the IMBH's orbital motion.

The angular momentum evolution of stars close to massive black holes (MBHs) is driven by secular torques. In contrast to two-body relaxation, where interactions between stars are incoherent, the resulting resonant relaxation (RR) process is characterized by coherence times of hundreds of orbital periods. In this paper, we show that all the statistical properties of RR can be reproduced in an autoregressive moving average (ARMA) model. We use the ARMA model, calibrated with extensive N-body simulations, to analyze the long-term evolution of stellar systems around MBHs with Monte Carlo simulations. We show that for a single-mass system in steady-state, a depression is carved out near an MBH as a result of tidal disruptions. Using Galactic center parameters, the extent of the depression is about 0.1 pc, of similar order to but less than the size of the observed "hole" in the distribution of bright late-type stars. We also find that the velocity vectors of stars around an MBH are locally not isotropic. In a second application, we evolve the highly eccentric orbits that result from the tidal disruption of binary stars, which are considered to be plausible precursors of the "S-stars" in the Galactic center. We find that RR predicts more highly eccentric (e > 0.9) S-star orbits than have been observed to date.

In this paper we revisit the "eccentric disc instability", an instability which occurs in coherently eccentric discs of stars orbiting massive black holes (MBHs) embedded in stellar clusters, which results in stars achieving either very high or low eccentricities. The preference for stars to attain higher or lower eccentricities depends significantly on the density distribution of the surrounding stellar cluster. Here we discuss its mechanism and the implications for the Galactic Centre, home to at least one circum-MBH stellar disc.

Two-body energy exchange between stars orbiting massive black holes (MBHs) leads to the formation of a power-law density distribution n(r)~r^(-a) that diverges towards the MBH. For a single mass population, a=7/4 and the flow of stars is much less than N(<r)/t_r (enclosed number of stars per relaxation time). This "zero-flow" solution is maintained for a multi-mass system for moderate mass ratios or systems where there are many heavy stars, and slopes of 3/2<a<2 are reached, with steeper slopes for the more massive stars. If the heavy stars are rare and massive however, the zero-flow limit breaks down and much steeper distributions are obtained. We discuss the physics driving mass-segregation with the use of Fokker-Planck calculations, and show that steady state is reached in 0.2-0.3 t_r. Since the relaxation time in the Galactic centre (GC) is t_r ~2-3 * 10^(10) yr, a cusp should form in less than a Hubble time. The absence of a visible cusp of old stars in the GC poses a challenge to these models, suggesting that processes other than two-body relaxation have played a role. We discuss astrophysical processes within the GC that depend crucially on the details of the stellar cusp.

We identify a new secular instability of eccentric stellar disks around supermassive black holes. We show that retrograde precession of the stellar orbits, due to the presence of a stellar cusp, induces coherent torques that amplify deviations of individual orbital eccentricities from the average, and thus drive all eccentricities away from their initial value. We investigate the instability using N-body simulations, and show that it can drive individual orbital eccentricities to significantly higher or lower values on the order of a precession time-scale. This physics is relevant for the Galactic center, where massive stars are likely to form in eccentric disks around the SgrA* black hole. We show that the observed bimodal eccentricity distribution of disk stars in the Galactic center is in good agreement with the distribution resulting from the eccentricity instability and demonstrate how the dynamical evolution of such a disk results in several of its stars acquiring high (1-e < 0.1) orbital eccentricity. Binary stars on such highly eccentric orbits would get tidally disrupted by the SgrA* black hole, possibly producing both S-stars near the black hole and high-velocity stars in the Galactic halo.