We perform the first self-consistent measurement of the rate of interactions between stellar tidal streams created by disrupting satellites and dark subhalos in a cosmological simulation of a Milky-Way-mass galaxy. Using a retagged version of the Aquarius A dark-matter-only simulation, we selected 18 streams of tagged star particles that appear thin at the present day and followed them from the point their progenitors accrete onto the main halo, recording in each snapshot the characteristics of all dark-matter subhalos passing within several distance thresholds of any tagged star particle in each stream. We considered distance thresholds corresponding to constant impact parameters (1, 2, and 5 kpc), as well as those proportional to the region of influence of each subhalo (one and two times its half-mass radius $r_{1/2}$). We then measured the age and present-day, phase-unwrapped length of each stream in order to compute the interaction rate in different mass bins and for different thresholds, and compared these to analytic predictions from the literature. We measure a median rate of $1.5^{+3.0}_{-1.1}\ (9.1^{+17.5}_{-7.1},\ 61.8^{+211}_{-40.6})$ interactions within 1 (2, 5) kpc of the stream per 10 kpc of stream length per 10 Gyr. Resolution effects (both time and particle number) affect these estimated rates by lowering them.

We employ high-resolution N-body simulations of isolated spiral galaxy models, from low-amplitude, multi-armed galaxies to Milky Way-like disks, to estimate the vertical action of ensembles of stars in an axisymmetric potential. In the multi-armed galaxy the low-amplitude arms represent tiny perturbations of the potential, hence the vertical action for a set of stars is conserved, although after several orbital periods of revolution the conservation degrades significantly. For a Milky Way-like galaxy with vigorous spiral activity and the formation of a bar, our results show that the potential is far from steady, implying that the action is not a constant of motion. Furthermore, because of the presence of high-amplitude arms and the bar, considerable in-plane and vertical heating occurs that forces stars to deviate from near-circular orbits, reducing the degree at which the actions are conserved for individual stars, in agreement with previous results, but also for ensembles of stars. If confirmed, this result has several implications, including the assertion that the thick disk of our Galaxy forms by radial migration of stars, under the assumption of the conservation of the action describing the vertical motion of stars.

We use numerical simulations to examine the effects of radial migration on the vertical structure of galaxy disks. The simulations follow three exponential disks of different mass but similar circular velocity, radial scalelength, and (constant) scale height. The disks develop different non-axisymmetric patterns, ranging from feeble, long-lived multiple arms to strong, rapidly-evolving few-armed spirals. These fluctuations induce radial migration through secular changes in the angular momentum of disk particles, mixing the disk radially and blurring pre-existing gradients. Migration affects primarily stars with small vertical excursions, regardless of spiral pattern. This "provenance bias" largely determines the vertical structure of migrating stars: inward migrators thin down as they move in, whereas outward migrators do not thicken up but rather preserve the disk scale height at destination. Migrators of equal birth radius thus develop a strong scale-height gradient, not by flaring out as commonly assumed, but by thinning down as they spread inward. Similar gradients have been observed for low-[$\alpha$/Fe] mono-abundance populations (MAPs) in the Galaxy but our results argue against interpreting them as a consequence of radial migration. This is because outward migration does not lead to thickening, implying that the maximum scale height of any population should reflect its value at birth. In contrast, Galactic MAPs have scale heights that increase monotonically outwards, reaching values that greatly exceed those at their presumed birth radii. Given the strong vertical bias affecting migration, a proper assessment of the importance of radial migration in the Galaxy should take carefully into account the strong radial dependence of the scale heights of the various stellar populations.

We use a set of high-resolution N-body simulations of the Galactic disk to study its interactions with the population of satellites predicted cosmologically. One simulation illustrates that multiple passages of massive satellites with different velocities through the disk generate a wobble, having the appearance of rings in face-on projections of the stellar disk. They also produce flares in the disk outer parts and gradually heat the disk through bending waves. A different numerical experiment shows that an individual satellite as massive as the Sagittarius dwarf galaxy passing through the disk will drive coupled horizontal and vertical oscillations of stars in underdense regions, with small significant associated heating. This experiment shows that vertical excursions of stars in these low-density regions can exceed 1 kpc in the Solar neighborhood, resembling the coherent vertical oscillations recently detected locally. They can also induce non-zero vertical streaming motions as large as 10-20 km s$^{-1}$, consistent with recent observations in the Galactic disk. This phenomenon appears as a local ring, with modest associated disk heating.

We compare the dark matter halos' structural parameters derived for four Milky Way dwarf spheroidal galaxies to those of subhalos found in cosmological $N$-body simulations. We confirm that estimates of the mass at a single fixed radius are fully consistent with the observations. However, when a second structural parameter such as the logarithmic slope of the dark halo density profile measured close to the half-light radius is included in the comparison, we find little to no overlap between the satellites and the subhalos. Typically the right mass subhalos have steeper profiles at these radii than measurements of the dSph suggest. Using energy arguments we explore if it is possible to solve this discrepancy by invoking baryonic effects. Assuming that feedback from supernovae can lead to a reshaping of the halos, we compute the required efficiency and find entirely plausible values for a significant fraction of the subhalos and even as low as 0.1%. This implies that care must be taken not to exaggerate the effect of supernovae feedback as this could make the halos too shallow. These results could be used to calibrate and possibly constrain feedback recipes in hydrodynamical simulations.

