The magnetised solar wind modulates the Galactic cosmic ray flux in the heliosphere up to rigidities as high as 40 GeV. In this work, we present a new and straightforward extension of the popular, but limited force-field model, thus providing a fast and robust method for phenomenological studies of Galactic cosmic rays. Our semi-analytical approach takes into account charge-sign dependent modulation due to drifts in the heliospheric magnetic field and has been validated via comparison to a fully numerical code. Our model nicely reproduces the time-dependent AMS-02 measurements and we find the strength of diffusion and drifts to be strongly correlated with the heliospheric tilt angle and magnitude of the magnetic field. We are able to predict the electron and positron fluxes beyond the range for which measurements by AMS-02 have been presented. We have made an example script for the semi-analytical model publicly available and we urge the community to adopt this approach for phenomenological studies.

In recent years, it has become possible to detect individual dark matter subhalos near images of strongly lensed extended background galaxies. Typically, only the most massive subhalos in the strong lensing region may be detected this way. In this work, we show that strong lenses may also be used to constrain the much more numerous population of lower mass subhalos that are too small to be detected individually. In particular, we show that the power spectrum of projected density fluctuations in galaxy halos can be measured using strong gravitational lensing. We develop the mathematical framework of power spectrum estimation, and test our method on mock observations. We use our results to determine the types of observations required to measure the substructure power spectrum with high significance. We predict that deep observations ($\sim10$ hours on a single target) with current facilities can measure this power spectrum at the $3\sigma$ level, with no apparent degeneracy with unknown clumpiness in the background source structure or fluctuations from detector noise. Upcoming ALMA measurements of strong lenses are capable of placing strong constraints on the abundance of dark matter subhalos and the underlying particle nature of dark matter.

We present an analysis of the effects of dissipational baryonic physics on the local dark matter (DM) distribution at the location of the Sun, with an emphasis on the consequences for direct detection experiments. Our work is based on a comparative analysis of two cosmological simulations of a Milky Way halo: Eris, a full hydrodynamic simulation, and ErisDark, its DM-only realization. We find that two distinct processes lead in Eris to a 30% enhancement of DM in the disk plane at the location of the Sun: the accretion and disruption of satellites resulting in a DM component with net angular momentum and the contraction of baryons pulling DM into the disk plane without forcing it to co-rotate. Owing to its particularly quiescent merger history for dark halos of Milky Way mass, the co-rotating dark disk in Eris is less massive than what has been suggested by previous work, contributing only 9% of the local DM density. The speed distribution in Eris is broadened and shifted to higher speeds compared to its DM-only twin simulation ErisDark. At high speeds $f(v)$ falls more steeply in Eris than in ErisDark or the SHM, easing the tension between recent results from the CDMS-II and XENON100 experiments. The non-Maxwellian aspects of $f(v)$ are still present, but much less pronounced in Eris than in DM-only runs. The weak dark disk increases the time-averaged scattering rate by only a few percent at low recoil energies. On the high velocity tail, however, the increase in typical speeds due to baryonic contraction results in strongly enhanced mean scattering rates compared to ErisDark, although they are still suppressed compared to the SHM. Similar trends are seen regarding the amplitude of the annual modulation, while the modulated fraction is increased compared to the SHM and decreased compared to ErisDark.

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.

We present a cosmological hydrodynamic simulation of the formation of dwarf galaxies at redshifts z>~2.5 using a physically-motivated model for H2-regulated star formation. Our simulation, performed using the Enzo code and reaching a peak resolution of 109 proper parsecs at z=2.5, extends the results of Kuhlen et al. (2012) to significantly lower redshifts. We show that a star formation prescription regulated by the local H2 abundance leads to the suppression of star formation in dwarf galaxy halos with M_h <~ 10^10 Msun and to a large population of gas-rich "dark galaxies" at z=2.5 with low star formation efficiencies and gas depletion timescales >20 Gyr. The fraction of dark galaxies is 60% at M_h ~ 10^10 Msun and increases rapidly with decreasing halo mass. Dark galaxies form late and their gaseous disks never reach the surface densities, > ~5700 Msun / pc^2 (Z/10^-3 Zsun)^(-0.88), that are required to build a substantial molecular fraction. Despite this large population of dark galaxies, we show that our H2-regulated simulation is consistent with both the observed luminosity function of galaxies and the cosmological mass density of neutral gas at z>~2.5. Moreover, our results provide a theoretical explanation for the recent detection in fluorescent Ly-alpha emission of gaseous systems at high redshift with little or no associated star formation. We further propose that H2-regulation may offer a fresh solution to a number of outstanding "dwarf galaxy problems" in LambdaCDM. In particular, H2-regulation leads galaxy formation to become effectively stochastic on mass scales of M_h ~ 10^10 Msun, and thus these massive dwarfs are not "too big to fail".

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way-mass galaxies. Zolotov et al. (2012) have suggested a correction to be applied to the central masses of dark matter-only satellites in order to mimic the affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites, and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this correction to the z=0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we find that the number of massive, luminous satellites around a Milky Way-mass galaxy is in agreement with the number of observed satellites around the Milky Way or M31. We conclude that baryonic processes have the potential to solve the missing satellites problem.

We show that the position of the central dark matter density peak may be expected to differ from the dynamical center of the Galaxy by several hundred parsec. In Eris, a high resolution cosmological hydrodynamics simulation of a realistic Milky-Way-analog disk galaxy, this offset is 300 - 400 pc (~3 gravitational softening lengths) after z=1. In its dissipationless dark-matter-only twin simulation ErisDark, as well as in the Via Lactea II and GHalo simulations, the offset remains below one softening length for most of its evolution. The growth of the DM offset coincides with a flattening of the central DM density profile in Eris inwards of ~1 kpc, and the direction from the dynamical center to the point of maximum DM density is correlated with the orientation of the stellar bar, suggesting a bar-halo interaction as a possible explanation. A dark matter density offset of several hundred parsec greatly affects expectations of the dark matter annihilation signals from the Galactic Center. It may also support a dark matter annihilation interpretation of recent reports by Weniger (2012) and Su & Finkbeiner (2012) of highly significant 130 GeV gamma-ray line emission from a region 1.5 degrees (~200 parsec projected) away from Sgr A* in the Galactic plane.

