Showing 1–44 of 44 results for author: Fouvry, J

Vector Resonant Relaxation and Statistical Closure Theory. I. Direct Interaction Approximation

Authors: Sofia Flores, Jean-Baptiste Fouvry

Abstract: Stars orbiting a supermassive black hole in the center of galaxies undergo very efficient diffusion in their orbital orientations. This is "Vector Resonant Relaxation", a diffusion process formally occurring on the unit sphere. Such a dynamics is intrinsically non-linear, stochastic, and correlated, hence bearing deep similarities with turbulence in fluid mechanics or plasma physics. In that conte… ▽ More Stars orbiting a supermassive black hole in the center of galaxies undergo very efficient diffusion in their orbital orientations. This is "Vector Resonant Relaxation", a diffusion process formally occurring on the unit sphere. Such a dynamics is intrinsically non-linear, stochastic, and correlated, hence bearing deep similarities with turbulence in fluid mechanics or plasma physics. In that context, we show how generic methods stemming from statistical closure theory, namely the celebrated "Martin-Siggia-Rose formalism", can be used to characterize the correlations describing the redistribution of orbital orientations. In particular, limiting ourselves to the leading order truncation in this closure scheme, the so-called "Direct Interaction Approximation", and placing ourselves in the limit of an isotropic distribution of orientations, we explicitly compare the associated prediction for the two-point correlation function with measures from numerical simulations. We discuss the successes and limitations of this approach and present possible future venues. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 15 pages, 7 figures. Submitted to APS

Kinetic Theory of Stellar Systems: A Tutorial

Authors: Chris Hamilton, Jean-Baptiste Fouvry

Abstract: Stellar systems - star clusters, galaxies, dark matter haloes, and so on - are ubiquitous characters in the evolutionary tale of our Universe. 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. We begin by introducing the basic phenomenology… ▽ More Stellar systems - star clusters, galaxies, dark matter haloes, and so on - are ubiquitous characters in the evolutionary tale of our Universe. 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. We begin by introducing the basic phenomenology of stellar systems, and explaining why and when we must develop a kinetic theory that transcends the traditional two-body relaxation picture of Chandrasekhar. We study the orbits that comprise stellar systems, how those orbits are modified by perturbations, how a system responds self-consistently to fluctuations in its gravitational potential, and how one can predict the long term fate of a stellar system in various dynamical regimes. Though our treatment is necessarily mathematical, we develop the formalism only to the extent that it facilitates real calculations. We give many examples throughout the text of the equations being applied to topics of major astrophysical importance. Furthermore, in the 1960s and 1970s the kinetic theory of stellar systems was a fledgling subject which developed in tandem with the kinetic theory of plasmas. However, the two fields have long since diverged. Yet once one has become fluent in both Plasmaish and Galacticese, and has a dictionary relating the two, one can pull ideas directly from one field to solve a problem in the other. Therefore, another aim of this tutorial article is to provide our plasma colleagues with a jargon-light understanding of the key properties of stellar systems, to point out the many direct analogies between stellar- and plasma-kinetic calculations, and ultimately to convince them that stellar dynamics and plasma kinetics are, in a deep and beautiful and useful sense, the same thing. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 66 Pages, 34 Figures. Invited tutorial article for Physics of Plasmas. Comments/criticism strongly encouraged, especially from non-specialists!

Non-resonant relaxation of rotating globular clusters

Authors: Kerwann Tep, Jean-Baptiste Fouvry, Christophe Pichon

Abstract: The long-term relaxation of rotating, spherically symmetric globular clusters is investigated through an extension of the orbit averaged Chandrasekhar non-resonant formalism. A comparison is made with the long-term evolution of the distribution function in action space, measured from averages of sets of $N$-body simulations up to core collapse. The impact of rotation on in-plane relaxation is foun… ▽ More The long-term relaxation of rotating, spherically symmetric globular clusters is investigated through an extension of the orbit averaged Chandrasekhar non-resonant formalism. A comparison is made with the long-term evolution of the distribution function in action space, measured from averages of sets of $N$-body simulations up to core collapse. The impact of rotation on in-plane relaxation is found to be weak. In addition, we observe a clear match between theoretical predictions and $N$-body measurements. For the class of rotating models considered, we find no strong gravo-gyro catastrophe accelerating core collapse. Both kinetic theory and simulations predict a reshuffling of orbital inclinations from overpopulated regions to underpopulated ones. This trend accelerates as the amount of rotation is increased. Yet, for orbits closer to the rotational plane, the non-resonant prediction does not reproduce numerical measurements. We argue that this mismatch stems from these orbits' coherent interactions, which are not captured by the non-resonant formalism that only addresses local deflections. △ Less

Submitted 22 April, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: 18 pages, 18 figures, submitted to A&A

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

Authors: Michael S. Petersen, Mathieu Roule, Jean-Baptiste Fouvry, Christophe Pichon, Kerwann Tep

Abstract: We present LinearResponse.jl, an efficient, versatile public library written in julia to compute the linear response of self-gravitating (3D spherically symmetric) stellar spheres and (2D axisymmetric razor-thin) discs. LinearResponse.jl can scan the whole complex frequency plane, probing unstable, neutral and (weakly) damped modes. Given a potential model and a distribution function, this numeric… ▽ More We present LinearResponse.jl, an efficient, versatile public library written in julia to compute the linear response of self-gravitating (3D spherically symmetric) stellar spheres and (2D axisymmetric razor-thin) discs. LinearResponse.jl can scan the whole complex frequency plane, probing unstable, neutral and (weakly) damped modes. Given a potential model and a distribution function, this numerical toolbox estimates the modal frequencies as well as the shapes of individual modes. The libraries are validated against a combination of previous results for the spherical isochrone model and Mestel discs, and new simulations for the spherical Plummer model. Beyond linear response theory, the realm of applications of LinearResponse.jl also extends to the kinetic theory of self-gravitating systems through a modular interface. △ Less

Submitted 17 November, 2023; originally announced November 2023.

Comments: Software available at https://github.com/michael-petersen/LinearResponse.jl

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

Authors: Ouassim Feliachi, Jean-Baptiste Fouvry

Abstract: 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. This result extends the classical ensemble-averaged kinetic theory given by the so-called inhomogeneous Landau equati… ▽ More 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. This result extends the classical ensemble-averaged kinetic theory given by the so-called inhomogeneous Landau equation, as it specifies the probability of typical and large dynamical fluctuations. We detail the main properties of the associated large deviation Hamiltonian, in particular how it complies with the system's conservation laws and possesses a gradient structure. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: 13 pages, 0 figures, submitted to APS

Kinetic blockings in long-range interacting inhomogeneous systems

Authors: Jean-Baptiste Fouvry, Mathieu Roule

Abstract: Long-range interacting systems unavoidably relax through Poisson shot noise fluctuations generated by their finite number of particles, $N$. When driven by two-body correlations, i.e. ${1/N}$ effects, this long-term evolution is described by the inhomogeneous Balescu-Lenard equation. Yet, in one-dimensional systems with a monotonic frequency profile and only subject to 1:1 resonances, this kinetic… ▽ More Long-range interacting systems unavoidably relax through Poisson shot noise fluctuations generated by their finite number of particles, $N$. When driven by two-body correlations, i.e. ${1/N}$ effects, this long-term evolution is described by the inhomogeneous Balescu-Lenard equation. Yet, in one-dimensional systems with a monotonic frequency profile and only subject to 1:1 resonances, this kinetic equation exactly vanishes: this is a first-order full kinetic blocking. These systems' long-term evolution is then driven by three-body correlations, i.e. ${1/N^2}$ effects. In the limit of dynamically hot systems, this is described by the inhomogeneous ${1/N^2}$ Landau equation. We investigate numerically the long-term evolution of systems for which this second kinetic equation also exactly vanishes: this a second-order bare kinetic blocking. We demonstrate that these systems relax through the "leaking" contributions of dressed three-body interactions that are neglected in the inhomogeneous ${1/N^2}$ Landau equation. Finally, we argue that these never-vanishing contributions prevent four-body correlations, i.e. ${1/N^{3}}$ effects, from ever being the main driver of relaxation. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 14 pages, 8 figures, submitted to APS

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

