Stellar migration, whether due to trapping by transient spirals (churning), or to scattering by non-axisymmetric perturbations, has been proposed to explain the presence of stars in outer disks. After a review of the basic theory, we present compelling, but not yet conclusive, evidence that churning has been important in the outer disks of galaxies with type II (down-bending) profiles, while scattering has produced the outer disks of type III (up-bending) galaxies. In contrast, field galaxies with type I (pure exponential) profiles appear to not have experienced substantial migration. We conclude by suggesting work that would improve our understanding of the origin of outer disks.

Recent analysis of the SDSS-III/APOGEE Data Release 12 stellar catalogue has revealed that the Milky Way's metallicity distribution function (MDF) changes shape as a function of radius, transitioning from being negatively skewed at small Galactocentric radii to positively skewed at large Galactocentric radii. Using a high resolution, N-body+SPH simulation, we show that the changing skewness arises from radial migration - metal-rich stars form in the inner disk and subsequently migrate to the metal-poorer outer disk. These migrated stars represent a large fraction (> 50%) of the stars in the outer disk; they populate the high metallicity tail of the MDFs and are, in general, more metal-rich than the surrounding outer disk gas. The simulation also reproduces another surprising APOGEE result: the spatially invariant high-[alpha/Fe] MDFs. This arises in the simulation from the migration of a population formed within a narrow range of radii (3.2+/-1.2 kpc) and time (8.8+/-0.6 Gyr ago), rather than from spatially extended star formation in a homogeneous medium at early times. These results point toward the crucial role radial migration has played in shaping our Milky Way.

