Robert J. J. Grand

The dynamics of stars around spiral arms

Spiral density wave theory attempts to describe the spiral pattern in spiral galaxies in terms of a long-lived wave structure with a constant pattern speed in order to avoid the winding dilemma. The pattern is consequently a rigidly rotating, long-lived feature. We run N-body simulations of a giant disc galaxy consisting of a pure stellar disc and a static dark matter halo, and find that the spiral arms are transient features whose pattern speeds decrease with radius in such a way that the pattern speed is almost equal to the rotation curve of the galaxy. We trace particle motion around the spiral arms. We show that particles from behind and in front of the spiral arm are drawn towards and join the arm. Particles move along the arm in the radial direction and we find a clear trend that they migrate towards the outer (inner) radii on the trailing (leading) side of the arm. Our simulations demonstrate that because the spiral arm feature is corotating, the particles continue to be accelerated (decelerated) by the spiral arm for long periods, which leads to strong and efficient migration, at all radii in the disc.

Dynamics of stars around spiral arms in an N-body/SPH simulated barred-spiral galaxy

We run N-body smoothed particle hydrodynamics (SPH) simulations of a Milky Way sized galaxy. The code takes into account hydrodynamics, self-gravity, star formation, supernova and stellar wind feedback, radiative cooling and metal enrichment. The simulated galaxy is a barred-spiral galaxy consisting of a stellar and gas disc, enveloped in a static dark matter halo. Similar to what is found in our pure N-body simulation of a non-barred galaxy in Grand et. al. (2012), we find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is similar to the rotation of star particles. Compared to the non-barred case, we find that the spiral arm pattern speed is slightly faster than the rotation speed of star particles: the bar appears to boost the pattern speed ahead of the rotational velocity. We trace particle motion around the spiral arms at different radii, and demonstrate that there are star particles that are drawn towards and join the arm from behind (in front of) the arm and migrate toward the outer (inner) regions of the disc until the arm disappears as a result of their transient nature. We see this migration over the entire radial range analysed, which is a consequence of the spiral arm rotating at similar speeds to star particles at all radii, which is inconsistent with the prediction of classical density wave theory. The bar does not prevent this systematic radial migration, which is shown to largely preserve circular orbits. We also demonstrate that there is no significant offset of different star forming tracers across the spiral arm, which is also inconsistent with the prediction of classical density wave theory.

Dynamics of stars around spiral arms in an n-body/sph simulated barred spiral galaxy

We run N-body smoothed particle hydrodynamics (SPH) simulations of a Milky Way sized galaxy. The code takes into account hydrodynamics, self-gravity, star formation, supernova and stellar wind feedback, radiative cooling and metal enrichment. The simulated galaxy is a barred-spiral galaxy consisting of a stellar and gas disc, enveloped in a static dark matter halo. Similar to what is found in our pure N-body simulation of a non-barred galaxy in Grand et. al. (2012), we find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is similar to the rotation of star particles. Compared to the non-barred case, we find that the spiral arm pattern speed is slightly faster than the rotation speed of star particles: the bar appears to boost the pattern speed ahead of the rotational velocity. We trace particle motion around the spiral arms at different radii, and demonstrate that there are star particles that are drawn towards and join the arm from behind (in front of) the arm and migrate toward the outer (inner) regions of the disc until the arm disappears as a result of their transient nature. We see this migration over the entire radial range analysed, which is a consequence of the spiral arm rotating at similar speeds to star particles at all radii, which is inconsistent with the prediction of classical density wave theory. The bar does not prevent this systematic radial migration, which is shown to largely preserve circular orbits. We also demonstrate that there is no significant offset of different star forming tracers across the spiral arm, which is also inconsistent with the prediction of classical density wave theory.

Vertical disc heating in Milky Way-sized galaxies in a cosmological context

Vertically extended, high velocity dispersion stellar distributions appear to be a ubiquitous feature of disc galaxies, and both internal and external mechanisms have been proposed to be the major driver of their formation. However, it is unclear to what extent each mechanism can generate such a distribution, which is likely to depend on the assembly history of the galaxy. To this end, we perform 16 high-resolution cosmological-zoom simulations of Milky Way-sized galaxies using the state-of-the-art cosmological magnetohydrodynamical code AREPO, and analyse the evolution of the vertical kinematics of the stellar disc in connection with various heating mechanisms. We find that the bar is the dominant heating mechanism in most cases, whereas spiral arms, radial migration and adiabatic heating from mid-plane density growth are all subdominant. The strongest source, though less prevalent than bars, originates from external perturbations from satellites/subhaloes of masses log10(M/M⊙) ≳ 10. However, in many simulations the orbits of newborn star particles become cooler with time, such that they dominate the shape of the age–velocity dispersion relation and overall vertical disc structure unless a strong external perturbation takes place.

