The Wheeler-DeWitt equation provides the probability distribution for the curvature perturbation, the gauge invariant quantum fluctuation of the inflaton. From this, we can find a tower of power spectrum which is not found in a perturbative approach. Since the power spectrum for the modes that cross the horizon contributes to the uncertainty in the classical inflaton displacement, we obtain new conditions for the eternal inflation. In the presence of the patch in the higher excitations, the bound on the slow-roll parameter allowing the eternal inflation is given by at most $\epsilon \lesssim (2n+1)(H/m_{\rm Pl})^2$ with $n$ integer. For large $n$, the bound on $\epsilon$ is relaxed such that the eternal inflation can take place with even larger value of $\epsilon$. While the second law of thermodynamics implies that $n=0$ state is preferred, we cannot ignore such large $n$ effect since the nonlinear interaction inducing transitions to the $n=0$ state is suppressed.

We study the Wigner function for the inflationary tensor perturbation defined in the real phase space. We compute explicitly the Wigner function including the contributions from the cubic self-interaction Hamiltonian of tensor perturbations. Then we argue that it is no longer an appropriate description for the probability distribution in the sense that quantum nature allows negativity around vanishing phase variables. This comes from the non-Gaussian wavefunction in the mixed state as a result of the non-linear interaction between super- and sub-horizon modes. We also show that this is related to the explicit infrared divergence in the Wigner function, in contrast to the trace of the density matrix.

We study the de Sitter (dS) swampland conjecture in the Dirac-Born-Infeld (DBI) inflation model. We obtain the dS swampland bound for the relativistic regime using Bousso's entropy bound argument and proper distance. It restricts $m_{\rm Pl}\nabla V/V$ by some positive constant depending on warping and the field range. In the specific case of the DBI model driven by the quadratic potential, the model-dependent backreaction argument is interpreted as a natural bound for the slow-roll parameter. This shows that quasi-dS spacetime in the DBI model is a result of tuning.

We study the quantum mechanical evolution of the tensor perturbations during inflation with non-linear tensor interactions. We first obtain the Lindblad terms generated by non-linear interactions by tracing out unobservable sub-horizon modes. Then we calculate explicitly the reduced density matrix for the super-horizon modes, and show that the probability of maintaining the unitarity of the squeezed state decreases in time. The decreased probability is transferred to other elements of the reduced density matrix including off-diagonal ones, so the evolution of the reduced density matrix describes the quantum-to-classical transition of the tensor perturbations. This is different from the classicality accomplished by the squeezed state, the suppression of the non-commutative effect, which is originated from the quadratic, linear interaction, and also maintains the unitarity. The quantum-to-classical transition occurs within 5 - 10 e-folds, faster than the curvature perturbation.

We study the consequences of spatial coordinate transformation in multi-field inflation. Among the spontaneously broken de Sitter isometries, only dilatation in the comoving gauge preserves the form of the metric and thus results in quantum-protected Slavnov-Taylor identities. We derive the corresponding consistency relations between correlation functions of cosmological perturbations in two different ways, by the connected and one-particle-irreducible Green's functions. The lowest-order consistency relations are explicitly given, and we find that even in multi-field inflation the consistency relations in the soft limit are independent of the detail of the matter sector.

It is well known that the memory effect in flat spacetime is parametrized by the BMS supertranslation. We investigate the relation between the memory effect and diffeomorphism in de Sitter spacetime. We find that gravitational memory is parametrized by a BMS-like supertranslation in the static patch of de Sitter spacetime. We also show a diffeomorphism that corresponds to gravitational memory in the Poincare/cosmological patch. Our method does not need to assume the separation between the source and the detector to be small compared with the Hubble radius, and can potentially be applicable to other FLRW universes, as well as "ordinary memory" mediated by massive messenger particles.

Thermal production of light dark matter with sub-GeV scale mass can be attributed to $3\rightarrow 2$ self-annihilation processes. We consider the thermal average for annihilation cross sections of dark matter at $3\rightarrow 2$ and general higher-order interactions. A correct thermal average for initial dark matter particles is important, in particular, for annihilation cross sections with overall velocity dependence and/or resonance poles. We apply our general results to benchmark models for SIMP dark matter and discuss the effects of the resonance pole in determining the relic density.

We develop the path integral formalism for studying cosmological perturbations in multi-field inflation, which is particularly well suited to study quantum theories with gauge symmetries such as diffeomorphism invariance. We formulate the gauge fixing conditions based on the Poisson brackets of the constraints, from which we derive two convenient gauges that are appropriate for multi-field inflation. We then adopt the in-in formalism to derive the most general expression for the power spectrum of the curvature perturbation including the corrections from the interactions of the curvature mode with other light degrees of freedom. We also discuss the contributions of the interactions to the bispectrum.

