LIGO has detected gravitational waves from massive binary black hole mergers. In order to explain the origin of such massive stellar-mass black holes, extreme metal poor stars including first stars have been invoked. However, black holes do not carry information of the metallicity. In order to check the metallicity dependence of the black hole formation, we focus on galactic black hole-main sequence binaries (BH-MSs). Using a binary population synthesis method, we find that $\gaia$ can detect $\sim200-400$ BH-MSs whose metallicity is $\zsun$ and $\sim70-400$ BH-MSs whose metallicity is $0.1\zsun$. With the spectroscopic observation on 4-m class telescopes, we can check the metallicity of BH-MSs. The metallicity dependence of the black hole formation might be checked by the astrometry and spectroscopic observations.

We study the frame dependence/independence of cosmological observables under disformal transformations, extending the previous results regarding conformal transformations, and provide the correspondence between Jordan-frame and Einstein-frame variables. We consider quantities such as the gravitational constant in the Newtonian limit, redshift, luminosity and angular diameter distances, as well as the distance-duality relation. Also, the Boltzmann equation, the observed specific intensity, and the adiabaticity condition are discussed. Since the electromagnetic action changes under disformal transformations, photons in the Einstein frame no longer propagate along null geodesics. As a result, several quantities of cosmological interest are modified. Nevertheless, we show that the redshift is invariant and the distance-duality relation (the relation between the luminosity distance and the angular diameter distance) still holds in general spacetimes even though the reciprocity relation (the relation between two geometrical distances) is modified.

We study the distortions of equilibrium spectra of relic neutrinos due to the interactions with electrons, positrons, and neutrinos in the early Universe. We solve the integro-differential kinetic equations for the neutrino density matrix, including three-flavor oscillations and finite temperature corrections from QED up to the next-to-leading order $\mathcal{O}(e^3)$ for the first time. In addition, the equivalent kinetic equations in the mass basis of neutrinos are directly solved, and we numerically evaluate the distortions of the neutrino spectra in the mass basis as well, which can be easily extrapolated into those for non-relativistic neutrinos in the current Universe. In both bases, we find the same value of the effective number of neutrinos, $N_{\rm eff} = 3.044$, which parameterizes the total neutrino energy density.

With an emphasis on improving the fidelity even in super-resolution regimes, new imaging techniques have been intensively developed over the last several years, which may provide substantial improvements to the interferometric observation of protoplanetary disks. In this study, sparse modeling (SpM) is applied for the first time to observational data sets taken by the Atacama Large Millimeter/submillimeter Array (ALMA). The two data sets used in this study were taken independently using different array configurations at Band 7 (330 GHz), targeting the protoplanetary disk around HD 142527; one in the shorter-baseline array configuration (~ 430 m), and the other in the longer-baseline array configuration (~ 1570 m). The image resolutions reconstructed from the two data sets are different by a factor of ~ 3. We confirm that the previously known disk structures appear on the images produced by both SpM and CLEAN at the standard beam size. The image reconstructed from the shorter-baseline data using the SpM matches that obtained with the longer-baseline data using CLEAN, achieving a super-resolution image from which a structure finer than the beam size can be reproduced. Our results demonstrate that on-going intensive development in the SpM imaging technique is beneficial to imaging with ALMA.

We propose a simple but novel cosmological scenario where both the Planck mass and the dark energy scale emerge from the same super-Hubble quantum fluctuations of a non-minimally coupled ultra-light scalar field during primordial inflation. The current cosmic and solar-system observations constrain the non-minimal coupling to be small.

We develop an effective-field-theory (EFT) framework for inflation with various symmetry breaking pattern. As a prototype, we formulate anisotropic inflation from the perspective of EFT and construct an effective action of the Nambu-Goldstone bosons for the broken time translation and rotation symmetries. We also calculate the statistical anisotropy in the scalar two-point correlation function for concise examples of the effective action.

In this paper, we fully investigate cosmological scenarios in which dark matter is non-minimally coupled to an extra scalar degree of freedom. The interaction is realized by means of conformal and disformal terms in the transformed gravitational metric. Considering linear perturbation theory, we show that the growth rate of dark matter differs from the uncoupled case and that the well-known Kaiser formula undergoes modification. As a result, redshift space distortion measurements cease to be a direct probe of the linear growth rate of total matter, since the distortion factor has an extra, coupling-dependent term. We study the effect of the coupling in three cosmological models, two conformally and one disformally coupled, and forecast the constraints on the coupling, and other cosmological parameters, from future galaxy surveys.

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary system containing the pulsar PSR B1259-63 orbiting around a Be star LS 2883 after the 2017 periastron passage. We detected radio continuum emission from the binary system in the millimeter/submillimeter wavelengths for the first time. At Band 3 (97 GHz), the flux 84 days after the periastron is almost the same as that 71 days after the periastron. Although the binary system showed intense GeV gamma-ray flares during our observations, the Band 3 flux did not indicate any time correlation with them. The Band 3 fluxes are consistent with an extrapolation of the radio spectrum at lower frequencies. Assuming that it is synchrotron emission, we constrain magnetic fields ($\lesssim 0.6$ G) and the high-energy cutoff of the electrons ($\gamma \gtrsim 360$). The flux at Band 7 (343 GHz) 69 days after the periastron shows a significant excess from the extrapolation of the radio spectrum at lower frequencies. The flux may be associated with the circumstellar disk around the Be star. We also present the results of Australian Telescope Compact Array (ATCA) observations at 94 GHz for the 2014 periastron passage, which show that the radio spectrum was relatively soft when the pulsar passed the disk.

We report our observations of HSC16aayt (SN 2016jiu), which was discovered by the Subaru/Hyper Suprime-Cam (HSC) transient survey conducted as part of Subaru Strategic Program (SSP). It shows very slow photometric evolution and its rise time is more than 100 days. The optical magnitude change in 400 days remains within 0.6 mag. Spectra of HSC16aayt show a strong narrow emission line and we classify it as a Type IIn supernova. The redshift of HSC16aayt is 0.6814 +/- 0.0002 from the spectra. Its host galaxy center is at 5 kpc from the supernova location and HSC16aayt might be another example of isolated Type IIn supernovae, although the possible existence of underlying star forming activity of the host galaxy at the supernova location is not excluded.

The presence of a light vector field coupled to a scalar field during inflation makes a distinct prediction: the observed correlation functions of the cosmic microwave background (CMB) become statistically anisotropic. We study the implications of the current bound on statistical anisotropy derived from the Planck 2013 CMB temperature data for such a model. The previous calculations based on the attractor solution indicate that the magnitude of anisotropy in the power spectrum is proportional to $N^2$, where $N$ is the number of $e$-folds of inflation counted from the end of inflation. In this paper, we show that the attractor solution is not compatible with the current bound, and derive new predictions using another branch of anisotropic inflation. In addition, we improve upon the calculation of the mode function of perturbations by including the leading-order slow-roll corrections. We find that the anisotropy is roughly proportional to $[2(\varepsilon_H+4\eta_H)/3-4(c-1)]^{-2}$, where $\varepsilon_H$ and $\eta_H$ are the usual slow-roll parameters and $c$ is the parameter in the model, regardless of the form of potential of an inflaton field. The bound from Planck implies that breaking of rotational invariance during inflation (characterized by the background homogeneous shear divided by the Hubble rate) is limited to be less than ${\cal O}(10^{-9})$. This bound is many orders of magnitude smaller than the amplitude of breaking of time translation invariance, which is observed to be ${\cal O}(10^{-2})$.

We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.

We discuss vacuum decays catalyzed by spherical and horizonless objects and show that an ultra compact object could catalyze a vacuum decay around it within the cosmological time. The catalytic effect of a horizonless compact object could be more efficient than that of a black hole since in this case there is no suppression of the decay rate due to the decrement of its Bekestein entropy. If there exists another minimum with AdS vacuum in the Higgs potential at a high energy scale, the abundance of compact objects such as monopoles, neutron stars, axion stars, oscillons, Q-balls, black hole remnants, gravastars and so on, could be severely constrained. We find that an efficient enhancement of nucleation rate occurs when the size of the compact object is comparable to its Schwarzschild radius and the bubble radius.

We propose a new class of higher derivative scalar-tensor theories without the Ostrogradsky's ghost instabilities. The construction of our theory is originally motivated by a scalar field with spacelike gradient, which enables us to fix a gauge in which the scalar field appears to be non-dynamical. We dub such a gauge as the spatial gauge. Though the scalar field loses its dynamics, the spatial gauge fixing breaks the time diffeomorphism invariance and thus excites a scalar mode in the gravity sector. We generalize this idea and construct a general class of scalar-tensor theories through a non-dynamical scalar field, which preserves only spatial covariance. We perform a Hamiltonian analysis and confirm that there are at most three (two tensors and one scalar) dynamical degrees of freedom, which ensures the absence of a degree of freedom due to higher derivatives. Our construction opens a new branch of scalar-tensor theories with higher derivatives.

We present Keck spectroscopic observations of three probable high redshift superluminous supernovae (SLSNe) from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), confirming redshifts of 1.851, 1.965 and 2.399. The host galaxies were selected for transient monitoring from multi-band photometric redshifts. The supernovae are detected during their rise, and the classically scheduled spectra are collected near maximum light. The restframe far-ultraviolet (FUV;$\sim$1000A-2500A) spectra include a significant host galaxy flux contribution and we compare our host galaxy subtracted spectra to UV-luminous SNe from the literature. While the signal-to-noise ratios of the spectra presented here are sufficient for redshift confirmation, supernova spectroscopic type confirmation remains inconclusive. The success of the first SHIZUCA Keck spectroscopic follow-up program demonstrates that campaigns such as SHIZUCA are capable of identifying high redshift SLSNe with sufficient accuracy, speed and depth for rapid, well-cadenced and informative follow-up.

We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper (Curtin et al. 2019). The supernovae with the confirmed redshifts of z ~ 2 have the rest ultraviolet peak magnitudes close to -21 mag and they are likely superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to the approximate event rate of ~ 900 +/- 520 Gpc-3 yr-1 with Poisson error, which is consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate the approximate superluminous supernova rates of ~ 400 +/- 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 +/- 500 Gpc-3 yr-1 (z ~ 4) with Poisson errors. Our initial results demonstrate the outstanding capability of Hyper Suprime-Cam to discover high-redshift supernovae.

We study twin Higgs models at non-zero temperature and discuss cosmological phase transitions as well as their implications on electroweak baryogenesis and gravitational waves. It is shown that the expectation value of the Higgs field at the critical temperature of the electroweak phase transition is much smaller than the critical temperature, which indicates two important facts: (i) the electroweak phase transition cannot be analyzed perturbatively (ii) the electroweak baryogenesis is hardly realized in the typical realizations of twin Higgs models. We also analyze the phase transition associated with the global symmetry breaking, through which the Standard Model Higgs is identified with one of the pseudo-Nambu-Goldstone bosons in terms of its linear realization, with and without supersymmetry. For this phase transition, we show that, only in the supersymmetric case, there are still some parameter spaces, in which the perturbative approach is validated and the phase transition is the first order. We find that the stochastic gravitational wave background is generated through this first order phase transition, but it is impossible to be detected by DECIGO or BBO in the linear realization and the decoupling limit. The detection of stochastic gravitational wave background with the feature of first order phase transition, therefore, will give strong constraints on twin Higgs models.