Authors: Jean-Baptiste Fouvry, María José Bustamante-Rosell, Aaron Zimmerman

Abstract: Stars evolving around a supermassive black hole see their orbital orientations diffuse efficiently, a process called "vector resonant relaxation". In particular, stars within the same disc, i.e. neighbors in orientations, will slowly diffuse away from one another through this stochastic process. We use jointly (i) detailed kinetic predictions for the efficiency of this dilution and (ii) the recent… ▽ More Stars evolving around a supermassive black hole see their orbital orientations diffuse efficiently, a process called "vector resonant relaxation". In particular, stars within the same disc, i.e. neighbors in orientations, will slowly diffuse away from one another through this stochastic process. We use jointly (i) detailed kinetic predictions for the efficiency of this dilution and (ii) the recent observation of a stellar disc around SgrA*, the supermassive black hole at the centre of the Milky-Way, to constrain SgrA*'s unobserved stellar cluster. Notably, we investigate quantitatively the impact of a population of intermediate mass black holes on the survivability of the stellar disc. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 10 pages, 10 figures, submitted to MNRAS

Milankovitch equations with spinors

Authors: Barnabás Deme, Jean-Baptiste Fouvry

Abstract: We investigate the use of spinors to describe the secular evolution of quasi-Keplerian systems. Evaluating their Poisson brackets, we show that the components of a properly-chosen spinor are canonical variables. We illustrate this formalism with a satellite's motion around an oblate body. We investigate the use of spinors to describe the secular evolution of quasi-Keplerian systems. Evaluating their Poisson brackets, we show that the components of a properly-chosen spinor are canonical variables. We illustrate this formalism with a satellite's motion around an oblate body. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: 12 pages, no figures

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

Authors: Jean-Baptiste Fouvry

Abstract: Long-range interacting systems irreversibly relax as a result of their finite number of particles, $N$. At order $1/N$, this process is described by the inhomogeneous Balescu--Lenard equation. Yet, this equation exactly vanishes in one-dimensional inhomogeneous systems with a monotonic frequency profile and sustaining only 1:1 resonances. In the limit where collective effects can be neglected, we… ▽ More Long-range interacting systems irreversibly relax as a result of their finite number of particles, $N$. At order $1/N$, this process is described by the inhomogeneous Balescu--Lenard equation. Yet, this equation exactly vanishes in one-dimensional inhomogeneous systems with a monotonic frequency profile and sustaining only 1:1 resonances. In the limit where collective effects can be neglected, we derive a closed and explicit $1/N^{2}$ collision operator for such systems. We detail its properties highlighting in particular how it satisfies an $H$-theorem for Boltzmann entropy. We also compare its predictions with direct $N$-body simulations. Finally, we exhibit a generic class of long-range interaction potentials for which this $1/N^{2}$ collision operator exactly vanishes. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: 9 pages, 1 figure, submitted to APS

Long-term relaxation of ${1D}$ self-gravitating systems

Authors: Mathieu Roule, Jean-Baptiste Fouvry, Christophe Pichon, Pierre-Henri Chavanis

Abstract: We investigate the long-term relaxation of one-dimensional (${1D}$) self-gravitating systems, using both kinetic theory and $N$-body simulations. We consider thermal and Plummer equilibria, with and without collective effects. All combinations are found to be in clear agreement with respect to the Balescu-Lenard and Landau predictions for the diffusion coefficients. Interestingly, collective effec… ▽ More We investigate the long-term relaxation of one-dimensional (${1D}$) self-gravitating systems, using both kinetic theory and $N$-body simulations. We consider thermal and Plummer equilibria, with and without collective effects. All combinations are found to be in clear agreement with respect to the Balescu-Lenard and Landau predictions for the diffusion coefficients. Interestingly, collective effects reduce the diffusion by a factor ${\sim 10}$. The predicted flux for Plummer equilibrium matches the measured one, which is a remarkable validation of kinetic theory. We also report on a situation of quasi kinetic blocking for the same equilibrium. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 12 pages, 11 figures, submitted to APS

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Non-resonant relaxation of anisotropic globular clusters

Authors: Kerwann Tep, Jean-Baptiste Fouvry, Christophe Pichon

Abstract: Globular clusters are dense stellar systems whose core slowly contracts under the effect of self-gravity. The rate of this process was recently found to be directly linked to the initial amount of velocity anisotropy: tangentially anisotropic clusters contract faster than radially anisotropic ones. Furthermore, initially anisotropic clusters are found to generically tend towards more isotropic dis… ▽ More Globular clusters are dense stellar systems whose core slowly contracts under the effect of self-gravity. The rate of this process was recently found to be directly linked to the initial amount of velocity anisotropy: tangentially anisotropic clusters contract faster than radially anisotropic ones. Furthermore, initially anisotropic clusters are found to generically tend towards more isotropic distributions during the onset of contraction. Chandrasekhar's "non-resonant" (NR) theory of diffusion describes this relaxation as being driven by a sequence of local two-body deflections along each star's orbit. We explicitly tailor this NR prediction to anisotropic clusters, and compare it with $N$-body realisations of Plummer spheres with varying degrees of anisotropy. The NR theory is shown to recover remarkably well the detailed shape of the orbital diffusion and the associated initial isotropisation, up to a global multiplicative prefactor which increases with anisotropy. Strikingly, a simple effective isotropic prescription provides almost as good a fit, as long as the cluster's anisotropy is not too strong. For these more extreme clusters, accounting for long-range resonant relaxation may be necessary to capture these clusters' long-term evolution. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 16 pages, 18 figures, submitted to MNRAS

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

Authors: Nathan Magnan, Jean-Baptiste Fouvry, Christophe Pichon, Pierre-Henri Chavanis

Abstract: Supermassive black holes dominate the gravitational potential in galactic nuclei. In these dense environments, stars follow nearly Keplerian orbits and see their orbital planes relax through the potential fluctuations generated by the stellar cluster itself. For typical astrophysical galactic nuclei, the most likely outcome of this vector resonant relaxation (VRR) is that the orbital planes of the… ▽ More Supermassive black holes dominate the gravitational potential in galactic nuclei. In these dense environments, stars follow nearly Keplerian orbits and see their orbital planes relax through the potential fluctuations generated by the stellar cluster itself. For typical astrophysical galactic nuclei, the most likely outcome of this vector resonant relaxation (VRR) is that the orbital planes of the most massive stars spontaneously self-align within a narrow disc. We present a maximum entropy method to systematically determine this long-term distribution of orientations and use it for a wide range of stellar orbital parameters and initial conditions. The heaviest stellar objects are found to live within a thin equatorial disk. The thickness of this disk depends on the stars' initial mass function, and on the geometry of the initial cluster. This work highlights a possible (indirect) novel method to constrain the distribution of intermediate mass black holes in galactic nuclei. △ Less

