Using $N$-body+SPH simulations of galaxies falling into a cluster, we study the evolution of their radial density profiles. When evolved in isolation, galaxies develop a type~II (down-bending) profile. In the cluster, the evolution of the profile depends on the minimum cluster-centric radius the galaxy reaches, which controls the degree of ram pressure stripping. If the galaxy falls to $\sim 50\%$ of the virial radius, then the profile remains type~II, but if the galaxy reaches down to $\sim 20\%$ of the virial radius the break weakens and the profile becomes more type~I like. The velocity dispersions are only slightly increased in the cluster simulations compared with the isolated galaxy; random motion therefore cannot be responsible for redistributing material sufficiently to cause the change in the profile type. Instead we find that the joint action of radial migration driven by tidally-induced spirals and the outside-in quenching of star formation due to ram pressure stripping alters the density profile. As a result, this model predicts a flattening of the age profiles amongst cluster lenticulars with type~I profiles, which can be observationally tested. %

We present the study of a set of N-body+SPH simulations of a Milky Way-like system produced by the radiative cooling of hot gas embedded in a dark matter halo. The galaxy and its gaseous halo evolve for 10 Gyr in isolation, which allows us to study how internal processes affect the evolution of the system. We show how the morphology, the kinematics and the evolution of the galaxy are affected by the input supernova feedback energy E$_{\rm SN}$, and we compare its properties with those of the Milky Way. Different values of E$_{\rm SN}$ do not significantly affect the star formation history of the system, but the disc of cold gas gets thicker and more turbulent as feedback increases. Our main result is that, for the highest value of E$_{\rm SN}$ considered, the galaxy shows a prominent layer of extra-planar cold (log(T)<4.3) gas extended up to a few kpc above the disc at column densities of $10^{19}$ cm$^{-2}$. The kinematics of this material is in agreement with that inferred for the HI halos of our Galaxy and NGC 891, although its mass is lower. Also, the location, the kinematics and the typical column densities of the hot (5.3<log(T)<5.7) gas are in good agreement with those determined from the O$_{\rm VI}$ absorption systems in the halo of the Milky Way and external galaxies. In contrast with the observations, however, gas at log(T)<5.3 is lacking in the circumgalactic region of our systems.

The radial density profiles of stellar galaxy discs can be well approximated as an exponential. Compared to this canonical form, however, the profiles in the majority of disc galaxies show downward or upward breaks at large radii. Currently, there is no coherent explanation in a galaxy formation context of the radial profile per se, along with the two types of profile breaks. Using a set of controlled hydrodynamic simulations of disc galaxy formation, we find a correlation between the host halo's initial angular momentum and the resulting radial profile of the stellar disc: galaxies that live in haloes with a low spin parameter {\lambda} <~ 0.03 show an up-bending break in their disc density profiles, while galaxies in haloes of higher angular momentum show a down-bending break. We find that the case of pure exponential profiles ({\lambda} ~ 0.035) coincides with the peak of the spin parameter distribution from cosmological simulations. Our simulations not only imply an explanation of the observed behaviours, but also suggest that the physical origin of this effect is related to the amount of radial redistribution of stellar mass, which is anti-correlated with {\lambda}.

Galaxy formation models and simulations rely on various feedback mechanisms to reproduce the observed baryonic scaling relations and galaxy morphologies. Although dwarf galaxy and giant elliptical properties can be explained using feedback from supernova and active galactic nuclei, Milky Way-sized galaxies still represent a challenge to current theories of galaxy formation. In this paper, we explore the possible role of feedback from stellar radiation in regulating the main properties of disk galaxies such as our own Milky Way. We have performed a suite of cosmological simulations of the same $\sim10^{12} {\rm M}_{\odot}$ halo selected based on its rather typical mass accretion history. We have implemented radiative feedback from young stars using a crude model of radiative transfer for ultraviolet (UV) and infrared (IR) radiation. However, the model is realistic enough such that the dust opacity plays a direct role in regulating the efficiency of our feedback mechanism. We have explored various models for the dust opacity, assuming different constant dust temperatures, as well as a varying dust temperature model. We find that while strong radiative feedback appears as a viable mechanism to regulate the stellar mass fraction in massive galaxies, it also prevents the formation of discs with reasonable morphologies. In models with strong stellar radiation feedback, stellar discs are systematically too thick while the gas disc morphology is completely destroyed due to the efficient mixing between the feedback-affected gas and its surroundings. At the resolution of our simulation suite, we find it impossible to preserve spiral disc morphology while at the same time expelling enough baryons to satisfy the abundance matching constraints.

