The data obtained by the recent modern sky surveys enable detailed studies of the stellar distribution in the multi-dimensional space spanned by spatial coordinates, velocity and metallicity, from ...the solar neighborhood all the way out to the outer Milky Way halo. While these results represent exciting observational breakthroughs, their interpretation is not simple. For example, traditional decomposition of the thin and thick disks predicts a strong correlation in metallicity and kinematics at \(\sim\)1 kpc from the Galactic plane; however, recent SDSS--based work has demonstrated an absence of this correlation for disk stars. Instead, the variation of the metallicity and rotational velocity distributions can be modeled using non--Gaussian functions that retain their shapes and only shift as the distance from the mid--plane increases. To fully contextualize these recent observational results, a detailed comparison with sophisticated numerical models is necessary. Modern simulations have sufficient resolution and physical detail to study the formation of stellar disks and spheroids over a large baseline of masses and cosmic ages. We discuss preliminary comparisons of various observed maps and N--body model predictions and find them encouraging. In particular, the N--body disk models of Roškar et al. \cite{Roskar 2008} reproduce a change of disk scale height reminiscent of thin/thick disk decomposition, as well as metallicity and rotational velocity gradients, while not inducing a correlation of the latter two quantities, in qualitative agreement with SDSS observations.
We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are ...taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydrogen dissociation stage. The effects of clump rotation, asymmetries, and radiative cooling are studied. Rotation provides support against fast collapse, but non-axisymmetric modes develop and efficiently transport angular momentum outward, forming a circumplanetary disc. This transport helps the clump reach the dynamical collapse phase, resulting from molecular hydrogen dissociation, on a thousand-year timescale, which is smaller than timescales predicted by some previous spherical 1D collapse models. Extrapolation to the threshold of the runaway hydrogen dissociation indicates that the collapse timescales can be shorter than inward migration timescales, suggesting that clumps could survive tidal disruption and deliver a proto-gas giant to distances of even a few AU from the central star.
The center of our disk galaxy, the Milky Way, is dominated by a boxy/peanut-shaped bulge. Numerous studies of the bulge based on stellar photometry have concluded that the bulge stars are exclusively ...old. The perceived lack of young stars in the bulge strongly constrains its likely formation scenarios, providing evidence that the bulge is a unique population that formed early and separately from the disk. However, recent studies of individual bulge stars using the microlensing technique have reported that they span a range of ages, emphasizing that the bulge may not be a monolithic structure. In this letter we demonstrate that the presence of young stars that are located predominantly near the plane is expected for a bulge that has formed from the disk via dynamical instabilities. Using an N-body+SPH simulation of a disk galaxy forming out of gas cooling inside a dark matter halo and forming stars, we find a qualitative agreement between our model and the observations of young metal-rich stars in the bulge. We are also able to partially resolve the apparent contradiction in the literature between results that argue for a purely old bulge population and those which show a population comprised of a range in ages; the key is where to look.
Abstract
We seek to understand the origin of radial migration in spiral galaxies by analysing in detail the structure and evolution of an idealized, isolated galactic disc. To understand the ...redistribution of stars, we characterize the time evolution of properties of spirals that spontaneously form in the disc. Our models unambiguously show that in such discs, single spirals are unlikely, but that a number of transient patterns may coexist in the disc. However, we also show that while spirals are transient in amplitude, at any given time the disc favours patterns of certain pattern speeds. Using several runs with different numerical parameters we show that the properties of spirals that occur spontaneously in the disc do not sensitively depend on resolution. The existence of multiple transient patterns has large implications for the orbits of stars in the disc, and we therefore examine the resonant scattering mechanisms that profoundly alter angular momenta of individual stars. We confirm that the corotation scattering mechanism described by Sellwood & Binney is responsible for the largest angular momentum changes in our simulations.
We present evidence that isolated growing discs, subject to internal spiral perturbations, thicken due to both heating and radial migration. We show this by demonstrating that the thickness and ...vertical velocity dispersions of coeval stars depend on their age as well as the change in their radii. While the disc thickens due to internal processes, we find that this induces only a minor amount of flaring. We further demonstrate the consequences of such thickening on the structural properties of stellar populations and find that they qualitatively agree with recent studies of the Milky Way disc.
The role of gas in the mass assembly at the nuclei of galaxies is still subject to some uncertainty. Stellar nuclear discs bridge the gap between the large-scale galaxy and the central massive ...objects that reside there. Using a high-resolution simulation of a galaxy forming out of gas cooling and settling into a disc, we study the formation and properties of nuclear discs. Gas, driven to the centre by a bar, settles into a rotating star-forming nuclear disc (ND). This ND is thinner, younger, kinematically cooler and more metal rich than the surrounding bar. The ND is elliptical and orthogonal to the bar. The complex kinematics in the region of the ND are a result of the superposition of older stars streaming along the bar and younger stars circulating within the ND. The signature of the ND is therefore subtle in the kinematics. Instead the ND stands out clearly in metallicity and age maps. We compare the model to the density and kinematics of real galaxies with NDs finding qualitative similarities. Our results suggest that gas dissipation is very important for forming nuclear structures.
Using high resolution, fully cosmological smoothed particle hydro-dynamical simulations of dwarf galaxies in a Lambda cold dark matter Universe, we show how baryons attain a final angular momentum ...distribution which allows pure disc galaxies to form. Blowing out substantial amounts of gas through supernovae and stellar winds, which is well supported observationally, is a key ingredient in forming bulgeless discs. We outline why galactic outflows preferentially remove low angular momentum material, and show that this is a natural result when structure forms in a cold dark matter cosmology. The driving factors are a) the mean angular momentum of accreted material increases with time, b) lower potentials at early times, c) the existence of an extended reservoir of high angular momentum gas which is not within star forming regions, meaning that only gas from the inner region (low angular momentum gas) is expelled and d) the tendency for outflows to follow the path of least resistance which is perpendicular to the disc. We also show that outflows are enhanced during mergers, thus expelling much of the gas which has lost its angular momentum during these events, and preventing the formation of "classical", merger driven bulges in low mass systems. Stars formed prior to such mergers form a diffuse, extended stellar halo component.
We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and ...kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding ~10 kpc and galactocentric distances exceeding ~20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix, we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of ~20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of ~20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio q super(Phi) sub(DM) = 0.7 + or - 0.1; this corresponds to an axis ratio q super(rho) sub(DM) ~ 0.4 + or - 0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.
Accurate phase-space coordinates (three components of position and velocity) of individual stars are rapidly becoming available with current and future resolved star surveys. These data will enable ...the computation of the full three-dimensional orbits of tens of thousands of stars in the Milky Way's stellar halo. We demonstrate that the analysis of stellar halo orbits in frequency space can be used to construct a 'frequency map' which provides a highly compact, yet intuitively informative way to represent the six-dimensional halo phase-space distribution function. This representation readily reveals the most important major orbit families in the halo, and the relative abundances of the different orbit families, which in turn reflect the shape and orientation of the dark matter halo relative to the disc. We demonstrate the value of frequency space orbit analysis by applying the method to halo orbits in a series of controlled simulations of disc galaxies. We show that the disc influences the shape of the inner halo making it nearly oblate, but the outer halo remains largely unaffected. Since the shape of the halo varies with radius, the frequency map provides a more versatile way to identify major and minor orbit families than traditional orbit classification schemes. Although the shape of the halo varies with radius, frequency maps of local samples of halo orbits confined to the inner halo contain most of the information about the global shape of the halo and its major orbit families. Frequency maps show that adiabatic growth of a disc traps halo orbits in numerous resonant orbit families (i.e. having commensurable frequencies). The locations and strengths of these resonant families are determined by both the global shape of the halo and its stellar distribution function. If a good estimate of the Galactic potential in the inner halo (within ∼ 50 kpc) is available, the appearance of strong, stable resonances in frequency maps of halo orbits will allow us to determine the degree of resonant trapping induced by the disc potential. We show that if the Galactic potential is not known exactly, a measure of the diffusion rate of a large sample of ∼ 104 halo orbits can help distinguish between the true potential and an incorrect potential. The orbital spectral analysis methods described in this paper provide a strong complementarity to existing methods for constraining the potential of the Milky Way halo and its stellar distribution function.
We present a detailed analysis of a disc galaxy forming in a high-resolution fully cosmological simulation to investigate the nature of the outer regions of discs and their relevance for the disc ...formation process. Specifically, we focus on the phenomenon of misaligned disc components and find that the outer disc warp is a consequence of the misalignment between the inner disc and the surrounding hot gaseous halo. As the infalling cold gas sinks towards the centre of the galaxy, it is strongly torqued by the hot gas halo. By the time the fresh gas reaches the central disc-forming region, its angular momentum is completely aligned with the spin of the hot gas halo. If the spin of the hot gas halo, in turn, is not aligned with that of the inner disc, a misaligned outer disc forms, comprised of newly accreted material. The inner and outer components are misaligned with each other because they respond differently to infalling substructure and accretion. The warped disc feeds the main gas disc due to viscous angular momentum losses, but small amounts of star formation in the warp itself form a low-metallicity thick disc. We show that observations of resolved stellar populations in warped galaxies in the local Universe could provide evidence for the presence of these processes and therefore indirectly reveal ongoing gas accretion and the existence of hot gas haloes.