We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate from the overall cluster dwarf population. UDGs in cluster environments form primarily from dwarf galaxies that experienced early cluster in-fall and subsequent quenching due to ram pressure. The ensuing dimming of these dwarf galaxies due to passive stellar evolution results in a population of very low surface brightness galaxies that are otherwise typical dwarfs. UDGs and non-UDGs alike are affected by tidal interactions with the cluster potential. Tidal stripping of dark matter, as well as mass loss from stellar evolution, results in the adiabatic expansion of stars, particularly in the lowest mass dwarfs. High mass dwarf galaxies show signatures of tidal heating while low mass dwarfs that survive until $z=0$ typically have not experienced such impulsive interactions. There is little difference between UDGs and non-UDGs in terms of their dark matter halos, stellar morphology, colors, and location within the cluster. In most respects cluster UDG and non-UDGs alike are similar to isolated dwarf galaxies, except for the fact that they are typically quenched.

Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relationship (AVR) of h277 is nearly identical to that of the analogous solar-neighborhood measurement in the Milky Way. A crucial element of this success is the simulation's dynamically cold multi-phase ISM, which allows young stars to form with a low velocity dispersion ($\sigma_{\mathrm{birth}}$$\sim 6 - 8\ \mathrm{km}\, \mathrm{s}^{-1}$) at late times. Older stars are born kinematically hotter (i.e., the disk settles over time in an "upside-down" formation scenario), and are subsequently heated after birth. The disk also grows "inside-out", and many of the older stars in the solar neighborhood at $z=0$ are present because of radial mixing. We demonstrate that the evolution of $\sigma_{\mathrm{birth}}$ in h277 can be explained by the same model used to describe the general decrease in velocity dispersion observed in disk galaxies from $z\sim 2-3$ to the present-day, in which the disk evolves in quasi-stable equilibrium and the ISM velocity dispersion decreases over time due to a decreasing gas fraction. Thus, our results tie together local observations of the Milky Way's AVR with observed kinematics of high $z$ disk galaxies.

We use the Romulus25 cosmological simulation volume to identify the largest-ever simulated sample of field ultra-diffuse galaxies (UDGs). At z=0, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. Although they tend to have slightly elevated HI masses, UDGs are consistent with the general isolated dwarf galaxy population and make up ~20% of all field galaxies with 10$^7$<M$_\star$/M$_\odot$<10$^{9}$. The properties of our UDGs are consistent with existing observations of field UDGs, but we predict that many isolated UDGs have been missed by current surveys. Despite their isolation at z=0, the UDGs in our sample are the products of major mergers at z>1. Mergers are no more common in UDG than non-UDG progenitors, but mergers that create UDGs tend to happen earlier, produce a temporary boost in spin, and cause star formation to migrate to the outskirts of galaxies, resulting in lower central star formation rates. The centers of the galaxies fade as their central stellar populations age, but their global star formation rates are maintained through bursts of star formation at larger radii. This formation channel is unique relative to other proposals for UDG formation in isolated galaxies, or for UDGs in groups and clusters, demonstrating that UDGs can potentially be formed through multiple mechanisms.

Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.

We report on high-resolution simulations that explore the orbital decay of massive black hole (MBH) pairs with masses between $10^5$ and $10^7 M_{\odot}$ embedded in a circumnuclear gas disk (CND). An adiabatic equation of state is adopted, with a range of adiabatic indices, which maintains a smooth flow. Mergers between MBHs in this mass range would be detectable by the upcoming Laser Inteferometer Space Antenna (LISA). The orbital evolution is followed from the CND scale ($100$~pc) down to separations of $0.1$--$0.01$~pc at which a circumbinary disk (CBD) could form. The decay is erratic and strongly dependent on the gas flow within the disk, that ultimately determines the net torques experienced by the sinking MBH. Overall, we can identify three different evolutionary stages: (i) an initially slow decay that leads to no significant change in the orbital angular momentum, resulting in some circularization; (ii) a fast migration phase in which the orbital angular momentum decreases rapidly; and (iii) a final, very slow decay phase, in which orbital angular momentum can even increase, and a CBD can form. The fast migration phase owes to disk-driven torques originating primarily from the co-orbital region of the secondary MBH, at a distance of 1--3 Hill radii. We find strong analogies with fast Type III migration for massive planets in protoplanetary disks. The CBD forms only when the decay rate becomes small enough to allow it enough time to carve a cavity around the primary MBH, at scales $\lesssim 1$~pc; when this happens, the MBH separation nearly stalls in our higher-resolution run. We suggest an empirically modified gap-opening criterion that takes into account such timescale effects as well as other deviations from standard assumptions made in the literature. [Abriged]

The dense environment of a galaxy cluster can radically transform the content of in-falling galaxies. Recent observations have found a significant population of active galactic nuclei (AGN) within "jellyfish galaxies," galaxies with trailing tails of gas and stars that indicate significant ram pressure stripping. The relationship between AGN and ram pressure stripping is not well understood. In this letter, we investigate the connection between AGN activity and ram pressure in a fully cosmological setting for the first time using the RomulusC simulation, one of the highest resolution simulations of a galaxy cluster to date. We find unambiguous morphological evidence for ram pressure stripping. For lower mass galaxies (with stellar masses < 10^9.5 solar masses) both star formation and black hole accretion are suppressed by ram pressure before they reach pericenter, whereas for more massive galaxies accretion onto the black hole is enhanced during pericentric passage. Our analysis also indicates that as long as the galaxy retains gas, AGN with higher Eddington ratios are more likely to be the found in galaxies experiencing higher ram pressure. We conclude that prior to quenching star formation, ram pressure triggers enhanced accretion onto the black hole, which then produces heating and outflows due to AGN feedback. AGN feedback may in turn serve to aid in the quenching of star formation in tandem with ram pressure.

The evolution of degenerate matter out of equilibrium is a topic of interest in fields such as condensed matter, nuclear and atomic physics, and increasingly cosmology, including inflaton physics prior to reheating. This follow-up paper extends the recent paper on the super-de Broglie dynamics of pure condensates of non-relativistic identical particles subject to non-local two-body interactions to the dynamics of mixed states. It is found that the two-body correlation function plays an increasingly dynamical role in these systems, driving the development of condensates and distributed phases alike. Examples of distribution and correlation evolution are presented, including instances of collapse, bound and unbound states, and phonons in the bulk. Potential applications are also discussed.

Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-faint regime, we present an improved treatment of the IMF. We implement a self-consistent, stochastically populated IMF in cosmological hydrodynamic simulations. We test our method using high-resolution simulations of a Milky Way halo, run to $z=6$, yielding a sample of nearly 100 galaxies. We also use an isolated dwarf galaxy to investigate the resulting systematic differences in galaxy properties. We find that a stochastic IMF in simulations makes feedback burstier, strengthening feedback, and quenching star formation earlier in small dwarf galaxies. For galaxies in halos with mass $\lesssim10^{8.5}$ M$_\odot$, a stochastic IMF typically leads to lower stellar mass compared to a continuous IMF, sometimes by more than an order of magnitude. We show that existing methods of ensuring discrete supernovae incorrectly determine the mass of the star particle and its associated feedback. This leads to overcooling of surrounding gas, with at least ${\sim}$10 per cent higher star formation and ${\sim}$30 per cent higher cold gas content. Going forward, to accurately model dwarf galaxies and compare to observations, it will be necessary to incorporate a stochastically populated IMF that samples the full spectrum of stellar masses.

The intracluster medium (ICM) is a multi-phase environment, dynamically regulated by Active Galactic Nuclei (AGN), the motions of galaxies through it, and mergers with other clusters. AGN as a central heating source are key to preventing runaway cooling flows, but their role in heating cores in a cosmological context is still poorly understood. The activity of the AGN is strongly linked to star formation, especially in the Brightest Cluster Galaxy (BCG), likely because both rely on cold phase gas. A self-consistent model for AGN and star formation in galaxy clusters thus requires cosmological context, higher resolution, and a careful modeling of cooling and heating balance. In this paper, we use the high-resolution hydrodynamical cosmological simulation of the RomulusC galaxy cluster to study in detail the role of AGN and a major, head-on merger in shaping the cluster core. The unprecedented resolution of the RomulusC simulation captures the multiphase structure of the ICM. The realistic large-scale outflows launched by very small-scale thermal injections, the improved modeling of turbulent diffusion and mixing, and the particle nature of the simulation allow us to carefully separate different heating channels. We show that AGN activity, while efficient at regulating star formation, is incapable of destroying a CC. Instead, that process is facilitated by a head-on, 1:8 mass ratio merger. The merger generates bulk and turbulent motions, which in turn mix high entropy gas generated by AGN and merger driven shocks, turbulent dissipation and sloshing of the ICM by infalling substructures. While central cooling times remain shorter than the Hubble time, restoring a CC is made more difficult by the reduced precipitation rates at larger radii, emphasizing that the AGN-ICM connection is truly a multi-scale problem.

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

We investigate the effects of massive black hole growth on the structural evolution of dwarf galaxies within the Romulus25 cosmological hydrodynamical simulation. We study a sample of 228 central, isolated dwarf galaxies with stellar masses $M_{star} < 10^{10} M_\odot$ and a central BH. We find that the local $M_{BH} - M_{star}$ relation exhibits a high degree of scatter below $M_{star} < 10^{10} M_\odot$, which we use to classify BHs as overmassive or undermassive relative to their host $M_{star}$. Overmassive BHs grow through a mixture of BH mergers and relatively high average accretion rates, while undermassive BHs grow slowly through accretion. We find that isolated dwarf galaxies that host overmassive BHs also follow different evolutionary tracks relative to their undermassive BH counterparts, building up their stars and dark matter earlier and experiencing star formation suppression starting around $z=2$. By $z=0.05$, overmassive BH hosts above $M_{star} > 10^{9} M_\odot$ are more likely to exhibit lower central stellar mass density, lower HI gas content, and lower star formation rates than their undermassive BH counterparts. Our results suggest that overmassive BHs in isolated galaxies above $M_{star} > 10^{9} M_\odot$ are capable of driving feedback, in many cases suppressing and even quenching star formation by late times.

A fraction of dwarf galaxies in the Virgo cluster contain disk features like bars and spiral arms. Using $N$-body simulations, we investigate the effects of tidal forces on the formation of such disk features in disk dwarf galaxies resembling VCC856. We consider 8 Cluster-Galaxy models in which disk dwarf galaxies with differing pericenter distance and spin orientation experience the tidal gravitational force of a Virgo-like NFW halo, and additional 8 Galaxy-Galaxy models in which two dwarf galaxies undergo tidal interactions with different strength. We find that the cluster tidal effect is moderate due to the small galaxy size, making the bars form earlier by $\sim1$--$1.5\Gyr$ compared to the cases in isolation. While the galactic halos significantly lose their mass within the virial radius due to the cluster tidal force, the mass of the stellar disks is nearly unchanged, suggesting that the inner regions of a disk-halo system is secured from the tidal force. The tidal forcing from either the cluster potential or a companion galaxy triggers the formation of two-armed spirals at early time before a bar develops. The tidally-driven arms decay and wind with time, suggesting that they are kinematic density waves. In terms of the strength and pitch angle, the faint arms in VCC856 are best matched with the arms in a marginally unstable galaxy produced by a distant tidal encounter with its neighbor $\sim0.85\Gyr$ ago.

We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development.

We investigate planetesimal accretion via a direct N-body simulation of an annulus at 1 AU orbiting a 1 $M_{\odot}$ star. The planetesimal ring, which initially contains N = $10^6$ bodies is evolved into the oligarchic growth phase. Unlike previous lower resolution studies, we find that the mass distribution of planetesimals develops a bump at intermediate mass after the oligarchs form. This feature marks a boundary between growth modes. The smallest planetesimals are packed tightly enough together to populate mean motion resonances with the oligarchs, which heats the small bodies, enhancing their growth. If we depopulate most of the resonances by decreasing the width of the annulus, this effect becomes weaker. To clearly demonstrate the dynamics driving these growth modes, we also examine the evolution of a planetary embryo embedded in an annulus of collisionless planetesimals. In this case, we find that the resonances push planetesimals away from the embryo, decreasing the surface density of the bodies adjacent to the embryo. This effect only occurs when the annulus is wide enough and the mass resolution of the planetesimals is fine enough to populate the resonances. The bump we observe in the mass distribution resembles the 100 km power law break seen in the size distribution of asteroid belt objects. Although the bump produced in our simulations occurs at a size larger than 100 km, we show that the bump location is sensitive to the initial planetesimal mass, which implies that this feature is potentially useful for constraining planetesimal formation models.

We study the link between supermassive black hole growth and the stellar mass assembly of their host galaxies in the state-of-the-art Romulus suite of simulations. The cosmological simulations Romulus25 and RomulusC employ innovative recipes for the seeding, accretion, and dynamics of black holes in the field and cluster environments respectively. We find that the black hole accretion rate traces the star formation rate among star-forming galaxies. This result holds for stellar masses between 10^8 and 10^12 solar masses, with a very weak dependence on host halo mass or redshift. The inferred relation between accretion rate and star formation rate does not appear to depend on environment, as no difference is seen in the cluster/proto-cluster volume compared to the field. A model including the star formation rate, the black hole-to-stellar mass ratio, and the cold gas fraction can explain about 70 per cent of all variations in the black hole accretion rate among star forming galaxies. Finally, bearing in mind the limited volume and resolution of these cosmological simulations, we find no evidence for a connection between black hole growth and galaxy mergers, on any timescale and at any redshift. Black holes and their galaxies assemble in tandem in these simulations, regardless of the larger-scale intergalactic environment, suggesting that black hole growth simply follows star formation on galactic scales.

The first paper in this series showed that QCD axion dark matter, as a highly-correlated Bose fluid, contains extra-classical physics on cosmological scales. The nature of the derived extra-classical physics is exchange-correlation interactions induced by the constraints of symmetric particle exchange and inter-axion correlations from self-gravitation. The paper also showed that the impact of extra-classical physics on early structure formation is marginal, as the exchange-correlation interaction is inherently non-linear. This paper continues the study of axion structure formation into the non-linear regime by now considering the case of full collapse and virialization. The N-body method is chosen to solve the collapse problem, and its algorithms are derived for a condensed Bose fluid. Simulations of isolated gravitational collapse are performed for both Bose and cold dark matter fluids using a prototype N-body code. Unique Bose structures are found to survive even the most violent collapses. Bose post-collapse features include dynamical changes to global structures, creation of new broad substructures, violations of classical binding energy conditions, and new fine structures. Effective models of the novel structures are constructed and possibilities for their discovery are discussed.

Quasar absorption-line studies in the ultraviolet (UV) can uniquely probe the nature of the multiphase cool-warm (10^4 < T < 10^6 K) gas in and around galaxy clusters, promising to provide unprecedented insights into 1) interactions between the circumgalactic medium (CGM) associated with infalling galaxies and the hot (T > 10^6 K) X-ray emitting intracluster medium (ICM), 2) the stripping of metal-rich gas from the CGM, and 3) a multiphase structure of the ICM with a wide range of temperatures and metallicities. In this work, we present results from a high-resolution simulation of a ~10^14 solar mass galaxy cluster to study the physical properties and observable signatures of this cool-warm gas in galaxy clusters. We show that the ICM becomes increasingly multiphased at large radii, with the cool-warm gas becoming dominant in cluster outskirts. The diffuse cool-warm gas also exhibits a wider range of metallicity than the hot X-ray emitting gas. We make predictions for the covering fractions of key absorption-line tracers, both in the ICM and in the CGM of cluster galaxies, typically observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope (HST). We further extract synthetic spectra to demonstrate the feasibility of detecting and characterizing the thermal, kinematic, and chemical composition of the cool-warm gas using H I, O VI, and C IV lines, and we predict an enhanced population of broad Ly-alpha absorbers tracing the warm gas. Lastly, we discuss future prospects of probing the multiphase structure of the ICM beyond HST.

Compact ellipticals (cEs) are outliers from the scaling relations of early-type galaxies, particularly the mass-metallicity relation which is an important outcome of feedback. The formation of such low-mass, but metal-rich and compact, objects is a long-standing puzzle. Using a pair of high-resolution N-body+gas simulations, we investigate the evolution of a gas-rich low-mass galaxy on a highly radial orbit around a massive host galaxy. As the infalling low-mass galaxy passes through the host's corona at supersonic speeds, its diffuse gas outskirts are stripped by ram pressure, as expected. However, the compactness increases rapidly because of bursty star formation in the gas tidally driven to the centre. The metal-rich gas produced by supernovae and stellar winds is confined by the ram pressure from the surrounding environment, leading to subsequent generations of stars being more metal-rich. After the gas is depleted, tidal interactions enhance the metallicity further via the stripping of weakly bound, old, and metal-poor stars, while the size of the satellite is changed only modestly. The outcome is a metal-rich cE that is consistent with observations. These results argue that classical cEs are neither the stripped remnants of much more massive galaxies nor the merger remnants of normal dwarfs. We present observable predictions that can be used to test our model.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

The existence of ultra-faint dwarf (UFD) galaxies highlights the need to push our theoretical understanding of galaxies to extremely low mass. We examine the formation of UFDs by twice running a fully cosmological simulations of dwarf galaxies, but varying star formation. One run uses a temperature-density threshold for star formation, while the other uses an H$_{2}$-based sub-grid star formation model. The total number of dwarf galaxies that forms is different by a factor of 2 between the two runs, but most of these are satellites, leading to a factor of 5 difference in the number of luminous UFD companions around more massive, isolated dwarfs. The first run yields a 47\% chance of finding a satellite around a M$_{halo}$ $\sim 10^{10}$ M$_{\odot}$ host, while the H$_2$ run predicts only a 16\% chance. Metallicity is the primary physical parameter that creates this difference. As metallicity decreases, the formation of H$_2$ is slowed and relegated to higher-density material. Thus, our H$_2$ run is unable to form many (and often, any) stars before reionization removes gas. These results emphasize that predictions for UFD properties made using hydrodynamic simulations, in particular regarding the frequency of satellites around dwarf galaxies, the slope of the stellar mass function at low masses, as well as the properties of ultra-faint galaxies occupying the smallest halos, are extremely sensitive to the subgrid physics of star formation contained within the simulation. However, upcoming discoveries of ultra-faint dwarfs will provide invaluable constraining power on the physics of the first star formation.

The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Way's thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies.

Star formation is observed to be strongly correlated to dense regions of molecular gas. Although the exact nature of the link between star formation and molecular hydrogen is still unclear, some have suggested that shielding of dense gas by dust grains is the key factor enabling the presence of both. We present a sub-grid model for use in galaxy formation simulations in which star formation is linked explicitly to local dust shielding. We developed and tested our shielding and star formation models using smoothed particle hydrodynamic simulations of solar and sub-solar metallicity isolated Milky Way-mass disk galaxies. We compared our dust shielding-based star formation model to two other star formation recipes that used gas temperature and H$_2$ fraction as star formation criteria. We further followed the evolution of a dwarf galaxy within a cosmological context using both the shielding and H$_2$-based star formation models. We find that the shielding-based model allows for star formation at higher temperatures and lower densities than a model in which star formation is tied directly to H$_2$ abundance, as requiring H$_2$ formation leads the gas to undergo additional gravitational collapse before star formation. However, the resulting galaxies are very similar for both the shielding and H$_2$-based star formation models, and both models reproduce the resolved Kennicutt-Schmidt law. Therefore, both star formation models appear viable in the context of galaxy formation simulations.

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies.

