Source author record

Jihad R. Touma

Jihad R. Touma appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
3topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2016arXiv

Stellar Dynamics around a Massive Black Hole I: Secular Collisionless Theory

We present a theory in three parts, of the secular dynamics of a (Keplerian) stellar system of mass $M$ orbiting a black hole of mass $M_\bullet \gg M$. Here we describe the collisionless dynamics; Papers II and III are on the (collisional) theory of Resonant Relaxation. The mass ratio, $\varepsilon = M/M_\bullet \ll 1$, is a natural small parameter implying a separation of time scales between the short Kepler orbital periods and the longer orbital precessional periods. The collisionless Boltzmann equation (CBE) for the stellar distribution function (DF) is averaged over the fast Kepler orbital phase using the method of multiple scales. The orbit-averaged system is described by a secular DF, $F$, in a reduced phase space. $F$ obeys a secular CBE that includes stellar self-gravity, general relativistic corrections up to 1.5 post-Newtonian order, and external sources varying over secular times. Secular dynamics, even with general time dependence, conserves the semi-major axis of every star. This additional integral of motion promotes extra regularity of the stellar orbits, and enables the construction of equilibria, $F_0$, through a secular Jeans theorem. A linearized secular CBE determines the response and stability of $F_0$. Spherical, non-rotating equilibria may support long-lived, warp-like distortions. We also prove that an axisymmetric, zero-thickness, flat disc is secularly stable to all in-plane perturbations, when its DF, $F_0$, is a monotonic function of the angular momentum at fixed energy.

preprint2016arXiv

Stellar Dynamics around a Massive Black Hole II: Resonant Relaxation

