We study scalar bubble collisions in first-order phase transitions focusing on the relativistic limit. We propose 'trapping equation' which describes the wall behavior after collision, and test it with numerical simulations in several setups. We also examine the energy dynamics after collision and discuss its implications to gravitational wave production.

We study the gravitational wave (GW) production in extremely strong first order phase transitions where the latent heat density dominates the plasma energy density, $\alpha \gtrsim 1$. In such transitions, bubbles develop extremely thin and relativistic fluid configurations, resulting in strong shock waves after collisions. We first propose a strategy to understand the GW production in such a system by separating the problem into the propagation part and the collision part. Focusing on the former, we next develop an effective theory for the propagation of the relativistic fluid shells. Using this effective theory, we finally calculate the expected duration of the relativistic fluid configurations and discuss its implications to the GW production.

We study gravitational-wave production from bubble dynamics (bubble collisions and sound waves) during a cosmic first-order phase transition with an analytic approach. We first propose modeling the system with the thin-wall approximation but without the envelope approximation often adopted in the literature, in order to take bubble propagation after collisions into account. The bubble walls in our setup are considered as modeling the scalar field configuration and/or the bulk motion of the fluid. We next write down analytic expressions for the gravitational-wave spectrum, and evaluate them with numerical methods. It is found that, in the long-lasting limit of the collided bubble walls, the spectrum grows from $\propto f^3$ to $\propto f^1$ in low frequencies, showing a significant enhancement compared to the one with the envelope approximation. It is also found that the spectrum saturates in the same limit, indicating a decrease in the correlation of the energy-momentum tensor at late times. We also discuss the implications of our results to gravitational-wave production both from bubble collisions (scalar dynamics) and sound waves (fluid dynamics).

We study I-balls/oscillons, which are long-lived, quasi-periodic, and spatially localized solutions in real scalar field theories. Contrary to the case of Q-balls, there is no evident conserved charge that stabilizes the localized configuration. Nevertheless, in many classical numerical simulations, it has been shown that they are extremely long-lived. In this paper, we clarify the reason for the longevity, and show how the exponential separation of time scales emerges dynamically. Those solutions are time-periodic with a typical frequency of a mass scale of a scalar field. This observation implies that they can be understood by the effective theory after integrating out relativistic modes. We find that the resulting effective theory has an approximate global U(1) symmetry reflecting an approximate number conservation in the non-relativistic regime. As a result, the profile of those solutions is obtained via the bounce method, just like Q-balls, as long as the breaking of the U(1) symmetry is small enough. We then discuss the decay processes of the I-ball/oscillon by the breaking of the U(1) symmetry, namely the production of relativistic modes via number violating processes. We show that the imaginary part is exponentially suppressed, which explains the extraordinary longevity of I-ball/oscillon. In addition, we find that there are some attractor behaviors during the evolution of I-ball/oscillon that further enhance the lifetime. The validity of our effective theory is confirmed by classical numerical simulations. Our formalism may also be useful to study condensates of ultra light bosonic dark matter, such as fuzzy dark matter, and axion stars, for instance.

We investigate the sensitivity of future space-based interferometers such as LISA and DECIGO to the parameters of new particle physics models which drive a first-order phase transition in the early Universe. We first perform a Fisher matrix analysis on the quantities characterizing the gravitational wave spectrum resulting from the phase transition, such as the peak frequency and amplitude. We next perform a Fisher analysis for the quantities which determine the properties of the phase transition, such as the latent heat and the time dependence of the bubble nucleation rate. Since these quantities are determined by the model parameters of the new physics, we can estimate the expected sensitivities to such parameters. We illustrate this point by taking three new physics models for example: (1) models with additional isospin singlet scalars (2) a model with an extra real Higgs singlet, and (3) a classically conformal $B-L$ model. We find that future gravitational wave observations play complementary roles to future collider experiments in pinning down the parameters of new physics models driving a first-order phase transition.

Coherent oscillation of axions or axion-like particles may give rise to long-lived clumps, called axion stars, because of the attractive gravitational force or its self-interaction. Such a kind of configuration has been extensively studied in the context of oscillons without the effect of gravity, and its stability can be understood by an approximate conservation of particle number in a non-relativistic effective field theory (EFT). We extend this analysis to the case with gravity to discuss the longevity of axion stars and clarify the EFT expansion scheme in terms of gradient energy and Newton's constant. Our EFT is useful to calculate the axion star configuration and its classical lifetime without any ad hoc assumption. In addition, we derive a simple stability condition against small perturbations. Finally, we discuss the consistency of other non-relativistic effective field theories proposed in the literature.

