A key uncertainty in galaxy evolution is the physics regulating star formation, ranging from small-scale processes related to the life-cycle of molecular clouds within galaxies to large-scale processes such as gas accretion onto galaxies. We study the imprint of such processes on the time-variability of star formation with an analytical approach tracking the gas mass of galaxies ("regulator model"). Specifically, we quantify the strength of the fluctuation in the star-formation rate (SFR) on different timescales, i.e. the power spectral density (PSD) of the star-formation history, and connect it to gas inflow and the life-cycle of molecular clouds. We show that in the general case the PSD of the SFR has three breaks, corresponding to the correlation time of the inflow rate, the equilibrium timescale of the gas reservoir of the galaxy, and the average lifetime of individual molecular clouds. On long and intermediate timescales (relative to the dynamical timescale of the galaxy), the PSD is typically set by the variability of the inflow rate and the interplay between outflows and gas depletion. On short timescales, the PSD shows an additional component related to the life-cycle of molecular clouds, which can be described by a damped random walk with a power-law slope of $\beta\approx2$ at high frequencies with a break near the average cloud lifetime. We discuss star-formation "burstiness" in a wide range of galaxy regimes, study the evolution of galaxies about the main sequence ridgeline, and explore the applicability of our method for understanding the star-formation process on cloud-scale from galaxy-integrated measurements.

We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

Physical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation of the IllustrisTNG project to explore the physical parameter space comprised of 8 resolved ISM properties in star-forming regions across a wide range of host galaxy mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$). We deconstruct our simulated galaxies into regions of 1 kpc in size, and measure the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region, and study them in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate weighted distributions of these properties, we show that stars primarily form in two distinct regimes, which are brought about by an underlying bi-component radial star formation rate profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time, and compare our findings with observations from the MaNGA IFU survey. Further, using principal component analysis, we characterise the aforementioned parameter space to reveal a high-degree of multicollinearity in relationships between ISM properties that drive the distribution of star formation at kiloparsec scales. Based on this, we find that a reduced three dimensional representation underpinned in essence by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.