Submitted 5 July, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 14 pages, 17 figures, accepted for publication in MNRAS

Journal ref: MNRAS 514, 3452 (2022)

Linear response theory and damped modes of stellar clusters

Authors: Jean-Baptiste Fouvry, Simon Prunet

Abstract: Because all stars contribute to its gravitational potential, stellar clusters amplify perturbations collectively. In the limit of small fluctuations, this is described through linear response theory, via the so-called response matrix. While the evaluation of this matrix is somewhat straightforward for unstable modes (i.e. with a positive growth rate), it requires a careful analytic continuation fo… ▽ More Because all stars contribute to its gravitational potential, stellar clusters amplify perturbations collectively. In the limit of small fluctuations, this is described through linear response theory, via the so-called response matrix. While the evaluation of this matrix is somewhat straightforward for unstable modes (i.e. with a positive growth rate), it requires a careful analytic continuation for damped modes (i.e. with a negative growth rate). We present a generic method to perform such a calculation in spherically symmetric stellar clusters. When applied to an isotropic isochrone cluster, we recover the presence of a low-frequency weakly damped $\ell = 1$ mode. We finally use a set of direct $N$-body simulations to test explicitly this prediction through the statistics of the correlated random walk undergone by a cluster's density centre. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 12 pages, 10 figures, submitted to MNRAS

Mapping the Galactic centre's dark cluster via Resonant Relaxation

Authors: Kerwann Tep, Jean-Baptiste Fouvry, Christophe Pichon, Gernot Heißel, Thibaut Paumard, Guy Perrin, Frederic Vincent

Abstract: Supermassive black holes in the centre of galaxies dominate the gravitational potential of their surrounding stellar clusters. In these dense environments, stars follow nearly Keplerian orbits, which get slowly distorted as a result of the potential fluctuations generated by the stellar cluster itself as a whole. In particular, stars undergo a rapid relaxation of their eccentricities through both… ▽ More Supermassive black holes in the centre of galaxies dominate the gravitational potential of their surrounding stellar clusters. In these dense environments, stars follow nearly Keplerian orbits, which get slowly distorted as a result of the potential fluctuations generated by the stellar cluster itself as a whole. In particular, stars undergo a rapid relaxation of their eccentricities through both resonant and non-resonant processes. An efficient implementation of the resonant diffusion coefficients allows for detailed and systematic explorations of the parameter space describing the properties of the stellar cluster. In conjunction with recent observations of the S-cluster orbiting SgrA*, this framework can be used to jointly constrain the distribution of the unresolved, old, background stellar cluster and the characteristics of a putative dark cluster. Specifically, we show how this can be used to estimate the typical mass and cuspide exponent of intermediate-mass black holes consistent with the relaxed state of the distribution of eccentricities in the observed S-cluster. This should prove useful in constraining super massive black hole formation scenarios. △ Less

Submitted 6 July, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 13 pages, 10 figures, submitted to MNRAS

Resonant and non-resonant relaxation of globular clusters

Authors: Jean-Baptiste Fouvry, Chris Hamilton, Simon Rozier, Christophe Pichon