We give a review on the SIMP paradigm and discuss a consistent model for SIMP dark mesons in the context of a dark QCD with flavor symmetry. The $Z'$-portal interaction is introduced being compatible with stable dark mesons and is responsible for making the SIMP dark mesons remain in kinetic equilibrium with the SM during the freeze-out process. The SIMP parameter space of the $Z'$ gauge boson can be probed by future collider and direct detection experiments.

The morphological types of 5840 galaxies were classified by a visual inspection of color images using the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) to produce a morphology catalog of a representative sample of local galaxies with $z<0.01$. The sample galaxies are almost complete for galaxies brighter than $r_{pet}=17.77$. Our classification system is basically the same as that of the Third Reference Catalog of Bright Galaxies with some simplification for giant galaxies. On the other hand, we distinguish fine features of dwarf elliptical-like galaxies to classify 5 subtypes: dwarf ellipticals (dE), blue-cored dwarf ellipticals (dE$_{bc}$), dwarf spheroidals (dSph), blue dwarf ellipticals (dE$_{blue}$), and dwarf lenticulars (dS0). In addition, we denote the presence of nucleation in dE, dSph, and dS0. Elliptical galaxies and lenticular galaxies contribute only $\sim1.5\%$ and $\sim4.9\%$ of the local galaxies, respectively, whereas spirals and irregulars contribute $\sim32.1\%$ and $\sim42.8\%$, respectively. The dE$_{blue}$ galaxies, which are recently found populations of galaxies, contribute a significant fraction to the dwarf galaxies. There seems to be structural difference between dSph and dE galaxies. The dSph galaxies are fainter and bluer with shallower surface brightness gradient than dE galaxies. They also have lower fraction of galaxies with small axis ratios ($b/a \lesssim 0.4$) than dE galaxies. The mean projected distance to the nearest neighbor galaxy is $\sim260$kpc. About $1\%$ of local galaxies have no neighbors with comparable luminosity within a projected distance of 2Mpc.

We study next-to-leading corrections to the effective action of the curvature perturbation obtained by integrating out the coupled heavy isocurvature perturbation. These corrections result from including higher order derivative operators, weighted by the mass scale of the heavy physics, in the effective theory expansion. We find that the correction terms are suppressed by the ratio of the Hubble parameter to the heavy mass scale. The corresponding corrections to the power spectrum of the curvature perturbation are presented for a simple illustrative example.

The QCD axion solving the strong CP problem may originate from antisymmetric tensor gauge fields in compactified string theory, with a decay constant around the GUT scale. Such possibility appears to be ruled out now by the detection of tensor modes by BICEP2 and the PLANCK constraints on isocurvature density perturbations. A more interesting and still viable possibility is that the string theoretic QCD axion is charged under an anomalous U(1)_A gauge symmetry. In such case, the axion decay constant can be much lower than the GUT scale if moduli are stabilized near the point of vanishing Fayet-Illiopoulos term, and U(1)_A-charged matter fields get a vacuum value far below the GUT scale due to a tachyonic SUSY breaking scalar mass. We examine the symmetry breaking pattern of such models during the inflationary epoch with the Hubble expansion rate 10^{14} GeV, and identify the range of the QCD axion decay constant, as well as the corresponding relic axion abundance, consistent with known cosmological constraints. In addition to the case that the PQ symmetry is restored during inflation, there are other viable scenarios, including that the PQ symmetry is broken during inflation at high scales around 10^{16}-10^{17} GeV due to a large Hubble-induced tachyonic scalar mass from the U(1)_A D-term, while the present axion scale is in the range 10^{9}-5\times 10^{13} GeV, where the present value larger than 10^{12} GeV requires a fine-tuning of the axion misalignment angle. We also discuss the implications of our results for the size of SUSY breaking soft masses.

I will talk on my recent works. Axino, related to the SUSY transformation of axion, can mix with Goldstino in principle. In this short talk, I would like to explain what is the axino mass and its plausible mass range. The axino mass is known to have a hierarchical mass structure depending on accidental symmetries. With only one axino, if G_A=0 where G=K+ 2ln|W|, we obtain axino mass= gravitino mass. For G_A nonzero, the axino mass depends on the details of the Kaehler potential. I also comment on the usefulness of a new parametrization of the CKM matrix.

Axino, related to the SUSY transformation of axion, can mix with goldstino in principle. This case is realized when some superfields carrying nonvanishing Peccei-Quinn charges develop both scalar VEVs and F-terms. In this case, we present a proper definition of axion and axino. With this definition, we present the QCD axino mass in the most general framework, including non-minimal K\"ahler potential. The axino mass is known to have a hierarchical mass structure depending on accidental symmetries. With only one axino, if $G_A=0$ where $G=K+\ln|W|^2$, we obtain $m_{\tilde a}=m_{3/2}$. For $G_A\ne 0$, the axino mass depends on the details of the K\"ahler potential. In the gauge mediation scenario, the gaugino mass is the dominant axino mass parameter. Therefore, we can take the theoretical QCD axino mass as a free parameter in the study of its cosmological effects, ranging from eV to multi-TeV scales, without a present knowledge on its ultraviolet completion.