We examine the effects of SMBH feedback on the CGM using a cosmological hydrodynamic simulation \citep[{\sc Romulus25};][]{Tremmel2017} and a set of four zoom-in `genetically modified' Milky Way-mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the circumgalactic medium (CGM), we show that \ion{O}{6} is a sensitive indicator of supermassive black hole (SMBH) feedback. First, we calculate the column densities of \ion{O}{6} in simulated Milky Way-mass galaxies and compare them with observations from the COS-Halos Survey. Our simulations show column densities of \ion{O}{6} in the CGM consistent with those of COS-Halos star forming and quenched galaxies. These results contrast with those from previous simulation studies which typically underproduce CGM column densities of \ion{O}{6}. We determine that a galaxy's star formation history and assembly record have little effect on the amount of \ion{O}{6} in its CGM. Instead, column densities of \ion{O}{6} are closely tied to galaxy halo mass and BH growth history. The set of zoom-in, genetically modified Milky Way-mass galaxies indicates that the SMBH drives highly metal-enriched material out into its host galaxy's halo which in turn elevates the column densities of \ion{O}{6} in the CGM.

Axions as dark matter is an increasingly important subject in astrophysics and cosmology. Experimental and observational searches are mounting across the mass spectrum of axion-like particles, many of which require detailed knowledge of axion structure over a wide range of scales. Current understanding of axion structures is far from complete, however, due largely to controversy in modeling the candidate's highly-degenerate state. The series Axion Structure Formation seeks to develop a consistent model of QCD axion dark matter dynamics that follows their highly-degenerate nature to the present using novel modeling techniques and sophisticated simulations. This inaugural paper presents the problem of describing many non-relativistic axions with minimal degrees of freedom and constructs a theory of axion infall for the limit of complete condensation. The derived model is shown to contain axion-specific dynamics not unlike the exchange-correlation influences experienced by identical fermions. Perturbative calculations are performed to explore the potential for imprints in early universe structures.

The majority of galactic baryons reside outside of the galactic disk in the diffuse gas known as the circumgalactic medium (CGM). While state-of-the art simulations excel at reproducing galactic disk properties, many struggle to drive strong galactic winds or to match the observed ionization structure of the CGM using only thermal supernova feedback. To remedy this, recent studies have invoked non-thermal cosmic ray (CR) stellar feedback prescriptions. However, numerical schemes of CR transport are still poorly constrained. We explore how the choice of CR transport affects the multiphase structure of the simulated CGM. We implement anisotropic CR physics in the astrophysical simulation code, {\sc Enzo} and simulate a suite of isolated disk galaxies with varying prescriptions for CR transport: isotropic diffusion, anisotropic diffusion, and streaming. We find that all three transport mechanisms result in strong, metal-rich outflows but differ in the temperature and ionization structure of their CGM. Isotropic diffusion results in a spatially uniform, warm CGM that underpredicts the column densities of low-ions. Anisotropic diffusion develops a reservoir of cool gas that extends further from the galactic center, but disperses rapidly with distance. CR streaming projects cool gas out to radii of 200 kpc, supporting a truly multiphase medium. In addition, we find that streaming is less sensitive to changes in constant parameter values like the CR injection fraction, transport velocity, and resolution than diffusion. We conclude that CR streaming is a more robust implementation of CR transport and motivate the need for detailed parameter studies of CR transport.

Recent discoveries of massive black holes (MBHs) in dwarf galaxies suggest that they may have a more common presence than once thought. Systematic searches are revealing more candidates, but this process could be accelerated by predictions from simulations. We perform a study of several high-resolution, cosmological, zoom-in simulations focusing on dwarf galaxies that host massive black holes at z = 0, with the aim of determining when the black holes are most observable. Larger dwarf galaxies are more likely to host MBHs than those of lower mass. About 50% of the MBHs in dwarfs are not centrally located, but rather are wandering within a few kpc of the galaxy center. The accretion luminosities of MBHs in dwarfs are low throughout cosmic time, rendering them extremely difficult to detect. However, the merger history of these MBHs is optimal for gravitational wave detection by LISA.

We analyze the metal accumulation in dwarf and spiral galaxies by following the history of metal enrichment and outflows in a suite of twenty high-resolution simulated galaxies. These simulations agree with the observed stellar and gas-phase mass-metallicity relation, an agreement that relies on large fractions of the produced metals escaping into the CGM. For instance, in galaxies with Mvir ~ 1e9.5 -- 1e10 solar masses, we find that about ~ 85% of the available metals are outside of the galactic disk at z = 0, although the fraction decreases to a little less than half in Milky Way-mass galaxies. In many cases, these metals are spread far beyond the virial radius. We analyze the metal deficit within the ISM and stars in the context of previous work tracking the inflow and outflow of baryons. Outflows are prevalent across the entire mass range, as is reaccretion. We find that between 40 and 80% of all metals removed from the galactic disk are later reaccreted. The outflows themselves are metal enriched relative to the ISM by a factor of 0.2 dex because of the correspondence between sites of metal enrichment and outflows. As a result, the metal mass loading factor scales as eta_metals \propto v_circ^-0.91, a somewhat shallower scaling than the total mass loading factor. We analyze the simulated galaxies within the context of analytic chemical evolution models by determining their net metal expulsion efficiencies, which encapsulate the rates of metal loss and reaccretion. We discuss these results in light of the inflow and outflow properties necessary for reproducing the mass-metallicity relation.

Systems of identical particles possessing non-local interactions are capable of exhibiting extra-classical properties beyond the characteristic quantum length scales. This letter derives the dynamics of such systems in the non-relativistic and degenerate limit, showing the effect of exchange symmetry and correlations on structure both in and out of equilibrium. Such descriptions may be crucial to understanding systems ranging from nuclei to dark matter. Appropriate limits for restoring the mean-field description are also discussed.

We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. The critical flux required to inhibit H$_2$ formation, $J_{\rm crit}$, is hotly debated, largely due to the uncertainties in the source radiation spectrum, H$_2$ self-shielding, and collisional dissociation rates. Here, we test the power of the direct collapse model in a self-consistent, time-dependant, non-uniform Lyman-Werner radiation field -- the first time such has been done in a cosmological volume -- using an updated version of the SPH+N-body tree code Gasoline with H$_2$ non-equilibrium abundance tracking, H$_2$ cooling, and a modern SPH implementation. We vary $J_{\rm crit} $ from $30$ to $10^3$ in units of $J_{21}$ to study how this parameter impacts the number of seed black holes and the type of galaxies which host them. We focus on black hole formation as a function of environment, halo mass, metallicity, and proximity of the Lyman-Werner source. Massive black hole seeds form more abundantly with lower $J_{\rm crit}$ thresholds, but regardless of $J_{\rm crit}$, these seeds typically form in halos that have recently begun star formation. Our results do not confirm the proposed atomic cooling halo pair scenario; rather black hole seeds predominantly form in low-metallicity pockets of halos which already host star formation.

Using a new, high-resolution cosmological hydrodynamic simulation of a Milky Way-type (MW-type) galaxy, we explore how a merger-rich assembly history affects the mass budget of the central supermassive black hole (SMBH). We examine a MW-mass halo at the present epoch whose evolution is characterized by several major mergers to isolate the importance of merger history on black hole accretion. This study is an extension of Bellovary et. al. 2013, which analyzed the accretion of high mass, high redshift galaxies and their central black holes, and found that the gas content of the central black hole reflects what is accreted by the host galaxy halo. In this study, we find that a merger-rich galaxy will have a central SMBH preferentially fed by merger gas. Moreover, we find that nearly 30$\%$ of the accreted mass budget of the SMBH enters the galaxy through the two major mergers in its history, which may account for the increase of merger-gas fueling the SMBH. Through an investigation of the angular momentum of the gas entering the host and its SMBH, we determine that merger gas enters the galaxy with lower angular momentum compared to smooth accretion, partially accounting for the preferential fueling witnessed in the SMBH. In addition, the presence of mergers, particularly major mergers, also helps funnel low angular momentum gas more readily to the center of the galaxy. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories.

We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as v_circ^(-2.2), with an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically about 1 Gyr, virtually independent of halo mass. Recycled material is re-accreted farther out in the disk and with typically about 2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most of cosmic star formation occurs.

Using an energy balance model with ice sheets, we examine the climate response of an Earth-like planet orbiting a G dwarf star and experiencing large orbital and obliquity variations. We find that ice caps couple strongly to the orbital forcing, leading to extreme ice ages. In contrast with previous studies, we find that such exo-Milankovitch cycles tend to impair habitability by inducing snowball states within the habitable zone. The large amplitude changes in obliquity and eccentricity cause the ice edge, the lowest latitude extent of the ice caps, to become unstable and grow to the equator. We apply an analytical theory of the ice edge latitude to show that obliquity is the primary driver of the instability. The thermal inertia of the ice sheets and the spectral energy distribution of the G dwarf star increase the sensitivity of the model to triggering runaway glaciation. Finally, we apply a machine learning algorithm to demonstrate how this technique can be used to extend the power of climate models. This work illustrates the importance of orbital evolution for habitability in dynamically rich planetary systems. We emphasize that as potentially habitable planets are discovered around G dwarfs, we need to consider orbital dynamics.

We present a self-consistent prediction from a large-scale cosmological simulation for the population of `wandering' supermassive black holes (SMBHs) of mass greater than $10^6$ M$_{\odot}$ on long-lived, kpc-scale orbits within Milky Way (MW)-mass galaxies. We extract a sample of MW-mass halos from the Romulus25 cosmological simulation (Tremmel et al. 2017), which is uniquely able to capture the orbital evolution of SMBHs during and following galaxy mergers. We predict that such halos, regardless of recent merger history or morphology, host an average of $5.1 \pm 3.3$ SMBHs, including their central black hole, within 10 kpc from the galactic center and an average of $12.2 \pm 8.4$ SMBHs total within their virial radius, not counting those in satellite halos. Wandering SMBHs exist within their host galaxies for several Gyrs, often accreted by their host halo in the early Universe. We find, with $>4\sigma$ significance, that wandering SMBHs are preferentially found outside of galactic disks.

We outline a mechanism that explains the observed lack of circumbinary planets (CBPs) via coupled stellar-tidal evolution of isolated binary stars. Tidal forces between low-mass, short-period binary stars on the pre-main sequence slow the stellar rotations, transferring rotational angular momentum to the orbit as the stars approach the tidally locked state. This transfer increases the binary orbital period, expanding the region of dynamical instability around the binary, and destabilizing CBPs that tend to preferentially orbit just beyond the initial dynamical stability limit. After the stars tidally lock, we find that angular momentum loss due to magnetic braking can significantly shrink the binary orbit, and hence the region of dynamical stability, over time impacting where surviving CBPs are observed relative to the boundary. We perform simulations over a wide range of parameter space and find that the expansion of the instability region occurs for most plausible initial conditions and that in some cases, the stability semi-major axis doubles from its initial value. We examine the dynamical and observable consequences of a CBP falling within the dynamical instability limit by running N-body simulations of circumbinary planetary systems and find that typically, at least one planet is ejected from the system. We apply our theory to the shortest period {\it Kepler} binary that possesses a CBP, Kepler-47, and find that its existence is consistent with our model. Under conservative assumptions, we find that coupled stellar-tidal evolution of pre-main sequence binary stars removes at least one close-in CBP in $87\%$ of multi-planet circumbinary systems.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We analyze the evolution of the potentially habitable planet Proxima Centauri b to identify environmental factors that affect its long-term habitability. We consider physical processes acting on size scales ranging from the galactic to the stellar system to the planet's core. We find that there is a significant probability that Proxima Centauri has had encounters with its companion stars, Alpha Centauri A and B, that are close enough to destabilize an extended planetary system. If the system has an additional planet, as suggested by the discovery data, then it may perturb planet b's eccentricity and inclination, possibly driving those parameters to non-zero values, even in the presence of strong tidal damping. We also model the internal evolution of the planet, evaluating the roles of different radiogenic abundances and tidal heating and find that magnetic field generation is likely for billions of years. We find that if planet b formed in situ, then it experienced 169 +/- 13 million years in a runaway greenhouse as the star contracted during its formation. This early phase could remove up to 5 times as much water as in the modern Earth's oceans, possibly producing a large abiotic oxygen atmosphere. On the other hand, if Proxima Centauri b formed with a substantial hydrogen atmosphere (0.01 - 1% of the planet's mass), then this envelope could have shielded the water long enough for it to be retained before being blown off itself. After modeling this wide range of processes we conclude that water retention during the host star's pre-main sequence phase is the biggest obstacle for Proxima b's habitability. These results are all obtained with a new software package called VPLANET.

The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.

How can scientists conclude with high confidence that an exoplanet hosts life? As telescopes come on line over the next 20 years that can directly observe photons from terrestrial exoplanets, this question will dictate the activities of many scientists across many fields. The expected data will be sparse and with low signal-to-noise, which will make disentangling biosignatures from abiotic features challenging. Our Earth is not just unique in that it hosts life, it is also the only terrestrial planet with direct observations of its interior through seismic waves, and compositional evolution through field and laboratory measurements. This extensive research reveals a planet born from collisions between worlds, followed by a complicated biogeochemical evolution. Exoplanet interiors, on the other hand, can only be constrained by the following observations: 1) photometric and spectroscopic analysis of the planet's atmosphere, 2) spectroscopic and photometric analysis of the host star, and 3) companion planet properties. From these (future) data, astrobiologists must generate plausible compositional and evolutionary models that constrain a potentially habitable exoplanet's internal properties and history, provide environmental context, and rule out geochemical explanations for any putative biosignatures. The goal of this white paper is to frame the role of geophysical and geochemical processes relevant to the search for life beyond the Solar System and to identify critical, but understudied, areas of future research.

We present the first self-consistent prediction for the distribution of formation timescales for close Supermassive Black Hole (SMBH) pairs following galaxy mergers. Using ROMULUS25, the first large-scale cosmological simulation to accurately track the orbital evolution of SMBHs within their host galaxies down to sub-kpc scales, we predict an average formation rate density of close SMBH pairs of 0.013 cMpc^-3 Gyr^-1. We find that it is relatively rare for galaxy mergers to result in the formation of close SMBH pairs with sub-kpc separation and those that do form are often the result of Gyrs of orbital evolution following the galaxy merger. The likelihood and timescale to form a close SMBH pair depends strongly on the mass ratio of the merging galaxies, as well as the presence of dense stellar cores. Low stellar mass ratio mergers with galaxies that lack a dense stellar core are more likely to become tidally disrupted and deposit their SMBH at large radii without any stellar core to aid in their orbital decay, resulting in a population of long-lived 'wandering' SMBHs. Conversely, SMBHs in galaxies that remain embedded within a stellar core form close pairs in much shorter timescales on average. This timescale is a crucial, though often ignored or very simplified, ingredient to models predicting SMBH mergers rates and the connection between SMBH and star formation activity.

The obliquity of the Earth, which controls our seasons, varies by only ~2.5 degrees over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth's ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities, high eccentricities, and dynamical variations can extend the outer edge of the habitable zone by preventing runaway glaciation (snowball states). We expand upon these studies by exploring the orbital dynamics with a semi-analytic model that allows us to map broad regions of parameter space. We find that in general, the largest drivers of obliquity variations are secular spin-orbit resonances. We show how the obliquity varies in several test cases, including Kepler-62 f, across a wide range of orbital and spin parameters. These obliquity variations, alongside orbital variations, will have a dramatic impact on the climates of such planets.

Within a Lambda Cold Dark Matter (LCDM) scenario, we use high resolution cosmological simulations spanning over four orders of magnitude in galaxy mass to understand the deficit of dwarf galaxies in observed velocity functions. We measure velocities in as similar a way as possible to observations, including generating mock HI data cubes for our simulated galaxies. We demonstrate that this apples-to-apples comparison yields an "observed" velocity function in agreement with observations, reconciling the large number of low-mass halos expected in a LCDM cosmological model with the low number of observed dwarfs at a given velocity. We then explore the source of the discrepancy between observations and theory, and conclude that the dearth of observed dwarf galaxies is primarily explained by two effects. The first effect is that galactic rotational velocities derived from the HI linewidth severely underestimate the maximum halo velocity. The second effect is that a large fraction of halos at the lowest masses are too faint to be detected by current galaxy surveys. We find that cored dark matter density profiles can contribute to the lower observed velocity of galaxies, but only for galaxies in which the velocity is measured interior to the size of the core (~3 kpc).

We describe the structure and implementation of a moving-mesh hydrodynamics solver in the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). While largely based on the algorithm described by Springel (2010) that is implemented in AREPO, our algorithm differs a few aspects. We describe our use of the Voronoi tessellation library, VORO++, to compute the Voronoi tessellation directly. We also incorporate some recent advances in gradient estimation and reconstruction that gives better accuracy in hydrodynamic solutions at minimal computational cost. We validate this module with a small battery of test problems against the smooth particle hydrodynamics solver included in ChaNGa. Finally, we study one example of a scientific problem involving the mergers of two main sequence stars and highlight the small quantitative differences between smooth particle and moving-mesh hydrodynamics. We close with a discussion of anticipated future improvements and advancements.

The methods in the Gasoline2 Smoothed Particle Hydrodynamics (SPH) code are described and tested. Gasoline2 is the most recent version of the Gasoline code for parallel hydrodynamics and gravity with identical hydrodynamics to the Changa code. As with other Modern SPH codes, we prevent sharp jumps in time steps, use upgraded kernels and larger neighbour numbers and employ local viscosity limiters. Unique features in Gasoline2 include its Geometric-Density-Average Force expression, explicit Turbulent Diffusion terms and Gradient-Based shock detection to limit artificial viscosity. This last feature allows Gasoline2 to completely avoid artificial viscosity in non-shocking compressive flows. We present a suite of tests demonstrating the value of these features with the same code configuration and parameter choices used for production simulations.

We present a novel implementation of supermassive black hole (SMBH) formation, dynamics, and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modeling of SMBHs in fully cosmological simulations, allowing for a more de- tailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies, and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multi-dimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization into a new set of high resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses, and global SF and SMBH accretion histories at high redshift that are consistent with observations. We discuss the importance of SMBH physics in shaping the evolution of massive galaxies and show how SMBH feedback is much more effective at regulating star formation compared to SNe feedback in this regime. Further, we show how each aspect of our SMBH model impacts this evolution compared to more common approaches. Finally, we present a science application of this scheme studying the properties and time evolution of an example dual AGN system, highlighting how our approach allows simulations to better study galaxy interactions and SMBH mergers in the context of galaxy-BH co-evolution.

Using a cosmological $N$-body simulation, we investigate the origin and distribution of stars in the intracluster light (ICL) of a Fornax-like cluster. In a dark matter only simulation we identify a halo which, at $z=0$, has $M_200 \simeq 4.1 \times 10^{13}M_{sun}$ and $r_{200} = 700kpc$, and replace infalling subhalos with models that include spheroid and disc components. As they fall into the cluster, the stars in some of these galaxies are stripped from their hosts, and form the ICL. We consider the separate contributions to the ICL from stars which originate in the haloes and the discs of the galaxies. We find that disc ICL stars are more centrally concentrated than halo ICL stars. The majority of the disc ICL stars are associated with one initially disc-dominated galaxy that falls to the centre of the cluster and is heavily disrupted, producing part of the cD galaxy. At radial distances greater than 200kpc, well beyond the stellar envelope of the cD galaxy, stars formerly from the stellar haloes of galaxies dominate the ICL. Therefore at large distances, the ICL population is dominated by older stars.

We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies, run with the N-body + SPH code, ChaNGA, at a high resolution of 60 pc. For M$_{halo} > 10^9$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter. However, the scatter in the SMHM relation increases dramatically for lower-mass halos. We find that some of this scatter is due to {\em dark dwarfs}, halos devoid of stars. However, even when only considering well-resolved halos that contain a stellar population, the scatter in stellar mass reaches nearly 1 dex for M$_{halo}$($z=0$) 10$^7$ M$_{\odot}$. Much of this scatter is due to including satellites of the dwarf galaxies that have had their halo masses reduced through tidal stripping. The fraction of dark dwarfs (those that contain no stars) increases substantially with decreasing halo mass. When these dark halos are considered, the true scatter in the SMHM at low masses is even larger. At the faintest end of the SMHM relation probed by our simulations, a galaxy cannot be assigned a unique halo mass based solely on its luminosity. We provide a formula to stochastically populate low-mass halos following our results. Our predicted large scatter at low halo masses increases the slope of the resulting stellar mass function on the ultra-faint dwarf galaxy scales currently being probed by such surveys as the Dark Energy Survey or the Hyper-Suprime Cam Subaru Strategic Program, and in the future by the Large Synoptic Survey Telescope.

The sources that reionized the universe are still unknown, but likely candidates are faint but numerous galaxies. In this paper we present results from running a high resolution, uniform volume simulation, the Vulcan, to predict the number densities of undetectable, faint galaxies and their escape fractions of ionizing radiation, $f_\mathrm{esc}$, during reionization. Our approach combines a high spatial resolution, a realistic treatment of feedback and hydro processes, a strict threshold for minimum number of resolution elements per galaxy, and a converged measurement of $f_\mathrm{esc}$. We calibrate our physical model using a novel approach to create realistic galaxies at z=0, so the simulation is predictive at high redshifts. With this approach we can (1) robustly predict the evolution of the galaxy UV luminosity function at faint magnitudes down to $M_\mathrm{UV}$~-15, two magnitudes fainter than observations, and (2) estimate $f_\mathrm{esc}$ over a large range of galaxy masses based on the detailed stellar and gas distributions in resolved galaxies. We find steep faint end slopes, implying high number densities of faint galaxies, and the dependence of $f_\mathrm{esc}$ on the UV magnitude of a galaxy, given by the power-law: log $f_\mathrm{esc} = (0.51 \pm 0.04)M_\mathrm{UV} + 7.3 \pm 0.8$, with the faint population having $f_\mathrm{esc}$~35%. Convolving the UV luminosity function with $f_\mathrm{esc}$($M_\mathrm{UV}$), we find an ionizing emissivity that is (1) dominated by the faintest galaxies and (2) reionizes the universe at the appropriate rate, consistent with observational constraints of the ionizing emissivity and the optical depth to the decoupling surface tau_es, without the need for additional sources of ionizing radiation.

Signal estimates for direct axion dark matter searches have used the isothermal sphere halo model for the last several decades. While insightful, the isothermal model does not capture effects from a halo's infall history nor the influence of baryonic matter, which has been shown to significantly influence a halo's inner structure. The high resolution of cavity axion detectors can make use of modern cosmological structure-formation simulations, which begin from realistic initial conditions, incorporate a wide range of baryonic physics, and are capable of resolving detailed structure. This letter uses a state-of-the-art cosmological N-body+Smoothed-Particle Hydrodynamics simulation to develop an improved signal model for axion cavity searches. Signal shapes from a class of galaxies encompassing the Milky Way are found to depart significantly from the isothermal sphere. A new signal model for axion detectors is proposed and projected sensitivity bounds on the Axion Dark Matter eXperiment data are presented.

The recent discoveries of circumbinary planets by $\it Kepler$ raise questions for contemporary planet formation models. Understanding how these planets form requires characterizing their formation environment, the circumbinary protoplanetary disc, and how the disc and binary interact and change as a result. The central binary excites resonances in the surrounding protoplanetary disc that drive evolution in both the binary orbital elements and in the disc. To probe how these interactions impact binary eccentricity and disc structure evolution, N-body smooth particle hydrodynamics (SPH) simulations of gaseous protoplanetary discs surrounding binaries based on Kepler 38 were run for $10^4$ binary periods for several initial binary eccentricities. We find that nearly circular binaries weakly couple to the disc via a parametric instability and excite disc eccentricity growth. Eccentric binaries strongly couple to the disc causing eccentricity growth for both the disc and binary. Discs around sufficiently eccentric binaries that strongly couple to the disc develop an $m = 1$ spiral wave launched from the 1:3 eccentric outer Lindblad resonance (EOLR) that corresponds to an alignment of gas particle longitude of periastrons. All systems display binary semimajor axis decay due to dissipation from the viscous disc.

We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $\Lambda$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.

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. %

Circumplanetary disks can be found around forming giant planets, regardless of whether core accretion or gravitational instability built the planet. We carried out state-of-the-art hydrodynamical simulations of the circumplanetary disks for both formation scenarios, using as similar initial conditions as possible to unveil possible intrinsic differences in the circumplanetary disk mass and temperature between the two formation mechanisms. We found that the circumplanetary disks mass linearly scales with the circumstellar disk mass. Therefore, in an equally massive protoplanetary disk, the circumplanetary disks formed in the disk instability model can be only a factor of eight more massive than their core-accretion counterparts. On the other hand, the bulk circumplanetary disk temperature differs by more than an order of magnitude between the two cases. The subdisks around planets formed by gravitational instability have a characteristic temperature below 100 K, while the core accretion circumplanetary disks are hot, with temperatures even greater than 1000 K when embedded in massive, optically thick protoplanetary disks. We explain how this difference can be understood as the natural result of the different formation mechanisms. We argue that the different temperatures should persist up to the point when a full-fledged gas giant forms via disk instability, hence our result provides a convenient criteria for observations to distinguish between the two main formation scenarios by measuring the bulk temperature in the planet vicinity.

Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here we use 1D coupled climate-photochemical models to generate self-consistent atmospheres for evolutionary scenarios predicted in our companion paper (Barnes et al., 2016). These include high-O2, high-CO2, and more Earth-like atmospheres, with either oxidizing or reducing compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and instrument models to explore our ability to discriminate between possible planetary states. These results are applicable not only to Proxima Cen b, but to other terrestrial planets orbiting M dwarfs. Thermal phase curves may provide the first constraint on the existence of an atmosphere, and JWST observations longward of 7 microns could characterize atmospheric heat transport and molecular composition. Detection of ocean glint is unlikely with JWST, but may be within the reach of larger aperture telescopes. Direct imaging spectra may detect O4, which is diagnostic of massive water loss and O2 retention, rather than a photosynthesis. Similarly, strong CO2 and CO bands at wavelengths shortward of 2.5 {\mu}m would indicate a CO2-dominated atmosphere. If the planet is habitable and volatile-rich, direct imaging will be the best means of detecting habitability. Earth-like planets with microbial biospheres may be identified by the presence of CH4 and either photosynthetically produced O2 or a hydrocarbon haze layer.

We investigate the conditions required for planet formation via gravitational instability (GI) and protoplanetary disk (PPD) fragmentation around M-dwarfs. Using a suite of 64 SPH simulations with $10^6$ particles, the parameter space of disk mass, temperature, and radius is explored, bracketing reasonable values based on theory and observation. Our model consists of an equilibrium, gaseous, and locally isothermal disk orbiting a central star of mass $M_*=M_{sol}/3$. Disks with a minimum Toomre $Q$ of $Q_{min} \lesssim 0.9$ will fragment and form gravitationally bound clumps. Some previous literature has found $Q_{min} < 1.3-1.5$ to be sufficient for fragmentation. Increasing disk height tends to stabilize disks, and when incorporated into $Q$ as $Q_{eff}\propto Q(H/R)^\alpha$ for $\alpha=0.18$ is sufficient to predict fragmentation. Some discrepancies in the literature regarding $Q_{crit}$ may be due to different methods of generating initial conditions (ICs). A series of 15 simulations demonstrates that perturbing ICs slightly out of equilibrium can cause disks to fragment for higher $Q$. Our method for generating ICs is presented in detail. We argue that GI likely plays a role in PPDs around M-dwarfs and that disk fragmentation at large radii is a plausible outcome for these disks.

