Search

We analyse how drag forces modify the orbits of objects moving through extended gaseous distributions. We consider how hydrodynamic (surface area) drag forces and dynamical friction (gravitational) drag forces drive the evolution of orbital eccentricity. While hydrodynamic drag forces cause eccentric orbits to become more circular, dynamical friction drag can cause orbits to become more eccentric. We develop a semi-analytic model that accurately predicts these changes by comparing the total work and torque applied to the orbit at periapse and apoapse. We use a toy model of a radial power-law density profile, $\rho \propto r^{-\gamma}$, to determine that there is a critical $\gamma = 3$ power index which separates the eccentricity evolution in dynamical friction: orbits become more eccentric for $\gamma < 3$ and circularize for $\gamma > 3$. We apply these findings to the infall of a Jupiter-like planet into the envelope of its host star. The hydrostatic envelopes of stars are defined by steep density gradients near the limb and shallower gradients in the interior. Under the influence of gaseous dynamical friction, an infalling object's orbit will first decrease in eccentricity, then increase. The critical separation that delineates these regimes is predicted by the local density slope and is linearly dependent on polytropic index. More broadly, our findings indicate that binary systems may routinely emerge from common envelope phases with non-zero eccentricities that were excited by the dynamical friction forces that drove their orbital tightening., Comment: 8 pages, 8 figures, Accepted for publication in MNRAS

Vector resonant relaxation (VRR) is known to be the fastest gravitational process that shapes the geometry of stellar orbits in nuclear star clusters. This leads to the realignment of the orbital planes on the corresponding VRR time scale $t_{\rm VRR}$ of a few million years, while the eccentricity $e$ and semimajor axis $a$ of the individual orbits are approximately conserved. The distribution of orbital inclinations reaches an internal equilibrium characterised by two conserved quantities, the total potential energy among stellar orbits, $E_{\rm tot}$, and the total angular momentum, $L_{\rm tot}$. On timescales longer than $t_\mathrm{VRR}$, the eccentricities and semimajor axes change slowly and the distribution of orbital inclinations are expected to evolve through a series of VRR equilibria. Using a Monte Carlo Markov Chain method, we determine the equilibrium distribution of orbital inclinations in the microcanonical ensemble with fixed $E_{\rm tot}$ and $L_{\rm tot}$ for isolated nuclear star clusters with a power-law distribution of $a$, $e$, and $m$, where $m$ is the stellar mass. We explore the possible equilibria for $9$ representative $E_{\rm tot}$--$L_{\rm tot}$ pairs that cover the possible parameter space. For all cases, the equilibria show anisotropic mass segregation where the distribution of more massive objects is more flattened than that for lighter objects. Given that stellar black holes are more massive than the average main sequence stars, these findings suggest that black holes reside in disc-like structures within nuclear star clusters for a wide range of initial conditions., Comment: Version accepted for publication in MNRAS

We investigate the dynamical evolution of an intermediate-mass black hole (IMBH) in a nuclear star cluster hosting a supermassive black hole (SMBH) and both a spherical and a flattened disk-like distribution of stellar-mass objects. We use a direct N-body (phiGPU) and an orbit-averaged (N-ring) numerical integrator to simulate the orbital evolution of stars and the IMBH. We find that the IMBH's orbit gradually aligns with the stellar disk if their mutual initial inclination is less than 90 degree. If it is larger than 90 degree, i.e. counterrotating, the IMBH does not align. Initially, the rate of orbital reorientation increases linearly with the ratio of the mass of the IMBH over the SMBH mass and it is orders of magnitude faster than ordinary (i.e. Chandrasekhar) dynamical friction, particularly for high SMBH masses. The semimajor axes of the IMBH and the stars are approximately conserved. This suggests that the alignment is predominantly driven by orbit-averaged gravitational torques of the stars, a process which may be called resonant dynamical friction. The stellar disk is warped by the IMBH, and ultimately increases its thickness. This process may offer a test for the viability of IMBH candidates in the Galactic Center. Resonant dynamical friction is not limited to IMBHs; any object much more massive than disk particles may ultimately align with the disk. This may have implications for the formation and evolution of black hole disks in dense stellar systems and gravitational wave source populations for LIGO, VIRGO, KAGRA, and LISA., Comment: 19 pages, 14 figures, 1 video, accepted by ApJ

We investigate the internal dynamics of anisotropic, rotating globular clusters with a multimass stellar population by performing new direct N-body simulations. In addition to the well-known radial mass segregation effect, where heavy stars and stellar remnants sink toward the center of the cluster, we find a mass segregation in the distribution of orbital inclinations as well. This newly discovered anisotropic mass segregation leads to the formation of a disk-like structure of massive objects near the equatorial plane of a rotating cluster. This result has important implications on the expected spatial distribution of black holes in globular clusters., Comment: 9 pages, 7 figures, 1 table

Gravitational torques among objects orbiting a supermassive black hole drive the rapid reorientation of orbital planes in nuclear star clusters (NSCs), a process known as vector resonant relaxation. In this Letter, we determine the statistical equilibrium of systems with a distribution of masses, semimajor axes, and eccentricities. We average the interaction over the apsidal precession time and construct a Monte Carlo Markov chain method to sample the microcanonical ensemble of the NSC. We examine the case of NSCs formed by 16 episodes of star formation or globular cluster infall. We find that the massive stars and stellar mass black holes form a warped disk, while low mass stars resemble a spherical distribution with a possible net rotation. This explains the origin of the clockwise disk in the Galactic center and predicts a population of black holes (BHs) embedded within this structure. The rate of mergers among massive stars, tidal disruption events of massive stars by BHs, and BH-BH mergers are highly increased in such disks. The first two may explain the origin of the observed G1 and G2 clouds, the latter may be important for gravitational wave detections with LIGO and VIRGO. More generally, black holes are expected to settle in disks in all dense spherical stellar systems assembled by mergers of smaller systems including globular clusters., Comment: 5 pages, 3 figures version published in Physical Review Letters

We introduce two novel time-dependent figures of merit for both online and offline optimizations of advanced gravitational-wave (GW) detector network operations with respect to (i) detecting continuous signals from known source locations and (ii) detecting GWs of neutron star binary coalescences from known local galaxies, which thereby have the highest potential for electromagnetic counterpart detection. For each of these scientific goals, we characterize an $N$-detector network, and all its $(N-1)$-detector subnetworks, to identify subnetworks and individual detectors (key contributors) that contribute the most to achieving the scientific goal. Our results show that aLIGO-Hanford is expected to be the key contributor in 2017 to the goal of detecting GWs from the Crab pulsar within the network of LIGO and Virgo detectors. For the same time period and for the same network, both LIGO detectors are key contributors to the goal of detecting GWs from the Vela pulsar, as well as to detecting signals from 10 high interest pulsars. Key contributors to detecting continuous GWs from the Galactic Center can only be identified for finite time intervals within each sidereal day with either the 3-detector network of the LIGO and Virgo detectors in 2017, or the 4-detector network of the LIGO, Virgo, and KAGRA detectors in 2019-2020. Characterization of the LIGO-Virgo detectors with respect to goal (ii) identified the two LIGO detectors as key contributors. Additionally, for all analyses, we identify time periods within a day when lock losses or scheduled service operations could result with the least amount of signal-to-noise or transient detection probability loss for a detector network., Comment: 16 pages, 5 figures