We present a review of the current state of the art of cosmological dark matter simulations, with particular emphasis on the implications for dark matter detection efforts and studies of dark energy. This review is intended both for particle physicists, who may find the cosmological simulation literature opaque or confusing, and for astro-physicists, who may not be familiar with the role of simulations for observational and experimental probes of dark matter and dark energy. Our work is complementary to the contribution by M. Baldi in this issue, which focuses on the treatment of dark energy and cosmic acceleration in dedicated N-body simulations. Truly massive dark matter-only simulations are being conducted on national supercomputing centers, employing from several billion to over half a trillion particles to simulate the formation and evolution of cosmologically representative volumes (cosmic scale) or to zoom in on individual halos (cluster and galactic scale). These simulations cost millions of core-hours, require tens to hundreds of terabytes of memory, and use up to petabytes of disk storage. The field is quite internationally diverse, with top simulations having been run in China, France, Germany, Korea, Spain, and the USA. Predictions from such simulations touch on almost every aspect of dark matter and dark energy studies, and we give a comprehensive overview of this connection. We also discuss the limitations of the cold and collisionless DM-only approach, and describe in some detail efforts to include different particle physics as well as baryonic physics in cosmological galaxy formation simulations, including a discussion of recent results highlighting how the distribution of dark matter in halos may be altered. We end with an outlook for the next decade, presenting our view of how the field can be expected to progress. (abridged)

We investigate how strong lensing of dusty, star-forming galaxies by foreground galaxies can be used as a probe of dark matter halo substructure. We find that spatially resolved spectroscopy of lensed sources allows dramatic improvements to measurements of lens parameters. In particular we find that modeling of the full, three-dimensional (angular position and radial velocity) data can significantly facilitate substructure detection, increasing the sensitivity of observables to lower mass subhalos. We carry out simulations of lensed dusty sources observed by early ALMA (Cycle 1) and use a Fisher matrix analysis to study the parameter degeneracies and mass detection limits of this method. We find that, even with conservative assumptions, it is possible to detect galactic dark matter subhalos of ~ 10^8 M_{\odot} with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG lenses, there is a ~ 55 % probability of detecting a substructure with M>10^8 M_{\odot} with more than 5 sigma detection significance in each lens, if the abundance of substructure is consistent with previous lensing results. The full ALMA array, with its significantly enhanced sensitivity and resolution, should improve these estimates considerably. Given the sample of ~100 lenses provided by surveys like the South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised to improve dramatically over the next few years.

Tidal stripping of dark matter from subhalos falling into the Milky Way produces narrow, cold tidal streams as well as more spatially extended "debris flows" in the form of shells, sheets, and plumes. Here we focus on the debris flow in the Via Lactea II simulation, and show that this incompletely phase-mixed material exhibits distinctive high velocity behavior. Unlike tidal streams, which may not necessarily intersect the Earth's location, debris flow is spatially uniform at 8 kpc and thus guaranteed to be present in the dark matter flux incident on direct detection experiments. At Earth-frame speeds greater than 450 km/s, debris flow comprises more than half of the dark matter at the Sun's location, and up to 80% at even higher speeds. Therefore, debris flow is most important for experiments that are particularly sensitive to the high speed tail of the dark matter distribution, such as searches for light or inelastic dark matter or experiments with directional sensitivity. We show that debris flow yields a distinctive recoil energy spectrum and a broadening of the distribution of incidence direction.

We distinguish between Local Group field galaxies which may have passed through the virial volume of the Milky Way, and those which have not, via a statistical compari- son against populations of dark matter haloes in the Via Lactea II (VLII) simulation with known orbital histories. Analysis of VLII provides expectations for this escaped population: they contribute 13 per cent of the galactic population between 300 and 1500 kpc from the Milky Way, and hence we anticipate that about 7 of the 54 known Local Group galaxies in that distance range are likely to be Milky Way escapees. These objects can be of any mass below that of the Milky Way, and they are expected to have positive radial velocities with respect to the Milky Way. Comparison of the radius-velocity distributions of VLII populations and measurements of Local Group galaxies presents a strong likelihood that Tucana, Cetus, NGC3109, SextansA, SextansB, Antlia, NGC6822, Phoenix, LeoT, and NGC185 have passed through the Milky Way. Most of these dwarfs have a lower HI mass fraction than the majority of dwarfs lying at similar distances to either the Milky Way or M31. Indeed, several of these galaxies - especially those with lower masses - contain signatures in their morphology, star formation history and/or gas content indicative of evolution seen in simulations of satellite/parent galactic interactions. Our results offer strong support for scenarios in which dwarfs of different types form a sequence in morphology and gas content, with evolution along the sequence being driven by interaction history.

We describe cosmological galaxy formation simulations with the adaptive mesh refinement code Enzo that incorporate a star formation prescription regulated by the local abundance of molecular hydrogen. We show that this H2-regulated prescription leads to a suppression of star formation in low mass halos (M_h < ~10^10 M_sun) at z>4, alleviating some of the dwarf galaxy problems faced by theoretical galaxy formation models. H2 regulation modifies the efficiency of star formation of cold gas directly, rather than indirectly reducing the cold gas content with "supernova feedback". We determine the local H2 abundance in our most refined grid cells (76 proper parsec in size at z=4) by applying the model of Krumholz, McKee, & Tumlinson, which is based on idealized 1D radiative transfer calculations of H2 formation-dissociation balance in ~100 pc atomic--molecular complexes. Our H2-regulated simulations are able to reproduce the empirical (albeit lower z) Kennicutt-Schmidt relation, including the low Sigma_gas cutoff due to the transition from atomic to molecular phase and the metallicity dependence thereof, without the use of an explicit density threshold in our star formation prescription. We compare the evolution of the luminosity function, stellar mass density, and star formation rate density from our simulations to recent observational determinations of the same at z=4-8 and find reasonable agreement between the two.

