We consider the effects of the injections of energetic photon and electron (or positron) on the big-bang nucleosynthesis. We study the photodissociation of light elements in the early Universe paying particular attention to the case that the injection energy is sub-GeV and derive upper bounds on the primordial abundances of the massive decaying particle as a function of its lifetime. We also discuss a solution of the $^7$Li problem in this framework.

We discuss the footprint of evaporation of primordial black holes (PBHs) on stochastic gravitational waves (GWs) induced by scalar perturbations. We consider the case where PBHs once dominated the Universe but eventually evaporated before the big bang nucleosynthesis. The reheating through the PBH evaporation could end with a sudden change in the equation of state of the Universe compared to the conventional reheating caused by particle decay. We show that this "sudden reheating" by the PBH evaporation enhances the induced GWs, whose amount depends on the length of the PBH-dominated era and the width of the PBH mass function. We explore the possibility to constrain the primordial abundance of the evaporating PBHs by observing the induced GWs. We find that the abundance parameter $\beta \gtrsim 10^{-5} \text{ - }10^{-8}$ for $\mathcal{O}(10^3 \text{ - } 10^5) \, \text{g}$ PBHs can be constrained by future GW observations if the width of the mass function is smaller than about a hundredth of the mass.

Q-balls are non-topological solitons whose classical stability is ensured by a global $U(1)$ charge. In particular, Q-balls are produced in the framework of the Affleck-Dine baryogenesis where $U(1)$ charge is the baryon number. Since some type of Q-balls is stable against quantum decay into nucleons, it can work as the dark matter. When the dark matter Q-balls are captured into the neutron star, they absorb the surrounding neutrons and grow to consume all the neutrons, which leads to a stringent constraint on the Q-ball dark matter. However, the Q-ball growth stops due to the $U(1)$ breaking A-term. We revisit the constraint by the neutron star for the gauge-mediation type and the new type Q-ball, and find the allowed parameter space, in which Q-balls work as the dark matter.

We consider the cosmological effects of sterile neutrinos with the masses of $150- 450$ MeV. The decay of sterile neutrinos changes the thermal history of the Universe and affects the energy density of radiation at the recombination and the big bang nucleosynthesis (BBN) results. We derive severe constraints on the parameters of sterile neutrinos from the primordial abundances of helium-4 and deuterium. We also find that in a particular model the constraints can be considerably relaxed by assuming a large lepton asymmetry in the active neutrinos. In this case, the consistent parameters result in $N_{\mathrm{eff}} \simeq 3.2- 3.4$ and can alleviate the Hubble tension.

We study the observational signatures from particle production of a $U(1)$ gauge field kinetically coupled to an inflaton. Regarding the form of gauge kinetic function, we consider the possibility that it becomes stabilized at a certain time, which makes the growing power of the gauge field evolve non-monotonically with a sharp transition. Remarkably, the copious production of the gauge field occurs on super-horizon scales at the late stage of inflation and perturbations are enhanced on the intermediate scales during inflation. We find that it can predict a bumpy shape of the curvature power spectrum which leads to the generation of primordial black holes as a dark matter after inflation. We also estimate two types of tensor modes sourced by the gauge field: the primordial gravitational waves generated during inflation and the induced gravitational waves provided by the enhanced curvature perturbation after inflation. We show that both of them are potentially testable with the future space-based gravitational wave interferometers.

Ultra-light axion-like particle (ULAP) is one of attractive candidates for cold dark matter. Because the de Broglie wavelength of ULAP with mass $\sim 10^{-22} {\rm eV}$ is $\mathcal{O}({\rm kpc})$, the suppression of the small scale structure by the uncertainty principle can solve the core-cusp problem. Frequently, ULAP is assumed to be uniformly distributed in the present universe. In typical ULAP potentials, however, strong self-resonance at the beginning of oscillation invokes the large fluctuations, which may cause the formation of the dense localized object, oscillon. % Such a dense object lives for a long time, it may affect the cosmological evolution. In this paper, we confirm the oscillon formation in a ULAP potential by numerical simulation and analytically derive its lifetime.

The axion causes two serious cosmological problems, domain wall and isocurvature perturbation problems. Linde pointed that the isocurvature perturbations are suppressed when the Peccei-Quinn (PQ) scalar field takes a large value $\sim M_{\text{pl}}$ (Planck scale) during inflation. In this case, however, the PQ field with large amplitude starts to oscillate after inflation and large fluctuations of the PQ field are produced through parametric resonance, which leads to the formation of domain walls. We consider a supersymmetric axion model and examine whether domain walls are formed by using lattice simulation. It is found that the domain wall problem does not appear in the SUSY axion model when the initial value of the PQ field is less than $10^{3}\times v$ where $v$ is the PQ symmetry breaking scale.

If the Peccei-Quinn symmetry is already broken during inflation, the decay constant $F_a$ of the axion can be in a wide region from $10^{11}$ GeV to $10^{18}$ GeV for the axion being the dominant dark matter. In this case, however, the axion causes the serious cosmological problem, isocurvature perturbation problem, which severely constrains the Hubble parameter during inflation. The constraint is relaxed when Peccei-Quinn scalar field takes a large value $\sim M_{p}$ (Planck scale) during inflation. In this letter, we point out that the allowed region of the decay constant $F_a$ is reduced to a rather narrow region for a given tensor-to-scalar ratio $r$ when Peccei-Quinn scalar field takes $\sim M_{p}$ during inflation. For example, if the ratio $r$ is determined as $r \gtrsim 10^{-3}$ in future measurements, we can predict $F_a \simeq (0.1-1.4)\times 10^{12}$ GeV for domain wall number $N_\text{DW}=6$.

I-ball/oscillon is a soliton-like oscillating configuration of a real scalar field which lasts for a long time. I-ball/oscillon is a minimum energy state for a given adiabatic invariant, and its approximate conservation guarantees the longevity. In this paper, we examine the stability of a special type of I-ball/oscillon, the "exact" I-ball/oscillon, whose adiabatic invariant is exactly conserved. We show that the exact I-ball/oscillon is stable in classical field theory, but not stable against small perturbations depending on the value of its adiabatic invariant. Accordingly, the exact I-ball/oscillon breaks up in the presence of the fluctuations with corresponding instability modes. We also confirm the fragileness of the exact I-ball/oscillon by the classical lattice simulation.

We study the supermassive black holes (SMBHs) observed in the galactic centers. Although the origin of SMBHs has not been well understood yet, previous studies suggest that seed black holes (BHs) with masses $10^{4-5}M_\odot$ exist at a high redshift ($z \sim 10$). We examine whether primordial black holes (PBHs) produced by inhomogeneous baryogenesis can explain those seed black holes. The inhomogeneous baryogenesis is realized in the modified Affleck-Dine mechanism. In this scenario, there is no stringent constraint from CMB $\mu$-distortion in contrast to the scenario where Gaussian fluctuations collapse into PBHs. It is found that the model can account for the origin of the seed BHs of the SMBHs.

I-balls/oscillons are long-lived and spatially localized solutions of real scalar fields. They are produced in various contexts of the early universe in, such as, the inflaton evolution and the axion evolution. However, their decay process has long been unclear. In this paper, we derive an analytic formula of the decay rate of the I-balls/oscillons within the classical field theory. In our approach, we calculate the Poynting vector of the perturbation around the I-ball/oscillon profile by solving a relativistic field equation, with which the decay rate of the I-ball/oscillon is obtained. We also perform a classical lattice simulation and confirm the validity of our analytical formula of the decay rate numerically.

The chemical separation of Li+ ions induced by a magnetic field during the hierarchical structure formation can reduce initial Li abundances in cosmic structures. It is shown that cosmological reionization of neutral Li atoms quickly completes as soon as the first star is formed. Since almost all Li is singly ionized during the main course of structure formation, it can efficiently separate from gravitationally collapsing neutral gas. The separation is more efficient in smaller structures which formed earlier. In the framework of the hierarchical structure formation, extremely metal-poor stars can have smaller Li abundances because of their earlier formations. It is found that the chemical separation by a magnetic field thus provides a reason that Li abundances in extremely metal-poor stars are lower than the Spite plateau and have a large dispersion as well as an explanation of the Spite plateau itself. In addition, the chemical separation scenario can explain Li abundances in NGC 6397 which are higher than the Spite plateau. Thus, Li abundances in metal-poor stars possibly keep information on the primordial magnetic field and the structure formation history.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

We study the effect of the nonlinear relation between density and curvature perturbations on the formation of PBHs. By calculating the variance and skewness of the density perturbation we derive the non-Gaussian property. As a criterion for PBH formation, the compaction function is used and it is found that larger curvature perturbations are required due to the nonlinear effect. We estimate the PBH abundance based on the Press-Schechter formalism with non-Gaussian probability density function during Radiation dominated era. It is found that the nonlinear effect slightly suppresses the PBH formation and the suppression is comparable to that expected if the primordial curvature perturbation would have the local form of non-Gaussianity with nonlinear parameter $f_\text{NL} \sim -1$.

The "swampland conjectures" have been recently suggested as a set of criteria to assess if effective field theories (EFTs) are consistent with a quantum gravity embedding. Such criteria, which restrict the behavior of scalar fields in the theory, have strong implications for cosmology in the early universe. As we demonstrate, they will also have direct consequences for formation of primordial black holes (PBHs) and dark matter (DM).

We study effects of long-lived massive particles, which decay during the big-bang nucleosynthesis (BBN) epoch, on the primordial abundances of light elements. Compared to the previous studies, (i) the reaction rates of the standard BBN reactions are updated, (ii) the most recent observational data of light element abundances and cosmological parameters are used, (iii) the effects of the interconversion of energetic nucleons at the time of inelastic scatterings with background nuclei are considered, and (iv) the effects of the hadronic shower induced by energetic high energy anti-nucleons are included. We compare the theoretical predictions on the primordial abundances of light elements with latest observational constraints, and derive upper bounds on relic abundance of the decaying particle as a function of its lifetime. We also apply our analysis to unstable gravitino, the superpartner of the graviton in supersymmetric theories, and obtain constraints on the reheating temperature after inflation.

We performed the detailed analysis of the primordial black hole (PBH) formation mechanism in an axionlike curvaton model with a coupling to inflaton. The phase direction of the complex scalar works as a curvaton and produces enough PBHs to explain the black hole binaries ($\sim 30 M_{\odot}$) observed in the LIGO-Virgo Collaboration or PBHs as dark matter (DM) ($\sim 10^{-12} M_{\odot}$). We examine whether our model satisfies the current constraints on the PBH mass spectrum, the curvature perturbation and the secondarily produced gravitational waves. We also take into account ambiguity about the choice of the window functions and effect of the non-Gaussianity.

We examine the oscillon formation in a recently proposed inflation model of the pure natural inflation, where the inflaton is an axion that couples to a strongly-coupled pure Yang-Mills theory. The plateau of the inflaton potential, which is favored by recent observations, drives the fragmentation of the inflaton and can produce spatially localized oscillons. We find that the oscillons are formed for F<~ O(0.1) M_{Pl}, with F the effective decay constant of the model. We also comment on observational implications of the oscillons.

We present results of new field-theoretic simulation of cosmological axion strings, which are eight times longer than previous ones. We have upgraded our simulation of physical strings in Hiramatsu et al. (2011) in terms of the number of grids as well as the suite of analysis methods. These improvements enable us to monitor a variety of quantities characterizing the dynamics of the physical string network for the longest term ever. Our extended simulations have revealed that global strings do not evolve according to the scaling solution but its scaling parameter, or the number of long strings per horizon, increases logarithmically in time. In addition, we have also found that the scaling parameter shows nontrivial dependence on the breaking scale of the Peccei-Quinn symmetry.

In this letter, we propose a novel scenario which simultaneously explains $\mathcal{O}(10)M_\odot$ primordial black holes (PBHs) and dark matter in the minimally supersymmetric standard model. Gravitational waves (GWs) events detected by LIGO-Virgo collaboration suggest an existence of black holes as heavy as $\sim 30M_\odot$. In our scenario, as seeds of the PBHs, we make use of the baryon number perturbations which are induced by the special type of Affleck-Dine mechanism. Furthermore, the scenario does not suffer from the stringent constraints from CMB $\mu$-distortion due to the Silk damping and pulsar timing. We find the scenario can explain not only the current GWs events consistently, but also dark matter abundance by the non-topological solitons formed after Affleck-Dine mechanism, called Q-balls.

The recent observations of the gravitational waves (GWs) by LIGO-Virgo collaboration infer the increasing possibility of the primordial black holes (PBHs). Recently it was pointed out that sufficient PBHs are produced by the Affleck-Dine mechanism where inhomogeneous baryogenesis takes place due to change of the Hubble induced mass during and after inflation and forms high baryon bubbles (HBBs). The produced HBBs have large density contrasts through the QCD phase transition or stable Q-ball formation, which leads to formation of the LIGO PBHs. Furthermore, in this model stringent constraints from CMB $\mu$-distortion and pulsar timing array (PTA) experiments are completely absent. In this paper, we study the model in full details based on gravity and gauge mediated supersymmetry breaking scenarios and show that the model can explain the current GWs events evading observational constraints.

We revise primordial black hole (PBH) formation in the axion-like curvaton model and investigate whether PBHs formed in this model can be the origin of the gravtitational wave (GW) signals detected by the Advanced LIGO. In this model, small-scale curvature perturbations with large amplitude are generated, which is essential for PBH formation. On the other hand, large curvature perturbations also become a source of primordial GWs by their second-order effects. Severe constraints are imposed on such GWs by pulsar timing array (PTA) experiments. We also check the consistency of the model with these constraints. In this analysis, it is important to take into account the effect of non-Gaussianity, which is generated easily in the curvaton model. We see that, if there are non-Gaussianities, the fixed amount of PBHs can be produced with a smaller amplitude of the primordial power spectrum.

We study the effect of large baryonic isocurvature perturbations on the abundance of deuterium (D) synthesized in big bang nucleosynthesis (BBN). We found that large baryonic isocurvature perturbations existing at the BBN epoch ($T\sim 0.1\,$MeV) change the D abundance by the second order effect, which, together with the recent precise D measurement, leads to a constraint on the amplitude of the power spectrum of the baryon isocurvature perturbations. The obtained constraint on the amplitude is $\lesssim 0.016 \,(2\sigma)$ for scale $k^{-1} \gtrsim 0.0025\,\text{pc}$. This gives the most stringent one for $0.1\,\text{Mpc}^{-1} \lesssim k \lesssim 4\times 10^8\, \text{Mpc}^{-1}$. We apply the BBN constraint to the relaxation leptogenesis scenario, where large baryon isocurvature perturbations are produced in the last $N_\text{last}$ $e$-fold of inflation, and we obtain a constraint on $N_\text{last}$.

Primordial black holes (PBHs) can be produced by the perturbations that exit the horizon during inflationary phase. While inflation models predict the power spectrum of the perturbations in Fourier space, the PBH abundance depends on the probability distribution function (PDF) of density perturbations in real space. In order to estimate the PBH abundance in a given inflation model, we must relate the power spectrum in Fourier space to the PDF in real space by coarse-graining the perturbations with a window function. However, there are uncertainties on what window function should be used, which could change the relation between the PBH abundance and the power spectrum. This is particularly important in considering PBHs with mass $30 M_\odot$ that account for the LIGO events because the required power spectrum is severely constrained by the observations. In this paper, we investigate how large influence the uncertainties on the choice of a window function have over the power spectrum required for LIGO PBHs. As a result, it is found that the uncertainties significantly affect the prediction for the stochastic gravitational waves (GWs) induced by the second order effect of the perturbations. In particular, the pulsar timing array constraints on the produced GWs could disappear for the real-space top-hat window function.

Primordial Black Hole (PBH) is one of the leading non-particle candidates for dark matter (DM). Although several observations severely constrain the amount of PBHs, it is recently pointed out that there is an uncertainty on the microlensing constraints below $\sim 10^{-10} M_\odot$ which was ignored originally but may weaken the constraints significantly. In this paper, facing this uncertainty, we investigate the possibility that PBHs can make up all DM in a broad mass spectrum. Moreover, we propose a concrete inflation model which can simultaneously produce PBHs for all DM in a broad mass spectrum around $\mathcal O(10^{-13}) M_\odot$ and PBHs for LIGO events in a sharp mass spectrum at $\mathcal O(10) M_\odot$.

The domain wall problem of the Peccei-Quinn mechanism can be solved if the Peccei-Quinn symmetry is explicitly broken by a small amount. Domain walls decay into axions, which may account for dark matter of the universe. This scheme is however strongly constrained by overproduction of axions unless the phase of the explicit breaking term is tuned. We investigate the case where the universe is matter-dominated around the temperature of the MeV scale and domain walls decay during this matter dominated epoch. We show how the viable parameter space is expanded.

LIGO-Virgo collaboration has found black holes as heavy as $M \sim 30M_\odot$ through the detections of the gravitational waves emitted during their mergers. Primordial black holes (PBHs) produced by inflation could be an origin of such events. While it is tempting to presume that these PBHs constitute all Dark Matter (DM), there exists a number of constraints for PBHs with $\mathcal{O} (10) M_\odot$ which contradict with the idea of PBHs as all DM. Also, it is known that weakly interacting massive particle (WIMP) that is a common DM candidate is almost impossible to coexist with PBHs. These observations motivate us to pursue another candidate of DM. In this paper, we assume that the string axion solving the strong CP problem makes up all DM, and discuss the coexistence of string axion DM and inflationary PBHs for LIGO events.