Dust grains of crystalline silicate, which are rarely presented in interstellar space, were found in cometary nuclei. These crystalline silicates are thought to have formed by annealing of amorphous silicate grains or direct condensation of gaseous materials near the Sun in the solar nebula, and incorporated into cometary nuclei in the cold comet-forming region after radial transportation of grains in the solar nebula. Abundances of the crystalline silicate dust grains were therefore expected to be smaller farther from the Sun. We aim to better understand the formation mechanism of minerals incorporated into comet 17P/Holmes based on its mineral abundances. To derive the mineral composition of comet 17P/Holmes, we applied a thermal emission model for cometary dust grains to mid-infrared spectra of comet 17P/Holmes taken with the Cooled Mid-Infrared Camera and Spectrograph (COMICS) mounted on the Subaru Telescope a few days later the great outburst in October 2007. The resulting mass fraction of crystalline silicate, f_ cry, and an olivine-to-pyroxene abundance ratio, f_OP, are f_cry = 0.31 +/- 0.03 and f_OP = 1.20 +0.16/-0.12, respectively. Based on a simple consideration of the mixing of dust grains originating in both the interstellar medium and solar nebula, the minerals of 17P/Holmes formed by non-equilibrium condensation. This result is consistent with theoretical and experimental predictions for vaporization and condensation of olivine in the solar nebula.

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.

We study the prospect of the Gaia satellite to identify black hole binary systems by detecting the orbital motion of the companion stars. Taking into account the initial mass function, mass transfer, common envelope phase, interstellar absorption and identifiability of black holes, we estimate the number of black hole binaries detected by Gaia and their distributions with respect to the black hole mass for several models with different parameters. We find that $\sim 200-1000$ black hole binaries will be detected by Gaia during its 5 years operation. The shape of distribution function of the black hole mass is affected most severely by the relation between the zero-age main sequence (ZAMS) stellar mass and the black hole mass in parameters adopted in this paper, which implies that black holes detected with Gaia enables us to constrain the mass relation.

As the first step to extend our understanding of higher-derivative theories, within the framework of analytic mechanics of point particles, we construct a ghost-free theory involving third-order time derivatives in Lagrangian. While eliminating linear momentum terms in the Hamiltonian is necessary and sufficient to kill the ghosts associated with higher derivatives for Lagrangian with at most second-order derivatives, we find that this is necessary but not sufficient for the Lagrangian with higher than second-order derivatives. We clarify a set of ghost-free conditions under which we show that the Hamiltonian is bounded, and that equations of motion are reducible into a second-order system.

The superhorizon (iso)curvature perturbations are conserved if the following conditions are satisfied: (i) (each) non adiabatic pressure perturbation is zero, (ii) the gradient terms are ignored, that is, at the leading order of the gradient expansion (iii) (each) total energy momentum tensor is conserved. We consider the case with the violation of the last two requirements and discuss the generation of secondary isocurvature perturbations during the late time universe. Second order gradient terms are not necessarily ignored even if we are interested in the long wavelength modes because of the convolutions which may pick products of short wavelength perturbations up. We then introduce second order conserved quantities on superhorizon scales under the conditions (i) and (iii) even in the presence of the gradient terms by employing the full second order cosmological perturbation theory. We also discuss the violation of the condition (iii), that is, the energy momentum tensor is conserved for the total system but not for each component fluid. As an example, we explicitly evaluate second order heat conduction between baryons and photons due to the weak Compton scattering, which dominates the period just before recombination. We show that such secondary effects can be recast into the isocurvature perturbations on superhorizon scales if the local type primordial non Gaussianity exists a priori.

We investigate inflation models in Jordan frame supergravity, in which an inflaton non-minimally couples to the scalar curvature. By imposing the condition that an inflaton would have the canonical kinetic term in the Jordan frame, we construct inflation models with asymptotically flat potential through pole inflation technique and discuss their relation to the models based on Einstein frame supergravity. We also show that the model proposed by Ferrara et al. has special position and the relation between the K\"ahler potential and the frame function is uniquely determined by requiring that scalars take the canonical kinetic terms in the Jordan frame and that a frame function consists only of a holomorphic term (and its anti-holomorphic counterpart) for symmetry breaking terms. Our case corresponds to relaxing the latter condition.

We perform a $z$-band survey for an optical counterpart of a binary neutron star coalescence GW170817 with Subaru/Hyper Suprime-Cam. Our untargeted transient search covers $23.6$ deg$^2$ corresponding to the $56.6\%$ credible region of GW170817 and reaches the $50\%$ completeness magnitude of $20.6$ mag on average. As a result, we find 60 candidates of extragalactic transients, including J-GEM17btc (a.k.a. SSS17a/DLT17ck). While J-GEM17btc is associated with NGC 4993 that is firmly located inside the 3D skymap of GW170817, the other 59 candidates do not have distance information in the GLADE v2 catalog or NASA/IPAC Extragalactic Database (NED). Among 59 candidates, 58 are located at the center of extended objects in the Pan-STARRS1 catalog, while one candidate has an offset. We present location, $z$-band apparent magnitude, and time variability of the candidates and evaluate the probabilities that they are located inside of the 3D skymap of GW170817. The probability for J-GEM17btc is $64\%$ being much higher than those for the other 59 candidates ($9.3\times10^{-3}-2.1\times10^{-1}\%$). Furthermore, the possibility, that at least one of the other 59 candidates is located within the 3D skymap, is only $3.2\%$. Therefore, we conclude that J-GEM17btc is the most-likely and distinguished candidate as the optical counterpart of GW170817.

We develop a method for identifying a compact object in binary systems with astrometric measurements and apply it to some binaries. Compact objects in some high-mass X-ray binaries and gamma-ray binaries are unknown, which is responsible for the fact that emission mechanisms in such systems have not yet confirmed. The accurate estimate of the mass of the compact object allows us to identify the compact object in such systems. Astrometric measurements are expected to enable us to estimate the masses of the compact objects in the binary systems via a determination of a binary orbit. We aim to evaluate the possibility of the identification of the compact objects for some binary systems. We then calculate probabilities that the compact object is correctly identified with astrometric observation (= confidence level) by taking into account a dependence of the orbital shape on orbital parameters and distributions of masses of white dwarfs, neutron stars, and black holes. We find that the astrometric measurements with the precision of 70 $\mu$as for $\gamma$~Cas allows us to identify the compact object at 99% confidence level if the compact object is a white dwarf with 0.6 M$_{\odot}$. In addition, we can identify the compact object with the precision of 10 $\mu$as at 97% or larger confidence level for LS~I~+61$^{\circ}$~303 and 99% or larger for HESS~J0632+057. These results imply that the astrometric measurements with the 10-$\mu$as precision level can realize the identification of compact objects for $\gamma$~Cas, LS~I~+61$^{\circ}$~303, and HESS~J0632+057.

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source.

We present an impact of coupling between dark matter and a scalar field, which might be responsible for dark energy, on measurements of redshift-space distortions. We point out that, in the presence of conformal and/or disformal coupling, linearized continuity and Euler equations for total matter fluid significantly deviate from the standard ones even in the sub-horizon scales. In such a case, a peculiar velocity of total matter field is determined not only by a logarithmic time derivative of its density perturbation but also by density perturbations for both dark matter and baryon, leading to a large modification of the physical interpretation of observed data obtained by measurements of redshift-space distortions. We reformulate galaxy two-point correlation function in the redshift space based on the modified continuity and Euler equations. We conclude from the resultant formula that the true value of the linear growth rate of large-scale structure cannot be necessarily constrained by single-redshift measurements of the redshift-space distortions, unless one observes the actual time-evolution of structure.

We analyse the performance of a gas time projection chamber (TPC) as a high-performance gamma-ray telescope and polarimeter in the e$^+$e$^-$ pair creation regime. We use data collected at a gamma-ray beam of known polarisation. The TPC provides two orthogonal projections $(x,z)$ and $(y,z)$ of the tracks induced by each conversion in the gas volume. We use a simple vertex finder in which vertices and pseudo-tracks exiting from them are identified. We study the various contributions to the single-photon angular resolution using Monte Carlo simulations and compare them with the experimental data and find that they are in excellent agreement. The distribution of the azimutal angle of pair conversions shows a bias due to the non-cylindrical-symmetric structure of the detector. This bias would average out for a long duration exposure on a space mission, but for this pencil-beam characterisation we have ensured its accurate simulation by a double systematics control scheme, data taking with the detector rotated at several angles with respect to the beam polarisation direction and systematics control with a non-polarised beam. We measure, for the first time, the polarisation asymmetry of a linearly polarised gamma-ray beam in the low energy pair creation regime. This sub-GeV energy range is critical for cosmic sources as their spectra are power laws which fall quickly as a function of increasing energy. This work could pave the way to extending polarised gamma-ray astronomy beyond the MeV energy regime.

We study co-existence system of both bosonic and fermionic degrees of freedom. For such system with up to first derivatives in Lagrangian, we find Ostrogradsky-type ghost-free condition in Hamiltonian analysis, which is found to be the same with requiring that the equations of motion of fermions are first-order in Lagrangian formulation. When fermionic degrees of freedom are present, uniqueness of time evolution is not guaranteed a priori because of the Grassmann property. We confirm that the additional condition, which is introduced to close Hamiltonian analysis, also ensures the uniqueness of the time evolution of system.

We revisit nonsingular cosmologies in which the limiting curvature hypothesis is realized. We study the cosmological perturbations of the theory and determine the general criteria for stability. For the simplest model, we find generic Ostrogradski instabilities unless the action contains the Weyl tensor squared with the appropriate coefficient. When considering two specific nonsingular cosmological scenarios(one inflationary and one genesis model), we find ghost and gradient instabilities throughout most of the cosmic evolution. Furthermore, we show that the theory is equivalent to a theory of gravity where the action is a general function of the Ricci and Gauss-Bonnet scalars, and this type of theory is known to suffer from instabilities in anisotropic backgrounds. This leads us to construct a new type of curvature invariant scalar function. We show that it does not have Ostrogradski instabilities, and it avoids ghost and gradient instabilities for most of the interesting background inflationary and genesis trajectories. We further show that it does not possess additional new degrees of freedom in an anisotropic spacetime. This opens the door for studying stable alternative nonsingular very early universe cosmologies.

We derive a direct correlation between the power spectrum and bispectrum of the primordial curvature perturbation in terms of the Goldstone mode based on the effective field theory approach to inflation. We show examples of correlated bispectra for the parametrized feature models presented by the Planck collaboration. We also discuss the consistency relation and the validity of our explicit correlation between the power spectrum and bispectrum.

In the context of classical mechanics, we study the conditions under which higher-order derivative theories can evade the so-called Ostrogradsky instability. More precisely, we consider general Lagrangians with second order time derivatives, of the form $L(\ddot\phi^a,\dot\phi^a,\phi^a;\dot q^i,q^i)$ with $a = 1,\cdots, n$ and $i = 1,\cdots, m$. For $n=1$, assuming that the $q^i$'s form a nondegenerate subsystem, we confirm that the degeneracy of the kinetic matrix eliminates the Ostrogradsky instability. The degeneracy implies, in the Hamiltonian formulation of the theory, the existence of a primary constraint, which generates a secondary constraint, thus eliminating the Ostrogradsky ghost. For $n>1$, we show that, in addition to the degeneracy of the kinetic matrix, one needs to impose extra conditions to ensure the presence of a sufficient number of secondary constraints that can eliminate all the Ostrogradsky ghosts. When these conditions that ensure the disappearance of the Ostrogradsky instability are satisfied, we show that the Euler-Lagrange equations, which involve a priori higher order derivatives, can be reduced to a second order system.