The "too big to fail" problem is revisited by studying the tidal evolution of populations of dwarf satellites with different density profiles. The high resolution cosmological $\rm \Lambda CDM$ "ErisMod" set of simulations is used. These simulations can model both the stellar and dark matter components of the satellites, and their evolution under the action of the tides of a MW-sized host halo at a force resolution better than 10 pc. The stronger tidal mass loss and re-shaping of the mass distribution induced in satellites with $\gamma=0.6$ dark matter density distributions, as those resulting from the effect of feedback in hydrodynamical simulations of dwarf galaxy formation, is sufficient to bring the circular velocity profiles in agreement with the kinematics of MW's dSphs. In contrast, in simulations in which the satellites retain cusps at $z=0$ there are several "massive failures" with circular velocities in excess of the observational constraints. Various sources of deviations in the conventionally adopted relation between the circular velocity at the half light radius and the one dimensional line-of-sight velocity dispersions are found. Such deviations are caused by the response of circular velocity profiles to tidal effects, which also varies depending on the initially assumed inner density profile, and by the complexity of the stellar kinematics, which include residual rotation and anisotropy. In addition tidal effects naturally induce large deviations in the stellar mass-halo mass relation for halo masses below $\rm 10^9 ~ M_{\odot}$, preventing any reliable application of the abundance matching technique to dwarf galaxy satellites.

The near universality of DM halo density profiles provided by N-body simulations has proven to be robust against changes in total mass density, power spectrum, and some forms of initial velocity dispersion. In this letter we study the effects of coherently spinning up an isolated DM-only progenitor on halo structure. Halos with spins within several standard deviations of the simulated mean ($\lambda \lesssim 0.20$) produce profiles with negligible deviations from the universal form. Only when the spin becomes quite large ($\lambda \gtrsim 0.20$) do departures become evident. The angular momentum distribution also exhibits a near universal form, which is also independent of halo spin up to $\lambda \lesssim 0.20$. A correlation between these epidemic profiles and the presence of a strong bar in the virialized halo is also observed. These bar structures bear resemblance to the radial orbit instability in the rotationless limit.

