A numerical study of spiral galaxies and dynamical effects of spiral arms

We use high-resolution numerical simulations of disc galaxies to study spiral arms and their dynamical influence on nearby stars. We find that in the numerical simulations the spiral arms are winding, 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 validate this for normal and barred-spirals, and demonstrate that there is no significant offset of different star-forming tracers across the spiral arm. We show from the dynamics of nearby star particles that star particles are drawn towards and join the spiral arm from behind (in front of) the arm and migrate towards the outer (inner) regions of the disc until the arm disappears. The resulting gain (loss) of angular momentum of star particles behind (in front of) the spiral arm is termed radial migration, and occurs over the entire radial range analyzed. A direct consequence of ubiquitous co-rotation radii is that the guiding centre radii of particles are changed while their degree of random motion is unchanged. We show that there are several types of migrator particles, as well as particles that do not migrate, that each have a different orbital evolution. We show that the orbital type depends on position in configuration and phase space, and show how the orbits can make up moving group features in velocity space. We investigate the correlation between the galactic shear rate and the pitch angle of both individual density waves and the apparent co-rotating spiral arms. We find that, in both cases, higher galactic shear rates produce more tightly wound spiral arms, in agreement with observation. We find also that winding spiral arms can naturally explain the scatter in the pitch angle-shear rate relation seen from observations

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

Recent numerical N-body simulations of spiral galaxies have shown that spiral arms in N-body simulations do not rotate rigidly as expected in classic density wave theory, but instead seem to rotate at a similar speed to the local rotation speed of the stellar disc material. This in turn yields winding, transient and recurrent spiral structure, whose co-rotating nature gives rise to changes in the angular momentum (radial migration) of star particles close to the spiral arm at many radii. From high resolution N-body simulations, we highlight the evolution of strongly migrating star particles (migrators) and star particles that do not migrate (non-migrators) around a spiral arm. We investigate the individual orbit histories of migrators and non-migrators and find that there are several types of migrator and non-migrator, each with unique radial evolution. We find the important quantities that affect the orbital evolution to be the radial and tangential velocity components in combination with the azimuthal distance to the spiral arm at the time the star particle begins to feel tangential force. We contrast each type of orbit to compare how these factors combine for migrators and non-migrators. We find that the positive (negative) migrators sustain a position behind (in front of) the spiral arm, and feel continuous tangential force as long as the spiral arm persists. This is because the positive (negative) migrators are close to the apocentre (pericentre) epicycle phase during their migration, and rotate slower (faster) than the co-rotating spiral arm. On the other hand, non-migrators stay close to the spiral arm, and pass or are passed by the spiral arm one or two times. Although they gain or lose the angular momentum when they are behind or in front of the spiral arm, their net angular momentum change becomes close to zero

Radial migration in numerical simulations of Milky-Way sized galaxies

We show that in rmN -body simulations of isolated spiral discs, spiral arms appear to transient, recurring features that co-rotate with the stellar disc stars at all radii. As a consequence, stars around the spiral arm continually feel a tangential force from the spiral and gain/lose angular momentum at all radii where spiral structure exists, without gaining significant amounts of random energy. We demonstrate that the ubiquitous radial migration in these simulations can be seen as outward (inward) systematic streaming motions along the trailing (leading) side of the spiral arms. We characterise these spiral induced peculiar motions and compare with those of the Milky Way obtained from APOGEE red clump data. We find that transient, co-rotating spiral arms are consistent with the data, in contrast with density wave-like spirals which are qualitatively inconsistent. In addition, we show that, in our simulations, radial migration does not change the radial metallicity gradient significantly, and broadens the metallicity distribution function at all radii, similar to some previous studies

Baryonic effects on the detectability of annihilation radiation from dark matter subhaloes around the Milky Way

We use six, high-resolution Λ-cold dark matter (ΛCDM) simulations of galaxy formation to study how emission from dark matter annihilation is affected by baryonic processes. These simulations produce isolated, disc-dominated galaxies with structure, stellar populations, and stellar and halo masses comparable to those of the Milky Way. They resolve dark matter structures with mass above ∼106 M⊙ and are each available in both full-physics and dark-matter-only versions. In the full-physics case, formation of the stellar galaxy enhances annihilation radiation from the dominant smooth component of the galactic halo by a factor of 3, and its central concentration increases substantially. In contrast, subhalo fluxes are reduced by almost an order of magnitude, partly because of changes in internal structure, partly because of increased tidal effects; they drop relative to the flux from the smooth halo by 1.5 orders of magnitude. The expected flux from the brightest Milky Way subhalo is four orders of magnitude below that from the smooth halo, making it very unlikely that any subhalo will be detected before robust detection of the inner Galaxy. We use recent simulations of halo structure across the full ΛCDM mass range to extrapolate to the smallest (Earth-mass) subhaloes, concluding, in contrast to earlier work, that the total annihilation flux from Milky Way subhaloes will be less than that from the smooth halo, as viewed both from the Sun and by a distant observer. Fermi-Large Area Telescope may marginally resolve annihilation radiation from the very brightest subhaloes, which, typically, will contain stars

Numerical simulations of bubble-induced star formation in dwarf irregular galaxies with a novel stellar feedback scheme

