In this work we investigate the systematic uncertainties that arise from the calculation of the peculiar velocity when estimating the Hubble constant ($H_0$) from gravitational wave standard sirens. We study the GW170817 event and the estimation of the peculiar velocity of its host galaxy, NGC 4993, when using Gaussian smoothing over nearby galaxies. NGC 4993 being a relatively nearby galaxy, at $\sim 40 \ {\rm Mpc}$ away, is subject to a significant effect of peculiar velocities. We demonstrate a direct dependence of the estimated peculiar velocity value on the choice of smoothing scale. We show that when not accounting for this systematic, a bias of $\sim 200 \ {\rm km \ s ^{-1}}$ in the peculiar velocity incurs a bias of $\sim 4 \ {\rm km \ s ^{-1} \ Mpc^{-1}}$ on the Hubble constant. We formulate a Bayesian model that accounts for the dependence of the peculiar velocity on the smoothing scale and by marginalising over this parameter we remove the need for a choice of smoothing scale. The proposed model yields $H_0 = 68.6 ^{+14.0}_{-8.5}~{\rm km\ s^{-1}\ Mpc^{-1}}$. We demonstrate that under this model a more robust unbiased estimate of the Hubble constant from nearby GW sources is obtained.

Recent works combining neutrino oscillation and cosmological data to determine the neutrino hierarchy found a range of odds in favour of the normal hierarchy. These results arise from differing approaches to incorporating prior knowledge about neutrinos. We develop a hierarchy-agnostic prior and show that the hierarchy cannot be conclusively determined with current data. The determination of the hierarchy is limited by the neutrino mass scale $\Sigma_{\nu}$ measurement. We obtain a target precision of $\sigma(\Sigma_{\nu}) = 0.014$ eV, necessary for conclusively establishing the normal hierarchy with future data.

We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative phase is approximately governed by a relativistic wave equation exhibiting true and false vacuum configurations. In previous work, a linear stability analysis identified exponentially growing short-wavelength modes driven by the time-dependent coupling. These modes threaten to destabilize the analog false vacuum. Here, we employ numerical simulations of the coupled Gross-Pitaevski equations (GPEs) to determine the non-linear evolution of these linearly unstable modes. We find that unless a physical mechanism modifies the GPE on short length scales, the analog false vacuum is indeed destabilized. We briefly discuss various physically expected corrections to the GPEs that may act to remove the exponentially unstable modes. To investigate the resulting dynamics in cases where such a removal mechanism exists, we implement a hard UV cutoff that excludes the unstable modes as a simple model for these corrections. We use this to study the range of phenomena arising from such a system. In particular, we show that by modulating the strength of the time-dependent coupling, it is possible to observe the crossover between a second and first order phase transition out of the false vacuum.

We introduce a new picture of vacuum decay which, in contrast to existing semiclassical techniques, provides a real-time description and does not rely on classically-forbidden tunneling paths. Using lattice simulations, we observe vacuum decay via bubble formation by generating realizations of vacuum fluctuations and evolving with the classical equations of motion. The decay rate obtained from an ensemble of simulations is in excellent agreement with existing techniques. Future applications include bubble correlation functions, fast decay rates, and decay of non-vacuum states.

It has been recently suggested that oscillons produced in the early universe from certain asymmetric potentials continue to emit gravitational waves for a number of $e$-folds of expansion after their formation, leading to potentially detectable gravitational wave signals. We revisit this claim by conducting a convergence study using graphics processing unit (GPU)-accelerated lattice simulations and show that numerical errors accumulated with time are significant in low-resolution scenarios, or in scenarios where the run-time causes the resolution to drop below the relevant scales in the problem. Our study determines that the dominant, growing high frequency peak of the gravitational wave signals in the fiducial "hill-top model" in [arXiv:1607.01314] is a numerical artifact. This finding prompts the need for a more careful analysis of the numerical validity of other similar results related to gravitational waves from oscillon dynamics.

Analog condensed matter systems present an exciting opportunity to simulate early Universe models in table-top experiments. We consider a recent proposal for an analog condensed matter experiment to simulate the relativistic quantum decay of the false vacuum. In the proposed experiment, two ultra-cold condensates are coupled via a time-varying radio-frequency field. The relative phase of the two condensates in this system is approximately described by a relativistic scalar field with a potential possessing a series of false and true vacuum local minima. If the system is set up in a false vacuum, it would then decay to a true vacuum via quantum mechanical tunnelling. Should such an experiment be realized, it would be possible to answer a number of open questions regarding non-perturbative phenomena in quantum field theory and early Universe cosmology. In this paper, we illustrate a possible obstruction: the time-varying coupling that is invoked to create a false vacuum for the long-wavelength modes of the condensate leads to a destabilization of shorter wavelength modes within the system via parametric resonance. We focus on an idealized setup in which the two condensates have identical properties and identical background densities. Describing the system by the coupled Gross-Pitaevskii equations (GPE), we use the machinery of Floquet theory to perform a linear stability analysis, calculating the wavenumber associated with the first instability band for a variety of experimental parameters. However, we demonstrate that, by tuning the frequency of the time-varying coupling, it may be possible to push the first instability band outside the validity of the GPE, where dissipative effects are expected to damp any instabilities. This provides a viable range of experimental parameters to perform analog experiments of false vacuum decay.