Recent observations have constrained the galaxy UV luminosity function up to z~10. However, these observations alone allow for a wide range of reionization scenarios due to uncertainties in the abundance of faint galaxies and the escape fraction of ionizing photons. We show that requiring continuity with post-reionization (z<6) measurements, where the Lya forest provides a complete probe of the cosmological emissivity of ionizing photons, significantly reduces the permitted parameter space. Models that are simultaneously consistent with the measured UV luminosity function, the Thomson optical depth to the CMB, and the Lya forest data require either: 1) extrapolation of the galaxy luminosity function down to very faint UV magnitudes M_lim ~ -10, corresponding roughly to the UV background suppression scale; 2) an increase of f_esc by a factor > ~10 from z=4 (where the best fit is 4%) to z=9; or 3) more likely, a hybrid solution in which undetected galaxies contribute significantly and f_esc increases more modestly. Models in which star formation is strongly suppressed in low-mass, reionization-epoch haloes of mass up to ~10^10 M_sun (e.g., owing to a metallicity dependence) are only allowed for extreme assumptions for the evolution of f_esc. However, variants of such models in which the suppression mass is reduced (e.g., assuming an earlier or higher metallicity floor) are in better agreement with the data. Concordance scenarios satisfying the available data predict a consistent redshift of 50% ionized fraction z_reion(50%) ~ 10. On the other hand, the duration of reionization is sensitive to the relative contribution of bright versus faint galaxies, with scenarios dominated by faint galaxies predicting a more extended reionization event. Scenarios relying too heavily on high-z dwarfs are disfavored by kinetic Sunyaev-Zeldovich measurements, which prefer a short reionization history.

We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with stars. The method is calibrated using the observed luminosity function of Milky Way satellites and the concentration of their stellar populations, and self-consistently follows the accretion and disruption of progenitor dwarfs and the build-up of the stellar halo in a cosmological "live host". Simple prescriptions for assigning stellar populations to collisionless particles are able to reproduce many properties of the observed Milky Way halo and its surviving dwarf satellites, like velocity dispersions, sizes, brightness profiles, metallicities, and spatial distribution. Our model predicts the existence of approximately 1,850 subhalos harboring "extremely faint" satellites (with mass-to-light ratios >5,000) lying beyond the Sloan Digital Sky Survey detection threshold. Of these, about 20 are "first galaxies", i.e. satellites that formed a stellar mass above 10 Msun before redshift 9. The ten most luminous satellites (L> 1e6 Lsun) in the simulation are hosted by subhalos with peak circular velocities today in the range V_max=10-40 km/s that have shed between 80% and 99% of their dark mass after being accreted at redshifts 1.7< z <4.6. The satellite maximum circular velocity and stellar line-of-sight velocity dispersion today follow the relation V_max=2.2 sigma_los. We apply a standard mass estimation algorithm based on Jeans modelling of the line-of-sight velocity dispersion profiles to the simulated dwarf spheroidals, and test the accuracy of this technique. The inner (within 300 pc) mass-luminosity relation for currently detectable satellites is nearly flat in our model, in qualitative agreement with the "common mass scale" found in Milky Way dwarfs. We do, however, predict a weak, but significant positive correlation for these objects: M_300 ~L^{0.088 \pm 0.024}.

A longstanding puzzle of fundamental importance in modern cosmology has been the origin of the nearly universal density profiles of dark matter halos found in N-body simulations -- the so-called NFW profile. We show how this behavior may be understood, simply, by applying adiabatic contraction to peaks of Gaussian random fields. We argue that dynamical friction acts to reduce enormously the effect of random scatter in the properties of initial peaks, providing a key simplification. We compare our model predictions with results of the ultra-high resolution Via Lactea-II N-body simulation, and find superb agreement. We show how our model may be used to predict the distribution of halo properties like concentration. Our results suggest that many of the basic properties of halo structure may be understood using extremely simple physics.

We present a study of the ability of the Fermi Gamma-ray Space Telescope to detect dark-matter annihilation signals from the Galactic subhalos predicted by the Via Lactea II N-body simulation. We implement an improved formalism for estimating the boost factor needed to account for the effect of dark-matter clumping on scales below the resolution of the simulation, and we incorporate a detailed Monte Carlo simulation of the response of the Fermi-LAT telescope, including a simulation of its all-sky observing mode integrated over a ten year mission. We find that for WIMP masses up to about 150 GeV in standard supersymmetric models with velocity-averaged cross section 3*10^-26 cm^3 s^-1, a few subhalos could be detectable with >5 standard deviations significance and would likely deviate significantly from the appearance of a point source.

We study the effects of substructure on the rate of dark-matter annihilation in the Galactic halo. We use an analytic model for substructure that can extend numerical simulation results to scales too small to be resolved by the simulations. We first calibrate the analytic model to numerical simulations, and then determine the annihilation boost factor, for standard WIMP models as well as those with Sommerfeld (or other) enhancements, as a function of Galactocentric radius in the Milky Way. We provide an estimate of the dependence of the gamma-ray intensity of WIMP annihilation as a function of angular distance from the Galactic center. This methodology, coupled with future numerical simulation results can be a powerful tool that can be used to constrain WIMP properties using Fermi all-sky data.

The velocity distribution function of dark matter particles is expected to show significant departures from a Maxwell-Boltzmann distribution. This can have profound effects on the predicted dark matter - nucleon scattering rates in direct detection experiments, especially for dark matter models in which the scattering is sensitive to the high velocity tail of the distribution, such as inelastic dark matter (iDM) or light (few GeV) dark matter (LDM), and for experiments that require high energy recoil events, such as many directionally sensitive experiments. Here we determine the velocity distribution functions from two of the highest resolution numerical simulations of Galactic dark matter structure (Via Lactea II and GHALO), and study the effects for these scenarios. For directional detection, we find that the observed departures from Maxwell-Boltzmann increase the contrast of the signal and change the typical direction of incoming DM particles. For iDM, the expected signals at direct detection experiments are changed dramatically: the annual modulation can be enhanced by more than a factor two, and the relative rates of DAMA compared to CDMS can change by an order of magnitude, while those compared to CRESST can change by a factor of two. The spectrum of the signal can also change dramatically, with many features arising due to substructure. For LDM the spectral effects are smaller, but changes do arise that improve the compatibility with existing experiments. We find that the phase of the modulation can depend upon energy, which would help discriminate against background should it be found.