Abstract: Globular clusters contain a finite number of stars. As a result, they inevitably undergo secular evolution (`relaxation') causing their mean distribution function (DF) to evolve on long timescales. On one hand, this long-term evolution may be interpreted as driven by the accumulation of local deflections along each star's mean field trajectory -- so-called `non-resonant relaxation'. On the other h… ▽ More Globular clusters contain a finite number of stars. As a result, they inevitably undergo secular evolution (`relaxation') causing their mean distribution function (DF) to evolve on long timescales. On one hand, this long-term evolution may be interpreted as driven by the accumulation of local deflections along each star's mean field trajectory -- so-called `non-resonant relaxation'. On the other hand, it can be thought of as driven by non-local, collectively dressed and resonant couplings between stellar orbits, a process termed `resonant relaxation'. In this paper we consider a model globular cluster represented by a spherical, isotropic isochrone DF, and compare in detail the predictions of both resonant and non-resonant relaxation theories against tailored direct $N$-body simulations. In the space of orbital actions (namely the radial action and total angular momentum), we find that both resonant and non-resonant theories predict the correct morphology for the secular evolution of the cluster's DF, although non-resonant theory over-estimates the amplitude of the relaxation rate by a factor ${\sim 2}$. We conclude that the secular relaxation of hot isotropic spherical clusters is not dominated by collectively amplified large-scale potential fluctuations, despite the existence of a strong ${\ell = 1}$ damped mode. Instead, collective amplification affects relaxation only marginally even on the largest scales. The predicted contributions to relaxation from smaller scale fluctuations are essentially the same from resonant and non-resonant theories. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Comments: 26 pages, 14 figures, submitted to MNRAS

Secular Dynamics around a Supermassive Black Hole via Multipole Expansion

Authors: Jean-Baptiste Fouvry, Walter Dehnen, Scott Tremaine, Ben Bar-Or

Abstract: In galactic nuclei, the gravitational potential is dominated by the central supermassive black hole, so stars follow quasi-Keplerian orbits. These orbits are distorted by gravitational forces from other stars, leading to long-term orbital relaxation. The direct numerical study of these processes is challenging because the fast orbital motion imposed by the central black hole requires very small ti… ▽ More In galactic nuclei, the gravitational potential is dominated by the central supermassive black hole, so stars follow quasi-Keplerian orbits. These orbits are distorted by gravitational forces from other stars, leading to long-term orbital relaxation. The direct numerical study of these processes is challenging because the fast orbital motion imposed by the central black hole requires very small timesteps. An alternative approach, pioneered by Gauss, is to use the secular approximation of smearing out the $N$ stars over their Keplerian orbits, using $K$ nodes along each orbit. In this study we propose three novel improvements to this method. First, we re-formulate the discretisation of the rates of change of the variables describing the orbital states to ensure that all conservation laws are exactly satisfied. Second, we replace the pairwise sum over nodes by a multipole expansion up to order $\ell_{\mathrm{max}}$, reducing the overall computational costs from $O(N^2K^2)$ to $O(NK\ell_{\mathrm{max}}^2)$. Finally, we show that the averaged dynamical system is equivalent to $2N$ interacting unit spin vectors and provide two time integrators: a second-order symplectic scheme and a fourth-order Lie-group Runge-Kutta method, both of which are straightforward to generalize to higher order. These new simulations recover the diffusion coefficients of stellar eccentricities obtained through analytical calculations of the secular dynamics. △ Less

Submitted 1 March, 2022; v1 submitted 3 November, 2020; originally announced November 2020.

Comments: 20 pages, 11 figures, submitted to AAS

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

Authors: Ben Bar-Or, Jean-Baptiste Fouvry, Scott Tremaine

Abstract: Fuzzy dark matter (FDM) is composed of ultra-light bosons having a de Broglie wavelength that is comparable to the size of the stellar component of galaxies at typical galactic velocities. FDM behaves like cold dark matter on large scales. However, on the scale of the de Broglie wavelength, an FDM halo exhibits density fluctuations that lead to relaxation, a process similar to the two-body relaxat… ▽ More Fuzzy dark matter (FDM) is composed of ultra-light bosons having a de Broglie wavelength that is comparable to the size of the stellar component of galaxies at typical galactic velocities. FDM behaves like cold dark matter on large scales. However, on the scale of the de Broglie wavelength, an FDM halo exhibits density fluctuations that lead to relaxation, a process similar to the two-body relaxation that occurs in classical gravitational N-body systems and is described by the Fokker-Planck equation. We derive the FDM analog of that kinetic equation, and solve it to find the evolution of the velocity distribution in a spatially homogeneous FDM halo. We also determine the dielectric function and the dispersion relation for linear waves in an FDM halo. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: 22 pages, 4 figures, submitted to AAS

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

Authors: Juan Giral Martínez, Jean-Baptiste Fouvry, Christophe Pichon

Abstract: We investigate the rate of orbital orientation dilution of young stellar clusters in the vicinity of supermassive black holes. Within the framework of vector resonant relaxation, we predict the time evolution of the two-point correlation function of the stellar orbital plane orientations as a function of their initial angular separation and diversity in orbital parameters (semi-major axis, eccentr… ▽ More We investigate the rate of orbital orientation dilution of young stellar clusters in the vicinity of supermassive black holes. Within the framework of vector resonant relaxation, we predict the time evolution of the two-point correlation function of the stellar orbital plane orientations as a function of their initial angular separation and diversity in orbital parameters (semi-major axis, eccentricity). As expected, the larger the spread in initial orientations and orbital parameters, the more efficient the dilution of a given set of co-eval stars, with a characteristic timescale set up by the coherence time of the background potential fluctuations. A Markovian prescription which matches numerical simulations allows us to efficiently probe the underlying kinematic properties of the unresolved nucleus when requesting consistency with a given dilution efficiency, imposed by the observed stellar disc within the one arcsecond of Sgr A*. As a proof of concept, we compute maps of constant dilution times as a function of the semi major axis cusp index and fraction of intermediate mass black holes in the old background stellar cluster. This computation suggests that vector resonant relaxation should prove useful in this context since it impacts orientations on timescales comparable to the stars' age. △ Less

Submitted 5 August, 2020; originally announced August 2020.

Comments: 20 pages, 16 figures, submitted to MNRAS

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

Authors: Jean-Baptiste Fouvry, Pierre-Henri Chavanis, Christophe Pichon

Abstract: Finite-$N$ effects unavoidably drive the long-term evolution of long-range interacting $N$-body systems. The Balescu-Lenard kinetic equation generically describes this process sourced by ${1/N}$ effects but this kinetic operator exactly vanishes by symmetry for one-dimensional homogeneous systems: such systems undergo a kinetic blocking and cannot relax as a whole at this order in ${1/N}$. It is t… ▽ More Finite-$N$ effects unavoidably drive the long-term evolution of long-range interacting $N$-body systems. The Balescu-Lenard kinetic equation generically describes this process sourced by ${1/N}$ effects but this kinetic operator exactly vanishes by symmetry for one-dimensional homogeneous systems: such systems undergo a kinetic blocking and cannot relax as a whole at this order in ${1/N}$. It is therefore only through the much weaker ${1/N^{2}}$ effects, sourced by three-body correlations, that these systems can relax, leading to a much slower evolution. In the limit where collective effects can be neglected, but for an arbitrary pairwise interaction potential, we derive a closed and explicit kinetic equation describing this very long-term evolution. We show how this kinetic equation satisfies an $H$-theorem while conserving particle number and energy, ensuring the unavoidable relaxation of the system towards the Boltzmann equilibrium distribution. Provided that the interaction is long-range, we also show how this equation cannot suffer from further kinetic blocking, i.e., the ${1/N^{2}}$ dynamics is always effective. Finally, we illustrate how this equation quantitatively matches measurements from direct $N$-body simulations. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 15 pages, 3 figures, submitted to Phys. Rev. E

Journal ref: Phys. Rev. E 102, 052110 (2020)

The Strength of the Dynamical Spiral Perturbation in the Galactic Disk

Authors: Anna-Christina Eilers, David W. Hogg, Hans-Walter Rix, Neige Frankel, Jason A. S. Hunt, Jean-Baptiste Fouvry, Tobias Buck

Abstract: The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\varphi)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, an… ▽ More The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\varphi)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, and find a good qualitative and quantitative match. Presuming that the amplitude of the gravitational potential perturbation is proportionate to that in the disk's surface mass density, we estimate the surface mass density amplitude to be $Σ_{\rm max} (R_{\odot})\approx 5.5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius when choosing a fixed pattern speed of $Ω_{\mathrm p}=12\,\rm km\,s^{-1}\,kpc^{-1}$. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 10\%$ at the solar radius, with an increases towards larger radii. Our model constrains the pitch angle of the dynamical spiral arms to be approximately $12^{\circ}$. △ Less

Submitted 3 August, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: Accepted for publication in ApJ

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

Authors: Jean-Baptiste Fouvry, Ben Bar-Or, Pierre-Henri Chavanis

Abstract: The long-term dynamics of long-range interacting $N$-body systems can generically be described by the Balescu-Lenard kinetic equation. However, for ${1D}$ homogeneous systems, this collision operator exactly vanishes by symmetry. These systems undergo a kinetic blocking, and cannot relax as a whole under ${1/N}$ resonant effects. As a result, these systems can only relax under ${1/N^{2}}$ effects,… ▽ More The long-term dynamics of long-range interacting $N$-body systems can generically be described by the Balescu-Lenard kinetic equation. However, for ${1D}$ homogeneous systems, this collision operator exactly vanishes by symmetry. These systems undergo a kinetic blocking, and cannot relax as a whole under ${1/N}$ resonant effects. As a result, these systems can only relax under ${1/N^{2}}$ effects, and their relaxation is drastically slowed down. In the context of the homogeneous Hamiltonian Mean Field model, we present a new, closed and explicit kinetic equation describing self-consistently the very long-term evolution of such systems, in the limit where collective effects can be neglected, i.e. for dynamically hot initial conditions. We show in particular how that kinetic equation satisfies an $H$-Theorem that guarantees the unavoidable relaxation to the Boltzmann equilibrium distribution. Finally, we illustrate how that kinetic equation quantitatively matches with the measurements from direct $N$-body simulations. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: 15 pages, 4 figures, submitted to PRE

Journal ref: Phys. Rev. E 100, 052142 (2019)

Mapping the stability of stellar rotating spheres via linear response theory

Authors: Simon Rozier, Jean-Baptiste Fouvry, Philip G. Breen, Anna Lisa Varri, Christophe Pichon, Douglas C. Heggie

Abstract: Rotation is ubiquitous in the Universe, and recent kinematic surveys have shown that early type galaxies and globular clusters are no exception. Yet the linear response of spheroidal rotating stellar systems has seldom been studied. This paper takes a step in this direction by considering the behaviour of spherically symmetric systems with differential rotation. Specifically, the stability of seve… ▽ More Rotation is ubiquitous in the Universe, and recent kinematic surveys have shown that early type galaxies and globular clusters are no exception. Yet the linear response of spheroidal rotating stellar systems has seldom been studied. This paper takes a step in this direction by considering the behaviour of spherically symmetric systems with differential rotation. Specifically, the stability of several sequences of Plummer spheres is investigated, in which the total angular momentum, as well as the degree and flavour of anisotropy in the velocity space are varied. To that end, the response matrix method is customised to spherical rotating equilibria. The shapes, pattern speeds and growth rates of the systems' unstable modes are computed. Detailed comparisons to appropriate N-body measurements are also presented. The marginal stability boundary is charted in the parameter space of velocity anisotropy and rotation rate. When rotation is introduced, two sequences of growing modes are identified corresponding to radially and tangentially-biased anisotropic spheres respectively. For radially anisotropic spheres, growing modes occur on two intersecting surfaces (in the parameter space of anisotropy and rotation), which correspond to fast and slow modes, depending on the net rotation rate. Generalised, approximate stability criteria are finally presented. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 18 pages, 20 figures. Submitted to MNRAS

Vector Resonant Relaxation of Stars around a Massive Black Hole

Authors: Jean-Baptiste Fouvry, Ben Bar-Or, Pierre-Henri Chavanis

Abstract: In the vicinity of a massive black hole, stars move on precessing Keplerian orbits. The mutual stochastic gravitational torques between the stellar orbits drive a rapid reorientation of their orbital planes, through a process called vector resonant relaxation. We derive, from first principles, the correlation of the potential fluctuations in such a system, and the statistical properties of random… ▽ More In the vicinity of a massive black hole, stars move on precessing Keplerian orbits. The mutual stochastic gravitational torques between the stellar orbits drive a rapid reorientation of their orbital planes, through a process called vector resonant relaxation. We derive, from first principles, the correlation of the potential fluctuations in such a system, and the statistical properties of random walks undergone by the stellar orbital orientations. We compare this new analytical approach with effective $N$-body simulations. We also provide a simple scheme to generate the random walk of a test star's orbital orientation using a stochastic equation of motion. We finally present quantitative estimations of this process for a nuclear stellar cluster such as the one of the Milky Way. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 23 pages, 10 figures, submitted to AAS journals

Instabilities in disc galaxies: from noise to grooves to spirals

Authors: Sven De Rijcke, Jean-Baptiste Fouvry, Christophe Pichon

Abstract: Using the linearized Boltzmann equation, we investigate how grooves carved in the phase space of a half-mass Mestel disc can trigger the vigorous growth of two-armed spiral eigenmodes. Such grooves result from the collisional dynamics of a disc subject to finite-N shot noise, as swing-amplified noise patterns push stars towards lower-angular momentum orbits at their inner Lindblad radius. Suppleme… ▽ More Using the linearized Boltzmann equation, we investigate how grooves carved in the phase space of a half-mass Mestel disc can trigger the vigorous growth of two-armed spiral eigenmodes. Such grooves result from the collisional dynamics of a disc subject to finite-N shot noise, as swing-amplified noise patterns push stars towards lower-angular momentum orbits at their inner Lindblad radius. Supplementing the linear theory with analytical arguments, we show that the dominant spiral mode is a cavity mode with reflections off the forbidden region around corotation and off the deepest groove. Other subdominant modes are identified as groove modes. We provide evidence that the depletion of near-circular orbits, and not the addition of radial orbits, is the crucial physical ingredient that causes these new eigenmodes. Thus, it is possible for an isolated, linearly stable stellar disc to spontaneously become linearly unstable via the self-induced formation of phase-space grooves through finite-N dynamics. These results may help explain the growth and maintenance of spiral patterns in real disc galaxies. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 11 pages, 7 figures, submitted to MNRAS

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

Authors: Sven De Rijcke, Jean-Baptiste Fouvry, Walter Dehnen

Abstract: Linear perturbation is used to investigate the effect of gravitational softening on the retrieved two-armed spiral eigenmodes of razor-thin stellar discs. We explore four softening kernels with different degrees of gravity bias, and with/without compact support (compact in the sense that they yield exactly Newtonian forces outside the softening kernel). These kernels are applied to two disc gala… ▽ More Linear perturbation is used to investigate the effect of gravitational softening on the retrieved two-armed spiral eigenmodes of razor-thin stellar discs. We explore four softening kernels with different degrees of gravity bias, and with/without compact support (compact in the sense that they yield exactly Newtonian forces outside the softening kernel). These kernels are applied to two disc galaxy models with well-known unsoftened unstable modes. We illustrate quantitatively the importance of a vanishing linear gravity bias to yield accurate frequency estimates of the unstable modes. As such, Plummer softening, while very popular amongst simulators, performs poorly in our tests. The best results, with excellent agreement between the softened and unsoftened mode properties, are obtained with softening kernels that have a reduced gravity bias, obtained by compensating for the sub-Newtonian forces at small interparticle distances with slightly super-Newtonian forces at radii near the softening length. We present examples of such kernels that, moreover, are analytically simple and computationally cheap. Finally, these results light the way to the construction of softening methods with even smaller gravity bias, although at the price of increasingly complex kernels. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 14 pages, 7 figures, submitted to MNRAS

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

Authors: Jean-Baptiste Fouvry, Ben Bar-Or, Pierre-Henri Chavanis

Abstract: We investigate the secular dynamics of long-range interacting particles moving on a sphere, in the limit of an axisymmetric mean field potential. We show that this system can be described by the general kinetic equation, the inhomogeneous Balescu--Lenard equation. We use this approach to compute long-term diffusion coefficients, that are compared with direct simulations. Finally, we show how the s… ▽ More We investigate the secular dynamics of long-range interacting particles moving on a sphere, in the limit of an axisymmetric mean field potential. We show that this system can be described by the general kinetic equation, the inhomogeneous Balescu--Lenard equation. We use this approach to compute long-term diffusion coefficients, that are compared with direct simulations. Finally, we show how the scaling of the system's relaxation rate with the number of particles fundamentally depends on the underlying frequency profile. This clarifies why systems with a monotonic profile undergo a kinetic blocking and cannot relax as a whole under ${1/N}$ resonant effects. Because of its general form, this framework can describe the dynamics of globally coupled classical Heisenberg spins, long-range couplings in liquid crystals, or the orbital inclination evolution of stars in nearly Keplerian systems. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: 8 pages, 6 figures, submitted to PRL

Journal ref: Phys. Rev. E 99, 032101 (2019)

Relaxation in a Fuzzy Dark Matter Halo

Authors: Ben Bar-Or, Jean-Baptiste Fouvry, Scott Tremaine

Abstract: Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $λ\sim 1 \,\mathrm{kpc}$ in typical galactic potentials. Such `fuzzy' dark matter (FDM) behaves like cold dark matter (CDM) on much larger scales than the de Broglie wavelength, but may resolve some of the challenges faced by CDM in explaining the properties of galaxies on small scales (… ▽ More Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $λ\sim 1 \,\mathrm{kpc}$ in typical galactic potentials. Such `fuzzy' dark matter (FDM) behaves like cold dark matter (CDM) on much larger scales than the de Broglie wavelength, but may resolve some of the challenges faced by CDM in explaining the properties of galaxies on small scales ($\lesssim 10\,\mathrm{kpc}$). Because of its wave nature, FDM exhibits stochastic density fluctuations on the scale of the de Broglie wavelength that never damp. The gravitational field from these fluctuations scatters stars and black holes, causing their orbits to diffuse through phase space. We show that this relaxation process can be analyzed quantitatively with the same tools used to analyze classical two-body relaxation in an $N$-body system, and can be described by treating the FDM fluctuations as quasiparticles, with effective mass $\sim 10^7 M_\odot {(1\,\mathrm{kpc}/r)}^2{(10^{-22}\,\mathrm{eV}/m_b)}^3$ in a galaxy with a constant circular speed of $200\,\mathrm{kms}$. This novel relaxation mechanism may stall the inspiral of supermassive black holes or globular clusters due to dynamical friction at radii of a few hundred pc, and can heat and expand the central regions of galaxies. These processes can be used to constrain the mass of the light bosons that might comprise FDM. △ Less