The redistribution of stars in galactic disks is an important aspect of disk galaxy evolution. Stars that efficiently migrate in such a way that does not also appreciably heat their orbits can drastically affect the stellar populations observed today and therefore influence constraints derived from such observations. Unfortunately, while the theoretical understanding of the migration process is becoming increasingly robust, there are currently few specific observable predictions. As a result, we do not yet have a clear handle on whether the process has been important for the Milky Way in the past or how to constrain it. I discuss some of the expected qualitative outcomes of migration as well as some current controversies.

We simulate an equal-mass merger of two Milky Way-size galaxy discs with moderate gas fractions at parsec-scale resolution including a new model for radiative cooling and heating in a multi-phase medium, as well as star formation and feedback from supernovae. The two discs initially have a $2.6\times10^6\mathrm{~M_{\odot}}$ supermassive black hole (SMBH) embedded in their centers. As the merger completes and the two galactic cores merge, the SMBHs form a a pair with a separation of a few hundred pc that gradually decays. Due to the stochastic nature of the system immediately following the merger, the orbital plane of the binary is significantly perturbed. Furthermore, owing to the strong starburst the gas from the central region is completely evacuated, requiring $\sim10$~Myr for a nuclear disc to rebuild. Most importantly, the clumpy nature of the interstellar medium has a major impact on the the dynamical evolution of the SMBH pair, which undergo gravitational encounters with massive gas clouds and stochastic torquing by both clouds and spiral modes in the disk. These effects combine to greatly delay the decay of the two SMBHs to separations of a few parsecs by nearly two orders of magnitude, $\sim 10^8$ yr, compared to previous work. In mergers of more gas-rich, clumpier galaxies at high redshift stochastic torques will be even more pronounced and potentially lead to stronger modulation of the orbital decay. This suggests that SMBH pairs at separations of several tens of parsecs should be relatively common at any redshift.

We present near-infrared integral field spectroscopy of the central kiloparsec of 17 nearby luminous and ultra-luminous infrared galaxies undergoing major mergers. These observations were taken with OSIRIS assisted by the Keck I and II Adaptive Optics systems, providing spatial resolutions of a few tens of parsecs. The resulting kinematic maps reveal gas disks in at least 16 out of 19 nuclei and stellar disks in 11 out of 11 nuclei observed in these galaxy merger systems. In our late-stages mergers, these disks are young (stellar ages $<30$ Myr) and likely formed as gas disks which became unstable to star formation during the merger. On average, these disks have effective radii of a few hundred parsecs, masses between $10^{8}$ and $10^{10} M_{Sun}$, and $v/\sigma$ between 1 and 5. These disks are similar to those created in high-resolution hydrodynamical simulations of gas-rich galaxy mergers, and favor short coalescence times for binary black holes. The few galaxies in our sample in earlier stages of mergers have disks which are larger ($r_{eff}\sim200-1800$ pc) and likely are remnants of the galactic disks that have not yet been completely disrupted by the merger.

