Research


You will find here a brief summary of the various projects I have been involved in.

Radial dependence of the WKB basis elements

Secular resonant dressed orbital diffusion - I. Method and WKB limit for tepid discs


We derived the secular diffusion equation of a self-gravitating collisionless system induced by external stochastic perturbations. In the case of a tepid galactic disc, relying on the WKB assumption that only tightly wound transient spirals are sustained by the disc, we obtained a simple quadrature for the diffusion coefficients, providing a straightforward understanding of the loci of maximal diffusion.

Prediction of the secular diffusion flux exhibiting a narrow ridge of resonant orbits

Secular resonant dressed orbital diffusion - II. Application to an isolated self-similar tepid galactic disc


We recovered the main orbital signatures of the secular evolution of an isolated self-gravitating stellar disc. The shot-noise-driven formation of narrow ridges of resonant orbits is recovered in the WKB limit of tightly wound transient spirals, in agreement with numerical simulations. This justifies the relevance of the dressed Fokker-Planck formalism in angle-action variables to describe the secular evolution of such systems.

Velocity flux in action space induced by Poisson shot noise

Self-gravity, Resonances, and Orbital Diffusion in Stellar Disks


Fluctuations in a stellar system's gravitational field cause the orbits of stars to evolve. The resulting evolution of the system can be computed with the orbit-averaged Fokker-Planck equation. We presented the formalism that enables one to compute the diffusion tensor from a given source of noise when the system's gravitational dynamical response to that noise is included. This formalism, which recovers the formations of narrow ridges of enhanced density in action space, appears as the ideal framework in which to study the long-term evolution of all kinds of stellar discs.

Illustration of a discrete resonant encounter between two stars, captured by the Balescu-Lenard equation. In the appropriate rotating frame, orbits are closed and resonate. In the long-run, this deforms orbits.

Secular diffusion in discrete self-gravitating tepid discs I : analytic solution in the tightly wound limit


We described the secular evolution of an infinitely thin isolated discrete self-gravitating stellar disc using the inhomogeneous Balescu-Lenard equation. Assuming that only tightly wound transient spirals are present in the disc, a WKB approximation provides a simple quadrature for the corresponding drift and diffusion coefficients. When applied to the secular evolution of an isolated stationary discrete self-gravitating Mestel disc, it predicts the formation of a ridge-like feature in action space, in agreement with simulations, but over-estimates the timescale involved in its appearance. Swing amplification is needed to resolve this discrepancy.

Comparison of the resonance ridges predicted via the Balescu-Lenard equation (top) and observed in the numerical simulations of Sellwood (2012) (bottom)

Secular diffusion in discrete self-gravitating tepid discs II: accounting for swing amplification via the matrix method


We investigated the secular evolution of an infinitely thin discrete self-gravitating stella disc using the inhomogeneous Balescu-Lenard equation in terms of angle-action variables. We implemented numerically the matrix method to capture the induced graviational polarisation. The position/shape of the induced resonant ridge are found to be in very good agreement with the numerical simulations, as well as the diffusion timescales. Quantitative comparisons with N-body simulations also yield consistent scalings with the number of particles.

Functional rewriting of the two first equations of the BBGKY hierarchy

Functional integral approach to the kinetic theory of inhomogeneous systems


We present a derivation of the kinetic equation describing the secular evolution of spatially inhomogeneous systems with long-range interactions, the so-called inhomogeneous Landau equation, by relying on a functional integral formalism.

Illustration of an axisymmetric system of 2D point vortices

Functional integral derivation of the kinetic equation of two-dimensional point vortices


We present a brief derivation of the kinetic equation describing the secular evolution of point vortices in two-dimensional hydrodynamics, by relying on a functional integral formalism.

Illustration of the typical dependence of the precession frequencies (induced by the cluster and relativistic corrections) as a function of the distance to the central BH

The secular evolution of discrete quasi-Keplerian systems. I. Kinetic theory of stellar clusters near black holes


We derive the kinetic equation that describes the secular evolution of a large set of particles orbiting a dominant massive object, such as stars bound to a supermassive black hole or a proto-planetary debris disc encircling a star. This degenerate Balescu-Lenard equation describes self-consistently the long-term evolution of the distribution of quasi-Keplerian orbits around the central object: it is the master equation that describes the secular effects of resonant relaxation.

Expression of the resonant dressed dynamical friction as given by the inhomogeneous Balescu-Lenard equation

Dressed diffusion and friction coefficients in inhomogeneous multi-component self-gravitating systems


We derive general self-consistent expressions for the coefficients of diffusion and dynamical friction in a stable, bound, multicomponent, self-gravitating, and inhomogeneous system.