Current $\gamma$-ray telescopes suffer from a gap in sensitivity in the energy range between 100keV and 100MeV, and no polarisation measurement has ever been done on cosmic sources above 1MeV. Past and present e$^+$e$^-$ pair telescopes are limited at lower energies by the multiple scattering of electrons in passive tungsten converter plates. This results in low angular resolution, and, consequently, a drop in sensitivity to point sources below 1GeV. The polarisation information, which is carried by the azimuthal angle of the conversion plane, is lost for the same reasons. HARPO (Hermetic ARgon POlarimeter) is an R\&D program to characterise the operation of a gaseous detector (a Time Projection Chamber or TPC) as a high angular-resolution and sensitivity telescope and polarimeter for $\gamma$ rays from cosmic sources. It represents a first step towards a future space instrument in the MeV-GeV range. We built and characterised a 30cm cubic demonstrator [SPIE 91441M], and put it in a polarised $\gamma$-ray beam at the NewSUBARU accelerator in Japan. Data were taken at photon energies from 1.74MeV to 74MeV, and with different polarisation configurations. We describe the experimental setup in beam. We then describe the software we developed to reconstruct the photon conversion events, with special focus on low energies. We also describe the thorough simulation of the detector used to compare results. Finally we will present the performance of the detector as extracted from this analysis and preliminary measurements of the polarisation asymmetry. This beam-test qualification of a gas TPC prototype in a $\gamma$-ray beam could open the way to high-performance $\gamma$-ray astronomy and polarimetry in the MeV-GeV energy range in the near future.

Minimal chaotic models of D-term inflation predicts too large primordial tensor perturbations. Although it can be made consistent with observations utilizing higher order terms in the K\"ahler potential, expansion is not controlled in the absence of symmetries. We comprehensively study the conditions of K\"ahler potential for D-term plateau-type potentials and discuss its symmetry. They include the alpha-attractor model with a massive vector supermultiplet and its generalization leading to pole inflation of arbitrary order. We extend the models so that it can describe Coulomb phase, gauge anomaly is cancelled, and fields other than inflaton are stabilized during inflation. We also point out a generic issue for large-field D-term inflation that the masses of the non-inflaton fields tend to exceed the Planck scale.

Current {\gamma}-ray telescopes based on photon conversions to electron-positron pairs, such as Fermi, use tungsten converters. They suffer of limited angular resolution at low energies, and their sensitivity drops below 1 GeV. The low multiple scattering in a gaseous detector gives access to higher angular resolution in the MeV-GeV range, and to the linear polarisation of the photons through the azimuthal angle of the electron-positron pair. HARPO is an R&D program to characterise the operation of a TPC (Time Projection Chamber) as a high angular-resolution and sensitivity telescope and polarimeter for {\gamma} rays from cosmic sources. It represents a first step towards a future space instrument. A 30 cm cubic TPC demonstrator was built, and filled with 2 bar argon-based gas. It was put in a polarised {\gamma}-ray beam at the NewSUBARU accelerator in Japan in November 2014. Data were taken at different photon energies from 1.7 MeV to 74 MeV, and with different polarisation configurations. The electronics setup is described, with an emphasis on the trigger system. The event reconstruction algorithm is quickly described, and preliminary measurements of the polarisation of 11 MeVphotons are shown.

The covariant St\"{u}ckelberg formalism is applied to bi-gravity in order to revisit the issue of absence of the Boulware-Deser (BD) ghost. We first confirm that the leading order action in the decoupling limit for helicity-2 modes of metrics and helicity-0 mode of St\"{u}ckelberg perturbations does not lead to higher time derivative in equations of motion, which suggests the absence of the BD ghost. By extending this method, we reconfirm that the BD ghost does not appear for arbitrary order of the perturbations at the decoupling limit in bi-gravity.

We investigate perturbations of a class of spherically symmetric solutions in massive gravity and bi-gravity. The background equations of motion for the particular class of solutions we are interested in reduce to a set of the Einstein equations with a cosmological constant. Thus, the solutions in this class include all the spherically symmetric solutions in general relativity, such as the Friedmann-Lema\^{i}tre-Robertson-Walker solution and the Schwarzschild (-de Sitter) solution, though the one-parameter family of two parameters of the theory admits such a class of solutions. We find that the equations of motion for the perturbations of this class of solutions also reduce to the perturbed Einstein equations at first and second order. Therefore, the stability of the solutions coincides with that of the corresponding solutions in general relativity. In particular, these solutions do not suffer from non-linear instabilities which often appear in the other cosmological solutions in massive gravity and bi-gravity.

It has been pointed out that the null energy condition can be violated stably in some non-canonical scalar-field theories. This allows us to consider the Galilean Genesis scenario in which the universe starts expanding from Minkowski spacetime and hence is free from the initial singularity. We use this scenario to study the early-time completion of inflation, pushing forward the recent idea of Pirtskhalava et al. We present a generic form of the Lagrangian governing the background and perturbation dynamics in the Genesis phase, the subsequent inflationary phase, and the graceful exit from inflation, as opposed to employing the effective field theory approach. Our Lagrangian belongs to a more general class of scalar-tensor theories than the Horndeski theory and Gleyzes-Langlois-Piazza-Vernizzi generalization, but still has the same number of the propagating degrees of freedom, and thus can avoid Ostrogradski instabilities. We investigate the generation and evolution of primordial perturbations in this scenario and show that one can indeed construct a stable model of inflation preceded by (generalized) Galilean Genesis.

The Horndeski theory is known as the most general scalar-tensor theory with second-order field equations. In this paper, we explore the bi-scalar extension of the Horndeski theory. Following Horndeski's approach, we determine all the possible terms appearing in the second-order field equations of the bi-scalar-tensor theory. We compare the field equations with those of the generalized multi-Galileons, and confirm that our theory contains new terms that are not included in the latter theory. We also discuss the construction of the Lagrangian leading to our most general field equations.

A time projection chamber (TPC) can be used to measure the polarization of gamma rays with excellent angular precision and sensitivity in the MeV-GeV energy range through the conversion of photons to e+e- pairs. The Hermetic ARgon POlarimeter (HARPO) prototype was built to demonstrate this concept. It was recently tested in the polarized photon beam at the NewSUBARU facility in Japan. We present this data-taking run, which demonstrated the excellent performance of the HARPO TPC.

We give new constraints on small-scale non-Gaussianity of primordial curvature perturbations by the use of anisotropies in acoustic reheating. Mixing of local thermal or local kinetic equilibrium systems with different temperatures yields a locally averaged temperature rise, which is proportional to the square of temperature perturbations damping in the photon diffusion scale. Such secondary temperature perturbations are indistinguishable from the standard temperature perturbations linearly coming from primordial curvature perturbations and hence should be subdominant compared to the standard ones. We show that small-scale higher order correlation functions (connected non-Gaussian and disconnected Gaussian parts) of primordial curvature perturbations can be probed by investigating auto power spectrum of the generated secondary perturbations and the cross power spectrum with the standard perturbations. This is simply because these power spectra come from higher order correlation functions of primordial curvature perturbations with non-linear parameters such as $f_{\rm NL}$ and $\tau_{\rm NL}$ since secondary temperature perturbations are second order effects. Thus, the observational results $l(l+1)C^{TT}_l\simeq 6\times 10^{-10}$ at large scales give a robust and universal upper bound on small-scale non-Gaussianities of primordial curvature perturbations.

The St\"uckelberg analysis of nonlinear massive gravity in the presence of a general fiducial metric is investigated. We develop a "covariant" formalism for the St\"uckelberg expansion by working with a local inertial frame, through which helicity modes can be characterized correctly. Within this covariant approach, an extended $\Lambda_3$ decoupling limit analysis can be consistently performed, which keeps $\bar{R}_{\mu\nu\rho\sigma}/m^2$ fixed with $\bar{R}_{\mu\nu\rho\sigma}$ the Riemann tensor of the fiducial metric. In this extended decoupling limit, the scalar mode $\pi$ acquires self-interactions due to the presence of the curvature of the fiducial metric. However, the equation of motion for $\pi$ remains of second order in derivatives, which extends the understanding of the absence of the Boulware Deser ghost in the case of a flat fiducial metric.

Shift-symmetry is essential to protect the flatness of the potential, even beyond the super-Planckian vacuum expectation value (VEV) for an inflaton field. The breaking of the shift-symmetry can yield potentials suitable for super-Planckian excursion of the inflaton. The aim of this paper is to illustrate that it is indeed possible to break the shift-symmetry dynamically within 4 dimensional supergravity prior to a long phase of inflation. Thanks to the shift-symmetry, the leading contribu- tion to the inflaton potential is free from the dangerous exponential factor even after its breaking, which is the main obstacle to realizing the super-Planckian inflation in supergravity. But, in our simple model, the resulting inflaton potential is a cosine type potential rather than the power-law one and it is difficult to realize a super-Planckian breaking scale unfortunately.

We newly introduce lift force into the aerodynamics of dust grains in the protoplanetary disk. Although many authors have so far investigated the effects of the drag force, gravitational force and electric force on the dust grains, the lift force has never been considered as a force exerted on the dust grains in the gas disk. If the grains are spinning and moving in the fluid, then the lift force is exerted on them. We show in this paper that the dust grains can be continuously spinning due to the frequent collisions so that the lift force continues to be exerted on them, which is valid in a certain parameter space where the grain size is larger than ~ 1 m and where the distance from the central star is larger than 1 AU for the minimum mass solar nebula. In addition, we estimate the effects of the force on the grain motion and obtain the result that the mean relative velocity between the grains due to the lift force is comparable to the gas velocity in the Kepler rotational frame when the Stokes number and lift-drag ratio are both ~ 1. This estimation is performed under the assumptions of the steady state and the isotropic spin angular momentum. We also estimate the mean relative velocity when the grains keep spinning and conclude that the lift force marginally affects the mean relative velocity in the minimum mass solar nebula. If there is a grain-concentrated part in the disk, the relative velocity due to the lift force may dominate there because of high collision rate.

We investigate whether the baryon asymmetry of the universe is explained in the framework of the supersymmetric extension of the Standard Model with R-parity violating interactions. It is shown that the Affleck-Dine mechanism naturally works via a trilinear interaction $LLE^c$, $LQD^c$, or $U^cD^cD^c$, if the magnitude of the coupling corresponding to the operator $\lambda$, $\lambda'$, or $\lambda''$ is sufficiently small. The formation of Q-balls and their subsequent evolution are also discussed. The present baryon asymmetry can be explained in the parameter region where R-parity is mildly violated $10^{-9}\lesssim\lambda,\lambda',\lambda''\lesssim 10^{-6}$ and the mass of the gravitino is relatively heavy $m_{3/2}\gtrsim 10^4\mathrm{GeV}$. On the other hand, it is difficult to explain the present baryon asymmetry for larger values of R-parity violating couplings $\lambda,\lambda',\lambda''\gtrsim 10^{-5}$, since Q-balls are likely to be destructed in the thermal environment and the primordial baryon number is washed away.

A moment approach for orbit determinations of astrometric binaries from astrometric observations alone has been recently studied for a low signal-to-noise ratio (Iwama et al. 2013, PASJ, 65, 2). With avoiding a direct use of the time-consuming Kepler equation, temporal information is taken into account to increase the accuracy of statistical moments. As numerical tests, 100 realizations are done and the mean and the standard deviation are also evaluated. For a semi-major axis, the difference between the mean of the recovered values and the true value decreases to less than a tenth in the case of $10000$ observed points. Therefore, the present moment approach works better than the previous one for the orbit determinations when one has a number of the observed points. The present approach is thus applicable to Cyg X-1.