We present a first-principles theory of Resonant Relaxation (RR) of a low mass stellar system orbiting a more massive black hole (MBH). We first extend the kinetic theory of Gilbert (1968) to include the Keplerian field of a black hole of mass $M_\bullet$. Specializing to a Keplerian stellar system of mass $M \ll M_\bullet$, we use the orbit-averaging method of Sridhar & Touma (2015; Paper I) to derive a kinetic equation for RR. This describes the collisional evolution of a system of $N \gg 1$ Gaussian Rings in a reduced 5-dim space, under the combined actions of self-gravity, 1 PN and 1.5 PN relativistic effects of the MBH and an arbitrary external potential. In general geometries RR is driven by both apsidal and nodal resonances, so the distinction between scalar-RR and vector-RR disappears. The system passes through a sequence of quasi-steady secular collisionless equilibria, driven by irreversible 2-Ring correlations that accrue through gravitational interactions, both direct and collective. This correlation function is related to a `wake function', which is the linear response of the system to the perturbation of a chosen Ring. The wake function is easier to appreciate, and satisfies a simpler equation, than the correlation function. We discuss general implications for the interplay of secular dynamics and non-equilibrium statistical mechanics in the evolution of Keplerian stellar systems toward secular thermodynamic equilibria, and set the stage for applications to the RR of axisymmetric discs in Paper III.

preprint2016arXiv

Stellar Dynamics around a Massive Black Hole III: Resonant Relaxation of Axisymmetric Discs

We study the Resonant Relaxation (RR) of an axisymmetric low mass (or Keplerian) stellar disc orbiting a more massive black hole (MBH). Our recent work on the general kinetic theory of RR is simplified in the standard manner by ignoring the effects of `gravitational polarization', and applied to a zero-thickness, flat, axisymmetric disc. The wake of a stellar orbit is expressed in terms of the angular momenta exchanged with other orbits, and used to derive a kinetic equation for RR under the combined actions of self-gravity, 1 PN and 1.5 PN relativistic effects of the MBH and an arbitrary external axisymmetric potential. This is a Fokker-Planck equation for the stellar distribution function (DF), wherein the diffusion coefficients are given self-consistently in terms of contributions from apsidal resonances between pairs of stellar orbits. The physical kinetics is studied for the two main cases of interest. (1) `Lossless' discs in which the MBH is not a sink of stars, and disc mass, angular momentum and energy are conserved: we prove that general H-functions can increase or decrease during RR, but the Boltzmann entropy is (essentially) unique in being a non-decreasing function of time. Therefore secular thermal equilibria are maximum entropy states, with DFs of the Boltzmann form; the two-Ring correlation function at equilibrium is computed. (2) Discs that lose stars to the MBH through an `empty loss-cone': we derive expressions for the MBH feeding rates of mass, angular momentum and energy in terms of the diffusive flux at the loss-cone boundary.

preprint2015arXiv

Migration into a Companion's Trap: Disruption of Multiplanet Systems in Binaries

Most exoplanetary systems in binary stars are of S--type, and consist of one or more planets orbiting a primary star with a wide binary stellar companion. Gravitational forcing of a single planet by a sufficiently inclined binary orbit can induce large amplitude oscillations of the planet's eccentricity and inclination through the Kozai-Lidov (KL) instability. KL cycling was invoked to explain: the large eccentricities of planetary orbits; the family of close--in hot Jupiters; and the retrograde planetary orbits in eccentric binary systems. However, several kinds of perturbations can quench the KL instability, by inducing fast periapse precessions which stabilize circular orbits of all inclinations: these could be a Jupiter--mass planet, a massive remnant disc or general relativistic precession. Indeed, mutual gravitational perturbations in multiplanet S--type systems can be strong enough to lend a certain dynamical rigidity to their orbital planes. Here we present a new and faster process that is driven by this very agent inhibiting KL cycling. Planetary perturbations enable secular oscillations of planetary eccentricities and inclinations, also called Laplace--Lagrange (LL) eigenmodes. Interactions with a remnant disc of planetesimals can make planets migrate, causing a drift of LL mode periods which can bring one or more LL modes into resonance with binary orbital motion. The results can be dramatic, ranging from excitation of large eccentricities and mutual inclinations to total disruption. Not requiring special physical or initial conditions, binary resonant driving is generic and could have profoundly altered the architecture of many S--type multiplanet systems. It can also weaken the multiplanet occurrence rate in wide binaries, and affect planet formation in close binaries.

preprint2012arXiv

Counter-rotating stellar discs around a massive black hole: self-consistent, time-dependent dynamics

We formulate the collisionless Boltzmann equation (CBE) for dense star clusters that lie within the radius of influence of a massive black hole in galactic nuclei. Our approach to these nearly Keplerian systems follows that of Sridhar and Touma (1999): Delaunay canonical variables are used to describe stellar orbits and we average over the fast Keplerian orbital phases. The stellar distribution function (DF) evolves on the longer time scale of precessional motions, whose dynamics is governed by a Hamiltonian, given by the orbit-averaged self-gravitational potential of the cluster. We specialize to razor-thin, planar discs and consider two counter-rotating ("$\pm$") populations of stars. To describe discs of small eccentricities, we expand the $\pm$ Hamiltonian to fourth order in the eccentricities, with coefficients that depend self-consistently on the $\pm$ DFs. We construct approximate $\pm$ dynamical invariants and use Jeans' theorem to construct time-dependent $\pm$ DFs, which are completely described by their centroid coordinates and shape matrices. When the centroid eccentricities are larger than the dispersion in eccentricities, the $\pm$ centroids obey a set of 4 autonomous ordinary differential equations. We show that these can be cast as a two-degree of freedom Hamiltonian system which is nonlinear, yet integrable. We study the linear instability of initially circular discs and derive a criterion for the counter-rotating instability. We then explore the rich nonlinear dynamics of counter-rotating discs, with focus on the variety of steadily precessing eccentric configurations that are allowed. The stability and properties of these configurations are studied as functions of parameters such as the disc mass ratios and angular momentum.

preprint2012arXiv

The Doubling of Stellar Black Hole Nuclei

It is strongly believed that Andromeda's double nucleus signals a disk of stars revolving around its central super-massive black hole on eccentric Keplerian orbits with nearly aligned apsides. A self-consistent stellar dynamical origin for such apparently long-lived alignment has so far been lacking, with indications that cluster self-gravity is capable of sustaining such lopsided configurations if and when stimulated by external perturbations. Here, we present results of N-body simulations which show unstable counter-rotating stellar clusters around super-massive black holes saturating into uniformly precessing lopsided nuclei. The double nucleus in our featured experiment decomposes naturally into a thick eccentric disk of apo-apse aligned stars which is embedded in a lighter triaxial cluster. The eccentric disk reproduces key features of Keplerian disk models of Andromeda's double nucleus; the triaxial cluster has a distinctive kinematic signature which is evident in HST observations of Andromeda's double nucleus, and has been difficult to reproduce with Keplerian disks alone. Our simulations demonstrate how the combination of eccentric disk and triaxial cluster arises naturally when a star cluster accreted over a pre-existing and counter-rotating disk of stars, drives disk and cluster into a mutually destabilizing dance. Such accretion events are inherent to standard galaxy formation scenarios. They are here shown to double stellar black hole nuclei as they feed them.