The Auriga Project: the properties and formation mechanisms of disc galaxies across cosmic time

We introduce a suite of 30 cosmological magneto-hydrodynamical zoom simulations of the formation of galaxies in isolated Milky Way mass dark haloes. These were carried out with the moving mesh code arepo, together with a comprehensive model for galaxy formation physics, including active galactic nuclei (AGN) feedback and magnetic fields, which produces realistic galaxy populations in large cosmological simulations. We demonstrate that our simulations reproduce a wide range of present-day observables, in particular, two-component disc-dominated galaxies with appropriate stellar masses, sizes, rotation curves, star formation rates and metallicities. We investigate the driving mechanisms that set present-day disc sizes/scalelengths, and find that they are related to the angular momentum of halo material. We show that the largest discs are produced by quiescent mergers that inspiral into the galaxy and deposit high-angular momentum material into the pre-existing disc, simultaneously increasing the spin of dark matter and gas in the halo. More violent mergers and strong AGN feedback play roles in limiting disc size by destroying pre-existing discs and by suppressing gas accretion on to the outer disc, respectively. The most important factor that leads to compact discs, however, is simply a low angular momentum for the halo. In these cases, AGN feedback plays an important role in limiting central star formation and the formation of a massive bulge.

Impact of radial migration on stellar and gas radial metallicity distribution

Radial migration is defined as the change in guiding centre radius of stars and gas caused by gains or losses of angular momentum that result from gravitational interaction with non-axisymmetric structure. This has been shown to have significant impact on the metallicity distribution in galactic discs, and therefore affects the interpretation of Galactic archeology. We use a simulation of a Milky Way-sized galaxy to examine the effect of radial migration on the star and gas radial metallicity distribution. We find that both the star and gas component show significant radial migration. The stellar radial metallicity gradient remains almost unchanged but the radial metallicity distribution of the stars is broadened to produce a greater dispersion at all radii. However, the metallicity dispersion of the gas remains narrow. We find that the main drivers of the gas metallicity distribution evolution are metal enrichment and mixing: more efficient metal enrichment in the inner region maintains a negative slope in the radial metallicity distribution, and the metal mixing ensures the tight relationship of the gas metallicity with the radius. The metallicity distribution function reproduces the trend in the age-metallicity relation found from observations for stars younger than 1.0 Gyr in the Milky Way.

Spiral arm pitch angle and galactic shear rate in N-body simulations of disc galaxies

Spiral galaxies are observed to exhibit a range of morphologies, in particular in the shape of spiral arms. A key diagnostic parameter is the pitch angle, which describes how tightly wound the spiral arms are. Observationally and analytically, a correlation between pitch angle and galactic shear rate has been detected. For the first time, we examine whether this effect is detected in N-body simulations by calculating and comparing pitch angles of both individual density waves and overall spiral structure in a suite of N-body simulations. We find that higher galactic shear rates produce more tightly wound spiral arms, both in individual mode patterns (density waves) and in the overall density enhancement. Although the mode pattern pitch angles by construction remain constant with time, the overall logarithmic spiral arm winds over time, which could help to explain the scatter in the relation between pitch angle versus shear seen from observations. The correlation between spiral arm pitch angle and galactic shear rate that we find in N-body simulations may also explain why late Hubble type of spiral galaxies tend to have more open arms.

A fully cosmological model of a Monoceros-like ring

We study the vertical structure of a stellar disk obtained from a fully cosmological high-resolution hydrodynamical simulation of the formation of a Milky Way-like galaxy. At the present day, the disks mean vertical height shows a well-defined and strong pattern, with amplitudes as large as 3 kpc in its outer regions. This pattern is the result of a satellite - host halo - disk interaction and reproduces, qualitatively, many of the observable properties of the Monoceros Ring. In particular we find disk material at the distance of Monoceros (

Gas and stellar motions and observational signatures of corotating spiral arms

R \sim

Orbits of radial migrators and non-migrators around a spiral arm in N-body simulations

12-16 kpc, galactocentric) extending far above the mid plane ( 30