We propose a new mechanism of generating the $\mu$ distortion in cosmic microwave background (CMB) originated from primordial gravitational waves. Such $\mu$ distortion is generated by the damping of the temperature anisotropies through the Thomson scattering, even on scales larger than that of Silk damping. This mechanism is in sharp contrast with that from the primordial curvature (scalar) perturbations, in which the temperature anisotropies mainly decay by Silk damping effects. We estimate the size of the $\mu$ distortion from the new mechanism, which can be used to constrain the amplitude of primordial gravitational waves on smaller scales independently from the CMB anisotropies, giving more wide-range constraint on their spectral index by combining the amplitude from the CMB anisotropies.

Quintessence, a scalar field model, has been proposed to account for the acceleration of the Universe at present. We discuss how accurately quintessence models are discriminated by future cosmological surveys, which include experiments of CMB, galaxy clustering, weak lensing, and the type Ia SNe surveys, by making use of the conventional parameterized dark energy models. We can see clear differences between the thawing and the freezing quintessence models at more than $1\sigma$ ($2\sigma$) confidence level as long as the present equation of state for quintessence is away from $-1$ as $w_X \gtrsim -0.95 (-0.90)$. However, it is found to be difficult to probe the effective mass squared for the potential in thawing models, whose signs are different between the quadratic and the cosine-type potentials. This fact may require us to invent a new estimator to distinguish quintessence models beyond the thawing and the freezing ones.

Effects of heavy fields on primordial spectra of curvature perturbations are discussed in inflationary models with a sudden turning trajectory. When heavy fields are excited after the sudden turn and oscillate around the bottom of the potential, the following two effects are generically induced: deformation of the inflationary background spacetime and conversion interactions between adiabatic and isocurvature perturbations, both of which can affect the primordial density perturbations. In this paper, we calculate primordial spectra in inflationary models with sudden turning potentials taking into account both of the two effects appropriately. We find that there are some non-trivial correlations between the two effects in the power spectrum and, as a consequence, the primordial scalar power spectrum has a peak around the scale exiting the horizon at the turn. Though both effects can induce parametric resonance amplifications, they are shown to be canceled out for the case with the canonical kinetic terms. The peak feature and the scale dependence of bispectra are also discussed.

We investigate the relation between the non-Gaussianities of the primordial perturbations and the sound speed of the tensor perturbations, that is, the propagation speed of the gravitational waves. We find that the sound speed of the tensor perturbations is directly related not to the auto-bispectrum of the tensor perturbations but to the cross-bispectrum of the primordial perturbations, especially, the scalar-tensor-tensor bispectrum. This result is in sharp contrast with the case of the scalar (curvature) perturbations, where their reduced sound speed enhances their auto-bispectrum. Our findings indicate that the scalar-tensor-tensor bispectrum can be a powerful tool to probe the sound speed of the tensor perturbations.

We investigate the structure of gravitational self-interactions of coherently oscillating axions in the general relativistic framework. A generic action for a massive scalar field in the Friedmann-Robertson-Walker background is first introduced based on the effective field theory approach to cosmological perturbations. Using the obtained setup, we evaluate the effective gravitational interaction of the massive scalar field, i.e. scalar quartic interactions mediated by metric perturbations. Applying the results to the system of dark matter axions, we estimate their self-interaction rate and discuss its implications for the axion Bose-Einstein condensate dark matter scenario. Leading contributions for the gravitational interactions of axions are given by the process mediated by the dynamical graviton field, which is essentially the Newtonian potential induced by fluctuations of the background fluids. We find that it leads to the same order of magnitude for the interaction rate of dark matter axions in the condensed regime, compared with the results of previous studies using the Newtonian approximation.

We derive the cosmological matching conditions for the homogeneous and isotropic background and for linear perturbations in Horndeski's most general second-order scalar-tensor theory. In general relativity, the matching is done in such a way that the extrinsic curvature is continuous across the transition hypersurface. This procedure is generalized so as to incorporate the mixing of scalar and gravity kinetic terms in the field equations of Horndeski's theory. Our matching conditions have a wide range of applications including the galilean genesis and the bounce scenarios, in which stable, null energy condition violating solutions play a central role. We demonstrate how our matching conditions are used in the galilean genesis scenario. In doing so, we extend the previous genesis models and provide a unified description of the theory admitting the solution that starts expanding from the Minkowski spacetime.

Higgs G-inflation is a Higgs inflation model with a generalized Galileon term added to the standard model Higgs field, which realizes inflation compatible with observations. Recently, it was claimed that the generalized Galileon term induces instabilities during the oscillation phase, and that the simplest Higgs G-inflation model inevitably suffers from this problem. In this paper, we extend the original Higgs G-inflation Lagrangian to a more general form, namely introducing a higher-order kinetic term and generalizing the form of the Galileon term, so that the Higgs field can oscillate after inflation without encountering instabilities. Moreover, it accommodates a large region of the n_s - r plane, most of which is consistent with current observations, leading us to expect the detection of B-mode polarization in the cosmic microwave background in the near future.

We provide the correspondence between the variables in the Jordan frame and those in the Einstein frame in scalar-tensor gravity and consider the frame-(in)dependence of the cosmological observables. In particular, we show that the cosmological observables/relations (redshift, luminosity distance, temperature anisotropies) are frame-independent. We also study the frame-dependence of curvature perturbations and find that the curvature perturbations are conformal invariant if the perturbation is adiabatic and the entropy perturbation between matter and the Brans-Dicke scalar is vanishing. The relation among various definitions of curvature perturbations in the both frames is also discussed, and the condition for the equivalence is clarified.

We revisit spectator field models including curvaton and modulated reheating scenarios, specifically focusing on their viability in the new Planck era, based on the derived expression for the spectral index in general spectator field models. Importantly, the recent Planck observations give strong preference to a red-tilted power spectrum, while the spectator field models tend to predict a scale-invariant one. This implies that, during inflation, either (i) the Hubble parameter varies significantly as in chaotic inflation, or (ii) a scalar potential for the spectator field has a relatively large negative curvature. Combined with the tight constraint on the non-Gaussianity, the Planck data provides us with rich implications for various spectator field models.

We propose a multi-field extension of (generalized) G-inflation, based on covariant multi-galileons and their generalization preserving second-order field equations. We compute the quadratic action for cosmological perturbations. By comparing the formulas for cosmological perturbations, it is highlighted that multi-field DBI galileon inflation is not included in the multi-field version of generalized G-inflation. Our result indicates that the generalized covariant multi-galileon theory is not the most general multi-scalar-tensor theory with second-order field equations.

We discuss implications of Planck results for models with local type non-Gaussianity. In light of the recent results of the Planck satellite, we constrain model parameters of several representative models and give the prediction of trispectrum, in particular, gNL. We also consider interesting possibilities that trispectrum appears as the first signature of the non-Gaussianities of the curvature perturbations, that is, fNL is small while gNL can be significantly large.

We apply the effective field theory approach to quasi-single field inflation, which contains an additional scalar field with Hubble scale mass other than inflaton. Based on the time-dependent spatial diffeomorphism, which is not broken by the time-dependent background evolution, the most generic action of quasi-single field inflation is constructed up to third order fluctuations. Using the obtained action, the effects of the additional massive scalar field on the primordial curvature perturbations are discussed. In particular, we calculate the power spectrum and discuss the momentum-dependence of three point functions in the squeezed limit for general settings of quasi-single field inflation. Our framework can be also applied to inflation models with heavy particles. We make a qualitative discussion on the effects of heavy particles during inflation and that of sudden turning trajectory in our framework.

We discuss the possibility that dark matter axions form a Bose-Einstein condensate (BEC) due to the gravitational self-interactions. The formation of BEC occurs in the condensed regime, where the transition rate between different momentum states is large compared to the energy exchanged in the transition. The time evolution of the quantum state occupation number of axions in the condensed regime is derived based on the in-in formalism. We recover the expression for the thermalization rate due to self interaction of the axion field, which was obtained in the other literature. It is also found that the leading order contributions for interactions between axions and other species vanish, which implies that the axion BEC does not give any significant modifications on standard cosmological parameters.

We compute the full bispectra, namely both auto- and cross- bispectra, of primordial curvature and tensor perturbations in the most general single-field inflation model whose scalar and gravitational equations of motion are of second order. The formulae in the limits of k-inflation and potential-driven inflation are also given. These expressions are useful for estimating the full bispectra of temperature and polarization anisotropies of the cosmic microwave background radiation.

We propose a model for the gamma-ray binary LS 5039 in which the X-ray emission is due to the inverse Compton (IC) process instead of the synchrotron radiation. Although the synchrotron model has been discussed in previous studies, it requires a strong magnetic field which leads to a severe suppression of the TeV gamma-ray flux in conflict with H.E.S.S. observations. In this paper, we calculate the IC emission by low energy electrons (\gamma_e \lesssim 10^3) in the Thomson regime. We find that IC emission of the low energy electrons can explain the X-ray flux and spectrum observed with Suzaku if the minimum Lorentz factor of injected electrons \gamma_min is around 10^3. In addition, we show that the Suzaku light curve is well reproduced if \gamma_min varies in proportion to the Fermi flux when the distribution function of injected electrons at higher energies is fixed. We conclude that the emission from LS 5039 is well explained by the model with the IC emission from electrons whose injection properties are dependent on the orbital phase. Since the X-ray flux is primarily determined by the total number of cooling electrons, this conclusion is rather robust, although some mismatches between the model and observations at the GeV band remain in the present formulation.

We discuss a supersymmetric version of DBI (Dirac-Born-Infeld) inflation, which is a typical inflation model in string cosmology. The supersymmetric DBI action together with a superpotential always leads to correction terms associated with the potential into the kinetic term, which drastically change the dynamics of DBI inflation. We find two significant features of supersymmetric DBI inflation. The first one is that ultra-relativistic motion is prohibited to cause inflation, which leads to order of unity sound velocity squared and hence small non-Gaussianities of primordial curvature perturbations. The second one is that the relation between the tensor-to-scalar ratio and the field variation is modified. Then, significant tensor-to-scalar ratio $r \gtrsim 0.01$ is possible even for sub-Planck variation of the field. These new features are in sharp contrast with those of the standard non-supersymmetric DBI inflation and hence have a lot of interest implications on upcoming observations of cosmic microwave background (CMB) anisotropies by the Planck satellite as well as direct detection experiments of gravitational waves like DECIGO and BBO.

We study Higgs inflation in the context of generalized G-inflation, i.e. the most general single-field inflation model with second-order field equations. The four variants of Higgs inflation proposed so far in the literature can be accommodated at one time in our framework. We also propose yet another class of Higgs inflation, the running Einstein inflation model, that can naturally arise from the generalized G-inflation framework. As a result, five Higgs inflation models in all should be discussed on an equal footing. Concise formulas for primordial fluctuations in these generalized Higgs inflation models are provided, which will be helpful to determine which model is favored from the future experiments and observations such as the Large Hadron Collider and the Planck satellite.