The coalescence of a massive black hole (MBH) binary leads to the gravitational-wave recoil of the system and its ejection from the galaxy core. We have carried out N-body simulations of the motion of a MBH = 3.7x10^6 Msun MBH remnant in the Via Lactea I simulation, a Milky Way sized dark matter halo. The black hole receives a recoil velocity of Vkick = 80, 120, 200, 300, and 400 km/s at redshift 1.5, and its orbit is followed for over 1 Gyr within a live host halo, subject only to gravity and dynamical friction against the dark matter background. We show that, owing to asphericities in the dark matter potential, the orbit of the MBH is hightly non-radial, resulting in a significantly increased decay timescale compared to a spherical halo. The simulations are used to construct a semi-analytic model of the motion of the MBH in a time-varying triaxial Navarro-Frenk-White dark matter halo plus a spherical stellar bulge, where the dynamical friction force is calculated directly from the velocity dispersion tensor. Such a model should offer a realistic picture of the dynamics of kicked MBHs in situations where gas drag, friction by disk stars, and the flattening of the central cusp by the returning black hole are all negligible effects. We find that MBHs ejected with initial recoil velocities Vkick > 500 km/s do not return to the host center within Hubble time. In a Milky Way-sized galaxy, a recoiling hole carrying a gaseous disk of initial mass ~MBH may shine as a quasar for a substantial fraction of its wandering phase. The long decay timescales of kicked MBHs predicted by this study may thus be favorable to the detection of off-nuclear quasar activity.

The unambiguous detection of Galactic dark matter annihilation would unravel one of the most outstanding puzzles in particle physics and cosmology. Recent observations have motivated models in which the annihilation rate is boosted by the Sommerfeld effect, a non-perturbative enhancement arising from a long range attractive force. Here we apply the Sommerfeld correction to Via Lactea II, a high resolution N-body simulation of a Milky-Way-size galaxy, to investigate the phase-space structure of the Galactic halo. We show that the annihilation luminosity from kinematically cold substructure can be enhanced by orders of magnitude relative to previous calculations, leading to the prediction of gamma-ray fluxes from up to hundreds of dark clumps that should be detectable by the Fermi satellite.

We study the local flux of electrons and positrons from annihilating Dark Matter (DM), and investigate how its spectrum depends on the choice of DM model and inhomogeneities in the DM distribution. Below a cutoff energy, the flux is expected to have a universal power-law form with an index n ~ -2. The cutoff energy and the behavior of the flux near the cutoff is model dependent. The dependence on the DM host halo profile may be significant at energies E < 100 GeV and leads to softening of the flux, n < -2. There may be additional features at high energies due to the presence of local clumps of DM, especially for models in which the Sommerfeld effect boosts subhalo luminosities. In general, the flux from a nearby clump gives rise to a harder spectrum of electrons and positrons, with an index n > -2. Using the Via Lactea II simulation, we estimate the probability of such subhalo effects in a generic Sommerfeld-enhanced model to be at least 4%, and possibly as high as 15% if subhalos below the simulation's resolution limit are accounted for. We discuss the consequences of these results for the interpretation of the ATIC, PAMELA, HESS, and Fermi data, as well as for future experiments.

Dark Matter annihilation holds great potential for directly probing the clumpiness of the Galactic halo that is one of the key predictions of the Cold Dark Matter paradigm of hierarchical structure formation. Here we review the gamma-ray signal arising from dark matter annihilation in the centers of Galactic subhalos. We consider both known Galactic dwarf satellite galaxies and dark clumps without a stellar component as potential sources. Utilizing the Via Lactea II numerical simulation, we estimate fluxes for 18 Galactic dwarf spheroidals with published central densities. The most promising source is Segue 1, followed by Ursa Major II, Ursa Minor, Draco, and Carina. We show that if any of the known Galactic satellites can be detected, then at least ten times more subhalos should be visible, with a significant fraction of them being dark clumps.

In order to attract and retain excellent researchers and diverse individuals in astrophysics, we recommend action be taken in several key areas impacting young scientists: (1) Maintain balance between large collaborations and individual projects through distribution of funding; encourage public releases of observational and simulation data for use by a broader community. (2) Improve the involvement of women, particularly at leading institutions. (3) Address the critical shortage of child care options and design reasonable profession-wide parental leave policies. (4) Streamline the job application and hiring process. We summarize our reasons for bringing these areas to the attention of the committee, and we suggest several practical steps that can be taken to address them.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The hierarchical theory of galaxy formation rests on the idea that smaller galactic structures merge to form the galaxies that we see today. The past decade has provided remarkable observational support for this scenario, driven in part by advances in spectroscopic instrumentation. Multi-object spectroscopy enabled the discovery of kinematically cold substructures around the Milky Way and M31 that are likely the debris of disrupting satellites. Improvements in high-resolution spectroscopy have produced key evidence that the abundance patterns of the Milky Way halo and its dwarf satellites can be explained by Galactic chemical evolution models based on hierarchical assembly. These breakthroughs have depended almost entirely on observations of nearby stars in the Milky Way and luminous red giant stars in M31 and Local Group dwarf satellites. In the next decade, extremely large telescopes will allow observations far down the luminosity function in the known dwarf galaxies, and they will enable observations of individual stars far out in the Galactic halo. The chemical abundance census now available for the Milky Way will become possible for our nearest neighbor, M31. Velocity dispersion measurements now available in M31 will become possible for systems beyond the Local Group such as Sculptor and M81 Group galaxies. Detailed studies of a greater number of individual stars in a greater number of spiral galaxies and their satellites will test hierarchical assembly in new ways because dynamical and chemical evolution models predict different outcomes for halos of different masses in different environments.

