We perform zoom-in cosmological simulations of a suite of dwarf galaxies, examining the impact of cosmic-rays generated by supernovae, including the effect of diffusion. We first look at the effect of varying the uncertain cosmic ray parameters by repeatedly simulating a single galaxy. Then we fix the comic ray model and simulate five dwarf systems with virial masses range from 8-30 $\times 10^{10}$ Msun. We find that including cosmic ray feedback (with diffusion) consistently leads to disk dominated systems with relatively flat rotation curves and constant star formation rates. In contrast, our purely thermal feedback case results in a hot stellar system and bursty star formation. The CR simulations very well match the observed baryonic Tully-Fisher relation, but have a lower gas fraction than in real systems. We also find that the dark matter cores of the CR feedback galaxies are cuspy, while the purely thermal feedback case results in a substantial core.

We explore the impact of cosmic rays (CRs) on cosmological adaptive-mesh refinement simulations of a forming 10^12 Msolar halo, focusing on the circumgalactic medium (CGM), and its resulting low-redshift structure and composition. In contrast to a run with star formation and energetic feedback but no CRs, the CR-inclusive runs feature a CGM substantially enriched with CRs and with metals to roughly 0.1 Zsolar, thanks to robust, persistent outflows from the disk. The CR-inclusive CGMs also feature more diffuse gas at lower temperatures, down to 10^4 K, than the non-CR run, with diffuse material often receiving a majority of its pressure support from the CR proton fluid. We compare to recent observations of the CGM of L ~ L* galaxies at low redshift, including UV absorption lines within background quasar spectra. The combination of metal-enriched, CR-driven winds and large swaths of CR pressure-supported, cooler diffuse gas leads to a CGM that provides a better match to data from COS-Halos (for HI, SiIV, CIII and OVI) than the non-CR run. We also compare our models to recent, preliminary observations of diffuse gamma-ray emission in local group halos. For our lower CR-diffusion runs with kappa_CR in the range 0.3 to 1 x 10^28 cm^2/s, the CR enriched CGM produces an inconsistently high level of gamma emission. But the model with a relatively high kappa_CR = 3 x 10^28 cm^2/s provided a gamma-ray luminosity consistent with the extra-galactic gamma-ray background observed by FERMI and roughly consistent with preliminary measures of the emission from M31's CGM.

Recent observations have constrained the orbit and structure of the Large Magellanic Cloud (LMC), implying a well-constrained pericentric passage about the Milky Way (MW) ~ 50 Myr ago. In this scenario, the LMC's gaseous disk has recently experienced maximal ram pressure stripping, suggesting the current extent of its HI disk directly probes the medium in which it is moving. From the observed stellar and HI distributions of the system we find evidence of a truncated gas profile along the windward ``leading edge' of the LMC disk, despite a far more extended stellar component. We explore the implications of this ram pressure stripping signature, using both analytic prescriptions and full 3-dimensional hydrodynamic simulations of the LMC. Our simulations subject the system to a headwind whose velocity components correspond directly to the recent orbital history of the LMC. We vary the density of this headwind, using a variety of sampled parameters for a Beta-profile for a theoretical MW circumgalactic medium (CGM), comparing the resulting HI morphology directly to observations of the LMC HI and stellar components. This model can match the radial extent of the LMC's leading (windward) edge only in scenarios where the MW CGM density at pericentric passage is n(R = 48.2 +/- 5 kpc) = 1.1 (+.44/-.45) x 1e-4 cm^-3. The implied pericentric density proves insensitive to both the broader CGM structure and temperature profile, thus providing a model-independent constraint on the local gas density. This result imposes an important constraint on the density profile of the MW's CGM, and thus the total baryon content of the MW. From our work, assuming a Beta-profile valid to ~ Rvir, we infer a total diffuse CGM mass M(300 kpc) = 2.6 +/- 1.4 x 1e10 Msun or approximately 15% of a 1e12 Msun MW's baryonic mass budget.