Submitted 25 January, 2019; v1 submitted 20 September, 2018; originally announced September 2018.

Comments: 23 pages, 7 figures, published on APJ

Journal ref: The Astrophysical Journal, 871:28, 2019

Revisiting relaxation in globular clusters

Authors: Chris Hamilton, Jean-Baptiste Fouvry, James Binney, Christophe Pichon

Abstract: The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the cluster's self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-parti… ▽ More The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the cluster's self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-particle scattering and is enhanced by self-gravity. The BL equation is adapted to a spherical system and used to estimate the flux through the action space of isochrone clusters with different velocity anisotropies. A range of fairly different secular behaviours is found depending on the fraction of radial orbits. Classical theory is also used to compute the corresponding classical fluxes. The BL and classical fluxes are very different because (a) the classical theory materially under-estimates the impact of large-scale collectively amplified fluctuations and (b) only the leading terms in an infinite sum for the BL flux are computed. A complete theory of cluster relaxation likely requires that the sum in the BL equation be decomposed into a sum over a finite number of small wavenumbers complemented by an integral over large wavenumbers analogous to classical theory. △ Less

Submitted 31 August, 2018; v1 submitted 29 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS. Replaced with updated manuscript including responses to referee's comments

Scalar Resonant Relaxation of Stars Around a Massive Black Hole