The final inspiral and coalescence of a black hole binary can produce highly beamed gravitational wave radiation. To conserve linear momentum, the black hole remnant can recoil with "kick" velocity as high as 4000 km/s. We present two sets of full N-body simulations of recoiling massive black holes (MBH) in high-resolution, non-axisymmetric potentials. The host to the first set of simulations is the main halo of the Via Lactea I simulation (Diemand et al. 2007). The nature of the resulting orbits is investigated through a numerical model where orbits are integrated assuming an evolving, triaxial NFW potential, and dynamical friction is calculated directly from the velocity dispersion along the major axes of the main halo of Via Lactea I. By comparing the triaxial case to a spherical model, we find that the wandering time spent by the MBH is significantly increased due to the asphericity of the halo. For kicks larger than 200 km/s, the remnant MBH does not return to the inner 200 pc within 1 Gyr, a timescale an order of magnitude larger than the upper limit of the estimated QSO lifetime. The second set of simulations is run using the outcome of a high-resolution gas-rich merger (Mayer et al. 2007) as host potential. In this case, a recoil velocity of 500 km/s cannot remove the MBH from the nuclear region.

We use the recently completed one billion particle Via Lactea II LambdaCDM simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo can not be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.

Our recently completed one billion particle Via Lactea II simulation of a Milky Way-sized dark matter halo resolves over 50,000 gravitationally bound clumps orbiting today within the virialized region of the main host. About 2,300 of these subhalos have one or more "progenitors" above 1e6 Msun at redshift 11, i.e. massive enough for their gas to have cooled via excitation of H2 and fragmented prior to the epoch of cosmic reionization. We count 4,500 such progenitors: if these were able to convert a fraction of their gas content into very metal-poor stars with a Salpeter initial mass function (IMF), they would be shining today with a visual magnitude M_V=6.7 per solar mass in stars. Assuming a universal baryon fraction, we show that mean star formation efficiencies as low as 0.1% in progenitors below 1e8 Msun would overproduce the abundance of the faint Galatic dwarf spheroidals observed by the Sloan Digital Sky Survey. Star formation at first light must have occurred either with an IMF lacking stars below 0.9 Msun, or was intrinsically very inefficient in small dark matter halos. If the latter, our results may be viewed as another hint of a minimum scale in galaxy formation.

It is a clear unique prediction of the cold dark matter paradigm of cosmological structure formation that galaxies form hierarchically and are embedded in massive, extended dark halos teeming with self-bound substructure or "subhalos". The amount and spatial distribution of subhalos around their host provide unique information and clues on the galaxy assembly process and the nature of the dark matter. Here we present results from the Via Lactea INCITE simulation, a one billion particle, one million cpu-hour simulation of the formation and evolution of a Galactic dark matter halo and its substructure population.

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

We present quantitative predictions for the detectability of individual Galactic dark matter subhalos in gamma-rays from dark matter pair annihilations in their centers. Our method is based on a hybrid approach, employing the highest resolution numerical simulations available (including the recently completed one billion particle Via Lactea II simulation) as well as analytical models for the extrapolation beyond the simulations' resolution limit. We include a self-consistent treatment of subhalo boost factors, motivated by our numerical results, and a realistic treatment of the expected backgrounds that individual subhalos must outshine. We show that for reasonable values of the dark matter particle physics parameters (M_X ~ 50 - 500 GeV and <sigma*v> ~ 10^-26 - 10^-25 cm^3/s) GLAST may very well discover a few, even up to several dozen, such subhalos, at 5 sigma significance, and some at more than 20 sigma. We predict that the majority of luminous sources would be resolved with GLAST's expected angular resolution. For most observer locations the angular distribution of detectable subhalos is consistent with a uniform distribution across the sky. The brightest subhalos tend to be massive (median Vmax of 24 km/s) and therefore likely hosts of dwarf galaxies, but many subhalos with Vmax as low as 5 km/s are also visible. Typically detectable subhalos are 20 - 40 kpc from the observer, and only a small fraction are closer than 10 kpc. The total number of observable subhalos has not yet converged in our simulations, and we estimate that we may be missing up to 3/4 of all detectable subhalos.

In cold dark matter cosmological models, structures form and grow by merging of smaller units. Numerical simulations have shown that such merging is incomplete; the inner cores of halos survive and orbit as "subhalos" within their hosts. Here we report a simulation that resolves such substructure even in the very inner regions of the Galactic halo. We find hundreds of very concentrated dark matter clumps surviving near the solar circle, as well as numerous cold streams. The simulation reveals the fractal nature of dark matter clustering: Isolated halos and subhalos contain the same relative amount of substructure and both have cuspy inner density profiles. The inner mass and phase-space densities of subhalos match those of recently discovered faint, dark matter-dominated dwarf satellite galaxies and the overall amount of substructure can explain the anomalous flux ratios seen in strong gravitational lenses. Subhalos boost gamma-ray production from dark matter annihilation, by factors of 4-15, relative to smooth galactic models. Local cosmic ray production is also enhanced, typically by a factor 1.4, but by more than a factor of ten in one percent of locations lying sufficiently close to a large subhalo. These estimates assume that gravitational effects of baryons on dark matter substructure are small.

We use the ``Via Lactea'' simulation to study the co-evolution of a Milky Way-size LambdaCDM halo and its subhalo population. While most of the host halo mass is accreted over the first 6 Gyr in a series of major mergers, the physical mass distribution [not M_vir(z)] remains practically constant since z=1. The same is true in a large sample of LambdaCDM galaxy halos. Subhalo mass loss peaks between the turnaround and virialization epochs of a given mass shell, and declines afterwards. 97% of the z=1 subhalos have a surviving bound remnant at the present epoch. The retained mass fraction is larger for initially lighter subhalos: satellites with maximum circular velocities Vmax=10 km/s at z=1 have today about 40% of their mass back then. At the first pericenter passage a larger average mass fraction is lost than during each following orbit. Tides remove mass in substructure from the outside in, leading to higher concentrations compared to field halos of the same mass. This effect, combined with the earlier formation epoch of the inner satellites, results in strongly increasing subhalo concentrations towards the Galactic center. We present individual evolutionary tracks and present-day properties of the likely hosts of the dwarf satellites around the Milky Way. The formation histories of ``field halos'' that lie today beyond the Via Lactea host are found to strongly depend on the density of their environment. This is caused by tidal mass loss that affects many field halos on eccentric orbits.