To study the star formation and feedback mechanism, we simulate the evolution of an isolated dwarf irregular galaxy (dIrr) in a fixed dark matter halo, similar in size to WLM. We use the new version of our original N-body/smoothed particle chemodynamics code, GCD+, which adopts improved hydrodynamics, metal diffusion between the gas particles and new modelling of star formation and stellar wind and supernovae (SNe) feedback. Comparing the simulations with and without stellar feedback effects, we demonstrate that the collisions of bubbles produced by strong feedback can induce star formation in a more widely spread area. We also demonstrate that the metallicity in star forming regions is kept low due to the mixing of the metal-rich bubbles and the metal-poor inter-stellar medium. Our simulations also suggest that the bubble-induced star formation leads to many counter-rotating stars. The bubble-induced star formation could be a dominant mechanism to maintain star formation in dIrrs, which is different from larger spiral galaxies where the non-axisymmetric structures, such as spiral arms, are a main driver of star formation

Dark matter annihilation and the Galactic Centre Excess

We compare the surface brightness profile and morphology of the Galactic Centre Excess (GCE) identified in wide-angle γ-ray maps from the Fermi-Large Area Telescope (LAT) to dark matter annihilation predictions derived from high-resolution Λ cold dark matter magnetohydrodynamic simulations of galaxy formation. These simulations produce isolated, disc-dominated galaxies with structure, stellar populations, gas content, and stellar and halo masses comparable to those of the Milky Way. For a specific choice of annihilation cross-section, they agree well with the Fermi-LAT data over the full observed angular range, 1°-15°, whereas their dark-matter-only counterparts, lacking any compression of the inner halo by the gravitational effects of the baryons, fail to predict emission as centrally concentrated as observed. These results provide additional support to the hypothesis that the GCE is produced by annihilating dark matter. If, however, it is produced by a different mechanism, they imply a strong upper limit on annihilation rates, which can be translated into upper limits on the expected γ-ray flux not only from the inner Galaxy, but also from any substructure, with or without stars, in the Galactic halo

Stellar Motion around Spiral Arms: Gaia Mock Data

We compare the stellar motion around a spiral arm created in two different scenarios, transient/co-rotating spiral arms and density-wave-like spiral arms. We generate Gaia mock data from snapshots of the simulations following these two scenarios using our stellar population code, SNAPDRAGONS, which takes into account dust extinction and the expected Gaia errors. We compare the observed rotation velocity around a spiral arm similar in position to the Perseus arm, and find that there is a clear difference in the velocity features around the spiral arm between the co-rotating spiral arm and the density-wave-like spiral arm. Our result demonstrates that the volume and accuracy of the Gaia data are sufficient to clearly distinguish these two scenarios of the spiral arms

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., 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 towards 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. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS

Impacts of a Flaring Star-forming Disc and Stellar Radial Mixing on the Vertical Metallicity Gradient

Using idealised N-body simulations of a Milky Way-sized disc galaxy, we qualitatively study how the metallicity distributions of the thin disc star particles are modified by the formation of the bar and spiral arm structures. The thin disc in our numerical experiments initially has a tight negative radial metallicity gradient and a constant vertical scale-height. We show that the radial mixing of stars drives a positive vertical metallicity gradient in the thin disc. On the other hand, if the initial thin disc is flared, with vertical scale-height increasing with galactocentric radius, the metal poor stars originally in the outer disc become dominant in regions above the disc plane at every radii. This process can drive a negative vertical metallicity gradient, which is consistent with the current observed trend. This model mimics a scenario where the star-forming thin disc was flared in the outer region at earlier epochs. Our numerical experiment with an initial flared disc predicts that the negative vertical metallicity gradient of the mono-age relatively young thin disc population should be steeper in the inner disc, and the radial metallicity gradient of the mono-age population should be shallower at greater heights above the disc plane. We also predict that the metallicity distribution function of mono-age young thin disc populations above the disc plane would be more positively skewed in the inner disc compared to the outer disc

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

Recent numerical N-body simulations of spiral galaxies have shown that spiral arms in N-body simulations do not rotate rigidly as expected in classic density wave theory, but instead seem to rotate at a similar speed to the local rotation speed of the stellar disc material. This in turn yields winding, transient and recurrent spiral structure, whose corotating nature gives rise to changes in the angular momentum (radial migration) of star particles close to the spiral arm at many radii. From high-resolution N-body simulations, we highlight the evolution of strongly migrating star particles (migrators) and star particles that do not migrate (non-migrators) around a spiral arm. We investigate the individual orbit histories of migrators and non-migrators and find that there are several types of migrator and non-migrator, each with unique radial evolution. We find the important quantities that affect the orbital evolution to be the radial and tangential velocity components in combination with the azimuthal distance to the spiral arm at the time the star particle begins to feel tangential force. We contrast each type of orbit to compare how these factors combine for migrators and non-migrators. We find that the positive (negative) migrators sustain a position behind (in front of) the spiral arm, and feel continuous tangential force as long as the spiral arm persists. This is because the positive (negative) migrators are close to the apocentre (pericentre) epicycle phase during their migration, and rotate slower (faster) than the corotating spiral arm. On the other hand, non-migrators stay close to the spiral arm, and pass or are passed by the spiral arm one or two times. Although they gain or lose the angular momentum when they are behind or in front of the spiral arm, their net angular momentum change becomes close to zero. We discuss also the long-term effects of radial migration on the radial metallicity distribution and radial angular momentum and mass profiles. © 2014 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society