The observational relation between the density of baryon and dark matter in the Universe, $\Omega_{\rm DM}/\Omega_B\simeq 5$, is one of the most difficult problems to solve in modern cosmology. We discuss a scenario that explains this relation by combining the asymmetric dark matter scenario and the spontaneous baryogenesis associated with the flat direction in the supersymmetric standard model. A part of baryon asymmetry is transferred to charge asymmetry $D$ that dark matter carries, if a symmetry violating interaction that works at high temperature breaks not only $B-L$ but also $D$ symmetries simultaneously. In this case, the present number density of baryon and dark matter can be same order if the symmetric part of dark matter annihilates sufficiently. Moreover, the baryon number density can be enhanced as compared to that of dark matter if another $B-L$ violating interaction is still in thermal equilibrium after the spontaneous genesis of dark matter, which accommodates a TeV scale asymmetric dark matter model.

We find new, simple cosmological solutions with flat, open, and closed spatial geometries, contrary to the previous wisdom that only the open model is allowed. The metric and the St\"{u}ckelberg fields are given explicitly, showing nontrivial configurations of the St\"{u}ckelberg in the usual Friedmann-Lema\^{i}tre-Robertson-Walker coordinates. The solutions exhibit self-acceleration, while being free from ghost instabilities. Our solutions can accommodate inhomogeneous dust collapse represented by the Lema\^{i}tre-Tolman-Bondi metric as well. Thus, our results can be used not only to describe homogeneous and isotropic cosmology but also to study gravitational collapse in massive gravity.

We explicitly show the fully non-linear equivalence of the $\delta$N and the covariant formalisms for the superhorizon curvature perturbations, which enables us to safely evaluate the non-Gaussian quantities of the curvature perturbation in either formalism. We also discuss isocurvature perturbations in the covariant formalism and clarify the relation between the fully non-linear evolution of the curvature covector and that of the curvature perturbation for multiple interacting fluids.

We show that the modulated reheating mechanism can naturally account for dark radiation, whose existence is hinted by recent observations of the cosmic microwave background radiation and the primordial Helium abundance. In this mechanism, the inflaton decay rate depends on a light modulus which acquires almost scale-invariant quantum fluctuations during inflation. We find that the light modulus is generically produced by the inflaton decay and therefore a prime candidate for the dark radiation. Interestingly, an almost scale-invariant power spectrum predicted in the modulated reheating mechanism gives a better fit to the observation in the presence of the extra radiation. We discuss the production mechanism of the light modulus in detail taking account of its associated isocurvature fluctuations. We also consider a case where the modulus becomes the dominant component of dark matter.

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

We study generalized Galileons as a framework to develop the most general single-field inflation models ever, Generalized G-inflation, containing yet further generalization of G-inflation, as well as previous examples such as k-inflation, extended inflation, and new Higgs inflation as special cases. We investigate the background and perturbation evolution in this model, calculating the most general quadratic actions for tensor and scalar cosmological perturbations to give the stability criteria and the power spectra of primordial fluctuations. It is pointed out in the Appendix that the Horndeski theory and the generalized Galileons are equivalent. In particular, even the non-minimal coupling to the Gauss-Bonnet term is included in the generalized Galileons in a non-trivial manner.

We completely clarify the feature of primordial non-Gaussianities of tensor perturbations in generalized G-inflation, i.e., the most general single-field inflation model with second order field equations. It is shown that the most general cubic action for the tensor perturbation (gravitational wave) $h_{ij}$ is composed only of two contributions, one with two spacial derivatives and the other with one time derivative on each $h_{ij}$. The former is essentially identical to the cubic term that appears in Einstein gravity and predicts a squeezed shape, while the latter newly appears in the presence of the kinetic coupling to the Einstein tensor and predicts an equilateral shape. Thus, only two shapes appear in the graviton bispectrum of the most general single-field inflation model, which could open a new clue to the identification of inflationary gravitational waves in observations of cosmic microwave background anisotropies as well as direct gravitational wave detection experiments.

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.

We present a comprehensive study of primordial fluctuations generated from G-inflation, in which the inflaton Lagrangian is of the form $K(\phi, X)-G(\phi, X)\Box\phi$ with $X=-(\partial\phi)^2/2$. The Lagrangian still gives rise to second-order gravitational and scalar field equations, and thus offers a more generic class of single-field inflation than ever studied, with a richer phenomenology. We compute the power spectrum and the bispectrum, and clarify how the non-Gaussian amplitude depends upon parameters such as the sound speed. In so doing we try to keep as great generality as possible, allowing for non slow-roll and deviation from the exact scale-invariance.

In this review, we discuss inflation models based on supergravity. After explaining the difficulties in realizing inflation in the context of supergravity, we show how to evade such difficulties. Depending on types of inflation, we give concrete examples, particularly paying attention to chaotic inflation because the ongoing experiments like Planck might detect the tensor perturbations in near future. We also discuss inflation models in Jordan frame supergravity, motivated by Higgs inflation.

Recently spatial as well as temporal variations of the fine structure constant alpha have been reported. We show that a "runaway domain wall", which arises for the scalar field potential without minima, can account for such variations simultaneously. The time variation is induced by a runaway potential and the spatial variation is induced by the formation of a domain wall. The model is consistent with the current cosmological data and can be tested by the future experiments to test the equivalence principle.

A new class of inflation models within the context of G-inflation is proposed, in which the standard model Higgs boson can act as an inflaton thanks to Galileon-like non-linear derivative interaction. The generated primordial density perturbation is shown to be consistent with the present observational data. We also make a general discussion on potential-driven G-inflation models, and find a new consistency relation between the tensor-to-scalar ratio $r$ and the tensor spectral index $n_T$, $r = -32 \sqrt{6}n_T / 9$, which is crucial in discriminating the present models from standard inflation with a canonical kinetic term.

We propose a new class of inflation model, G-inflation, which has a Galileon-like nonlinear derivative interaction of the form $G(\phi, (\nabla\phi)^2)\Box\phi$ in the Lagrangian with the resultant equations of motion being of second order. It is shown that (almost) scale-invariant curvature fluctuations can be generated even in the exactly de Sitter background and that the tensor-to-scalar ratio can take a significantly larger value than in the standard inflation models, violating the standard consistency relation. Furthermore, violation of the null energy condition can occur without any instabilities. As a result, the spectral index of tensor modes can be blue, which makes it easier to observe quantum gravitational waves from inflation by the planned gravitational-wave experiments such as LISA and DECIGO as well as by the upcoming CMB experiments such as Planck and CMBpol.

We study time evolution of the $Q$ ball in thermal logarithmic potential using lattice simulations. As the temperature decreases due to the cosmic expansion, the thermal logarithmic term in the potential is eventually overcome by a mass term, and we confirm that the $Q$ ball transforms from the thick-wall type to the thin-wall type for a positive coefficient of radiative corrections to the mass term, as recently suggested. Moreover, we find that the $Q$ ball finally "melts down" when the $Q$-ball solution disappears. We also discuss the effects of this phenomenon on the detectability of gravitational waves from the $Q$-ball formation.

We study the detectability of the gravitational waves (GWs) from the Q-ball formation associated with the Affleck-Dine (AD) mechanism, taking into account both of the dilution effect due to Q-ball domination and of finite temperature effects. The AD mechanism predicts the formation of non-topological solitons, Q-balls, from which GWs are generated. Q-balls with large conserved charge $Q$ can produce a large amount of GWs. On the other hand, the decay rate of such Q-balls is so small that they may dominate the energy density of the universe, which implies that GWs are significantly diluted and that their frequencies are redshifted during Q-ball dominated era. Thus, the detectability of the GWs associated with the formation of Q-balls is determined by these two competing effects. We find that there is a finite but small parameter region where such GWs may be detected by future detectors such as DECIGO or BBO, only in the case when the thermal logarithmic potential dominates the potential of the AD field. Otherwise GWs from Q-balls would not be detectable even by these futuristic detectors: $\Omega_{\rm GW}^0<10^{-21}$. Unfortunately, for such parameter region the present baryon asymmetry of the universe can hardly be explained unless one fine-tunes A-terms in the potential. However the detection of such a GW background may give us an information about the early universe, for example, it may suggest that the flat directions with $B-L=0$ are favored.

We classify models generating large local-type non-Gaussianity into some categories by using some "consistency relations" among the non-linearity parameters f_{NL}^{local}, \tau_{NL}^{local} and g_{NL}^{local}, which characterize the size of bispectrum for the former and trispectrum for the later two. Then we discuss how one can discriminate models of large local-type non-Gaussianity with such relations. We first classify the models by using the ratio of \tau_{NL}^{local}/(6f_{NL}^{local}/5)^2, which is unity for "single-source" models and deviates from unity for "multi-source" ones. We can make a further classification of models in each category by utilizing the relation between f_{NL}^{local} and g_{NL}^{local}. Our classification suggests that observations of trispectrum would be very helpful to distinguish models of large non-Gaussianity and may reveal the generation mechanism of primordial fluctuations.

We show that current cosmic acceleration can be explained by an almost massless scalar field experiencing quantum fluctuations during primordial inflation. Provided its mass does not exceed the Hubble parameter today, this field has been frozen during the cosmological ages to start dominating the universe only recently. By using supernovae data, completed with baryonic acoustic oscillations from galaxy surveys and cosmic microwave background anisotropies, we infer the energy scale of primordial inflation to be around a few TeV, which implies a negligible tensor-to-scalar ratio of the primordial fluctuations. Moreover, our model suggests that inflation lasted for an extremely long period. Dark energy could therefore be a natural consequence of cosmic inflation close to the electroweak energy scale.

The emission mechanism of the gamma-ray binary LS5039 in energy bands of TeV, GeV, and X-ray is investigated. Observed light curves in LS5039 show that TeV and GeV fluxes anticorrelate and TeV and X-ray fluxes correlate. However, such correlated variations have not been explained yet reasonably at this stage. Assuming that relativistic electrons are injected constantly at the location of the compact object as a point source, and that they lose energy only by the inverse Compton (IC) process, we calculate gamma-ray spectra and light curves by the Monte Carlo method, including the full electromagnetic cascade process. Moreover, we calculated X-ray spectra and light curves by using the resultant electron distribution. As a result, we are able to reproduce qualitatively spectra and light curves observed by HESS, Fermi, and Suzaku for the inclination angle i = 30 dig and the index of injected electron distribution p = 2.5. We conclude that TeV-GeV anticorrelation is due to anisotropic IC scattering and anisotropic gamma-gamma absorption, and that TeV-X correlation is due to the dependence of IC cooling time on orbital phases. In addition, the constraint on the inclination angle implies that the compact object in LS5039 is a black hole.

In Horava's theory of gravity, Lorentz symmetry is broken in exchange for renormalizability, but the original theory has been argued to be plagued with problems associated with a new scalar mode stemming from the very breaking of Lorentz symmetry. Recently, Blas, Pujolas, and Sibiryakov have proposed a healthy extension of Horava gravity, in which the behavior of the scalar mode is improved. In this paper, we study scalar modes of cosmological perturbations in extended Horava gravity. The evolution of metric and density perturbations is addressed analytically and numerically. It is shown that for vanishing non-adiabatic pressure of matter the large scale evolution of cosmological perturbations converges to that described by a single constant, $\zeta$, which is an analog of a curvature perturbation on the uniform-density slicing commonly used in usual gravitational theories. The subsequent evolution is thus determined completely by the value of $\zeta$.

We derive the slow-roll conditions for a non-minimally coupled scalar field (extended quintessence) during the radiation/matter dominated era extending our previous results for thawing quintessence. We find that the ratio $\ddot\phi/3H\dot\phi$ becomes constant but negative, in sharp contrast to the ratio for the minimally coupled scalar field. We also find that the functional form of the equation of state of the scalar field asymptotically approaches that of the minimally coupled thawing quintessence.