Using cosmological galaxy formation simulations from the MaGICC project, spanning more than three magnitudes in stellar mass (~10^7-3x10^{10} Msun), we trace the baryonic cycle of infalling gas from the virial radius through to its participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and CGM. We derive the distributions of time between gas crossing the virial radius and being accreted to the star forming region (which allows mixing within the corona), as well as the time between gas being accreted to the star forming region and then forming stars (which allows mixing within the disc). Significant numbers of stars are formed from gas that cycles back through the hot halo after first accreting to the star forming region. Gas entering high mass galaxies is pre-enriched in low mass proto-galaxies prior to entering the virial radius of the central progenitor, with only small amounts of primordial gas accreted, even at high redshift (z~5). After entering the virial radius, significant further enrichment occurs prior to the accretion of the gas to the star forming region, with gas that is feeding the star forming region surpassing 0.1Z by z=0. Mixing with halo gas, itself enriched via galactic fountains, is thus crucial in determining the metallicity at which gas is accreted to the disc. The lowest mass simulation (Mvir~2x10^{10}Msun, with M*~10^7Msun), by contrast, accretes primordial gas through the virial radius and onto the disc at all times. Much like the classical analytical solutions to the `G-dwarf problem', overproduction of low-metallicity stars is ameliorated by the inefficiency of star formation. Finally, gas outflow/metal removal rates from star forming regions as a function of galactic mass are presented.

The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback.

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.

Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy mergers with two codes: the Smoothed Particle Hydrodynamics code Gasoline, and the Adaptive Mesh Refinement code Ramses. We seek to identify predictions of these models that are robust despite differences in hydrodynamic methods and implementations of sub-grid physics. We find that the general behavior is consistent between codes. Black hole accretion is minimal while the galaxies are well-separated (and even as they "fly-by" within 10 kpc at first pericenter). At late stages, when the galaxies pass within a few kpc, tidal torques drive nuclear gas inflow that triggers bursts of black hole accretion accompanied by star formation. We also note quantitative discrepancies that are model-dependent: our Ramses simulations show less star formation and black hole growth, and a smoother gas distribution with larger clumps and filaments, than our Gasoline simulations. We attribute these differences primarily to the sub-grid models for black hole fueling and feedback and gas thermodynamics. The main conclusion is that differences exist quantitatively between codes, and this should be kept in mind when making comparisons with observations. However, reassuringly, both codes capture the same dynamical behaviors in terms of triggering of black hole accretion, star formation, and black hole dynamics.

We show that the cool gas masses of galactic discs reach a steady state that lasts many Gyr after their last major merger in cosmological hydrodynamic simulations. The mass of disc gas, M$_{\rm gas}$, depends upon a galaxy halo's spin and virial mass, but not upon stellar feedback. Halos with low spin have high star formation efficiency and lower disc gas mass. Similarly, lower stellar feedback leads to more star formation so the gas mass ends up nearly the same irregardless of stellar feedback strength. Even considering spin, the M$_{\rm gas}$ relation with halo mass, M$_{200}$ only shows a factor of 3 scatter. The M$_{\rm gas}$--M$_{200}$ relation show a break at M$_{200}$=$2\times10^{11}$ M$_\odot$ that corresponds to an observed break in the M$_{\rm gas}$--M$_\star$ relation. The constant disc mass stems from a shared halo gas density profile in all the simulated galaxies. In their outer regions, the profiles are isothermal. Where the profile rises above $n=10^{-3}$ cm$^{-3}$, the gas readily cools and the profile steepens. Inside the disc, rotation supports gas with a flatter density profile except where supernova explosions disrupt the disc. Energy injection from stellar feedback also provides pressure support to the halo gas to prevent runaway cooling flows. The resulting constant gas mass makes simpler models for galaxy formation possible, either using a "bathtub" model for star formation rates or when modeling chemical evolution.

We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $\Lambda$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that naturally drives potential fluctuations able to create dark matter cores. Remarkably, SIDM fails to significantly lower the central dark matter density at halo peak velocities V$_{max}$ $<$ 30 Km/s. This is due to the fact that the central regions of very low--mass field halos have relatively low central velocity dispersion and densities, leading to time scales for SIDM collisions greater than a Hubble time. CDM halos with V$_{max}$ $<$ 30 km/s have inefficient star formation, and hence weak supernova feedback. At a fixed 2 cm2/g SIDM cross section, the DM content of very low mass CDM and SIDM halos differs by no more than a factor of two within 100-200pc. At larger halo masses ($\sim$ 10$^{10}$ solar masses), the introduction of baryonic processes creates field dwarf galaxies with dark matter cores and central DM$+$baryon distributions that are effectively indistinguishable between CDM and SIDM. Both models are in broad agreement with observed Local Group field galaxies across the range of masses explored. To significantly differentiate SIDM from CDM at the scale of faint dwarf galaxies, a velocity dependent cross section that rapidly increases to values larger than 2 $cm^{2}/g$ for halos with V$_{max}$ < 25-30 Km/s needs to be introduced.

We introduce a force correction term to better model the dynamical friction (DF) experienced by a supermassive black hole (SMBH) as it orbits within its host galaxy. This new approach accurately follows the orbital decay of a SMBH and drastically improves over commonly used advection methods. The force correction introduced here naturally scales with the force resolution of the simulation and converges as resolution is increased. In controlled experiments we show how the orbital decay of the SMBH closely follows analytical predictions when particle masses are significantly smaller than that of the SMBH. In a cosmological simulation of the assembly of a small galaxy, we show how our method allows for realistic black hole orbits. This approach overcomes the limitations of the advection scheme, where black holes are rapidly and artificially pushed toward the halo center and then forced to merge, regardless of their orbits. We find that SMBHs from merging dwarf galaxies can spend significant time away from the center of the remnant galaxy. Improving the modeling of SMBH orbital decay will help in making robust predictions of the growth, detectability, and merger rates of SMBHs, especially at low galaxy masses or at high redshift.

We present the first simulations of tidal stirring of dwarf galaxies in the Local Group carried out in a cosmological context. We use the ErisDARK simulation of a MW-sized galaxy to identify some of the most massive subhalos ($M_{vir} > 10^8 M_{\odot}$) that fall into the main host before $z=2$. Subhalos are replaced before infall with high-resolution models of dwarf galaxies comprising a faint stellar disk embedded in a dark matter halo. The set of models contains cuspy halos as well as halos with "cored" profiles (with asymptotic inner slope $\gamma = 0.6$). The simulations are then run to $z=0$ with as many as 54 million particles and resolution as small as $\sim 4$ pc using the N-Body code ChaNGa. The stellar components of all satellites are significantly affected by tidal stirring, losing stellar mass and undergoing a morphological transformation towards a pressure supported spheroidal system. However, while some remnants with cuspy halos maintain significant rotational flattening and disk-like features, all the shallow halo models achieve $v/\sigma < 0.5$ and round shapes typical of dSph satellites of the MW and M31. Mass loss is also enhanced in the latter, and remnants can reach luminosities and velocity dispersions as low as those of Ultra Faint Dwarfs (UFDs). We argue that cuspy progenitors must be the exception rather than the rule among satellites of the MW since all the MW and M31 satellites in the luminosity range of our remnants are dSphs, a result matched only in the simulation with "cored" models.

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}.

We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.

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 use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

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 Upsilon Andromedae system is the first exoplanetary system to have the relative inclination of two planets' orbital planes directly measured, and therefore offers our first window into the 3-dimensional configurations of planetary systems. We present, for the first time, full 3-dimensional, dynamically stable configurations for the 3 planets of the system consistent with all observational constraints. While the outer 2 planets, c and d, are inclined by about 30 degrees, the inner planet's orbital plane has not been detected. We use N-body simulations to search for stable 3-planet configurations that are consistent with the combined radial velocity and astrometric solution. We find that only 10 trials out of 1000 are robustly stable on 100 Myr timescales, or about 8 billion orbits of planet b. Planet b's orbit must lie near the invariable plane of planets c and d, but can be either prograde or retrograde. These solutions predict b's mass is in the range 2 - 9 $M_{Jup}$ and has an inclination angle from the sky plane of less than 25 degrees. Combined with brightness variations in the combined star/planet light curve ("phase curve"), our results imply that planet b's radius is about 1.8 $R_{Jup}$, relatively large for a planet of its age. However, the eccentricity of b in several of our stable solutions reaches values greater than 0.1, generating upwards of $10^{19}$ watts in the interior of the planet via tidal dissipation, possibly inflating the radius to an amount consistent with phase curve observations.

ChaNGa is an N-body cosmology simulation application implemented using Charm++. In this paper, we present the parallel design of ChaNGa and address many challenges arising due to the high dynamic ranges of clustered datasets. We focus on optimizations based on adaptive techniques for scaling to more than 128K cores. We demonstrate strong scaling on up to 512K cores of Blue Waters evolving 12 and 24 billion particles. We also show strong scaling of highly clustered datasets on up to 128K cores.

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 explore the impact of obliquity variations on planetary habitability in hypothetical systems with high mutual inclination. We show that large amplitude, high frequency obliquity oscillations on Earth-like exoplanets can suppress the ice-albedo feedback, increasing the outer edge of the habitable zone. We restrict our exploration to hypothetical systems consisting of a solar-mass star, an Earth-mass planet at 1 AU, and 1 or 2 larger planets. We verify that these systems are stable for $10^8$ years with N-body simulations, and calculate the obliquity variations induced by the orbital evolution of the Earth-mass planet and a torque from the host star. We run a simplified energy balance model on the terrestrial planet to assess surface temperature and ice coverage on the planet's surface, and we calculate differences in the outer edge of the habitable zone for planets with rapid obliquity variations. For each hypothetical system, we calculate the outer edge of habitability for two conditions: 1) the full evolution of the planetary spin and orbit, and 2) the eccentricity and obliquity fixed at their average values. We recover previous results that higher values of fixed obliquity and eccentricity expand the habitable zone, but also find that obliquity oscillations further expand habitable orbits in all cases. Terrestrial planets near the outer edge of the habitable zone may be more likely to support life in systems that induce rapid obliquity oscillations as opposed to fixed-spin planets. Such planets may be the easiest to directly characterize with space-borne telescopes.

We present a detailed comparison between the photometric properties of the bulges of two simulated galaxies and those of a uniform sample of observed galaxies. This analysis shows that the simulated galaxies have bulges with realistic surface brightnesses for their sizes and magnitude. These two field disc galaxies have rotational velocities ~ 100 km/s and were integrated to a redshift of zero in a fully cosmological Lambda cold dark matter context as part of high-resolution smoothed particle hydrodynamic simulations. We performed bulge-disc decompositions of the galaxies using artificial observations, in order to conduct a fair comparison to observations. We also dynamically decomposed the galaxies and compared the star formation histories of the bulges to those of the entire galaxies. These star formation histories showed that the bulges were primarily formed before z = 1 and during periods of rapid star formation. Both galaxies have large amounts of early star formation, which is likely related to the relatively high bulge-to-disc ratios also measured for them. Unlike almost all previous cosmological simulations, the realistically concentrated bulges of these galaxies do not lead to unphysically high rotational velocities, causing them to naturally lie along the observed Tully-Fisher relation.

We compare the central mass distribution of galaxies simulated with three different models of the interstellar medium (ISM) with increasing complexity: primordial (H+He) cooling down to 10^4K, additional cooling via metal lines and to lower temperatures, and molecular hydrogen (H_2) with shielding of atomic and molecular hydrogen, in addition to metal line cooling. In order to analyze the effect of these models, we follow the evolution of four field galaxies with V_peak < 120 km/s to a redshift of zero using high-resolution Smoothed Particle Hydrodynamic simulations in a fully cosmological LCDM context. The spiral galaxies produced in simulations with either primordial cooling or H_2 physics have realistic, rising rotation curves. In contrast, the simulations with metal line cooling and otherwise similar feedback and star formation produced galaxies with the peaked rotation curves typical of most previous LCDM simulations of spiral galaxies. The less-massive bulges and non-peaked rotation curves in the galaxies simulated with primordial cooling or H_2 are linked to changes in the angular momentum distribution of the baryons. These galaxies had smaller amounts of low-angular momentum baryons because of increased gas loss from stellar feedback. When there is only primordial cooling, the star forming gas is hotter and the feedback-heated gas cools more slowly than when metal line cooling is included and so requires less energy to be expelled. When H_2 is included, the accompanying shielding produces large amounts of clumpy, cold gas where H_2 forms. Star formation in clumpy gas results in more concentrated supernova feedback and greater efficiency of mass loss. The higher feedback efficiency causes a decrease of low-angular momentum material. (abridged)

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

Using a set of zoomed-in cosmological simulations of high-redshift progenitors of massive galaxies, we isolate and trace the history of gas that is accreted by central supermassive black holes. We determine the origins of the accreted gas, in terms of whether it entered the galaxy during a merger event or was smoothly accreted. Furthermore, we designate whether the smoothly accreted gas is accreted via a cold flow or is shocked upon entry into the halo. For moderate-mass (10^6 - 10^7 Msun) black holes at z ~ 4, there is a preference to accrete cold flow gas than gas of shocked or merger origin. However, this result is a consequence of the fact that the entire galaxy has a higher fraction of gas from cold flows. In general, each black hole tends to accrete the same fractions of smooth- and merger-accreted gas as is contained in its host galaxy, suggesting that once gas enters a halo it becomes well-mixed, and its origins are erased. We find that the angular momentum of the gas upon halo entry is a more important factor; black holes preferentially accrete gas that had low angular momentum when it entered the galaxy, regardless of whether it was accreted smoothly or through mergers.

We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly $10^{12}$ solar masses. Using a kinematic decomposition of these galaxies into present-day bulge and disk components, we find that the typical pressure of the star-forming ISM in the present-day bulge is higher than that in the present-day disk by an order of magnitude. We also find that pressure of the star-forming ISM at high redshift is on average, higher than ISM pressures at low redshift. This explains the why the bulge forms at higher pressures: the disk assembles at lower redshift, when the ISM is lower pressure and the bulge forms at high redshift, when the ISM is at higher pressure. If ISM pressure and IMF variation are tied together as suggested in studies like \cite{Conroy2012}, these results could indicate a time-dependent IMF in Milky-Way like systems, as well as a different IMF in the bulge and the disk.

We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar profiles. These effects are most evident for the 1 keV candidate but almost disappear for $m_{\mathrm{WDM}}>2$ keV. All simulations form similar stellar discs independent of WDM particle mass. In particular, the disc scale length does not change when WDM is considered. The reduced amount of star formation in the case of 1 keV particles is due to the effects of WDM on merging satellites which are on average less concentrated and less gas rich. The altered satellites cause a reduced starburst during mergers because they trigger weaker disc instabilities in the main galaxy. Nevertheless we show that disc galaxy evolution is much more sensitive to stellar feedback than it is to WDM candidate mass. Overall we find that WDM, especially when restricted to current observational constraints ($m_{\mathrm{WDM}}>2$ keV), has a minor impact on disc galaxy formation.

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation.

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.

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 study the influence of outer Solar System architecture on the structural evolution of the Oort Cloud (OC) and the flux of Earth-crossing comets. In particular, we seek to quantify the role of the giant planets as "planetary protectors". To do so, we have run simulations in each of four different planetary mass configurations to understand the significance of each of the giant planets. Because the outer planets modify the structure of the OC throughout its formation, we integrate each simulation over the full age of the Solar System. Over this time, we follow the evolution of cometary orbits from their starting point in the protoplanetary disk to their injection into the OC to their possible re-entry into the inner planetary region. We find that the overall structure of the OC, including the location of boundaries and the relative number of comets in the inner and outer parts, does not change significantly between configurations; however, as planetary mass decreases, the trapping efficiency (TE) of comets into the OC and the flux of comets into the observable region increases. We determine that those comets that evolve onto Earth-crossing orbits come primarily from the inner OC but show no preference for initial protoplanetary disk location. We also find that systems that have at least a Saturn-mass object are effective at deflecting possible Earth-crossing comets but the difference in flux between systems with and without such a planet is less than an order of magnitude. We conclude by discussing the individual roles of the planets and the implications of incorporating more realistic planetary accretion and migration scenarios into simulations, particularly on existing discrepancies between low TE and the mass of the protoplanetary disk and on determining the structural boundaries of the OC.

The advent of high-resolution imaging of galaxy clusters using the Sunyaev-Zel'dovich Effect (SZE) provides a unique probe of the astrophysics of the intracluster medium (ICM) out to high redshifts. To investigate the effects of cluster mergers on resolved SZE images, we present a high-resolution cosmological simulation of a 1.5E15 M_sun adiabatic cluster using the TreeSPH code ChaNGa. This massive cluster undergoes a 10:3:1 ratio triple merger accompanied by a dramatic rise in its integrated Compton-Y, peaking at z = 0.05. By modeling the thermal SZE (tSZ) and kinetic SZE (kSZ) spectral distortions of the Cosmic Microwave Background (CMB) at this redshift with relativistic corrections, we produce various mock images of the cluster at frequencies and resolutions achievable with current high-resolution SZE instruments. The two gravitationally-bound merging subclusters account for 10% and 1% of the main cluster's integrated Compton-Y, and have extended merger shock features in the background ICM visible in our mock images. We show that along certain projections and at specific frequencies, the kSZ CMB intensity distortion can dominate over the tSZ due to the large line of sight velocities of the subcluster gas and the unique frequency-dependence of these effects. We estimate that a one-velocity assumption in estimation of line of sight velocities of the merging subclusters from the kSZ induces a bias of ~10%. This velocity bias is small relative to other sources of uncertainty in observations, partially due to helpful bulk motions in the background ICM induced by the merger. Our results show that high-resolution SZE observations, which have recently detected strong kSZ signals in subclusters of merging systems, can robustly probe the dynamical as well as the thermal state of the ICM.

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

We introduce the Making Galaxies in a Cosmological Context (MaGICC) program of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter study of galaxy formation simulations of an L* galaxy that uses early stellar feedback combined with supernova feedback to match the stellar mass--halo mass relationship. While supernova feedback alone can reduce star formation enough to match the stellar mass--halo mass relationship, the galaxy forms too many stars before z=2 to match the evolution seen using abundance matching. Our early stellar feedback is purely thermal and thus operates like a UV ionization source as well as providing some additional pressure from the radiation of massive, young stars. The early feedback heats gas to >10^6 K before cooling to 10^4 K. The pressure from this hot gas creates a more extended disk and prevents more star formation prior to z=1 than supernovae feedback alone. The resulting disk galaxy has a flat rotation curve, an exponential surface brightness profile, and matches a wide range of disk scaling relationships. The disk forms from the inside-out with an increasing exponential scale length as the galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match observational results at low and high redshifts.

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Ly\alpha absorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.

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.

We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the polar disc. We also note that there is a twist in the shape profile, where the innermost 3 kpc of the halo flattens in the direction perpendicular to the old disc, and then aligns with the polar disc out until the virial radius. This result is then compared to the halo shape inferred from the circular velocities of the two discs. We also use the temporal information of the simulation to track the system's evolution, and identify the processes which give rise to this unusual galaxy type. We confirm the proposal that the polar disc galaxy is the result of the last major merger, where the angular moment of the interaction is orthogonal to the angle of the infalling gas. This merger is followed by the resumption of coherent gas infall. We emphasise that the disc is rapidly restored after the major merger and that after this event the galaxy begins to tilt. A significant proportion of the infalling gas comes from filaments. This infalling gas from the filament gives the gas its angular momentum, and, in the case of the polar disc galaxy, the direction of the gas filament does not change before or after the last major merger.

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).

Motivated by the observed connection between molecular hydrogen (H2) and star formation, we present a method for tracking the non-equilibrium abundance and cooling processes of H2 and H2-based star formation in Smoothed Particle Hydrodynamic simulations. The local abundances of H2 are calculated by integrating over the hydrogen chemical network. This calculation includes the gas-phase and dust grain formation of H2, shielding of HI and H2, and photodissociation of H2 by Lyman-Werner radiation from nearby stellar populations. Because this model does not assume equilibrium abundances, it is particularly well suited for simulations that model low-metallicity environments, such as dwarf galaxies and the early Universe. We further introduce an explicit link between star formation and local H2 abundance. This link limits star formation to "star-forming regions," represented by areas with abundant H2. With this implementation, we determine the effect of H2 on star formation in a cosmological simulation of a dwarf galaxy. This simulation is the first cosmological simulation with non-equilibrium H2 abundances to be integrated to a redshift of zero or to include efficient SN feedback. We find that our simulations are consistent with the observed Tully-Fisher, global Kennicutt-Schmidt, and resolved Kennicutt-Schmidt relations. We find that the inclusion of shielding of both the atomic and molecular hydrogen and, to a lesser extent, the additional cooling from H2 at temperatures between 200 and 5000 K increases the amount of cold gas in the galaxies. The changes to the ISM result in an increased amount of cold, dense gas in the disk of the galaxy and the formation of a clumpier interstellar media (ISM). The explicit link between star formation and H2 and the clumpier ISM results in a bluer galaxy with a greater spatial distribution of star formation at a redshift of zero. (abridged)

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

We use the same physical model to simulate four galaxies that match the relation between stellar and total mass, over a mass range that includes the vast majority of disc galaxies. The resultant galaxies, part of the Making Galaxies in a Cosmological Context (MaGICC) program, also match observed relations between luminosity, rotation velocity, size, colour, star formation rate, HI mass, baryonic mass, and metallicity. Radiation from massive stars and supernova energy regulate star formation and drive outflows, balancing the complex interplay between cooling gas, star formation, large scale outflows, and recycling of gas in a manner which correctly scales with the mass of the galaxy. Outflows also play a key role in simulating galaxies with exponential surface brightness profiles, flat rotation curves and dark matter cores. Our study implies that large scale outflows are the primary driver of the dependence of disc galaxy properties on mass. We show that the amount of outflows invoked in our model is required to meet the constraints provided by observations of OVI absorption lines in the circum-galactic-media of local galaxies.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u,g,r,i,z) photometry is used to determine the comets' colors, sizes, surface brightness profiles, and rates of dust production in terms of the Af{\rho} formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(< H) \propto 10(0.49\pm0.05)H for H < 18, with evidence of a much shallower fit N(< H) \propto 10(0.19\pm0.03)H for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 \pm 0.05 in g - r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.

We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies.

A clear prediction of the Cold Dark Matter model is the existence of cuspy dark matter halo density profiles on all mass scales. This is not in agreement with the observed rotation curves of spiral galaxies, challenging on small scales the otherwise successful CDM paradigm. In this work we employ high resolution cosmological hydro-dynamical simulations to study the effects of dissipative processes on the inner distribution of dark matter in Milky-Way like objects (M~1e12 Msun). Our simulations include supernova feedback, and the effects of the radiation pressure of massive stars before they explode as supernovae. The increased stellar feedback results in the expansion of the dark matter halo instead of contraction with respect to N-body simulations. Baryons are able to erase the dark matter cuspy distribution creating a flat, cored, dark matter density profile in the central several kpc of a massive Milky-Way like halo. The profile is well fit by a Burkert profile, with fitting parameters consistent with the observations. In addition, we obtain flat rotation curves as well as extended, exponential stellar disk profiles. While the stellar disk we obtain is still partially too thick to resemble the MW thin disk, this pilot study shows that there is enough energy available in the baryonic component to alter the dark matter distribution even in massive disc galaxies, providing a possible solution to the long standing problem of cusps vs. cores.

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).

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.

We investigate the formation of the first massive black holes in high redshift galaxies, with the goal of providing insights to which galaxies do or do not host massive black holes. We adopt a novel approach to forming seed black holes in galaxy halos in cosmological SPH+N-body simulations. The formation of massive black hole seeds is dictated directly by the local gas density, temperature, and metallicity, and motivated by physical models of massive black hole formation. We explore seed black hole populations as a function of halo mass and redshift, and examine how varying the efficiency of massive black hole seed formation affects the relationship between black holes and their hosts. Seed black holes tend to form in halos with mass between 10^7 and 10^9 Msun, and the formation rate is suppressed around z = 5 due to the diffusion of metals throughout the intergalactic medium. We find that the time of massive black hole formation and the occupation fraction of black holes are a function of the host halo mass. By z = 5, halos with mass M_halo > 3 x 10^9 Msun host massive black holes regardless of the efficiency of seed formation, while the occupation fraction for smaller halos increases with black hole formation efficiency. Our simulations explain why massive black holes are found in some bulgeless and dwarf galaxies, but we also predict that their occurrence becomes rarer and rarer in low-mass systems.

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.

We perform a suite of high-resolution smoothed particle hydrodynamics simulations to investigate the orbital decay and mass evolution of massive black hole (MBH) pairs down to scales of ~30 pc during minor mergers of disk galaxies. Our simulation set includes star formation and accretion onto the MBHs, as well as feedback from both processes. We consider 1:10 merger events starting at z~3, with MBH masses in the sensitivity window of the Laser Interferometer Space Antenna, and we follow the coupling between the merger dynamics and the evolution of the MBH mass ratio until the satellite galaxy is tidally disrupted. While the more massive MBH accretes in most cases as if the galaxy were in isolation, the satellite MBH may undergo distinct episodes of enhanced accretion, owing to strong tidal torques acting on its host galaxy and to orbital circularization inside the disk of the primary galaxy. As a consequence, the initial 1:10 mass ratio of the MBHs changes by the time the satellite is disrupted. Depending on the initial fraction of cold gas in the galactic disks and the geometry of the encounter, the mass ratios of the MBH pairs at the time of satellite disruption can stay unchanged or become as large as 1:2. Remarkably, the efficiency of MBH orbital decay correlates with the final mass ratio of the pair itself: MBH pairs that increase significantly their mass ratio are also expected to inspiral more promptly down to nuclear-scale separations. These findings indicate that the mass ratios of MBH pairs in galactic nuclei do not necessarily trace the mass ratios of their merging host galaxies, but are determined by the complex interplay between gas accretion and merger dynamics.

We investigate the evolution of grains composed of an ice shell surrounding an olivine core as they pass through a spiral shock in a protoplanetary disk. We use published three-dimensional radiation-hydrodynamics simulations of massive self-gravitating protoplanetary disks to extract the thermodynamics of spiral shocks in the region between 10 and 20 AU from the central star. As the density wave passes, it heats the grains, causing them to lose their ice shell and resulting in a lowering of the grain opacity. In addition, since grains of different sizes will have slightly different temperatures, there will be a migration of ice from the hotter grains to the cooler ones. The rate of migration depends on the temperature of the background gas, so a grain distribution that is effectively stable for low temperatures, can undergo an irreversible change in opacity if the gas is temporarily heated to above $\sim 150$\,K. We find that the opacity can drop more, and at a significantly faster rate throughout the spiral shocks relative to the prediction of standard dust grains models adopted in hydrodynamical calculations of protoplanetary disks. This would lead to faster gas cooling within spiral arms. We discuss the implications of our results on the susceptibility of disk fragmentation into sub-stellar objects at distances of a few tens of astronomical units.

We present observations and dynamical models of the stellar nuclear clusters (NCs) at the centres of NGC 4244 and M33. We then compare these to an extensive set of simulations testing the importance of purely stellar dynamical mergers on the formation and growth of NCs. Mergers of star clusters are able to produce a wide variety of observed properties, including densities, structural scaling relations, shapes (including the presence of young discs) and even rapid rotation. Nonetheless, difficulties remain, most notably that the second order kinematic moment V_rms = (V^2 + sigma^2)^(1/2) of the models is too centrally peaked to match observations. This can be remedied by the merger of star clusters onto a pre-existing nuclear disc, but the line-of-sight velocity V is still more slowly rising than in NGC 4244. Our results therefore suggest that purely stellar dynamical mergers cannot form NCs, and that gas dissipation is a necessary ingredient for at least ~50% of a NC's mass. The negative vertical anisotropy found in NGC 4244 however requires at least 10% of the mass to be accreted as stars, since gas dissipation and in situ star formation leads to positive vertical anisotropy.

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.

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 perform SPH+N-body cosmological simulations of massive disk galaxies, including a formalism for black hole seed formation and growth, and find that satellite galaxies containing supermassive black hole seeds are often stripped as they merge with the primary galaxy. These events naturally create a population of "wandering" black holes that are the remnants of stripped satellite cores; galaxies like the Milky Way may host 5 -- 15 of these objects within their halos. The satellites that harbor black hole seeds are comparable to Local Group dwarf galaxies such as the Small and Large Magellanic Clouds; these galaxies are promising candidates to host nearby intermediate mass black holes. Provided that these wandering black holes retain a gaseous accretion disk from their host dwarf galaxy, they give a physical explanation for the origin and observed properties of some recently discovered off-nuclear ultraluminous X-ray sources such as HLX-1.

It has previously been shown that varying the numerical timestep during a symplectic orbital integration leads to a random walk in energy and angular momentum, destroying the phase space-conserving property of symplectic integrators. Here we show that when altering the timestep symplectic correctors can be used to reduce this error to a negligible level. Furthermore, these correctors can also be employed to avoid a large error introduction when changing the Hamiltonian's partitioning. We have constructed a numerical integrator using this technique that is nearly as accurate as widely used fixed-step routines. In addition, our algorithm is drastically faster for integrations of highly eccentricitic, large semimajor axis orbits, such as those found in the Oort Cloud.

We evaluate the orbital evolution and several plausible origins scenarios for the mutually inclined orbits of Upsilon Andromedae c and d. These two planets have orbital elements that oscillate with large amplitudes and lie close to the stability boundary. This configuration, and in particular the observed mutual inclination, demands an explanation. The planetary system may be influenced by a nearby low-mass star, Upsilon And B, which could perturb the planetary orbits, but we find it cannot modify two coplanar orbits into the observed mutual inclination of ~30 deg. However, it could incite ejections or collisions between planetary companions that subsequently raise the mutual inclination to >30 deg. Our simulated systems with large mutual inclinations tend to be further from the stability boundary than Upsilon And, but we are able to produce similar systems. We conclude that scattering is a plausible mechanism to explain the observed orbits of Upsilon And c and d, but we cannot determine whether the scattering was caused by instabilities among the planets themselves or by perturbations from Upsilon And B. We also develop a procedure to quantitatively compare numerous properties of the observed system to our numerical models. Although we only implement this procedure to Upsilon And, it may be applied to any exoplanetary system.

A sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the Lambda CDM cosmology. Artificial images in the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface brightness profiles as observations would, and avoiding any assumption that light must follow mass as previous models have assumed. We demonstrate that these simulated disks are an excellent match to the observed magnitude - size relation for both local disks, and for disks at z=1 in the magnitude/mass range of overlap. We disentangle the evolution seen in the population as a whole from the evolution of individual disk galaxies. In agreement with observations, our simulated disks undergo roughly 1.5 magnitudes/arcsec^2 of surface brightness dimming since z=1. We find evidence that evolution in the magnitude - size plane varies by mass, such that galaxies with M* > 10^9 M_sun undergo more evolution in size than luminosity, while dwarf galaxies tend to evolve potentially more in luminosity. The disks grow in such a way as to stay on roughly the same stellar mass - size relation with time. Finally, due to an evolving stellar mass - SFR relation, a galaxy at a given stellar mass (or size) at z=1 will reside in a more massive halo and have a higher SFR, and thus a higher luminosity, than a counterpart of the same stellar mass at z=0.

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

We examine the effect of mass and force resolution on a specific star formation (SF) recipe using a set of N-body/Smooth Particle Hydrodynamic simulations of isolated galaxies. Our simulations span halo masses from 10^9 to 10^13 solar masses, more than four orders of magnitude in mass resolution, and two orders of magnitude in the gravitational softening length, epsilon, representing the force resolution. We examine the total global star formation rate, the star formation history, and the quantity of stellar feedback and compare the disk structure of the galaxies. Based on our analysis, we recommend using at least 10^4 particles each for the dark matter and gas component and a force resolution of epsilon approximately equal to 10^-3 R_vir when studying global SF and feedback. When the spatial distribution of stars is important, the number of gas and dark matter particles must be increased to at least 10^5 of each. Low mass resolution simulations with fixed softening lengths show particularly weak stellar disks due to two-body heating. While decreasing spatial resolution in low mass resolution simulations limits two-body effects, density and potential gradients cannot be sustained. Regardless of the softening, low-mass resolution simulations contain fewer high density regions where SF may occur. Galaxies of approximately 10^10 solar masses display unique sensitivity to both mass and force resolution. This mass of galaxy has a shallow potential and is on the verge of forming a disk. The combination of these factors give this galaxy the potential for strong gas outflows driven by supernova feedback and make it particularly sensitive to any changes to the simulation parameters.

We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$\times10^{11}$ M$_\odot$ to 2$\times10^{12}$ M$_\odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.

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.

We summarize recent developments in the use of spectral methods for analyzing large numbers of orbits in N-body simulations to obtain insights into the global phase space structure of dark matter halos. The fundamental frequencies of oscillation of orbits can be used to understand the physical mechanism by which the shapes of dark matter halos evolve in response to the growth of central baryonic components. Halos change shape primarily because individual orbits change their shapes adiabatically in response to the growth of a baryonic component, with those at small radii become preferentially rounder. Chaotic scattering of orbits occurs only when the central point mass is very compact and is equally effective for centrophobic long-axis tube orbits as it is for centrophilic box orbits.

Using spectral methods, we analyse the orbital structure of prolate/triaxial dark matter (DM) halos in N-body simulations to understand the processes that drive the evolution of shapes of DM halos and elliptical galaxies in which central masses are grown. A longstanding issue is whether the change in the shapes of DM halos is the result of chaotic scattering of box orbits, or whether they change shape adiabatically in response to the evolving galactic potential. We use orbital frequencies to classify orbits, to quantify orbital shapes, and to identify resonant orbits and chaotic orbits. The frequency-based method overcomes the limitations of Lyapunov exponents which are sensitive to numerical discreteness effects. Regardless of the distribution of the baryonic component, the shape of a DM halo changes primarily due to changes in the shapes of individual orbits within a given family. Orbits with small pericentric radii are more likely to change both their orbital type and shape than orbits with large pericentric radii. Whether the evolution is regular (and reversible) or chaotic (and irreversible), depends primarily on the radial distribution of the baryonic component. The growth of an extended baryonic component of any shape results in a regular rather than chaotic change in orbital populations. In contrast the growth of a massive and compact central component results in chaotic scattering of a significant fraction of both box and long-axis tube orbits. The growth of a disk causes a significant fraction of halo particles to become trapped by major global orbital resonances. Despite the fact that shape of a DM halo is always quite oblate following the growth of a central baryonic component, a significant fraction of its orbit population has characteristics of its triaxial or prolate progenitor (ABRIDGED).

Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This method is implemented in the parallel N-body code, PKDGRAV and tested on some simple orbits. The method demonstrates a lack of secular changes in orbital elements, making it a very useful technique for integrating Hill's equations over many dynamical times. Furthermore, the method allows for efficient collision searching using linear extrapolation of particle positions.

We present numerical simulations to model the production of observable long-period comets (LPCs) from the Oort Cloud, a vast reservoir of icy bodies surrounding the Sun. We show that inner Oort Cloud objects can penetrate Jupiter's orbit via a largely unexplored dynamical pathway, and they are an important, if not the dominant, source of known LPCs. We use this LPC production to place observationally motivated constraints on the population and mass of the inner Oort Cloud, which are consistent with giant planet formation theory. These constraints indicate that only one comet shower producing late Eocene bombardment levels has likely occurred since the Cambrian Explosion, making these phenomena an improbable cause of additional extinction events.

In recent work, we have shown that it is possible to link quantitatively many aspects of damped Lyman alpha (DLA) absorbers in the spectra of quasars to high resolution simulations of galaxy formation. Using runs from the same series of hydrodynamic numerical studies, we consider the expected properties of Lyman alpha absorbers seen in the spectra of high redshift (z>2) gamma ray burst afterglows (GRB-DLAs). If GRBs are associated with the death of massive stars, their afterglows provide insights into otherwise unprobed regions of protogalactic objects, but detailed physical interpretations are currently embryonic. We find that median impact parameters (measured from the potential minimum) are approximately 1 kpc for GRBs compared with 4 kpc for QSO-DLAs. However, an equally important difference is that GRB-DLAs are predominantly associated with halos of mass 10^10<M_vir/M_sol<10^12, an order of magnitude larger than the hosts of QSO-DLAs. Accordingly, there are differences in the stellar properties of hosts. Our simulations accurately predict the form of the GRB-DLA HI column density distribution, producing quantitative agreement for N_HI>10^19 cm^-2, but they somewhat underpredict the incidence of low column densities N_HI<10^19 cm^-2. Line-of-sight neutral gas metallicities predicted by our simulations (10^-2 < Z/Z_sol < 1) are consistent with the modest observational constraints. Because of large internal dispersions in gas metallicities, this agreement is not significantly compromised by imposing a cut-off on the metallicity of stars able to launch GRBs (Z_star<Z_sol/3), confounding claims that the observed metallicity of GRB-DLAs poses a challenge to current GRB models. (Abridged.)

We trace the evolution of central galaxies in three ~10^13 M_sun galaxy groups simulated at high resolution in cosmological hydrodynamical simulations. The evolution in the group potential leads, at z=0, to central galaxies that are massive, gas-poor early-type systems supported by stellar velocity dispersion resembling either elliptical or S0 galaxies. Their z~2-2.5 main progenitors are massive M* ~ 3-10 x 10^10 M_sun, star forming (20-60 M_sun/yr) galaxies which host substantial reservoirs of cold gas (~5 x 10^9 M_sun) in extended gas disks. Our simulations thus show that star forming galaxies observed at z~2 are likely the main progenitors of central galaxies in galaxy groups at z=0. Their central stellar densities stay approximately constant from z~1.5 down to z=0. Instead, the galaxies grow inside-out, by acquiring a stellar envelope outside the innermost ~2 kpc. Consequently the density within the effective radius decreases by up to two orders of magnitude. Both major and minor mergers contribute to most of the mass accreted outside the effective radius and thus drive the evolution of the half-mass radii. In one of the three simulated groups the short central cooling time leads to a dramatic rejuvenation of the central group galaxy at z<1, affecting its morphology, kinematics and colors. This episode is eventually terminated by a group-group merger. Our simulations demonstrate that, in galaxy groups, the interplay between halo mass assembly, galaxy merging and gas accretion has a substantial influence on the star formation histories and z=0 morphologies of central galaxies.[Abridged]

For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario.

We investigate the formation of the stellar halos of four simulated disk galaxies using high resolution, cosmological SPH + N-Body simulations. These simulations include a self-consistent treatment of all the major physical processes involved in galaxy formation. The simulated galaxies presented here each have a total mass of ~10^12 M_sun, but span a range of merger histories. These simulations allow us to study the competing importance of in-situ star formation (stars formed in the primary galaxy) and accretion of stars from subhalos in the building of stellar halos in a LambdaCDM universe. All four simulated galaxies are surrounded by a stellar halo, whose inner regions (r < 20 kpc) contain both accreted stars, and an in-situ stellar population. The outer regions of the galaxies' halos were assembled through pure accretion and disruption of satellites. Most of the in-situ halo stars formed at high redshift out of smoothly accreted cold gas in the inner 1 kpc of the galaxies' potential wells, possibly as part of their primordial disks. These stars were displaced from their central locations into the halos through a succession of major mergers. We find that the two galaxies with recently quiescent merger histories have a higher fraction of in-situ stars (~20-50%) in their inner halos than the two galaxies with many recent mergers (~5-10% in-situ fraction). Observational studies concentrating on stellar populations in the inner halo of the Milky Way will be the most affected by the presence of in-situ stars with halo kinematics, as we find that their existence in the inner few tens of kpc is a generic feature of galaxy formation.

Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z.

We review the problem of the formation of terrestrial planets, with particular emphasis on the interaction of dynamical and geochemical models. The lifetime of gas around stars in the process of formation is limited to a few million years based on astronomical observations, while isotopic dating of meteorites and the Earth-Moon system suggest that perhaps 50-100 million years were required for the assembly of the Earth. Therefore, much of the growth of the terrestrial planets in our own system is presumed to have taken place under largely gas-free conditions, and the physics of terrestrial planet formation is dominated by gravitational interactions and collisions. The earliest phase of terrestrial-planet formation involve the growth of km-sized or larger planetesimals from dust grains, followed by the accumulations of these planetesimals into ~100 lunar- to Mars-mass bodies that are initially gravitationally isolated from one-another in a swarm of smaller planetesimals, but eventually grow to the point of significantly perturbing one-another. The mutual perturbations between the embryos, combined with gravitational stirring by Jupiter, lead to orbital crossings and collisions that drive the growth to Earth-sized planets on a timescale of 10-100 million years. Numerical treatment of this process has focussed on the use of symplectic integrators which can rapidly integrate the thousands of gravitationally-interacting bodies necessary to accurately model planetary growth. While the general nature of the terrestrial planets--their sizes and orbital parameters--seem to be broadly reproduced by the models, there are still some outstanding dynamical issues. One of these is the presence of an embryo-sized body, Mars, in our system in place of the more massive objects that simulations tend to yield. [Abridged]

The ACS Nearby Galaxy Survey Treasury (ANGST) is a systematic survey to establish a legacy of uniform multi-color photometry of resolved stars for a volume-limited sample of nearby galaxies (D<4 Mpc). The survey volume encompasses 69 galaxies in diverse environments, including close pairs, small & large groups, filaments, and truly isolated regions. The galaxies include a nearly complete range of morphological types spanning a factor of ~10^4 in luminosity and star formation rate. The survey data consists of images taken with ACS on HST, supplemented with archival data and new WFPC2 imaging taken after the failure of ACS. Survey images include wide field tilings covering the full radial extent of each galaxy, and single deep pointings in uncrowded regions of the most massive galaxies in the volume. The new wide field imaging in ANGST reaches median 50% completenesses of m_F475W=28.0 mag, m_F606W=27.3 mag, and m_F814W=27.3 mag, several magnitudes below the tip of the red giant branch (TRGB). The deep fields reach magnitudes sufficient to fully resolve the structure in the red clump. The resulting photometric catalogs are publicly accessible and contain over 34 million photometric measurements of >14 million stars. In this paper we present the details of the sample selection, imaging, data reduction, and the resulting photometric catalogs, along with an analysis of the photometric uncertainties (systematic and random), for both the ACS and WFPC2 imaging. We also present uniformly derived relative distances measured from the apparent magnitude of the TRGB.

We examine the pairing process of supermassive black holes (SMBHs) down to scales of 20-100 pc using a set of N-body/SPH simulations of binary mergers of disk galaxies with mass ratios of 1:4 and 1:10. Our numerical experiments are designed to represent merger events occurring at various cosmic epochs. The initial conditions of the encounters are consistent with the LambdaCDM paradigm of structure formation, and the simulations include the effects of radiative cooling, star formation, and supernovae feedback. We find that the pairing of SMBHs depends sensitively on the amount of baryonic mass preserved in the center of the companion galaxies during the last phases of the merger. In particular, due to the combination of gasdynamics and star formation, we find that a pair of SMBHs can form in 1:10 minor mergers provided that galaxies are relatively gas-rich (gas fractions of 30% of the disk mass) and that the mergers occur at relatively high redshift (z~3), when dynamical friction timescales are shorter. Since 1:10 mergers are most common events during the assembly of galaxies, and mergers are more frequent at high redshift when galaxies are also more gas-rich, our results have positive implications for future gravitational wave experiments such as the Laser Interferometer Space Antenna.

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 solar mass star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10^5 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.

Stellar population studies show that low mass galaxies in all environments exhibit stellar halos that are older and more spherically distributed than the main body of the galaxy. In some cases, there is a significant intermediate age component that extends beyond the young disk. We examine a suite of Smoothed Particle Hydrodynamic (SPH) simulations and find that elevated early star formation activity combined with supernova feedback can produce an extended stellar distribution that resembles these halos for model galaxies ranging from $v_{200}$ = 15 km s$^{-1}$ to 35 km s$^{-1}$, without the need for accretion of subhalos.

We report the discovery of a minor planet (2006 SQ372) on an orbit with a perihelion of 24 AU and a semimajor axis of 796 AU. Dynamical simulations show that this is a transient orbit and is unstable on a timescale of 200 Myrs. Falling near the upper semimajor axis range of the scattered disk and the lower semimajor axis range of the Oort Cloud, previous membership in either class is possible. By modeling the production of similar orbits from the Oort Cloud as well as from the scattered disk, we find that the Oort Cloud produces 16 times as many objects on SQ372-like orbits as the scattered disk. Given this result, we believe this to be the most distant long-period comet ever discovered. Furthermore, our simulation results also indicate that 2000 OO67 has had a similar dynamical history. Unaffected by the "Jupiter-Saturn Barrier," these two objects are most likely long-period comets from the inner Oort Cloud.

We use high resolution cosmological hydrodynamical simulations to demonstrate that cold flow gas accretion, particularly along filaments, modifies the standard picture of gas accretion and cooling onto galaxy disks. In the standard picture, all gas is initially heated to the virial temperature of the galaxy as it enters the virial radius. Low mass galaxies are instead dominated by accretion of gas that stays well below the virial temperature, and even when a hot halo is able to develop in more massive galaxies there exist dense filaments that penetrate inside of the virial radius and deliver cold gas to the central galaxy. For galaxies up to ~L*, this cold accretion gas is responsible for the star formation in the disk at all times to the present. Even for galaxies at higher masses, cold flows dominate the growth of the disk at early times. Within this modified picture, galaxies are able to accrete a large mass of cold gas, with lower initial gas temperatures leading to shorter cooling times to reach the disk. Although star formation in the disk is mitigated by supernovae feedback, the short cooling times allow for the growth of stellar disks at higher redshifts than predicted by the standard model.

We investigate the population of z=3 damped Lyman alpha absorbers (DLAs) in a recent series of high resolution galaxy formation simulations. The simulations are of interest because they form at z=0 some of the most realistic disk galaxies to date. No free parameters are available in our study: these have been fixed by physical and z=0 observational constraints, and thus our study provides a genuine consistency test. The precise role of DLAs in galaxy formation remains in debate, but they provide a number of strong constraints on the nature of our simulated bound systems at z=3 because of their coupled information on neutral H I densities, kinematics, metallicity and estimates of star formation activity. Our results, without any parameter-tuning, closely match the observed incidence rate and column density distributions of DLAs. Our simulations are the first to reproduce the distribution of metallicities (with a median of Z_{DLA} = Z_{solar}/20) without invoking observationally unsupported mechanisms such as dust biasing. This is especially encouraging given that these simulations have previously been shown to have a realistic 0<z<2 stellar mass-metallicity relation. Additionally, we see a strong positive correlation between sightline metallicity and low-ion velocity width, the normalization and slope of which comes close to matching recent observational results. However, we somewhat underestimate the number of observed high velocity width systems; the severity of this disagreement is comparable to other recent DLA focused studies. By z=0 the majority of the z=3 neutral gas forming the DLAs has been converted into stars, in agreement with rough physical expectations. [Abridged]

The ACS Nearby Galaxy Survey Treasury (ANGST) is a large Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) treasury program to obtain resolved stellar photometry for a volume-limited sample of galaxies out to 4 Mpc. As part of this program, we have obtained deep ACS imaging of a field in the outer disk of the large spiral galaxy M81. The field contains the outskirts of a spiral arm as well as an area containing no current star formation. Our imaging results in a color-magnitude diagram (CMD) reaching to F814W = 28.8 and F606W = 29.5, one magnitude fainter than the red clump. Through detailed modeling of the full CMD, we quantify the age and metallicity distribution of the stellar populations contained in the field. The mean metallicity in the field is -1<[M/H]<0 and only a small fraction of stars have ages <~1 Gyr. The results show that most of the stars in this outer disk field were formed by z~1 and that the arm structure at this radius has a lifetime of >~100 Myr. We discuss the measured evolution of the M81 disk in the context of surveys of high-redshift disk galaxies and deep stellar photometry of other nearby galaxies. All of these indicate that massive spiral disks are mostly formed by z~1 and that they have experienced rapid metal enrichment.

We study dynamics of bars in models of disk galaxies embeded in realistic dark matter halos. We find that disk thickness plays an important, if not dominant, role in the evolution and structure of the bars. We also make extensive numerical tests of different N-body codes used to study bar dynamics. Models with thick disks typically used in this type of modeling (height-to-length ratio hz/Rd=0.2) produce slowly rotating, and very long, bars. In contrast, more realistic thin disks with the same parameters as in our Galaxy (hz/Rd= 0.1) produce bars with normal length Rbar approx R_d, which rotate quickly with the ratio of the corotation radius to the bar radius 1.2-1.4 compatible with observations. Bars in these models do not show a tendency to slow down, and may lose as little as 2-3 percent of their angular momentum due to dynamical friction with the dark matter over cosmological time. We attribute the differences between the models to a combined effect of high phase-space density and smaller Jeans mass in the thin disk models, which result in the formation of a dense central bulge. Special attention is paid to numerical effects such as the accuracy of orbital integration, force and mass resolution. Using three N-body codes -- Gadget, ART, and Pkdgrav -- we find that numerical effects are very important and, if not carefully treated, may produce incorrect and misleading results. Once the simulations are performed with sufficiently small time-steps and with adequate force and mass resolution, all the codes produce nearly the same results: we do not find any systematic deviations between the results obtained with TREE codes (Gadget and Pkdgrav) and with the Adaptive-Mesh-Refinement (ART) code.

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.

For a decade, N-body simulations have revealed a nearly universal dark matter density profile, which appears to be robust to changes in the overall density of the universe and the underlying power spectrum. Despite its universality, the physical origin of this profile has not yet been well understood. Semi--analytic models by Barnes et al. (2005) have suggested that the density structure of dark matter halos is determined by the onset of the radial orbit instability (ROI). We have tested this hypothesis using N-body simulations of collapsing dark matter halos with a variety of initial conditions. For dynamically cold initial conditions, the resulting halo structures are triaxial in shape, due to the mild aspect of the instability. We examine how variations in initial velocity dispersion affect the onset of the instability, and find that an isotropic velocity dispersion can suppress the ROI entirely, while a purely radial dispersion does not. The quantity sigma^2/vc^2 is a criterion for instability, where regions with sigma^2/vc^2 <~1 become triaxial due to the ROI or other perturbations. We also find that the radial orbit instability sets a scale length at which the velocity dispersion changes rapidly from isotropic to radially anisotropic. This scale length is proportional to the radius at which the density profile changes shape, as is the case in the semi--analytic models; however, the coefficient of proportionality is different by a factor of ~2.5. We conclude that the radial orbit instability is likely to be a key physical mechanism responsible for the nearly universal profiles of simulated dark matter halos.

Cold dark matter cosmogony predicts triaxial dark matter halos, whereas observations find quite round halos. This is most likely due to the condensation of baryons leading to rounder halos. We examine the halo phase space distribution basis for such shape changes. Triaxial halos are supported by box orbits, which pass arbitrarily close to the density center. The decrease in triaxiality caused by baryons is thought to be due to the scattering of these orbits. We test this hypothesis with simulations of disks grown inside triaxial halos. After the disks are grown we check whether the phase space structure has changed by evaporating the disks and comparing the initial and final states. While the halos are substantially rounder when the disk is at full mass, their final shape after the disk is evaporated is not much different from the initial. Likewise, the halo becomes (more) radially anisotropic when the disk is grown, but the final anisotropy is consistent with the initial. Only if the baryons are unreasonably compact or massive does the halo change irreversibly. We show that the character of individual orbits is not generally changed by the growing mass. Thus the central condensation of baryons does not destroy enough box orbits to cause the shape change. Rather, box orbits merely become rounder along with the global potential. However, if angular momentum is transferred to the halo, either via satellites or via bars, a large irreversible change in the halo distribution occurs. The ability of satellites to alter the phase space distribution of the halo is of particular concern to galaxy formation simulations since halo triaxiality can profoundly influence the evolution of disks.

Using effective temperature and metallicity derived from SDSS spectra for ~60,000 F and G type main sequence stars (0.2<g-r<0.6), we develop polynomial models for estimating these parameters from the SDSS u-g and g-r colors. We apply this method to SDSS photometric data for about 2 million F/G stars and measure the unbiased metallicity distribution for a complete volume-limited sample of stars at distances between 500 pc and 8 kpc. The metallicity distribution can be exquisitely modeled using two components with a spatially varying number ratio, that correspond to disk and halo. The two components also possess the kinematics expected for disk and halo stars. The metallicity of the halo component is spatially invariant, while the median disk metallicity smoothly decreases with distance from the Galactic plane from -0.6 at 500 pc to -0.8 beyond several kpc. The absence of a correlation between metallicity and kinematics for disk stars is in a conflict with the traditional decomposition in terms of thin and thick disks. We detect coherent substructures in the kinematics--metallicity space, such as the Monoceros stream, which rotates faster than the LSR, and has a median metallicity of [Fe/H]=-0.96, with an rms scatter of only ~0.15 dex. We extrapolate our results to the performance expected from the Large Synoptic Survey Telescope (LSST) and estimate that the LSST will obtain metallicity measurements accurate to 0.2 dex or better, with proper motion measurements accurate to ~0.2 mas/yr, for about 200 million F/G dwarf stars within a distance limit of ~100 kpc (g<23.5). [abridged]

We study the influence of an open cluster environment on the formation and current structure of the Oort cloud. To do this, we have run 19 different simulations of the formation of the Oort Cloud for 4.5 Gyrs. In each simulation, the solar system spends its first 100 Myrs in a different open cluster environment before transitioning to its current field environment. We find that, compared to forming in the field environment, the inner Oort Cloud is preferentially loaded with comets while the Sun resides in the open cluster and that most of this material remains locked in the interior of the cloud for the next 4.4 Gyrs. In addition, the outer Oort Cloud trapping efficiencies we observe in our simulations are lower than previous formation models by about a factor of 2, possibly implying an even more massive early planetesimal disk. Furthermore, some of our simulations reproduce the orbits of observed extended scattered disk objects, which may serve as an observational constraint on the Sun's early environment. Depending on the particular open cluster environment, the properties of the inner Oort Cloud and extended scattered disk can vary widely. On the other hand, the outer portions of the Oort Cloud in each of our simulations are all similar.

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.

Polar Ring Galaxies, such as NGC4650A, are a class of galaxy which have two kinematically distinct components that are inclined by almost 90 degrees to each other. These striking galaxies challenge our understanding of how galaxies form; the origin of their distinct components has remained uncertain, and the subject of much debate. We use high-resolution cosmological simulations of galaxy formation to show that Polar Ring Galaxies are simply an extreme example of the angular moment misalignment that occurs during the hierarchical structure formation characteristic of Cold Dark Matter cosmology. In our model, Polar Ring Galaxies form through the continuous accretion of gas whose angular momentum is misaligned with the central galaxy.

Supermassive black holes (SMBHs) are a ubiquitous component of the nuclei of galaxies. It is normally assumed that, following the merger of two massive galaxies, a SMBH binary will form, shrink due to stellar or gas dynamical processes and ultimately coalesce by emitting a burst of gravitational waves. However, so far it has not been possible to show how two SMBHs bind during a galaxy merger with gas due to the difficulty of modeling a wide range of spatial scales. Here we report hydrodynamical simulations that track the formation of a SMBH binary down to scales of a few light years following the collision between two spiral galaxies. A massive, turbulent nuclear gaseous disk arises as a result of the galaxy merger. The black holes form an eccentric binary in the disk in less than a million years as a result of the gravitational drag from the gas rather than from the stars.

We simulate the collapse of isolated dwarf galaxies using SPH + N-Body simulations including a physically motivated description of the effects of supernova feedback. As the gas collapses and stars form, the supernova feedback disrupts enough gas to temporarily quench star formation. The gas flows outward into a hot halo, where it cools until star formation can continue once more and the cycle repeats. The star formation histories of isolated Local Group dwarf galaxies exhibit similar episodic bursts of star formation. We examine the mass dependence of the stellar velocity dispersions and find that they are no less than half the velocity of the halos measured at the virial radius.

We construct mock galaxy catalogues to analyse clustering properties of a Lambda cold dark matter (LCDM) universe within a cosmological dark matter simulation of sufficient resolution to resolve structure down to the scale of dwarfs. We show that there is a strong age-clustering correlation for objects likely to host luminous galaxies, which includes the satellite halo (subhalo) population. Older mock galaxies are significantly more clustered in our catalog, which consists of satellite haloes as well as the central peaks of discrete haloes, selected solely by peak circular velocity. This age dependence is caused mainly by the age-clustering relation for discrete haloes, recently found by Gao et al., acting mostly on field members, combined with the tendency for older mock galaxies to lie within groups and clusters, where galaxy clustering is enhanced. Our results suggest that the clustering age dependence is manifested in real galaxies. At small scales (less than ~5 Mpc/h), the very simple assumption that galaxy colour depends solely on halo age is inconsistent with the strength of the observed clustering colour trends, where red galaxies become increasingly more clustered than blue galaxies toward smaller scales, suggesting that luminosity weighted galaxy ages do not closely trace the assembly epoch of their dark matter hosts. The age dependence is present but is much weaker for satellite haloes lying within groups and clusters than for the global population.

We model the formation, evolution and astrophysical effects of dark compact Scalar Miniclusters (``ScaMs''). These objects arise when a scalar field, with an axion-like or Higgs-like potential, undergoes a second order phase transition below the QCD scale. Such a scalar field may couple too weakly to the standard model to be detectable directly through particle interactions, but may still be detectable by gravitational effects, such as lensing and baryon accretion by large, gravitationally bound miniclusters. The masses of these objects are shown to be constrained by the Ly$\alpha$ power spectrum to be less than $\sim 10^4 M_\odot$, but they may be as light as classical axion miniclusters, of the order of $10^{-12} M_\odot$. We simulate the formation and nonlinear gravitational collapse of these objects around matter-radiation equality using an N-body code, estimate their gravitational lensing properties, and assess the feasibility of studying them using current and future lensing experiments. Future MACHO-type variability surveys of many background sources can reveal either high-amplification, strong lensing events, or measure density profiles directly via weak-lensing variability, depending on ScaM parameters and survey depth. However, ScaMs, due to their low internal densities, are unlikely to be responsible for apparent MACHO events already detected in the Galactic halo. A simple estimate is made of parameters that would give rise to early structure formation; in principle, early stellar collapse could be triggered by ScaMs as early as recombination, and significantly affect cosmic reionization.

We use the suite of simulations presented in Poole et al (2006) to examine global X-ray and Sunyaev-Zel'dovich (SZ) observables for systems of merging relaxed X-ray clusters. In all of the scaling relations generated from the properties we have studied, we observe a generic evolution (in good qualitative agreement with previous authors): a rapid transient roughly along the mass scaling relations, a subsequent slow drift across the scatter until virialization, followed by a slow evolution along and up the mass scaling relations as cooling recovers in the cluster cores. However, this drift is not sufficient to account for the observed scatter in the scaling relations. We also study the effects of mergers on several theoretical temperature measures of the intracluster medium: emission weighted measures (T_ew), the spectroscopic-like measure proposed by Mazzotta et al (2004; T_sl) and plasma model fits to the integrated spectrum of the system (T_spec). We find that T_sl tracks T_spec for the entire duration of our mergers, illustrating that it remains a good tool for observational comparison even for highly disturbed systems. Furthermore, the transient temperature increases produced during first and second pericentric passage are 15-40% larger for T_ew than for T_sl or T_spec. This suggests that the effects of transient temperature increases on sigma_8 and Omega_M derived by Randall et al (2002) are over estimated. Lastly, we examine the X-ray SZ proxy proposed by Kravtsov et al (2006) and find that the tight mass scaling relation they predict remains secure through the entire duration of a merger event, independent of projection effects.

We discuss several results made possible by accurate SDSS astrometric measurements in a large sky area, with emphasis on asteroids and stellar proper motions obtained by comparing POSS and SDSS. SDSS has observed over 200,000 moving objects in five photometric bands, corresponding to about two orders of magnitude increase over previous multi--color surveys. These data were used to extend the measurement of asteroid size distribution to a smaller size limit, to demonstrate that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space, and to discover a correlation between asteroid age and colors. A preliminary analysis of SDSS-POSS proper motions for about 1 million M dwarf stars demonstrates that, in the 0.1-1 kpc distance range, the rotational velocity and its dispersion for disk stars increase with the distance from the Galactic plane.

The water content and habitability of terrestrial planets are determined during their final assembly, from perhaps a hundred 1000-km "planetary embryos" and a swarm of billions of 1-10 km "planetesimals." During this process, we assume that water-rich material is accreted by terrestrial planets via impacts of water-rich bodies that originate in the outer asteroid region. We present analysis of water delivery and planetary habitability in five high-resolution simulations containing about ten times more particles than in previous simulations (Raymond et al 2006a, Icarus, 183, 265-282). These simulations formed 15 terrestrial planets from 0.4 to 2.6 Earth masses, including five planets in the habitable zone. Every planet from each simulation accreted at least the Earth's current water budget; most accreted several times that amount (assuming no impact depletion). Each planet accreted at least five water-rich embryos and planetesimals from past 2.5 AU; most accreted 10-20 water-rich bodies. We present a new model for water delivery to terrestrial planets in dynamically calm systems, with low-eccentricity or low-mass giant planets -- such systems may be very common in the Galaxy. We suggest that water is accreted in comparable amounts from a few planetary embryos in a "hit or miss" way and from millions of planetesimals in a statistically robust process. Variations in water content are likely to be caused by fluctuations in the number of water-rich embryos accreted, as well as from systematic effects such as planetary mass and location, and giant planet properties.

We examine the origin and evolution of the mass-metallicity relationship (MZR, M-Z) for galaxies using high resolution cosmological SPH + N-Body simulations that include a physically motivated description of supernovae feedback and subsequent metal enrichment. We discriminate between two sources that may contribute to the origin of the MZR: 1) metal and baryon loss due to gas outflow, or 2) inefficient star formation at the lowest galaxy masses. Our simulated galaxies reproduce the observed MZR in shape and normalization both at z=0 and z=2. We find that baryon loss occurs due to UV heating before star formation turns on in galaxies with M_baryon < 10^8 M_sun, but that some gas loss due to supernovae induced winds is required to subsequently reproduce the low effective chemical yield observed in low mass galaxies. Despite this, we show that low star formation efficiencies, regulated by supernovae feedback, are primarily responsible for the lower metallicities of low mass galaxies and the overall M-Z trend. We find that the shape of the MZR is relatively constant with redshift, but that its normalization increases with time. Simulations with no energy feedback from supernovae overproduce metals at low galaxy masses by rapidly transforming a large fraction of their gas into stars. Despite the fact that our low mass galaxies have lost a majority of their baryons, they are still the most gas rich objects in our simulations due to their low star formation efficiencies.

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N--Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.

We present an analysis of star formation and feedback recipes appropriate for galactic smoothed particle hydrodynamics simulations. Using an isolated Milky Way-like galaxy, we constrain these recipes based on well-established observational results. Our star formation recipe is based on that of Katz (1992) with the additional inclusion of physically motivated supernova feedback recipes. We propose a new feedback recipe in which type II supernovae are modelled using an analytical treatment of blastwaves. With this feedback mechanism and a tuning of other star formation parameters, the star formation in our isolated Milky Way-like galaxy is constant and follows the slope and normalisation of the observed Schmidt law. In addition, we reproduce the low density cutoff and filamentary structure of star formation observed in disk galaxies. Our final recipe will enable better comparison of cosmological N-body simulations with observations.

We report on the analysis of a suite of SPH simulations (incorporating cooling and star formation) of mergers involving idealised X-ray clusters whose initial conditions resemble relaxed clusters with cool compact cores observed by Chandra and XMM. The simulations sample the most interesting, theoretically plausible, range of impact parameters and progenitor mass ratios. We find that all mergers evolve via a common progression. We illustrate this progression in the projected gas density, X-ray surface brightness, SZ, temperature, and gas entropy maps. Several different classes of transient ``cold front''-like features can arise over the course of a merger. We find that all of these classes are present in Chandra and XMM observations of merging systems and propose a naming scheme for these features: ``comet-like'' tails, bridges, plumes, streams and edges. In none of the cases considered do the initial cool compact cores of the primary and the secondary get destroyed during the course of the mergers. We quantify the evolving morphology of our mergers using centroid variance, power ratios and offset between the X-ray and the projected mass maps. We find that the centroid variance best captures the dynamical state of the cluster. Placing the system at z=0.1, we find that all easily identified observable traces of the secondary disappear from a simulated 50 ks Chandra image following the second pericentric passage. The system, however, takes approximately 2 additional Gyrs to relax and virialize. Temperature fluctuations at the level of 20% can persist in the final systems well past the point of virialization, suggesting that that the existence of temperature fluctuations, in and of themselves, do not necessarily indicate a disturbed or unrelaxed system.

Two well studied dwarf galaxies -- NGC 3109 and NGC 6822 -- present some of the strongest observational support for a flat core at the center of galactic dark matter (DM) halos. We use detailed cosmologically motivated numerical models to investigate the systematics and the accuracy of recovering parameters of the galaxies. Some of our models match the observed structure of the two galaxies remarkably well. Our analysis shows that the rotation curves of these two galaxies are instead quite compatible with their DM halos having steep cuspy density profiles. The rotation curves in our models are measured using standard observational techniques. The models reproduce the rotation curves of both galaxies, the disk surface brightness profiles as well as the profile of isophotal ellipticity and position angle. The models are centrally dominated by baryons; however, the dark matter component is globally dominant. The simulated disk mass is marginally consistent with a stellar mass-to-light ratio in agreement with the observed colors. We show that non-circular motions combined with gas pressure support and projection effects results in a large underestimation of the circular velocity in the central $\sim 1$ kpc region, creating the illusion of a constant density core. Although the systematic effects mentioned above are stronger in barred systems, they are also present in axisymetric disks. Our results strongly suggest that there is no contradiction between the observed rotation curves in dwarf galaxies and the cuspy central dark matter density profiles predicted by Cold Dark Matter models.

We present a comprehensive series of $N$-body as well as $N$-body + SPH simulations to study the secular evolution of the structure of disk galaxies. Our simulations are organized in a hierarchy of increasing complexity, ranging from rigid-halo collisionless simulations to fully live simulations with gas and star formation. Comparisons between the different types of simulations allow us to isolate the role of various physical mechanisms. We focus on the evolution of systems expected in a LCDM universe. Our goal is to examine which structural properties of disk galaxies may result from secular evolution rather than from direct hierarchical assembly. In the vertical direction, we find that various mechanisms can lead to heating. The strongest heating occurs during the vertical buckling instability of a bar. Among the consequences of this instability is the formation of peanut-shaped bulges which produce clear kinematic signatures when observed face-on. We find that bars are robust structures that are not destroyed by buckling. They can be destroyed instead by a central mass concentration but we find that this mass needs to be a large fraction of the total mass of the disk. We then study the evolution of stellar surface density profiles showing how angular momentum redistribution leads to increasing central densities and disk scale lengths and to profile breaks at large radii. The breaks in these simulations are in excellent agreement with observed breaks, even when the evolution is purely collisionless. Disk scale-lengths increase even when the total disk angular momentum is conserved; thus mapping halo angular momenta to scale-lengths is non-trivial. [Abridged]

We report on the results of the first 3D SPH simulation of massive, gravitationally unstable protoplanetary disks with radiative transfer. We adopt a flux-limited diffusion scheme justified by the high opacity of most of the disk. The optically thin surface of the disk cools as a blackbody. The disks grow slowly in mass starting from a Toomre-stable initial condition to the point at which they become marginally unstable. We find that gravitationally bound clumps with masses close to a Jupiter mass can arise. Fragmentation appears to be driven by vertical convective-like motions capable of transporting the heat from the disk midplane to its surface on a timescale of only about 40 years at 10 AU. A larger or smaller cooling efficiency of the disk at the optically thin surface can promote or stifle fragmentation by affecting the vertical temperature profile, which determines whether convection can happen or not, and by regulating the accretion flow from optically thin regions towards overdense regions. We also find that the chances of fragmentation increase for a higher mean molecular weight $\mu$ since compressional heating is reduced. Around a star with mass $1 M_{\odot}$ only disks with $\mu \ge 2.4$, as expected for gas with a metallicity comparable to solar or higher, fragment. This suggests that disk instability, like core-accretion, should be more effective in forming gas giants at higher gas metallicities, consistent with the observed correlation between metallicity of the planet-hosting stars and frequency of planets

Protoplanetary gas disks are likely to experience gravitational instabilites (GI's) during some phase of their evolution. Density perturbations in an unstable disk grow on a dynamic time scale into spiral arms that produce efficient outward transfer of angular momentum and inward transfer of mass through gravitational torques. In a cool disk with rapid enough cooling, the spiral arms in an unstable disk form self-gravitating clumps. Whether gas giant protoplanets can form by such a disk instability process is the primary question addressed by this review. We discuss the wide range of calculations undertaken by ourselves and others using various numerical techniques, and we report preliminary results from a large multi-code collaboration. Additional topics include -- triggering mechanisms for GI's, disk heating and cooling, orbital survival of dense clumps, interactions of solids with GI-driven waves and shocks, and hybrid scenarios where GI's facilitate core accretion. The review ends with a discussion of how well disk instability and core accretion fare in meeting observational constraints.

We present the results of the first multi-scale N-Body+SPH simulations of merging galaxies containing central supermassive black holes (SMBHs) and having a spatial resolution of only a few parsecs. Strong gas inflows associated with equal-mass mergers produce non-axisymmetric nuclear disks with masses of order $10^9 M_{\odot}$, resolved by about $10^6$ SPH particles. Such disks have sizes of several hundred parsecs but most of their mass is concentrated within less than $50$ pc. We find that a close SMBH pair forms after the merger. The separation of the two SMBHs then shrinks further owing to dynamical friction against the predominantly gaseous background. The orbits of the SMBHs decay down to the minimum resolvable scale in a few million years for an ambient gas temperature and density typical of a region undergoing a starburst. These results suggest the initial conditions necessary for the eventual coalescence of the two holes arise naturally from the merging of two equal-mass galaxies whose structure and orbits are consistent with the predictions of the $\Lambda$CDM model. Our findings have important implications for planned gravitational wave detection experiments such as {\it LISA}.

The final stage in the formation of terrestrial planets consists of the accumulation of ~1000-km ``planetary embryos'' and a swarm of billions of 1-10 km ``planetesimals.'' During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ~0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: 1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. 2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. 3) If planetary embryos form relatively slowly, following the models of Kokubo & Ida, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. 4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.

We used fully cosmological, high resolution N-body+SPH simulations to follow the formation of disk galaxies with a rotational velocity between 140 and 280 Km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation (SF), the effects of a uniform UV background and a physically motivated description of feedback from supernovae (SN). Feedback parameters have been chosen to match the star formation rate and interstellar medium (ISM) properties of local galaxies. In cosmological simulations galaxies formed rotationally supported disks with realistic exponential scale lengths and fall on the I-band and baryonic Tully Fisher relations. The combination of UV background and SN feedback drastically reduced the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Feedback delays SF in small galaxies and more massive ones contain older stellar populations. Here we focus on the SF and feedback implementations. We also briefly discuss how high mass and force resolution and a realistic description of SF and feedback are important ingredients to match the observed properties of galaxies.

We use cosmological LCDM numerical simulations to model the evolution of the substructure population in sixteen dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10^11 to 10^14 Msun at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation: dn/dv=(1/8)(v_cmax/v_cmax,host)^-4. Halo to halo variance in the VDF is a factor of roughly two to four. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for r >~ 0.3r_vir, where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease toward the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias toward more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power law cosmologies of P ~ k^n, and demonstrate that a steeper spectral index, n, results in significantly less substructure.

The ``minimum-mass solar nebula'' (MMSN) model estimates the surface density distribution of the protoplanetary disk by assuming the planets to have formed in situ. However, significant radial migration of the giant planets likely occurred in the Solar system, implying a distortion in the values derived by the MMSN method. The true density profiles of protoplanetary disks is therefore uncertain. Here we present results of simulations of late-stage terrestrial accretion, each starting from a disk of planetary embryos. We assume a power-law surface density profile that varies with heliocentric distance r as r^-alpha, and vary alpha between 1/2 and 5/2 (alpha = 3/2 for the MMSN model). We find that for steeper profiles (higher values of alpha), the terrestrial planets (i) are more numerous, (ii) form more quickly, (iii) form closer to the star, (iv) are more massive, (v) have higher iron contents, and (vi) have lower water contents. However, the possibility of forming potentially habitable planets does not appear to vary strongly with alpha.

We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html

We analyze the effect of dissipation on the orbital evolution of supermassive black holes (SMBHs) using high-resolution self-consistent gasdynamical simulations of binary equal- and unequal-mass mergers of disk galaxies. The galaxy models are consistent with the LCDM paradigm of structure formation and the simulations include the effects of radiative cooling and star formation. We find that equal-mass mergers always lead to the formation of a close SMBH pair at the center of the remnant with separations limited solely by the adopted force resolution of ~ 100 pc. Instead, the final SMBH separation in unequal-mass mergers depends sensitively on how the central structure of the merging galaxies is modified by dissipation. In the absence of dissipation, the satellite galaxy can be entirely disrupted before the merger is completed leaving its SMBH wandering at a distance too far from the center of the remnant for the formation of a close pair. In contrast, we show that gas cooling facilitates the pairing process by increasing the resilience of the companion galaxy to tidal disruption. Moreover, we demonstrate that merging disk galaxies constructed to obey the M(BH)-sigma relation, move relative to it depending on whether they undergo a dissipational or collisionless merger, regardless of the mass ratio of the merging systems. Collisionless simulations reveal that remnants tend to move away from the mean relation highlighting the role of gas-poor mergers as a possible source of scatter. In dissipational mergers, the interplay between strong gas inflows associated with the formation of massive nuclear disks and the consumption of gas by star formation provides the necessary fuel to the SMBHs and allows the merger remnants to satisfy the relation.

Neutralino dark matter leads to the formation of numerous earth mass dark matter haloes at redshifts z\approx 60 (Diemand et al. 2005). These abundant CDM micro-haloes have cuspy density profiles that can easily withstand the Galactic tidal field at the solar radius. Zhao, Taylor, Silk & Hooper (astro-ph/0502049) concluded that ``...the majority of dark matter substructures with masses \sim 10^{-6}M_o will be tidally disrupted due to interactions with stars in the Galactic halo''. However these authors assumed a halo density of stars that is at least an order of magnitude higher than observed. We show that the appropriate application of the impulse approximation is to the regime of multiple encounters, not single disruptive events as adopted by Zhao et al., which leads to a survival time of several Hubble times. Therefore we do not expect the tidal heating by Galactic stars to affect the abundance of micro-haloes. Disk crossing will cause some mass loss but the central cores are likely to survive and could be detected as gamma-ray sources with proper motions of several arc minutes per year.

We use a multi million particle N-body + SPH simulation to follow the formation of a rich galaxy cluster in a Lambda+CDM cosmology, with the goal of understanding the origin and properties of intracluster stars. The simulation includes gas cooling, star formation, the effects of a uniform UVB and feedback from supernovae. Halos that host galaxies as faint as M_R = -19.0 are resolved by this simulation, which includes 85% of the total galaxy luminosity in a rich cluster. We find that the accumulation of intracluster light (ICL) is an ongoing process, linked to infall and stripping events. The unbound star fraction increases with time and is 20% at z = 0, consistent with observations of galaxy clusters. The surface brightness profile of the cD shows an excess compared to a de Vaucouleur profile near 200 kpc, which is also consistent with observations. Both massive and small galaxies contribute substantially to the formation of the ICL, with stars stripped preferentially from the outer parts of their stellar distributions. Simulated observations of planetary nebulae (PNe) show significant substructure in velocity space. Despite this, individual intracluster PNe might be useful mass tracers if more than 5 fields at a range of radii have measured line-of-sight velocities, where an accurate mass calculation depends more on the number of fields than the number of PNe measured per field. However, the orbits of IC stars are more anisotropic than those of galaxies or dark matter, which leads to a systematic underestimate of cluster mass relative to that calculated with galaxies, if not accounted for in dynamical models. Overall, the properties of ICL formed in a hierarchical scenario are in good agreement with current observations. (Abridged)

We use N-body hydrodynamical simulations to study the structure of disks in triaxial potentials resembling CDM halos. Our analysis focuses on the accuracy of the dark mass distribution inferred from rotation curves derived from simulated long-slit spectra. We consider a massless disk embedded in a halo with axis ratios of 0.5:0.6:1.0 and with its rotation axis aligned with the minor axis of the halo. Closed orbits for the gaseous particles deviate from coplanar circular symmetry, resulting in a variety of long-slit rotation curve shapes, depending on the orientation of the disk relative to the line of sight. Rotation curves may thus differ significantly from the spherically-averaged circular velocity profile of the dark matter halo. "Solid-body" rotation curves--typically interpreted as a signature of a constant density core in the dark matter distribution--are obtained about 25% of the time for random orientations although the dark matter follows the cuspy density profile proposed by Navarro, Frenk & White (NFW). We conclude that the discrepancies reported between the shape of the rotation curve of low surface brightness galaxies and the structure of CDM halos may be resolved once the complex effects of halo triaxiality on the dynamics of the gas component is properly taken into account.

`Hot jupiters,' giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form (Armitage 2003). We present results of simulations of terrestrial planet formation in the presence of hot jupiters, broadly defined as having orbital radii <= 0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g. 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with hot jupiters.