We show that most particle and subhalo orbits in simulated cosmological cold dark matter halos are surprisingly regular and periodic: The phase space structure of the outer halo regions shares some of the properties of the classical self-similar secondary infall model. Some of the outer branches are clearly visible in the radial velocity - radius plane at certain epochs. However, they are severely broadened in realistic, triaxial halos with non-radial, clumpy, mass accretion. This prevents the formation of high density caustics: Even in the best cases there are only broad, very small (<10 percent) enhancements in the spherical density profile. Larger fluctuations in rho(r) caused by massive satellites are common. Infall caustics are therefore too weak to affect lensing or dark matter annihilation experiments. Their detection is extremely challenging, as it requires a large number of accurate tracer positions and radial velocities in the outer halo. The stellar halo of the Milky Way is probably the only target where this could become feasible in the future.

The Via Lactea simulation of the dark matter halo of the Milky Way predicts the existence of many thousands of bound subhalos distributed approximately with equal mass per decade of mass. Here we show that: a) a similar steeply rising subhalo mass function is also present at redshift 0.5 in an elliptical-sized halo simulated with comparable resolution in a different cosmology. Compared to Via Lactea, this run produces nearly a factor of two more subhalos with large circular velocities; b) the fraction of Via Lactea mass brought in by subhalos that have a surviving bound remnant today with present-day peak circular velocity Vmax>2 km/s (>10 km/s) is 45% (30%); c) because of tidal mass loss, the number of subhalos surviving today that reached a peak circular velocity of >10 km/s throughout their lifetime exceeds half a thousand, five times larger than their present-day abundance and more than twenty times larger than the number of known satellites of the Milky Way; e) unless the circular velocity profiles of Galactic satellites peak at values significantly higher that expected from the stellar line-of-sight velocity dispersion, only about one in five subhalos with Vmax>20 km/s today must be housing a luminous dwarf; f) small dark matter clumps appear to be relatively inefficient at forming stars even well beyond the virial radius; g) the observed Milky Way satellites appear to follow the overall dark matter distribution of Via Lactea, while the largest simulated subhalos today are found preferentially at larger radii; h) subhalos have central densities that increase with Vmax and reach 0.1-0.3 Msun/pc3 comparable to the central densities inferred in dwarf spheroidals with core radii >250 pc.

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

During the epoch of reionization, the formation of radiation sources is accompanied by the growth of a H- photodissociating flux. We estimate the impact of this flux on the formation of molecular hydrogen and cooling in the first galaxies, assuming different types of radiation sources (e.g. Pop II and Pop III stars, miniquasars). We find that H- photodissociation reduces the formation of H2 molecules by a factor of ~1+1000k_s*x/(f_esc*delta), where x is the mean ionized fraction in the IGM, f_esc is the fraction of ionizing photons that escape from their progenitor halos, delta is the local gas overdensity and k_s is an order unity constant which depends on the type of radiation source. By the time a significant fraction of the universe becomes ionized, H- photodissociation may significantly reduce the H2 abundance and, with it, the primordial star formation rate, delaying the progress of reionization.

We present a study of the shapes, orientations, and alignments of Galactic dark matter subhalos in the ``Via Lactea'' simulation of a Milky Way-size LCDM host halo. Whereas isolated dark matter halos tend to be prolate, subhalos are predominantly triaxial. Overall subhalos are more spherical than the host halo, with minor to major and intermediate to major axis ratios of 0.68 and 0.83, respectively. Like isolated halos, subhalos tend to be less spherical in their central regions. The principal axis ratios are independent of subhalo mass, when the shapes are measured within a physical scale like r_Vmax, the radius of the peak of the circular velocity curve. Subhalos tend to be slightly more spherical closer to the host halo center. The spatial distribution of the subhalos traces the prolate shape of the host halo when they are selected by the largest V_max they ever had, i.e. before they experienced strong tidal mass loss. The subhalos' orientation is not random: the major axis tends to align with the direction towards the host halo center. This alignment disappears for halos beyond 3 r_200 and is more pronounced when the shapes are measured in the outer regions of the subhalos. The radial alignment is preserved during a subhalo's orbit and they become elongated during pericenter passage, indicating that the alignment is likely caused by the host halo's tidal forces. These tidal interactions with the host halo act to make subhalos rounder over time.

We discuss the possibility of GLAST detecting gamma-rays from the annihilation of neutralino dark matter in the Galactic halo. We have used "Via Lactea", currently the highest resolution simulation of Galactic cold dark matter substructure, to quantify the contribution of subhalos to the annihilation signal. We present a simulated allsky map of the expected gamma-ray counts from dark matter annihilation, assuming standard values of particle mass and cross section. In this case GLAST should be able to detect the Galactic center and several individual subhalos.

We present initial results from ``Via Lactea'', the highest resolution simulation to date of Galactic CDM substructure. It follows the formation of a Milky Way-size halo with Mvir=1.8x10^12 Msun in a WMAP 3-year cosmology, using 234 million particles. Over 10,000 subhalos can be identified at z=0: Their cumulative mass function is well-fit by N(>Msub)= 0.0064 (Msub/Mvir)^(-1) down to Msun=4x10^6 Msun. The total mass fraction in subhalos is 5.3%, while the fraction of surface mass density in substructure within a projected distance of 10 kpc from the halo center is 0.3%. Because of the significant contribution from the smallest resolved subhalos, these fractions have not converged yet. Sub-substructure is apparent in all the larger satellites, and a few dark matter lumps are resolved even in the solar vicinity. The number of dark satellites with peak circular velocities above 10 km/s (5 km/s) is 124 (812): of these, 5 (26) are found within 0.1 Rvir, a region that appeared practically smooth in previous simulations. The neutralino self-annihilation gamma-ray emission from dark matter clumps is approximately constant per subhalo mass decade. Therefore, while in our run the contribution of substructure to the gamma-ray luminosity of the Galactic halo amounts to only 40% of the total spherically-averaged smooth signal, we expect this fraction to grow significantly as resolution is increased further. An all-sky map of the expected annihilation gamma-ray flux reaching a fiducial observer at 8 kpc from the Galactic center shows that at the current resolution a small number of subhalos start to be bright enough to be visible against the background from the smooth density field surrounding the observer.