Cosmic inflation, a period of accelerated expansion in the early universe, can give rise to large amplitude ultra-large scale inhomogeneities on distance scales comparable to or larger than the observable universe. The cosmic microwave background (CMB) anisotropy on the largest angular scales is sensitive to such inhomogeneities and can be used to constrain the presence of ultra-large scale structure (ULSS). We numerically evolve nonlinear inhomogeneities present at the beginning of inflation in full General Relativity to assess the CMB quadrupole constraint on the amplitude of the initial fluctuations and the size of the observable universe relative to a length scale characterizing the ULSS. To obtain a statistically significant number of simulations, we adopt a toy model in which inhomogeneities are injected along a preferred direction. We compute the likelihood function for the CMB quadrupole including both ULSS and the standard quantum fluctuations produced during inflation. We compute the posterior given the observed CMB quadrupole, finding that when including gravitational nonlinearities, ULSS curvature perturbations of order unity are allowed by the data, even on length scales not too much larger than the size of the observable universe. Our results illustrate the utility and importance of numerical relativity for constraining early universe cosmology.

We study collisions between nearly planar domain walls including the effects of small initial nonplanar fluctuations. These perturbations represent the small fluctuations that must exist in a quantum treatment of the problem. In a previous paper, we demonstrated that at the linear level a subset of these fluctuations experience parametric amplification as a result of their coupling to the planar symmetric background. Here we study the full three-dimensional nonlinear dynamics using lattice simulations, including both the early time regime when the fluctuations are well described by linear perturbation theory as well as the subsequent stage of fully nonlinear evolution. We find that the nonplanar fluctuations have a dramatic effect on the overall evolution of the system. Specifically, once these fluctuations begin to interact nonlinearly the split into a planar symmetric part of the field and the nonplanar fluctuations loses its utility. At this point the colliding domain walls dissolve, with the endpoint of this being the creation of a population of oscillons in the collision region. The original (nearly) planar symmetry has been completely destroyed at this point and an accurate study of the system requires the full three-dimensional simulation.

We study collisions between pairs of bubbles nucleated in an ambient false vacuum. For the first time, we include the effects of small initial (quantum) fluctuations around the instanton profiles describing the most likely initial bubble profile. Past studies of this problem neglect these fluctuations and work under the assumption that the collisions posess an exact SO(2,1) symmetry. We use three-dimensional lattice simulations to demonstrate that for double-well potentials, small initial perturbations to this symmetry can be amplified as the system evolves. Initially the amplification is well-described by linear perturbation theory around the SO(2,1) background, but the onset of strong nonlinearities amongst the fluctuations quickly leads to a drastic breaking of the original SO(2,1) symmetry and the production of oscillons in the collision region. We explore several single-field models, and we find it is hard to both realize inflation inside of a bubble and produce oscillons in a collision. Finally, we extend our results to a simple two-field model. The additional freedom allowed by the second field allows us to construct viable inflationary models that allow oscillon production in collisions. The breaking of the SO(2,1) symmetry allows for a new class of observational signatures from bubble collisions that do not posess azimuthal symmetry, including the production of gravitational waves which cannot be supported by an SO(2,1) spacetime.

This is the first paper in a series where we study collisions of nucleated bubbles taking into account the effects of small initial (quantum) fluctuations in a fully 3+1-dimensional setting. In this paper, we consider the evolution of linear fluctuations around highly symmetric though inhomogeneous backgrounds. We demonstrate that a large degree of asymmetry develops over time from tiny fluctuations superposed upon planar and SO(2,1) symmetric backgrounds. These fluctuations arise from zero-point vacuum oscillations, so excluding them by enforcing a spatial symmetry is inconsistent in a quantum treatment. We consider the limit of two colliding planar walls, with fluctuation mode functions characterized by the wavenumber transverse to the collision direction and a longitudinal shape along the collision direction $x$, which we solve for. Initially, the fluctuations obey a linear wave equation with a time- and space-dependent mass $m_{eff}(x,t)$. When the walls collide multiple times, $m_{eff}$ oscillates in time. We use Floquet theory to study the fluctuations and generalize techniques familiar from preheating to the case with many coupled degrees of freedom. This inhomogeneous case has bands of unstable transverse wavenumbers $k_\perp$ with exponentially growing mode functions. From the detailed spatial structure of the mode functions in $x$, we identify both broad and narrow parametric resonance generalizations of the homogeneous $m_{eff}(t)$ case of preheating. The unstable $k_\perp$ modes are longitudinally localized, yet can be described as quasiparticles in the Bogoliubov sense. We define an effective occupation number to show they are created in bursts for the case of well-defined collisions in the background. The transverse-longitudinal coupling accompanying nonlinearity radically breaks this localized particle description, with nonseparable 3D modes arising.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

We study in detail (p)reheating after multi-field inflation models with a particular focus on N-flation. We consider a variety of different couplings between the inflatons and the matter sector, including both quartic and trilinear interactions with a light scalar field. We show that the presence of multiple oscillating inflatons makes parametric resonance inefficient in the case of the quartic interactions. Moreover, perturbative processes do not permit a complete decay of the inflaton for this coupling. In order to recover the hot big bang, we must instead consider trilinear couplings. In this case we show that strong nonperturbative preheating is possible via multi-field tachyonic resonance. In addition, late-time perturbative effects do permit a complete decay of the condensate. We also study the production of gauge fields for several prototype couplings, finding similar results to the trilinear scalar coupling. During the course of our analysis we develop the mathematical theory of the quasi-periodic Mathieu equation, the multi-field generalization of the Floquet theory familiar from preheating after single field inflation. We also elaborate on the theory of perturbative decays of a classical inflaton condensate, which is applicable in single-field models also.