We study gravitational-wave production from bubble collisions in a cosmic first-order phase transition, focusing on the possibility of model separation by the bubble nucleation rate dependence of the resulting gravitational-wave spectrum. By using the method of relating the spectrum with the two-point correlator of the energy-momentum tensor $\left< T(x)T(y) \right>$, we first write down analytic expressions for the spectrum with a Gaussian correction to the commonly used nucleation rate, $\Gamma \propto e^{\beta t}\rightarrow e^{\beta t-\gamma^2t^2}$, under the thin-wall and envelope approximations. Then we quantitatively investigate how the spectrum changes with the size of the Gaussian correction. It is found that the spectral shape shows ${\mathcal O}(10)\%$ deviation from $\Gamma \propto e^{\beta t}$ case for some physically motivated scenarios. We also briefly discuss detector sensitivities required to distinguish different spectral shapes.

Cosmic-ray anti-deuterium and anti-helium have long been suggested as probes of dark matter, as their secondary astrophysical production was thought extremely scarce. But how does one actually predict the secondary flux? Anti-nuclei are dominantly produced in pp collisions, where laboratory cross section data is lacking. We make a new attempt at tackling this problem by appealing to a scaling law of nuclear coalescence with the physical volume of the hadronic emission region. The same volume is probed by Hanbury Brown-Twiss (HBT) two-particle correlations. We demonstrate the consistency of the scaling law with systems ranging from central and off-axis AA collisions to pA collisions, spanning 3 orders of magnitude in coalescence yield. Extending the volume scaling to the pp system, HBT data allows us to make a new estimate of coalescence, that we test against preliminary ALICE pp data. For anti-helium the resulting cross section is 1-2 orders of magnitude higher than earlier estimates. The astrophysical secondary flux of anti-helium could be within reach of a five-year exposure of AMS02.

Secondary astrophysical production of $e^+$ and $\bar{p}$ cosmic rays is considered. Inclusive $\pi$, $K$, and $\bar{p}$ production cross sections in $pp$ collisions at large $\sqrt{s}$ are parametrised using recent experimental data at LHC energies. The astrophysical production rate ratio $Q_{e^+}/Q_{\bar{p}}$ is calculated for an input cosmic ray proton flux consistent with local measurements. At $10<E<100$$\sim$GeV the cosmic ray flux ratio $J_{e^+}/J_{\bar{p}}$ measured by AMS02 falls below the production rate ratio by about 50\%, while at high energy $E>100$$\sim$GeV the measured flux ratio coincides with the production rate ratio of the secondary source.

We consider gravitational wave production by bubble collisions during a cosmological first-order phase transition. In the literature, such spectra have been estimated by simulating the bubble dynamics, under so-called thin-wall and envelope approximations in a flat background metric. However, we show that, within these assumptions, the gravitational wave spectrum can be estimated in an analytic way. Our estimation is based on the observation that the two-point correlator of the energy-momentum tensor $\langle T(x)T(y)\rangle$ can be expressed analytically under these assumptions. Though the final expressions for the spectrum contain a few integrations that cannot be calculated explicitly, we can easily estimate it numerically. As a result, it is found that the most of the contributions to the spectrum come from single-bubble contribution to the correlator, and in addition the fall-off of the spectrum at high frequencies is found to be proportional to $f^{-1}$. We also provide fitting formulae for the spectrum.

We study the cosmological history of the classical conformal $B-L$ gauge extension of the standard model, in which the physical scales are generated via the Coleman-Weinberg-type symmetry breaking. Especially, we consider the thermal phase transition of the U$(1)_{B-L}$ symmetry in the early universe and resulting gravitational-wave production. Due to the classical conformal invariance, the phase transition tends to be a first-order one with ultra-supercooling, which enhances the strength of the produced gravitational waves. We show that, requiring (1) U$(1)_{B-L}$ is broken after the reheating, (2) the $B-L$ gauge coupling does not blow up below the Planck scale, (3) the thermal phase transition completes in almost all the patches in the universe, the gravitational wave spectrum can be as large as $\Omega_{\rm GW} \sim 10^{-8}$ at the frequency $f \sim 0.01$-$1$Hz for some model parameters, and a vast parameter region can be tested by future interferometer experiments such as eLISA, LISA, BBO and DECIGO.

