The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\phi)$ 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 $\Sigma_{\rm max} (R_{\odot})= 14\pm5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 20\%$ at the solar radius, with an increase towards larger Galactocentric radii. Our model constrains the pitch angle of the {\it dynamical} spiral arms to be $12.7^{\circ}\pm1.2^{\circ}$.

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.

Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $\lambda \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.

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.

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.

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.

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.

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.

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 $\eta$-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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.