Authors: Ben Bar-Or, Jean-Baptiste Fouvry

Abstract: In nuclear star clusters, the potential is governed by the central massive black hole, so that stars move on nearly Keplerian orbits and the total potential is almost stationary in time. Yet, the deviations of the potential from the Keplerian one, due to the enclosed stellar mass and general relativity, will cause the stellar orbits to precess. Moreover, as a result of the finite number of stars,… ▽ More In nuclear star clusters, the potential is governed by the central massive black hole, so that stars move on nearly Keplerian orbits and the total potential is almost stationary in time. Yet, the deviations of the potential from the Keplerian one, due to the enclosed stellar mass and general relativity, will cause the stellar orbits to precess. Moreover, as a result of the finite number of stars, small deviations of the potential from spherical symmetry induce residual torques that can change the stars' angular momentum faster than the standard two-body relaxation. The combination of these two effects drives a stochastic evolution of orbital angular momentum, a process named "resonant relaxation". Owing to recent developments in the description of the relaxation of self-gravitating systems, we can now fully describe scalar resonant relaxation (relaxation of the magnitude of the angular momentum) as a diffusion process. In this framework, the potential fluctuations due to the complex orbital motion of the stars are described by a random correlated noise, whose statistical properties are fully characterized by the stars' mean field motion. On long timescales, the cluster can be regarded as a diffusive system, whose diffusion coefficients depend explicitly on the mean field stellar distribution through the properties of the noise. We show here, for the first time, how the diffusion coefficients of scalar resonant relaxation, for a spherically symmetric system, can be fully calculated from first principles, without any free parameters. We also provide an open source code that evaluates these diffusion coefficients numerically. △ Less

Submitted 24 February, 2018; originally announced February 2018.

Comments: 8 pages, 4 figures

Relaxation in self-gravitating systems

Authors: Jean-Baptiste Fouvry, Ben Bar-Or

Abstract: The long timescale evolution of a self-gravitating system is generically driven by two-body encounters. In many cases, the motion of the particles is primarily governed by the mean field potential. When this potential is integrable, particles move on nearly fixed orbits, which can be described in terms of angle-action variables. The mean field potential drives fast orbital motions (angles) whose a… ▽ More The long timescale evolution of a self-gravitating system is generically driven by two-body encounters. In many cases, the motion of the particles is primarily governed by the mean field potential. When this potential is integrable, particles move on nearly fixed orbits, which can be described in terms of angle-action variables. The mean field potential drives fast orbital motions (angles) whose associated orbits (actions) are adiabatically conserved on short dynamical timescales. The long-term stochastic evolution of the actions is driven by the potential fluctuations around the mean field and in particular by "resonant two-body encounters", for which the angular frequencies of two particles are in resonance. We show that the stochastic gravitational fluctuations acting on the particles can generically be described by a correlated Gaussian noise. Using this approach, the so-called $η$-formalism, we derive a diffusion equation for the actions in the test particle limit. We show that in the appropriate limits, this diffusion equation is equivalent to the inhomogeneous Balescu-Lenard and Landau equations. This approach provides a new view of the resonant diffusion processes associated with long-term orbital distortions. Finally, by investigating the example of the Hamiltonian Mean Field Model, we show how the present method generically allows for alternative calculations of the long-term diffusion coefficients in inhomogeneous systems. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Comments: 28 pages, 3 figures, submitted to MNRAS

Distribution functions for resonantly trapped orbits in the Galactic disc

Authors: Giacomo Monari, Benoit Famaey, Jean-Baptiste Fouvry, James Binney

Abstract: The present-day response of a Galactic disc stellar population to a non-axisymmetric perturbation of the potential has previously been computed through perturbation theory within the phase-space coordinates of the unperturbed axisymmetric system. Such an Eulerian linearized treatment however leads to singularities at resonances, which prevent quantitative comparisons with data. Here, we manage to… ▽ More The present-day response of a Galactic disc stellar population to a non-axisymmetric perturbation of the potential has previously been computed through perturbation theory within the phase-space coordinates of the unperturbed axisymmetric system. Such an Eulerian linearized treatment however leads to singularities at resonances, which prevent quantitative comparisons with data. Here, we manage to capture the behaviour of the distribution function (DF) at a resonance in a Lagrangian approach, by averaging the Hamiltonian over fast angle variables and re-expressing the DF in terms of a new set of canonical actions and angles variables valid in the resonant region. We then follow the prescription of Binney (2016), assigning to the resonant DF the time average along the orbits of the axisymmetric DF expressed in the new set of actions and angles. This boils down to phase-mixing the DF in terms of the new angles, such that the DF for trapped orbits only depends on the new set of actions. This opens the way to quantitatively fitting the effects of the bar and spirals to Gaia data in terms of distribution functions in action space. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 8 pages, 7 figures, accepted for publication in MNRAS

Dressed diffusion and friction coefficients in inhomogeneous multicomponent self-gravitating systems

Authors: Jean Heyvaerts, Jean-Baptiste Fouvry, Pierre-Henri Chavanis, Christophe Pichon

Abstract: General self-consistent expressions for the coefficients of diffusion and dynamical friction in a stable, bound, multicomponent self-gravitating and inhomogeneous system are derived. They account for the detailed dynamics of the colliding particles and their self-consistent dressing by collective gravitational interactions. The associated Fokker-Planck equation is shown to be fully consistent with… ▽ More General self-consistent expressions for the coefficients of diffusion and dynamical friction in a stable, bound, multicomponent self-gravitating and inhomogeneous system are derived. They account for the detailed dynamics of the colliding particles and their self-consistent dressing by collective gravitational interactions. The associated Fokker-Planck equation is shown to be fully consistent with the corresponding inhomogeneous Balescu-Lenard equation and, in the weak self-gravitating limit, to the inhomogeneous Landau equation. Hence it provides an alternative derivation to both and demonstrates their equivalence. The corresponding stochastic Langevin equations are presented: they can be a practical alternative to numerically solving the inhomogeneous Fokker-Planck and Balescu-Lenard equations. The present formalism allows for a self-consistent description of the secular evolution of different populations covering a spectrum of masses, with a proper accounting of the induced secular mass segregation, which should be of interest to various astrophysical contexts, from galactic centers to protostellar discs. △ Less

Submitted 19 June, 2017; originally announced June 2017.

Comments: 27 pages, 1 figure

Journal ref: MNRAS (2017) 469 (4): 4193-4220

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