We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [\alpha/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances.

We present evidence that isolated growing discs, subject to internal spiral perturbations, thicken due to both heating \emph{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.

We seek to understand the origin of radial migration in spiral galaxies by analyzing in detail the structure and evolution of an idealized, isolated galactic disk. To understand the redistribution of stars, we characterize the time-evolution of properties of spirals that spontaneously form in the disk. Our models unambiguously show that in such disks, single spirals are unlikely, but that a number of transient patterns may coexist in the disk. However, we also show that while spirals are transient in amplitude, at any given time the disk favors patterns of certain pattern speeds. Using several runs with different numerical parameters we show that the properties of spirals that occur spontaneously in the disk do not sensitively depend on resolution. The existence of multiple transient patterns has large implications for the orbits of stars in the disk, 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 (2002) is responsible for the largest angular momentum changes in our simulations.

Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- components which can be assigned to a thin stellar disc, thick disk, and a low mass stellar halo via a chemical decomposition. The thin and thick disc populations so selected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young (<6.6 Gyr), possess low velocity dispersion ({\sigma}U,V,W = 41, 31, 25 km/s), high [Fe/H], and low [O/Fe]. The thick disc stars are old (6.6<age<9.8 Gyrs), lag the thin disc by \sim21 km/s, possess higher velocity dispersion ({\sigma}U,V,W = 49, 44, 35 km/s), relatively low [Fe/H] and high [O/Fe]. The halo component comprises less than 4% of stars in the "solar annulus" of the simulation, has low metallicity, a velocity ellipsoid defined by ({\sigma}U,V,W = 62, 46, 45 km/s) and is formed primarily in-situ during an early merger epoch. Gas-rich mergers during this epoch play a major role in fuelling the formation of the old disc stars (the thick disc). This is consistent with studies which show that cold accretion is the main source of a disc galaxy's baryons. Our simulation initially forms a relatively short (scalelength \sim1.7 kpc at z=1) and kinematically hot disc, primarily from gas accreted during the galaxy's merger epoch. Far from being a competing formation scenario, migration is crucial for reconciling the short, hot, discs which form at high redshift in {\Lambda}CDM, with the properties of the thick disc at z=0. The thick disc, as defined by its abundances maintains its relatively short scale-length at z = 0 (2.31 kpc) compared with the total disc scale-length of 2.73 kpc. The inside-out nature of disc growth is imprinted the evolution of abundances such that the metal poor {\alpha}-young population has a larger scale-length (4.07 kpc) than the more chemically evolved metal rich {\alpha}-young population (2.74 kpc).

Within a fully cosmological hydrodynamical simulation, we form a galaxy which rotates at 140 km/s, and is characterised by two loose spiral arms and a bar, indicative of a Hubble Type SBc/d galaxy. We show that our simulated galaxy has no classical bulge, with a pure disc profile at z=1, well after the major merging activity has ended. A long-lived bar subsequently forms, resulting in the formation of a secularly-formed "pseudo" bulge, with the final bulge-to-total light ratio B/T=0.21. We show that the majority of gas which loses angular momentum and falls to the central region of the galaxy during the merging epoch is blown back into the hot halo, with much of it returning later to form stars in the disc. We propose that this mechanism of redistribution of angular momentum via a galactic fountain, when coupled with the results from our previous study which showed why gas outflows are biased to have low angular momentum, can solve the angular momentum/bulgeless disc problem of the cold dark matter paradigm.

Recent numerical simulations have demonstrated that the Sun's dynamical history within the Milky Way may be much more complex than that suggested by its current low peculiar velocity. In particular, the Sun may have radially migrated through the galactic disk by up to 5-6 kpc. This has important ramifications for the structure of the Oort Cloud, as it means that the solar system may have experienced tidal and stellar perturbations that were significantly different from its current local galactic environment. To characterize the effects of solar migration within the Milky Way, we use direct numerical simulations to model the formation of an Oort Cloud around stars that end up on solar-type orbits in a galactic-scale simulation of a Milky Way-like disk formation. Surprisingly, our simulations indicate that Sedna's orbit may belong to the classical Oort Cloud. Contrary to previous understanding, we show that field star encounters play a pivotal role in setting the Oort Cloud's extreme inner edge, and due to their stochastic nature this inner edge sometimes extends to Sedna's orbit. The Sun's galactic migration heightens the chance of powerful stellar passages, and Sedna production occurs around ~20-30% of the solar-like stars we study. Considering the entire Oort Cloud, we find its median distance depends on the minimum galactocentric distance attained during the Sun's orbital history. The inner edge also shows a similar dependence but with increased scatter due to the effects of powerful stellar encounters. Both of these Oort Cloud parameters can vary by an order of magnitude and are usually overestimated by an Oort Cloud formation model that assumes a fixed galactic environment. In addition, the amount of material trapped in outer Oort Cloud orbits can be extremely low and may present difficulties for models of Oort Cloud formation and long-period comet production.

Recent theoretical work suggests that it may be common for stars in the disks of spiral galaxies to migrate radially across significant distances in the disk. Such migrations are a result of resonant scattering with spiral arms and move the guiding centers of the stars while preserving the circularities of their orbits. Migration can therefore efficiently mix stars in all parts of the Galactic disk. We are rapidly approaching an important confluence of theory and observation, where we may soon be able to uncover signatures of such processes in our own Milky Way. The resolution and robustness of the physical modeling in simulations has improved drastically, while observational datasets are increasing in depth and astrometric accuracy. Here, we discuss the results from our idealized N-body/SPH simulations of disk formation and evolution, emphasizing specifically the observational consequences of stellar migration on the solar neighborhood and the vertical structure of the disk. We demonstrate that radial mixing of stars is a crucial dynamical process that we must try to understand if we are to draw significant conclusions about our Galactic history from the properties of stars in our vicinity.

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 toward 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 halos.

Like the majority of spiral galaxies, NGC 6155 exhibits an exponential surface brightness profile that steepens significantly at large radii. Using the VIRUS-P IFU spectrograph, we have gathered spatially resolved spectra of the system. Modifying the GANDALF spectral fitting routine for use on the complex stellar populations found in spirals, we find that the average stellar ages increase significantly beyond the profile break radius. This result is in good agreement with recent simulations that predict the outskirts of disk galaxies are populated through stellar migration. With the ability to bin multiple fibers, we are able to measure stellar population ages down to mu_V~24 mag/sq arcsec.

Stars in disks of spiral galaxies are usually assumed to remain roughly at their birth radii. This assumption is built into decades of modelling of the evolution of stellar populations in our own Galaxy and in external systems. We present results from self-consistent high-resolution $N$-body + Smooth Particle Hydrodynamics simulations of disk formation, in which stars migrate across significant galactocentric distances due to resonant scattering with transient spiral arms, while preserving their circular orbits. We investigate the implications of such migrations for observed stellar populations. Radial migration provides an explanation for the observed flatness and spread in the age-metallicity relation and the relative lack of metal poor stars in the solar neighborhood. The presence of radial migration also prompts rethinking of interpretations of extra-galactic stellar population data, especially for determinations of star formation histories.

The outer disks of galaxies present a unique laboratory for studying the process of disk formation. A considerable fraction of observed disks exhibit a break in their surface brightness profiles. The ubiquity of these features points to a crucial aspect of disk formation which must be explained. Recent theoretical work suggests that such breaks are related to significant amounts of radial migration. We discuss the current observational evidence which supports this picture.

We have performed a high mass and force resolution simulation of an idealized galaxy forming from dissipational collapse of gas embedded in a spherical dark matter halo. The simulation includes star formation and effects of stellar feedback. In our simulation a stellar disk forms with a surface density profile consisting of an inner exponential breaking to a steeper outer exponential. The break forms early on and persists throughout the evolution, moving outward as more gas is able to cool and add mass to the disk. The parameters of the break are in excellent agreement with observations. The break corresponds with a rapid drop in the star formation rate associated with a drop in the cooled gas surface density, but the outer exponential is populated by stars that were scattered outward on nearly circular orbits from the inner disk by spiral arms. The resulting profile and its associated break are therefore a consequence of the interplay between a radial star formation cutoff and redistribution of stellar mass by secular processes. A consequence of such evolution is a sharp change in the radial mean stellar age profile at the break radius.