In the non-relativistic theory of gravity recently proposed by Horava, the Hamiltonian constraint is not satisfied locally at each point in space. The absence of the local Hamiltonian constraint allows the system to have an extra dark-matter-like component as an integration constant. We discuss consequences of this fact in the context of cosmological perturbations, paying a particular attention to the large scale evolution of the curvature perturbation. The curvature perturbation is defined in a gauge invariant manner with this ``dark matter'' taken into account. We then clarify the conditions under which the curvature perturbation is conserved on large scales. This is done by using the evolution equations.

We investigate primordial non-Gaussianity and dark matter isocurvature fluctuations in the modulated reheating and the curvaton scenarios. In these scenarios, large non-Gaussianity can be generated, on the other hand, depending on how dark matter is produced, too large isocurvature fluctuations can also arise, which is inconsistent with current observations. In this paper, we study this issue in a mixed scenario where the curvature fluctuations can also be produced from the inflaton fluctuations as well as those from a light scalar field such as the modulus and the curvaton. We show that primordial fluctuations can be highly non-Gaussian without conflicting the current constraint on isocurvature fluctuations for such mixed scenarios. However, if the constraint on isocurvature fluctuations becomes severer as expected by the Planck satellite, $f_{\rm NL}$, a nonlinearity parameter for adiabatic fluctuations, should be very small as $f_{\rm NL} \lesssim 3$, which would give interesting implications for the generation mechanism of dark matter. Non-Gaussianity from isocurvature fluctuations is also discussed in these scenarios.

We investigate density fluctuations in a scenario with gravitino dark matter in the framework of modulated reheating, which is known to generate large non-Gaussianity. We show that gravitino dark matter is disfavored if primordial fluctuations have large local-type non-Gaussianity in this framework. We also briefly discuss the case with the curvaton mechanism and some other possible dark matter scenarios.

We show that adiabatic, super-Hubble, and almost scale invariant density fluctuations are produced by cosmic strings in a contracting universe. An essential point is that isocurvature perturbations produced by topological defects such as cosmic strings on super-Hubble scales lead to a source term which seeds the growth of curvature fluctuations on these scales. Once the symmetry has been restored at high temperatures, the isocurvature seeds disappear, and the fluctuations evolve as adiabatic ones in the expanding phase. Thus, cosmic strings may be resurrected as a mechanism for generating the primordial density fluctuations observed today.

In this paper, we put constraints on neutrino properties such as mass $m_{\nu}$ and degeneracy parameters $\xi_i$ from WMAP5 data and light element abundances by using a Markov chain Monte Carlo (MCMC) approach. In order to take consistently into account the effects of the degeneracy parameters, we run the Big Bang Nucleosynthesis code for each value of $\xi_i$ and the other cosmological parameters to estimate the Helium abundance, which is then used to calculate CMB anisotropy spectra instead of treating it as a free parameter. We find that the constraint on $m_{\nu}$ is fairly robust and does not vary very much even if the lepton asymmetry is allowed, and is given by $\sum m_\nu < 1.3 \rm eV$ ($95 % \rm C.L.$).

We study the effect of primordial isocurvature perturbations on non-Gaussian properties of CMB temperature anisotropies. We consider generic forms of the non-linearity of isocurvature perturbations which can be applied to a wide range of theoretical models. We derive analytical expressions for the bispectrum and the Minkowski Functionals for CMB temperature fluctuations to describe the non-Gaussianity from isocurvature perturbations. We find that the isocurvature non-Gaussianity in the quadratic isocurvature model, where the isocurvature perturbation S is written as a quadratic function of the Gaussian variable sigma, S=sigma^2-<sigma^2>, can give the same signal-to-noise as f_NL=30 even if we impose the current observational limit on the fraction of isocurvature perturbations contained in the primordial power spectrum alpha. We give constraints on isocurvature non-Gaussianity from Minkowski Functionals using WMAP 5-year data. We do not find a significant signal of the isocurvature non-Gaussianity. For the quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature fraction alpha<0.070 (95% CL) for a scale invariant spectrum which is comparable to the limit obtained from the power spectrum.

In modern terahertz (THz) sensing and imaging spectroscopy, water is considered a nemesis to be avoided due to strong absorption in the THz frequency range. Here we report the first experimental demonstration and theoretical implications of using femtosecond laser pulses to generate intense broadband THz emission from water vapor. When we focused an intense laser pulse in water vapor contained in a gas cell or injected from a gas jet nozzle, an extraordinarily strong THz field from optically excited water vapor is observed. Water vapor has more than 50% greater THz generation efficiency than dry nitrogen. It had previously been assumed that the nonlinear generation of THz waves in this manner primarily involves a free-electron plasma, but we show that the molecular structure plays an essential role in the process. In particular, we found that THz wave generation from H2O vapor is significantly stronger than that from D2O vapor. Vibronic activities of water cluster ions, occurring naturally in water vapor, may possibly contribute to the observed isotope effect along with rovibrational contributions from the predominant monomers.

As an extension of our previous study, we derive slow-roll conditions for multiple scalar fields which are non-minimally coupled with gravity and for generalized gravity theories of the form $f(\phi,R)$. We provide simple formulae of the spectral indices of scalar/tensor perturbations in terms of the slow-roll parameters.

A simple model is discussed to give rise to successful reheating in chaotic D-term inflation with a quadratic inflaton potential, introducing a trilinear coupling in the K\"ahler potential. Leptogenesis through the inflaton decay is also discussed in this model.

We derive slow-roll conditions for a scalar field which is non-minimally coupled with gravity in a consistent manner and express spectral indices of scalar/tensor perturbations in terms of the slow-roll parameters. The conformal invariance of the curvature perturbation is proved without linear approximations. Rapid-roll conditions are also derived, and the relation with the slow-roll conditions is discussed.

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.

We study non-Gaussianity, the spectral index of primordial scalar fluctuations and tensor modes in models where fluctuations from the inflaton and the curvaton can both contribute to the present cosmic density fluctuations. Even though simple single-field inflation models generate only tiny non-Gaussianity, if we consider such a mixed scenario, large non-Gaussianity can be produced. Furthermore, we study the inflationary parameters such as the spectral index and the tensor-to-scalar ratio in this kind of models and discuss in what cases models predict the spectral index and tensor modes allowed by the current data while generating large non-Gaussianity, which may have many implications for model-buildings of the inflationary universe.

We investigate non-Gaussianity in the modulated reheating scenario where fluctuations of the decay rate of the inflaton generate adiabatic perturbations, paying particular attention to the non-linearity parameters $f_{\rm NL}, \tau_{\rm NL}$ and $g_{\rm NL}$ as well as the scalar spectral index and tensor-to-scalar ratio which characterize the nature of the primordial power spectrum. We also take into account the pre-existing adiabatic perturbations produced from the inflaton fluctuations. It has been known that the non-linearity between the curvature perturbations and the fluctuations of the decay rate can yield non-Gaussianity at the level of $f_{\rm NL} \sim \mathcal{O}(1)$, but we find that the non-linearity between the decay rate and the modulus field which determines the decay rate can generate much greater non-Gaussianity. We also discuss a consistency relation among non-linearity parameters which holds in the scenario and find that the modulated reheating yields a different one from that of the curvaton model. In particular, they both can yield a large positive $f_{\rm NL}$ but with a different sign of $g_{\rm NL}$. This provides a possibility to discriminate these two competitive models by looking at the sign of $g_{\rm NL}$. Furthermore, we work on some concrete inflation models and investigate in what cases models predict the spectral index and the tensor-to-scalar ratio allowed by the current data while generating large non-Gaussianity, which may have many implications for model-buildings of the inflationary universe.

We present a new model of D-term dominated chaotic inflation in supergravity. The F-flat direction present in this model is lifted by the dominant D-term, which leads to chaotic inflation and subsequent reheating. No cosmic string is formed after inflation because the U(1) gauge symmetry is broken during inflation. The leptogenesis scenario via the inflaton decay in our D-term chaotic inflation scenario is also discussed.

We investigate the non-Gaussianity of primordial curvature perturbation in the modulated reheating scenario where the primordial perturbation is generated due to the spacial fluctuation of the inflaton decay rate to radiation. We use the $\delta N$ formalism to evaluate the trispectrum of curvature perturbation as well as its bispectrum. We give expressions for three non-linear parameters $f_{NL}, \tau_{NL}$ and $g_{NL}$ in the modulated reheating scenario. If the intrinsic non-Gaussianity of scalar field fluctuations and third derivative of the decay rate with respect to scalar fields are negligibly small, $g_{NL}$ has at least the same order of magnitude as $f_{NL}$. We also give general inequality between $f_{NL}$ and $\tau_{NL}$ which is true for other inflationary scenarios as long as primordial non-Gaussianity comes from super-horizon evolution.

Models where the dark energy is a scalar field with a non-standard Dirac-Born-Infeld (DBI) kinetic term are investigated. Scaling solutions are studied and proven to be attractors. The corresponding shape of the brane tension and of the potential is also determined and found to be, as in the standard case, either exponentials or power-law of the DBI field. In these scenarios, in contrast to the standard situation, the vacuum expectation value of the field at small redshifts can be small in comparison to the Planck mass which could be an advantage from the model building point of view. This situation arises when the present-day value of the Lorentz factor is large, this property being per se interesting. Serious shortcomings are also present such as the fact that, for simple potentials, the equation of state appears to be too far from the observational favored value -1. Another problem is that, although simple stringy-inspired models precisely lead to the power-law shape that has been shown to possess a tracking behavior, the power index turns out to have the wrong sign. Possible solutions to these issues are discussed.

The hard X-ray detector (HXD) on board the X-ray satellite Suzaku is designed to have a good timing capability with a 61 $\mu$s time resolution. In addition to detailed descriptions of the HXD timing system, results of in-orbit timing calibration and performance of the HXD are summarized. The relative accuracy of time measurements of the HXD event was confirmed to have an accuracy of $1.9\times 10^{-9}$ s s$^{-1}$ per day, and the absolute timing was confirmed to be accurate to 360 $\mu$s or better. The results were achieved mainly through observations of the Crab pulsar, including simultaneous ones with RXTE, INTEGRAL, and Swift.

Even though the chaotic inflation is one of the most popular inflation models for its simple dynamics and compelling resolutions to the initial condition problems, its realization in supergravity has been considered a challenging task. We discuss how the chaotic inflation dominated by the D-term can be induced in supergravity, which would give a new perspective on the inflation model building in supergravity.

We calculate the signal-to-noise ratio (SNR) of the stochastic gravitational-wave background in an extreme case that its spectrum has a sharp falloff with its amplitude close to the detection threshold. Such a spectral feature is a characteristic imprint of the change in the number of relativistic degrees of freedom on the stochastic background generated during inflation in the early Universe. We find that, although SNR is maximal with the correct template which is proportional to the assumed real spectrum, its sensitivity to the shape of template is fairly weak indicating that a simple power-law template is sufficient to detect the signature.

We show how leptogenesis can occur at the TeV scale with neutrinos that possess almost purely Dirac masses and negligible Majorana mass contributions as a consequence of the small wavefunction overlap in a warped fifth dimension. Lepton number violation at the Planck scale is introduced via a Majorana mass term on the Planck brane. Such a Majorana mass operator leads to the small mass splitting of otherwise degenerate Kaluza-Klein excited states on the TeV brane. This tiny mass splitting can compensate for the small Yukawa couplings to give a CP asymmetry large enough to produce the sufficient baryon asymmetry from the decay of the nearly degenerate neutrino Kaluza-Klein states. In this way the standard baryogenesis via leptogenesis scenario can naturally occur at the TeV scale without the need for a high mass scale.