Norm of the collisional diffusion flux given by the thickened WKB limit of the Balescu-Lenard equation. The presence of an enhanced diffusion flux in the inner region of the disc is compatible with the self-induced formation of vertical resonant ridges.

Resonant thickening of self-gravitating discs: imposed or self-induced orbital diffusion in the tightly wound limit


We investigate the secular thickening of a self-gravitating stellar galactic disc using the dressed collisionless Fokker-Planck equation and the inhomogeneous multi-component Balescu-Lenard equation, by relying on a generalised thickened WKB approximation. When applied to a tepid stable tapered disc perturbed by shot noise, these two frameworks predict the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but over-estimates the timescale involved in their appearance. Swing amplification is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs.

Response of the Galactic's DF to a bar perturbation, in the region of trapping.

Distribution functions for resonantly trapped orbits in the Galactic disc


We show how to compute the response of the DF of a galactic disc to a bar-like perturbation in the region of orbit-trapping, where the traditional perturbation theory fails.

Illustration in physical space of the mass segregation of a two-component quasi-Keplerian disc.

The secular evolution of discrete quasi-Keplerian systems. II. Application to a multi-mass axisymmetric disc around a supermassive black hole


Relying on Gauss' method, we compute the drift and diffusion coefficients characterising the properties of the resonant relaxation of a razor-thin quasi-Keplerian disc.

Phase-space structure of the diffusion coefficients of scalar resonant relaxation.

Scalar Resonant Relaxation of Stars around a Massive Black Hole


We describe from first principles the process of scalar resonant relaxation (relaxation of the norm of the angular momentum) in galactic nuclei.

Rate of change a globular cluster's DF (blue increasing, red decreasing) in action space as predicted by the inhomogeneous Balescu-Lenard equation.

Revisiting relaxation in globular clusters


We adapt the Balescu-Lenard equation to spherical systems such as globular clusters. We use it to estimate the flux through action space, and compare it with classical theory.

Diffusion coefficient in a self-interacting HMF bath, obtained through different methods.

Relaxation in self-gravitating systems


We show how the stochastic η-formalism allows for the recovery of the Balescu-Lenard and Landau equations. This approach provides a new view of the resonant diffusion processes associated with long-term orbital distortions.

Inspiral of a massive object in a FDM halo.

Relaxation in a Fuzzy Dark Matter Halo


By treating the FDM fluctuations as quasiparticles, we show how the relaxation in a FDM halo can be analysed quantitatively with the same tools used to analyse classical two-body relaxation in an N-body system.

Unstable growth rate in a Mestel disc for different softening kernels, as predicted by linear theory.

How gravitational softening affects galaxy stability. I. Linear mode analysis of disc galaxies


We use linear pertubation theory to investigate the effect of graviational softening on the spiral eigenmodes of razor-thin stellar discs.

Second-order diffusion coefficients predicted by the Balescu-Lenard equation and compared with direct simulations.

Secular dynamics of long-range interacting particles on a sphere in the axisymmetric limit


We investigate the secular dynamics of long-range interacting particles moving on a sphere, in the limit of an axisymmetric mean field potential.

Unstable spiral mode a grooved stellar disc, as predicted by linear theory.

Instabilities in disc galaxies: from noise to grooves to spirals


Using linear response theory, we show how grooves in action space (spontaneously generated via finite-N effects) can lead unstable modes in razor-thin stellar discs.

Illustration of an unstable mode of a rotating stellar cluster.

Mapping the stability of stellar rotating spheres via linear response theory


We investigate the linear stability of rotating stellar clusters with varying velocity anisotropies.

Illustration of the stochastic dynamics of stellar orientations in a galactic nuclei.

Vector Resonant Relaxation of Stars around a Massive Black Hole


We derive from first principles the statistical properties of random walks undergone by stellar orbital orientations in galactic nuclei.

Diffusion flux for various temperatures (Q), compared with kinetic theory.

Kinetic theory of 1D homogeneous long-range interacting systems sourced by 1/N^2 effects


We present a 1/N^2 kinetic theory of the HMF model in the limit where collective effects can be neglected.

Kinematic spiral signature within the Milky Way's disc.

The Strength of the Dynamical Spiral Perturbation in the Galactic Disk


We characterise the amplitude of the spiral perturbation in the Milky Way.

Illustration of the (slow) separation of two neighbour stars with similar initial orientations.

Young stellar cluster dilution near supermassive black holes: the impact of Vector Resonant Relaxation on neighbour separation


We develop the theory of "neighbor separation" for the vector resonant relaxation in galactic nuclei.

Diffusion flux in a 1D homogeneous model, numerical measurements vs. kinetic predictions.

Kinetic theory of one-dimensional homogeneous long-range interacting systems with an arbitrary potential of interaction