We review evidence that the census of Milky Way satellites similar to those known may be incomplete at low latitude due to obscuration and in the outer halo due to a decreasing sensitivity to dwarf satellites with distance. We evaluate the possible impact that incompleteness has on comparisons with substructure models by estimating corrections to the known number of dwarfs using empirical and theoretical models. If we assume that the true distribution of Milky Way satellites is uniform with latitude, then we estimate a 33% incompleteness in the total number of dwarfs due to obscuration at low latitude. Similarly, if we suppose that the radial distribution of Milky Way satellites matches that of M31, or that of the oldest sub-halos or the most massive sub-halos in a simulation, we estimate a total number of Milky Way dwarfs ranging from 1 -- 3 times the known population. Although the true level of incompleteness is quite uncertain, the fact that our extrapolations yield average total numbers of MW dwarfs that are realistically 1.5 -- 2 times the known population, shows that incompleteness needs to be taken seriously when comparing to models of dwarf galaxy formation. Interestingly, the radial distribution of the oldest sub-halos in a Lambda+CDM simulation of a Milky Way-like galaxy possess a close match to the observed distribution of M31's satellites, which suggests that reionization may be an important factor controlling the observability of sub-halos. We also assess the prospects for a new SDSS search for Milky Way satellites to constrain the possible incompleteness in the outer halo.

