Distribution functions for resonantly trapped orbits in the Galactic disc Giacomo Monari, Benoit Famaey, Jean-Baptiste Fouvry, James Binney
To cite this version:
Giacomo Monari, Benoit Famaey, Jean-Baptiste Fouvry, James Binney. Distribution functions for resonantly trapped orbits in the Galactic disc. Monthly Notices of the Royal Astronomical Society, Oxford University Press (OUP): Policy P - Oxford Open Option A, 2017, 471 (4), pp.4314-4322. 10.1093/mnras/stx1825. hal-03146696
HAL Id: hal-03146696 https://hal.archives-ouvertes.fr/hal-03146696 Submitted on 1 Mar 2021
HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. MNRAS 471, 4314–4322 (2017) doi:10.1093/mnras/stx1825 Advance Access publication 2017 July 20
Distribution functions for resonantly trapped orbits in the Galactic disc
Giacomo Monari,1‹ Benoit Famaey,2 Jean-Baptiste Fouvry3† and James Binney4 1The Oskar Klein Centre for Cosmoparticle Physics, Department of Physics, Stockholm University, AlbaNova, SE-10691 Stockholm, Sweden 2Observatoire astronomique de Strasbourg, Universite´ de Strasbourg, CNRS UMR 7550, 11 rue de l’Universite,´ F-67000 Strasbourg, France 3Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540, USA 4Rudolf Peierls Centre for Theoretical Physics, Keble Road, Oxford OX1 3NP, UK
Accepted 2017 July 17. Received 2017 July 17; in original form 2017 June 12 Downloaded from https://academic.oup.com/mnras/article/471/4/4314/3980207 by guest on 01 March 2021
ABSTRACT The present-day response of a Galactic disc stellar population to a non-axisymmetric pertur- bation 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 compar- isons 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, 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 depends only 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 DFs in action space. Key words: Galaxy: disc – Galaxy: kinematics and dynamics – Galaxy: structure.
that the phase-space position is 2π-periodic in them (e.g. Binney & 1 INTRODUCTION Tremaine 2008). In the absence of perturbations, these angles evolve The optimal exploitation of the next data releases of the Gaia linearly with time, θ(t) = θ 0 + t,where ( J) ≡ ∂H0/∂ J is the mission (Gaia Collaboration et al. 2016) will necessarily in- vector of fundamental orbital frequencies. In an equilibrium con- volve the construction of a fully dynamical model of the Milky figuration, the angle coordinates of stars are phase-mixed on or- Way. Rather than trying to construct a quixotic full ab initio bital tori that are defined by the actions J alone, and the phase- 3 hydrodynamical model of the Galaxy, which would never be space density of stars f0( J)d J corresponds to the number of able to reproduce all the details of the Gaia data, a promising stars dN in a given infinitesimal action range divided by (2π)3. approach is to rather construct a multicomponent phase-space In an axisymmetric configuration, the action variables can simply distribution function (DF) representing each stellar population be chosen to be the radial, azimuthal and vertical actions, respec- as well as dark matter, and to compute the potential that tively. By constructing DFs depending on these action variables, these populations jointly generate (e.g. Binney & Piffl 2015). the current best axisymmetric models of the Milky Way have been To do so, one can make use of the Jeans theorem, constraining the constructed (Cole & Binney 2017). DF of an equilibrium configuration to depend only on three integrals The Milky Way is, however, not axisymmetric: It harbours both of motion. Choosing three integrals of motion that have canonically a bar (de Vaucouleurs 1964; Binney et al. 1991; Binney, Ger- conjugate variables allows us to express the Hamiltonian H0 in its hard & Spergel 1997; Wegg, Gerhard & Portail 2015; Monari simplest form, i.e. depending only on these three integrals. Such et al. 2017a,b) and spiral arms, the exact number, dynamics and integrals are called the ‘action variables’ J, and are new gener- nature of which are still under debate (Sellwood & Carlberg 2014; alized momenta having the dimension of velocity times distance, Grand et al. 2015). Whilst Trick et al. (2017) showed that spiral while their dimensionless canonically conjugate variables are called arms might not affect the axisymmetric fit, the combined effects of the ‘angle variables’ θ, because they are usually normalized such spiral arms and the central bar of the Milky Way are clearly impor- tant observationally (e.g. McMillan 2013;Bovyetal.2015). Hence, non-axisymmetric DFs are needed to pin down the present structure E-mail: [email protected] of the non-axisymmetric components of the potential, which have † Hubble Fellow. enormous importance as drivers of the secular evolution of the disc
C 2017 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society Distribution functions for trapped orbits 4315 Downloaded from https://academic.oup.com/mnras/article/471/4/4314/3980207 by guest on 01 March 2021
Figure 1. Two orbits in the Galactic potential described in this work, shown in the reference frame corotating with the bar. Here, the pattern speed is chosen to be b = 1.8 0. The red points correspond to the position of the relative apocentres. The thick horizontal line represents the long axis of the bar. The bar rotates counterclockwise. Left-hand panel: an orbit trapped to the outer Lindblad resonance, where the angle φmax (dashed lines) represents the maximum extension of the apocentres excursions. Right-hand panel: a circulating orbit.
(Fouvry, Binney & Pichon 2015; Fouvry, Pichon & Prunet 2015; average over the fast angles (Arnold 1978). Under this averaged Aumer, Binney & Schonrich¨ 2016; Aumer & Binney 2017). Hamiltonian, the slow variables have the dynamics of a pendulum. A recent step (Monari, Famaey & Siebert 2016) has been to In Section 4, we introduce the pendulum’s angles and actions. In derive from perturbation theory explicit DFs for present-day snap- Section 5, we discuss how to build the DF using the newly intro- shots of the disc as a function of the actions and angles of the duced pendulum angles and actions, both inside and outside the unperturbed axisymmetric system. This work, which is an Eulerian zone of trapping at resonances. In Section 6, we present the form of approach to the problem posed by non-axisymmetry, has allowed us the DFs in velocity space, in cases of astrophysical interest related to probe the effect of stationary spiral arms in three spatial dimen- to the Galactic bar. We conclude in Section 7. sions, away from the main resonances. In particular, the moments of the perturbed DF describe ‘breathing’ modes of the Galactic disc in perfect accordance with simulations (Monari et al. 2016). 2 THE BAR AND TRAPPED ORBITS Such a breathing mode might actually have been detected in the extended solar neighbourhood (Williams et al. 2013), but with a Let us consider orbits in the Galactic plane, and let (R, φ)bethe larger amplitude, perhaps because the spiral arms are transient. Galactocentric radius and azimuth. The logarithmic potential, cor- = Although such an Eulerian treatment has also been used to gain responding to a flat circular velocity curve vc(R) v0, is a rough but qualitative insights into the effects of non-axisymmetries near res- simple representation of the potential of the Galaxy. In a formula, onances (Monari et al. 2017a), no quantitative assessments can (R) = v2 ln(R/R ), (1) be made with such an approach because the linear treatment di- 0 0 verges at resonances (the problem of small divisors; Binney & with R0 the distance of the Sun from the Galactic centre. Motion Tremaine 2008). in this planar axisymmetric potential admits actions Jφ ,whichis In the present contribution, we solve this problem by developing simply the angular momentum about the Galactic Centre, and JR, the Lagrangian approach to the impact of non-axisymmetries at res- which quantifies the extent of radial oscillations. onances. The basic idea is to model the deformation of the orbital tori Let the axisymmetric potential be perturbed by a non- outside of the trapping region, and to construct new tori, complete axisymmetric component with a new system of angle–action variables, within the trapping − = im(φ bt) region (Kaasalainen 1994). Finally, following Binney (2016), we 1(R,φ,t) Re a(R)e , (2) populate the new tori by phase-averaging the unperturbed DF over where is the pattern speed. We will specialize to the bar adopted the new tori. b by Dehnen (2000) and Monari et al. (2017a), so we set m = 2and In Section 2, we present some examples of trapped and untrapped adopt orbits. In Section 3, we summarize the standard approach to a reso- nance, namely to make a canonical transformation to fast and slow R5 2 (R/R )−3 R ≥ R , (R) =−α 0 0 × b b (3) angles and actions, and to replace the real Hamiltonian with its a b 3 − 3 3Rb 2 (R/Rb) R MNRAS 471, 4314–4322 (2017) 4316 G. Monari et al. Downloaded from https://academic.oup.com/mnras/article/471/4/4314/3980207 by guest on 01 March 2021 Figure 2. The slow variables θs + π (left-hand panel) and Js (right-hand panel) for the librating orbit of Fig. 1. The blue lines correspond to the numerically integrated orbit and the orange lines to the pendulum approximation. where ≡ (R ) is the circular frequency at the solar radius, 2 = 0 0 Here the guiding radius Rg( Jφ )isdefinedbyRg (Rg) Jφ , the ec- R is the length of the bar and α represents the maximum ratio b b centricity by e(J ,J ) = 2J /(κR2)andγ = 2 (R )/κ(R ).2 between the radial force contributed by the bar and the axisymmetric R φ R g g g background at the Sun (see also Monari, Famaey & Siebert 2015; In the epicyclic approximation, the radial and azimuthal fre- −1 quencies of an orbit of actions (JR, Jφ )are R = κ(Rg)and Monari et al. 2016). We further set R0 = 8.2 kpc, v0 = 243 km s , φ = (Rg) + [dκ(Rg)/dJφ ]JR, respectively. Rb = 3.5 kpc and α = 0.01. At a resonance, the orbital frequencies R and φ satisfy Fig. 1 shows for b = 1.8 0 two orbits in the reference frame corotating with the bar. The red points in this plot correspond to the l R + m( φ − b) = 0. (6) relative apocentres of the orbit.1 In the left-hand panel of Fig. 1, the azimuths of the apocentre points do not cover the whole [0, 2π] The three main resonances are the corotation resonance (l = 0), = =− range, but rather oscillate within the interval [−φmax, φmax]. The and the outer (l 1) and inner (l 1) Lindblad resonances. To orbit is said to be ‘trapped’ at a resonance, in this case the outer capture the behaviour of the slow- and fast-varying motions near the Lindblad resonance (see below), and the angle of the apocentres, resonances, one makes a canonical transformation of coordinates the ‘precession angle’, is said to ‘librate’. For comparison, the right- defined by the following time-dependent generating function of hand panel of Fig. 1 shows an orbit for which the precession angle type 2 (e.g. Weinberg 1994): coversthewhole[0, 2π] range. This orbit is said to ‘circulate’. In S(θ ,θ ,J ,J ,t) = lθ + m θ − t J + θ J . (7) the following, we give a quantitative description of these two types R φ s f R φ b s R f of orbits, using perturbation theory. The new angles and actions (θ f, θ s, Jf, Js) are then related to the old ones by 3 REDUCTION TO A PENDULUM θ = lθ + m(θ − t),J = mJ , s R φ b φ s (8) θ = θ ,J= lJ + J . We will hereafter work within the epicyclic approximation f R R s f (Binney & Tremaine 2008), in which radial oscillations are har- By taking the time-derivative of θ s in the unperturbed system, and monic with angular frequency κ(R), so the radial action JR = ER/κ, by the definition of the resonance in equation (6), one finds that the where ER = E − Ec is the energy of these oscillations (Ec being evolution of the new ‘slow’ angle θ s is indeed slow near a resonance. the energy of a circular orbit with the same angular momentum Jφ ). Given that along nearly circular orbits θ φ φ, from equation (8), we Using the formulae reported in Dehnen (1999) and Monari et al. can understand why θ s represents the azimuth of the apocentres and (2017a), we can then relate the coordinates of the trapped orbit to pericentres3 of the orbit in the frame of reference that corotates with the angles θ R and θ φ of the unperturbed system. We can also rewrite the bar: At θ R = 0(θR =−π), we are at the pericentre (apocentre) the perturbing potential in the Galactic plane in terms of actions and of the orbit and θ s (θs + π)ism times the star’s azimuth θ φ − bt angles as a Fourier series in the frame of reference corotating with the orbit. ⎧ ⎫ In the new canonical coordinates of equation (8), the motion in ⎨ 1 ⎬ + − the perturbed system is described by the following Hamiltonian (J ,J ,θ ,θ ) = Re c ei[jθR m(θφ bt)] , (4) 1 R φ R φ ⎩ jm ⎭ (also called the Jacobi integral): j=−1 ⎧ ⎫ ⎨ 1 ⎬ with − + = + i[(j l)θf θs] − H H0 Re ⎩ cjme ⎭ m bJs, (9) m j=−1 c (J ,J ) ≡ δ + δ| | sgn(j)γe (R , 0) jm R φ j0 j 1 2 a g where H0 is the Hamiltonian of the unperturbed axisymmetric Rg ∂a − δ| | e (R , 0). (5) system and the cjm coefficients are the Fourier coefficients from j 1 2 ∂R g equation (5), expressed as functions of the actions (Jf, Js), thanks 1 While in an axisymmetric potential the apocentres of the orbits in the √ √ Galactic plane are always at the same distance R from the Galactic Centre, 2 In the potential from equation (1), κ = 2 ,soγ = 2. this is not the case in a non-axisymmetric potential like the one used in this 3 According to the convention of Dehnen (1999), in this work the angle work. Therefore, the plotted points correspond to relative apocentres. θ R = 0 at the pericentre, and θR = π at the apocentre. MNRAS 471, 4314–4322 (2017) Distribution functions for trapped orbits 4317 to the canonical transformation from equation (8). Since θ f evolves much faster than θ s,weaverageH along θ f (the averaging princi- ple, e.g. Arnold 1978; Weinberg 1994; Binney & Tremaine 2008), to obtain iθs H = H0(Jf,Js) − m bJs + Re clm(Jf ,Js)e . (10) Since J˙f =−∂H/∂θf = 0, Jf is an integral of motion, and for each Jf, the motion of every orbit can be described purely in the (θ s, Js) plane. For each value of Jf, let us then define Js,res as the value of Js where