Large surveys have shown that red galaxies are preferentially aligned with their halos while blue galaxies have a more isotropic distribution. Since halos generally align with their filaments this introduces a bias in the measurement of the cosmic shear from weak lensing. It is therefore vitally important to understand why this difference arises. We explore the stability of different disc orientations within triaxial halos. We show that, in the absence of gas, the disc orientation is most stable when its spin is along the minor axis of the halo. Instead when gas cools onto a disc it is able to form in almost arbitrary orientation, including off the main planes of the halo (but avoiding an orientation perpendicular to the halo's intermediate axis). Substructure helps gasless galaxies reach alignment with the halo faster, but have less effect on galaxies when gas is cooling onto the disc. Our results provide a novel and natural interpretation for why red, gas poor galaxies are preferentially aligned with their halo, while blue, star-forming, galaxies have nearly random orientations, without requiring a connection between galaxies' current star formation rate and their merger history.

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 (2013), 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 qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-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.

We present novel 3D multi-scale SPH simulations of gas-rich galaxy mergers between the most massive galaxies at $z \sim 8 - 10$, designed to scrutinize the direct collapse formation scenario for massive black hole seeds proposed in \citet{mayer+10}. The simulations achieve a resolution of 0.1 pc, and include both metallicity-dependent optically-thin cooling and a model for thermal balance at high optical depth. We consider different formulations of the SPH hydrodynamical equations, including thermal and metal diffusion. When the two merging galaxy cores collide, gas infall produces a compact, optically thick nuclear disk with densities exceeding $10^{-10}$ g cm$^3$. The disk rapidly accretes higher angular momentum gas from its surroundings reaching $\sim 5$ pc and a mass of $\gtrsim 10^9$ $M_{\odot}$ in only a few $10^4$ yr. Outside $\gtrsim 2$ pc it fragments into massive clumps. Instead, supersonic turbulence prevents fragmentation in the inner parsec region, which remains warm ($\sim 3000-6000$ K) and develops strong non-axisymmetric modes that cause prominent radial gas inflows ($> 10^4$ $M_{\odot}$ yr$^{-1}$), forming an ultra-dense massive disky core. Angular momentum transport by non-axisymmetric modes should continue below our spatial resolution limit, quickly turning the disky core into a supermassive protostar which can collapse directly into a massive black hole of mass $10^8-10^9$ $M_{\odot}$ via the relativistic radial instability. Such a "cold direct collapse"' explains naturally the early emergence of high-z QSOs. Its telltale signature would be a burst of gravitational waves in the frequency range $10^{-4} - 10^{-1}$ Hz, possibly detectable by the planned eLISA interferometer.

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.

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.

We study the dynamics of massive black hole pairs in clumpy gaseous circumnuclear disks. We track the orbital decay of the light, secondary black hole $M_{\bullet2}$ orbiting around the more massive primary at the center of the disk, using $N$-body/smoothed particle hydrodynamic simulations. We find that the gravitational interaction of $M_{\bullet2}$ with massive clumps $M_{\rm cl}$ erratically perturbs the otherwise smooth orbital decay. In close encounters with massive clumps, gravitational slingshots can kick the secondary black hole out of the disk plane. The black hole moving on an inclined orbit then experiences the weaker dynamical friction of the stellar background, resulting in a longer orbital decay timescale. Interactions between clumps can also favor orbital decay when the black hole is captured by a massive clump which is segregating toward the center of the disk. The stochastic behavior of the black hole orbit emerges mainly when the ratio $M_{\bullet2}/M_{\rm cl}$ falls below unity, with decay timescales ranging from $\sim1$ to $\sim50$ Myr. This suggests that describing the cold clumpy phase of the inter-stellar medium in self-consistent simulations of galaxy mergers, albeit so far neglected, is important to predict the black hole dynamics in galaxy merger remnants.

Models of the Sagittarius Stream have consistently found that the Milky Way disc is oriented such that its short axis is along the intermediate axis of the triaxial dark matter halo. We attempt to build models of disc galaxies in such an `intermediate-axis orientation'. We do this with three models. In the first two cases we simply rigidly grow a disc in a triaxial halo such that the disc ends up perpendicular to the global intermediate axis. We also attempt to coax a disc to form in an intermediate-axis orientation by producing a gas+dark matter triaxial system with gas angular momentum about the intermediate axis. In all cases we fail to produce systems which remain with stellar angular momentum aligned with the halo's intermediate axis, even when the disc's potential flattens the inner halo such that the disc is everywhere perpendicular to the halo's local minor axis. For one of these unstable simulations we show that the potential is even rounder than the models of the Milky Way potential in the region probed by the Sagittarius Stream. We conclude that the Milky Way's disc is very unlikely to be in an intermediate axis orientation. However we find that a disc can persist off one of the principal planes of the potential. We propose that the disc of the Milky Way must be tilted relative to the principal axes of the dark matter halo. Direct confirmation of this prediction would constitute a critical test of Modified Newtonian Dynamics.

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.

We analyzed the radial surface brightness profile of the spiral galaxy NGC 7793 using HST/ACS images from the GHOSTS survey and a new HST/WFC3 image across the disk break. We used the photometry of resolved stars to select distinct populations covering a wide range of stellar ages. We found breaks in the radial profiles of all stellar populations at 280" (~5.1 kpc). Beyond this disk break, the profiles become steeper for younger populations. This same trend is seen in numerical simulations where the outer disk is formed almost entirely by radial migration. We also found that the older stars of NGC 7793 extend significantly farther than the underlying HI disk. They are thus unlikely to have formed entirely at their current radii, unless the gas disk was substantially larger in the past. These observations thus provide evidence for substantial stellar radial migration in late-type disks.

We use the Mitchell Spectrograph (formerly VIRUS-P) to observe 12 nearby disk galaxies. We successfully measure ages in the outer disk in six systems. In three cases (NGC 2684, NGC 6155, and NGC 7437), we find that a downward break in the disk surface brightness profile corresponds with a change in the dominant stellar population with the interior being dominated by active star formation and the exterior having older stellar populations that are best-fit with star formation histories that decline with time. The observed increase in average stellar ages beyond a profile break is similar to theoretical models that predict surface brightness breaks are caused by stellar migration, with the outer disk being populated from scattered old interior stars. In three more cases (IC 1132, NGC 4904, and NGC 6691), we find no significant change in the stellar population as one crosses the break radius. In these galaxies, both the inner and outer disks are dominated by active star formation and younger stellar populations. While radial migration can contribute to the stellar populations beyond the break, it appears more than one mechanism is required to explain all of our observed stellar profile breaks.

The phase space coordinates of individual halo stars obtained by Galactic surveys enable the computation of their full 3-dimensional orbits. Spectral analysis of halo orbits can be used to construct "frequency maps" which provide a compact representation of the 6-dimensional phase space distribution function. Frequency maps identify important major orbit families, and the orbital abundances reflect the shape and orientation of the dark matter halo relative to the disk. We apply spectral analysis to halo orbits in a series of controlled simulations of disk galaxies. Although the shape of the simulated 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. Quiescent or adiabatic disk formation results in significant trapping of halo orbits in resonant orbit families (i.e. orbits with commensurable frequencies). If a good estimate of the Galactic potential in the inner halo (within ~50kpc) 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 disk potential. The locations and strengths of these resonant families are determined both by the global shape of the halo and its distribution function. Identification of such resonances in the Milky Way's stellar halo would therefore provide evidence of an extended period of adiabatic disk growth. If the Galactic potential is not known exactly, a measure of the diffusion rate of large sample of 10^4 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 (ABRIDGED).

We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.

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.

The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular orbits, as described by Sellwood & Binney (2002). As stars move outwards they populate a thickened component. Such stars found at the present time in the solar neighborhood formed early in the disk's history at smaller radii where stars are more metal-poor and alpha-enhanced, leading to the properties observed for thick disk stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to an apparent separation in the other property as observed. This scenario is supported by the SDSS observation that stars in the transition region do not show any correlation between rotational velocity and metallicity. The good qualitative agreement between our simulation and observations in the Milky Way hint that the thick disk may be a ubiquitous galaxy feature generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies.

We present new Hubble Space Telescope (HST) observations of NGC 300 taken as part of the ACS Nearby Galaxy Survey Treasury (ANGST). Individual stars are resolved in these images down to an absolute magnitude of M_F814W = 1.0 (below the red clump). We determine the star formation history of the galaxy in 6 radial bins by comparing our observed color-magnitude diagrams (CMDs) with synthetic CMDs based on theoretical isochrones. We find that the stellar disk out to 5.4 kpc is primarily old, in contrast with the outwardly similar galaxy M33. We determine the scale length as a function of age and find evidence for inside-out growth of the stellar disk: the scale length has increased from 1.1 +/- 0.1 kpc 10 Gyr ago to 1.3 +/- 0.1 kpc at present, indicating a buildup in the fraction of young stars at larger radii. As the scale length of M33 has recently been shown to have increased much more dramatically with time, our results demonstrate that two galaxies with similar sizes and morphologies can have very different histories. With an N-body simulation of a galaxy designed to be similar to NGC 300, we determine that the effects of radial migration should be minimal. We trace the metallicity gradient as a function of time and find a present day metallicity gradient consistent with that seen in previous studies. Consistent results are obtained from archival images covering the same radial extent but differing in placement and filter combination.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present resolved stellar photometry of NGC 2976 obtained with the Advanced Camera for Surveys (ACS) as part of the ACS Nearby Galaxy Survey Treasury (ANGST) program. The data cover the radial extent of the major axis of the disk out to 6 kpc, or ~6 scale lengths. The outer disk was imaged to a depth of M_F606W ~ 1, and an inner field was imaged to the crowding limit at a depth of M_F606W ~ -1. Through detailed analysis and modeling of these CMDs we have reconstructed the star formation history of the stellar populations currently residing in these portions of the galaxy, finding similar ancient populations at all radii but significantly different young populations at increasing radii. In particular, outside of the well-measured break in the disk surface brightness profile, the age of the youngest population increases with distance from the galaxy center, suggesting that star formation is shutting down from the outside-in. We use our measured star formation history, along with H I surface density measurements, to reconstruct the surface density profile of the disk during previous epochs. Comparisons between the recovered star formation rates and reconstructed gas densities at previous epochs are consistent with star formation following the Schmidt law during the past 0.5 Gyrs, but with a drop in star formation efficiency at low gas densities, as seen in local galaxies at the present day. The current rate and gas density suggest that rapid star formation in NGC 2976 is currently in the process of ceasing from the outside-in due to gas depletion. This process of outer disk gas depletion and inner disk star formation was likely triggered by an interaction with the core of the M81 group >~1 Gyr ago that stripped the gas from the galaxy halo and/or triggered gas inflow from the outer disk toward the galaxy center.

We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk.

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\v{s}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.