We study the radial migration of stars driven by recurring multi-arm spiral features in an exponential disk embedded in a dark matter halo. The spiral perturbations redistribute angular momentum within the disk and lead to substantial radial displacements of individual stars, in a manner that largely preserves the circularity of their orbits and that results, after 5 Gyr (~40 full rotations at the disk scalelength), in little radial heating and no appreciable changes to the vertical or radial structure of the disk. Our results clarify a number of issues related to the spatial distribution and kinematics of migrators. In particular, we find that migrators are a heavily biased subset of stars with preferentially low vertical velocity dispersions. This "provenance bias" for migrators is not surprising in hindsight, for stars with small vertical excursions spend more time near the disk plane and thus respond more readily to non-axisymmetric perturbations. We also find that the vertical velocity dispersion of outward migrators always decreases, whereas the opposite holds for inward migrators. To first order, newly arrived migrators simply replace stars that have migrated off to other radii, thus inheriting the vertical bias of the latter. Extreme migrators might therefore be recognized, if present, by the unexpectedly small amplitude of their vertical excursions. Our results show that migration, understood as changes in angular momentum that preserve circularity, can affect strongly the thin disk, but cast doubts on models that envision the Galactic thick disk as a relic of radial migration.

We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21cm absorption spectra from the 21cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the HI line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes.

We analyze the Aquarius simulations to characterize the shape of dark matter halos with peak circular velocity in the range 8<Vmax<200 km/s, and perform a convergence study using the various Aquarius resolution levels. For the converged objects, we determine the principal axis (a<b<c) of the normalized inertia tensor as a function of radius. We find that the triaxiality of field halos is an increasing function of halo mass, so that the smallest halos in our sample are ~40-50% rounder than Milky Way-like objects at the radius where the circular velocity peaks, rmax. We find that the distribution of subhalo axis ratios is consistent with that of field halos of comparable Vmax. Inner and outer contours within each object are well aligned, with the major axis preferentially pointing in the radial direction for subhalos closest to the center of their host halo. We also analyze the dynamical structure of subhalos likely to host luminous satellites comparable to the classical dwarf spheroidals in the Local Group. These halos have axis ratios that increase with radius, and which are mildly triaxial with <b/a>~0.75 and <c/a>~0.60 at r~1 kpc. Their velocity ellipsoid become strongly tangentially biased in the outskirts as a consequence of tidal stripping.

We propose a new model for the dark matter halo of the Milky Way that fits the properties of the stellar stream associated with the Sagittarius dwarf galaxy. Our dark halo is oblate with q_z = 0.9 for r < 10 kpc, and can be made to follow the Law & Majewski model at larger radii. However, we find that the dynamical perturbations induced by the Large Magellanic Cloud on the orbit of Sgr cannot be neglected when modeling its streams. When taken into account, this leads us to constrain the Galaxy's outer halo shape to have minor-to-major axis ratio (c/a)_\Phi = 0.8 and intermediate-to-major axis ratio (b/a)_\Phi = 0.9, in good agreement with cosmological expectations.

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.

We present a new analysis of the Aquarius simulations done in combination with a semi-analytic galaxy formation model. Our goal is to establish whether the subhalos present in LCDM simulations of Milky Way-like systems could host the dwarf spheroidal (dSph) satellites of our Galaxy. Our analysis shows that, contrary to what has been assumed in most previous work, the mass profiles of subhalos are generally not well fit by NFW models but that Einasto profiles are preferred. We find that for shape parameters alpha = 0.2 - 0.5 and Vmax = 10 - 30 km/s there is very good correspondence with the observational constraints obtained for the nine brightest dSph of the Milky Way. However, to explain the internal dynamics of these systems as well as the number of objects of a given circular velocity the total mass of the Milky Way should be ~ 8x10^11 Msun, a value that is in agreement with many recent determinations, and at the low mass end of the range explored by the Aquarius simulations. Our simulations show important scatter in the number of bright satellites, even when the Aquarius Milky Way-like hosts are scaled to a common mass, and we find no evidence for a missing population of massive subhalos in the Galaxy. This conclusion is also supported when we examine the dynamics of the satellites of M31.

One of the strongest predictions of the LambdaCDM cosmological model is the presence of dark satellites orbiting all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and demonstrate that these encounters can be dramatic. Although mergers with M_sat > M_d are not very common, because of the lower baryonic content they occur much more frequently on the dwarf scale than for L_*-galaxies. As an example, we present a numerical simulation of a 20% (virial) mass ratio merger between a dark satellite and a disky dwarf (akin to the Fornax dwarf galaxy in luminosity) that shows that the merger remnant has a spheroidal morphology. We conclude that perturbations by dark satellites provide a plausible path for the formation of dSph systems and also could trigger starbursts in gas rich dwarf galaxies. Therefore the transition from disky to the often amorphous, irregular, or spheroidal morphologies of dwarfs could be a natural consequence of the dynamical heating of hitherto unobservable dark satellites.

We use the high resolution cosmological N-body simulations from the Aquarius project to investigate in detail the mechanisms that determine the shape of Milky Way-type dark matter haloes. We find that, when measured at the instantaneous virial radius, the shape of individual haloes changes with time, evolving from a typically prolate configuration at early stages to a more triaxial/oblate geometry at the present day. This evolution in halo shape correlates well with the distribution of the infalling material: prolate configurations arise when haloes are fed through narrow filaments, which characterizes the early epochs of halo assembly, whereas triaxial/oblate configurations result as the accretion turns more isotropic at later times. Interestingly, at redshift z=0, clear imprints of the past history of each halo are recorded in their shapes at different radii, which also exhibit a variation from prolate in the inner regions to triaxial/oblate in the outskirts. Provided that the Aquarius haloes are fair representatives of Milky Way-like 10^12 Msun objects, we conclude that the shape of such dark matter haloes is a complex, time-dependent property, with each radial shell retaining memory of the conditions at the time of collapse.