The Andromeda galaxy (M31) contains a box/peanut bulge (BPB) entangled with a classical bulge (CB) requiring a triaxial modelling to determine the dynamics, stellar and dark matter mass. We construct made-to-measure models fitting new VIRUS-W IFU bulge stellar kinematic observations, the IRAC-3.6$\mu$m photometry, and the disc's HI rotation curve. We explore the parameter space for the 3.6$\mu$m mass-to-light ratio $(\Upsilon_{3.6})$, the bar pattern speed ($\Omega_p$), and the dark matter mass in the composite bulge ($M^B_{DM}$) within 3.2kpc. Considering Einasto dark matter profiles, we find the best models for $\Upsilon_{3.6}=0.72\pm0.02\,M_\odot/L_\odot$, $M^B_{DM}=1.2^{+0.2}_{-0.4}\times10^{10}M_\odot$ and $\Omega_p=40\pm5\,km/s/kpc$. These models have a dynamical bulge mass of $M_{dyn}^B=4.25^{+0.10}_{-0.29}\times10^{10}M_{\odot}$ including a stellar mass of $M^B=3.09^{+0.10}_{-0.12}\times10^{10}M_\odot$(73%), of which the CB has $M^{CB}=1.18^{+0.06}_{-0.07}\times10^{10}M_\odot$(28%) and the BPB $M^{BPB}=1.91\pm0.06\times10^{10}M_\odot$(45%). We also explore models with NFW haloes finding that, while the Einasto models better fit the stellar kinematics, the obtained parameters agree within the errors. The $M^B_{DM}$ values agree with adiabatically contracted cosmological NFW haloes with M31's virial mass and radius. The best model has two bulge components with completely different kinematics that only together successfully reproduce the observations ($\mu_{3.6},\upsilon_{los},\sigma_{los},h3,h4$). The modelling includes dust absorption which reproduces the observed kinematic asymmetries. Our results provide new constraints for the early formation of M31 given the lower mass found for the classical bulge and the shallow dark matter profile, as well as the secular evolution of M31 implied by the bar and its resonant interactions with the classical bulge, stellar halo and disc.
The Hercules stream is a group of comoving stars in the solar neighbourhood, which can potentially be explained as a signature of either the outer Lindblad resonance (OLR) of a fast Galactic bar or the corotation resonance (CR) of a slower bar. In either case, the feature should be present over a large area of the disc. With the recent commissioning of the APOGEE-2 Southern spectrograph we can search for the Hercules stream at (l, b) = (270°, 0), a direction in which the Hercules stream, if caused by the bar's OLR, would be strong enough to be detected using only the line-of-sight velocities. We clearly detect a narrow, Hercules-like feature in the data that can be traced from the solar neighbourhood to a distance of about 4 kpc. The detected feature matches well the line-of-sight velocity distribution from the fast-bar (OLR) model. Confronting the data with a model where the Hercules stream is caused by the CR of a slower bar leads to a poorer match, as the corotation model does not predict clearly separated modes, possibly because the slow-bar model is too hot.
Andromeda is our nearest neighbouring disk galaxy and a prime target for detailed modelling of the evolutionary processes that shape galaxies. We analyse the nature of M31's triaxial bulge with an extensive set of N-body models, which include Box/Peanut (B/P) bulges as well as initial classical bulges (ICBs). Comparing with IRAC 3.6$\mu m$ data, only one model matches simultaneously all the morphological properties of M31's bulge, and requires an ICB and a B/P bulge with 1/3 and 2/3 of the total bulge mass respectively. We find that our pure B/P bulge models do not show concentrations high enough to match the S\'ersic index ($n$) and the effective radius of M31's bulge. Instead, the best model requires an ICB component with mass $M^{\rm ICB}=1.1\times10^{10}{\rm M_{\odot}}$ and three-dimensional half-mass radius $r_{\rm half}^{\rm ICB}$=0.53 kpc (140 arcsec). The B/P bulge component has a mass of $M^{\rm B/P}=2.2\times10^{10}{\rm M_{\odot}}$ and a half-mass radius of $r_{\rm half}^{\rm B/P}$=1.3 kpc (340 arcsec). The model's B/P bulge extends to $r^{\rm B/P}$=3.2 kpc (840 arcsec) in the plane of the disk, as does M31's bulge. In this composite bulge model, the ICB component explains the velocity dispersion drop observed in the centre within $R<$190 pc (50 arcsec), while the B/P bulge component reproduces the observed rapid rotation and the kinematic twist of the observed zero velocity line. This model's pattern speed is $\Omega_p$=38 km/s/kpc, placing corotation at $r_{\rm cor}$=5.8 kpc (1500 arcsec). The outer Lindblad resonance (OLR) is then at $r_{\rm OLR}$=10.4kpc, near the 10kpc-ring of M31, suggesting that this structure may be related to the bar's OLR. By comparison with an earlier snapshot, we estimate that M31's thin bar extends to $r_{\rm bar}^{\rm thin}\sim$4.0 kpc (1000 arcsec) in the disk plane, and in projection extends to $R_{\rm bar}^{\rm thin}\sim$2.3 kpc (600 arcsec).
We investigated the relationship between the kinematics and mass of young (<3 × 108 yr) white dwarfs using proper motions. Our sample is taken from the colour-selected catalogues of the Sloan Digital Sky Survey and the Palomar–Green Survey, both of which have spectroscopic temperature and gravity determinations. We find that the dispersion decreases with increasing white dwarf mass. This can be explained as a result of less scattering by objects in the Galactic disc during the shorter lifetime of their more massive progenitors. A direct result of this is that white dwarfs with high mass have a reduced scale height, and hence their local density is enhanced over their less massive counterparts. In addition, we have investigated whether the kinematics of the highest mass white dwarfs (>0.95 M⊙) are consistent with the expected relative contributions of single star evolution and mergers. We find that the kinematics are consistent with the majority of high-mass white dwarfs being formed through single star evolution.
We consider the formation of extreme mass-ratio inspirals (EMRIs) sourced from a stellar cusp centred on a primary supermassive black hole (SMBH) and perturbed by an inspiraling less massive secondary SMBH. The problem is approached numerically, assuming the stars are non-interacting over these short time-scales and performing an ensemble of restricted three-body integrations. From these simulations, we see that not only can EMRIs be produced during this process, but the dynamics are also quite rich. In particular, most of the EMRIs are produced through a process akin to the Kozai–Lidov mechanism, but with strong effects due to the non-Keplerian stellar potential, general relativity and non-secular oscillations in the angular momentum on the orbital time-scale of the binary SMBH system.
We construct a large set of dynamical models of the galactic bulge, bar and inner disc using the made-to-measure method. Our models are constrained to match the red clump giant density from a combination of the VVV, UKIDSS and 2MASS infrared surveys together with stellar kinematics in the bulge from the BRAVA and OGLE surveys, and in the entire bar region from the ARGOS Survey. We are able to recover the bar pattern speed and the stellar and dark matter mass distributions in the bar region, thus recovering the entire galactic effective potential. We find a bar pattern speed of 39.0 ± 3.5 km s− 1 kpc− 1, placing the bar corotation radius at 6.1 ± 0.5 kpc and making the Milky Way bar a typical fast rotator. We evaluate the stellar mass of the long bar and bulge structure to be Mbar/bulge = 1.88 ± 0.12 × 1010 M⊙, larger than the mass of disc in the bar region, Minner disc = 1.29 ± 0.12 × 1010 M⊙. The total dynamical mass in the bulge volume is 1.85 ± 0.05 × 1010 M⊙. Thanks to more extended kinematic data sets and recent measurement of the bulge initial mass function, our models have a low dark matter fraction in the bulge of 17 ± 2 per cent. We find a dark matter density profile which flattens to a shallow cusp or core in the bulge region. Finally, we find dynamical evidence for an extra central mass of ∼ 0.2 × 1010 M⊙, probably in a nuclear disc or discy pseudo-bulge.
We present the first self-consistent chemodynamical model fitted to reproduce data for the galactic bulge, bar and inner disc. We extend the Made-to-Measure method to an augmented phase-space including the metallicity of stars, and show its first application to the bar region of the Milky Way. Using data from the ARGOS and APOGEE (DR12) surveys, we adapt the recent dynamical model from Portail et al. to reproduce the observed spatial and kinematic variations as a function of metallicity, thus allowing the detailed study of the 3D density distributions, kinematics and orbital structure of stars in different metallicity bins. We find that metal-rich stars with [Fe/H] ≥ −0.5 are strongly barred and have dynamical properties that are consistent with a common disc origin. Metal-poor stars with [Fe/H] ≤ −0.5 show strong kinematic variations with metallicity, indicating varying contributions from the underlying stellar populations. Outside the central kpc, metal-poor stars are found to have the density and kinematics of a thick disc while in the inner kpc, evidence for an extra concentration of metal-poor stars is found. Finally, the combined orbit distributions of all metallicities in the model naturally reproduce the observed vertex deviations in the bulge. This paper demonstrates the power of Made-to-Measure chemodynamical models, that when extended to other chemical dimensions will be very powerful tools to maximize the information obtained from large spectroscopic surveys such as APOGEE, GALAH and MOONS.
We derive new constraints on the mass, rotation, orbit structure and statistical parallax of the Galactic old nuclear star cluster (NSC) and the mass of the supermassive black hole. We combine star counts and kinematic data from Fritz et al (2014), including 2'500 line-of-sight velocities and 10'000 proper motions. We show that the difference between the proper motion dispersions sigma_l and sigma_b cannot be explained by rotation, but is a consequence of the flattening of the NSC. We fit the surface density distribution of stars in the central 1000" by a spheroidal cluster with scale ~100" and a much larger nuclear disk component. We compute the two-integral distribution function f(E,Lz) for this density model, and add rotation self-consistently. We find that: (i) The orbit structure of the f(E,Lz) gives an excellent match to the observed velocity dispersion profiles as well as the proper motion and line-of-sight velocity histograms, including the double-peak in the v_l-histograms. (ii) This requires an axial ratio of q= 0.73+-0.04 for r<70" consistent with our determination from star counts. (iii) The NSC is approximately described by an isotropic rotator model. (iv) Using the corresponding Jeans equations to fit the proper motion and line-of-sight velocity dispersions, we obtain best estimates for the NSC mass, black hole mass, and distance M*(r<100")=(8.94+-0.31|stat+-0.9|syst)x10^6Msun, Mbh=(3.86+-0.14|stat+-0.4|syst)x10^6Msun, and R0=8.27+-0.09|stat+-0.1|syst kpc, where the systematic errors estimate additional uncertainties in the dynamical modeling. (v) The combination of the cluster dynamics with the S-star orbits around Sgr A* strongly reduces the degeneracy between black hole mass and Galactic centre distance present in previous S-star studies. A joint statistical analysis with the results of Gillessen et al (2009) gives Mbh=(4.23+-0.14)x10^6Msun and R0=8.33+-0.11kpc.
ABSTRACT We study the late-time evolution of the central regions of two Milky Way (MW)-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, Rb, measurements fluctuate on a dynamical time-scale by up to 100 per cent, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15 per cent, correlating with Rb. The Tremaine–Weinberg method estimates of the bars’ instantaneous pattern speeds show variations around the mean of up to $\sim \!20{{\ \rm per\ cent}}$, typically anticorrelating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ∼60–200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50 per cent of bar measurements in MW-mass external galaxies, the bar lengths of SBab-type galaxies are overestimated by $\sim \!15{{\ \rm per\ cent}}$ and those of SBbc types by $\sim \!55{{\ \rm per\ cent}}$. Consequently, bars longer than their corotation radius reported in the literature, dubbed ‘ultrafast bars’, may simply correspond to the largest biases. Given that the Scutum–Centaurus arm is likely connected to the near half of the MW bar, recent direct measurements may be overestimating its length by 1–1.5 kpc, while its present pattern speed may be 5–10 $\rm km\ s^{-1}\ kpc^{-1}$ smaller than its time-averaged value.
Abstract We have derived absolute proper motions of the entire Galactic bulge region from VVV Infrared Astrometric Catalogue (VIRAC) and Gaia. We present these both as integrated on-sky maps and, after isolating standard candle red clump (RC) stars, as a function of distance using RC magnitude as a proxy. These data provide a new global, 3D view of the Milky Way barred bulge kinematics. We find a gradient in the mean longitudinal proper motion, $\langle \mu _ l^\star \rangle $, between the different sides of the bar, which is sensitive to the bar pattern speed. The split RC has distinct proper motions and is colder than other stars at similar distance. The proper motion correlation map has a quadrupole pattern in all magnitude slices showing no evidence for a separate, more axisymmetric inner bulge component. The line-of-sight integrated kinematic maps show a high central velocity dispersion surrounded by a more asymmetric dispersion profile. $\sigma _{\mu _l} / \sigma _{\mu _b}$ is smallest, ≈1.1, near the minor axis and reaches ≈1.4 near the disc plane. The integrated $\langle \mu_b\rangle$ pattern signals a superposition of bar rotation and internal streaming motion, with the near part shrinking in latitude and the far part expanding. To understand and interpret these remarkable data, we compare to a made-to-measure barred dynamical model, folding in the VIRAC selection function to construct mock maps. We find that our model of the barred bulge, with a pattern speed of 37.5 $\mathrm{ \mathrm{ km \, s^{-1}} \, kpc^{-1} }$, is able to reproduce all observed features impressively well. Dynamical models like this will be key to unlocking the full potential of these data.