We use high resolution 3D SPH simulations to study the evolution of self-gravitating binary protoplanetary disks. Heating by shocks and cooling are included. We consider different orbital separations and masses of the disks and central stars. Isolated massive disks ($M \sim 0.1 M_{{\o}dot}$) fragment into protoplanets as a result of gravitational instability for cooling times comparable to the orbital time. Fragmentation does not occur in binary systems with a separation of about 60 AU. This is because efficient heating owing to strong tidally induced spiral shocks damps any overdensity. The resulting temperatures, above 200 K, would vaporize water ice in the outer disk, posing a problem even for the other model of giant planet formation, core-accretion. Light disks ($M \sim 0.01 M_{\odot}$) do not fragment but remain cold because their low self-gravity inhibits strong shocks. Core accretion would not be hampered in the latter. At separations of about 120 AU the efficiency of fragmentation by disk instability rises and approaches that in isolated systems. If disk instability is the main formation mechanism for giant planets, on going surveys targeting binary systems should find considerably fewer planets in systems with separations below 100 AU.

We have undertaken a thorough dynamical investigation of five extrasolar planetary systems using extensive numerical experiments. The systems Gl 777 A, HD 72659, Gl 614, 47 Uma and HD 4208 were examined concerning the question of whether they could host terrestrial like planets in their habitable zones (=HZ). First we investigated the mean motion resonances between fictitious terrestrial planets and the existing gas giants in these five extrasolar systems. Then a fine grid of initial conditions for a potential terrestrial planet within the HZ was chosen for each system, from which the stability of orbits was then assessed by direct integrations over a time interval of 1 million years. The computations were carried out using a Lie-series integration method with an adaptive step size control. This integration method achieves machine precision accuracy in a highly efficient and robust way, requiring no special adjustments when the orbits have large eccentricities. The stability of orbits was examined with a determination of the Renyi entropy, estimated from recurrence plots, and with a more straight forward method based on the maximum eccentricity achieved by the planet over the 1 million year integration. Additionally, the eccentricity is an indication of the habitability of a terrestrial planet in the HZ; any value of e>0.2 produces a significant temperature difference on a planet's surface between apoapse and periapse. The results for possible stable orbits for terrestrial planets in habitable zones for the five systems are summarized as follows: for Gl 777 A nearly the entire HZ is stable, for 47 Uma, HD 72659 and HD 4208 terrestrial planets can survive for a sufficiently long time, while for Gl 614 our results exclude terrestrial planets moving in stable orbits within the HZ.

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We compare results from a LambdaCDM simulation with a LambdaWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force resolution in both the gas and dark matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. The stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of that of the dark halo in the LambdaCDM and LambdaWDM models respectively. The LambdaWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the LambdaCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. (abridged)

We present a detailed observational strategy for finding Near Earth Objects (NEOs) with the Sloan Digital Sky Survey (SDSS) telescope. We investigate strategies in normal, unbinned mode as well as binning the CCDs 2x2 or 3x3, which affects the sky coverage rate and the limiting apparent magnitude. We present results from 1 month, 3 year and 10 year simulations of such surveys. For each cadence and binning mode, we evaluate the possibility of achieving the Spaceguard goal of detecting 90% of 1 km NEOs (absolute magnitude H <= 18 for an albedo of 0.1). We find that an unbinned survey is most effective at detecting H <= 20 NEOs in our sample. However, a 3x3 binned survey reaches the Spaceguard Goal after only seven years of operation. As the proposed large survey telescopes (PanStarss; LSST) are at least 5-10 years from operation, an SDSS NEO survey could make a significant contribution to the detection and photometric characterization of the NEO population.

We present results of numerical simulations which examine the dynamical stability of known planetary systems, a star with two or more planets. First we vary the initial conditions of each system based on observational data. We then determine regions of phase space which produce stable planetary configurations. For each system we perform 1000 ~1 million year integrations. We examine upsilon And, HD83443, GJ876, HD82943, 47UMa, HD168443, and the solar system (SS). We find that the resonant systems, 2 planets in a first order mean motion resonance, (HD82943 and GJ876) have very narrow zones of stability. The interacting systems, not in first order resonance, but able to perturb each other (upsilon And, 47UMa, and SS) have broad regions of stability. The separated systems, 2 planets beyond 10:1 resonance, (we only examine HD83443 and HD168443) are fully stable. Furthermore we find that the best fits to the interacting and resonant systems place them very close to unstable regions. The boundary in phase space between stability and instability depends strongly on the eccentricities, and (if applicable) the proximity of the system to perfect resonance. In addition to million year integrations, we also examined stability on ~100 million year timescales. For each system we ran ~10 long term simulations, and find that the Keplerian fits to these systems all contain configurations which may be regular on this timescale.

We investigate the mass profile of LambdaCDM halos using a suite of numerical simulations spanning five decades in halo mass, from dwarf galaxies to rich galaxy clusters. Our analysis confirms the proposal of Navarro, Frenk & White (NFW) that the shape of LambdaCDM halo mass profiles differs strongly from a power law and depends little on mass. The logarithmic slope of the spherically-averaged density profile, as measured by beta=-dln(rho)/dln(r), decreases monotonically towards the center and becomes shallower than isothermal (beta<2) inside a characteristic radius, r_{-2}. Although the fitting formula proposed by NFW provides a reasonably good approximation to the density and circular velocity profiles of individual halos, systematic deviations from the best NFW fits are also noticeable. Inside r_{-2}, the profile of simulated halos gets shallower with radius more gradually than predicted and, as a result, NFW fits tend to underestimate the dark matter density in these regions. This discrepancy has been interpreted as indicating a steeply divergent cusp, but our results suggest a different interpretation. We use the density and enclosed mass at our innermost resolved radii to place strong constraints on beta_{0}: density cusps as steep as r^{-1.5} are inconsistent with most of our simulations, although beta_{0}=1 is still consistent with our data. Our density profiles show no sign of converging to a well-defined asymptotic inner power law. We propose a simple formula that reproduces the radial dependence of the slope better than the NFW profile, and so may minimize errors when extrapolating our results inward to radii not yet reliably probed by numerical simulations.

We present three-dimensional self-gravitating smoothed-particle hydrodynamics (SPH) simulations of an isothermal gaseous disc interacting with an embedded planet. Discs of varying stability are simulated with planets ranging from 10 Earth-masses to 2 Jupiter-masses. The SPH technique provides the large dynamic range needed to accurately capture the large scale behavior of the disc as well as the small scale interaction of the planet with surrounding material. Most runs used 10^5 gas particles, giving us the spatial resolution required to observe the formation of planets. We find four regions in parameter space: low-mass planets undergo Type I migration; higher-mass planets can form a gap; the gravitational instability mode of planet formation in marginally stable discs can be triggered by embedded planets; discs that are completely unstable can fragment to form many planets. The disc stability is the most important factor in determing which interaction a system will exhibit. For the stable disc cases, our migration and accretion time-scales are shorter and scale differently than previously suggested.

We present three-dimensional SPH calculations of giant planets embedded in gaseous disks. Our findings are collected into a map of parameter space, exhibiting four distinct regions: Type I migration, gap formation, triggered formation of more planets, and wholly unstable disks. For Type I migration of the planet due to secular interactions with the disk material, the migration rate depends linearly on the disk mass, and is independent of the initial planet mass. For more massive disks, the planet can disturb the disk strongly enough to trigger the collapse of gas into additional giant planets. When additional planets form, their interaction as point masses dominates the subsequent behavior of the system. This mechanism allows for the rapid formation of Jupiter-mass and higher planets. Migration due to interaction with the disk can significantly change the orbits of giant planets in gas disks.

We present galaxy circular velocity functions (GCVFs) for 34 low redshift (z < 0.15) clusters identified in the Sloan Digital Sky Survey (SDSS), for fifteen clusters drawn from dark matter simulations of hierarchical structure growth in a LambdaCDM cosmology, and for ~29,000 SDSS field galaxies. We find that the observed and simulated cluster GCVFs take the form of a power law. The cumulative GCVFs of the simulated clusters are very similar across a wide range of cluster masses, provided individual subhalo circular velocities are scaled by the circular velocities of the parent cluster. Once all sources of error are accounted for, the intrinsic scatter in the cumulative, scaled observed cluster GCVF is consistent with the simulations. The slope of the observed cluster GCVF is -2.5, independent of cluster velocity dispersion. The average slope of the simulated GCVFs is somewhat steeper, although formally consistent given the errors. Using our highest resolution cluster, we find that the effects of baryons on galaxy rotation curves is to flatten the simulated cluster GCVF into better agreement with observations. Finally, we find that the field GCVF deviates significantly from a power law, being flatter than the cluster GCVF at circular velocities less than 200 km/s, and steeper at circular velocities greater than 200 km/s.

We present the results of high resolution SPH simulations of the evolution of gravitationally unstable protoplanetary disks. We report on calculations in which the disk is evolved using a locally isothermal or adiabatic equation of state (with shock heating), and also on new simulations in which cooling and heating by radiation are explicitly modeled. We find that disks with a minimum Toomre parameter $< 1.4$ fragment into several gravitationally bound protoplanets with masses from below to a few Jupiter masses. This is confirmed also in runs where the disk is given a quiet start, growing gradually in mass over several orbital times. A cooling time comparable to the orbital time is needed to achieve fragmentation, for disk masses in the range $0.08-0.1 M_{\odot}$. After about 30 orbital times, merging between the bound condensations always leads to 2-3 protoplanets on quite eccentric orbits.

We carry out a large set of very high resolution, three dimensional smoothed particle hydrodynamics (SPH) simulations describing the evolution of gravitationally unstable gaseous protoplanetary disks. We consider a broad range of initial disk parameters. Disk masses out to 20 AU range from 0.075 to 0.125 $M_{\odot}$, roughly consistent with the high-end of the mass distribution inferred for disks around T Tauri stars. The initial disks span minimum $Q$ parameters between 0.8 and 2, with most models being around $\sim 1.4$. The disks are evolved assuming either a locally isothermal equation of state or an adiabatic equation of state with varying $\gamma$. When overdensities above a specific threshold appear as a result of gravitational instability in a locally isothermal calculation, the equation of state is switched to adiabatic to account for the increased optical depth. We show that when a disk has a minimum $Q$ parameter less than 1.4 strong trailing spiral instabilities, typically three or four armed modes, form and grow until fragmentation occurs along the arms after about 5 mean disk orbital times. The resulting clumps contract quickly to densities several orders of magnitude higher than the initial disk density, and the densest of them survive even under adiabatic conditions. These clumps are stable to tidal disruption and merge quickly, leaving 2-3 protoplanets on fairly eccentric orbits (the mean eccentricity being around 0.2) with masses between 0.7 and more than $7 M_{Jup}$, well in agreement with those of detected extrasolar planets.after $\sim 10^3$ years. Fragmentation is not strongly dependent on whether the disk starts from a marginally unstable state or gradually achieves it, as shown by a test where the disk gradually achieves the critical $Q$ by growing in mass.