The string theory predicts many light fields called moduli and axions, which cause a cosmological problem due to the overproduction of their coherent oscillation after inflation. One of the prominent solutions is an adiabatic suppression mechanism, which, however, is non-trivial to achieve in the case of axions because it necessitates a large effective mass term which decreases as a function of time. The purpose of this paper is twofold. First, we provide an analytic method to calculate the cosmological abundance of coherent oscillation in a general situation under the adiabatic suppression mechanism. Secondly, we apply our method to some concrete examples, including the one where a string axion acquires a large effective mass due to the Witten effect in the presence of hidden monopoles.

We study the dynamics of the Affleck-Dine field after inflation in more detail. After inflation, the Affleck-Dine field inevitably oscillates around the potential minimum. This oscillation is hard to decay and can cause accidental suppression of the consequential baryon asymmetry. This suppression is most effective for the model with non-renormalizable superpotential $W_{AD}\sim \Phi^4$ ($\Phi$: Affleck-Dine field). It is found that the Affleck-Dine leptogenesis in high-scale inflation, which suffers from serious gravitino overproduction, becomes workable owing to this effect.

Primordial black holes (PBHs) are one of the candidates to explain the gravitational wave (GW) signals observed by the LIGO detectors. Among several phenomena in the early Universe, cosmic inflation is a major example to generate PBHs from large primordial density perturbations. In this paper, we discuss the possibility to interpret the observed GW events as mergers of PBHs which are produced by cosmic inflation. The primordial curvature perturbation should be large enough to produce a sizable amount of PBHs and thus we have several other probes to test this scenario. We point out that the current pulsar timing array (PTA) experiments already put severe constraints on GWs generated via the second-order effects, and that the observation of the cosmic microwave background (CMB) puts severe restriction on its $\mu$ distortion. In particular, it is found that the scalar power spectrum should have a very sharp peak at $k \sim 10^{6}$ Mpc$^{-1}$ to fulfill the required abundance of PBHs while evading constraints from the PTA experiments together with the $\mu$ distortion. We propose a mechanism which can realize such a sharp peak. In the future, simple inflation models that generate PBHs via almost Gaussian fluctuations could be probed/excluded.

Following a new microlensing constraint on primordial black holes (PBHs) with $\sim10^{20}$--$10^{28}\,\mathrm{g}$~[1], we revisit the idea of PBH as all Dark Matter (DM). We have shown that the updated observational constraints suggest the viable mass function for PBHs as all DM to have a peak at $\simeq 10^{20}\,\mathrm{g}$ with a small width $\sigma \lesssim 0.1$, by imposing observational constraints on an extended mass function in a proper way. We have also provided an inflation model that successfully generates PBHs as all DM fulfilling this requirement.

Observations of cosmic infrared background (CIB) radiation exhibit significant fluctuations on small angular scales. A number of explanations have been put forth, but there is currently no consensus on the origin of these large fluctuations. We consider the possibility that small-scale fluctuations in matter-antimatter asymmetry could lead to variations in star formation rates which are responsible for the CIB fluctuations. We show that the recently proposed Higgs relaxation leptogenesis mechanism can produce such small-scale baryonic isocurvature perturbations which can explain the observed excess in the CIB fluctuations.

We revisit the compatibility between the chaotic inflation, which provides a natural solution to the initial condition problem, and the metastable electroweak vacuum, which is suggested by the results of LHC and the current mass measurements of top quark and Higgs boson. It is known that the chaotic inflation poses a threat to the stability of the electroweak vacuum because it easily generates large Higgs fluctuations during inflation or preheating and triggers the catastrophic vacuum decay. In this paper, we propose a simple cosmological solution in which the vacuum is stabilized during chaotic inflation, preheating and after that. This simple solution naturally predicts the formation of primordial black holes. We find interesting parameter regions where the present dark matter density is provided by them. Also, the thermal leptogenesis can be accommodated in our scenario.

We propose a novel scenario to produce abundant primordial black holes (PBHs) in new inflation which is a second phase of a double inflation in the supergravity frame work. In our model, some preinflation phase before the new inflation is assumed and it would be responsible for the primordial curvature perturbations on the cosmic microwave background scale, while the new inflation produces only the small scale perturbations. Our new inflation model has linear, quadratic, and cubic terms in its potential and PBH production corresponds with its flat inflection point. The linear term can be interpreted to come from a supersymmetry-breaking sector, and with this assumption, the vanishing cosmological constant condition after inflation and the flatness condition for the inflection point can be consistently satisfied.

We revisit the constraints on the small scale density perturbations ($10^4\,\mathrm{Mpc}^{-1}\lesssim k \lesssim10^5\,\mathrm{Mpc}^{-1}$) from the modification of the freeze-out value of the neutron-proton ratio at big-bang nucleosynthesis era. Around the freeze-out temperature $T\sim 0.5\,\mathrm{MeV}$, the universe can be divided into several local patches which have different temperatures since any perturbation which enters the horizon after the neutrino decoupling has not diffused yet. Taking account of this situation, we calculate the freeze-out value in detail. We find that the small scale perturbations decrease the n-p ratio in contrast to previous works. With use of the latest observed $^4$He abundance, we obtain the constraint on the power spectrum of the curvature perturbations as $\Delta^2_\mathcal{R}\lesssim 0.018$ on $10^4\,\mathrm{Mpc}^{-1}\lesssim k \lesssim 10^5\,\mathrm{Mpc}^{-1}$.

We studied the possibility whether the massive primordial black holes (PBHs) surviving today can be produced in hybrid inflation. Though it is of great interest since such PBHs can be the candidate for dark matter or seeds of the supermassive black holes in galaxies, there have not been quantitatively complete works yet because of the non-perturbative behavior around the critical point of hybrid inflation. Therefore, combining the stochastic and $\delta N$ formalism, we numerically calculated the curvature perturbations in a non-perturbative way and found, without any specific assumption of the types of hybrid inflation, PBHs are rather overproduced when the waterfall phase of hybrid inflation continues so long that the PBH scale is well enlarged and the corresponding PBH mass becomes sizable enough.

We investigate how well we can constrain the mass of light gravitino m_3/2 by using future observations of 21 cm line fluctuations such as Square Kilometre Array (SKA) and Omniscope. Models with light gravitino with the mass m_3/2 < O(10) eV are quite interesting because they are free from the cosmological gravitino problem and consistent with many baryogenesis/leptogenesis scenarios. We evaluate expected constraints on the mass of light gravitino from the observations of 21 cm line, and show that the observations are quite useful to prove the mass. If the gravitino mass is m_3/2 = 1 eV, we found expected 1 sigma errors on m_3/2 are sigma(m_3/2) = 0.25 eV (SKA phase 1), 0.16 eV (SKA phase 2) and 0.067 eV (Omniscope) in combination with Planck + Simons Array + DESI (BAO) + H_0. Additionally, we also discuss detectability of the effective number of neutrino species by varying the effective number of neutrino species for light gravitino N_3/2 and constraints on the mass of light gravitino in the presence of massive neutrinos. We show that 21 cm line observations can detect the nonzero value of N_3/2, and allow us to distinguish the effects of the light gravitino from those of massive neutrino.

We update the constraint on the dark matter annihilation cross section by using the recent measurements of the CMB anisotropy by the Planck satellite. We fully calculate the cascade of dark matter annihilation products and their effects on ionization, heating and excitation of the hydrogen, hence do not rely on any assumption on the energy fractions that cause these effects.