Authors: Jean-Baptiste Fouvry, Christophe Pichon, Pierre-Henri Chavanis

Abstract: The drift and diffusion coefficients of the inhomogeneous multi-mass degenerate Landau equation are computed to describe the self-induced resonant relaxation of a discrete self-gravitating quasi-Keplerian razor-thin axisymmetric disc orbiting a massive black hole while relying on Gauss' method. For a disc-like configuration in our Galactic centre, secular diffusion induces an adiabatic distortion… ▽ More The drift and diffusion coefficients of the inhomogeneous multi-mass degenerate Landau equation are computed to describe the self-induced resonant relaxation of a discrete self-gravitating quasi-Keplerian razor-thin axisymmetric disc orbiting a massive black hole while relying on Gauss' method. For a disc-like configuration in our Galactic centre, secular diffusion induces an adiabatic distortion of orbits. When considering a disc composed of multiple masses similarly distributed, the population of lighter stars will gain eccentricity, driving it closer to the central black hole provided the distribution function increases with angular momentum. The quenching of the diffusion of a test star in the vicinity of the black hole due to the divergence of the relativistic precessions (the "Schwarzschild barrier") is correctly recovered by the kinetic equation. The dual stochastic Langevin formulation yields consistent results and provides a versatile framework in which to incorporate other stochastic processes. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: 16 pages, 14 figures, submitted to A&A

Journal ref: A&A 609, A38 (2018)

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

Authors: Jean-Baptiste Fouvry, Christophe Pichon, John Magorrian

Abstract: 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. Because the particles move in a quasi-Keplerian potential, their orbits can be approximated by ellipses whose orientations remain fixed over many dynamical times. The… ▽ More 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. Because the particles move in a quasi-Keplerian potential, their orbits can be approximated by ellipses whose orientations remain fixed over many dynamical times. The kinetic equation is obtained by simply averaging the BBGKY equations over the fast angle that describes motion along these ellipses. This so-called Balescu-Lenard equation describes self-consistently the long-term evolution of the distribution of quasi-Keplerian orbits around the central object: it models the diffusion and drift of their actions, induced through their mutual resonant interaction. Hence, it is the master equation that describes the secular effects of resonant relaxation. We show how it captures the phenonema of mass segregation and of the relativistic Schwarzschild barrier recently discovered in $N$-body simulations. △ Less

Submitted 22 February, 2017; v1 submitted 17 June, 2016; originally announced June 2016.

Comments: 24 pages, 3 figures

Journal ref: A&A 598, A71 (2017)

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

Authors: Jean-Baptiste Fouvry, Pierre-Henri Chavanis, Christophe Pichon

Abstract: 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. We start from Liouville's equation which describes the exact dynamics of a two-dimensional system of point vortices. At the order ${1/N}$, the evolution of the system is characterised by the first two equations of… ▽ More 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. We start from Liouville's equation which describes the exact dynamics of a two-dimensional system of point vortices. At the order ${1/N}$, the evolution of the system is characterised by the first two equations of the BBGKY hierarchy involving the system's 1-body distribution function and its 1-body correlation function. Thanks to the introduction of auxiliary fields, these two evolution constraints may be rewritten as a functional integral. When functionally integrated over the 2-body correlation function, this rewriting leads to a new constraint coupling the 1-body distribution function and the two auxiliary fields. Once inverted, this constraint provides, through a new route, the closed non-linear kinetic equation satisfied by the 1-body distribution function. Such a method sheds new lights on the origin of these kinetic equations complementing the traditional derivation methods. △ Less

Submitted 25 November, 2017; v1 submitted 13 June, 2016; originally announced June 2016.

Comments: 7 pages, 1 figure

Journal ref: Physics Letters A 380, 2589 (2016)

Secular diffusion in discrete self-gravitating tepid discs III. Resonant thickening in the tightly wound limit

Authors: Jean-Baptiste Fouvry, Christophe Pichon, Pierre-Henri Chavanis

Abstract: The secular thickening of a discrete self-gravitating galactic disc is investigated using the inhomogeneous multi-component Balescu-Lenard equation. The thick WKB limit for the diffusion and drift coefficients is found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields a simple double quadrature for the drift… ▽ More The secular thickening of a discrete self-gravitating galactic disc is investigated using the inhomogeneous multi-component Balescu-Lenard equation. The thick WKB limit for the diffusion and drift coefficients is found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields a simple double quadrature for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc induced by the effects of a finite number of particles. When applied to a tepid stable tapered disc, the Balescu-Lenard formalism predicts 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 amplication is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. The joint evolution of a population of giant molecular clouds within the disc may accelerate the secular disc's thickening induced by finite${-N}$ effects, but the observed number of clouds in the Milky Way does not seem to be sufficient to explain its thick disc. △ Less

Submitted 24 August, 2017; v1 submitted 11 May, 2016; originally announced May 2016.

Comments: This paper has been withdrawn by the authors. This paper has been superseded by arXiv:1605.03379

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

Authors: Jean-Baptiste Fouvry, Christophe Pichon, Pierre-Henri Chavanis, Laura Monk

Abstract: The secular thickening of a self-gravitating stellar galactic disc is investigated using the dressed collisionless Fokker-Planck equation and the inhomogeneous multicomponent Balescu-Lenard equation. The thick WKB limits for the diffusion fluxes are found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields simp… ▽ More The secular thickening of a self-gravitating stellar galactic disc is investigated using the dressed collisionless Fokker-Planck equation and the inhomogeneous multicomponent Balescu-Lenard equation. The thick WKB limits for the diffusion fluxes are found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields simple quadratures for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc, induced by fluctuations either external or due to the finite number of particles. These thick limits also offer a consistent derivation of a thick disc Toomre parameter, which is shown to be exponentially boosted by the ratio of the vertical to radial scale heights. Dressed potential fluctuations within the disc statistically induce a vertical bending of a subset of resonant orbits, triggering the corresponding increase in vertical velocity dispersion. When applied to a tepid stable tapered disc perturbed by shot noise, these two frameworks reproduce qualitatively the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but overestimates the time-scale involved in their appearance. Swing amplification is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. Other sources of thickening are also investigated, such as fading sequences of slowing bars, or the joint evolution of a population of giant molecular clouds within the disc. △ Less

Submitted 22 August, 2017; v1 submitted 11 May, 2016; originally announced May 2016.

Comments: 31 pages, 19 figures

Journal ref: MNRAS (2017) 471 (3): 2642-2673

Functional integral approach to the kinetic theory of inhomogeneous systems

Authors: Jean-Baptiste Fouvry, Pierre-Henri Chavanis, Christophe Pichon

Abstract: 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. We start from the BBGKY hierarchy derived from the Liouville equation. At the order ${1/N}$, where $N$ is the number of particles, the evolution of the system is… ▽ More 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. We start from the BBGKY hierarchy derived from the Liouville equation. At the order ${1/N}$, where $N$ is the number of particles, the evolution of the system is characterised by its 1-body distribution function and its 2-body correlation function. Introducing associated auxiliary fields, the evolution of these quantities may be rewritten as a traditional functional integral. By functionally integrating over the 2-body autocorrelation, one obtains a new constraint connecting the 1-body DF and the auxiliary fields. When inverted, this constraint allows us to obtain the closed non-linear kinetic equation satisfied by the 1-body distribution function. This derivation provides an alternative to previous methods, either based on the direct resolution of the truncated BBGKY hierarchy or on the Klimontovich equation. It may turn out to be fruitful to derive more accurate kinetic equations, e.g., accounting for collective effects, or higher order correlation terms. △ Less

