No Arabic abstract
We describe a method to extract resonant orbits from N-body simulations exploiting the fact that they close in a frame rotating with a constant pattern speed. Our method is applied to the N-body simulation of the Milky Way by Shen et al. (2010). This simulation hosts a massive bar, which drives strong resonances and persistent angular momentum exchange. Resonant orbits are found throughout the disc, both close to the bar itself and out to the very edges of the disc. Using Fourier spectrograms, we demonstrate that the bar is driving kinematic substructure even in the very outer parts of the disc. We identify two major orbit families in the outskirts of the disc that make significant contributions to the kinematic landscape, namely the m:l = 3:-2 and 1:-1 families resonating with the pattern speed of the bar. A mechanism is described that produces bimodal distributions of Galactocentric radial velocities at selected azimuths in the outer disc. It occurs as a result of the temporal coherence of particles on the 3:-2 resonant orbits, which causes them to arrive simultaneously at pericentre or apocentre. This resonant clumping, due to the in-phase motion of the particles through their epicycle, leads to both inward and outward moving groups which belong to the same orbital family and consequently produce bimodal radial velocity distributions. This is a possible explanation of the bimodal velocity distributions observed towards the Galactic anti-Centre by Liu et al. (2012). Another consequence is that transient overdensities appear and dissipate (in a symmetric fashion) on timescales equal to the their epicyclic period resulting in a periodic pulsing of the discs surface density.
We examine an ensemble of 48 simulated clusters to determine the effects of small-scale density fluctuations and large-scale substructure on X-ray measurements of the intracluster medium (ICM) mass. We measure RMS density fluctuations in the ICM which can be characterized by a mean mass-weighted clumping factor C = <rho^2>/<rho>^2 between 1.3 and 1.4 within a density contrast of 500 times the critical density. These fluctuations arise from the cluster history of accretion shocks and major mergers, and their presence enhances the clusters luminosity relative to the smooth case. We expect, therefore, that ICM mass measurements utilizing models which assume uniform density at a given radius carry a bias of order sqrt(C) = 1.16. We verify this result by performing ICM mass measurements on X-ray images of the simulations and finding the expected level of bias. The varied cluster morphologies in our ensemble also allow us to investigate the effects of departures from spherical symmetry on our measurements. We find that the presence of large-scale substructure does not further bias the resulting gas mass unless it is pronounced enough to produce a second peak in the image of at least 1% the maximum surface brightness. We analyze the subset of images with no secondary peaks and find a bias of 9% and a Gaussian random error of 4% in the derived mass.
Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark-matter halo, which has important consequences for dark-matter direct-detection experiments.
We have measured the amount of kinematic substructure in the Galactic halo using the final data set from the Spaghetti project, a pencil-beam high latitude sky survey. Our sample contains 101 photometrically selected and spectroscopically confirmed giants with accurate distance, radial velocity and metallicity information. We have developed a new clustering estimator: the 4distance measure, which when applied to our data set leads to the identification of 1 group and 7 pairs of clumped stars. The group, with 6 members, can confidently be matched to tidal debris of the Sagittarius dwarf galaxy. Two pairs match the properties of known Virgo structures. Using models of the disruption of Sagittarius in Galactic potentials with different degrees of dark halo flattening, we show that this favors a spherical or prolate halo shape, as demonstrated by Newberg et al. (2007) using SDSS data. One additional pair can be linked to older Sagittarius debris. We find that 20% of the stars in the Spaghetti data set are in substructures. From comparison with random data sets we derive a very conservative lower limit of 10% to the amount of substructure in the halo. However, comparison to numerical simulations shows that our results are also consistent with a halo entirely built up from disrupted satellites, provided the dominating features are relatively broad due to early merging or relatively heavy progenitor satellites.
Although thick stellar discs are detected in nearly all edge-on disc galaxies, their formation scenarios still remain a matter of debate. Due to observational difficulties, there is a lack of information about their stellar populations. Using the Russian 6-m telescope BTA we collected deep spectra of thick discs in three edge-on S0-a disc galaxies located in different environments: NGC4111 in a dense group, NGC4710 in the Virgo cluster, and NGC5422 in a sparse group. We see intermediate age (4-5 Gyr) metal rich ([Fe/H] $sim$ -0.2 - 0.0 dex) stellar populations in NGC4111 and NGC4710. On the other hand, NGC5422 does not harbour young stars, its disc is thick and old (10 Gyr), without evidence for a second component, and its $alpha$-element abundance suggests a 1.5-2 Gyr long formation epoch implying its formation at high redshift. Our results suggest the diversity of thick disc formation scenarios.
We present simulations of isolated disc galaxies in a realistic environment performed with the Tree-SPMHD-Code Gadget-3. Our simulations include a spherical circum-galactic medium (CGM) surrounding the galactic disc, motivated by observations and the results of cosmological simulations. We present three galactic models with different halo masses between 10e10 Msol and 10e12 Msol, and for each we use two different approaches to seed the magnetic field, as well as a control simulation without a magnetic field. We find that the amplification of the magnetic field in the centre of the disc leads to a biconical magnetic outflow of gas that magnetizes the CGM. This biconical magnetic outflow reduces the star formation rate (SFR) of the galaxy by roughly 40 percent compared to the simulations without magnetic fields. As the key aspect of our simulations, we find that small scale turbulent motion of the gas in the disc leads to the amplification of the magnetic field up to tens of 10e-6 G, as long as the magnetic field strength is low. For stronger magnetic fields turbulent motion does not lead to significant amplification but is replaced by an alpha-omega dynamo. The occurance of a small scale turbulent dynamo becomes apparent through the magnetic power spectrum and analysis of the field lines curvature. In accordance with recent observations we find an anti-correlation between the spiral structure in the gas density and in the magnetic field due to a diffusion term added to the induction equation.