We use a set of high-resolution cosmological N-body simulations to investigate the inner mass profile of galaxy-sized cold dark matter (CDM) halos. These simulations extend the thorough numerical convergence study presented in Paper I of this series (Power et al. 2003), and demonstrate that the mass profile of CDM halos can be robustly estimated beyond a minimum converged radius of order r_conv ~ 1 kpc/h in our highest resolution runs. The density profiles of simulated halos become progressively shallow from the virial radius inwards, and show no sign of approaching a well-defined power-law behaviour near the centre. At r_conv, the logarithmic slope of the density profile is steeper than the asymptotic \rho \propto r^-1 expected from the formula proposed by Navarro, Frenk, and White (1996), but significantly shallower than the steeply divergent \rho \propto r^-1.5 cusp proposed by Moore et al. (1999). We perform a direct comparison of the spherically-averaged dark matter circular velocity (V_c) profiles with rotation curves of low surface brightness (LSB) galaxies from the samples of de Blok et al. (2001), de Blok and Bosma (2002), and Swaters et al. (2003). Most (about two-thirds) LSB galaxies in this dataset are roughly consistent with CDM halo V_c profiles. However, about one third of LSBs in these samples feature a sharp transition between the rising and flat part of the rotation curve that is not seen in the V_c profiles of CDM halos. This discrepancy has been interpreted as excluding the presence of cusps, but we argue that it might simply reflect the difference between circular velocity and gas rotation speed likely to arise in gaseous disks embedded within realistic, triaxial CDM halos.

We use a high resolution $\Lambda$CDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen, in a 50 $h^{-1}$Mpc volume containing 80 million particles. Our low redshift results allow us to probe low $\sigma$ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth and Tormen mass function provides an excellent match to all of our data except for redshifts of ten and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 percent for the $10^{10}-10^{11} \msun$ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth and Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high redshift (z $\simeq$ 6) QSOs, which are thought to be hosted by haloes at 5$\sigma$ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50%. We test the sensitivity of our results to force accuracy, starting redshift, and halo finding algorithm.

