Theoretical astrophysics and plasma physics at RPC

Abstract:

We use the second data releases of the European Space AgencyGaia astrometric survey and the high-resolution Galactic Archaeology with HERMES (GALAH) spectroscopic survey to analyse the structure of our Galaxy’s disc components. With GALAH, we separate the α-rich and α-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zVz phase plane, for both Vϕ and VR, and recover the remarkable ‘phase spiral’ discovered by Gaia. We identify the anticipated quadrupole signature in zVz of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the α-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zVz phase plane lasts ≳2Gyr⁠, but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited ≲0.5Gyr ago by an object like Sgr with total mass ∼3 × 1010 M⊙ (stripped down from ∼5 × 1010 M⊙ when it first entered the halo) passing through the plane.

Abstract:

We present a method for fitting orbit-superposition models to the kinematics of discrete stellar systems when the available stellar sample has been filtered by a known selection function. The fitting method can be applied to any model in which the distribution function is represented as a linear superposition of basis elements with unknown weights. As an example, we apply it to Fritz et al.'s kinematics of the innermost regions of the Milky Way's nuclear stellar cluster. Assuming spherical symmetry, our models fit a black hole of mass $M_\bullet=(3.76\pm0.22)\times10^6\,M_\odot$, surrounded by an extended mass $M_\star=(6.57\pm0.54)\times10^6\,M_\odot$ within $4\,\pc$. Within $1\,\pc$ the best-fitting mass models have an approximate power-law density cusp $\rho\propto r^{-\gamma}$ with $\gamma=1.3\pm0.3$. We carry out an extensive investigation of how our modelling assumptions might bias these estimates: $M_\bullet$ is the most robust parameter and $\gamma$ the least. Internally the best-fitting models have broadly isotropic orbit distributions, apart from a bias towards circular orbits between 0.1 and 0.3 parsec.