Horizon Run 5 (HR5) is a cosmological hydrodynamics simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1 kpc. This enormous dynamic range allows us to simultaneously capture the physics of the cosmic web on very large scales and account for the formation and evolution of dwarf galaxies on much smaller scales. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of 1049x114x114 Mpc^3. The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, reionization, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and active galactic nuclei (AGN) feedback in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OpenMP version of the RAMSES code, specifically targeted for modern many-core many thread parallel architectures. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder to analyse the large outputs of HR5. The simulation successfully reproduces many observations, such as the cosmic star formation history, connectivity of galaxy distribution and stellar mass functions. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one should carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.

A new detector, the Fermilab Holometer, consists of separate yet identical 39-meter Michelson interferometers. Strain sensitivity achieved is better than $10^{-21} /{\sqrt{\rm{Hz}}}$ between 1 to 13 MHz from a 130-hr dataset. This measurement exceeds the sensitivity and frequency range made from previous high frequency gravitational wave experiments by many orders of magnitude. Constraints are placed on a stochastic background at 382 Hz resolution. The 3$\sigma$ upper limit on $\Omega_{\rm{GW}}$, the gravitational wave energy density normalized to the closure density, ranges from $5.6 \times 10^{12}$ at 1 MHz to $8.4 \times 10^{15}$ at 13 MHz. Another result from the same dataset is a search for nearby primordial black hole binaries (PBHB). There are no detectable monochromatic PBHBs in the mass range $0.83$ - $3.5 \times 10^{21}$g between the earth and the moon. Projections for a chirp search with the same dataset increases the mass range to $0.59 - 2.5 \times 10^{25}$g and distances out to Jupiter. This result presents a new method for placing limits on a poorly constrained mass range of primordial black holes. Additionally, solar system searches for PBHBs place limits on their contribution to the total dark matter fraction.

Suppression of the scalar power spectrum on large scales is one way to reconcile the tension between Planck and BICEP2 data. This suppression can occur by introducing a phase transition from the fast-roll phase to the slow-roll phase in a single field inflation model. In this paper we consider a deformed single field inflation model in terms of three SO(3) symmetric moduli fields. We find that spatially linear solutions for the moduli fields induces a phase transition during the early stage of the inflation and the suppression of scalar power spectrum at large scale perturbation modes.

We investigate the tensor and the scalar perturbations in the symmetric bouncing universe driven by one ordinary field and its Lee-Wick partner field which is a ghost. We obtain the even- and the odd-mode functions of the tensor perturbation in the matter-dominated regime. The tensor perturbation grows in time during the contracting phase of the Universe, and decays during the expanding phase. The power spectrum for the tensor perturbation is evaluated and the spectral index is given by $n_{\rm T} =6$. We add the analysis on the scalar perturbation by inspecting the even- and the odd-mode functions in the matter-dominated regime, which was studied numerically in our previous work. We conclude that the comoving curvature by the scalar perturbation is constant in the super-horizon scale and starts to decay in the far sub-horizon scale while the Universe expands.

We investigate the scalar perturbation in the Lee-Wick bouncing universe driven by an ordinary scalar field plus a ghost field. We consider only a symmetric evolution of the universe and the scalar fields about the bouncing point. The gauge invariant Sasaki-Mukhanov variable is numerically solved in the spatially flat gauge. We find a new form of the initial perturbation growing during the contracting phase. After the bouncing, this growing mode stabilizes to a constant mode which is responsible for the late-time power spectrum.

We investigate the decay of an inhomogeneous D1-brane wrapped on a $S^1$ with an electric field. The model that we consider consists of an array of tachyon kink and anti-kink with a constant electric flux. Beginning with an initially static configuration, we numerically evolve the tachyon field with some perturbations under a fixed boundary condition at diametrically opposite points on the circle $S^1$. When the electric flux is smaller than the critical value, the tachyon kink becomes unstable; the tachyon field rolls down the potential, and the lower dimensional D0- and $\bar {\rm D}0$-brane become thin, which resembles the caustic formation known for this type of the system in the literature. For the supercritical values of the electric flux, the tachyon kink remains stable.

Vacuum structure and global cosmic strings are analyzed in the effective theory of self-interacting O(2) scalar fields on (3+1)-manifolds with conical singularities. In the context of one-loop effective action computed by heat-kernel methods with $\zeta$-function regularization, we find an inhomogeneous vacuum of minimum energy and suggest some reason why low-energy global strings are likely to be generated at the conical singularities.