We present results from 42 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary "embryos" (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 42 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 "habitable" planets between 0.9 and 1.1 AU. These planets range from dry worlds to "water worlds" with 100+ oceans of water (1 ocean = 1.5x10^24 g), and vary in mass between 0.23 and 3.85 Earth masses. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call "super embryos," 1-2 Earth mass protoplanets which can serve as the accretion seeds for 2+ Earth mass planets with large water contents.

We use a high resolution collisionless simulation of a Virgo--like cluster in a $\Lambda$CDM cosmology to determine the velocity and clustering properties of the diffuse stellar component in the intracluster region at the present epoch. The simulated cluster builds up hierarchically and tidal interactions between member galaxies and the cluster potential produce a diffuse stellar component free-flying in the intracluster medium. Here we adopt an empirical scheme to identify tracers of the stellar component in the simulation and hence study its properties. We find that at $z=0$ the intracluster stellar light is mostly unrelaxed in velocity space and clustered in structures whose typical clustering radii are about 50 kpc at R=400--500 kpc from the cluster center, and predict the radial velocity distribution expected in spectroscopic follow-up surveys. Finally, we compare the spatial clustering in the simulation with the properties of the Virgo intracluster stellar population, as traced by ongoing intracluster planetary nebulae surveys in Virgo. The preliminary results indicate a substantial agreement with the observed clustering properties of the diffuse stellar population in Virgo.

The key algorithms and features of the Gasoline code for parallel hydrodynamics with self-gravity are described. Gasoline is an extension of the efficient Pkdgrav parallel N-body code using smoothed particle hydrodynamics. Accuracy measurements, performance analysis and tests of the code are presented. Recent successful Gasoline applications are summarized. These cover a diverse set of areas in astrophysics including galaxy clusters, galaxy formation and gas-giant planets. Future directions for gasdynamical simulations in astrophysics and code development strategies for tackling cutting edge problems are discussed.

We study the evolution of cold, gravitationally unstable protoplanetary gaseous disks performing 3D SPH simulations with up to a million particles on large parallel machines. We show that self-gravitating protoplanets can form in disks with initial $Q_{min} \sim 1.4$ provided that disks are evolved with a locally isothermal equation of state until localized regions reach a sufficiently high density. Lower resolution simulations carried out for as long as a thousand years show that systems of giant protoplanets with masses and orbits comparable to those of the observed extrasolar planets arise naturally.

We present a comprehensive set of convergence tests which explore the role of various numerical parameters on the equilibrium structure of a simulated dark matter halo. We report results obtained with two independent, state-of-the-art, multi-stepping, parallel N-body codes: PKDGRAV and GADGET. We find that convergent mass profiles can be obtained for suitable choices of the gravitational softening, timestep, force accuracy, initial redshift, and particle number. For softenings chosen so that particle discreteness effects are negligible, convergence in the circular velocity is obtained at radii where the following conditions are satisfied: (i) the timestep is much shorter than the local orbital timescale; (ii) accelerations do not exceed a characteristic acceleration imprinted by the gravitational softening; and (iii) enough particles are enclosed so that the collisional relaxation timescale is longer than the age of the universe. The most stringent requirement for convergence is typically that imposed on the particle number by the collisional relaxation criterion, which implies that in order to estimate accurate circular velocities at radii where the density contrast may reach $\sim 10^6$, the region must enclose of order 3000 particles (or more than a few times $10^6$ within the virial radius). Applying these criteria to a galaxy-sized $\Lambda$CDM halo, we find that the spherically-averaged density profile becomes progressively shallower from the virial radius inwards, reaching a logarithmic slope shallower than -1.2 at the innermost resolved point, $r \sim 0.005 r_{200}$, with little evidence for convergence to a power-law behaviour in the inner regions.

The evolution of gravitationally unstable protoplanetary gaseous disks has been studied with the use of three-dimensional smoothed particle hydrodynamics simulations with unprecedented resolution. We have considered disks with initial masses and temperature profiles consistent with those inferred for the protosolar nebula and for other protoplanetary disks. We show that long-lasting, self-gravitating protoplanets arise after a few disk orbital periods if cooling is efficient enough to maintain the temperature close to 50 K. The resulting bodies have masses and orbital eccentricities similar to those of detected extrasolar planets.

We investigate the evolution of substructure in cold dark matter halos using N-body simulations of tidal stripping of substructure halos (subhalos) within a static host potential. We find that halos modeled following the Navarro, Frenk & White (NFW) mass profile lose mass continuously due to tides from the massive host, leading to the total disruption of satellite halos with small tidal radii. The structure of stripped NFW halos depends mainly on the fraction of mass lost, and can be expressed in terms of a simple correction to the original NFW profile. We apply these results to substructure in the Milky Way, and conclude that the dark matter halos surrounding its dwarf spheroidal (dSph) satellites have circular velocity curves that peak well beyond the luminous radius at velocities significantly higher than expected from the stellar velocity dispersion. Our modeling suggests that the true tidal radii of dSphs lie well beyond the putative tidal cutoff observed in the surface brightness profile, suggesting that the latter are not really tidal in origin but rather features in the light profile of limited dynamical relevance. For Draco, in particular, our modeling implies that its tidal radius is much larger than derived by Irwin & Hatzidimitriou (1995), lending support to the interpretation of recent Sloan survey data by Odenkirchen et al. (2001). Similarly, our model suggests that Carina's halo has a peak circular velocity of ~55 km/s, which may help explain how this small galaxy has managed to retain enough gas to undergo several bursts of star formation. Our results imply a close correspondence between the most massive subhalos expected in a CDM universe and the known satellites of the Milky Way, and suggest that only subhalos with peak circular velocities below 35 km/s lack readily detectable luminous counterparts.

Using a catalog of 147,986 galaxy redshifts and fluxes from the Sloan Digital Sky Survey (SDSS) we measure the galaxy luminosity density at z=0.1 in five optical bandpasses corresponding to the SDSS bandpasses shifted to match their restframe shape at z=0.1. We denote the bands {0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}, with \lambda_{eff} = [3216, 4240, 5595, 6792, 8111] Angstroms respectively. We use a maximum likelihood method which allows for a general form for the shape of the luminosity function, simple luminosity and number evolution, incorporates flux uncertainties, and accounts for the flux limits of the survey. We find luminosity densities at z=0.1 in absolute AB magnitudes in a Mpc^3 of [-14.10 \pm 0.15, -15.18 \pm 0.03, -15.90 \pm 0.03, -16.24 \pm 0.03, -16.56 \pm 0.02] in [{0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}], respectively, using \Omega_0 =0.3, \Omega_\Lambda=0.7, and h=1, and using Petrosian magnitudes. Similar results are obtained using Sersic model magnitudes, suggesting that flux from outside the Petrosian apertures is not a major correction. In the {0.1}{r} band, the best fit Schechter function to our results has \phi_\ast = (1.49 \pm 0.04) \times 10^{-2} h^3 Mpc^{-3}, M_\ast - 5\log_{10} h = -20.44 \pm 0.01, and \alpha = -1.05\pm 0.01. In solar luminosities, the luminosity density in {0.1}{r} is (1.84 \pm 0.04) h 10^8 L_{{0.1}{r},\odot} Mpc^{-3}. Our results are consistent with other estimates of the luminosity density, from the Two-degree Field Galaxy Redshift Survey and the Millenium Galaxy Catalog. They represent a substantial change (\sim 0.5 mag) from earlier SDSS luminosity density results based on commissioning data, almost entirely because of the inclusion of evolution in the luminosity function model.

We announce the first public release of the SDSS Moving Object Catalog, with SDSS observations for 58,117 asteroids. The catalog lists astrometric and photometric data for moving objects observed prior to Dec 15, 2001, and also includes orbital elements for 10,592 previously known objects. We analyze the correlation between the orbital parameters and optical colors for the known objects, and confirm that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space. Their distinctive optical colors indicate that the variations in chemical composition within a family are much smaller than the compositional differences between families, and strongly support earlier suggestions that asteroids belonging to a particular family have a common origin.

We discuss optical colors of 10,592 asteroids with known orbits selected from a sample of 58,000 moving objects observed by the Sloan Digital Sky Survey (SDSS). This is more than ten times larger sample that includes both orbital parameters and multi-band photometric measurements than previously available. We confirm that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space. In particular, we demonstrate that the three major asteroid families (Eos, Koronis, and Themis), together with the Vesta family, represent four main asteroid color types. Their distinctive optical colors indicate that the variations in chemical composition within a family are much smaller than the compositional differences between families, and strongly support earlier suggestions that asteroids belonging to a particular family have a common origin. We estimate that over 90% of asteroids belong to families.

We present a set of gas-dynamical simulations of galaxy groups and clusters aimed at exploring the effect of non-gravitational heating. We use GASOLINE, a parallel Tree+SPH code, to simulate the formation of four cosmic halos with temperature 0.5<T<8 keV. Non-gravitational heating is implemented in two different ways: (1) by imposing a minimum entropy floor at a given redshift, 1<z<5; (2) by gradually heating gas, proportionally to the SN rate expected from semi-analytical modeling of galaxy formation. Our main results are the following. (a) An extra heating energy of about 1 keV per gas particle is required to reproduce the observed Lx-T relation, independent of whether it is provided so as to create an entropy floor of 50-100 keV cm^2, or is modulated in redshift; our SN feedback recipe provides only 1/3 keV/part. (b) The M-T relation is almost unaffected by non-gravitational heating and follows the M T^{3/2} scaling, with a normalization ~40% higher than observed, independent of the heating scheme. The inclusion of cooling in a run of a small group has the effects of increasing T_ew by ~30%, possibly reconciling simulated and observed M-T relations, and of decreasing Lx by ~40%. In spite of the inclusion of SN feedback energy, almost 40% of the gas becomes cold, in excess of current observational estimates. (abridged)

We positionally correlate asteroids from existing catalogs with a sample of $\about$18,000 asteroids detected by the Sloan Digital Sky Survey (SDSS, Ivezi\'{c} {\em et al.} 2001). We find 2641 unique matches, which represent the largest sample of asteroids with both accurate multi-color photometry and known orbital parameters. The matched objects are predominantly bright, and demonstrate that the SDSS photometric pipeline recovers \about90% of the known asteroids in the observed region. For the recovered asteroids we find a large offset (\about 0.4 mag) between Johnson V magnitudes derived from SDSS photometry and the predicted catalog-based visual magnitudes. This offset varies with the asteroid color from 0.34 mag for blue asteroids to 0.44 mag for red asteroids, and is probably caused by the use of unfiltered CCD observations in the majority of recent asteroid surveys. This systematic photometric error leads to an overestimate of the number of asteroids brighter than a given absolute magnitude limit by a factor of \about 1.7. The distribution of the matched asteroids in orbital parameter space indicates strong color segregation. We confirm that some families are dominated by a single asteroid type (e.g. the Koronis family by red asteroids and the Themis family by blue asteroids), while others appear to be a mixture of blue and red objects (e.g. the Nysa/Polana family). Asteroids with the bluest \iz colors, which can be associated with the Vesta family, show particularly striking localization in orbital parameter space.

We examine recent suggestions that substructure in cold dark matter (CDM) halos may be in conflict with the presence of thin, dynamically fragile stellar disks. N-body simulations of an isolated disk/bulge/halo model of the Milky Way that includes several hundred dark matter satellites with masses, densities and orbits derived from high-resolution cosmological CDM simulations indicate that substructure at $z=0$ plays only a minor dynamical role in the heating of the disk over several Gyrs. This is because the orbits of satellites in present-day CDM halos seldom take them near the disk, where their tidal effects are greatest. Unless the effects of substructure are very different at earlier times, our models suggest that substructure might not preclude virialized CDM halos from being acceptable hosts of thin stellar disks like that of the Milky Way.

Recent observations suggest that Carina and other nearby dwarf spheroidal galaxies are surrounded by unbound stars tidally stripped by the Milky Way. We run high-resolution N-Body simulations of dwarf galaxies orbiting within the Milky Way halo to determine if such observations can be explained with dark matter potentials as those implied by current structure formation models. We show that tidal forces acting on dwarfs with constant density cores or with cuspy profiles having a low concentration parameter ($c < 5$) lead to flat outer stellar density profiles like that of Carina for a variety of orbital configurations. On the contrary, it is more difficult to remove stars from cuspy dark matter halos with concentrations as high as predicted by CDM models at the mass scale of dwarf galaxies ($c \simgt 10$) and the data can only be reproduced assuming nearly radial orbits. Our simulations show that Carina is losing mass at a fractional rate $< 0.1$ Gyr$^{-1}$ and its mass-to-light ratio could be inflated by at most a factor of 2 due to unbound stars projected along the line of sight. We follow the evolution of the tidal debris within a triaxial clumpy cold dark matter Milky Way halo which causes differential precession and small scale heating of the stellar streams. This renders their use as a dynamical tracer of the Galactic potential practically useless, but does provide a novel test of the nature of the dark matter. Models with warm dark matter (WDM) or fluid dark matter (FDM) produce dwarf halos with lower central densities than CDM and would be consistent with the observed tidal tails even for orbits with eccentricities as low as indicated by current data on nearby dwarf spheroidals. Galactic halos in FDM are smooth and spherical and would be favored by the detection of coherent streams.

We present a new algorithm (PINOCCHIO, PINpointing Orbit-Crossing Collapsed HIerarchical objects) to predict accurately the formation and evolution of individual dark matter haloes in a given realization of an initial linear density field. Compared with the halo population formed in a large (360^3 particles) collisionless simulation of a CDM universe, our method is able to predict to better than 10 per cent statistical quantities such as the mass function, two-point correlation function and progenitor mass function of the haloes. Masses of individual haloes are estimated accurately as well, with errors typically of order 30 per cent in the mass range well resolved by the numerical simulation. These results show that the hierarchical formation of dark matter haloes can be accurately predicted using local approximations to the dynamics when the correlations in the initial density field are properly taken into account. The approach allows one to automatically generate a large ensemble of accurate merging histories of haloes with complete knowledge of their spatial distribution. The construction of the full merger tree for a 256^3 realisation requires a few hours of CPU-time on a personal computer, orders of magnitude faster than the corresponding N-body simulation would take, and without needing any extensive post-processing. The technique can be efficiently used, for instance, for generating the input for galaxy formation modeling.

We present results from high-resolution Tree+SPH simulations of galaxy clusters and groups, aimed at studying the effect of non-gravitational heating on the entropy of the ICM. We simulate three systems, having emission-weighted temperature T=0.6,1 and 3 keV, with spatial resolution better than 1% of the virial radius. We consider the effect of different prescriptions for non-gravitational ICM heating, such as SN energy feedback, as predicted by semi-analytical models of galaxy formation, and two different minimum entropy floors, S_fl=50 and 100 keV cm^2, imposed at z=3. Simulations with only gravitational heating nicely reproduce predictions from self-similar ICM models, while extra heating is shown to break the self-similarity, by a degree which depends on total injected energy and on cluster mass. We use observational results on the excess entropy in central regions of galaxy systems, to constrain the amount of extra-heating required. We find that setting the entropy floor S_fl=50 keV cm^2, which corresponds to an extra heating energy of about 1 keV per particle, is able to reproduce the observed excess of ICM entropy.

We use a cosmological simulation of the Local Group to make quantitative and speculative predictions for direct detection experiments. Cold dark matter (CDM) halos form via a complex series of mergers, accretion events and violent relaxation which precludes the formation of significant caustic features predicted by axially symmetric collapse. The halo density profiles are combined with observational constraints on the galactic mass distribution to constrain the local density of cold dark matter to lie in the range 0.18 <~ rho_CDM(R_solar)/GeV cm^-3 <~ 0.30. In velocity space, coherent streams of dark matter from tidally disrupted halos fill the halo and provide a tracer of the merging hierarchy. The particle velocities within triaxial CDM halos cannot be approximated by a simple Maxwellian distribution and is radially biased at the solar position. The detailed phase space structure within the solar system will depend on the early merger history of the progenitor halos and the importance of major mergers over accretion dominated growth. We follow the formation of a ``Draco'' sized dSph halo of mass 10^8M_solar with several million particles and high force accuracy. Its internal structure and substructure resembles that of galactic or cluster mass halos: the density profile has a singular central cusp and it contains thousands of sub-halos orbiting within its virial radius demonstrating a self-similar nature to collisionless dark matter sub-clustering. The singular cores of substructure halos always survive complete tidal disruption although mass loss is continuous and rapid. Extrapolating wildly to earth mass halos with velocity dispersion of 1 m s^-1 (roughly equal to the free streaming scale for neutralinos) we find that most of the dark matter may remain attached to bound subhalos. (Abridged)

We discuss measurements of the properties of about 10,000 asteroids detected in 500 deg2 of sky in the Sloan Digital Sky Survey (SDSS) commissioning data. The moving objects are detected in the magnitude range 14 < r < 21.5, with a baseline of 5 minutes. Extensive tests show that the sample is at least 98% complete, with the contamination rate of less than 3%. We find that the size distribution of asteroids resembles a broken power-law, independent of the heliocentric distance: D^{-2.3} for 0.4 km < D < 5 km, and D^{-4} for 5 km < D < 40 km. As a consequence of this break, the number of asteroids with r < 21.5 is ten times smaller than predicted by extrapolating the power-law relation observed for brighter asteroids (r < 18). The observed counts imply that there are about 530,000 objects with D>1 km in the asteroid belt, or about four times less than previous estimates. The distribution of main belt asteroids in the 4-dimensional SDSS color space is bimodal, and the two groups can be associated with S (rocky) and C (carbonaceous) asteroids. A strong bimodality is also seen in the heliocentric distribution of asteroids and suggests the existence of two distinct belts: the inner rocky belt, about 1 AU wide (FWHM) and centered at R~2.8 AU, and the outer carbonaceous belt, about 0.5 AU wide and centered at R~3.2 AU. The colors of Hungarias, Mars crossers, and near-Earth objects are more similar to the C-type than to S-type asteroids, suggesting that they originate in the outer belt. (abridged).

We present results from high-resolution N-Body/SPH simulations of rotationally supported dwarf irregular galaxies moving on bound orbits in the massive dark matter halo of the Milky Way.The dwarf models span a range in disk surface density and the masses and sizes of their dark halos are consistent with the predictions of cold dark matter cosmogonies. We show that the strong tidal field of the Milky Way determines severe mass loss in their halos and disks and induces bar and bending instabilities that transform low surface brightness dwarfs (LSBs) into dwarf spheroidals (dSphs) and high surface brightness dwarfs (HSBs) into dwarf ellipticals (dEs) in less than 10 Gyr. The final central velocity dispersions of the remnants are in the range 8-30 km/s and their final $v/\sigma$ falls to values $< 0.5$, matching well the kinematics of early-type dwarfs. The transformation requires the orbital time of the dwarf to be $\simlt 3-4$ Gyr, which implies a halo as massive and extended as predicted by hierarchical models of galaxy formation to explain the origin of even the farthest dSph satellites of the Milky Way, Leo I and Leo II. Only dwarfs with central dark matter densities as high as those of Draco and Ursa Minor can survive for 10 Gyr in the proximity of the Milky Way: this is naturally achieved within hierarchical models, where the densest objects should have small orbital times due to their early formation epochs. Part of the gas is stripped and part is funneled to the center due to the bar, generating one strong burst of star formation in HSBs and smaller, multiple bursts in LSBs. Extended low-surface brightness stellar and gaseous streams originate from LSBs and, when projected along the line of sight, can lead to overestimate the mass-to-light ratio of the bound remnant by a factor $\simlt 2$,

An all high-latitude sky survey for cool carbon giant stars in the Galactic halo has revealed 75 such stars, of which the majority are new detections. Of these, more than half are clustered on a Great Circle on the sky which intersects the center of Sagittarius dwarf galaxy (Sgr) and is parallel to its proper motion vector, while many of the remainder are outlying Magellanic Cloud C-stars. A pole-count analysis of the carbon star distribution clearly indicates that the Great Circle stream we have isolated is statistically significant, being a 5-6 sigma over-density. These two arguments strongly support our conclusion that a large fraction of the Halo carbon stars originated in Sgr. The stream orbits the Galaxy between the present location of Sgr, 16 kpc from the Galactic center, and the most distant stream carbon star, at ~60 kpc. It follows neither a polar nor a Galactic plane orbit, so that a large range in both Galactic R and z distances are probed. That the stream is observed as a Great Circle indicates that the Galaxy does not exert a significant torque upon the stream, so the Galactic potential must be nearly spherical in the regions probed by the stream. We present N-body experiments simulating this disruption process as a function of the distribution of mass in the Galactic halo. A likelihood analysis shows that, in the Galactocentric distance range 16 kpc < R < 60 kpc, the dark halo is most likely almost spherical. We rule out, at high confidence levels, the possibility that the Halo is significantly oblate, with isodensity contours of aspect q_m < 0.7. This result is quite unexpected and contests currently popular galaxy formation models. (Abridged)

N-Body/SPH simulations are used to study the evolution of dwarf irregular galaxies (dIrrs) entering the dark matter halo of the Milky Way or M31 on plunging orbits. We propose a new dynamical mechanism driving the evolution of gas rich, rotationally supported dIrrs, mostly found at the outskirts of the Local Group (LG), into gas free, pressure supported dwarf spheroidals (dSphs) or dwarf ellipticals (dEs), observed to cluster around the two giant spirals. The initial model galaxies are exponential disks embedded in massive dark matter halos and reproduce nearby dIrrs. Repeated tidal shocks at the pericenter of their orbit partially strip their halo and disk and trigger dynamical instabilities that dramatically reshape their stellar component. After only 2-3 orbits low surface brightness (LSB) dIrrs are transformed into dSphs, while high surface brightness (HSB) dIrrs evolve into dEs. This evolutionary mechanism naturally leads to the morphology-density relation observed for LG dwarfs. Dwarfs surrounded by very dense dark matter halos, like the archetypical dIrr GR8, are turned into Draco or Ursa Minor, the faintest and most dark matter dominated among LG dSphs. If disks include a gaseous component, this is both tidally stripped and consumed in periodic bursts of star formation. The resulting star formation histories are in good qualitative agreement with those derived using HST color-magnitude diagrams for local dSphs.

Can N-body simulations reliably determine the structural properties of dark matter halos? Focussing on a Virgo-sized galaxy cluster, we increase the resolution of current ``high resolution simulations'' by almost an order of magnitude to examine the convergence of the important physical quantities. We have 4 million particles within the cluster and force resolution 0.5 kpc/h (0.05% of the virial radius). The central density profile has a logarithmic slope of -1.5, as found in lower resolution studies of the same halo, indicating that the profile has converged to the ``physical'' limit down to scales of a few kpc. Also the abundance of substructure is consistent with that derived from lower resolution runs; on the scales explored, the mass and circular velocity functions are close to power laws of exponents ~ -1.9 and -4. Overmerging appears to be globally unimportant for suhalos with circular velocities > 100 km/s. We can trace most of the cluster progenitors from z=3 to the present; the central object (the dark matter analog of a cD galaxy)is assembled between z=3 and 1 from the merging of a dozen halos with v_circ \sim 300 km/s. The mean circular velocity of the subhalos decreases by ~ 20% over 5 billion years, due to tidal mass loss. The velocity dispersions of halos and dark matter globally agree within 10%, but the halos are spatially anti-biased, and, in the very central region of the cluster, they show positive velocity bias; however, this effect appears to depend on numerical resolution.

Because of the high eccentricities (~0.3) of two of the possible planets about the star Upsilon Andromeda, the stability of the system requires careful study. We present results of 1000 numerical simulations which explore the orbital parameter space as constrained by the observations. The orbital parameters of each planet are chosen from a Gaussian error distribution, and the resulting configuration is integrated for 1,000,000 years. We find that 84% of these integrations are stable. Configurations in which the eccentricity of the third planet is <0.3 are always stable, but when the eccentricity is >0.45 the system is always unstable, typically producing a close encounter between the second and third planets. A similar exercise with the gas giants in our own Solar System sampled with the same error distribution was performed. Approximately 81% of these simulations were stable for 10^6 years.

We combine semi-analytical methods with a ultra-high resolution simulation of a galaxy cluster (of mass 2.3 10^14h-1Msolar, and 4 10^6 particles within its virial radius) formed in a standard CDM universe to study the spatial distribution and orbital properties of the present-day descendents of Lyman Break Galaxies (LBGs). At the present time only five (out of 12) of halos containing LBGs survive as separate entities inside the cluster virial radius. Their circular velocities are in the range 200 - 550 km/sec. Seven halos merged together to form the central object at the very center of the cluster. Using semi-analytical modeling of galaxy evolution we show that descendents of halos containing LBGs now host giant elliptical galaxies. Galaxy orbits are radial, with a pericenter to apocenter ratio of about 1:5. The orbital eccentricities of LBGs descendents are statistically indistinguishable from those of the average galaxy population inside the cluster, suggesting that the orbits of these galaxies are not significantly affected by dynamical friction decay after the formation of the cluster's main body. In this cluster, possibly due to its early formation time, the descendents of LBGs are contained within the central 60% of the cluster virial radius and have an orbital velocity dispersion lower than the global galaxy population, originating a mild luminosity segregation for the brightest cluster members. Mass estimates based only on LBGs descendents (especially including the central cD) reflect this bias in space and velocity and underestimate the total mass of this well virialized cluster by up to a factor of two compared to estimates using at least 20 cluster members.

We present the analysis of a suite of simulations of a Virgo mass galaxy cluster. Undertaken within the framework of standard cold dark matter cosmology, these simulations were performed at differing resolutions and with increasingly complex physical processes, with the goal of identifying the effects of each on the evolution of the cluster. We focus on the cluster at the present epoch and examine properties including the radial distributions of density, temperature, entropy and velocity. We also map `observable' projected properties such as the surface mass density, X-ray surface brightness and SZ signature. We identify significant differences between the simulations, which highlights the need for caution when comparing numerical simulations to observations of galaxy clusters. While resolution affects the inner density profile in dark matter simulations, the addition of a gaseous component, especially one that cools and forms stars, affects the entire cluster. We conclude that both resolution and included physical processes play an important role in simulating the formation and evolution of galaxy clusters. Therefore, physical inferences drawn from simulations that do not include a gaseous component that can cool and form stars present a poor representation of reality. (Abridged)

There is increasing evidence that many km-sized bodies in the Solar System are piles of rubble bound together by gravity. We present results from a project to map the parameter space of collisions between km-sized spherical rubble piles. The results will assist in parameterization of collision outcomes for Solar System formation models and give insight into fragmentation scaling laws. We use a direct numerical method to evolve the positions and velocities of the rubble pile particles under the constraints of gravity and physical collisions. We test the dependence of the collision outcomes on impact parameter and speed, impactor spin, mass ratio, and coefficient of restitution. Speeds are kept low (< 10 m/s, appropriate for dynamically cool systems such as the primordial disk during early planet formation) so that the maximum strain on the component material does not exceed the crushing strength. We compare our results with analytic estimates and hydrocode simulations. Off-axis collisions can result in fast-spinning elongated remnants or contact binaries while fast collisions result in smaller fragments overall. Clumping of debris escaping from the remnant can occur, leading to the formation of smaller rubble piles. In the cases we tested, less than 2% of the system mass ends up orbiting the remnant. Initial spin can reduce or enhance collision outcomes, depending on the relative orientation of the spin and orbital angular momenta. We derive a relationship between impact speed and angle for critical dispersal of mass in the system. We find that our rubble piles are relatively easy to disperse, even at low impact speed, suggesting that greater dissipation is required if rubble piles are the true progenitors of protoplanets.

We use numerical simulations to examine the substructure within galactic and cluster mass halos that form within a hierarchical universe. Clusters are easily reproduced with a steep mass spectrum of thousands of substructure clumps that closely matches observations. However, the survival of dark matter substructure also occurs on galactic scales, leading to the remarkable result that galaxy halos appear as scaled versions of galaxy clusters. The model predicts that the virialised extent of the Milky Way's halo should contain about 500 satellites with circular velocities larger than Draco and Ursa-Minor i.e. bound masses > 10^8Mo and tidally limited sizes > kpc. The substructure clumps are on orbits that take a large fraction of them through the stellar disk leading to significant resonant and impulsive heating. Their abundance and singular density profiles has important implications for the existence of old thin disks, cold stellar streams, gravitational lensing and indirect/direct detection experiments.

We study the evolution of galaxy satellites with high resolution N-body simulations. Satellites are modeled as replicas of typical low and high surface brightness galaxies (LSBs and HSBs). Encounters on high eccentricity orbits (as typical in hierarchical models of galaxy formation) strip LSBs of most of their stars and tend to decrease their surface brightness. On the contrary, bar instability in HSBs leads to substantial loss of angular momentum of the stellar component and to an increase of central surface brightness. In both cases the remnant resembles a spheroidal galaxy with an exponential surface brightness profile. A simple modeling of color evolution and interactions driven star formation gives M/L ratios for the remnants that are roughly consistent with observations. These results suggest an evolutionary scenario for the dwarf galaxies in our Local Group, faint dSphs being the descendents of LSBs and brighter dSphs/dEs being the final state of HSB satellites.

We show that a universe dominated by cold dark matter fails to reproduce the rotation curves of dark matter dominated galaxies, one of the key problems that it was designed to resolve. We perform numerical simulations of the formation of dark matter halos, each containing \gsim 10^6 particles and resolved to 0.003 times the virial radius, allowing an accurate comparison with rotation curve data. A good fit to both galactic and cluster sized halos can be achieved using the density profile rho(r) \propto [(r/r_s)^1.5(1+(r/r_s)^1.5)]^-1, where r_s is a scale radius. This profile has a steeper asymptotic slope, rho(r) \propto r^-1.5, and a sharper turnover than found by lower resolution studies. The central structure of relaxed halos that form within a hierarchical universe has a remarkably small scatter (unrelaxed halos would not host disks). We compare the results with a sample of dark matter dominated, low surface brightness (LSB) galaxies with circular velocities in the range 100-300 km/s. The rotation curves of disks within cold dark matter halos rise too steeply to match these data which require a constant mass density in the central regions. The same conclusion is reached if we compare the scale free shape of observed rotation curves with the simulation data. It is important to confirm these results using stellar rather than HI rotation curves for LSB galaxies. We test the effects of introducing a cut-off in the power spectrum that may occur in a universe dominated by warm dark matter. In this case halos form by a monolithic collapse but the final density profile hardly changes, demonstrating that the merger history does not play a role in determining the halo structure.

We follow the evolution of disk galaxies within a cluster that forms hierarchically in a standard cold dark matter N-body simulation. At a redshift z=0.5 we select several dark matter halos that have quiet merger histories and are about to enter the newly forming cluster environment. The halos are replaced with equilibrium high resolution model spirals that are constructed to represent luminous examples of low surface brightness (LSB) and high surface brightness (HSB) galaxies. Whilst the models have the same total luminosity, \sim L_*, they have very different internal mass profiles, core radii and disk scale lengths, however they all lie at the same place on the Tully-Fisher relation. Due to their ``soft'' central potentials, LSB galaxies evolve dramatically under the influence of rapid encounters with dark matter halos, substructure and strong tidal shocks from the global cluster potential - galaxy harassment. As much as 90% of the LSB disk stars are tidally stripped and congregate in large diffuse tails that trace the orbital path of the galaxy and form the diffuse intra-cluster light. The bound stellar remnants closely resemble the dwarf spheroidals (dE's) that populate nearby clusters, with large scale lengths and low central surface brightness.

We follow the evolution of disk galaxies within a cluster that forms hierarchically in a cold dark matter N-body simulation. At a redshift z=0.5 we select several dark matter halos that have quiet merger histories and are about to enter the newly forming cluster environment. The halos are replaced with equilibrium high resolution model spirals that are constructed to represent examples of low surface brightness (LSB) and high surface brightness (HSB) galaxies. Varying the disk and halo structural parameters reveals that the response of a spiral galaxy to tidal encounters depends primarily on the potential depth of the mass distribution and the disk scale length. LSB galaxies, characterised by slowly rising rotation curves and large scale lengths, evolve dramatically under the influence of rapid encounters with substructure and strong tidal shocks from the global cluster potential --- galaxy harassment. We find that up to 90% of their stars are tidally stripped and congregate in large diffuse tails that trace the orbital path of the galaxy and form the diffuse intra-cluster light. The bound stellar remnants closely resemble the dwarf spheroidals (dE's) that populate nearby clusters. HSB galaxies are stable to the chaos of cluster formation and tidal encounters. These disks lie well within the tidally limited dark matter halos and their potentials are more concentrated. Although very few stars are stripped, the scale height of the disks increases substantially and no spiral features remain, therefore we speculate that these galaxies would be identified as S0 galaxies in present day clusters.

We analyse parallel N-body simulations of three Cold Dark Matter (CDM) universes to study the abundance and clustering of galaxy clusters. The simulations cover a volume comparable to the forthcoming SDSS. We are able to make robust measurements of cluster properties to a redshift larger than unity. We extract halos using two independent, public domain group finders (FOF & HOP) and find consistent results. The correlation function of clusters is in very good agreement with a simple analytic prescription based upon a Lagrangian biasing scheme developed by Mo & White (1996) and the Press-Schechter (PS) formalism for the mass function. The R_0--D_c relation for the open CDM model is in good agreement with the results from the APM Cluster Survey. The SCDM universe shows a robust deviation in the shape and evolution of the mass function when compared with that predicted by the PS formalism. Critical models with a low sigma_8 normalization or small shape parameter Gamma show an excess of massive clusters compared with the PS prediction. When cluster normalized, the SCDM universe at z =1 contains 10 times more clusters with temperatures greater than 7keV, compared with the PS prediction. The agreement between the analytic and N-body mass functions of SCDM can be improved if the value of the delta_c (the extrapolated linear theory threshold for collapse) is revised to be $ delta_c(z) = 1.685[(0.7/sigma_8)(1+z)]^{-0.125}. Our best estimate for the amplitude of fluctuations inferred from the local cluster abundance for SCDM is sigma_{8} = 0.5 \pm 0.04. However, the discrepancy between the temperature function predicted in a critical density universe and that observed at z=0.33 (Henry et al. 1998) remains. (abridged)

Higher order statistics are investigated in ($\Omega$)CDM universes by analyzing $500\mpc$ high resolution tree N-body simulations with both $\Omega = 1$, and $\Omega < 1$. The amplitudes of the N-point correlation functions are calculated from moments of counts-in-cells determined by a pair of new algorithms especially developed for large simulations. This approach enables massive oversampling with $\simeq 10^{9-14}$ cells for accurate determination of factorial moments from up to 47 million particles in the scale range of $8 \kpc - 125\mpc$. Thorough investigation shows that there are three scale ranges in the simulations: $\ge 8\mpc$, weakly non-linear regime, where perturbation theory applies with utmost precision, $1\mpc - 8\mpc$, the non-linear plateau, and finally $\le 1\mpc$, a regime where dynamical discreteness effects dominate the higher order statistics. In the physically relevant range of $1\mpc-125\mpc$ the results i) confirm the validity of perturbation theory in the weakly non-linear regime, ii) establish the existence of a plateau in the highly non-linear regime similar to the one observed in scale free simulations iii) show extended perturbation theory to be an excellent approximation for the non-linear regime iv) find the time dependence of the $S_N$'s to be negligible in both regimes v) in comparison with similar measurements in the EDSGC survey, strongly support $\Omega < 1$ with no biasing vi) show that the formulae of Szapudi & Colombi (1996) provide a good approximation for errors on higher order statistics measured in N-body simulations.

The discovery of a population of young galaxies at an epoch when the universe was about one tenth of its current age has shed new light on the question of when and how galaxies formed. Within the context of popular models this is the population of primeval galaxies that built themselves up to the size of present--day galaxies through the process of repeated mergers. But the recent detection of a large concentration of these primeval galaxies (Steidel et al. 1997) appears to be incompatible with hierarchical clustering models, which generally predict that clusters of this size are fully formed later in time. Here we use a combination of two powerful theoretical techniques --semi-analytic modelling and N-body simulations-- to show that such large concentrations should be quite common in a universe dominated by cold dark matter, and that they are the progenitors of the rich galaxy clusters seen today. We predict the clustering properties of primeval galaxies which should, when compared with data that will be collected in the near future, test our understanding of galaxy formation within the framework of a universe dominated by cold dark matter.

We examine the properties of dark matter halos within a rich galaxy cluster using a high resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 halos with circular velocities larger than 80 kms within the cluster's virial radius of 2 Mpc. This enables an unprecedented study of the statistical properties of a large sample of dark matter halos evolving in a dense environment. The cumulative fraction of mass attached to these halos varies from 0% at 200 kpc, to 13% at the virial radius. Even at this resolution the overmerging problem persists; halos that pass within 200 kpc of the cluster center are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor halos. The median ratio of apocentric to pericentric radii is 6:1; the orbital distribution is close to isotropic, circular orbits are rare, radial orbits are common. The orbits of halos are unbiased with respect to both position within the cluster and with the orbits of the smooth dark matter background and no velocity bias is detected. The tidal radii of surviving halos are generally well-fit using the simple analytic prediction applied to their orbital pericenters. Halos within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other halos that cause significant mass loss - ``galaxy harassment''. Mergers between halos do not occur inside the clusters virial radius.

We examine the effects of mass resolution and force softening on the density profiles of cold dark matter halos that form within cosmological N-body simulations. As we increase the mass and force resolution, we resolve progenitor halos that collapse at higher redshifts and have very high densities. At our highest resolution we have nearly 3 million particles within the virial radius, several orders of magnitude more than previously used and we can resolve more than one thousand surviving dark matter halos within this single virialised system. The halo profiles become steeper in the central regions and we may not have achieved convergence to a unique slope within the inner 10% of the virialised region. Results from two very high resolution halo simulations yield steep inner density profiles, $\rho(r)\sim r^{-1.4}$. The abundance and properties of arcs formed within this potential will be different from calculations based on lower resolution simulations. The kinematics of disks within such a steep potential may prove problematic for the CDM model when compared with the observed properties of halos on galactic scales.

We investigate the internal structure of cold dark matter halos using high resolution N-body simulations. As the mass and force resolution are increased, halo density profiles become steeper, asymptoting to a slope of $\sim -4/3$ in the central regions and may not have converged to a unique result. At our highest resolution we have nearly 3 million particles within the virial radius, $R_{200}$, and force softening that is $\sim 0.2% R_{200}$. This resolution has also allowed us to resolve a large part of the overmerging problem - we find over 1000 surviving dark halos orbiting within a single cluster potential. These data have given us unprecedented insights into the dynamics and structure of ``halos within halos'', allowing us for the first time, to compare the distribution of dark matter with observations of galaxies in clusters.

Using cosmological N-body simulations of critical (SCDM) and open (Omega_0=0.3, OCDM) cold dark matter models we have identified dark matter halos which are associated with the progenitors of present day bright early type galaxies. By following their merging history, we show how early type galaxies that formed within massive halos at redshift around 3 are now preferentially residing in clusters and groups. On the other hand, those that formed through major merging events between redshift 1 and the present have not yet been accreted into larger, virialized structures. This result is in agreement with analytical predictions in hierarchical clustering models. CDM models are able to explain both the ancient and uniform population of ellipticals that dominates in clusters together with the more recent and heterogeneous population of field ellipticals. Predictions for the comoving number density of bright early-type galaxies are given, and are shown to be consistent with the observed luminosity function. We predict that the number density of interacting bright binary galaxies, from which the field population of ellipticals may have originated, is proportional to (1+z)**4.2+-0.28 and (1+z)**2.5+-0.42 in SCDM and OCDM respectively. This result is consistent with previous analytical estimates and is discussed together with recent observational constraints.

Leapfrog integration has been the method of choice in N-body simulations owing to its low computational cost for a symplectic integrator with second order accuracy. We introduce a new leapfrog integrator that allows for variable timesteps for each particle in large N-body simulations. Tests with single particles in fixed potentials show that it behaves as a symplectic integrator. We then examine the results of both standard leapfrog and our temporally adaptive leapfrog on full N-body integrations of clusters and large scale structure establishing accuracy criteria for both methods. The adaptive method shows significant speed-ups over single step integrations---but the integrator no longer appears to be symplectic or, in the case of large scale structure simulations, accurate. This loss of accuracy appears to be caused by the way that the timestep is chosen, not by the integrator itself. We present a related integration technique that does retain sufficient accuracy. Although it is not symplectic, it is apparently better than previous implementations and is our current integrator of choice for large astrophysical simulations. We also note that the standard leapfrog difference equations used in cosmological N-body integrations in comoving coordinates are not symplectic. We derive an implementation of leapfrog that is in comoving canonical coordinates to correct for this deficiency.

The dynamics of the Local Group and its environment provide a unique challenge to cosmological models. The velocity field within 5h-1 Mpc of the Local Group (LG) is extremely ``cold''. The deviation from a pure Hubble flow, characterized by the observed radial peculiar velocity dispersion, is measured to be about 60km/s. We compare the local velocity field with similarly defined regions extracted from N-body simulations of Universes dominated by cold dark matter (CDM). This test is able to strongly discriminate between models that have different mean mass densities. We find that neither the Omega=1 (SCDM) nor Omega=0.3 (OCDM) cold dark matter models can produce a single candidate Local Group that is embedded in a region with such small peculiar velocities. For these models, we measure velocity dispersions between 500-700km/s and 150-300km/s respectively, more than twice the observed value. Although both CDM models fail to produce environments similar to those of our Local Group on a scale of a few Mpc, they can give rise to many binary systems that have similar orbital properties as the Milky Way--Andromeda system. The local, gravitationally induced bias of halos in the CDM ``Local Group'' environment, if defined within a sphere of 10 Mpc around each Local Group is about 1.5, independent of Omega. No biasing scheme could reconcile the measured velocity dispersions around Local Groups with the observed one. Identification of binary systems using a halo finder (named Skid (http://www-hpcc.astro.washington.edu/tools/DENMAX for a public version)) based on local density maxima instead of a simple linking algorithm, gives a much more complete sample. We show that a standard ``friend of friends'' algorithm would miss 40% of the LG candidates present in the simulations.

We test the precision with which weak lensing data can provide characteristic cluster mass profiles within Cold Dark Matter scenarios. Using a parallel treecode to simulate volumes as large as 500h-1 Mpc with good resolution, we generate samples of large clusters within a standard CDM model and an open CDM model with Omega_o=0.3. We mock high-quality lensing data by including realistic errors, selecting cluster samples based on velocity dispersion, and fitting profiles within a realistic range in radius. We find that a sample of ten clusters can determine logarithmic profile slopes with 1 sigma errors of about 7%. Increasing the sample size to twenty brings this error down to less than 5%, but this is still insufficient to distinguish the two models. However, measures of cluster profiles obtained with weak lensing do place strong constraints for general CDM--like models of structure formation, and we discuss the optimal strategy for obtaining data samples to use for this purpose.

It has been argued that a UV photoionizing background radiation field suppresses the formation of dwarf galaxies, and may even inhibit the formation of larger galaxies. In order to test this, we present gas-dynamical simulations of the formation of small objects in a CDM universe with and without a photoionizing background. The objects are selected from a collisionless simulation at a redshift of 2.4, and rerun at higher resolution including the effects of gas dynamics and using a hierarchical grid of particles. Five objects, each with a circular speed of 46 km/sec are simulated. The presence of the photoionizing background has only a small effect on the amount of gas that collapses in these objects, reducing the amount of cold collapsed gas by at most 30%. Analysis of the smaller objects found in the higher resolution simulation indicates that the photoionizing background only significantly affects the formation of objects with a virialized halo mass less than 10^9 soalr masses and circular speeds less than 23 km/sec. However, the ionization balance is greatly changed by the presence of the background radiation field. Typical lines of sight through the objects have 4 orders of magnitude less neutral hydrogen column density when the photoionizing background is included.

We investigate the injection spectrum of ultra-high-energy (>10^{15} eV) cosmic rays under the hypotheses that (1) these cosmic rays are protons and (2) the sources of these cosmic rays are extra-galactic and are homogeneously distributed in space. The most puzzling aspect of the observed ultra-high-energy cosmic ray spectrum is the apparent nonexistence of a ``Greisen cut-off'' at about 10^{19.8} eV. We show that this fact could be naturally explained if most (or all) of the cosmic rays presently observed above about 10^{19.6} eV were initially injected with energy above the Greisen cut-off. However, we find that the injection of cosmic rays above the Greisen cut-off cannot account for the observed flux below about 10^{19.6} eV unless the injection rate of these particles was enormously higher in the past, as would be the case if the injection resulted from the decay of an ultra-massive particle with lifetime of order 10^{9} yr. Even with such a rapid source evolution, the observed cosmic ray spectrum below about 10^{18.5} eV cannot be explained in this manner. However, we show that a 1/E^3 injection spectrum can account for the observed spectrum below 10^{18.5} eV}, with the steepening observed by the Fly's Eye group between 10^{17.6} eV and 10^{18.5} eV being very naturally explained by e+ - e- production effects. This latter fact lends support to the hypothesis that the cosmic rays in this energy regime are protons. However, due to e+ - e- production effects, a 1/E^3 injection spectrum cannot account for the observed flux above about 10^{18.5} eV.