Submitted 25 November, 2017; v1 submitted 8 April, 2016; originally announced April 2016.

Comments: 13 pages

Journal ref: Physica A 459, 117 (2016)

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

Authors: Jean-Baptiste Fouvry, Christophe Pichon, John Magorrian, Pierre-Henri Chavanis

Abstract: The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is investigated using the inhomogeneous Balescu-Lenard equation expressed in terms of angle-action variables. The matrix method is implemented numerically in order to model the induced gravitational polarization. Special care is taken to account for the amplification of potential fluctuati… ▽ More The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is investigated using the inhomogeneous Balescu-Lenard equation expressed in terms of angle-action variables. The matrix method is implemented numerically in order to model the induced gravitational polarization. Special care is taken to account for the amplification of potential fluctuations of mutually resonant orbits and the unwinding of the induced swing amplified transients. Quantitative comparisons with ${N-}$body simulations yield consistent scalings with the number of particles and with the self-gravity of the disc: the fewer particles and the colder the disc, the faster the secular evolution. Secular evolution is driven by resonances, but does not depend on the initial phases of the disc. For a Mestel disc with ${Q \sim 1.5}$, the polarization cloud around each star boosts up its secular effect by a factor of the order of a thousand or more, promoting accordingly the dynamical relevance of self-induced collisional secular evolution. The position and shape of the induced resonant ridge are found to be in very good agreement with the prediction of the Balescu-Lenard equation, which scales with the square of the susceptibility of the disc. In astrophysics, the inhomogeneous Balescu-Lenard equation may describe the secular diffusion of giant molecular clouds in galactic discs, the secular migration and segregation of planetesimals in proto-planetary discs, or even the long-term evolution of population of stars within the Galactic centre. It could be used as a valuable check of the accuracy of ${N-}$body integrators over secular timescales. △ Less

Submitted 24 July, 2015; originally announced July 2015.

Comments: 21 pages, 19 figures

Journal ref: A&A 584, A129 (2015)

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

Authors: Jean-Baptiste Fouvry, Christophe Pichon, Pierre-Henri Chavanis

Abstract: The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is described using the inhomogeneous Balescu-Lenard equation. Assuming that only tightly wound transient spirals are present in the disc, a WKB approximation provides a simple and tractable quadrature for the corresponding drift and diffusion coefficients. It provides insight into the phys… ▽ More The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is described using the inhomogeneous Balescu-Lenard equation. Assuming that only tightly wound transient spirals are present in the disc, a WKB approximation provides a simple and tractable quadrature for the corresponding drift and diffusion coefficients. It provides insight into the physical processes at work during the secular diffusion of a self-gravitating discrete disc and makes quantitative predictions on the initial variations of the distribution function in action space. When applied to the secular evolution of an isolated stationary self-gravitating Mestel disc, this formalism predicts initially the importance of the corotation resonance in the inner regions of the disc leading to a regime involving radial migration and heating. It predicts in particular 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 likely to resolve this discrepancy. In astrophysics, the inhomogeneous Balescu-Lenard equation and its WKB limit may also describe the secular diffusion of giant molecular clouds in galactic discs, the secular migration and segregation of planetesimals in proto-planetary discs, or even the long-term evolution of population of stars within the Galactic center. △ Less

Submitted 20 April, 2015; originally announced April 2015.

Comments: 22 pages, 12 figures

Journal ref: A&A 581, A139 (2015)

Self-gravity, resonances and orbital diffusion in stellar discs

Authors: Jean-Baptiste Fouvry, James Binney, Christophe Pichon

Abstract: 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 once the diffusion tensor is known. We present the formalism that enables one to compute the diffusion tensor from a given source of noise in the gravitational field when the system's dynamical response to t… ▽ More 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 once the diffusion tensor is known. We present the formalism that enables one to compute the diffusion tensor from a given source of noise in the gravitational field when the system's dynamical response to that noise is included. In the case of a cool stellar disc we are able to reduce the computation of the diffusion tensor to a one-dimensional integral. We implement this formula for a tapered Mestel disc that is exposed to shot noise and find that we are able to explain analytically the principal features of a numerical simulation of such a disc. In particular the formation of narrow ridges of enhanced density in action space is recovered. As the disc's value of Toomre's $Q$ is reduced and the disc becomes more responsive, there is a transition from a regime of heating in the inner regions of the disc through the inner Lindblad resonance to one of radial migration of near-circular orbits via the corotation resonance in the intermediate regions of the disc. The formalism developed here provides the ideal framework in which to study the long-term evolution of all kinds of stellar discs. △ Less

Submitted 19 April, 2015; originally announced April 2015.

Comments: 11 pages, 7 figures

Journal ref: ApJ 806, 117 (2015)

Secular resonant dressed orbital diffusion II : application to an isolated self similar tepid galactic disc

Authors: Jean-Baptiste Fouvry, Christophe Pichon

Abstract: The main orbital signatures of the secular evolution of an isolated self-gravitating stellar Mestel disc are recovered using a dressed Fokker-Planck formalism in angle-action variables. The shot-noise-driven formation of narrow ridges of resonant orbits is recovered in the WKB limit of tightly wound transient spirals, for a tepid Toomre-stable tapered disc. The relative effect of the bulge, the ha… ▽ More The main orbital signatures of the secular evolution of an isolated self-gravitating stellar Mestel disc are recovered using a dressed Fokker-Planck formalism in angle-action variables. The shot-noise-driven formation of narrow ridges of resonant orbits is recovered in the WKB limit of tightly wound transient spirals, for a tepid Toomre-stable tapered disc. The relative effect of the bulge, the halo, the disc temperature and the spectral properties of the shot noise are investigated in turn. For such galactic discs all elements seem to impact the locus and direction of the ridge. For instance, when the halo mass is decreased, we observe a transition between a regime of heating in the inner regions of the disc through the inner Lindblad resonance to a regime of radial migration of quasi-circular orbits via the corotation resonance in the outer part of the disc. The dressed secular formalism captures both the nature of collisionless systems (via their natural frequencies and susceptibility), and their nurture via the structure of the external perturbing power spectrum. Hence it provides the ideal framework in which to study their long term evolution. △ Less

Submitted 19 April, 2015; originally announced April 2015.

Comments: 15 pages, 11 figures

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 449 (2): 1982-1995

Secular resonant dressed orbital diffusion I : method and WKB limit for tepid discs

Authors: Jean-Baptiste Fouvry, Christophe Pichon, Simon Prunet

Abstract: The equation describing the secular diffusion of a self-gravitating collisionless system induced by an exterior perturbation is derived while assuming that the timescale corresponding to secular evolution is much larger than that corresponding to the natural frequencies of the system. Its two dimensional formulation for a tepid galactic disc is also derived using the epicyclic approximation. Its W… ▽ More The equation describing the secular diffusion of a self-gravitating collisionless system induced by an exterior perturbation is derived while assuming that the timescale corresponding to secular evolution is much larger than that corresponding to the natural frequencies of the system. Its two dimensional formulation for a tepid galactic disc is also derived using the epicyclic approximation. Its WKB limit is found while assuming that only tightly wound transient spirals are sustained by the disc. It yields a simple quadrature for the diffusion coefficients which provides a straightforward understanding of the loci of maximal diffusion within the disc. △ Less

Submitted 19 April, 2015; originally announced April 2015.

Comments: 15 pages, 1 figure

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 449 (2): 1967-1981