We show that axinos, which are dominantly generated by the decay of the next-to-lightest supersymmetric particles produced from the leptonic $Q$-ball ($L$-ball), become warm dark matter suitable for the solution of the missing satellite problem and the cusp problem. In addition, $\Omega_b - \Omega_{DM}$ coincidence is naturally explained in this scenario.

We consider electroweak baryogenesis mediated by embedded domain walls. Embedded domain walls originating from a symmetry breaking phase transition are stabilized by thermal plasma effects, so that the electroweak symmetry is unbroken in their cores. For this reason, the cosmological evolution of such domain walls can generate a sufficiently large baryon asymmetry, irrespective of the order of the electroweak phase transition. For embedded domain walls, the condition that the energy of the universe not be dominated by the energy of the domain walls is relaxed significantly, and it is shown to be compatible with our scenario of electroweak baryogenesis.

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.

Recent astrophysical observations indicate that the proton-electron mass ratio and the fine structure constant have gone through nontrivial time evolution. We discuss their time variation in the context of a dilaton runaway scenario with gauge coupling unification at the string scale $M_{\rm s}$. We show that the choice of adjustable parameters allows them to fit the same order magnitude of both variations and their (opposite) signs in such a scenario.

Effects of neutrino free streaming is evaluated on the primordial spectrum of gravitational radiation taking both neutrino chemical potential and masses into account. The former or the lepton asymmetry induces two competitive effects, namely, to increase anisotropic pressure, which damps the gravitational wave more, and to delay the matter-radiation equality time, which reduces the damping. The latter effect is more prominent and a large lepton asymmetry would reduce the damping. We may thereby be able to measure the magnitude of lepton asymmetry from the primordial gravitational wave spectrum.

We study the effect of free streaming in a universe with cosmic strings with time-varying tension as well as with constant tension. Although current cosmological observations suggest that fluctuation seeded by cosmic strings cannot be the primary source of cosmic density fluctuation, some contributions from them are still allowed. Since cosmic strings actively produce isocurvature fluctuation, the damping of small scale structure via free streaming by dark matter particles with large velocity dispersion at the epoch of radiation-matter equality is less efficient than that in models with conventional adiabatic fluctuation. We discuss its implications to the constraints on the properties of particles such as massive neutrinos and warm dark matter.

Any one-field inflation is actually realized in a multifield configuration because the inflaton must have couplings with other fields to reheat the universe and is coupled to all other fields at least gravitationally. In all single inflaton models, it is explicitly or implicitly assumed that the heavier fields are stuck to their potential minima during inflation, which are time-dependent in general. We present a formalism to calculate curvature perturbations in such a time-dependent background and show that the proper expression can be obtained using a single-field analysis with a reduced potential in which all these heavy fields are situated at their respective, time-dependent minima. Our results provide a firm ground on the conventional calculation.

We investigate cosmological evolution and implications of cosmic strings with time-dependent tension. We derive basic equations of time development of the correlation length and the velocity of such strings, based on the one scale model. Then, we find that, in the case where the tension depends on some power of the cosmic time, cosmic strings with time-dependent tension goes into the scaling solution if the power is lower than a critical value. We also discuss cosmic microwave background anisotropy and matter power spectra produced by these strings. The constraints on their tensions from the Wilkinson microwave anisotropy probe (WMAP) three year data and Sloan digital sky survey (SDSS) data are also given.

Cosmological issues of the gravitino production by the decay of a heavy scalar field $X$ are examined, assuming that the damped coherent oscillation of the scalar once dominates the energy of the universe. The coupling of the scalar field to a gravitino pair is estimated both in spontaneous and explicit supersymmetry breaking scenarios, with the result that it is proportional to the vacuum expectation value of the scalar field in general. Cosmological constraints depend on whether the gravitino is stable or not, and we study each case separately. For the unstable gravitino with $M_{3/2} \sim$ 100GeV--10TeV, we obtain not only the upper bound, but also the lower bound on the reheating temperature after the $X$ decay, in order to retain the success of the big-bang nucleosynthesis. It is also shown that it severely constrains the decay rate into the gravitino pair. For the stable gravitino, similar but less stringent bounds are obtained to escape the overclosure by the gravitinos produced at the $X$ decay. The requirement that the free-streaming effect of such gravitinos should not suppress the cosmic structures at small scales eliminates some regions in the parameter space, but still leaves a new window of the gravitino warm dark matter. Implications of these results to inflation models are discussed. In particular, it is shown that modular inflation will face serious cosmological difficulty when the gravitino is unstable, whereas it can escape the constraints for the stable gravitino. A similar argument offers a solution to the cosmological moduli problem, in which the moduli is relatively heavy while the gravitino is light.

We propose a new scenario of non-thermal production of neutralino cold dark matter, in which the overproduction problem of lightest supersymmetric particles (LSPs) in the standard thermal history is naturally solved. The mechanism requires a heavy modulus field which decays mainly to ordinary particles releasing large entropy to dilute gravitinos produced just after inflation and thermal relics of LSPs. Significant amount of gravitinos are also pair-produced at the decay, which subsequently decay into the neutralinos. We identify the regions of the parameter space in which the requisite abundance of the neutralino dark matter is obtained without spoiling the big-bang nucleosynthesis by injection of hadronic showers from gravitino decay. The neutralino abundance obtained in this mechanism is insensitive to the details of the superparticle mass spectrum, unlike the standard thermal abundance. We also briefly mention the testability of the scenario in future experiments.

We discuss the cosmological evolution of cosmic strings with time-dependent tension. We show that, in the case that the tension changes as a power of time, the cosmic string network obeys the scaling solution: the characteristic scale of the string network grows with the time. But due to the time dependence of the tension, the ratio of the energy density of infinite strings to that of the background universe is {\it not} necessarily constant.

The brane world scenario is a new approach to resolve the problem on how to compactify the higher dimensional spacetime to our 4-dimensional world. One of the remarkable features of this scenario is the higher dimensional effects in classical gravitational interactions at short distances. Due to this feature, there are black string solutions in our 4-dimensional world. In this paper, assuming the simplest model of complex minimally coupled scalar field with the local U(1) symmetry, we show a possibility of black-string formation by merging processes of type I long cosmic strings in our 4-dimensional world. No fine tuning for the parameters in the model might be necessary.

It is shown that in a smooth hybrid inflation model in supergravity adiabatic fluctuations with a running spectral index with $\ns >1$ on a large scale and $\ns <1$ on a smaller scale can be naturally generated, as favored by the first-year data of WMAP. It is due to the balance between the nonrenormalizable term in the superpotential and the supergravity effect. However, since smooth hybrid inflation does not last long enough to reproduce the central value of observation, we invoke new inflation after the first inflation. Its initial condition is set dynamically during smooth hybrid inflation and the spectrum of fluctuations generated in this regime can have an appropriate shape to realize early star formation as found by WMAP. Hence two new features of WMAP observations are theoretically explained in a unified manner.

We study the cosmological consequences of generic scalar fields like moduli which decay only through gravitationally suppressed interactions. We consider a new production mechanism of gravitinos from moduli decay, which might be more effective than previously known mechanisms, and calculate the final gravitino-to-entropy ratio to compare with the constraints imposed by successful big bang nucleosynthesis (BBN) etc., taking possible hadronic decays of gravitinos into account. We find the modulus mass smaller than $\sim 10^4$ TeV is excluded. On the other hand, inflation models with high reheating temperatures $T_{R,\rm inf} \sim 10^{16}$ GeV can be compatible with BBN thanks to the late-time entropy production from the moduli decay if model parameters are appropriately chosen.

We discuss a spontaneous baryogenesis mechanism in flat directions. After identifying the Nambu-Goldstone mode which derivatively couples to the associated $U$(1) current and rotates due to the A-term, we show that spontaneous baryogenesis can be naturally realized in the context of the flat direction. As applications, we discuss two scenarios of baryogenesis in detail. One is baryogenesis in a flat direction with a vanishing $B-L$ charge, especially, with neither baryon nor lepton charge, which was recently proposed by Chiba and the present authors. The other is a baryogenesis scenario compatible with a large lepton asymmetry.

The recent high precision maps of cosmic microwave anisotropies combined with measurements of the galaxy power spectrum from new large-scale redshift surveys have allowed stringent bounds on the sum of the neutrino masses to be placed. The past analyses, however, have implicitly assumed that the spectrum of primordial density fluctuations is adiabatic and coherent, as predicted in the simplest models of inflation. In this paper, we show that the limits hold even if the assumption on the primordial power spectrum is relaxed to allow for a contribution of nonadiabatic, incoherent fluctuations such as would be predicted by topological defects.

We present a new mechanism of spontaneous baryogenesis. Usually such mechanisms require a derivative coupling between a scalar field and baryon current. In our model, the scalar field instead derivatively couples to a current associated with some global symmetry $U(1)_Q$. In this case, baryogenesis is still possible provided that an interaction exists, which violates the baryon and $U(1)_Q$ symmetries simultaneously. As a concrete example, we discuss baryogenesis in a flat direction with neither baryon nor lepton charge.

We propose an inflation model in supergravity, in which chaotic and hybrid inflation occurs successively, followed by new inflation. During hybrid inflation, adiabatic fluctuations with a running spectral index with $\ns >1$ on a large scale and $\ns <1$ on a smaller scale are generated, as favored by recent results of the first year Wilkinson Microwave Anisotropy Probe. The initial condition of new inflation is also set dynamically during hybrid inflation, and its duration and the amplitude of density fluctuations take appropriate values to help early star formation to realize early reionization.

We propose a new scenario to simultaneously explain a large lepton asymmetry and a small baryon asymmetry. We consider a rolling scalar field and its derivative coupling to the lepton number current. The presence of an effective nonzero time derivative of the scalar field leads to $CPT$ violation so that the lepton asymmetry can be generated even in thermal equilibrium as pointed out by Cohen and Kaplan. In this model, the lepton asymmetry varies with time. In particular, we consider the case where it grows with time. The final lepton asymmetry is determined by the decoupling of the lepton number violating interaction and can be as large as order unity. On the other hand, if the decoupling takes place after the electroweak phase transition, a small baryon asymmetry is obtained from the small lepton asymmetry at that time through sphaleron effects. We construct a model in which a rolling scalar field is identified with a quintessence field.

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 discuss spontaneous baryogenesis in the warm inflation scenario. In contrast with standard inflation models, radiation always exists in the warm inflation scenario, and the inflaton must be directly coupled to it. Also, the transition to the post-inflationary radiation dominated phase is smooth and the entropy is not significantly increased at the end of the period of inflation. In addition, after the period of warm inflation ends, the inflaton does not oscillate coherently but slowly rolls. We show that as a consequence of these features of warm inflation, the scenario can well accommodate the spontaneous baryogenesis mechanism, provided that the decoupling temperature $T_{D}$ of the baryon or the $B-L$ violating interactions is higher than the temperature of radiation during the late stages of inflation.

We clarify the quantitative nature of the cosmological evolution of a global string network, that is, the energy density, peculiar velocity, velocity squared, Lorentz factor, formation rate of loops, and emission rate of Nambu-Goldstone bosons, based on a new type of numerical simulation of scalar fields in Eulerian meshes. We give a detailed explanation of a method to extract the above-mentioned quantities to characterize string evolution by analyzing scalar fields in Eulerian meshes from a Lagrangian viewpoint. We confirm our previous claim that the number of long strings per horizon volume in a global string network is smaller than in the case of a local string network by a factor of $\sim$10 even under cosmological situations, and its reason is clarified.

