Dr Tjarda Boekholt

Abstract:

<jats:title>ABSTRACT</jats:title> <jats:p>Motivated by the increasing number of detections of merging black holes by LIGO-VIRGO-KAGRA, black hole (BH) binary mergers in the discs of active galactic nuclei (AGNs) is investigated as a possible merger channel. In this pathway, BH encounters in the gas disc form mutually bound BH binary systems through interaction with the gas in the disc and subsequently inspiral through gravitational torques induced by the local gas. To determine the feasibility of this merger pathway, we present the first three-dimensional global hydrodynamic simulations of the formation and evolution of a stellar-mass BH binaries AGN discs with three different AGN disc masses and five different initial radial separations. These 15 simulations show binary capture of prograde and retrograde binaries can be successful in a range of disc densities including cases well below that of a standard radiatively efficient alpha disc, identifying that the majority of these captured binaries are then subsequently hardened by the surrounding gas. The eccentricity evolution depends strongly on the orbital rotation where prograde binaries are governed by gravitational torques form their circumbinary mini disc, with eccentricities being damped, while for retrograde binaries the eccentricities are excited to &amp;gt;∼ 0.9 by accretion torques. In two cases, retrograde binaries ultimately undergo a close periapsis passage which results in a merger via gravitational waves after only a few thousand binary orbits. Thus, the merger time-scale can be far shorter than the AGN disc lifetime. These simulations support an efficient AGN disc merger pathway for BHs.</jats:p>

Abstract:

We investigate the effects of an oscillating gas filament on the dynamics of its embedded stellar clusters. Motivated by recent observational constraints, we model the host gas filament as a cylindrically symmetrical potential, and the star cluster as a Plummer sphere. In the model, the motion of the filament will produce star ejections from the cluster, leaving star cluster remnants that can be classified into four categories: (a) filament-associated clusters, which retain most of their particles (stars) inside the cluster and inside the filament; (b) destroyed clusters, where almost no stars are left inside the filament, and there is no surviving bound cluster; (c) ejected clusters, that leave almost no particles in the filament, since the cluster leaves the gas filament; and (d) transition clusters, corresponding to those clusters that remain in the filament, but that lose a significant fraction of particles due to ejections induced by filament oscillation. Our numerical investigation predicts that the Orion Nebula cluster is in the process of being ejected, after which it will most likely disperse into the field. This scenario is consistent with observations which indicate that the Orion Nebula cluster is expanding, and somewhat displaced from the integral-shaped filament ridgeline.

Abstract:

Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars, and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent and upcoming observing missions of exoplanets and stars in the galaxy help to provide constraints on the physics of tidal dissipation. It is timely to develop new N-body codes, which allow for experimentation with various tidal models and numerical implementations. We present the open-source N-body code TIDYMESS, which stands for ‘TIdal DYnamics of Multibody ExtraSolar Systems’. This code implements a Creep deformation law for the bodies, parametrized by their fluid Love numbers and fluid relaxation times. Due to tidal and centrifugal deformations, we approximate the general shape of a body to be an ellipsoid. We calculate the associated gravitational field to quadruple order, from which we derive the gravitational accelerations and torques. The equations of motion for the orbits, spins and deformations are integrated directly using a fourth-order integration method based on a symplectic composition. We implement a novel integration method for the deformations, which allows for a time-step solely dependent on the orbits, and not on the spin periods or fluid relaxation times. This feature greatly speeds up the calculations, while also improving the consistency when comparing different tidal regimes. We demonstrate the capabilities and performance of TIDYMESS, particularly in the niche regime of parameter space where orbits are chaotic and tides become non-linear.