In the course of doing astronomy, one often encounters plots of densities, for example probability densities, flux densities, and mass functions. Quite frequently the ordinate of these diagrams is plotted logarithmically to accommodate a large dynamic range. In this situation, I argue that it is critical to adjust the density appropriately, rather than simply setting the x-scale to `log' in your favorite plotting code. I will demonstrate the basic issue with a pedagogical example, then mention a few common plots where this may arise, and finally some possible exceptions to the rule.

We explore what may be learned by close encounters between extrasolar minor bodies like `Oumuamua and the Sun. These encounters may yield strong constraints on the bulk composition and possible origin of `Oumuamua-like objects. We find that such objects collide with the Sun once every 30 years, while about 2 pass within the orbit of Mercury each year. We identify preferred orientations for the orbits of extrasolar objects and point out known Solar System bodies with these orientations. We conclude using a simple Bayesian analysis that about one of these objects is extrasolar in origin, even if we cannot tell which.

Almost by definition brown dwarfs are objects with masses below the hydrogen burning limit, around $0.07\ M_\odot$. Below this mass, objects never reach a steady state where they can fuse hydrogen. Here we demonstrate, in contrast to this traditional view, that brown dwarfs with masses greater than the hydrogen burning limit may in principle exist in the universe. These objects, which we term "overmassive brown dwarfs" form a continuous sequence with traditional brown dwarfs in any property (mass, effective temperature, radius, luminosity). To form an overmassive brown dwarf, mass must be added sufficiently slowly to a sufficiently old traditional brown dwarf below the hydrogen burning limit. We identify mass transfer in binary brown dwarf systems via Roche lobe overflow driven by gravitational waves to be the most plausible mechanism to produce the bulk of the putative overmassive brown dwarf population.

Cold clouds in hot media are quickly crushed, shredded, and then accelerated as a result of their interaction with the background gas. The persistence of cold clouds moving at substantial velocities in harsh environments is a common yet puzzling feature of many astrophysical systems, from quasar absorption lines probing galactic halos to clouds of dust passing near Sgr $A^*$. Here we run a set of idealized numerical experiments, subjecting a line of cold clouds at a series of mutual separations to a hot background wind. We find that this stream of clouds is able to shield itself from hydrodynamic destruction by accelerating the hot background material, creating a protective layer of co-moving gas. We write down a simple diffusion equation that reproduces the behavior of the simulations, and we discuss the implications for cosmological gas accretion and G2.

We present a flexible, detailed model for the evolution of galactic discs in a cosmological context since $z\approx 4$, including a physically-motivated model for radial transport of gas and stars within galactic discs. This expansion beyond traditional semi-analytic models that do not include radial structure, or include only a prescribed radial structure, enables us to study the internal structure of disc galaxies and the processes that drive it. In order to efficiently explore the large parameter space allowed by this model, we construct a neural network-based emulator that can quickly return a reasonable approximation for many observables we can extract from the model, e.g. the star formation rate or the half mass stellar radius, at different redshifts. We employ the emulator to constrain the model parameters with Bayesian inference by comparing its predictions to 11 observed galaxy scaling relations at a variety of redshifts. The constrained models agree well with observations, both those used to fit the data and those not included in the fitting procedure. These models will be useful theoretical tools for understanding the increasingly detailed observational datasets from IFUs.

Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties.

We introduce a new model for the structure and evolution of the gas in galactic discs. In the model the gas is in vertical pressure and energy balance. Star formation feedback injects energy and momentum, and non-axisymmetric torques prevent the gas from becoming more than marginally gravitationally unstable. From these assumptions we derive the relationship between galaxies' bulk properties (gas surface density, stellar content, and rotation curve) and their star formation rates, gas velocity dispersions, and rates of radial inflow. We show that the turbulence in discs can be powered primarily by star formation feedback, radial transport, or a combination of the two. In contrast to models that omit either radial transport or star formation feedback, the predictions of this model yield excellent agreement with a wide range of observations, including the star formation law measured in both spatially resolved and unresolved data, the correlation between galaxies' star formation rates and velocity dispersions, and observed rates of radial inflow. The agreement holds across a wide range of galaxy mass and type, from local dwarfs to extreme starbursts to high-redshifts discs. We apply the model to galaxies on the star-forming main sequence, and show that it predicts a transition from mostly gravity-driven turbulence at high redshift to star formation-driven turbulence at low redshift. This transition, and the changes in mass transport rates that it produces, naturally explain why galaxy bulges tend to form at high redshift and discs at lower redshift, and why galaxies tend to quench inside-out.

Black holes growing via the accretion of gas emit radiation that can photoevaporate the atmospheres of nearby planets. Here we couple planetary structural evolution models of sub-Neptune mass planets to the growth of the Milky way's central supermassive black-hole, Sgr A$^*$ and investigate how planetary evolution is influenced by quasar activity. We find that, out to ${\sim} 20$ pc from Sgr A$^*$, the XUV flux emitted during its quasar phase can remove several percent of a planet's H/He envelope by mass; in many cases, this removal results in bare rocky cores, many of which situated in the habitable zones (HZs) of G-type stars. The erosion of sub-Neptune sized planets may be one of the most prevalent channels by which terrestrial super-Earths are created near the Galactic Center. As such, the planet population demographics may be quite different close to Sgr A$^*$ than in the Galaxy's outskirts. The high stellar densities in this region (about seven orders of magnitude greater than the solar neighborhood) imply that the distance between neighboring rocky worlds is short ($500-5000$~AU). The proximity between potentially habitable terrestrial planets may enable the onset of widespread interstellar panspermia near the nuclei of galaxies. More generally, we predict these phenomena to be ubiquitous for planets in nuclear star clusters and ultra-compact dwarfs. Globular clusters, on the other hand, are less affected by the black holes.

Planetary atmospheres are subject to mass loss through a variety of mechanisms including irradiation by XUV photons from their host star. Here we explore the consequences of XUV irradiation by supermassive black holes as they grow by the accretion of gas in galactic nuclei. Based on the mass distribution of stars in galactic bulges and disks and the luminosity history of individual black holes, we estimate the probability distribution function of XUV fluences as a function of galaxy halo mass, redshift, and stellar component. We find that about 50% of all planets in the universe may lose the equivalent of a Martian atmosphere, 10% may lose an Earth's atmosphere, and 0.2% may lose the mass of Earth's oceans. The fractions are appreciably higher in the spheroidal components of galaxies, and depend strongly on galaxy mass, but only weakly on redshift.

Self-gravity and stellar feedback are capable of driving turbulence and transporting mass and angular momentum in disk galaxies, but the balance between them is not well understood. In the previous paper in this series, we showed that gravity alone can drive turbulence in galactic disks, regulate their Toomre $Q$ parameters to $\sim$ 1, and transport mass inwards at a rate sufficient to fuel star formation in the centers of present-day galaxies. In this paper we extend our models to include the effects of star formation feedback. We show that feedback suppresses galaxies' star formation rates by a factor of $\sim$ 5 and leads to the formation of a multi-phase atomic and molecular ISM. Both the star formation rate and the phase balance produced in our simulations agree well with observations of nearby spirals. After our galaxies reach steady state, we find that the inclusion of feedback actually lowers the gas velocity dispersion slightly compared to the case of pure self-gravity, and also slightly reduces the rate of inward mass transport. Nevertheless, we find that, even with feedback included, our galactic disks self-regulate to $Q$ $\sim$ 1, and transport mass inwards at a rate sufficient to supply a substantial fraction of the inner disk star formation. We argue that gravitational instability is therefore likely to be the dominant source of turbulence and transport in galactic disks, and that it is responsible for fueling star formation in the inner parts of galactic disks over cosmological times.

Photoelectric heating has long been recognized as the primary source of heating for the neutral interstellar medium. Simulations of spiral galaxies found some indication that photoelectric heating could suppress star formation. However, simulations that include photoelectric heating have typically found that it has little effect on the rate of star formation in either spiral galaxies or dwarfs suggesting that supernovae and not photoelectric heating are responsible for setting the star formation law in galaxies. This result is in tension with recent work indicating that a star formation law that depends on galaxy metallicity, as expected for photoelectric heating but not for supernovae, reproduces the present-day galaxy population better than a metallicity-independent one. Here we report a series of simulations of dwarf galaxies, where the effects of both photoelectric heating and supernovae are expected to be strongest. We simultaneously include space- and time-dependent photoelectric heating, and we resolve the Sedov phase of every supernova blast wave, allowing us to make a direct measurement of the relative importance of momentum injection by supernovae and dust heating by far ultraviolet (FUV) photons in suppressing star formation. We find that supernovae are unable to account for the long observed gas depletion times in dwarf galaxies. Instead, ordinary photoelectric heating is the dominant means by which dwarf galaxies regulate their star formation rate at any given time, suppressing the star formation rate by more than an order of magnitude relative to simulations with only supernovae.

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

The role of gravitational instability-driven turbulence in determining the structure and evolution of disk galaxies, and the extent to which gravity rather than feedback can explain galaxy properties, remains an open question. To address it, we present high resolution adaptive mesh refinement simulations of Milky Way-like isolated disk galaxies, including realistic heating and cooling rates and a physically motivated prescription for star formation, but no form of star formation feedback. After an initial transient, our galaxies reach a state of fully-nonlinear gravitational instability. In this state, gravity drives turbulence and radial inflow. Despite the lack of feedback, the gas in our galaxy models shows substantial turbulent velocity dispersions, indicating that gravitational instability alone may be able to power the velocity dispersions observed in nearby disk galaxies on 100 pc scales. Moreover, the rate of mass transport produced by this turbulence approaches $\sim 1$ $M_\odot$ yr$^{-1}$ for Milky Way-like conditions, sufficient to fully fuel star formation in the inner disks of galaxies. In a companion paper we add feedback to our models, and use the comparison between the two cases to understand what galaxy properties depend sensitively on feedback, and which can be understood as the product of gravity alone. All of the code, initial conditions, and simulation data for our model are publicly available.

Spinning planetesimals in a gaseous protoplanetary disk may experience a hydrodynamical force perpendicular to their relative velocities. We examine the effect this force has on the dynamics of these objects using analytical arguments based on a simple laminar disk model and numerical integrations of the equations of motion for individual grains. We focus in particular on meter-sized boulders traditionally expected to spiral in to the central star in as little as 100 years from 1 A.U. We find that there are plausible scenarios in which this force extends the lifetime of these solids in the disk by a factor of several. More importantly the velocities induced by the Magnus force can prevent the formation of planetesimals via gravitational instability in the inner disk if the size of the dust particles is larger than of order 10 cm. We find that the fastest growing linear modes of the streaming instability may still grow despite the diffusive effect of the Magnus force, but it remains to be seen how the Magnus force will alter the non-linear evolution of these instabilities.

Metal production in galaxies traces star formation, and is highly concentrated toward the centers of galactic discs. This suggests that galaxies should have inhomogeneous metal distributions with strong radial gradients, but observations of present-day galaxies show only shallow gradients with little azimuthal variation, implying the existence of a redistribution mechanism. We study the role of gravitational instability-driven turbulence as a mixing mechanism by simulating an isolated galactic disc at high resolution, including metal fields treated as passive scalars. Since any cylindrical field can be decomposed into a sum of Fourier-Bessel basis functions, we set up initial metal fields characterized by these functions and study how different modes mix. We find both shear and turbulence contribute to mixing, but the mixing strongly depends on the symmetries of the mode. Non-axisymmetric modes have decay times smaller than the galactic orbital period because shear winds them up to small spatial scales, where they are erased by turbulence. The decay timescales for axisymmetric modes are much greater, though for all but the largest-scale inhomogeneities the mixing timescale is still short enough to erase chemical inhomogeneities over cosmological times. These different timescales provide an explanation for why galaxies retain metallicity gradients while there is almost no variation at a fixed radius. Moreover, the comparatively long timescales required for mixing axisymmetric modes may explain the greater diversity of metallicity gradients observed in high redshift galaxies as compared to local ones: these systems have not yet reached equilibrium between metal production and diffusion.

The evolution of thin axisymmetric viscous accretion disks is a classic problem in astrophysics. While models based on this simplified geometry provide only approximations to the true processes of instability-driven mass and angular momentum transport, their simplicity makes them invaluable tools for both semi-analytic modeling and simulations of long-term evolution where two- or three-dimensional calculations are too computationally costly. Despite the utility of these models, the only publicly-available frameworks for simulating them are rather specialized and non-general. Here we describe a highly flexible, general numerical method for simulating viscous thin disks with arbitrary rotation curves, viscosities, boundary conditions, grid spacings, equations of state, and rates of gain or loss of mass (e.g., through winds) and energy (e.g., through radiation). Our method is based on a conservative, finite-volume, second-order accurate discretization of the equations, which we solve using an unconditionally-stable implicit scheme. We implement Anderson acceleration to speed convergence of the scheme, and show that this leads to factor of $\sim 5$ speed gains over non-accelerated methods in realistic problems, though the amount of speedup is highly problem-dependent. We have implemented our method in the new code Viscous Accretion Disk Evolution Resource (VADER), which is freely available for download from https://bitbucket.org/krumholz/vader/ under the terms of the GNU General Public License.

We present a simple toy model to understand what sets the scatter in star formation and metallicity of galaxies at fixed mass. The scatter ultimately arises from the intrinsic scatter in the accretion rate, but may be substantially reduced depending on the timescale on which the accretion varies compared to the timescale on which the galaxy loses gas mass. This model naturally produces an anti-correlation between star formation and metallicity at a fixed mass, the basis of the fundamental metallicity relation. We show that observational constraints on the scatter in galaxy scaling relations can be translated into constraints on the galaxy-to-galaxy variation in the mass loading factor, and the timescales and magnitude of stochastic accretion onto star-forming galaxies. We find a remarkably small scatter in the mass loading factor, < 0.1 dex, and that the scatter in accretion rates is smaller than expected from N-body simulations.

Over the past 10 Gyr, star-forming galaxies have changed dramatically, from clumpy and gas rich, to rather quiescent stellar-dominated disks with specific star formation rates lower by factors of a few tens. We present a general theoretical model for how this transition occurs, and what physical processes drive it, making use of 1D axisymmetric thin disk simulations with an improved version of the Gravitational Instability-Dominated Galaxy Evolution Tool (GIDGET) code. We show that at every radius galaxies tend to be in a slowly evolving equilibrium state wherein new accretion is balanced by star formation, galactic winds, and radial transport of gas through the disk by gravitational instability (GI) -driven torques. The gas surface density profile is determined by which of these terms are in balance at a given radius - direct accretion is balanced by star formation and galactic winds near galactic centers, and by transport at larger radii. We predict that galaxies undergo a smooth transition from a violent disk instability phase to secular evolution. This model provides a natural explanation for the high velocity dispersions and large clumps in z~2 galaxies, the growth and subsequent quenching of bulges, and features of the neutral gas profiles of local spiral galaxies.

Observations of disk galaxies at z~2 have demonstrated that turbulence driven by gravitational instability can dominate the energetics of the disk. We present a 1D simulation code, which we have made publicly available, that economically evolves these galaxies from z~2 to z~0 on a single CPU in a matter of minutes, tracking column density, metallicity, and velocity dispersions of gaseous and multiple stellar components. We include an H$_2$ regulated star formation law and the effects of stellar heating by transient spiral structure. We use this code to demonstrate a possible explanation for the existence of a thin and thick disk stellar population and the age-velocity dispersion correlation of stars in the solar neighborhood: the high velocity dispersion of gas in disks at z~2 decreases along with the cosmological accretion rate, while at lower redshift, the dynamically colder gas forms the low velocity dispersion stars of the thin disk.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

Cosmological constraints from the cosmic microwave background (CMB) anisotropies rely on accurate theoretical calculations of the cosmic recombination history. Recent work has emphasized the importance of radiative transfer calculations due to the high optical depth in the HI Lyman lines. Transfer in the Lyman-alpha line is dominated by true emission and absorption, Hubble expansion, and resonant scattering. Resonant scattering causes photons to diffuse in frequency due to random kicks from the thermal velocities of hydrogen atoms, and also to drift toward lower frequencies due to energy loss via atomic recoil. Past analyses of Lyman-alpha transfer during the recombination era have either considered a subset of these processes, ignored time dependence, or incorrectly assumed identical emission and absorption profiles. We present here a fully time-dependent radiative transfer calculation of the Lyman-alpha line including all of these processes, and compare it to previous results that ignored the resonant scattering. We find a faster recombination due to recoil enhancement of the Lyman-alpha escape rate, leading to a reduction in the free electron density of 0.45% at z=900. This results in an increase in the small-scale CMB power spectrum that is negligible for the current data but will be a 0.9 sigma correction for Planck. We discuss the reasons why we find a smaller correction than some other recent computations.

We present first evidence for the so-called Head-Tail asymmetry signature of neutron-induced nuclear recoil tracks at energies down to 1.5 keV/amu using the 1m^3 DRIFT-IIc dark matter detector. This regime is appropriate for recoils induced by Weakly Interacting Massive Particle (WIMPs) but one where the differential ionization is poorly understood. We show that the distribution of recoil energies and directions induced here by Cf-252 neutrons matches well that expected from massive WIMPs. The results open a powerful new means of searching for a galactic signature from WIMPs.

The DRIFT collaboration utilizes low pressure gaseous detectors to search for WIMP dark matter with directional signatures. A 252Cf neutron source was placed on each of the principal axes of a DRIFT detector in order to test its ability to measure directional signatures from the three components of very low energy (~keV/amu) recoil ranges. A high trigger threshold and the event selection procedure ensured that only sulfur recoils were analyzed. Sulfur recoils produced in the CS2 target gas by the 252Cf source closely match those expected from massive WIMP induced sulfur recoils. For each orientation of the source a directional signal from the range components was observed, indicating that the detector is directional along all 3 axes. An analysis of these results yields an optimal orientation for DRIFT detectors when searching for a directional signature from WIMPs. Additional energy dependent information is provided to aid in understanding this effect.