We investigate the dynamical impact of cosmic rays in cosmological simulations of galaxy formation using adaptive-mesh refinement simulations of a $10^{12}$ solar mass halo. In agreement with previous work, a run with only our standard thermal energy feedback model results in a massive spheroid and unrealistically peaked rotation curves. However, the addition of a simple two-fluid model for cosmic rays drastically changes the morphology of the forming disk. We include an isotropic diffusive term and a source term tied to star formation due to (unresolved) supernova-driven shocks. Over a wide range of diffusion coefficients, the CRs generate thin, extended disks with a significantly more realistic (although still not flat) rotation curve. We find that the diffusion of CRs is key to this process, as they escape dense star forming clumps and drive outflows within the more diffuse ISM.

Galactic-scale winds are a generic feature of massive galaxies with high star formation rates across a broad range of redshifts. Despite their importance, a detailed physical understanding of what drives these mass-loaded global flows has remained elusive. In this paper, we explore the dynamical impact of cosmic rays by performing the first three-dimensional, adaptive mesh refinement simulations of an isolated starbursting galaxy that includes a basic model for the production, dynamics and diffusion of galactic cosmic rays. We find that including cosmic rays naturally leads to robust, massive, bipolar outflows from our 10^12 Msun halo, with a mass-loading factor Mout/SFR = 0.3 for our fiducial run. Other reasonable parameter choices led to mass-loading factors above unity. The wind is multiphase and is accelerated to velocities well in excess of the escape velocity. We employ a two-fluid model for the thermal gas and relativistic CR plasma and model a range of physics relevant to galaxy formation, including radiative cooling, shocks, self-gravity, star formation, supernovae feedback into both the thermal and CR gas, and isotropic CR diffusion. Injecting cosmic rays into star-forming regions can provide significant pressure support for the interstellar medium, suppressing star formation and thickening the disk. We find that CR diffusion plays a central role in driving superwinds, rapidly transferring long-lived CRs from the highest density regions of the disk to the ISM at large, where their pressure gradient can smoothly accelerate the gas out of the disk.

Our local Hubble volume might be contained within a bubble that nucleated in a false vacuum with only two large spatial dimensions. We study bubble collisions in this scenario and find that they generate gravity waves, which are made possible in this context by the reduced symmetry of the global geometry. These gravity waves would produce B-mode polarization in the cosmic microwave background, which could in principle dominate over the inflationary background.

In the context of eternal inflation, cosmological predictions depend on the choice of measure to regulate the diverging spacetime volume. The spectrum of inflationary perturbations is no exception, as we demonstrate by comparing the predictions of the fat geodesic and causal patch measures. To highlight the effect of the measure---as opposed to any effects related to a possible landscape of vacua---we take the cosmological model, including the model of inflation, to be fixed. We also condition on the average CMB temperature accompanying the measurement. Both measures predict a 1-point expectation value for the gauge-invariant Newtonian potential, which takes the form of a (scale-dependent) monopole, in addition to a related contribution to the 3-point correlation function, with the detailed form of these quantities differing between the measures. However, for both measures both effects are well within cosmic variance. Our results make clear the theoretical relevance of the measure, and at the same time validate the standard inflationary predictions in the context of eternal inflation.

To make predictions for an eternally inflating "multiverse", one must adopt a procedure for regulating its divergent spacetime volume. Recently, a new test of such spacetime measures has emerged: normal observers - who evolve in pocket universes cooling from hot big bang conditions - must not be vastly outnumbered by "Boltzmann brains" - freak observers that pop in and out of existence as a result of rare quantum fluctuations. If the Boltzmann brains prevail, then a randomly chosen observer would be overwhelmingly likely to be surrounded by an empty world, where all but vacuum energy has redshifted away, rather than the rich structure that we observe. Using the scale-factor cutoff measure, we calculate the ratio of Boltzmann brains to normal observers. We find the ratio to be finite, and give an expression for it in terms of Boltzmann brain nucleation rates and vacuum decay rates. We discuss the conditions that these rates must obey for the ratio to be acceptable, and we discuss estimates of the rates under a variety of assumptions.