We consider a cosmologically consistent scenario with a heavy Polonyi field. The Polonyi field with a mass of ${\cal O}(100){\rm\,TeV}$ decays before the Big-Bang Nucleosynthesis (BBN) and avoids the severe constraint from the BBN. However, the abundance of the Lightest Supersymmetric Particle (LSP) produced from the decay often exceeds the observed dark matter density. In our scenario, the dark matter density is obtained by the LSP abundance with an aid of entropy production, and baryon asymmetry is generated by the Affleck-Dine mechanism. We show that the observed baryon-to-dark matter ratio of ${\cal O}(0.1\mathchar`-1)$ is naturally explained in sequestering models with a QCD axion.

We provide an anthropic reason that the supersymmetry breaking scale is much higher than the electroweak scale as indicated by the null result of collider experiments and observed 125 GeV Higgs boson. We focus on a new inflation model as a typical low-scale inflation model that may be expected in the string landscape. In this model, the R-symmetry is broken at the minimum of the inflaton potential and its breaking scale is related to the reheating temperature. Once we admit that the anthropic principle requires thermal leptogenesis, we obtain a lower bound on gravitino mass, which is related to R-symmetry breaking scale. This scenario and resulting gravitino mass predict the consistent amplitude of density perturbations. We also find that string axions and saxions are consistently implemented in this scenario.

We study the Witten effect of hidden monopoles on the QCD axion dynamics, and show that its abundance as well as isocurvature perturbations can be significantly suppressed if there is a sufficient amount of hidden monopoles. When the hidden monopoles make up a significant fraction of dark matter, the Witten effect suppresses the abundance of axion with the decay constant smaller than $10^{12}$ GeV. The cosmological domain wall problem of the QCD axion can also be avoided, relaxing the upper bound on the decay constant when the Peccei-Quinn symmetry is spontaneously broken after inflation.

We investigate the Q-ball decay into the axino dark matter in the gauge-mediated supersymmetry breaking. In our scenario, the Q ball decays mainly into nucleons and partially into axinos to account for the baryon asymmetry and the dark matter of the universe simultaneously. The Q ball decays well before the big bang nucleosynthesis so that it is not affected by the decay. The decay into the supersymmetric particles of the minimal supersymmetric standard model is kinematically prohibited until the very end of the decay, and we could safely make their abundances small enough for the successful big bang nucleosynthesis. We show the regions of axino model parameters and the Q-ball parameters which realize this scenario.

The simplest two-field completion of natural inflation has a regime in which both fields are active and in which its predictions are within the Planck 1-$\sigma$ confidence contour. We show this for the original model of natural inflation, in which inflation is achieved through the explicit breaking of a U(1) symmetry. We consider the case in which the mass coming from explicit breaking of this symmetry is comparable to that from spontaneous breaking, which we show is consistent with a hierarchy between the corresponding energy scales. While both masses are comparable when the observable modes left the horizon, the mass hierarchy is restored in the last e-foldings of inflation, rendering the predictions consistent with the isocurvature bounds. For completeness, we also study the predictions for the case in which there is a large hierarchy of masses and an initial period of inflation driven by the (heavy) radial field.

The next-to-minimal supersymmetric standard model predicts the formation of domain walls due to the spontaneous breaking of the discrete $Z_3$-symmetry at the electroweak phase transition, and they collapse before the epoch of big bang nucleosynthesis if there exists a small bias term in the potential which explicitly breaks the discrete symmetry. Signatures of gravitational waves produced from these unstable domain walls are estimated and their parameter dependence is investigated. It is shown that the amplitude of gravitational waves becomes generically large in the decoupling limit, and that their frequency is low enough to be probed in future pulsar timing observations.

We study the effects of dark-matter annihilation during the epoch of big-bang nucleosynthesis on the primordial abundances of light elements. We improve the calculation of the light-element abundances by taking into account the effects of anti-nucleons emitted by the annihilation of dark matter and the interconversion reactions of neutron and proton at inelastic scatterings of energetic nucleons. Comparing the theoretical prediction of the primordial light-element abundances with the latest observational constraints, we derive upper bounds on the dark-matter pair-annihilation cross section. Implication to some of particle-physics models are also discussed.

Thermal inflation is an attractive idea to dilute cosmic density of unwanted particles such as moduli fields which cause cosmological difficulties. However, it also dilutes preexisting baryon asymmetry and some viable baryogenesis is necessary for a cosmologically consistent scenario. We investigate whether the Affleck-Dine mechanism can produce baryon asymmetry enough to survive after the dilution in gauge-mediated SUSY breaking models. Flat directions except for $LH_u$ flat direction cannot provide such huge baryon number because of Q-ball formation. We show that although the $LH_u$ flat direction is special in terms of having $\mu$-term which prevents Q-ball formation, it cannot explain the observed baryon asymmetry either.

The cosmological scenario where the Peccei-Quinn symmetry is broken after inflation is investigated. In this scenario, topological defects such as strings and domain walls produce a large number of axions, which contribute to the cold dark matter of the universe. The previous estimations of the cold dark matter abundance are updated and refined based on the field-theoretic simulations with improved grid sizes. The possible uncertainties originated in the numerical calculations are also discussed. It is found that axions can be responsible for the cold dark matter in the mass range $m_a=(0.8-1.3)\times 10^{-4}\mathrm{eV}$ for the models with the domain wall number $N_{\rm DW}=1$, and $m_a\approx\mathcal{O}(10^{-4}-10^{-2})\mathrm{eV}$ with a mild tuning of parameters for the models with $N_{\rm DW}>1$. Such higher mass ranges can be probed in future experimental studies.

Real scalar fields are known to fragment into spatially localized and long-lived solitons called oscillons or $I$-balls. We prove the adiabatic invariance of the oscillons/$I$-balls for a potential that allows periodic motion even in the presence of non-negligible spatial gradient energy. We show that such potential is uniquely determined to be the quadratic one with a logarithmic correction, for which the oscillons/$I$-balls are absolutely stable. For slightly different forms of the scalar potential dominated by the quadratic one, the oscillons/$I$-balls are only quasi-stable, because the adiabatic charge is only approximately conserved. We check the conservation of the adiabatic charge of the $I$-balls in numerical simulation by slowly varying the coefficient of logarithmic corrections. This unambiguously shows that the longevity of oscillons/$I$-balls is due to the adiabatic invariance.

The Peccei-Quinn mechanism suffers from the problem of the isocurvature perturbations. The isocurvature perturbations are suppressed if the Peccei-Quinn breaking scale is large during inflation. The oscillation of the Peccei-Quinn breaking field after inflation, however, leads to the formation of domain walls due to the parametric resonance effect. In this paper, we discuss the evolution of the Peccei-Quinn breaking field after inflation in detail, and propose a model where the parametric resonance is ineffective and hence domain walls are not formed. We also discuss consistency of our model with supersymmetric theory.

We point out that the prediction of the minimal chaotic inflation model is altered if a scalar field takes a large field value close to the Planck scale during inflation due to a negative Hubble induced mass. In particular, we show that the inflaton potential is effectively suppressed at a large inflaton field value in the presence of such a scalar field. The scalar field may be identified with the standard model Higgs field or flat directions in supersymmetric theory. With such spontaneous suppression, we find that the minimal chaotic inflation model, especially the model with a quadratic potential, is consistent with recent observations of the cosmic microwave background fluctuation without modifying the inflation model itself.

We propose a QCD axion model where the origin of PQ symmetry and suppression of axion isocurvature perturbations are explained by introducing an extra dimension. Each extra quark-antiquark pair lives on branes separately to suppress PQ breaking operators. The size of the extra dimension changes after inflation due to an interaction between inflaton and a bulk scalar field, which implies that the PQ symmetry can be drastically broken during inflation to suppress undesirable axion isocurvature fluctuations.

A coherently oscillating real scalar field with potential shallower than quadratic one fragments into spherical objects called I-balls. We study the I-ball formation for logarithmic potential which appears in many cosmological models. We perform lattice simulations and find that the I-balls are formed when the potential becomes dominated by the quadratic term. Furthermore, we estimate the I-ball profile assuming that the adiabatic invariant is conserved during formation and obtain the result that agrees to the numerical simulations.

We study Q-ball dark matter in gauge-mediated supersymmetry breaking, and seek the possibility of detection in the IceCube experiment. We find that the Q balls would be the dark matter in the parameter region different from that for gravitino dark matter. In particular, the Q ball is a good dark matter candidate for low reheating temperature, which may be suitable for the Affleck-Dine baryogenesis and/or nonthermal leptogenesis. Dark matter Q balls are detectable by IceCube-like experiments in the future, which is the peculiar feature compared to the case of gravitino dark matter.

It is known that after Affleck-Dine baryogenesis, spatial inhomogeneities of Affleck-Dine field grow into non-topological solitons called Q-balls. In gauge mediated SUSY breaking models, sufficiently large Q-balls with baryon charge are stable while Q-balls with lepton charge can always decay into leptons. For a Q-ball that carries nonzero $B$ and $L$ charges, the difference between the baryonic component and the leptonic component in decay rate may induce nonzero electric charge on the Q-ball. This implies that charged Q-ball, also called gauged Q-ball, may emerge in our universe. In this paper, we investigate two complex scalar fields, a baryonic scalar field and a leptonic one, in an Abelian gauge theory. We find stable solutions of gauged Q-balls for different baryon and lepton charges. Those solutions shows that a Coulomb potential arises and the Q-ball becomes electrically charged as expected. It is energetically favored that some amount of leptonic component decays, but there is an upper bound on its amount due to the Coulomb force. The baryonic decay also becomes possible by virtue of electrical repulsion and we find the condition to suppress it so that the charged Q-balls can survive in the universe.

We propose a QCD axion model that avoids the cosmological domain wall problem, introducing a global SU(3)_f family symmetry to which we embed the unwanted PQ discrete symmetry. The spontaneous breaking of SU(3)_f and PQ symmetry predicts eight NG bosons as well as axion, all of which contribute to dark radiation in the Universe. The derivation from the standard model prediction of dark radiation can be observed by future observations of CMB fluctuations. Our model also predicts a sizable exotic kaon decay rate, which is marginally consistent with the present collider data and would be tested by future collider experiments.

We investigate the scenario that one flat direction creates baryon asymmetry of the unverse, while Q balls from another direction can be the dark matter in the gauge-mediated supersymmetry breaking for high-scale inflation. Isocurvature fluctuations are suppressed by the fact that the Affleck-Dine field stays at around the Planck scale during inflation. We find that the dark matter Q balls can be detected in IceCube-like experiments in the future.

We investigate the Affleck-Dine baryogenesis after D-term inflation with a positive Hubble-induced mass term for a B-L flat direction. It stays at a large field value during D-term inflation, and just after inflation ends it starts to oscillate around the origin of the potential due to the positive Hubble-induced mass term. The phase direction is kicked by higher-dimensional Kahler potentials to generate the B-L asymmetry. The scenario predicts nonzero baryonic isocurvature perturbations, which would be detected by future observations of CMB fluctuations. We also provide a D-term inflation model which naturally explain the coincidence of the energy density of baryon and dark matter.

We consider an initial condition problem in a nearly quadratic chaotic inflation model in supergravity. We introduce shift symmetry breaking not only in the superpotential but also in the Kahler potential. In this model the inflaton potential is nearly quadratic for inflaton field values around the Planck scale, but deviates from the quadratic one for larger field values. As a result, the prediction on the tensor-to-scalar ratio can be smaller than that of a purely quadratic model. Due to the shift symmetry breaking in the Kahler potential, the inflaton potential becomes steep for large inflaton field values, which may prevent inflation from naturally taking place in a closed universe. We estimate an upper bound on the magnitude of the shift symmetry breaking so that inflation takes place before a closed universe with a Planck length size collapses, which yields a lower bound on the tensor-to-scalar ratio, $r \mathop{}_{\textstyle \sim}^{\textstyle >} 0.1$.

During the structure formation, charged and neutral chemical species may have separated from each other at the gravitational contraction in primordial magnetic field (PMF). A gradient in the PMF in a direction perpendicular to the field direction leads to the Lorentz force on the charged species. Resultantly, an ambipolar diffusion occurs, and charged species can move differently from neutral species, which collapses gravitationally during the structure formation. We assume a gravitational contraction of neutral matter in a spherically symmetric structure, and calculate fluid motions of charged and neutral species. It is shown that the charged fluid, i.e., proton, electron and $^7$Li$^+$, can significantly decouple from the neutral fluid depending on the field amplitude. The charged species can, therefore, escape from the gravitational collapse. We take the structure mass, the epoch of the gravitational collapse, and the comoving Lorenz force as parameters. We then identify a parameter region for an effective chemical separation. This type of chemical separation can reduce the abundance ratio of Li/H in early structures because of inefficient contraction of $^7$Li$^+$ ion. Therefore, it may explain Li abundances of Galactic metal-poor stars which are smaller than the prediction in standard big bang nucleosynthesis model. Amplitudes of the PMFs are controlled by a magneto-hydrodynamic turbulence. The upper limit on the field amplitude derived from the turbulence effect is close to the value required for the chemical separation.

Matter isocurvature perturbations are strictly constrained from cosmic microwave background observations. We study a sneutrino curvaton model where both cold dark matter (CDM)/baryon isocurvature perturbations are generated. In our model, total matter isocurvature perturbations are reduced since the CDM/baryon isocurvature perturbations compensate for each other. We show that this model can not only avoid the stringent observational constraints but also suppress temperature anisotropies on large scales, which leads to improved agreement with observations.

We investigate the degeneracy of the isocurvature perturbations and the primordial gravitational waves, by using recent observations of the cosmic microwave background (CMB) reported by Planck and BICEP2 collaborations. We show that the tension in the bound on the tensor-to-scalar ratio $r$ between Planck and BICEP2 can be resolved by introducing the anti-correlated isocurvature perturbations. Quantitatively, we find that with the anti-correlated isocurvature perturbations the constraints on $r$ from Planck alone and BICEP2 results can be consistent at 68 % C.L.

In our previous paper, we have proposed a new algorithm to calculate the power spectrum of the curvature perturbations generated in inflationary universe with use of the stochastic approach. Since this algorithm does not need the perturbative expansion with respect to the inflaton fields on super-horizon scale, it works even in highly stochastic cases. For example, when the curvature perturbations are very large or the non-Gaussianities of the curvature perturbations are sizable, the perturbative expansion may break down but our algorithm enables to calculate the curvature perturbations. We apply it to two well-known inflation models, chaotic and hybrid inflation, in this paper. Especially for hybrid inflation, while the potential is very flat around the critical point and the standard perturbative computation is problematic, we successfully calculate the curvature perturbations.

The discovery of the primordial B-mode polarisation by the BICEP2 experiment indicates inflation with a relatively high energy scale. Taking this indication into account, we propose consistent scenarios to account for the observed baryon and dark matter densities in gravity and gauge mediated supersymmetry breaking models. The baryon asymmetry is explained by the Afflck-Dine mechanism, which requires relatively low reheating temperature to avoid a sizable baryonic isocurvature perturbation. The low reheating temperature then requires non-thermal production of dark matter to account for the correct relic density of dark matter. Our scenarios can account for the observations of baryon and dark matter density in gravity and gauge mediation and predict some parameters, including the mass of dark matter.

We investigate that the two types of the Q balls explain the baryon asymmetry and the dark matter of the universe in the gauge-mediated supersymmetry breaking. The gauge-mediation type Q balls of one flat direction produce baryon asymmetry, while the new type Q balls of another flat direction become the dark matter. We show that the dark matter new type Q balls are free from the neutron star constraint. n=5 gauge mediation type and n=6 new type Q balls are displayed as an example, where the potential is lifted by the superpotential Phi^n. These dark matter Q balls may be detected by future observations, such as in advanced IceCube-like observations.

We comprehensively explore the quadratic curvaton models in the chaotic inflation. In the light of the BICEP2 result $r \approx 0.2$, all model parameters and relevant observables are computed. It is found the curvaton field value is constrained into a narrow range, $\sigma_* = \mathcal{O}(10^{-2}$-$10^{-1})$ and the running of the spectral index is $n_s' \gtrsim -10^{-3}$. We show that if the curvaton is added, the models are heavily degenerated on the $n_s$ - $r$ plane. However, introducing a new plane, the degeneracy can be resolved. To distinguish the curvaton models, precise measurements of not only $r$ but also $n_s'$ and the tensor tilt $n_T$ are required.

If inflation scale is high, light scalars acquire large quantum fluctuations during inflation.If sufficiently long-lived, they will give rise to CDM isocurvature perturbations, which are highly constrained by the Planck data. Focusing on string axions as such light scalars, we show that thermal inflation can provide a sufficiently large entropy production to dilute the CDM isocurvature perturbations. Importantly, efficient dilution is possible for the string axions, because effectively no secondary coherent oscillations are induced at the end of thermal inflation, in contrast to the moduli fields. We also study the viability of the axion dark matter with mass of about 7 keV as the origin of the 3.5 keV X-ray line excess, in the presence of large entropy production.

We show that the baryon-dark matter coincidence problem is solved in the CMSSM. The baryons and dark matter are generated simultaneously through the late-time decay of non-topological solitons, Q-balls, which are formed after the Affleck-Dine baryogenesis. A certain relation between the universal scalar mass, $m_0$, and the universal gaugino mass, $M_{1/2}$, is required to solve the coincidence problem, marginally depending on the other CMSSM parameters, and the result is consistent with the observation of the 126-GeV Higgs boson.

Recently, BICEP2 has reported the large tensor-to-scalar ratio $r=0.2^{+0.07}_{-0.05}$ from the observation of the cosmic microwave background (CMB) B-mode at degree-scales. Since tensor modes induce not only CMB B-mode but also the temperature fluctuations on large scales, to realize the consistent temperature fluctuations with the Planck result we should consider suppression of scalar perturbations on corresponding large scales. To realize such a suppression, we consider anti-correlated iso-curvature perturbations which could be realized in the simple curvaton model.

We estimate dark matter density for the Universe with a reheating temperature smaller than the mass of dark matter, assuming dark matter to be a weakly interacting massive particle. During the reheating process, an inflaton decays and releases high-energy particles, which are scattered inelastically by the thermal plasma and emit many particles. Dark matters are produced through these inelastic scattering processes and pair creation processes by high-energy particles. We properly take account of the Landau-Pomeranchuk-Migdal effect on inelastic processes and show that the resultant energy density of dark matter is much larger than that estimated in the literature and can be consistent with that observed when the mass of dark matter is larger than $O(100)$ GeV.

We revisit the production of gravitational waves from unstable domain walls analyzing their spectrum by the use of field theoretic lattice simulations with grid size $1024^3$, which is larger than the previous study. We have recognized that there exists an error in the code used in the previous study, and the correction of the error leads to the suppression of the spectrum of gravitational waves at high frequencies. The peak of the spectrum is located at the scale corresponding to the Hubble radius at the time of the decay of domain walls, and its amplitude is consistent with the naive estimation based on the quadrupole formula. Using the numerical results, the magnitude and the peak frequency of gravitational waves at the present time are estimated. It is shown that for some choices of parameters the signal of gravitational waves is strong enough to be probed in the future gravitational wave experiments.

I-balls/oscillons are long-lived spatially localized lumps of a scalar field which may be formed after inflation. In the scalar field theory with monomial potential nearly and shallower than quadratic, which is motivated by chaotic inflationary models and supersymmetric theories, the scalar field configuration of I-balls is approximately Gaussian. If the I-ball interacts with another scalar field, the I-ball eventually decays into radiation. Recently, it was pointed out that the decay rate of I-balls increases exponentially by the effects of Bose enhancement under some conditions and a non-perturbative method to compute the exponential growth rate has been derived. In this paper, we apply the method to the Gaussian-type I-ball in 3+1 dimensions assuming spherical symmetry, and calculate the partial decay rates into partial waves, labelled by the angular momentum of daughter particles. We reveal the conditions that the I-ball decays exponentially, which are found to depend on the mass and angular momentum of daughter particles and also be affected by the quantum uncertainty in the momentum of daughter particles.

We investigate the consistency of a scenario in which the baryon asymmetry, dark matters, as well as the cosmic density perturbation are generated simultaneously through the evaporation of primordial black holes (PBHs). This scenario can explain the coincidence of the dark matter and the baryon density of the universe, and is free from the isocurvature perturbation problem. We show that this scenario predicts the masses of PBHs, right-handed neutrinos and dark matters, the Hubble scale during inflation, the non-gaussianity and the running of the spectral index. We also discuss the testability of the scenario by detecting high frequency gravitational waves from PBHs.

We present a novel mechanism to generate the cosmic perturbation from evaporation of primordial black holes. A mass of a field is fluctuated if it is given by a vacuum expectation value of a light scalar field because of the quantum fluctuation during inflation. The fluctuated mass causes variations of the evaporation time of the primordial black holes. Therefore provided the primordial black holes dominate the universe when they evaporate, primordial cosmic perturbations are generated. We find that the amplitude of the large scale curvature perturbation generated in this scenario can be consistent with the observed value. Interestingly, our mechanism works even if all fields which are responsible for inflation and the generation of the cosmic perturbation are decoupled from the visible sector except for the gravitational interaction. An implication to the running spectral index is also discussed.

Axion models have two serious cosmological problems, domain wall and isocurvature perturbation problems. In order to solve these problems we investigate the Linde's model in which the field value of the Peccei-Quinn (PQ) scalar is large during inflation. In this model the fluctuations of the PQ field grow after inflation through the parametric resonance and stable axionic strings may be produced, which results in the domain wall problem. We study formation of axionic strings using lattice simulations. It is found that in chaotic inflation the axion model is free from both the domain wall and the isocurvature perturbation problems if the initial misalignment angle \theta_a is smaller than O(10^{-2}). Furthermore, axions can also account for the dark matter for the breaking scale v \simeq 10^{12-16} GeV and the Hubble parameter during inflation H_{inf} \lesssim 10^{11-12} GeV in general inflation models.

The presence of a hot dark matter component has been hinted at 3 sigma by a combination of the results from different cosmological observations. We examine a possibility that pseudo Nambu-Goldstone bosons account for both hot and cold dark matter components. We show that the QCD axions can do the job for the axion decay constant f_a < O(10^10) GeV, if they are produced by the saxion decay and the domain wall annihilation. We also investigate the cases of thermal QCD axions, pseudo Nambu-Goldstone bosons coupled to the standard model sector through the Higgs portal, and axions produced by modulus decay.

We reinvestigate the scenario that the amount of the baryons and the gravitino dark matter is naturally explained by the decay of the Q balls in the gauge-mediated SUSY breaking. Equipped by the more correct decay rates into gravitinos and baryons recently derived, we find that the scenario with the direct production of the gravitino dark matter from the Q-ball decay works naturally.

We calculate the production of the gravitational waves from a double inflation model with lattice simulations. Between the two inflationary stages, gravitational waves with a characteristic frequency are produced by fluctuations of the scalar fields enhanced through parametric resonance. The wavelength of the produced gravitational waves gets extra redshift during the second inflationary stage and it can be in the observable range for the direct gravitational wave detectors. It is found that there is a possibility for the produced gravitational waves to be detected in the planned experiments.

We propose a new approach for calculating the curvature perturbations produced during inflation in the stochastic formalism. In our formalism, the fluctuations of the e-foldings are directly calculated without perturbatively expanding the inflaton field and they are connected to the curvature perturbations by the $\delta N$ formalism. The result automatically includes the contributions of the higher order perturbations because we solve the equation of motion non-perturbatively. In this paper, we analytically prove that our result (the power spectrum and the nonlinearity parameter) is consistent with the standard result in single field slow-roll inflation. We also describe the algorithm for numerical calculations of the curvature perturbations in more general inflation models.

We investigate the gravitational wave background induced by the first order scalar perturbations in the curvaton models. We consider the quadratic and axion-like curvaton potential which can generate the blue-tilted power spectrum of curvature perturbations on small scales and derive the maximal amount of gravitational wave background today. We find the power spectrum of the induced gravitational wave background has a characteristic peak at the frequency corresponding to the scale reentering the horizon at the curvaton decay, in the case where the curvaton does not dominate the energy density of the Universe. We also find the enhancement of the amount of the gravitational waves in the case where the curvaton dominates the energy density of the Universe. Such induced gravitational waves would be detectable by the future space-based gravitational wave detectors or pulsar timing observations.

The observed baryon and dark matter densities are equal up to a factor of 5. This observation indicates that the baryon asymmetry and dark matter have the same origin. The Affleck-Dine baryogenesis is one of the most promising mechanisms in this context. Q balls, which are often formed in the early Universe associated with the Affleck-Dine baryogenesis, decay both into supersymmetric particles and into quarks. Recently, it was pointed out that annihilation of squarks into quarks gives a dominant contribution to the Q-ball decay rate and the branching ratio of Q-ball decay into supersymmetric particles changes from the previous estimate. In this paper, the scenario of baryon and dark matter cogenesis from Q ball in gravity mediation is revisited in respect of the improved Q-ball decay rates. It is found that the successful cogenesis takes place when a wino with mass 0.3-1 TeV is dark matter.

It is known that supersymmetric hybrid inflation model may require severe tunings on the initial condition for large gravitino mass of order 100 - 1000 TeV due to the constant term in the superpotential. We propose a modified hybrid inflation model, where the constant term is suppressed during inflation and generated after inflation by replacing a constant term with dynamical field. In this modified model, successful inflation consistent with large gravitino mass takes place without severe tunings on the initial condition. Constraint from cosmic strings is also relaxed.

We study the CMB constraint on non-Gaussianity in CDM isocurvature perturbations. Non-Gaussian isocurvature perturbations can be produced in various models at the very early stage of the Universe. Since the isocurvature perturbations little affect the structure formation at late times, CMB is the best probe of isocurvature non-Gaussianity at least in the near future. In this paper, we focus on uncorrelated isocurvature perturbations and constrain their non-Gaussianity. For this purpose, we employ several state-of-art techniques for the analysis of CMB data and simulation. We use the WMAP 7 year data of temperature anisotropy. When the adiabatic perturbations are assumed to be Gaussian, we obtained a constraint on the isocurvature non-Gaussianity alpha^2 f_{NL}^{(ISO)}=40+-66 for the scale invariant isocurvature power spectrum, where alpha is the ratio of the power spectrum of isocurvature perturbations to that of the adiabatic ones. When we assume that the adiabatic perturbations can also be non-Gaussian, we obtain f_{NL}=38+-24 and alpha^2 f_{NL}^{(ISO)}=-8+-72. We also discuss implications our results for the axion CDM isocurvature model.

We evaluate the effective mass of a scalar field phi coupled to thermal plasma through Planck-suppressed interactions. We find it useful to rescale the coupled fields so that all the phi-dependences are absorbed into the yukawa and gauge couplings, which allows us to read off the leading order contributions to the effective mass \tilde m_{phi} from the 2-loop free energy calculated with the rescaled couplings. We give an analytical expression for \tilde m_{phi} at a sufficiently high temperature in the case where phi is coupled to the MSSM chiral superfields through non-minimal Kahler potential. We find that \tilde m_{phi}^2 is about 10^{-3} H^2 \sim 10^{-2} H^2 for typical parameter sets, where H is the Hubble expansion rate in the radiation-dominated era.

We argue that the existence of the cold dark matter is explained by primordial black holes.We show that a significant number of primordial black holes can be formed in an axion-like curvaton model, in which the highly blue-tilted power spectrum of primordial curvature perturbations is achieved.It is found that the produced black holes with masses $\sim 10^{20} -10^{38} \mathrm{g}$ account for the present cold dark matter.We also argue the possibility of forming the primordial black holes with mass $\sim 10^5 M_{\odot}$ as seeds of the supermassive black holes.

Spectroscopic observations of metal-poor stars have indicated possible 6Li abundances that are much larger than the primordial abundance predicted in the standard big bang nucleosynthesis model. Possible mechanisms of 6Li production in metal-poor stars include pregalactic and cosmological cosmic ray (CR) nucleosynthesis and nucleosynthesis by flare-accelerated nuclides. We study 9Be production via two-step alpha-fusion reactions of CR or flare-accelerated 3,4He through 6He and 6,7Li, in pregalactic structure, intergalactic medium, and stellar surfaces. We solve transfer equations of CR or flare particles and calculate nuclear yields of 6He, 6,7Li, and 9Be taking account of probabilities of processing 6He and 6,7Li into 9Be via fusions with alpha particles. Yield ratios, i.e., 9Be/6Li, are then calculated for the CR and flare nucleosynthesis models. We suggest that the future observations of 9Be in metal-poor stars may find enhanced abundances originating from metal-poor CR or flare activities.

We study non-Gaussianity of density perturbations generated by an axionic curvaton, focusing on the case that the curvaton sits near the hilltop of the potential during inflation. Such hilltop curvatons can generate a red-tilted density perturbation spectrum without invoking large-field inflation. We show that, even when the curvaton dominates the Universe, the non-Gaussianity parameter fNL is positive and mildly increases towards the hilltop of the curvaton potential, and that fNL = O(10) is a general and robust prediction of such hilltop axionic curvatons. In particular, we find that the non-Gaussianity parameter is bounded as fNL <~ 30 - 40 for a range of the scalar spectral index, ns = 0.94 - 0.99, and that fNL = 20 - 40 is realized for the curvaton mass m_\sigma = 10 - 10^6 GeV and the decay constant f = 10^{12} - 10^{17} GeV. One of the plausible candidates for the axionic curvaton is an imaginary component of a modulus field with mass of order 10 - 100 TeV and decay constant of 10^{16} - 10^{17} GeV. We also discuss extreme cases where the curvaton drives a second inflation and find that fNL is typically smaller compared to non-inflating cases.

We investigate the cosmological constraints on axion models where the domain wall number is greater than one. In these models, multiple domain walls attached to strings are formed, and they survive for a long time. Their annihilation occurs due to the effects of explicit symmetry breaking term which might be raised by Planck-scale physics. We perform three-dimensional lattice simulations and compute the spectra of axions and gravitational waves produced by long-lived domain walls. Using the numerical results, we estimated relic density of axions and gravitational waves. We find that the existence of long-lived domain walls leads to the overproduction of cold dark matter axions, while the density of gravitational waves is too small to observe at the present time. Combining the results with other observational constraints, we find that the whole parameter region of models are excluded unless an unacceptable fine-tuning exists.

We show that the smooth hybrid inflation is naturally realized in a framework of supersymmetric axion model. Identifying the Peccei-Quinn scalar fields as a part of the infaton sector, successful inflation takes place reproducing the amplitude and spectral index of the curvature perturbation observed by WMAP. A relatively large axion isocurvature perturbation and its non-Gaussianity are predicted in our model. The saxion coherent oscillation has a large amplitude and dominates the Universe. The subsequent decay of the saxion produces huge amount of entropy, which dilutes unwanted relics. Winos, the lightest supersymmetric particles in this scenario, are produced non-thermally in the decay and account for dark matter.

We review recent developments on axion cosmology. Topics include : axion cold dark matter, axions from topological defects, axion isocurvature perturbation and its non-Gaussianity and axino/saxion cosmology in supersymmetric axion model.

We study CMB constraints on non-Gaussianity from isocurvature perturbations of general types. Specifically, we study CDM/neutrino isocurvature perturbations which are uncorrelated or totally correlated with adiabatic ones. Using the data from the WMAP 7-year observation at V and W bands, we obtained optimal constraints on the nonlinearity parameters of adiabatic and isocurvature perturbations. Our result shows that primordial perturbations are consistent with Gaussian ones at around 2 sigma level for above mentioned isocurvature modes.

We propose a curvaton model in which the initial condition of the curvaton oscillation is determined by its attractor behavior during inflation. Assuming a chaotic inflation model, we find that the initial condition determined by the attractor behavior is appropriate to generate a sizable non-Gaussianity contribution to the curvature perturbation, which will be tested in the foreseeable future. Implications on the thermal history of the universe and on particle physics models are also discussed.

The recently observed mass ~ 125 GeV for the Higgs boson suggests a high-energy scale SUSY breaking, above O(10) TeV. It is, however, very puzzling why nature chooses such a high energy scale for the SUSY breaking, if the SUSY is a solution to the hierarchy problem. We show that the pure gravity mediation provides us with a possible solution to this puzzle if the topological inflation is the last inflation in the early universe. We briefly discuss a chaotic inflation model in which a similar solution can be obtained.

We evaluate the effective mass of a scalar field which interacts with visible sector via Planck-suppressed coupling in supergravity framework. We focus on the radiation-dominated (RD) era after inflation and the contribution from a fermionic field in the thermal bath. We find that, in RD era, the fermion kinetic term gives the effective mass of the order of Hubble scale to the scalar field.

We analyze the spectrum of axions radiated from collapse of domain walls, which have received less attention in the literature. The evolution of topological defects related to the axion models is investigated by performing field-theoretic lattice simulations. We simulate the whole process of evolution of the defects, including the formation of global strings, the formation of domain walls and the annihilation of the defects due to the tension of walls. The spectrum of radiated axions has a peak at the low frequency, which implies that axions produced by the collapse of domain walls are not highly relativistic. We revisit the relic abundance of cold dark matter axions and find that the contribution from the decay of defects can be comparable with the contribution from strings. This result leads to a more severe upper bound on the axion decay constant.

The Affleck-Dine mechanism, which is one of the most attractive candidates for the baryogenesis in supersymmetric theories, often predicts the existence of baryonic Q balls in the early universe. In this scenario, there is a possibility to explain the observed baryon-to-dark matter ratio because Q balls decay into supersymmetric particles as well as into quarks. If the gravitino mass is small compared to the typical interaction energy, the longitudinal component of the gravitino behaves like the massless goldstino. We numerically calculate the goldstino production rates from Q balls in the leading semi-classical approximation without using large radius limit or effective coupling. We also calculate the quark production rates from Q balls in the Yukawa theory with a massive fermion. In deriving the decay rate we also take into account the scalar field configuration of the Q ball. These results are applied to a realistic model in the gauge-mediated supersymmetry breaking and yield the branching ratio of the Q ball decay into the gravitino. We obtain the branching ratio much smaller than the one estimated in the previous analysis.

We revisit the time evolution of a flat and non-flat direction system during inflation. In order to take into account quantum noises in the analysis, we base on stochastic formalism and solve coupled Langevin equations numerically. We focus on a class of models in which tree-level Hubble-induced mass is not generated. Although the non-flat directions can block the growth of the flat direction's variance in principle, the blocking effects are suppressed by the effective masses of the non-flat directions. We find that the fate of the flat direction during inflation is determined by one-loop radiative corrections and non-renormalizable terms as usually considered, if we remove the zero-point fluctuation from the noise terms.

We study non-Gaussian properties of the isocurvature perturbations in the dark radiation, which consists of the active neutrinos and extra light species, if exist. We first derive expressions for the bispectra of primordial perturbations which are mixtures of curvature and dark radiation isocurvature perturbations. We also discuss CMB bispectra produced in our model and forecast CMB constraints on the nonlinearity parameters based on the Fisher matrix analysis. Some concrete particle physics motivated models are presented in which large isocurvature perturbations in extra light species and/or the neutrino density isocurvature perturbations as well as their non-Gaussianities may be generated. Thus detections of non-Gaussianity in the dark radiation isocurvature perturbation will give us an opportunity to identify the origin of extra light species and lepton asymmetry.

A new upper limit on the amplitude of primordial magnetic field (PMF) is derived by a comparison between a calculation of elemental abundances in big bang nucleosynthesis (BBN) model and the latest observational constraints on the abundances. Updated nuclear reaction rates are adopted in the calculation. Effects of PMF on the abundances are consistently taken into account in the numerical calculation with the precise formulation of changes in physical variables. We find that abundances of 3He and 6Li increase while that of 7Li decreases when the PMF amplitude increases, in the case of the baryon-to-photon ratio determined from the measurement of cosmic microwave background radiation. We derive a constraint on the present amplitude of PMF, i.e., B(0)<1.5 micro G [corresponding to the amplitude less than 2.0x10^{11} G at BBN temperature of T=10^9 K] based on the rigorous calculation.

We reconsider the effective mass of a scalar field which interact with visible sector via Planck-suppressed coupling in supergravity framework. We focus on the radiation-dominated (RD) era after inflation. In this era, the effective mass is given by thermal average of interaction terms. To make our analysis clear, we rely on Kadanoff-Baym equations to evaluate the thermal average. We find that, in RD era, a scalar field acquires the effective mass of the order of $H$.

Supermassive black holes exist in the centers of galaxies, including Milky Way, but there is no compelling theory of their formation. Furthermore, observations of quasars imply that supermassive black holes have already existed at some very high redshifts, suggesting the possibility of their primordial origin. In a class of well-motivated models, inflationary epoch could include two or more periods of inflation dominated by different scalar fields. The transition between such periods of inflation could enhance the spectrum of density perturbations on some specific scale, which could lead to formation of primordial black holes with a very narrow range of masses of the order of 10^5 solar masses. These primordial black holes could have provided the requisite seeds for the observed population of supermassive black holes.

Recent cosmological observations, including measurements of the CMB anisotropy and the primordial helium abundance, indicate the existence of an extra radiation component in the Universe beyond the standard three neutrino species. In this paper we explore the possibility that the extra radiation has isocurvatrue fluctuations. A general formalism to evaluate isocurvature perturbations in the extra radiation is provided in the mixed inflaton-curvaton system, where the extra radiation is produced by the decay of both scalar fields. We also derive constraints on the abundance of the extra radiation and the amount of its isocurvature perturbation. Current observational data favors the existence of an extra radiation component, but does not indicate its having isocurvature perturbation. These constraints are applied to some particle physics motivated models. If future observations detect isocurvature perturbations in the extra radiation, it will give us a hint to the origin of the extra radiation.

We investigate density perturbations sourced by a curvaton with a generic energy potential. The key feature of a curvaton potential which deviates from a quadratic is that the curvaton experiences a non-uniform onset of its oscillation. This sources additional contributions to the resulting density perturbations, and we especially find that the non-Gaussianity parameter f_NL can become large with either sign no matter whether the curvaton dominates or subdominates the universe when it decays. Such non-quadratic curvaton potentials are required in order to produce a red-tilted density perturbation spectrum (without invoking large-field inflation), and are also motivated by explicit curvaton models based on microscopic physics. We further apply our generic results to the case where the curvaton is a pseudo-Nambu-Goldstone (NG) boson of a broken U(1) symmetry, and show that the resulting density perturbations are strongly enhanced towards the hilltop region of the energy potential, accompanied by a mild increase of the non-Gaussianity. Such hilltop NG curvatons can produce observationally suggested density perturbations under wide ranges of inflation/reheating scales, and further predict the non-Gaussianity of the density perturbations to lie within the range 10 <~ f_NL <~ 30.

We discuss how one can discriminate models with cold dark matter (CDM) and baryon isocurvature fluctuations. Although current observations such as cosmic microwave background (CMB) can severely constrain the fraction of such isocurvature modes in the total density fluctuations, CMB cannot differentiate CDM and baryon modes by the shape of its power spectra. However, the evolution of CDM and baryon density fluctuations are different for each model, thus it would be possible to discriminate those isocurvature modes by extracting information on the fluctuations of CDM/baryon itself. We discuss that observations of 21 cm fluctuations can in principle differentiate these modes and demonstrate to what extent we can distinguish them with future 21 cm surveys. We show that, when the isocurvature mode has a large blue-tilted initial spectrum, 21 cm survey can clearly probe the difference.

We study the cosmological evolution of domain walls bounded by strings which arise naturally in axion models. If we introduce a bias in the potential, walls become metastable and finally disappear. We perform two dimensional lattice simulations of domain wall networks and estimate the decay rate of domain walls. By using the numerical results, we give a constraint for the bias parameter and the Peccei-Quinn scale. We also discuss the possibility to probe axion models by direct detection of gravitational waves produced by domain walls.

In this paper, following the previous study, we evaluate the spectrum of gravitational wave background generated by domain walls which are produced if some discrete symmetry is spontaneously broken in the early universe. We apply two different methods to calculate the gravitational wave spectrum: One is to calculate the gravitational wave spectrum directly from numerical simulations, and another is to calculate it indirectly by estimating the unequal time anisotropic stress power spectrum of the scalar field. Both analysises indicate that the slope of the spectrum changes at two characteristic frequencies corresponding to the Hubble radius at the decay of domain walls and the width of domain walls, and that the spectrum between these two characteristic frequencies becomes flat or slightly red tilted. The second method enables us to evaluate the GW spectrum semi-analytically for the frequencies which can not be resolved in the finite box lattice simulations, but relies on the assumptions for the unequal time correlations of the source.

We study the excitation of fine structure levels of C I, C II and O I by ultraviolet (UV) photons around strong UV sources which are also ionizing sources of the cosmological reionization at redshift of $\sim 10$. The evolutions of ionized regions around a point source are calculated by solving rate equations for non-equilibrium chemistry. Signals of UV photons through the fine structure lines are considered to be stronger at locations of more abundant chemical species of C I, C II and O I. Such environments would be realized where strong fluxes of non-ionizing UV line photons available for the pumping up of fine structure levels exist, and simultaneously ionizing UV photons are effectively shielded by dense H I regions. Signals from H I regions of moderately large densities induced by redshifted UV photons emitted at the point sources are found to be dominantly large over those of others. We discuss the detectability of the signals, and show that signals from idealized environments will be possibly detected by radio observations with next-generation arrays to come after the Atacama Large Millimeter/submillimeter Array (ALMA).

We investigate the Q-ball decay in the gauge-mediated SUSY breaking. Q balls decay mainly into nucleons, and partially into gravitinos, while they are kinematically forbidden to decay into sparticles which would be cosmologically harmful. This is achieved by the Q-ball charge small enough to be unstable for the decay, and large enough to be protected kinematically from unwanted decay channel. We can then have right amounts of the baryon asymmetry and the dark matter of the universe, evading any astrophysical and cosmological observational constraints such as the big bang nucleosynthesis, which has not been treated properly in the literatures.

We show that the hybrid inflation is naturally realized in the framework of a supersymmetric axion model, which is consistent with the WMAP observation if the Peccei-Quinn symmetry breaking scale is around 10^{15}GeV. By solving the post inflationary scalar dynamics, it is found that the scalar partner of the axion, saxion, oscillates with large amplitude and its decay produces a huge entropy and dilutes the axion. As a result, the axion coherent oscillation can be the dominant component of the dark matter in the Universe. Cosmological gravitino and axino problems are solved.

Cosmological evolution of axionic string network is analyzed in terms of field-theoretic simulations in a box of 512^3 grids, which are the largest ever, using a new and more efficient identification scheme of global strings. The scaling parameter is found to be \xi=0.87 +- 0.14 in agreement with previous results. The energy spectrum is calculated precisely using a pseudo power spectrum estimator which significantly reduces the error in the mean reciprocal comoving momentum. The resultant constraint on the axion decay constant leads to f_a <= 3*10^11 GeV. We also discuss implications for the early Universe.

In supersymmetric theories, cosmic strings produced in the early Universe often have a width of TeV scale, while the tension is much larger. In a scaling regime, an infinite cosmic string releases significant fraction of its energy in the form of string loops. These thick string loops lose their energies efficiently by particle emissions, and hence it may have effects on cosmological observations. We study cosmological implications of string loops with TeV scale width in detail and derive constraints on the tension of the string. Implications on future gravitational wave detectors are also discussed.

We derive cosmological constraints on the annihilation cross section of dark matter with velocity-dependent structure, motivated by annihilating dark matter models through Sommerfeld or Breit-Wigner enhancement mechanisms. In models with annihilation cross section increasing with decreasing dark matter velocity, big-bang nucleosynthesis and cosmic microwave background give stringent constraints.

We identify reactions which destroy 7Be and 7Li during big bang nucleosynthesis (BBN) in the scenario of BBN catalyzed by a long-lived sub-strongly interacting massive particle (sub-SIMP or X particle). The destruction associated with non radiative X captures of the nuclei can be realized only if the interaction strength between an X particle and a nucleon is properly weaker than that between two nucleons to a degree depending on the mass of X. Binding energies of nuclei to an X particle are estimated taking the mass and the interaction strength to nuclei of the X as input parameters. Nuclear reaction rates associated with the X are estimated naively, and adopted in calculating evolutions of nuclear abundances. We suggest that the 7Li problem, which might be associated with as-yet-unrecognized particle processes operating during BBN, can be solved if the X particle interacts with nuclei strongly enough to drive 7Be destruction but not strongly enough to form a bound state with 4He of relative angular momentum L=1. Justifications of this scenario by rigorous calculations of reaction rates using quantum mechanical many-body models are highly desirable since this result involves many significant uncertainties.

We show that double inflation is naturally realized in K\"ahler moduli inflation, which is caused by moduli associated with string compactification. We find that there is a small coupling between the two inflatons which leads to amplification of perturbations through parametric resonance in the intermediate stage of double inflation. This results in the appearance of a peak in the power spectrum of the primordial curvature perturbation. We numerically calculate the power spectrum and show that the power spectrum can have a peak on observationally interesing scales. We also compute the TT-spectrum of CMB based on the power spectrum with a peak and see that it better fits WMAP 7-years data.

We show that a supersymmetric axion model naturally induces a hybrid inflation with the waterfall field identified as a Peccei-Quinn scalar. The Peccei-Quinn scale is predicted to be around 10^{15}GeV for reproducing the large-scale density perturbation of the Universe. After the built-in late-time entropy-production process, the axion becomes a dark matter candidate. Several cosmological implications are discussed.

We investigate the effect of the stochastic gravitational wave (GW) background produced by kinks on infinite cosmic strings, whose spectrum was derived in our previous work, on the B-mode power spectrum of the cosmic microwave background (CMB) anisotropy. We find that the B-mode polarization due to kinks is comparable to that induced by the motion of the string network and hence the contribution of GWs from kinks is important for estimating the B-mode power spectrum originating from cosmic strings. If the tension of cosmic strings \mu is large enough i.e., G\mu >~ 10^{-8}, B-mode polarization induced by cosmic strings can be detected by future CMB experiments.

We study Q-ball formation in the expanding universe on 1D, 2D and 3D lattice simulations. We obtain detailed Q-ball charge distributions, and find that the distribution is peaked at Q^{3D}_{peak} \simeq 1.9\times 10^{-2} (|\Phi_{in}|/m)^2, which is greater than the existing result by about 60%. Based on the numerical simulations, we discuss how the Q-ball formation proceeds. Also we make a comment on possible deviation of the charge distributions from what was conjectured in the past.

We study the production of gravitational waves from cosmic domain walls created during phase transition in the early universe. We investigate the process of formation and evolution of domain walls by running three dimensional lattice simulations. If we introduce an approximate discrete symmetry, walls become metastable and finally disappear. We calculate the spectrum of gravitational waves produced by collapsing metastable domain walls. Extrapolating the numerical results, we find the signal of gravitational waves produced by domain walls whose energy scale is around 10^10-10^12GeV will be observable in the next generation gravitational wave interferometers.

Gravitational waves emitted by kinks on infinite strings are investigated using detailed estimations of the kink distribution on infinite strings. We find that gravitational waves from kinks can be detected by future pulsar timing experiments such as SKA for an appropriate value of the the string tension, if the typical size of string loops is much smaller than the horizon at their formation. Moreover, the gravitational wave spectrum depends on the thermal history of the Universe and hence it can be used as a probe into the early evolution of the Universe.

We argue that a primordial black hole is a natural and unique candidate for all dark matter. We show that, in a smooth-hybrid new double inflation model, a right amount of the primordial black holes, with a sharply-defined mass, can be produced at the end of the smooth-hybrid regime, through preheating. We first consider masses < 10^(-7)M_sun which are allowed by all the previous constraints. We next discuss much heavier mass 10^5 M_sun hinted at by entropy, and galactic size evolution, arguments. Effects on the running of the scalar spectral index are computed.

We investigate constraints on power spectra of the primordial curvature and tensor perturbations with priors based on single-field slow-roll inflation models. We stochastically draw the Hubble slow-roll parameters and generate the primordial power spectra using the inflationary flow equations. Using data from recent observations of CMB and several measurements of geometrical distances in the late Universe, Bayesian parameter estimation and model selection are performed for models that have separate priors on the slow-roll parameters. The same analysis is also performed adopting the standard parameterization of the primordial power spectra. We confirmed that the scale-invariant Harrison-Zel'dovich spectrum is disfavored with increased significance from previous studies. While current observations appear to be optimally modeled with some simple models of single-field slow-roll inflation, data is not enough constraining to distinguish these models.

We study the light-element abundances in supersymmetric model where the right-handed sneutrino is the lightest superparticle (LSP), assuming that the neutrino masses are purely Dirac-type. In such a scenario, the lightest superparticle in the minimal supersymmetric standard model sector (which we call MSSM-LSP) becomes long-lived, and thermal relic MSSM-LSP may decay after the big-bang nucleosynthesis starts. We calculate the light-element abundances including non-standard nuclear reactions induced by the MSSM-LSP decay, and derive constraints on the scenario of right-handed sneutrino LSP.

We consider primordial fluctuations in thermal inflation scenario. Since the thermal inflation drives about 10 $e$-folds after the standard inflation, the time of horizon-exit during inflation corresponding to the present observational scale shifts toward the end of inflation. It generally makes the primordial power spectrum more deviated from a scale-invariant one and hence renders some models inconsistent with observations. We present a mechanism of generating the primordial curvature perturbation at the end of thermal inflation utilizing a fluctuating coupling of a flaton field with the fields in thermal bath. We show that, by adopting the mechanism, some inflation models can be liberated even in the presence of the thermal inflation. We also discuss non-Gaussianity in the mechanism and show that large non-Gaussianity can be generated in this scenario.

Non-Gaussianity may exist in the CDM isocurvature perturbation. We provide general expressions for the bispectrum and trispectrum of both adiabatic and isocurvature pertubations. We apply our result to the QCD axion case, and found a consistency relation between the coefficients of the bispectrum and trispectrum : tau_{NL}^(iso)~10^3 [f_{NL}^(iso)]^{4/3}, if the axion is dominantly produced by quantum fluctuation. Thus future observations of the trispectrum, as well as the bispectrum, will be important for understanding the origin of the CDM and baryon asymmetry.

We investigate the possibilities of constraining the light gravitino mass m_{3/2} from future cosmic microwave background (CMB) surveys. A model with light gravitino with the mass m_{3/2}<O(10) eV is of great interest since it is free from the cosmological gravitino problem and, in addition, can be compatible with many baryogenesis/leptogenesis scenarios such as the thermal leptogenesis. We show that the lensing of CMB anisotropies can be a good probe for m_{3/2} and obtain an expected constraint on m_{3/2} from precise measurements of lensing potential in the future CMB surveys, such as the PolarBeaR and CMBpol experiments. If the gravitino mass is m_{3/2}=1 eV, we will obtain the constraint for the gravitino mass as m_{3/2} < 3.2 eV (95% C.L.) for the case with Planck+PolarBeaR combined and m_{3/2}=1.04^{+0.22}_{-0.26} eV (68% C.L.) for CMBpol. The issue of Bayesian model selection is also discussed.

We construct an axion model for generating isocurvature fluctuations with blue spectrum, n_{iso}=2-4, which is suggested by recent analyses of admixture of adiabatic and isocurvature perturbations with independent spectral indices, n_{ad} \ne n_{iso}. The distinctive feature of the model is that the spectrum is blue at large scales while scale invariant at small scales. This is naturally realized by the dynamics of the Peccei-Quinn scalar field.

We study the effects of dark matter annihilation during and after the cosmic recombination epoch on the cosmic microwave background anisotropy, taking into account the detailed energy deposition of the annihilation products. It is found that a fairly stringent constraint on the annihilation cross section is imposed for TeV scale dark matter masses from WMAP observations.

The paper by Sato and Kobayashi in 1977 studied the cosmological effects of a massive neutrino and obtained constraints on its properties. This paper initiated many studies to use cosmology as a laboratory of particle physics or to use particle physics to explore the very early universe.

We study the annihilating dark matter contribution to the extra-galactic diffuse gamma-ray background spectrum, motivated by the recent observations of cosmic-ray positron/electron anomalies. The observed diffuse gamma-ray flux provides stringent constraint on dark matter models and we present upper bounds on the annihilation cross section of the dark matter. It is found that for the case of cored dark matter halo profile, the diffuse gamma-rays give more stringent bound compared with gamma-rays from the Galactic center. The Fermi satellite will make the bound stronger.

Recent measurements of cosmic-ray electron and positron fluxes by PAMELA and ATIC experiments may indicate the existence of annihilating dark matter with large annihilation cross section. We show that the dark matter annihilation in the big-bang nucleosynthesis epoch affects the light element abundances, and it gives stringent constraints on such annihilating dark matter scenarios for the case of hadronic annihilation. Constraints on leptonically annihilating dark matter models are less severer.

The positron excess observed by the PAMELA experiment may come from dark matter annihilation, if the annihilation cross section is large enough. We show that the dark matter annihilation scenarios to explain the positron excess may also be compatible with the discrepancy of the cosmic lithium abundances between theory and observations. The wino-like neutralino in the supersymmetric standard model is a good example for it. This scenario may be confirmed by Fermi satellite experiment.

The excess of cosmic-ray electron and positron fluxes measured by the PAMELA satellite and ATIC balloon experiments may be interpreted as the signals of the dark matter annihilation or decay into leptons. In this letter we show that the dark matter annihilation/decay which reproduces the electron/positron excess may yield a significant amount of high-energy neutrinos from the Galactic center. In the case, future kilometer-square size experiments may confirm such a scenario, or even the Super-Kamiokande results already put constraints on some dark matter models.

Light scalars may be ubiquitous in nature, and their quantum fluctuations can produce large non-Gaussianity in the cosmic microwave background temperature anisotropy. The non-Gaussianity may be accompanied with a small admixture of isocurvature perturbations, which often have correlations with the curvature perturbations. We present a general method to calculate the non-Gaussianity in the adiabatic and isocurvature perturbations with and without correlations, and see how it works in several explicit examples. We also show that they leave distinct signatures on the bispectrum of the cosmic microwave background temperature fluctuations.

We study non-Gaussianity induced by a pseudo Nambu-Goldstone boson with a cosine-type scalar potential. We focus on how the non-Gaussianity is affected when the pseudo Nambu-Goldstone boson rolls down from near the top of the scalar potential where the deviation from a quadratic potential is large. We find that the resultant non-Gaussianity is similar to that obtained in the quadratic potential, if the pseudo Nambu-Goldstone boson accounts for the curvature perturbation; the non-Gaussianity is enhanced, otherwise.

We study a scenario that large non-Gaussianity arises from the baryon asymmetry of the Universe. There are baryogenesis scenarios containing a light scalar field, which may result in baryonic isocurvature perturbations with some amount of non-Gaussianity. As an explicit example we consider the Affleck-Dine mechanism and show that a flat direction of the supersymmeteric standard model can generate large non-Gaussianity in the curvature perturbations, satisfying the observational constraints on the baryonic isocurvature perturbations. The sign of a non-linearity parameter, f_{NL}, is negative, if the Affleck-Dine mechanism accounts for the observed baryon asymmetry; otherwise it can be either positive or negative.

We develop a formalism to study non-Gaussianity in both curvature and isocurvature perturbations. It is shown that non-Gaussianity in the isocurvature perturbation between dark matter and photons leaves distinct signatures in the CMB temperature fluctuations, which may be confirmed in future experiments, or possibly, even in the currently available observational data. As an explicit example, we consider the QCD axion and show that it can actually induce sizable non-Gaussianity for the inflationary scale, H_{inf} = O(10^9 - 10^{11})GeV.

Formation of primordial black holes (PBHs) on astrophysical mass scales is a natural consequence of inflationary cosmology if the primordial perturbation spectrum has a large and negative running of the spectral index as observationally inferred today, because double inflation is required to explain it and fluctuations on some astrophysical scales are enhanced in the field oscillation regime in between. It is argued that PBHs thus produced can serve as intermediate-mass black holes (IMBHs) which act as the observed ultraluminous X-ray sources (ULXs) by choosing appropriate values of the model parameters in their natural ranges. Our scenario can be observationally tested in near future because the mass of PBHs is uniquely determined once we specify the values of the spectral index and its running on large scales.

A model of inhomogeneous baryogenesis based on the Affleck and Dine mechanism is described. A simple coupling of the scalar baryon field to the inflaton allows for formation of astronomically significant bubbles with a large baryon (or antibaryon) asymmetry. During the farther evolution these domains form compact stellar-like objects, or lower density clouds, or primordial black holes of different size. According to the scenario, such high baryonic number objects occupy relatively small fraction of space but despite that they may significantly contribute to the cosmological mass density. For some values of parameters the model allows the possibility the whole dark matter in the universe to be baryonic. Furthermore, the model allows the existence of the antibaryonic B-bubbles, i.e. a significant fraction of the mass density in the universe can be in the form of the compact antimatter objects (e.g. anti-stars).

We consider energy deposition of high energy electrons and photons in universe. We carry out detailed calculations of fractions of the initial energy of the injected electron or photon which are used to heat, ionize and excite background plasma in the early universe for various ionization states and redshifts.

We derive big-bang nucleosynthesis (BBN) constraints on both unstable and stable gravitino taking account of recent progresses in theoretical study of the BBN processes as well as observations of primordial light-element abundances. In the case of unstable gravitino, we set the upper limit on the reheating temperature assuming that the primordial gravitinos are mainly produced by the scattering processes of thermal particles. For stable gravitino, we consider Bino, stau and sneutrino as the next-to-the-lightest supersymmetric particle and obtain constraints on their properties. Compared with the previous works, we improved the following points: (i) we use the most recent observational data, (ii) for gravitino production, we include contribution of the longitudinal component, and (iii) for the case with unstable long-lived stau, we estimate the bound-state effect of stau accurately by solving the Boltzmann equation.

We investigate the modifications to predictions for the abundances of light elements from standard Big-Bang nucleosynthesis when exotic late-decaying particles with lifetimes exceeding ~1 sec are prominent in the early Universe. Utilising a model-independent analysis of the properties of these long-lived particles, we identify the parameter space associated with models that are consistent with all observational data and hence resolve the much discussed discrepancies between observations and theoretical predictions for the abundances of Li^7 and Li^6.

We investigate cosmological constraints on primordial isocurvature and tensor perturbations, using recent observations of the cosmic microwave background and the large scale structure. We find that present observations are consistent with purely adiabatic initial conditions for the structure formation under any priors on correlations of isocurvature modes, and upper limits on the contribution of isocurvature and tensor perturbations are presented. We also apply the obtained constraints to some specific theoretical models, axion isocurvature perturbation models and curvaton models, and give some implications for theoretical models.

We present models of decaying particles for increasing the effective number of neutrinos N_nu after big bang nucleosynthesis but before the structure formation begins. We point out that our scenario not only solves the discrepancy between the constraints on N_nu from these two epochs, but also provides a possible answer to deeper inconsistency in the estimation of the matter power spectrum amplitude at small scales, represented by sigma_8, between the WMAP and some small scale matter power measurements such as the Lyman-alpha forest and weak lensing. We consider (a) saxion decay into two axions; (b) gravitino decay into axino and axion; (c) Dirac right-handed sneutrino decay into gravitino and right-handed neutrino.

We derive general constraints on the relic abundances of a long-lived particle which mainly decays into a neutrino (and something else) at cosmological time scales. Such an exotic particle may show up in various particle-physics models based on physics beyond the standard model. The constraints are obtained from big-bang nucleosynthesis, cosmic microwave background and diffuse neutrino and photon fluxes, depending on the lifetime and the electromagnetic and hadronic branching ratios.

We derive constraint on the effective number of neutrino species N_nu from the cosmic microwave background power spectrum of the WMAP and galaxy clustering power spectrum of the SDSS luminous red galaxies (LRGs). Using these two latest data sets of CMB and galaxy clustering alone, we obtain the limit 0.9 < N_nu < 8.2 (95% C.L.) for the power-law LambdaCDM flat universe, with no external prior. The lower limit corresponds to the lower bound on the reheating temperature of the universe T_R > 2 MeV.

We numerically investigate density perturbations generated in the smooth hybrid new inflation model, a kind of double inflation model that is designed to reproduce the running spectral index suggested by the WMAP results. We confirm that this model provides the running spectral index within 1sigma range of the three year WMAP result. In addition, we find a sharp and strong peak on the spectrum of primordial curvature perturbation at small scales. This originates from amplification of fluctuation in the first inflaton fields due to parametric resonance, which takes place in the oscillatory phase between two inflationary regime. Formation probability of primordial black holes (PBHs) is discussed as a consequence of such peak.

We consider constraints on long-lived charged scalar leptons \tilde{l}^\pm in supersymmetric models, where gravitino is the lightest superparticle. We study the decay and hadronization processes of \tilde{l}^\pm. We also take into account the significant enhancement of Li6 production due to the formation of the bound-state (He4 \tilde{l}^-); for this purpose, we use the reaction rate given by the most recent calculation based on coupled-channel method.

We study the reionization histories where ionizing UV photons are emitted from decaying particles, in addition to usual contributions from stars and quasars, taking account of the fact that the universe is not fully ionized until z = 6 as observed by Sloan Digital Sky Survey. Likelihood analysis of the three-year data from the WMAP (Wilkinson Microwave Anisotropy Probe) severely constrains the decaying particle scenario.In particular, the decaying particle with relatively short lifetime is not favored by the polarization data.

We point out that supergravity effects enable the inflaton to decay into all matter fields, including the visible and the supersymmetry breaking sectors, once the inflaton acquires a non-vanishing vacuum expectation value. The new decay processes have great impacts on cosmology; the reheating temperature is bounded below; the gravitinos are produced by the inflaton decay in a broad class of the dynamical supersymmetry breaking models. We derive the bounds on the inflaton mass and the vacuum expectation value, which severely constrain high-scale inflations such as the hybrid and chaotic inflation models.

A reanalysis is made for the helium abundance determination for the Izotov-Thuan (2004) spectroscopic sample of extragalactic H II regions. We find that the effect of underlying stellar absorption of the He I lines, which is more important for metal poor systems, affects significantly the inferred primordial helium abundance $Y_p$ obtained in the zero metallicity limit and the slope of linear extrapolation, $dY/dZ$. This brings $Y_p$ from $0.234\pm0.004$ to $0.250\pm 0.004$ and $dY/dZ=4.7\pm 1.0$ to $1.1\pm 1.4$. Conservatively, this indicates the importance of the proper understanding of underlying stellar absorption for accurate determinations of the primordial helium abundance to the error of $\delta Y_p\simeq 0.002-0.004$.

We show that the gravitino-overproduction problem is prevalent among inflation models in supergravity. An inflaton field \phi generically acquires (effective) non-vanishing auxiliary field G_\phi, if the K\"ahler potential is non-minimal. The inflaton field then decays into a pair of the gravitinos. We extensively study the cosmological constraints on G_\phi for a wide range of the gravitino mass. For many inflation models we explicitly estimate G_\phi, and show that the gravitino-overproduction problem severely constrains the inflation models, unless such an interaction as K = \kappa/2 |\phi|^2 z^2 + h.c. is suppressed (here z is the field responsible for the supersymmetry breaking). We find that many of them are already excluded or on the verge of, if \kappa \sim O(1).

We study the Affleck-Dine mechanism with various types of the Kahler potential, and investigate whether or not the Affleck-Dine field could acquire a large VEV as an initial condition for successful baryogenesis. In addition to a negative Hubble-induced mass term, we find that large enough Hubble-induced A-terms could also develop the minimum at large amplitude of the field. Therefore, the Affleck-Dine mechanism works for broader classes of the theories.

We investigate how the nature of dark energy affects the determination of the curvature of the universe from recent observations. For this purpose, we consider the constraints on the matter and dark energy density using observations of type Ia supernovae, baryon acoustic oscillation peak and cosmic microwave background with several types of dark energy equation of state. Although it is usually said that the combination of current observations favors a flat universe, we found that a relatively large parameter space allows the universe to be open for a particular model of dark energy. We also discuss what kind of dark energy model or prior allow a non-flat universe.

We study constraints on time variation of the fine structure constant alpha from cosmic microwave background (CMB) taking into account simultaneous change in alpha and the electron mass m_e which might be implied in unification theories. We obtain the constraints -0.097 < Delta alpha/alpha < 0.034 at 95% C.L. using WMAP data only, and -0.042 < Delta alpha/alpha < 0.026 combining with the constraint on the Hubble parameter by the HST Hubble Key Project. These are improved by 15% compared with constraints assuming only alpha varies. We discuss other relations between variations in alpha and m_e but we do not find evidence for varying alpha.

We derive an upper limit on the neutrino mass from the WMAP three-year data alone by employing a deterministic minimisation method based on a grid search in multi-parameter space. Assuming the flat LCDM model with power-law adiabatic perturbations, we find the sum of m_nu < 2.0 eV in agreement with the result of the WMAP team. This result, the limit being nearly the same as that from the WMAP first-year data, means that the fundamental limit on the neutrino mass obtainable from the cosmic microwave background alone is already nearly met, as anticipated from the previous analysis. We also clarify the role of the polarisation data in deriving the limit on the neutrino mass.

Most of the inflation models end up with non-vanishing vacuum expectation values of the inflaton fields \phi in the true vacuum, which induce, in general, nonvanishing auxiliary field G_\phi for the inflaton potential in supergravity. We show that the presence of nonzero G_\phi gives rise to inflaton decay into a pair of the gravitinos and are thereby severely constrained by cosmology especially if the gravitino is unstable and its mass is in a range of O(100) GeV \sim O(10) TeV. For several inflation models, we explicitly calculate the values of G_\phi and find that most of them are excluded or on the verge of being excluded for the gravitino mass in that range. We conclude that an inflation model with vanishing G_\phi, typically realized in a chaotic inflation, is favored in a sense that it naturally avoids the potential gravitino overproduction problem.

We obtain the spectrum of gamma ray emissions from the moduli whose decay into $e^+ e^-$ accounts for the 511 keV line observed by SPI/INTERGRAL. The moduli emit gamma rays through internal bremsstrahlung, and also decay directly into two gammas via tree and/or one-loop diagrams. We show that the internal bremsstahlung constrains the mass of the moduli below $\sim 40$ MeV model-independently. On the other hand, the flux of two gammas directly decayed from the moduli through one loop diagrams will exceed the observed galactic diffuse gamma-ray background if the moduli mass exceeds $\sim 20$ MeV in the typical situation. Moreover, forthcoming analysis of SPI data in the range of 1-8 MeV may detect the line emisson with the energy half the moduli mass in the near future, which confirms the decaying moduli scenario.

The effects of radiatively decaying, long-lived particles on big-bang nucleosynthesis (BBN) are discussed. If high-energy photons are emitted after BBN, they may change the abundances of the light elements through photodissociation processes, which may result in a significant discrepancy between the BBN theory and observation. We calculate the abundances of the light elements, including the effects of photodissociation induced by a radiatively decaying particle, but neglecting the hadronic branching ratio. Using these calculated abundances, we derive a constraint on such particles by comparing our theoretical results with observations. Taking into account the recent controversies regarding the observations of the light-element abundances, we derive constraints for various combinations of the measurements. We also discuss several models which predict such radiatively decaying particles, and we derive constraints on such models.

We present the results of a systematic study using ASCA of the ionization state for six ``mixed-morphology'' supernova emnants (MMSNRs): IC 443, W49B, W28, W44, 3C391, and Kes 27. MMSNRs show centrally filled thermal X-ray emission, which contrasts to shell-like radio morphology, a set of haracteristics at odds with the standard model of SNR evolution (e.g., the Sedov model). We have therefore studied the evolution of the MMSNRs from the ionization conditions inferred from the X-ray spectra, independent of X-ray morphology. We find highly ionized plasmas approaching ionization equilibrium in all the mmsnrs. The degree of ionization is systematically higher than the plasma usually seen in shell-like SNRs. Radial temperature gradients are also observed in five remnants, with cooler plasma toward the limb. In IC 443 and W49B, we find a plasma structure consistent with shell-like SNRs, suggesting that at least some MMSNRs have experienced similar evolution to shell-like SNRs. In addition to the results above, we have discovered an ``overionized'' ionization state in W49B, in addition to that previously found in IC 443. Thermal conduction can cause the hot interior plasma to become overionized by reducing the temperature and density gradients, leading to an interior density increase and temperature decrease. Therefore, we suggest that the ``center-filled'' X-ray morphology develops as the result of thermal conduction, and should arise in all SNRs. This is consistent with the results that MMSNRs are near collisional ionization equilibrium since the conduction timescale is roughly similar to the ionization timescale. Hence, we conclude that MMSNRs are those that have evolved over$\sim10^4$ yr. We call this phase as the ``conduction phase.''

We study how the oscillations of the neutrinos affect their thermalization process during the reheating period with temperature O(1) MeV in the early universe. We follow the evolution of the neutrino density matrices and investigate how the predictions of big bang nucleosynthesis vary with the reheating temperature. For the reheating temperature of several MeV, we find that including the oscillations makes different predictions, especially for $^4$He abundance. Also, the effects on the lower bound of the reheating temperature from cosmological observations are discussed.

We show that the $e^++e^-$ decay of a light scalar boson of mass 1-10 MeV may account for the fluxes of 511 keV gamma ray observed by SPI/INTEGRAL. We argue that candidates of such a light scalar boson is one of the string moduli or a scalar partner of the axion in a supersymmetric theory.

Q-ball is a non-topological soliton whose stability is ensured by global U(1) symmetry. We study a Q-ball which arises in the Affleck-Dine mechanism for baryogenesis and consider its possible instability due to U(1) breaking term ($A$-term) indispensable for successful baryogenesis. It is found that the instability destroys the Q-ball if its growth rate exceeds inverse of the typical relaxation time scale of the Q-ball. However, the instability is not so strong as it obstructs the cosmological formation of the Q-balls.

We study the big-bang nucleosynthesis (BBN) with the long-lived exotic particle, called X. If the lifetime of X is longer than \sim 0.1 sec, its decay may cause non-thermal nuclear reactions during or after the BBN, altering the predictions of the standard BBN scenario. We pay particular attention to its hadronic decay modes and calculate the primordial abundances of the light elements. Using the result, we derive constraints on the primordial abundance of X. Compared to the previous studies, we have improved the following points in our analysis: The JETSET 7.4 Monte Carlo event generator is used to calculate the spectrum of hadrons produced by the decay of X; The evolution of the hadronic shower is studied taking account of the details of the energy-loss processes of the nuclei in the thermal bath; We have used the most recent observational constraints on the primordial abundances of the light elements; In order to estimate the uncertainties, we have performed the Monte Carlo simulation which includes the experimental errors of the cross sections and transfered energies. We will see that the non-thermal productions of D, He3, He4 and Li6 provide stringent upper bounds on the primordial abundance of late-decaying particle, in particular when the hadronic branching ratio of X is sizable. We apply our results to the gravitino problem, and obtain upper bound on the reheating temperature after inflation.

It is shown that a subelectronvolt upper limit can be derived on the neutrino mass from the CMB data alone in the Lambda CDM model with the power-law adiabatic perturbations, without the aid of any other cosmological data. Assuming the flatness of the universe, the constraint we can derive from the current WMAP observations is \sum m_nu < 2.0 eV at the 95% confidence level for the sum over three species of neutrinos (m_nu < 0.66 eV for the degenerate neutrinos). This constraint modifies little even if we abandon the flatness assumption for the spatial curvature. We argue that it would be difficult to improve the limit much beyond \sum m_nu \lesssim 1.5 eV using only the CMB data, even if their statistics are substantially improved. However, a significant improvement of the limit is possible if an external input is introduced that constrains the Hubble constant from below. The parameter correlation and the mechanism of CMB perturbations that give rise to the limit on the neutrino mass are also elucidated.

It is shown that late-time decay of domain walls can dilute unwanted relics such as moduli, if the universe was dominated by frustrated domain walls with tension $\sigma = (1 - 100 TeV)^3$. Since energy density of the frustrated domain walls decreases as slow as the inverse of the scale factor, an overclosure limit on the axion decay constant $f_a$ is also considerably relaxed. In fact $f_a$ can be as large as the Planck scale, which may enable us to naturally implement the QCD axion in the string scheme. Furthermore, in contrast to thermal inflation models, the Affleck-Dine baryogenesis can generate enough asymmetry to explain the present baryon abundance, even in the presence of late-time entropy production.

We discuss cosmic microwave background anisotropies in models with an unconventional primordial power spectrum. In particular, we consider an initial power spectrum with some ``spiky'' corrections. Interestingly, such a primordial power spectrum generates ``wavy'' structure in the CMB angular power spectrum.

We study the reionization scenario in which ionizing UV photons emitted from decaying particle, in addition to usual contributions from stars and quasars, ionize the universe. It is found that the scenario is consistent with both the first year data of the Wilkinson Microwave Anisotropy Probe and the fact that the universe is not fully ionized until z \sim 6 as observed by Sloan Digital Sky Survey. Likelihood analysis revealed that rather broad parameter space can be chosen. This scenario will be discriminated by future observations, especially by the EE polarization power spectrum of cosmic microwave background radiation.

We study the big-bang nucleosynthesis (BBN) scenario with late-decaying exotic particles with lifetime longer than $\sim 1$ sec. With a late-decaying particle in the early universe, predictions of the standard BBN scenario can be significantly altered. Therefore, we derive constraints on its primordial abundance. We pay particular attention to hadronic decay modes of such particles. We see that the non-thermal production process of D, ${\rm ^{3}He}$ and ${\rm ^6Li}$ provides a stringent upper bound on the primordial abundance oflate-decaying particles with hadronic branching ratio.

Within the standard big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) framework, the baryon density measured by the Wilkinson Microwave Anisotropy Probe (WMAP) or the primordial D abundance is much higher than the one measured by the $^4$He or $^7$Li abundances. To solve the discrepancy, we propose a scenario in which additional baryons appear after BBN. We show that simply adding the baryons can not be a solution but the existence of a large lepton asymmetry before BBN makes the scenario successful. These extra baryons and leptons, in addition to the initial baryons which exist before the BBN, can be all produced from $Q$-balls.

We study dynamics of flat directions in the minimal supersymmetric standard model taking account of its constituent fields. It is found that there exist new instabilities due to the D-term potential and the nature of these instabilities depends on the eccentricity of the orbit. For a roughly circular orbit, it is similar to narrow-band resonance described by the Mathieu equation. For an elliptic orbit, the instabilities exhibit tachyonic nature. In the narrow-band regime, we show that the existence of the instabilities is crucial to the formation of $Q$-balls. The $Q$-ball formation proceeds through quasi-stable state called ``$I$-ball''. The transition from $I$-balls to $Q$-balls is made efficient by the $D$-term instability and during this process some fraction of the charge of the $Q$-ball is emitted. This discovery may revive the scenario which relates the baryon number and dark matter density of the universe. Furthermore, the tachyonic $D$-term instability causes the drastic decay of the flat direction moving in an orbit with relatively large eccentricity. Thus the evolution of flat directions is completely altered by the appearance of this instability.

We consider spherical collapse in the Randall-Sundrum type II model and estimate the critical over density for black hole formation in the radiation dominated era. It is found that when (density)$^{2}$-term is dominant in the modified Friedmann equation the critical density is smaller than in the standard cosmology, which implies that PBHs are more easily produced in the Randall-Sundrum model.

Flat directions in the minimal supersymmetric standard model are known to deform into non-topological solitons, Q-balls, which generally possess both baryon and lepton asymmetries. We investigate how Q-balls evolve if some of the constituent fields of the flat direction decay into light species. It is found that the Q-balls takes a new configuration whose energy per charge slightly increases due to the decay. Specifically, we show that all the stable Q-balls eventually transform into pure B-balls via the decay into neutrinos.

We calculate primordial abundances of light elements produced during big bang nucleosynthesis when the fine structure constant and/or the cosmic expansion rate take non-standard values. We compare them with the recent values of observed D, He4 and Li7 abundances, which show slight inconsistency among themselves in the standard big bang nucleosynthesis scenario. This inconsistency is not solved by considering either a varying fine structure constant or a non-standard expansion rate separately but solutions are found by their simultaneous existence.

The stability of cosmological solutions in the recently suggested specific mechanism of dynamical compensation of vacuum energy is studied. It is found that the solutions in the original version lead to cosmological singularity which could be reached in final (and short) time. A modification of the interaction of the compensating field with gravity is suggested which allows to escape such singularity. It is shown that generic cosmological solution in this model tends to the Friedmann expansion regime even starting from initially large vacuum energy.

We show that UV photons produced by decaying particles can partially reionize the universe and explain the large optical depth observed by Wilkinson Microwave Anisotropy Probe. Together with UV fluxes from early formed stars and quasars, it is possible that the universe is fully ionized at z \lesssim 6 and partially ionized at z \gtrsim 6 as observed by Sloan Digital Sky Survey for large parameter space of the decaying particle. This scenario will be discriminated by future observations, especially by the EE polarization power spectrum of cosmic microwave background radiation.

We study the curvaton scenario in supersymmetric framework paying particular attention to the fact that scalar fields are inevitably complex in supersymmetric theories. If there are more than one scalar fields associated with the curvaton mechanism, isocurvature (entropy) fluctuations between those fields in general arise, which may significantly affect the properties of the cosmic density fluctuations. We examine several candidates for the curvaton in the supersymmetric framework, such as moduli fields, Affleck-Dine field, $F$- and $D$-flat directions, and right-handed sneutrino. We estimate how the isocurvature fluctuations generated in each case affect the cosmic microwave background angular power spectrum. With the use of the recent observational result of the WMAP, stringent constraints on the models are derived and, in particular, it is seen that large fraction of the parameter space is excluded if the Affleck-Dine field plays the role of the curvaton field. Natural and well-motivated candidates of the curvaton are also listed.

We show that the recently suggested explanations of cosmic acceleration by the modification of gravity at small curvature suffer violent instabilities and strongly disagree with the known properties of gravitational interactions.

The recent observation of the cosmic microwave background anisotropy by the WMAP confirmed that the lower multipoles are considerably suppressed. From the standpoint of the cosmic variance, it is nothing but a statistical accident. Alternatively, one can attribute the deficit of fluctuation on the large scale to the cosmic history, which might be explained in the context of the inflationary physics. In this paper, we show that it is possible to explain the suppressed lower multipoles in the hybrid new inflation model.

We propose a model with a compensating scalar field whose back reaction to the cosmological curvature cancels possible vacuum energy density down to the terms of the order of the time dependent critical energy density. Thus the model simultaneously solves the mystery of the compensation of vacuum energy with the accuracy of 120 orders of magnitude and explains existence of the observed dark energy. At an early stage the suggested cosmological model might experience exponential expansion without an additional inflaton field.

We reconsider the evolution of strongly degenerate neutrinos in the early universe. Our chief concern is the validity of the entropy conservation after the neutrino annihilation process has frozen out (so that the establishment of chemical equilibrium is not trivial). We argue that the entropy indeed conserves because elastic scattering keeps the neutrino and antineutrino distribution functions in the equilibrium form and the sum of their chemical potential keeps zero even after the neutrino annihilation freeze-out. We also simulate the evolution of the degenerate neutrino spectrum to support the argument. We conclude that the change in the neutrino degeneracy parameter when the relativistic degrees of freedom in the universe decreases is calculated using the entropy conservation and the lepton number conservation without worrying about at what temperature the neutrino annihilation process freezes out.

We show that the earliest possible reionisation of the Universe is approximately at $z\simeq 13.5$ and the optical depth is $\tau\simeq 0.17$ in the conventionally accepted $\Lambda$ cold dark matter ($\Lambda$CDM) model with adiabatic fluctuations of a flat spectrum normalised to the cosmic microwave background. This is consistent with the reionisation found by WMAP (the apparently earlier reionisation epoch of the WMAP is ascribed to the adoption of the instantaneous reionisation approximation), i.e., the WMAP result is realised only if reionisation of the universe takes place nearly at the maximal efficiency in the $\Lambda$CDM model.

The first year Wilkinson Microwave Anisotropy Probe data favors primordial adiabatic fluctuation with a running spectral index with $\ns >1$ on a large scale and $\ns <1$ on a smaller scale. The model building of inflation that predicts perturbations with such a spectrum is a challenge, because most models predict fluctuations with either $\ns >1$ (hybrid inflation) or $\ns <1$ (new, chaotic, or topological inflation). We give a sensible particle physics model in supergravity that accommodates the desired running of the spectral index using double inflation.

We investigate whether the Affleck-Dine (AD) mechanism works when the contribution of the two-loop thermal correction to the potential is negative in the gauge-mediated supersymmetry breaking models. The AD field is trapped far away from the origin by the negative thermal correction for a long time until the temperature of the universe becomes low enough. The most striking feature is that the Hubble parameter becomes much smaller than the mass scale of the radial component of the AD field, during the trap. Then, the amplitude of the AD field decreases so slowly that the baryon number is not fixed even after the onset of radial oscillation. The resultant baryon asymmetry crucially depends on whether the Hubble parameter, $H$, is larger than the mass scale of the phase component of the AD field, $M_\theta$, at the beginning of oscillation. If $H < M_\theta$ holds, the formation of Q balls plays an essential role to determine the baryon number, which is found to be washed out due to the nonlinear dynamics of Q-ball formation. On the other hand, if $H > M_\theta$ holds, it is found that the dynamics of Q-ball formation does not affect the baryon asymmetry, and that it is possible to generate the right amount of the baryon asymmetry.

We study effects of the cosmic string in the D-term inflation model on the cosmic microwave background (CMB) anisotropy. In the D-term inflation model, gauged cosmic string is usually formed, which may significantly affect the CMB anisotropy. We see that its study imposes important constraint. In order to realize the minimal model of D-term inflation, we see that the coupling constant \lambda in the superpotential generating the inflaton potential should be significantly small, \lambda \lesssim O(10^{-4}-10^{-5}). We also discuss that a consistent scenario of the D-term inflation with \lambda \sim 1 can be constructed by adopting the curvaton mechanism where the cosmic density fluctuations are generated by a primordial fluctuation of a late-decaying scalar field other than the inflaton.

We discuss the implication of the observation of supernova relic neutrinos on the study of the star formation rate (SFR) in galaxies. The limit recently obtained at Super-Kamiokande (SK) is already marginally significant: The SFR we derived $\psi(t_0)<0.040M_{\odot}$yr$^{-1}$Mpc$^{-3}$(at a 90% CL) is about twice the SFR estimated from radio observations, and five times the rate from H$\alpha$ allowing for uncertainties in the model supernova neutrino flux.

We find that there exists a soliton-like solution ``I-ball'' in theories of a real scalar field if the scalar potential satisfies appropriate conditions. Although the I-ball does not have any topological or global U(1) charges, its stability is ensured by the adiabatic invariance for the oscillating field.

We re-investigate the evolution of the strongly degenerate neutrinos in the early universe. With the larger degeneracy, the neutrino number freezes at higher temperatures because the neutrino annihilation rate decreases. We consider very large degeneracy so large that the neutrino number freezes before events in which the particle degrees of freedom in the universe decrease (e.g. the muon annihilation and the quark-hadron phase transition). In such a case, the degeneracy by the time of nucleosynthesis becomes smaller than the initial degeneracy. We calculate how much it decreases from the initial value on the basis of the conservation of the neutrino number and the total entropy. We found a large drop in the degeneracy but it is not large enough to affect the current constraints on the neutrino degeneracy from BBN and CMBR.

We propose a scenario which can explain large lepton asymmetry and small baryon asymmetry simultaneously. Large lepton asymmetry is generated through Affleck-Dine (AD) mechanism and almost all the produced lepton numbers are absorbed into Q-balls (L-balls). If the lifetime of the L-balls is longer than the onset of electroweak phase transition but shorter than the epoch of big bang nucleosynthesis (BBN), the large lepton asymmetry in the L-balls is protected from sphaleron effects. On the other hand, small (negative) lepton numbers are evaporated from the L-balls due to thermal effects, which are converted into the observed small baryon asymmetry by virtue of sphaleron effects. Large and positive lepton asymmetry of electron type is often requested from BBN. In our scenario, choosing an appropriate flat direction in the minimal supersymmetric standard model (MSSM), we can produce positive lepton asymmetry of electron type but totally negative lepton asymmetry.

We propose a topological inflation model in supergravity. In this model, the vacuum expectation value (VEV) of a scalar field takes a value much larger than the gravitational scale $M_{G} \simeq 2.4 \times 10^{18}$ GeV, which is large enough to cause topological inflation. On the other hand, expansions of the K\"ahler potential and the superpotential beyond the gravitational scale are validated by the introduction of a Nambu-Goldstone-like shift symmetry. Thus, topological inflation inevitably takes place in our model.

We consider effects of the isocurvature perturbation in the framework of the tracker-type quintessence models. During the inflation, fluctuations in the amplitude of the quintessence field are generated, which provide isocurvature component of the cosmic density perturbation. Contrary to the conventional notion, it is shown that effects of the isocurvature fluctuation may become sizable in some cases, and in particular, the cosmic microwave background angular power spectrum may be significantly enhanced due to the effect of the isocurvature mode. Such an enhancement may be detectable in on-going and future experiments.

We present the results of X-ray spatial and spectral studies of the ``mixed-morphology'' supernova remnant IC 443 using ASCA. IC 443 has a center-filled image in X-ray band, contrasting with the shell-like appearance in radio and optical bands. The overall X-ray emission is thermal, not from a synchrotron nebula. ASCA observed IC 443 three times, covering the whole remnant. From the image analysis, we found that the softness-ratio map reveals a shell-like structure. At the same time, its spectra require two (1.0 keV and 0.2 keV) plasma components; the emission of the 0.2 keV plasma is stronger in the region near the shell than the center. These results can be explained by a simple model that IC 443 has a hot (1.0 keV) interior surrounded by a cool (0.2 keV) outer shell. From the emission measures, we infer that the 0.2 keV plasma is denser than the 1.0 keV plasma, suggesting pressure equilibrium between the two. In addition, we found that the ionization temperature of sulfur, obtained from H-like K$\alpha$ to He-like K$\alpha$ intensity ratio, is 1.5 keV, significantly higher than the gas temperature of 1.0 keV suggested from the continuum spectrum. The same can be concluded for silicon. Neither an additional, hotter plasma component nor a multi-temperature plasma successfully accounts for this ratio, and we conclude that the 1.0 keV plasma is overionized. This is the first time that overionized gas has been detected in a SNR. For the gas to become overionized in the absence of a photoionizing flux, it must cool faster than the ions recombine. Thermal conduction from the 1.0 keV plasma to the 0.2 keV one could cause the 1.0 keV plasma to become overionized, which is plausible within an old (3$\times10^4$ yr) SNR.

We study the Cosmic Microwave Background (CMB) anisotropies produced by cosine-type quintessence models. In our analysis, effects of the adiabatic and isocurvature fluctuations are both taken into account. For purely adiabatic fluctuations with scale invariant spectrum, we obtain a stringent constraint on the model parameters using the CMB data from COBE, BOOMERanG and MAXIMA. Furthermore, it is shown that isocurvature fluctuations have significant effects on the CMB angular power spectrum at low multipoles in some parameter space, which may be detectable in future satellite experiments. Such a signal may be used to test the cosine-type quintessence models.

We comprehensively investigate a chaotic inflation model proposed recently in the framework of supergravity. In this model, the form of K\"ahler potential is determined by a symmetry, that is, the Nambu-Goldstone-like shift symmetry, which guarantees the absence of the exponential factor in the potential for the inflaton field. Though we need the introduction of small parameters, the smallness of the parameters is justified also by symmetries. That is, the zero limit of the small parameters recovers symmetries, which is natural in the 't Hooft's sense. The leptogenesis scenario via the inflaton decay in this chaotic inflation model is also discussed. We find that the lepton asymmetry enough to explain the present baryon number density is produced for low reheating temperatures avoiding the overproduction of gravitinos.

We consider the effects on big bang nucleosynthesis (BBN) of the radiative decay of a long-lived massive particle. If high-energy photons are emitted after the BBN epoch ($t \sim 1 - 10^3$ sec), they may change the abundances of the light elements through photodissociation processes, which may result in a significant discrepancy between standard BBN and observation. Taking into account recent observational and theoretical developments in this field, we revise our previous study constraining the abundance of the radiatively-decaying particles. In particular, on the theoretical side, it was recently claimed that the non-thermal production of $^6$Li, which is caused by the photodissociation of $\hefour$, most severely constrains the abundance of the radiatively-decaying particle. We will see, however, it is premature to emphasize the importance of the non-thermal production of $^6$Li because (i) the theoretical computation of the $^6$Li abundance has large uncertainty due to the lack of the precise understanding of the $^6$Li production cross section, and (ii) the observational data of $^6$Li abundance has large errors.

We investigate a topological inflation model in supergravity. By means of numerical simulations, it is confirmed that topological inflation can take place in supergravity. We also show that the condition for successful inflation depends not only on the vacuum-expectation value (VEV) of inflaton field but also on the form of its K\"ahler potential. In fact, it is found that the required VEV of the inflaton $\phi$ can be as small as $<\phi > \simeq 1 \times M_G$, where $M_{G}$ is the gravitational scale.

We propose a chaotic inflation model in supergravity. In the model the K\"ahler potential has a Nambu-Goldstone-like shift symmetry of the inflaton chiral multiplet which ensures the flatness of the inflaton potential beyond the Planck scale. We show that the chaotic inflation naturally takes place by introducing a small breaking term of the shift symmetry in the superpotential. This may open a new branch of model building for inflationary universe in the framework of supergravity.

One of the biggest mystery of the recent observation of cosmic microwave background anisotropies by the boomerang team is insignificance of the second acoustic peak in the angular power spectrum. It is very difficult to explain such a low amplitude without assuming the higher baryon density than predicted by the Big Bang Nucleosynthesis (BBN). Employing the double inflation model in supergravity, we show that the density fluctuations produced by this inflation model can produce a sufficient low second acoustic peak with the standard value of the baryon density from BBN. It is shown that these density fluctuations are also consistent with the observations of cluster abundances and galaxy distributions.

We study a double inflation model (a hybrid inflation + a new inflation) in supergravity and discuss the formation of primordial black holes (PBHs) with mass \sim 10^{-20}-10^{5}M_{\odot}. We find that in a wide range of parameter space, we obtain PBHs which amount to \Omega \simeq 1, i.e., PBH dark matter. Also, we find a set of inflation parameters which produces PBHs evaporating now. Those PBHs may be responsible for antiproton fluxes observed by the BESS experiment.

The cosmological implication of a double inflation model with hybrid + new inflations in supergravity is studied. The hybrid inflation drives an inflaton for new inflation close to the origin through supergravity effects and new inflation naturally occurs. If the total e-fold number of new inflation is smaller than $\sim 60$, both inflations produce cosmologically relevant density fluctuations. Both cluster abundances and galaxy distributions provide strong constraints on the parameters in the double inflation model assuming $\Omega_0=1$ standard cold dark matter scenario. The future satellite experiments to measure the angular power spectrum of the cosmic microwave background will make a precise determination of the model parameters possible.

We present new observations of two anomalous X-ray pulsars (AXP) 4U 0142+61 and 1E 1048.1-5937 made in 1998 with the ASCA. The energy spectra of these two AXPs are found to consist of two components, a power-law and a blackbody emission from the neutron star surface. These observations, when compared to earlier ASCA observations in 1994 show remarkable stability in the intensity, spectral shape and pulse profile. However, we find that the spin-down rate in 1E 1048.1-5937 is not constant. In this source, we have clearly identified three epochs with spin-down rates different from each other and the average value. This has very strong implications for the magnetar hypothesis of AXPs. We also note that the spin-down rate and its variations in 1E 1048.1-5937 are much larger than what can normally be produced by an accretion disk with very low mass accretion rate corresponding to its low X-ray luminosity.

The late-time entropy production by the massive particle decay induces the various cosmological effects in the early epoch and modify the standard scenario. We investigate the thermalization process of the neutrinos after the entropy production by solving the Boltzmann equations numerically. We find that if the large entropy are produced at t $\sim$ 1 sec, the neutrinos are not thermalized very well and do not have the perfect Fermi-Dirac distribution. Then the freeze-out value of the neutron to proton ratio is altered considerably and the produced light elements, especially He4, are drastically changed. Comparing with the observational light element abundances, we find that $T_R$ < 0.7 MeV is excluded at 95 % C.L. We also study the case in which the massive particle has a decay mode into hadrons. Then we find that $T_R$ should be a little higher, i.e. $T_R$ > 2.5 MeV - 4 MeV, for the hadronic branching ratio $B_h = 10^{-2} - 1$. Possible influence of late-time entropy production on the large scale structure formation and temperature anisotropies of cosmic microwave background is studied. It is expected that the future satellite experiments (MAP and PLANCK) to measure anisotropies of cosmic microwave background radiation temperature can detect the vestige of the late-time entropy production as a modification of the effective number of the neutrino species $N_{\nu}^{\rm eff}$.

We evolve the network of global strings in the matter-dominated universe by means of numerical simulations. The existence of the scaling solution is confirmed as in the radiation-dominated universe but the scaling parameter $\xi$ takes a slightly smaller value, $\xi \simeq 0.6 \pm 0.1$, which is defined as $\xi = \rho_{s} t^{2} / \mu$ with $\rho_{s}$ the energy density of global strings and $\mu$ the string tension per unit length. The change of $\xi$ from the radiation to the matter-dominated universe is consistent with that obtained by Albrecht and Turok by use of the one-scale model. We also study the loop distribution function and find that it can be well fitted with that predicted by the one-scale model, where the number density $n_{l}(t)$ of the loop with the length $l$ is given by $n_{l}(t) = \nu/[t^2 (l + \kappa t)^2]$ with $\nu \sim 0.040$ and $\kappa \sim 0.48$. Thus, the evolution of the global string network in the matter-dominated universe can be well described by the one-scale model as in the radiation-dominated universe.

We study the initial value problem of new inflation models in the scheme of supergravity. We find that an inflaton generally acquires a mass of the order of the Hubble parameter via gravitationally suppressed interactions with particles in the thermal bath of the universe. This additional mass term stabilizes dynamically the inflaton at the local maximum of the potential, which offers the initial condition of new inflation.

The axion is the most attractive candidate to solve the strong CP problem in QCD. If it exists, the inflationary universe produces axion fluctuations which are mixtures of isocurvature and adiabatic fluctuations in general. We investigate how large isocurvature fluctuations are allowed or favored in order to explain observations of the large scale structure of the present universe. Generic flat universe models with mixed (isocurvature + adiabatic) density fluctuations are studied. It is found that the observations are consistent with the mixed fluctuation model if the ratio \alpha of the power spectrum of isocurvature fluctuations to that of adiabatic fluctuations is less than \sim 0.1. In particular, the mixed fluctuation model with \alpha \sim 0.05, total matter density \Omega_{0} =0.4, and Hubble parameter H_0=70 km/s/Mpc gives a very good fit to the observational data. Since the height of the acoustic peak in the angular power spectrum of the cosmic microwave background (CMB) radiation decreases significantly when the isocurvature fluctuations are present, the mixed fluctuation model can be tested in future satellite experiments. Ratios of the amplitude at the peak location to that at the COBE normalization scale for various models are given. Furthermore, we also obtain the amplitude of isocurvature fluctuations as a function of axion parameters and the Hubble parameter during the inflation. We discuss the axion fluctuations in some realistic inflation models and find that a significant amount of the isocurvature fluctuations are naturally produced.

Cosmological evolution of axionic strings is investigated by numerically solving field equations of a complex scalar field in 3+1 dimensions. It is shown that the global strings relax to the scaling solution with a significantly smaller number density than the case of local strings. The power spectrum of axions radiated from them is calculated from the simulation data, which is found to be highly peaked around the Hubble scale, and a more accurate constraint on the Peccei-Quinn breaking scale is obtained.

We investigate cosmological effects concerning the late-time entropy production due to the decay of non-relativistic massive particles. The thermalization process of neutrinos after the entropy production is properly solved by using the Boltzmann equation. If a large entropy production takes place at late time t$\simeq$ 1 sec, it is found that a large fraction of neutrinos cannot be thermalized. This fact loosens the tight constraint on the reheating temperature T_R from the big bang nucleosynthesis and T_R could be as low as 0.5 MeV. The influence on the large scale structure formation and cosmic microwave background anisotropies is also discussed.

We investigate cosmological moduli problem in the gauge-mediated supersymmetry breaking (GMSB). A mini-inflation (oscillating inflation) takes place when a scalar field corresponding to the flat direction in GMSB oscillates along the logarithmic potential induced by the gauge-mediation mechanism. It is shown that this oscillating inflation can sufficiently dilute the relic abundance of the string moduli for some ranges of parameters in the GMSB models.

We study a double inflation model (a preinflation + a new inflation) in supergravity and discuss the formation of primordial black holes which may be identified with massive compact halo objects (MACHOs) observed in the halo of our galaxy. The preinflation drives an inflaton for the new inflation close to the origin through supergravity effects and the new inflation naturally occurs. If the total e-fold number of the new inflation is smaller than $\sim 60$, both inflations produce cosmologically relevant density fluctuations. If the coherent inflaton oscillation after the preinflation continues until the beginning of the new inflation, density fluctuations on small cosmological scales can be set suitably large to produce black holes MACHOs of masses $\sim 1 M_{\odot}$ in a wide region of parameter space in the double inflation model.

We propose a topological inflation model in the framework of supergravity with $R$ invariance. This topological inflation model is not only free from the initial value problem of the inflaton field but also gives low reheating temperature which is favored in supergravity since the overproduction of gravitinos is avoided. Furthermore, the predicted spectrum of the density fluctuations is generally tilted, which will be tested by future observations on CMB anisotropies and large scale structure of the universe.

The M-theory, the strong-coupling heterotic string theory, presents various interesting new phenomenologies. The M-theory bulk axion is one of these. The decay constant in this context is estimated as $F_a\simeq 10^{16}$ GeV. Direct searches for the M-theory axion seem impossible because of the large decay constant. However, we point out that large isocurvature fluctuations of the M-theory axion are obtained in a hybrid inflation model, which will most likely be detectable in future satellite experiments on anisotropies of cosmic microwave background radiation.

We consider the formation of topological defects after inflation. In order to take into account the effects of the rescattering of fluctuations, we integrate the classical equation that describes the evolution of a complex scalar field on the two-dimensional lattice with a slab symmetry. The growth of fluctuations during preheating is found not to be enough for defect formation, and rather a long stage of the rescattering of fluctuations after preheating is necessary. We conclude that the topological defects are not formed if the breaking scale \eta is lager than \sim (2 - 3)\times 10^{16} GeV.

We simulate the formation and the evolution of global strings taking into account the expansion of the universe and the concomitant change of the effective potential, that is, the change from the restoration stage of the global {\it U}(1)-symmetry to the broken stage. Starting from the thermal equilibrium state, we run simulations for 100 sets of different initial conditions in 2+1 dimensions. As a result, we find deviation from scale invariance. This is because the energy loss mechanism is pair annihilation rather than intercommutation, so that a string and an anti-string attract each other under a logarithmic potential.

We examine the $X$-ray spectrum from the decay of the dark-matter moduli with mass $\sim {\cal O}(100)$keV, in particular, paying attention to the line spectrum from the moduli trapped in the halo of our galaxy. It is found that with the energy resolution of the current experiments ($\sim 10$%) the line intensity is about twice stronger than that of the continuum spectrum from the moduli that spread in the whole universe. Therefore, in the future experiments with higher energy resolutions it may be possible to detect such line photons. We also investigate the $\gamma$-ray spectrum emitted from the decay of the multi-GeV moduli. It is shown that the emitted photons may form MeV-bump in the $\gamma$-ray spectrum. We also find that if the modulus mass is of the order of 10 GeV, the emitted photons at the peak of the continuum spectrum loses their energy by the scattering and the shape of the spectrum is significantly changed, which makes the constraint weaker than that obtained in the previous works.

A generic class of string theories predicts the existence of light moduli fields, and they are expected to have masses $m_\phi$ comparable to the gravitino mass $m_{3/2}$ which is in a range of $10^{-2}$keV--1GeV in gauge-mediated supersymmetry breaking theories. Such light fields with weak interactions suppressed by the Planck scale can not avoid some stringent cosmological constraints, that is, they suffer from `cosmological moduli problems'. We show that all the gravitino mass region $10^{-2}$keV $\lesssim m_{3/2} \lesssim$ 1GeV is excluded by the constraints even if we incorporate a late-time mini-inflation (thermal inflation). However, a modification of the original thermal inflation model enables the region $10^{-2}$keV $\lesssim m_{3/2} \lesssim$ 500keV to survive the constraints. It is also stressed that the moduli can be dark matter in our universe for the mass region $10^{-2}$keV $\lesssim m_\phi \lesssim$ 100keV.

Comparing Big Bang Nucleosynthesis predictions with the light element abundances inferred from observational data, we can obtain the strong constraints on some neutrino properties, e.g. number of neutrino species, mass, lifetime. Recently the deuterium abundances were measured in high red-shift QSO absorption systems. It is expected that they are close to the primordial values, however, two groups have reported inconsistent values which are different in one order of magnitude. In this paper we show how we can constrain on $\tau$ neutrino mass and its lifetime in each case when we adopt either high or low deuterium data. We find that if $0.01 \sec \lesssim \tau_{\nutau} \lesssim 1 \sec$ and $10\mev \lesssim m_{\nutau} \lesssim 24\mev$, the theoretical predictions agree with the low D/H abundances. On the other hand if we adopt the high D/H abundances, we obtain the upper bound of $\tau$ neutrino mass, $m_{\nutau}\lesssim 20 \mev$.

We study the dynamics of a scalar field \Phi with the potential g(|\Phi|^2-\eta^2)^2/2 (g=self-coupling constant, \eta=symmetry breaking scale) after inflation and make clear whether topological defects can ever be formed during preheating. In particular, we pay attention to GUT defects (\eta \sim 10^{15}GeV - 10^{17}GeV), and consider three types of fluctuations. The first one is produced due to parametric resonance, the second is due to the negative curvature of the potential, and the last is created during inflation. We search for the parameter region that nonthermal fluctuations of the scalar field produced through the parametric resonant decay of its homogeneous part do not lead to defect formation. We find that this region is rather wide, and the GUT defects are not produced after inflation. This fact shows that the positiveness of the effective mass square of the field and production of large fluctuations whose amplitude is as large as that of homogeneous mode are not enough conditions for full symmetry restoration.

The dilaton field in string theories (if exists) is expected to have a mass of the order of the gravitino mass $m_{3/2}$ which is in a range of $10^{-2}$keV--1GeV in gauge-mediated supersymmetry breaking models. If it is the case, the cosmic energy density of coherent dilaton oscillation easily exceeds the critical density of the present universe. We show that even if this problem is solved by a late-time entropy production (thermal inflation) a stringent constraint on the energy density of the dilaton oscillation is derived from experimental upperbounds on the cosmic X($\gamma$)-ray backgrounds. This excludes an interesting mass region, $500keV \lesssim m_{3/2} \lesssim 1GeV$, in gauge-mediated supersymmetry breaking models.

We study the domain wall problem of an axion in chaotic inflation models. We show that the production of domain walls does not occurs if the Peccei-Quinn scalar has a flat potential and its breaking $F_{a}$ is larger than $\sim 10^{15}$GeV. We find that too large isocurvature fluctuations are produced for such high $F_{a}$. In order for those isocurvature fluctuations to be consistent with observations of the large scale structure of the universe, the self-coupling constant $g$ and the breaking scale of the Peccei-Quinn scalar should be $g \simeq (1-2.8)\times 10^{-13}$ and $F_{a} \simeq (0.6-1.5)\times 10^{15}$GeV, respectively. In particular the value of self-coupling constant is almost the same as that required for a chaotic inflation, which strongly suggests that the Peccei-Quinn scalar itself is an inflaton.

Recently it is reported that there is the discrepancy between big bang nucleosynthesis theory and observations (BBN crisis). We show that BBN predictions agree with the primordial abundances of light elements, He4, D, He3 and Li7 inferred from the observational data if an electron neutrino has a net chemical potential xi_{nu_e} due to lepton asymmetry. We estimate that xi_{nu_e} = 0.043^{+0.040}_{-0.040} (95% C.L.) and Omega_bh^2 = 0.015^{+0.006}_{-0.003} (95% C.L.).

Banks and Dine have recently shown that the M theory naturally accommodates the Peccei-Quinn axion. Since the decay constant $F_a$ of the axion is large as $F_a \simeq 10^{15}-10^{16}$GeV, the halo axion is hardly detected in coming axion-search experiments. However, we show that isocurvature fluctuations of the M-theory axion produced at the inflationary epoch are most likely detectable in future satellite experiments on anisotropies of the cosmic microwave background radiation.

We discuss a new mechanism which can solve the crisis in standard big-bang nucleosynthesis. A long-lived particle $X$ ($10^4 sec \lesssim \tau_X \lesssim 10^6 sec$) which decays into photon(s) will induce cascade photons, and destroy significant amounts of D and $^3$He without destroying $^4$He or too much $^7$Li. We numerically investigate this process and derive a constraint on the properties of $X$ such that the theoretical values of the primordial light-element abundances agree with observation. We also present some candidates for the unstable particle $X$.

We argue that an extension of the Minimal Supersymmetric Standard Model that gives rise to viable thermal inflation, and so does not suffer from a Polonyi/moduli problem, should contain right-handed neutrinos which acquire their masses due to the vacuum expectation value of the flaton that drives thermal inflation. This strongly disfavours SO(10) Grand Unified Theories. The $\mu$-term of the MSSM should also arise due to the vev of the flaton. With the extra assumption that $ m_L^2 - m_{H_u}^2 < 0 $, but of course $ m_L^2 - m_{H_u}^2 + |\mu|^2 > 0 $, we show that a complicated Affleck-Dine type of baryogenesis employing an $LH_u$ $D$-flat direction can naturally generate the baryon asymmetry of the Universe.

The $ee\gamma\gamma +\rlap/E_T$ event observed by the CDF at Fermilab is naturally explained by dynamically supersymmetry breaking models and suggests the presence of the light gravitino which can be a warm dark matter. We consider large scale structure of the universe in the worm dark matter model and find that the warm dark matter plays almost the same role in the formation of the large scale structure as a cold dark matter if its mass is about $0.5$keV. We also study the Ly~$\alpha$ absorption systems which are presumed to be galaxies at high redshifts and show that the baryon density in the damped Ly~$\alpha$ absorption systems predicted by the warm dark matter model is quite consistent with the present observation.

In the chaotic inflation models, quantum fluctuations for axion fields lead to the overproduction of domain walls and too large isocurvature fluctuations which is inconsistent with the observations of cosmic microwave background anisotropies. These problems are solved by assuming a very flat potential for the Peccei-Quinn scalar. As the simplest possibility, we consider a model where the Peccei-Quinn scalar is an inflaton itself and show that the isocurvature fluctuations can be comparable with the adiabatic ones. We investigate cosmological implications in the case that both adiabatic and isocurvature fluctuations exist and find that the amplitude of the matter spectrum becomes smaller than that for the pure adiabatic case. This leads to relatively high bias parameter ($b \simeq 2$) which is favoured by the current observations.

We have reexamine the effect of entropy production on the cosmic axion density and find that the Peccei-Quinn scale $F_a$ larger than about $10^{15}$~GeV is not allowed even if large entropy is produced by the decays of coherent oscillations or non-relativistic massive particles. We stress that this result is independent of the details of models for the decaying particles.

We study the Polonyi problem in the framework of no-scale type supergravity models. We show that the lightest superparticle (LSP) produced in the decay of the Polonyi field may contribute too much to the present density of the universe. By requiring that LSP should not overclose the universe, we obtain a stringent constraint on the reheating temperature after the decay of the Polonyi field. We calculate the LSP density with physical parameters obtained by solving renormalization group equations in the minimal supersymmetric SU(5) model and find that the reheating temperature should be greater than about 100MeV which corresponds to $O$(100)TeV of the Polonyi mass.

We investigate the electromagnetic cascade initiated by injection of very high energy photons in the early Universe and calculate the cascade spectrum by solving a set of Boltzmann equations numerically. In the calculation we take account of Compton scattering off the background electrons and pair creation off the background nucleons as well as photon-photon processes and inverse Compton scattering. We also apply our cascade spectrum to the big bang nucleosynthesis with photo-dissociation processes due to heavy unstable particles and obtain the constraint on their lifetime and abundance.

Gravitino produced in the inflationary universe are studied. When the gravitino decays into a neutrino and a sneutrino, the emitted high energy neutrinos scatter off the background neutrinos and produce charged leptons (mainly electrons and positrons), which cause the electro-magnetic cascades and produce many soft photons. We obtain the spectra of the high energy neutrinos as well as the spectrum of the high energy photon by integrating a set of Boltzmann equations. Requiring these photons should not alter the abundances of the light elements (D, $^3$He, $^4$He) in the universe, we can set the stringent upperbound on the reheating temperature after the inflation. We find that $T_R \lesssim (10^{10}-10^{12})$GeV for $m_{3/2}\sim (100\GEV - 1\TEV)$, which is more stringent than the constraints in the previous works.

We examine the collapse of an axion domain wall bounded by an axionic string. It is found that the collapse proceeds quickly and axion domain walls disappear. However axions are emitted in the collapse and its energy density increases during radiation dominated era and contributes significantly to the present mass density of the universe. In particular the axion emitted from the wall can account for the dark matter in the universe for $F_a\gsim 10^{10}\GeV$.

Gravitino production and decay in the inflationary universe are reexanimed. Assuming that gravitino mainly decays into photon and photino, we have calculated the upperbound of the reheating temperature. Compared to previous works, we have essentially improved the following two points; (i) the helicity $\pm\frac{3}{2}$ gravitino production cross sections are calculated by using the full relevant terms in the supergravity lagrangian, and (ii) high energy photon spectrum is obtained by solving the Boltzmann equations numerically. Photo-dissociation of the light elements (D, T, $^3$He, $^4$He) leads to the most stringent upperbound of the reheating temperature, which is given by ($10^{6}$--$10^{9}$)GeV for the gravitino mass 100GeV--1TeV. On the other hand, requiring that the present mass density of photino should be smaller than the critical density, we find that the reheating temperature have to be smaller than ($10^{11}$--$10^{12}$)GeV for the photino mass (10--100)GeV, irrespectively of the gravitino mass. The effect of other decay channel is also considered.

We investigate reheating of the universe by early formation of stars and quasars in the hierarchical clustering scheme of cold dark matter scenario, with perturbation fluctuations normalized by the COBE data. It is found that ionizing uv flux from OB stars with the abundance given by the standard initial mass function is strong enough to ionize the universe from $z\approx 30$ to the present epoch, if 1--2\% of the collapsed baryons go into stars. This lessens significantly the CMB anisotropies at a small angular scale. Reionization also increases the Jeans mass to $M_{\rm luminous}\approx 10^9M_\odot$ for $z\lsim 10$, which leads to a cut-off of the luminosity function of normal galaxies on a faint side. A strong uv flux is expected at $z\approx 2-5$, and the null result of the Gunn-Peterson test is naturally explained. Early star formation also results in production of heavy elements, and the observational metal abundance sets a strong constraint on the photon energy injection into the intergalactic space.