We establish a method to trace the Lagrangian evolution of extended objects consisting of a multicomponent scalar field in terms of a numerical calculation of field equations in three dimensional Eulerian meshes. We apply our method to the cosmological evolution of global strings and evaluate the energy density, peculiar velocity, Lorentz factor, formation rate of loops, and emission rate of Nambu-Goldstone (NG) bosons. We confirm the scaling behavior with a number of long strings per horizon volume smaller than the case of local strings by a factor of $\sim$ 10. The strategy and the method established here are applicable to a variety of fields in physics.

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 investigate the scaling property of global monopoles in the expanding universe. By directly solving the equations of motion for scalar fields, we follow the time development of the number density of global monopoles in the radiation dominated (RD) universe and the matter dominated (MD) universe. It is confirmed that the global monopole network relaxes into the scaling regime and the number per hubble volume is a constant irrespective of the cosmic time. The number density $n(t)$ of global monopoles is given by $n(t) \simeq (0.43\pm0.07) / t^{3}$ during the RD era and $n(t) \simeq (0.25\pm0.05) / t^{3}$ during the MD era. We also examine the peculiar velocity $v$ of global monopoles. For this purpose, we establish a method to measure the peculiar velocity by use of only the local quantities of the scalar fields. It is found that $v \sim (1.0 \pm 0.3)$ during the RD era and $v \sim (0.8 \pm 0.3)$ during the MD era. By use of it, a more accurate analytic estimate for the number density of global monopoles is obtained.

In this talk we propose a natural double inflation model in supergravity. Chaotic inflation first takes place by virtue of the Nambu-Goldstone-like shift symmetry. During chaotic inflation, an initial value of second inflation (new inflation) is set, which is adequately far from the local maximum of the potential due to the small linear term in the K\"ahler potential. Then, primordial fluctuations within the present horizon scale may be produced during both inflations. Primordial fluctuations responsible for anisotropies of the cosmic microwave background radiation and the large scale structure are produced during chaotic inflation, while fluctuations on smaller scales are produced during new inflation. Because of the peculiar nature of new inflation, they can become as large as $10^{-1}$-$10^{-2}$, which may lead to the formation of primordial black holes.

We investigate the cosmological evolution of global monopoles in the radiation dominated (RD) and matter dominated (MD) universes by numerically solving field equations of scalar fields. It is shown that the global monopole network relaxes into the scaling regime, unlike the gauge monopole network. The number density of global monopoles is given by $n(t) \simeq (0.43\pm0.07) / t^{3}$ during the RD era and $n(t) \simeq (0.25\pm0.05) / t^{3}$ during the MD era. Thus, we have confirmed that density fluctuations produced by global monopoles become scale invariant and are given by $\delta \rho \sim 7.2(5.0) \sigma^{2} / t^{2} $ during the RD (MD) era, where $\sigma$ is the breaking scale of the symmetry.

We investigate the recently proposed natural double inflation model in supergravity. Chaotic inflation first takes place by virtue of the Nambu-Goldstone-like shift symmetry. During chaotic inflation, an initial value of second inflation (new inflation) is set, which is adequately far from the local maximum of the potential due to the small linear term in the K\"ahler potential. Then, primordial fluctuations within the present horizon scale may be produced during both inflations. Primordial fluctuations responsible for anisotropies of the cosmic microwave background radiation and the large scale structure are produced during chaotic inflation, while fluctuations on smaller scales are produced during new inflation. Because of the peculiar nature of new inflation, they can become as large as $10^{-1}$-$10^{-2}$, which may lead to the formation of primordial black holes.

We propose a natural double inflation model in supergravity. In this model, chaotic inflation first takes place by virtue of the Nambu-Goldstone-like shift symmetry, which guarantees the absence of the exponential factor in the potential for the inflaton field. During chaotic inflation, an initial value of the second inflation (new inflation) is set. In this model, the initial value of new inflation can be adequately far from the local maximum of the potential for new inflation due to the small linear term of the inflaton in the K\"ahler potential. Therefore, the primordial fluctuations within the present horizon scale may be attributed to both inflations; that is, the first chaotic inflation produces the primordial fluctuations on the large cosmological scales while the second new inflation on the smaller scales. The successive decay of the inflaton for new inflation leads to a reheating temperature low enough to avoid the overproduction of gravitinos in a wide range of the gravitino mass.

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 propose a self-consistent scenario of new inflation in supergravity. Chaotic inflation first takes place around the Planck scale, which solves the longevity problem, namely, why the universe can live much beyond the Planck time, and also gives an adequate initial condition for new inflation. Then, new inflation lasts long enough to generate primordial fluctuations for the large scale structure, which generally has a tilted spectrum with the spectral index $n_{s} < 1$. The successive decay of the inflaton leads to the reheating temperature low enough to avoid the overproduction of gravitinos in a wide range of the gravitino mass.

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.

Recently, a novel class of models for inflation has been found in which the inflationary dynamics is driven solely by (non-canonical) kinetic terms rather than by potential terms. As an obvious extension, we show that a scalar field with non-canonical kinetic terms alone behaves like an energy component which is time-varying and has negative pressure presently, i.e. quintessence. We present a model which has a constant equation of state, that is, a ``kinetic'' counterpart of the Ratra-Peebles model of a quintessence field with a potential term. We make clear the structure of the phase plane and show that the quintessential solution is a late-time attractor. We also give a model for the ``phantom'' component which has an equation of state with $w=p/\rho <-1$.

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 present a new type of one-field models for open inflation utilizing a nonminimally coupled scalar field with polynomial potentials in which the Coleman-de Luccia instanton does exist and slow-roll inflation after the bubble nucleation is realized.

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.

We investigate the evolution of the global string network in the radiation dominated universe by use of numerical simulations in 3+1 dimensions. We find that the global string network settles down to the scaling regime where the energy density of global strings, $\rho_{s}$, is given by $\rho_{s} = \xi \mu / t^2$ with $\mu$ the string tension per unit length and the scaling parameter, $\xi \sim (0.9-1.3)$, irrespective of the cosmic time. We also find that the loop distribution function can be fitted with that predicted by the so-called one scale model. Concretely, the number density, $n_{l}(t)$, of the loop with the length, $l$, is given by $n_{l}(t) = \nu/[t^{3/2} (l + \kappa t)^{5/2}]$ where $\nu \sim 0.0865$ and $\kappa$ is related with the Nambu-Goldstone(NG) boson radiation power from global strings, $P$, as $P = \kappa \mu$ with $\kappa \sim 0.535$. Therefore, the loop production function also scales and the typical scale of produced loops is nearly the horizon distance. Thus, the evolution of the global string network in the radiation dominated universe can be well described by the one scale model in contrast with that of the local string network.

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 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 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.

A simple hadronic axion model is proposed in the framework of gauge-mediated supersymmetry breaking. Dynamics of Peccei-Quinn symmetry breaking is governed by supersymmetry breaking effects and the Peccei-Quinn breaking scale $f_{PQ}$ is inversely proportional to the gravitino mass. The gravitino mass range which corresponds to the axion window $f_{PQ} \simeq 10^{9}$ GeV -- $10^{13}$ GeV lies in the region predicted by gauge-mediated supersymmetry breaking models. The model is also shown to be cosmologically viable.

We present a formalism to calculate the probability distribution function of a scalar field coarse-grained over some spatial scales with a Gaussian filter at finite temperature. As an application, we investigate the role of subcritical fluctuations in the electroweak phase transition in the minimal standard model. It is concluded that the universe was in a mixed state of true and false vacua already at the critical temperature.

We analyse the formation of cosmic structures in models where the dark matter is dominated by light gravitinos with mass of $ 100$ eV -- 1 keV, as predicted by gauge-mediated supersymmetry (SUSY) breaking models. After evaluating the number of degrees of freedom at the gravitinos decoupling ($g_*$), we compute the transfer function for matter fluctuations and show that gravitinos behave like warm dark matter (WDM) with free-streaming scale comparable to the galaxy mass scale. We consider different low-density variants of the WDM model, both with and without cosmological constant, and compare the predictions on the abundances of neutral hydrogen within high-redshift damped Ly--$\alpha$ systems and on the number density of local galaxy clusters with the corresponding observational constraints. We find that none of the models satisfies both constraints at the same time, unless a rather small $\Omega_0$ value ($\mincir 0.4$) and a rather large Hubble parameter ($\magcir 0.9$) is assumed. Furthermore, in a model with warm + hot dark matter, with hot component provided by massive neutrinos, the strong suppression of fluctuation on scales of $\sim 1\hm$ precludes the formation of high-redshift objects, when the low--$z$ cluster abundance is required. We conclude that all different variants of a light gravitino DM dominated model show strong difficulties for what concerns cosmic structure formation. This gives a severe cosmological constraint on the gauge-mediated SUSY breaking scheme.

We investigate whether the universe was homogeneously in the false vacuum state at the critical temperature of a weakly first-order phase transition such as the electroweak phase transition in terms of a series of numerical simulations of a phenomenological Langevin equation, whose noise term is derived from the effective action but the dissipative term is set so that the fluctuation-dissipation relation is met. The correlation function of the noise terms given by a non-equilibrium field theory has a distinct feature if it originates from interactions with a boson or with a fermion. The spatial correlation function of noises from a massless boson damps with a power-law, while the fermionic noises always damp exponentially above the inverse-temperature scale. In the simulation with one-loop effective potential of the Higgs field, the latter turns out to be more effective to disturb the homogeneous field configuration. Since noises of the both types are present in the electroweak phase transition, our results suggest that conventional picture of a phase transition, namely, nucleation of critical bubbles in a homogeneous background does not apply or the one-loop approximation breaks down in the standard model.

In theories with a gauge-mediated mechanism of supersymmetry breaking the gravitino is likely to be the lightest superparticle and, hence, a candidate for dark matter. We show that the decay of the next-to-lightest superparticle into a gravitino can yield a non-thermal population of gravitinos which behave as a hot dark matter component. Together with the warm component, which is provided by the population of gravitinos of thermal origin, they can give rise to viable schemes of mixed dark matter. This realization has some specific and testable features both in particle physics and astrophysics. We outline under which conditions the mechanism remains viable even when R parity is broken.

Thermal history after inflation is studied in a chaotic inflation model with supersymmetric couplings of the inflaton to matter fields. Time evolution equation is solved in a formalism that incorporates both the back reaction of particle production and the cosmological expansion. The effect of the parametric resonance gives rise to a rapid initial phase of the inflaton decay followed by a slow stage of the Born term decay. Thermalization takes place immediately after the first explosive stage for a medium strength of the coupling among created particles. As an application we calculate time evolution of the gravitino abundance that is produced by ordinary particles directly created from the inflaton decay, which typically results in much more enhanced yield than what a naive estimate based on the Born term would suggest.

Large amplitude oscillation of cosmic field that may occur right after inflation and in the decay process of weakly interacting fields gives rise to violent particle production via the parametric resonance. In the large amplitude limit the problem of back reaction against the field oscillation is solved and the energy spectrum of created particles is determined in a semi-classical approximation. For large enough coupling or large enough amplitude the resulting energy spectrum is broadly distributed, implying larger production of high energy particles than what a simple estimate of the reheating temperature due to the Born formula would suggest.