We describe a new concept for a MEMS-based active spatial filter for astronomical spectroscopy. The goal of this device is to allow the use of a diffraction-limited spectrometer on a seeing limited observation at improved throughput over a comparable seeing-limited spectrometer, thus reducing the size and cost of the spectrometer by a factor proportional to r0/D (For the case of a 10 meter telescope this size reduction will be approximately a factor of 25 to 50). We use a fiber-based integral field unit (IFU) that incorporates an active MEMS mirror array to feed an astronomical spectrograph. A fast camera is used in parallel to sense speckle images at a spatial resolution of lambda/D and at a temporal frequency greater than that of atmospheric fluctuations. The MEMS mirror-array is used as an active shutter to feed speckle images above a preset intensity threshold to the spectrometer, thereby increasing the signal-to-noise ratio (SNR) of the spectrogram. Preliminary calculations suggests an SNR improvement of a factor of about 1.4. Computer simulations have shown an SNR improvement of 1.1, but have not yet fully explored the parameter space.

Earth-mass ``microhalos'' may be the first objects to virialize in the early universe. Their ability to survive the hierarchical clustering process as substructure in the larger halos that form subsequently has implications for dark matter detection experiments. We present a large N-body simulation of early substructure in a supersymmetric cold dark matter (SUSY-CDM) scenario characterized by an exponential cutoff in the power spectrum at M_c=10^-6 Msun. The simulation resolves a 0.014 Msun parent ``SUSY'' halo at z=75 with 14 million particles. On these scales the effective index of the power spectrum approaches -3, and a range of mass scales collapse almost simultaneously. Compared to a z=0 galaxy cluster substructure within our SUSY host is less evident both in phase-space and in physical space, and it is less resistant against tidal disruption. As the universe expands by a factor of 1.3, we find that between 20 and 40 percent of well-resolved SUSY substructure is destroyed, compared to only about 1 percent in the low-redshift cluster. Nevertheless SUSY substructure is just as abundant as in z=0 galaxy clusters, i.e. the normalized mass and circular velocity functions are very similar. The dark matter self-annihilation gamma-ray luminosity from bound subhalos and other deviations from a smooth spherical configuration is at least comparable to the spherically-averaged signal in the SUSY host, and at least three times larger than the spherically-averaged signal in the cluster host. Such components must be taken into account when estimating the total cosmological extragalactic gamma-ray annihilation background.(ABRIDGED)

We use numerical hydrodynamical simulations of early structure formation in a LCDM universe to investigate the spin temperature and 21cm brightness of the diffuse intergalactic medium (IGM) prior to the epoch of cosmic reionization, at z<20. In the absence of any radiative heating, collisions between neutral hydrogen atoms can efficiently decouple the spin temperature from the CMB only in dense minihalos and filaments. Shock heated gas shines in emission, surrounded by cooler gas visible in absorption. In the case of a warm, mostly neutral IGM, produced here by X-ray emission from an early miniquasar, the 21cm signal is strongly enhanced. Even slightly overdense filaments now shine in emission against the CMB, possibly allowing future radio arrays to probe the distribution of neutral hydrogen before reionization.

The thermonuclear runaway that culminates in the explosion of a Chandrasekhar mass white dwarf as a Type Ia supernova begins centuries before the star actually explodes. Here, using a 3D anelastic code, we examine numerically the convective flow during the last minute of that runaway, a time that is crucial in determining just where and how often the supernova ignites. We find that the overall convective flow is dipolar, with the higher temperature fluctuations in an outbound flow preferentially on one side of the star. Taken at face value, this suggests an asymmetric ignition that may well persist in the geometry of the final explosion. However, we also find that even a moderate amount of rotation tends to fracture this dipole flow, making ignition over a broader region more likely. Though our calculations lack the resolution to study the flow at astrophysically relevant Rayleigh numbers, we also speculate that the observed dipolar flow will become less organized as the viscosity becomes very small. Motion within the dipole flow shows evidence of turbulence, suggesting that only geometrically large fluctuations (~1 km) will persist to ignite the runaway. We also examine the probability density function for the temperature fluctuations, finding evidence for a Gaussian, rather than exponential distribution, which suggests that ignition sparks may be strongly spatially clustered.