We develop the 1/N^2 kinetic theory of 1D homogeneous long-range interacting systems.

Illustration of the velocity in a fuzzy halo as predicted by kinetic theory and measured in numerical simulations.

Relaxation in a Fuzzy Dark Matter Halo. II. Self-consistent kinetic equations


We derive a kinetic equation to describe the self-consistent relaxation of fuzzy dark matter halos.

Small-scale Coulomb logarithm as recovered from the RR kinetic theory.

Resonant and non-resonant relaxation of globular clusters


We compute the diffusion rate in action space using both resonant and non-resonant kinetic theories. We compare with direct N-body measurements

RR and NR eccentricity diffusion coefficients in a galactic nucleus.

Mapping the Galactic centre's dark cluster via Resonant Relaxation


We use scalar resonant relaxation and the observed distribution of the S-stars' eccentricities to constrain the stellar content around SgrA*.

Diffusion coefficients in eccentricities as predicted by kinetic theory and measured in numerical simulations.

Secular Dynamics around a Supermassive Black Hole via Multipole Expansion


Within the secular approximation of smearing out stars along their underlying Keplerian orbits, we detail how a multipole expansion of the pairwise interaction between the stars yields an efficient numerical code to investigate the long-term evolution of stellar orbital parameters.

Dispersion function in complex frequency for the isotropic isochrone cluster.

Linear response theory and damped modes of stellar clusters


We present a new method to perform explicitly Landau's prescription in stellar systems, and compute the associated response matrix.

Time evolution of the clusters' core radius for various velocity anisotropies.

Non-resonant relaxation of anisotropic globular clusters


We consider globular clusters with anisotropic velocity distributions. We compare non-resonant kinetic theory with direct N-body measurements.

Illustration of VRR in a galactic nucleus.

Orbital alignment and mass segregation in galactic nuclei via vector resonant relaxation


We implement a maximum entropy method to predict the thermodynamical equilibrium distribution of stellar orientations around a supermassive BH.

Energy diffusion coefficients in the 1D self-gravitating system.

Long-term relaxation of 1D self-gravitating systems


We implement explicitly the (inhomogeneous) Balescu-Lenard equation in 1D self-gravitating systems. We perform quantitative comparisons with direct N-body measurements.

Diffusion flux in an inhomogeneous kinetically blocked spin system.

Kinetic theory of one-dimensional inhomogeneous long-range interacting N-body systems at order 1/N^2 without collective effects


In the limit where collective effects can be neglected, we derive a closed and explicit 1/N^2 collision operator for 1D inhomogeneous systems.

Diffusion processes around a supermassive BH.

Astrophysics with the Laser Interferometer Space Antenna


Review summarising the astrophysical processes relevant for the upcoming LISA interferometer.

Numerical integration of Kozai-Lidov oscillations using spinors

Milankovitch equations with spinors


We investigate the use of spinors to describe the secular evolution of quasi-Keplerian systems.

Illustration of the dilution of a stellar disc through vector resonant relaxation.

Constraining intermediate-mass black holes from the stellar disc of SgrA*


We use jointly kinetic predictions for the efficiency of dilution through vector resonant relaxation along with the recent observation of a stellar disc around SgrA*, to constrain the possible presence of intermediate mass black holes around SgrA*.

Dependence of the relaxation time as a function of the number of particles and the strength of the pairwise interaction.

Kinetic blockings in long-range interacting systems


We investigate the long-term relaxation of systems submitted to a second-order bare kinetic blocking.

Radial orbit instability in globular clusters.

Predicting the linear response of self-gravitating stellar spheres and discs with LinearResponse.jl


We present LinearResponse.jl, an efficient, versatile public library written in julia to compute the linear response of self-gravitating stellar spheres and razor-thin discs.

Large deviation principle for the system's empirical angle-averaged distribution function.

Dynamical large deviations for long-range interacting inhomogeneous systems without collective effects


We consider the long-term evolution of an inhomogeneous long-range interacting N -body system. Placing ourselves in the dynamically hot limit, i.e. neglecting collective effects, we derive a large deviation principle for the system's empirical angle-averaged distribution function.

Relaxation rate for various velocity anisotropies and strengths of rotation.

Non-resonant relaxation of rotating globular clusters


The long-term relaxation of rotating globular clusters is investigated through an extension of the orbit averaged Chandrasekhar non-resonant formalism.

Schematic illustration of typical orbits in (a) cold disks and (b) hot spheres.

Kinetic Theory of Stellar Systems: A Tutorial


This tutorial article is an introduction to the collective dynamical evolution of the very large numbers of stars and/or other self-gravitating objects that comprise such systems, i.e. their kinetic theory.