In a wide class of new physics models, there exist scalar fields that obtain vacuum expectation values of high energy scales. We study the possibility that the standard model Higgs field has experienced first order phase transition at the high energy scale due to the couplings with these scalar fields. We estimate the amount of gravitational waves produced by the phase transition, and discuss observational consequences.

Cosmological dynamics of scalar field with a monomial potential $\phi^{n}$ with a general background equation of state is revisited. It is known that if $n$ is smaller than a critical value, the scalar field exhibits a coherent oscillation and if $n$ is larger it obeys a scaling solution without oscillation. We study in detail the case where $n$ is equal to the critical value, and find a peculiar scalar dynamics which is neither oscillating nor scaling solution, and we call it a pseudo scaling solution. We also discuss cosmological implications of a pseudo scaling scalar dynamics, such as the curvature perturbation and the domain wall problem.

If a real scalar field is dominated by non-relativistic modes, then it approximately conserves its particle number and obeys an equation that governs a complex scalar field theory with a conserved global U(1) symmetry. From this fact, it is shown that the I-ball (oscillon) can be naturally understood as a projection (e.g., real part) of the non-relativistic Q-ball solution. In particular, we clarify that the stability of the I-ball is guaranteed by the U(1) symmetry in the corresponding complex scalar field theory as long as the non-relativistic condition holds. We also discuss the longevity of I-ball from the perspective of the complex scalar field in terms of U(1) charge violating processes.

We revisit the cosmological aspects of axion models. In the high-scale inflation scenario, the Peccei-Quinn (PQ) symmetry is likely to be restored during/after inflation. If the curvature of the PQ scalar potential at the origin is smaller than its vacuum expectation value; for instance in a class of SUSY axion models, thermal inflation happens before the radial component of the PQ scalar (saxion) relaxes into the global minimum of the potential and the decay of saxion coherent oscillation would produce too much axion dark radiation. In this paper, we study how to avoid the overproduction of axion dark radiation with some concrete examples. We show that, by taking account of the finite-temperature dissipation effect appropriately, the overproduction constraint can be relaxed since the PQ scalar can take part in the thermal plasma again even after the PQ phase transition. We also show that it can be further relaxed owing to the late time decay of another heavy CP-odd scalar, if it is present.

We show that the local type non-Gaussianity in a class of curvaton models is suppressed, i.e. the non-linearity parameters f_{NL} and those related with higher order statistics can be at most O(1), even if the curvaton energy density is subdominant at the decay. This situation is naturally realized in a very simple curvaton potential with quadratic term plus quartic term.

We revisit the dynamics of the curvaton in detail taking account of effects from thermal environment, effective potential and decay/dissipation rate for general field values and couplings. We also consider the curvature perturbation generated through combinations of various effects: large scale modulation of the oscillation epoch, the effective dissipation rate and the timing at which the equation of state changes. In particular, we find that it tends to be difficult to explain the observed curvature perturbation by the curvaton mechanism without producing too large non-Gaussianity if the curvaton energy density is dissipated through thermal effects. Therefore, the interaction between the curvaton and light elements in the thermal plasma should be suppressed in order for the curvaton to survive the thermal dissipation.

The evolution of a coherently oscillating scalar field with $Z_2$ symmetry is studied in detail. We calculate the dissipation rate of the scalar field based on the closed time path formalism. Consequently, it is shown that the energy density of the coherent oscillation can be efficiently dissipated if the coupling constant is larger than the critical value, even though the scalar particle is stable due to the $Z_2$ symmetry.

We study in detail the dynamics of a scalar field in thermal bath with symmetry breaking potential. In particular, we focus on the process of trapping of a scalar field at an enhanced symmetry point through the thermal/non-thermal particle production, taking into account the interactions of produced particles with the standard model particles. As an explicit example, we revisit the saxion dynamics with an initial amplitude much larger than the Peccei-Quinn scale and show that the saxion trapping phenomenon happens for the most cases and it often leads to thermal inflation. We also study the saxion dynamics after thermal inflation, and it is shown that thermal dissipation effect on the saxion can relax the axion overproduction problem from the saxion decay.