We investigate the environmental impact of the first active galactic nuclei that may have formed ~150 Myr after the big bang in low-mass ~10^6 Msun minihaloes. Using Enzo, an adaptive-mesh refinement cosmological hydrodynamics code, we carry out three-dimensional simulations of the radiative feedback from `miniquasars' powered by intermediate-mass black holes. We follow the non-equilibrium multispecies chemistry of primordial gas in the presence of a point source of X-ray radiation, which starts shining in a rare high-sigma peak at z=21 and emits a power-law spectrum in the 0.2-10 keV range. We find that, after one Salpeter time-scale, the miniquasar has heated up the simulation box to a volume-averaged temperature of 2800 K. The mean electron and H2 fractions are now 0.03 and 4e-5: the latter is 20 times larger than the primordial value, and will delay the buildup of a uniform UV photodissociating background. The net effect of the X-rays is to reduce gas clumping in the IGM by as much as a factor of 3. While the suppression of baryonic infall lowers the gas mass fraction at overdensities delta in the range 20-2000, enhanced molecular cooling increases the amount of dense material at delta>2000. In many haloes within the proximity of our miniquasar the H2-boosting effect of X-rays is too weak to overcome heating, and the cold and dense gas mass actually decreases. We find little evidence for an entropy floor in gas at intermediate densities preventing gas contraction and H2 formation. Overall, the radiative feedback from X-rays enhances gas cooling in lower-sigma peaks that are far away from the initial site of star formation, thus decreasing the clustering bias of the early pregalactic population, but does not appear to dramatically reverse or promote the collapse of pregalactic clouds as a whole. (abridged)

We investigate the effect of dark energy on the density profiles of dark matter haloes with a suite of cosmological N-body simulations and use our results to test analytic models. We consider constant equation of state models, and allow both w>-1 and w<-1. Using five simulations with w ranging from -1.5 to -0.5, and with more than ~1600 well-resolved haloes each, we show that the halo concentration model of Bullock et al. (2001) accurately predicts the median concentrations of haloes over the range of w, halo masses, and redshifts that we are capable of probing. We find that the Bullock et al. (2001) model works best when halo masses and concentrations are defined relative to an outer radius set by a cosmology-dependent virial overdensity. For a fixed power spectrum normalization and fixed-mass haloes, larger values of w lead to higher concentrations and higher halo central densities, both because collapse occurs earlier and because haloes have higher virial densities. While precise predictions of halo densities are quite sensitive to various uncertainties, we make broad comparisons to galaxy rotation curve data. At fixed power spectrum normalization (fixed sigma_8), w>-1 quintessence models seem to exacerbate the central density problem relative to the standard w=-1 model. Meanwhile w<-1 models help to reduce the apparent discrepancy. We confirm that the Jenkins et al. (2001) halo mass function provides an excellent approximation to the abundance of haloes in our simulations and extend its region of validity to include models with w<-1.

Gravitational lensing can amplify the apparent brightness of distant sources. Images that are highly magnified are often part of multiply-imaged systems, but we consider the possibility of having large magnifications without additional detectable images. In rare but non-negligible situations, lensing can produce a singly highly magnified image; this phenomenon is mainly associated with massive cluster-scale halos (>~1e13.5 Msun). Alternatively, lensing can produce multiply-imaged systems in which the extra images are either unresolved or too faint to be detectable. This phenomenon is dominated by galaxies and lower-mass halos (<~1e12 Msun), and is very sensitive to the inner density profile of the halos. Although we study the general problem, we customize our calculations to four quasars at redshift z~6 in the Sloan Digital Sky Survey (SDSS), for which Richards et al. (2004) have ruled out the presence of extra images down to an image splitting of 0.3" and a flux ratio of f=0.01. We predict that 9-29% of all z~6 quasars that are magnified by a factor of mu>10 would lack detectable extra images, with 5-10% being true singly-imaged systems. The maximum of 29% is reached only in the unlikely event that all low-mass (<~1e10 Msun) halos have highly concentrated (isothermal) profiles. In more realistic models where dwarf halos have flatter (NFW) inner profiles, the maximum probability is ~10%. We conclude that the probability that all four SDSS quasars are magnified by a factor of 10 is <~1e-4. The only escape from this conclusion is if there are many (>10) multiply-imaged z~6 quasars in the SDSS database that have not yet been identified, which seems unlikely. In other words, lensing cannot explain the brightnesses of the z~6 quasars, and models that invoke lensing to avoid having billion-Msun black holes in the young universe are not viable.

We study lens statistics in flat, low-density universes with different equations of state $w=p_Q/\rho_Q$ for the dark energy component. Dark energy modifies the distance-redshift relation and the mass function of dark matter halos leading to changes in the lensing optical depth as a function of image separation. Those effects must, however, be distinguished from effects associated with the structure of dark matter halos. Baryonic cooling causes galaxy-mass halos to have different central density profiles than group- and cluster-mass halos, which causes the distribution of normal arcsecond-scale lenses to differ from the distribution of ``wide-separation'' ($\Delta\theta \gtrsim 4\arcsec$) lenses. Fortunately, the various parameters related to cosmology and halo structure have very different effects on the overall image separation distribution: (1) the abundance of wide-separation lenses is exremely sensitive (by orders of magnitude) to the distribution of ``concentration'' parameters for massive halos modeled with the Navarro-Frenk-White profile; (2) the transition between normal and wide-separation lenses depends mainly on the mass scale where baryonic cooling ceases to be efficient; and (3) dark energy has effects at all image separation scales. While current lens samples cannot usefully constrain all of the parameters, ongoing and future imaging surveys should discover hundreds or thousands of lenses and make it possible to disentangle the various effects and constrain all of the parameters simultaneously. (abridged)

The observable properties of a Type Ia supernova are sensitive to how the nuclear runaway ignites in a Chandrasekhar mass white dwarf - at a single point at its center, off-center, or at multiple points and times. We present a simple analytic model for the runaway based upon a combination of stellar mixing-length theory and recent advances in understanding Rayleigh-Benard convection. The convective flow just prior to runaway is likely to have a strong dipolar component, though higher multipoles may contribute appreciably at the very high Rayleigh number (10$^{25}$) appropriate to the white dwarf core. A likely outcome is multi-point ignition with an exponentially increasing number of ignition points during the few tenths of a second that it takes the runaway to develop. The first sparks ignite approximately 150 - 200 km off center, followed by ignition at smaller radii. Rotation may be important to break the dipole asymmetry of the ignition and give a healthy explosion.

The development of primordial inhomogeneities into the non-linear regime and the formation of the first astrophysical objects within dark matter halos mark the transition from a simple, neutral, cooling universe -- described by just a few parameters -- to a messy ionized one -- the realm of radiative, hydrodynamic, and star formation processes. The recent measurement by the WMAP satellite of a large optical depth to electron scattering implies that this transition must have begun very early, and that the universe was reionized at redshift z_ion=17\pm 5. It is an early generation of extremely metal-poor massive stars and/or `seed' accreting black holes in subgalactic halos that may have generated the ultraviolet radiation and mechanical energy that reheated and reionized most of the hydrogen in the cosmos. The detailed thermal, ionization, and chemical enrichment history of the universe during the crucial formative stages around z=10-20 depends on the power-spectrum of density fluctuations on small scales, the stellar initial mass function and star formation efficiency, a complex network of poorly understood `feedback' mechanisms, and remains one of the crucial missing links in galaxy formation and evolution studies.

After briefly describing the anelastic approximation and Glatzmaier's code, we present results from our preliminary studies of core convection during the hydrogen burning phase of a 15 solar mass star, as well as our most recent results concerning convection in the oxygen shell of a 25 solar mass star.