Our universe may have formed via bubble nucleation in an eternally-inflating background. Furthermore, the background may have a compact dimension--the modulus of which tunnels out of a metastable minimum during bubble nucleation--which subsequently grows to become one of our three large spatial dimensions. Then the reduced symmetry of the background is equivalent to anisotropic initial conditions in our bubble universe. We compute the inflationary spectrum in such a scenario and, as a first step toward understanding the effects of anisotropy, project it onto spherical harmonics. The resulting spectrum exhibits anomalous multipole correlations, their relative amplitude set by the present curvature parameter, which extend to arbitrarily large multipole moments. This raises the possibility of future detection, if slow-roll inflation does not last too long within our bubble. A full understanding of the observational signal must account for the effects of background anisotropy on photon free streaming, and is left to future work.

Our universe may be contained in one among a diverging number of bubbles that nucleate within an eternally inflating multiverse. A promising measure to regulate the diverging spacetime volume of such a multiverse is the scale-factor cutoff, one feature of which is bubbles are not rewarded for having a longer duration of slow-roll inflation. Thus, depending on the landscape distribution of the number of e-folds of inflation among bubbles like ours, we might hope to measure spacetime curvature. We study a recently proposed cartoon model of inflation in the landscape and find a reasonable chance (about ten percent) that the curvature in our universe is well above the value expected from cosmic variance. Anthropic selection does not strongly select for curvature as small as is observed (relative somewhat larger values), meaning the observational bound on curvature can be used to rule out landscape models that typically give too little inflation.

Recently a mechanism was proposed whereby the primordial density perturbations are generated at the end of inflation. We continue the analysis of the proposed model of this mechanism and calculate the maximum extent to which the density perturbations produced via this model can dominate over those of the standard inflationary paradigm. In addition, we provide a straightforward variation of this model which allows for greater amplification of the density perturbations. Finally, we show that a variation in the implementation of the original model results in significant non-gaussianities in the resulting spectrum of density perturbations. The level of non-gaussianities can be made to saturate the current observational bound.

Arguably a major success of the landscape picture is the prediction of a small, non-zero vacuum energy density. The details of this prediction depends in part on how the diverging spacetime volume of the multiverse is regulated, a question that remains unresolved. One proposal, the causal diamond measure, has demonstrated many phenomenological successes, including predicting a distribution of positive vacuum energy densities in good agreement with observation. In the string landscape, however, the vacuum energy density is expected to take positive and negative values. We find the causal diamond measure gives a poor fit to observation in such a landscape -- in particular, 99.6% of observers in galaxies seemingly just like ours measure a vacuum energy density smaller than we do, most of them measuring it to be negative.

It is well known that anthropic selection from a landscape with a flat prior distribution of cosmological constant Lambda gives a reasonable fit to observation. However, a realistic model of the multiverse has a physical volume that diverges with time, and the predicted distribution of Lambda depends on how the spacetime volume is regulated. We study a simple model of the multiverse with probabilities regulated by a scale-factor cutoff, and calculate the resulting distribution, considering both positive and negative values of Lambda. The results are in good agreement with observation. In particular, the scale-factor cutoff strongly suppresses the probability for values of Lambda that are more than about ten times the observed value. We also discuss several qualitative features of the scale-factor cutoff, including aspects of the distributions of the curvature parameter Omega and the primordial density contrast Q.

It is possible that the scale of gravity, parameterized by the apparent Planck mass, may obtain different values within different universes in an encompassing multiverse. We investigate the range over which the Planck mass may scan while still satisfying anthropic constraints. The window for anthropically allowed values of the Planck mass may have important consequences for landscape predictions. For example, if the likelihood to observe some value of the Planck mass is weighted by the inflationary expansion factors of the universes that contain that value, then it appears extremely unlikely to observe the value of the Planck mass that is measured within our universe. This is another example of the runaway inflation problem discussed in recent literature. We also show that the window for the Planck mass significantly weakens the anthropic constraint upon the cosmological constant when both are allowed to vary over a landscape.

Fluctuations in the mass and decay rate of a heavy particle which for some period dominates the energy density of the universe are known to lead to adiabatic density perturbations. We show that generically the annihilation cross section of the same particle also receives fluctuations, which leads to entropy perturbations at freeze-out. If the particle comes to dominate the energy density of the universe and subsequently decays, this leads to an additional source of adiabatic density perturbations. On the other hand, non-adiabatic density perturbations result when the particle does not decay but contributes to the observed dark matter.