Late-time cosmology in the extended cuscuton theory is studied, in which gravity is modified while one still has no extra dynamical degrees of freedom other than two tensor modes. We present a simple example admitting analytic solutions for the cosmological background evolution that mimics $\Lambda$CDM cosmology. We argue that the extended cuscuton as dark energy can be constrained, like usual scalar-tensor theories, by the growth history of matter density perturbations and the time variation of Newton's constant.

We show that an axionlike field coupled to a new confining gauge group can generate the primordial density perturbation in the post-inflation universe. The axion decay constant and strong coupling scale are uniquely determined by observations of the density perturbation, which further suggest a temporal deconfinement of the gauge group after inflation. The resulting temperature-dependent axion potential, together with its periodic nature, gives rise to the red-tilted density perturbation, with a positive local-type non-Gaussianity of order unity.

A kinetic mixing between the axion and the inflaton allows for a production of axion dark matter even if the inflationary Hubble scale is smaller than the zero-temperature axion mass. We analyze the axion dynamics in this recently discovered "inflaxion" framework, and present a new cosmological scenario where the axion drifts away from its vacuum during the reheating epoch, giving rise to the observed dark matter abundance. We discuss the implications for both the QCD axion and axion-like particles.

It has been found that the primordial non-Gaussianity of the curvature perturbation in the case of non-Bunch-Davies initial states can be enhanced compared with those in the case of the Bunch-Davies one due to the interactions among the perturbations on subhorizon scales. The purpose of the present paper is to investigate whether tensor non-Gaussianities can also be enhanced or not by the same mechanism. We consider general gravity theory in the presence of an inflaton, and evaluate the tensor auto-bispectrum and the cross-bispectrum involving one tensor and two scalar modes with the non-Bunch-Davies initial states for tensor modes. The crucial difference from the case of the scalar auto-bispectrum is that the tensor three-point function vanishes at the flattened momentum triangles. We point out that the cross-bispectrum can potentially be enhanced at non-trivial triangle shapes due to the non-Bunch-Davies initial states.

We explore the possibility of a blue-tilted gravitational wave spectrum from potential-driven slow-roll inflation in the Horndeski theory. In Kamada et al. (2012), it was claimed that a blue gravitational wave spectrum cannot be obtained from stable potential-driven slow-roll inflation within the Horndeski framework. However, it has been demonstrated that the spectrum of primordial gravitational waves can be blue in inflation with the Gauss-Bonnet term, where the potential term is dominant and slow-roll conditions as well as the stability conditions are satisfied. To fill in this gap, we clarify where the discrepancy is coming from. We extend the formulation of Kamada et al. (2012) and show that a blue gravitational wave spectrum can certainly be generated from stable slow-roll inflation if some of the conditions previously imposed on the form of the free functions in the Lagrangian are relaxed.

In the quiet regions on the solar surface, turbulent convective motions of granulation play an important role in creating small-scale magnetic structures, as well as in energy injection into the upper atmosphere. The turbulent nature of granulation can be studied using spectral line profiles, especially line broadening, which contains information on the flow field smaller than the spatial resolution of an instrument. Moreover, the Doppler velocity gradient along a line-of-sight (LOS) causes line broadening as well. However, the quantitative relationship between velocity gradient and line broadening has not been understood well. In this study, we perform bisector analyses using the spectral profiles obtained using the Spectro-Polarimeter of the Hinode/Solar Optical Telescope to investigate the relationship of line broadening and bisector velocities with the granulation flows. The results indicate that line broadening has a positive correlation with the Doppler velocity gradients along the LOS. We found excessive line broadening in fading granules, that cannot be explained only by the LOS velocity gradient, although the velocity gradient is enhanced in the process of fading. If this excessive line broadening is attributed to small-scale turbulent motions, the averaged turbulent velocity is obtained as 0.9 km/s.

We study relativistic stars in Hordenski theories that evade the gravitational wave constraints and exhibit the Vainshtein mechanism, focusing on a model based on the cubic Galileon Lagrangian. We derive the scalar field profile for static spherically symmetric objects in asymptotically de Sitter space-time with a linear time dependence. The exterior solution matches to the black hole solution found in the literature. Due to the Vainshtein mechanism, the stellar structure is indistinguishable from that of General Relativity with the same central density as long as the radius of the star is shorter than the Vainshtein radius. On the other hand, the scalar field is not suppressed beyond the Vainshtein radius. These solutions have an additional integration constant in addition to the mass of the star.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

We revisit and clarify the gauge dependence of gravitational waves generated at second order from scalar perturbations. In a universe dominated by a perfect fluid with a constant equation-of-state parameter $w$, we compute the energy density of such induced gravitational waves in the Newtonian, comoving, and uniform curvature gauges. Huge differences are found between the Newtonian and comoving gauge results for any $w \,(\ge 0)$. This is always caused by the perturbation of the shift vector. Interestingly, the Newtonian and uniform curvature gauge calculations give the same energy density for $w>0$. In the case of $w=0$, the uniform curvature gauge result differs only by a factor from that of the comoving gauge, but deviates significantly from that of the Newtonian gauge. Our calculation is done analytically for $w=0$ and $w=1/3$, and our result is consistent with the previous numerical one.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

It has been pointed out that matter bounce cosmology driven by a k-essence field cannot satisfy simultaneously the observational bounds on the tensor-to-scalar ratio and non-Gaussianity of the curvature perturbation. In this paper, we show that this is not the case in more general scalar-tensor theories. To do so, we evaluate the power spectra and the bispectra of scalar and tensor perturbations on a general contracting background in the Horndeski theory. We then discuss how one can discriminate contracting models from inflation based on non-Gaussian signatures of tensor perturbations.

A new mechanism for producing axion dark matter is proposed. By invoking low-scale inflation and a kinetic mixing between the axion and the inflaton, it is shown that the axion is driven to a field point slightly displaced from the potential minimum, which can give rise to the observed dark matter abundance. In this framework, different combinations of the axion and inflaton fields play various cosmological roles, including generating the cosmological perturbations, reheating the universe, and serving as dark matter. The kinetic mixing also relates the dark matter lifetime with the reheating temperature. The mechanism tames axions that would otherwise overdominate the universe, and thus opens up new windows in the axion parameter space, including decay constants at the GUT scale and higher.

It is usually assumed that when Weyl invariance is unbroken in the electromagnetic sector, the energy density of primordial magnetic fields will redshift as radiation. Here we show that primordial magnetic fields do not exhibit radiation-like redshifting in the presence of stronger electric fields, as a consequence of Faraday's law of induction. In particular for the standard Maxwell theory, magnetic fields on super-horizon scales can redshift as $B^2 \propto a^{-6} H^{-2}$, instead of the usually assumed $a^{-4}$. Taking into account this effect for inflationary magnetogenesis can correct previous estimates of the magnetic field strength by up to 37 orders of magnitude. This opens new possibilities for inflationary magnetogenesis, and as an example we propose a scenario where femto-Gauss intergalactic magnetic fields are created on Mpc scales, with high-scale inflation producing observable primordial gravitational waves, and reheating happening at low temperatures.

This article is intended to review the recent developments in the Horndeski theory and its generalization, which provide us with a systematic understanding of scalar-tensor theories of gravity as well as a powerful tool to explore astrophysics and cosmology beyond general relativity. This review covers the generalized Galileons, (the rediscovery of) the Horndeski theory, cosmological perturbations in the Horndeski theory, cosmology with a violation of the null energy condition, degenerate higher-order scalar-tensor theories and their status after GW170817, the Vainshtein mechanism in the Horndeski theory and beyond, and hairy black hole solutions.

We consider a subclass of degenerate higher-order scalar-tensor (DHOST) theories in which gravitational waves propagate at the speed of light and do not decay into scalar fluctuations. The screening mechanism in DHOST theories evading these two gravitational wave constraints operates very differently from that in generic DHOST theories. We derive a spherically symmetric solution in the presence of nonrelativistic matter. General relativity is recovered in the vacuum exterior region provided that functions in the Lagrangian satisfy a certain condition, implying that fine-tuning is required. Gravity in the matter interior exhibits novel features: although the gravitational potentials still obey the standard inverse power law, the effective gravitational constant is different from its exterior value, and the two metric potentials do not coincide. We discuss possible observational constraints on this subclass of DHOST theories, and argue that the tightest bound comes from the Hulse-Taylor pulsar.

We investigate the potential of cosmological observations, such as galaxy surveys, for constraining degenerate higher-order scalar-tensor (DHOST) theories, focusing in particular on the linear growth of the matter density fluctuations. We develop a formalism to describe the evolution of the matter density fluctuations during the matter dominated era and in the early stage of the dark energy dominated era in DHOST theories, and give an approximate expression for the gravitational growth index in terms of several parameters characterizing the theory and the background solution under consideration. By employing the current observational constraints on the growth index, we obtain a new constraint on a parameter space of DHOST theories. Combining our result with other constraints obtained from the Newtonian stellar structure, we show that the degeneracy between the effective parameters of DHOST theories can be broken without using the Hulse-Taylor pulsar constraint.

We perform a fully relativistic analysis of even-parity linear perturbations around a static and spherically symmetric solution in the most general scalar-tensor theory with second-order field equations. This paper is a sequel to Kobayashi {\em et al.} (2012), in which the linear perturbation analysis for the odd-parity modes is presented. Expanding the Horndeski action to second order in perturbations and eliminating auxiliary variables, we derive the quadratic action for even-parity perturbations written solely in terms of two dynamical variables. The two perturbations can be interpreted as the gravitational and scalar waves. Correspondingly, we obtain two conditions to evade ghosts and two conditions for the absence of gradient instabilities. Only one in each pair of conditions yields a new stability criterion, as the conditions derived from the stability of the gravitational-wave degree of freedom coincide with those in the odd-parity sector. Similarly, the propagation speed of one of the two modes is the same as that for the odd-parity mode, while the other differs in general from them. Our result is applicable to all the theories of gravitation with an extra single scalar degree of freedom such as the Brans-Dicke theory, $f(R)$ models, and Galileon gravity.

The JVO ALMA WebQL web service - available through the JVO ALMA FITS archive - has been upgraded to include legacy data from other telescopes, for example Nobeyama NRO45M in Japan. The updated server software has been renamed FITSWebQL. In addition, a standalone desktop version supporting Linux, macOS and Windows 10 Linux Subsystem (Bash on Windows) is also available for download from http://jvo.nao.ac.jp/~chris/ . The FITSWebQL server enables viewing of even 100GB-large FITS files in a web browser running on a PC with a limited amount of RAM. Users can interactively zoom-in to selected areas of interest with the corresponding frequency spectrum being calculated on the server in near real-time. The client (a browser) is a JavaScript application built on WebSockets, HTML5, WebGL and SVG. There are many challenges when providing a web browser-based real-time FITS data cube preview service over high-latency low-bandwidth network connections. The upgraded version tries to overcome the latency issue by predicting user mouse movements with a Kalman Filter in order to speculatively deliver the real-time spectrum data at a point where the user is likely to be looking at. The new version also allows one to view multiple FITS files simultaneously in an RGB composite mode (NRO45M FUGIN only), where each dataset is assigned one RGB channel to form a colour image. Spectra from multiple FITS cubes are shown together too. The paper briefly describes main features of FITSWebQL. We also touch on some of the recent developments, such as an experimental switch from C/C++ to Rust (see https://www.rust-lang.org/) for improved stability, better memory management and fearless concurrency, or attempts to display FITS data cubes in the form of interactive on-demand video streams in a web browser.

Among single-field scalar-tensor theories, there is a special class called "cuscuton," which is represented as some limiting case of k-essence in general relativity. This theory has a remarkable feature that the number of propagating degrees of freedom is only two in the unitary gauge in contrast to ordinary scalar-tensor theories with three degrees of freedom. We specify a general class of theories with the same property as the cuscuton in the context of the beyond Horndeski theory, which we dub as the extended cuscuton. We also study cosmological perturbations in the presence of matter in these extended cuscuton theories.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

The quantum effective action for the electromagnetic field in an expanding universe has an anomalous dependence on the scale factor of the metric arising from virtual charged particles in the loops. It has been argued that this Weyl anomaly of quantum electrodynamics sources cosmological magnetic fields in the early universe. We examine this long-standing claim by using the effective action beyond the weak gravitational field limit which has recently been determined. We introduce a general criteria for assessing the quantumness of field fluctuations, and show that the Weyl anomaly is not able to convert vacuum fluctuations of the gauge field into classical fluctuations. We conclude that there is no production of coherent magnetic fields in the universe from the Weyl anomaly of quantum electrodynamics, irrespective of the number of massless charged particles in the theory.

We consider the cosmological dynamics of a scalar field in a potential with multiple troughs and peaks. We show that the dynamics of the scalar field will evolve from light dark matter-like behaviour (such as that of a light axion) to a combination of heavy dark matter-like and dark energy-like behaviour. We discuss the phenomenology of such a model, explaining how it can give rise to the cosmological constant, as well as how it can decouple the dark sector densities between the time of recombination and today, for both the homogeneous background and perturbations. The final form of the dark matter is axion-like, but with abundance and primordial isocurvature modes taking very different values from traditional, axionic, dark matter.

The Horndeski scalar-tensor theory and its recent extensions allow nonlinear derivative interactions of the scalar degree of freedom. We study the matter bispectrum of large scale structure as a probe of these modified gravity theories, focusing in particular on the effect of the terms that newly appear in the so-called "beyond Horndeski" theories. We derive the second-order solution for the matter density perturbations and find that the interactions beyond Horndeski lead to a new time-dependent coefficient in the second-order kernel which differs in general from the standard value of general relativity and the Horndeski theory. This can deform the matter bispectrum at the folded triangle configurations ($k_1+k_2=k_3$), while it is never possible within the Horndeski theory.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We study relativistic stars in degenerate higher-order scalar-tensor theories that evade the constraint on the speed of gravitational waves imposed by GW170817. It is shown that the exterior metric is given by the usual Schwarzschild solution if the lower order Horndeski terms are ignored in the Lagrangian and a shift symmetry is assumed. However, this class of theories exhibits partial breaking of Vainshtein screening in the stellar interior and thus modifies the structure of a star. Employing a simple concrete model, we show that for high-density stars the mass-radius relation is altered significantly even if the parameters are chosen so that only a tiny correction is expected in the Newtonian regime. We also find that, depending on the parameters, there is a maximum central density above which solutions cease to exist.

We propose a new model of single-field D-term inflation in supergravity, where the inflation is driven by a single modulus field which transforms non-linearly under the U(1) gauge symmetry. One of the notable features of our modulus D-term inflation scenario is that the global U(1) remains unbroken in the vacuum and hence our model is not plagued by the cosmic string problem which can exclude most of the conventional D-term inflation models proposed so far due to the CMB observations.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

We investigate constraints on scalar dark matter (DM) by analyzing the Lyman-alpha forest, which probes structure formation at medium and small scales, and also by studying its cosmological consequences at high and low redshift. For scalar DM that constitutes more than 30% of the total DM density, we obtain a lower limit m >~ 10^{-21} eV for the mass of scalar DM. This implies an upper limit on the initial field displacement (or the decay constant for an axion-like field) of phi <~ 10^{16} GeV. We also derive limits on the energy scale of cosmic inflation and establish an upper bound on the tensor-to-scalar ratio of r < 10^{-3} in the presence of scalar DM. Furthermore, we show that there is very little room for ultralight scalar DM to solve the "small-scale crisis" of cold DM without spoiling the Lyman-alpha forest results. The constraints presented in this paper can be used for testing generic theories that contain light scalar fields.

We propose a novel class of degenerate higher-order scalar-tensor theories as an extension of mimetic gravity. By performing a noninvertible conformal transformation on "seed" scalar-tensor theories which may be nondegenerate, we can generate a large class of theories with at most three physical degrees of freedom. We identify a general seed theory for which this is possible. Cosmological perturbations in these extended mimetic theories are also studied. It is shown that either of tensor or scalar perturbations is plagued with gradient instabilities, except for a special case where the scalar perturbations are presumably strongly coupled, or otherwise there appear ghost instabilities.

We perform a fully relativistic analysis of odd-type linear perturbations around a static and spherically symmetric solution in the most general scalar-tensor theory with second-order field equations in four-dimensional spacetime. It is shown that, as in the case of general relativity, the quadratic action for the perturbations reduces to the one having only a single dynamical variable, from which concise formulas for no-ghost and no-gradient instability conditions are derived. Our result is applicable to all the theories of gravity with an extra scalar degree of freedom. We demonstrate how the generic formulas can be applied to some particular examples such as the Brans-Dicke theory, $f(R)$ models, and Galileon gravity.

In string-inspired supergravity theory, K\"{a}hler metric of chiral matter fields often has a pole. Such K\"{a}hler metric is interesting from the viewpoint of the framework of the pole inflation, where the scalar potential can be stretched out to be flat around the pole for a canonically normalized field and inflation can be realized. However, when K\"{a}hler metric has a pole, the scalar potential can also have a pole at the same point in supergravity theory. We study such supergravity models with the pole, and provide numerical analysis of inflationary dynamics and resultant density perturbation. In contrast with usual pole inflation models, inflation in this supergravity based model occurs not on the pole but region apart from the pole. We show that the existence of the pole in the scalar potential is crucial nevertheless. We also examine attractor behavior of our model.

We investigate cosmological constraints on the original relaxion scenario proposed by Graham, Kaplan and Rajendran. We first discuss the appropriate sign choice of the terms in the scalar potential, when the QCD axion is the relaxion with a relaxion-inflaton coupling proposed in the original paper. We next derive the cosmologically consistent ranges of the mass and a coupling of the relaxion for both the QCD relaxion and non-QCD relaxion. The mass range is obtained by $10^{-5}$ eV $\ll m_{\phi} \lesssim 10^4$ eV. We also find that a strong correlation between the Hubble parameter at the relaxion stabilization and the scale $\Lambda$ of non-QCD strong dynamics, which generates the non-perturbative relaxion cosine potential. For a higher relaxion mass, a large scale $\Lambda$ becomes available. However, for its lower mass, $\Lambda$ should be small and constructing such a particle physics model is challenging.

We study the stability of a recently proposed model of scalar-field matter called mimetic dark matter or imperfect dark matter. It has been known that mimetic matter with higher derivative terms suffers from gradient instabilities in scalar perturbations. To seek for an instability-free extension of imperfect dark matter, we develop an effective theory of cosmological perturbations subject to the constraint on the scalar field's kinetic term. This is done by using the unifying framework of general scalar-tensor theories based on the ADM formalism. We demonstrate that it is indeed possible to construct a model of imperfect dark matter which is free from ghost and gradient instabilities. As a side remark, we also show that mimetic $F({\cal R})$ theory is plagued with the Ostrogradsky instability.

We propose a new model of the D-term hybrid inflation in the framework of supergravity. Although our model introduces, analogously to the conventional D-term inflation, the inflaton and a pair of scalar fields charged under a $U(1)$ gauge symmetry, we study the logarithmic and exponential dependence on the inflaton field, respectively, for the K\"ahler and superpotential. This results in a characteristic one-loop scalar potential consisting of linear and exponential terms, which realizes the small-field inflation dominated by the Fayet-Iliopoulos term. With the reasonable values for the coupling coefficients and, in particular, with the $U(1)$ gauge coupling constant comparable to that of the Standard Model, our D-term inflation model can solve the notorious problems in the conventional D-term inflation, namely, the CMB constraints on the spectral index and the generation of cosmic strings.

We investigate the axion inflation model derived by poly-instanton effects in type II superstring theories. Poly-instanton effects are instanton effects corrected by another instanton and it can generate the modulus-axion potential with the double exponential function. Although the axion has a period of small value, this potential can have a flat region because its derivatives are exponentially suppressed by non-perturbative effects. From the view point of the cosmic inflation, such potential is interesting. In this paper, we numerically study the possibilities for realizing the cosmic inflation. We also study their spectral index and other cosmological observables, numerically.

We propose a new type of axion inflation with complex structure moduli in the framework of type IIB superstring theory compactified on Calabi-Yau manifold. The inflaton is identified as the axion for the complex structure moduli whose potential is originating from instantonic corrections appearing through the period vector of mirror Calabi-Yau manifold. The axionic shift symmetry is broken down to the discrete one by the inclusion of instantonic correction and certain three-from fluxes. Our proposed inflation scenario is compatible with K\"ahler moduli stabilization. We also study a typical reheating temperature in the case of complex structure moduli inflation.

An invertible field transformation is such that the old field variables correspond one-to-one to the new variables. As such, one may think that two systems that are related by an invertible transformation are physically equivalent. However, if the transformation depends on field derivatives, the equivalence between the two systems is nontrivial due to the appearance of higher derivative terms in the equations of motion. To address this problem, we prove the following theorem on the relation between an invertible transformation and Euler-Lagrange equations: If the field transformation is invertible, then any solution of the original set of Euler-Lagrange equations is mapped to a solution of the new set of Euler-Lagrange equations, and vice versa. We also present applications of the theorem to scalar-tensor theories.

String theory generically predicts the coupling between the Affleck-Dine field and axion field through higher-dimensional operators. We thus explore the Affleck-Dine baryogenesis on an axion background. It turns out that the axion oscillation produces an enough amount of baryon asymmetry of the Universe just after the inflation, even without a soft supersymmetry-breaking $A$-term. This baryogenesis scenario is applicable to the string axion inflation.

We study general class of small-field axion inflations which are the mixture of polynomial and sinusoidal functions suggested by the natural and axion monodromy inflations. The axion decay constants leading to the successful axion inflations are severely constrained in order not to spoil the Big-Bang nucleosynthesis and overproduce the isocurvature perturbation originating from the QCD axion. We in turn find that the cosmologically favorable axion decay constants are typically of order the grand unification scale or the string scale which is consistent with the prediction of closed string axions.

It has been pointed out that non-singular cosmological solutions in second-order scalar-tensor theories generically suffer from gradient instabilities. We extend this no-go result to second-order gravitational theories with an arbitrary number of interacting scalar fields. Our proof follows directly from the action of generalized multi-Galileons, and thus is different from and complementary to that based on the effective field theory approach. Several new terms for generalized multi-Galileons on a flat background were proposed recently. We find a covariant completion of them and confirm that they do not participate in the no-go argument.

We present a new scenario for generating the baryon asymmetry of the universe that is induced by a Nambu-Goldstone (NG) boson. The shift symmetry naturally controls the operators in the theory, while allowing the NG boson to couple to the spacetime geometry as well as to the baryons. The cosmological background thus sources a coherent motion of the NG boson, which leads to baryogenesis. Good candidates of the baryon-generating NG boson are the QCD axion and axion-like fields. In these cases the axion induces baryogenesis in the early universe, and can also serve as dark matter in the late universe.

We propose a novel branch of the Galilean Genesis scenario as an alternative to inflation, in which the universe starts expanding from Minkowski in the asymptotic past with a gross violation of the null energy condition (NEC). This variant, described by several functions and parameters within the Horndeski scalar-tensor theory, shares the same background dynamics with the existing Genesis models, but the nature of primordial quantum fluctuations is quite distinct. In some cases, tensor perturbations grow on superhorizon scales. The tensor power spectrum can be red, blue, or scale invariant, depending on the model, while scalar perturbations are nearly scale invariant. This is in sharp contrast to typical NEC-violating cosmologies, in which a blue tensor tilt is generated. Though the primordial tensor and scalar spectra are both nearly scale invariant as in the inflationary scenario, the consistency relation in our variant of Galilean Genesis is non-standard.

If there are a plethora of axions in nature, they may have a complicated potential and create an axion landscape. We study a possibility that one of the axions is so light that it is cosmologically stable, explaining the observed dark matter density. In particular we focus on a case in which two (or more) shift-symmetry breaking terms conspire to make the axion sufficiently light at the potential minimum. In this case the axion has a flat-bottomed potential. In contrast to the case in which a single cosine term dominates the potential, the axion abundance as well as its isocurvature perturbations are significantly suppressed. This allows an axion with a rather large mass to serve as dark matter without fine-tuning of the initial misalignment, and further makes higher-scale inflation to be consistent with the scenario.

We propose a new class of spontaneous baryogenesis models that does not produce baryon isocurvature perturbations. The baryon chemical potential in these models is independent of the field value of the baryon-generating scalar, hence the scalar field fluctuations are blocked from propagating into the baryon isocurvature. We demonstrate this mechanism in simple examples where spontaneous baryogenesis is driven by a non-canonical scalar field. The suppression of the baryon isocurvature allows spontaneous baryogenesis to be compatible even with high-scale inflation.

The conventional slow-roll approximation is broken in the so-called "ultra slow-roll" models of inflation, for which the inflaton potential is exactly (or extremely) flat. The interesting nature of (canonical) ultra slow-roll inflation is that the curvature perturbation grows on superhorizon scales, but has a scale-invariant power spectrum. We study the ultra slow-roll inflationary dynamics in the presence of non-canonical kinetic terms of the scalar field, namely ultra slow-roll G-inflation. We compute the evolution of the curvature perturbation and show that the primordial power spectrum follows a broken power law with an oscillation feature. It is demonstrated that this could explain the lack of large-scale power in the cosmic microwave background temperature anisotropies. We also point out that the violation of the null energy condition is prohibited in ultra slow-roll G-inflation and hence a blue tensor tilt is impossible as long as inflation is driven by the potential. This statement is, however, not true if the energy density is dominated by the kinetic energy of the scalar field.

Galilean genesis is an alternative to inflation, in which the universe starts expanding from Minkowski with the stable violation of the null energy condition. In this paper, we discuss how the early universe is reheated through the gravitational particle production at the trasition from the genesis phase to the subsequent phase where the kinetic energy of the scalar field is dominant. We then study the consequences of gravitational reheating after Galilean genesis on the spectrum of primordial gravitational waves. The resultant spectrum is strongly blue, and at high frequencies $\Omega_{\rm gw}\propto f^3$ in terms of the energy density per unit logarithmic frequency. Though this cannot be detected in existing detctors, the amplitude can be as large as $\Omega_{\rm gw}\sim 10^{-12}$ at $f\sim 100\,$MHz, providing a future test of the genesis scenario. The analysis is performed within the framework of generalized Galilean genesis based on the Horndeski theory, which enables us to derive generic formulas.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

A QCD axion with a time-dependent decay constant has been known to be able to accommodate high-scale inflation without producing topological defects or too large isocurvature perturbations on CMB scales. We point out that a dynamical decay constant also has the effect of enhancing the small-scale axion isocurvature perturbations. The enhanced axion perturbations can even exceed the periodicity of the axion potential, and thus lead to the formation of axionic domain walls. Unlike the well-studied axionic walls, the walls produced from the enhanced perturbations are not bounded by cosmic strings, and thus would overclose the universe independently of the number of degenerate vacua along the axion potential.

We investigate cosmological aspects of spontaneous baryogenesis driven by a scalar field, and present general constraints that are independent of the particle physics model. The relevant constraints are obtained by studying the backreaction of the produced baryons on the scalar field, the cosmological expansion history after baryogenesis, and the baryon isocurvature perturbations. We show that cosmological considerations alone provide powerful constraints, especially for the minimal scenario with a quadratic scalar potential. Intriguingly, we find that for a given inflation scale, the other parameters including the reheat temperature, decoupling temperature of the baryon violating interactions, and the mass and decay constant of the scalar are restricted to lie within ranges of at most a few orders of magnitude. We also discuss possible extensions to the minimal setup, and propose two ideas for evading constraints on isocurvature perturbations: one is to suppress the baryon isocurvature with nonquadratic scalar potentials, another is to compensate the baryon isocurvature with cold dark matter isocurvature by making the scalar survive until the present.

The null energy condition can be violated stably in generalized Galileon theories, which gives rise to the possibilities of healthy non-singular cosmologies. However, it has been reported that in many cases cosmological solutions are plagued with instabilities or have some pathologies somewhere in the whole history of the universe. Recently, this was shown to be generically true in a certain subclass of the Horndeski theory. In this short paper, we extend this no-go argument to the full Horndeski theory, and show that non-singular models (with flat spatial sections) in general suffer either from gradient instabilities or some kind of pathology in the tensor sector. This implies that one must go beyond the Horndeski theory to implement healthy non-singular cosmologies.

We study an axion inflation model recently proposed within the framework of type IIB superstring theory, where we pay a particular attention to a sub-Planckian axion decay constant. Our axion potential can lead to the small field inflation with a small tensor-to-scalar ratio, and a typical reheating temperature can be as low as GeV.

We present a real-time supernova neutrino burst monitor at Super-Kamiokande (SK). Detecting supernova explosions by neutrinos in real time is crucial for giving a clear picture of the explosion mechanism. Since the neutrinos are expected to come earlier than light, a fast broadcasting of the detection may give astronomers a chance to make electromagnetic radiation observations of the explosions right at the onset. The role of the monitor includes a fast announcement of the neutrino burst detection to the world and a determination of the supernova direction. We present the online neutrino burst detection system and studies of the direction determination accuracy based on simulations at SK.

Modular invariance is a striking symmetry in string theory, which may keep stringy corrections under control. In this paper, we investigate a phenomenological consequence of the modular invariance, assuming that this symmetry is preserved as well as in a four dimensional (4D) low energy effective field theory. As a concrete setup, we consider a modulus field $T$ whose contribution in the 4D effective field theory remains invariant under the modular transformation and study inflation drived by $T$. The modular invariance restricts a possible form of the scalar potenntial. As a result, large field models of inflation are hardly realized. Meanwhile, a small field model of inflation can be still accomodated in this restricted setup. The scalar potential traced during the slow-roll inflation mimics the hilltop potential $V_{ht}$, but it also has a non-negligible deviation from $V_{ht}$. Detecting the primordial gravitational waves predicted in this model is rather challenging. Yet, we argue that it may be still possible to falsify this model by combining the information in the reheating process which can be determined self-completely in this setup.

We study the moduli-induced gravitino problem within the framework of the phenomenologically attractive mirage mediations. The huge amount of gravitino generated by the moduli decay can be successfully diluted by introducing an extra light modulus field which does not induce the supersymmetry breaking. Since the lifetime of extra modulus field becomes longer than usually considered modulus field, our proposed mechanism is applied to both the low- and high-scale supersymmetry breaking scenarios. We also point out that such an extra modulus field appears in the flux compactification of type II string theory.

We present generic bounds on magnetic fields produced from cosmic inflation. By investigating field bounds on the vector potential, we constrain both the quantum mechanical production of magnetic fields and their classical growth in a model independent way. For classical growth, we show that only if the reheating temperature is as low as T_{reh} <~ 10^2 MeV can magnetic fields of 10^{-15} G be produced on Mpc scales in the present universe. For purely quantum mechanical scenarios, even stronger constraints are derived. Our bounds on classical and quantum mechanical scenarios apply to generic theories of inflationary magnetogenesis with a two-derivative time kinetic term for the vector potential. In both cases, the magnetic field strength is limited by the gravitational back-reaction of the electric fields that are produced simultaneously. As an example of quantum mechanical scenarios, we construct vector field theories whose time diffeomorphisms are spontaneously broken, and explore magnetic field generation in theories with a variable speed of light. Transitions of quantum vector field fluctuations into classical fluctuations are also analyzed in the examples.

We clarify the features of primordial non-Gaussianities of tensor perturbations in Gao's unifying framework of scalar-tensor theories. The general Lagrangian is given in terms of the ADM variables so that the framework maintains spatial covariance and includes the Horndeski theory and Gleyzes-Langlois-Piazza-Vernizzi (GLPV) generalization as specific cases. It is shown that the GLPV generalization does not give rise to any new terms in the cubic action compared to the case of the Horndeski theory, but four new terms appear in more general theories beyond GLPV. We compute the tensor 3-point correlation functions analytically by treating the modification to the dispersion relation as a perturbation. The relative change in the 3-point functions due to the modified dispersion relation is only mildly configuration-dependent. When the effect of the modified dispersion relation is small, there is only a single cubic term generating squeezed non-Gaussianity, which is the only term present in general relativity. The corresponding non-Gaussian amplitude has a fixed and universal feature, and hence offers a "consistency relation" for primordial tensor modes in a quite wide class of single-field inflation models. All the other cubic interactions are found to give peaks at equilateral shapes.

The JVO ALMA Archive provides users one of the easiest ways to access the ALMA archival data. The users can have a quick look at a 3 or 4-dimensional data cube without downloading multiple huge tarballs from a science portal of ALMA Regional Centers (ARCs). Since we just synchronize all datasets with those of ARCs, the metadata are identical to the upstream, including ``target name'' for each dataset. The name is not necessarily a common one like NGC numbers, but sometimes one of sequential numbers assigned in an observation proposal. Compilation of these artificial names into astronomical ones could provide users more flexible and powerful search interfaces; for instance, with the knowledge of the redshift for each source, the users can easily find the datasets which observed their interested emission/absorption lines at not the observer frame but the rest frame, fitting well with theoretical studies. To implement this functionality, cross-identification of all the sources in our archive with those in some other astronomical databases such as NED and SIMBAD is required. We developed a tiny Java application named ``Blade Runner'' for this purpose. The program works as a crawler for both the JVO ALMA Archive and SIMBAD, storing all information onto a SQLite-based database file; this portable design enables us to communicate results to each other even under different computing environments. In this paper, we introduce its software design and our recent work on the application, and report a preliminary result on the source identification in our archive.

We address the question of how one can modify the inflationary tensor spectrum without changing at all the successful predictions on the curvature perturbation. We show that this is indeed possible, and determine the two quadratic curvature corrections that are free from instabilities and affect only the tensor sector at the level of linear cosmological perturbations. Both of the two corrections can reduce the tensor amplitude, though one of them generates large non-Gaussianity of the curvature perturbation. It turns out that the other one corresponds to so-called Lorentz-violating Weyl gravity. In this latter case one can obtain as small as 65% of the standard tensor amplitude. Utilizing this effect we demonstrate that even power-law inflation can be within the 2$\sigma$ contour of the Planck results.

We consider domain walls in the $Z_3$ symmetric NMSSM. The spontaneous $Z_3$ discrete symmetry breaking produces domain walls, and the stable domain walls are problematic. Thus, we assume the $Z_3$ symmetry is slightly but explicitly broken and the domain walls decay. Such a decay causes a large late-time entropy production. We study its cosmological implications on unwanted relics such as moduli, gravitino, LSP and axion.

We carry out a detailed study of the orbital dynamics and structural evolution of over 6000 subhalos in the Via Lactea II simulation, from infall to present. By analyzing subhalos with masses down to m = 4e5 Msun, we find that lower mass subhalos, which are not strongly affected by dynamical friction, exhibit behaviors qualitatively different from those found previously for more massive ones. Furthermore, there is a clear trend of subhalos that fell into the host earlier being less concentrated. We show that the concentration at infall characterizes various aspects of subhalo evolution. In particular, tidal effects truncate the growth of less concentrated subhalos at larger distances from the host; subhalos with smaller concentrations have larger infall radii. The concentration at infall is further shown to be a determining factor for the subsequent mass loss of subhalos within the host, and also for the evolution of their internal structure in the v_max-r_max plane. Our findings raise the prospects of using the concentration to predict the tidal evolution of subhalos, which will be useful for obtaining analytic models of galaxy formation, as well as for near field cosmology.

We study the inflation scenarios, in the framework of superstring theory, where the inflaton is an axion producing the adiabatic curvature perturbations while there exists another light axion producing the isocurvature perturbations. We discuss how the non-trivial couplings among string axions can generically arise, and calculate the consequent cross-correlations between the adiabatic and isocurvature modes through concrete examples. Based on the Planck analysis on the generally correlated isocurvature perturbations, we show that there is a preference for the existence of the correlated isocurvature modes for the axion monodromy inflation while the natural inflation disfavors such isocurvature modes.

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.

Axions in the Peccei-Quinn (PQ) mechanism provide a promising solution to the strong CP problem in the standard model of particle physics. Coherently generated PQ scalar fields could dominate the energy density in the early Universe and decay into relativistic axions, which would conflict with the current dark radiation constraints. We study the possibility that a thermal inflation driven by a $U(1)$ gauged Higgs field dilutes such axions. A well-motivated extra gauged $U(1)$ would be the local $B-L$ symmetry. We also discuss the implication for the case of $U(1)_{B-L}$ and an available baryogenesis mechanism in such cosmology.

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.

Recently it was argued that gravity with the squire of the Ricci tensor can be stabilized by adding constraints to the theory. This was so far demonstrated for fluctuations on the Minkowski/de Sitter background. We show that the same scheme works equally well for removing Ostrogradski's ghost from fluctuations on a cosmological background in generic $f(R,R_{\mu\nu}^2,C_{\mu\nu\rho\sigma}^2)$-type theories of gravity. We also derive the general formula for the spectrum of primordial tensor perturbations from the stabilized theory.

Recent detection of $B$-mode polarization induced from tensor perturbations by the BICEP2 experiment implies so-called large field inflation, where an inflaton field takes super-Planckian expectation value during inflation, at a high energy scale. We show however, if another inflation follows hybrid inflation, the hybrid inflation can generate a large tensor perturbation with not super-Plankian but Planckian field value. This scenario would relax the tension between BICEP2 and Planck concerning the tensor-to-scalar ratio, because a negative large running can also be obtained for a certain number of e-fold of the hybrid inflation. A natural interpretation of a large gravitational wave mode with or without the scalar spectral running might be multiple inflation in the early Universe.

We propose a mechanism for the natural inflation with and without modulation in the framework of type IIB string theory on toroidal orientifold or orbifold. We explicitly construct the stabilization potential of complex structure, dilaton and K\"ahler moduli, where one of the imaginary component of complex structure moduli becomes light which is identified as the inflaton. The inflaton potential is generated by the gaugino-condensation term which receives the one-loop threshold corrections determined by the field value of complex structure moduli and the axion decay constant of inflaton is enhanced by the inverse of one-loop factor. We also find the threshold corrections can also induce the modulations to the original scalar potential for the natural inflation. Depending on these modulations, we can predict several sizes of tensor-to-scalar ratio as well as the other cosmological observables reported by WMAP, Planck and/or BICEP2 collaborations.

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.

We provide a unified description of the hemispherical asymmetry in the cosmic microwave background generated by the mechanism proposed by Erickcek, Kamionkowski, and Carroll, using a delta N formalism that consistently accounts for the asymmetry-generating mode throughout. We derive a general form for the power spectrum which explicitly exhibits the broken translational invariance. This can be directly compared to cosmic microwave background observables, including the observed quadrupole and fNL values, automatically incorporating the Grishchuk--Zel'dovich effect. Our calculation unifies and extends previous calculations in the literature, in particular giving the full dependence of observables on the phase of our location in the super-horizon mode that generates the asymmetry. We demonstrate how the apparently different results obtained by previous authors arise as different limiting cases. We confirm the existence of non-linear contributions to the microwave background quadrupole from the super-horizon mode identified by Erickcek et al. and further explored by Kanno et al., and show that those contributions are always significant in parameter regimes capable of explaining the observed asymmetry. We indicate example parameter values capable of explaining the observed power asymmetry without violating other observational bounds.

The galilean genesis scenario is an alternative to inflation in which the universe starts expanding from Minkowski in the asymptotic past by violating the null energy condition stably. Several concrete models of galilean genesis have been constructed so far within the context of galileon-type scalar-field theories. We give a generic, unified description of the galilean genesis scenario in terms of the Horndeski theory, i.e., the most general scalar-tensor theory with second-order field equations. In doing so we generalize the previous models to have a new parameter (denoted by {\alpha}) which results in controlling the evolution of the Hubble rate. The background dynamics is investigated to show that the generalized galilean genesis solution is an attractor, similarly to the original model. We also study the nature of primordial perturbations in the generalized galilean genesis scenario. In all the models described by our generalized genesis Lagrangian, amplification of tensor perturbations does not occur as opposed to what happens in quasi-de Sitter inflation. We show that the spectral index of curvature perturbations is determined solely from the parameter {\alpha} and does not depend on the other details of the model. In contrast to the original model, a nearly scale-invariant spectrum of curvature perturbations is obtained for a specific choice of {\alpha}.

The Horndeski theory of gravity is known as the most general scalar-tensor theory with second-order field equations. Recently, it was demonstrated by Gleyzes et al. that the Horndeski theory can further be generalized in such a way that although field equations are of third order, the number of propagating degrees of freedom remains the same. We study small-scale gravity in the generalized Horndeski theory, focusing in particular on an impact of the new derivative interaction beyond Horndeski on the Vainshtein screening mechanism. In the absence of the quintic galileon term and its generalization, we show that the new interaction does not change the qualitative behavior of gravity outside and near the source: the two metric potentials coincide, $\Phi = \Psi \;(\sim r^{-1})$, while the gravitational coupling is given by the cosmological one and hence is time-dependent in general. We find, however, that the gravitational field inside the source shows a novel behavior due to the interaction beyond Horndeski: the gravitational attraction is not determined solely from the enclosed mass and two potentials do not coincide, indicating breaking of the screening mechanism.

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.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

We study parameter space in the axion monodromy inflation corrected by dynamically generated terms involving with the axion. The potential has the linear potential with multiple sinusoidal functions, which play a role in generating modulations. We show that this potential leads both to a large tensor-to-scalar ratio $r_T \sim 0.16$ and to a large negative running of spectral index $\alpha_s \sim - (0.02 -0.03)$. To realize these results, a small hierarchy among dynamical scales is required whereas the decay constants in sinusoidal functions remain sub-Planckian in this model. We discuss also reheating process after the inflation in a bottom-up approach.

We study the axion monodromy inflation with a non-perturbatively generated sinusoidal term. The potential form is a mixture between the natural inflation and the axion monodromy inflation potentials. The sinusoidal term is subdominant in the potential, but leaves significant effects on the resultant fluctuation generated during inflation. A larger tensor-to-scalar ratio can be obtained in our model. We study two scenarios, single inflation scenario and the double inflation scenario. In the first scenario, the axion monodromy inflation with a sufficient number of e-fold generates a larger tensor-to-scalar ratio about $0.1 - 0.15$ but also a tiny running of spectral index. In the second scenario of double inflation, axion monodromy inflation is its first stage and, we assume another inflation follows. In this case, our model can realize a larger tensor-to-scalar ratio and a large negative running of spectral index simultaneously.

We investigate pair creation by an electric field in four-dimensional de Sitter space. The expectation value of the induced current is computed, using the method of adiabatic regularization. Under strong electric fields the behavior of the current is similar to that in flat space, while under weak electric fields the current becomes inversely proportional to the mass squared of the charged field. Thus we find that the de Sitter space obtains a large conductivity under weak electric fields in the presence of a charged field with a tiny mass. We then apply the results to constrain electromagnetic fields in the early universe. In particular, we study cosmological scenarios for generating large-scale magnetic fields during the inflationary era. Electric fields generated along with the magnetic fields can induce sufficiently large conductivity to terminate the phase of magnetogenesis. For inflationary magnetogenesis models with a modified Maxwell kinetic term, the generated magnetic fields cannot exceed 10^{-30} G on Mpc scales in the present epoch, when a charged field carrying an elementary charge with mass of order the Hubble scale or smaller exists in the Lagrangian. Similar constraints from the Schwinger effect apply for other magnetogenesis mechanisms.

We propose the natural inflation from the heterotic string theory on "Swiss-Cheese" Calabi-Yau manifold with multiple $U(1)$ magnetic fluxes. Such multiple $U(1)$ magnetic fluxes stabilize the same number of the linear combination of the universal axion and K\"ahler axions and one of the K\"ahler axions is identified as the inflaton. This axion decay constant can be determined by the size of one-loop corrections to the gauge kinetic function of the hidden gauge groups, which leads effectively to the trans-Planckian axion decay constant consistent with the WMAP, Planck and/or BICEP2 data. During the inflation, the real parts of the moduli are also stabilized by employing the nature of the "Swiss-Cheese" Calabi-Yau manifold.

We revisit large field inflation models with modulations in light of the recent discovery of the primordial B-mode polarization by the BICEP2 experiment, which, when combined with the Planck + WP + highL data, gives a strong hint for additional suppression of the CMB temperature fluctuations at small scales. Such a suppression can be explained by a running spectral index. In fact, it was pointed out by two of the present authors (TK and FT) that the existence of both tensor mode perturbations and a sizable running of the spectral index is a natural outcome of large inflation models with modulations. We find that this holds also in the recently proposed multi-natural inflation, in which the inflaton potential consists of multiple sinusoidal functions and therefore the modulations are a built-in feature.

We explore cosmological magnetogenesis in the post-inflationary universe, when the inflaton oscillates around its potential minimum and the universe is effectively dominated by cold matter. During this epoch prior to reheating, large-scale magnetic fields can be significantly produced by the cosmological background. By considering magnetogenesis both during and after inflation, we demonstrate that magnetic fields stronger than 10^{-15} G can be generated on Mpc scales without having strong couplings in the theory, or producing too large electric fields that would dominate the universe.

Large field inflation models are favored by the recent BICEP2 that has detected gravitational wave modes generated during inflation. We study general large field inflation models for which the potential contains (constant) quadratic and quartic terms of inflaton field. We show, in this framework, those inflation models can generate the fluctuation with the tensor-to-scalar ratio of $0.2$ as well as the scalar spectral index of $0.96$: those are very close to the center value of the tensor-to-scalar ratio reported by BICEP2 as well as Planck. Finally, we briefly discuss the particle physics model building of inflation.

The breaking of $U(1)_R$ symmetry plays a crucial role in modeling the breaking of supersymmetry (SUSY). In the models that possess both SUSY preserving and SUSY breaking vacua, tube-like cosmic strings called R-tubes, whose surfaces are constituted by domain walls interpolating a false and a true vacuum with some winding numbers, can exist. Their (in)stability can strongly constrain SUSY breaking models theirselves. In the present study, we investigate the dynamical (in)stability of two colliding metastable tube-like strings by field-theoretic simulations. From them, we find that the strings become unstable, depending on the relative collision angle and speed of two strings, and the false vacuum is eventually filled out by the true vacuum owing to rapid expansion of the strings or unstable bubbles created as remnants of the collision.

In the modulated reheating scenario, the field that drives inflation has a spatially varying decay rate, and the resulting inhomogeneous reheating process generates adiabatic perturbations. We examine the statistical properties of the density perturbations generated in this scenario. Unlike earlier analyses, we include the dynamics of the field that determines the inflaton decay rate. We show that the dynamics of this modulus field can significantly alter the amplitude of the power spectrum and the bispectrum, even if the modulus field has a simple potential and its effective mass is smaller than the Hubble rate. In some cases, the evolution of the modulus amplifies the non-Gaussianity of the perturbations to levels that are excluded by recent observations of the cosmic microwave background. Therefore, a proper treatment of the modulus dynamics is required to accurately calculate the statistical properties of the perturbations generated by modulated reheating.

We study the spontaneous R-symmetry breaking model and investigate the cosmological constraints on this model due to the pseudo Nambu-Goldstone boson, R-axion. We consider the R-axion which has relatively heavy mass in order to complement our previous work. In this regime, model parameters, R-axions mass and R-symmetry breaking scale, are constrained by Big Bang Nucleosynthesis and overproduction of the gravitino produced from R-axion decay and thermal plasma. We find that the allowed parameter space is very small for high reheating temperature. For low reheating temperature, the U(1)_R breaking scale f_a is constrained as f_a<10^{12-14}GeV regardless of the value of R-axion mass.

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.

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

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 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 revisit the isocurvature density perturbations induced by quantum fluctuations of the axion field by extending a recently developed analytic method and approximations to a time-dependent scalar potential, which enables us to follow the evolution of the axion until it starts to oscillate. We find that, as the initial misalignment angle approaches the hilltop of the potential, the isocurvature perturbations become significantly enhanced, while the non-Gaussianity parameter increases slowly but surely. As a result, the isocurvature constraint on the inflation scale is tightened as Hinf < O(100) GeV for the axion decay constant fa < 10^10 GeV, near the smaller end of the axion dark matter window. We also derive useful formulae for the power spectrum and non-Gaussianity of the isocurvature perturbations.

Cosmic inflation driven by branes wrapping the extra dimensions involves Kaluza-Klein (KK) degrees of freedom in addition to the zero-mode position of the brane which plays the role of the inflaton. As the wrapped brane passes by localized sources or features along its inflationary trajectory in the extra dimensional space, the KK modes along the wrapped direction are excited and start to oscillate during inflation. We show that the oscillating KK modes induce parametric resonance for the curvature perturbations, generating sharp signals in the perturbation spectrum. The effective four dimensional picture is a theory where the inflaton couples to the heavy KK modes. The Nambu-Goto action of the brane sources couplings between the inflaton kinetic terms and the KK modes, which trigger significant resonant amplification of the curvature perturbations. We find that the strong resonant effects are localized to narrow wave number ranges, producing spikes in the perturbation spectrum. Investigation of such resonant signals opens up the possibility of probing the extra dimensional space through cosmological observations.

We show that density perturbations seeded by the inflaton can be suppressed when having additional light degrees of freedom contributing to the production of perturbations. The inflaton fluctuations affect the light field dynamics by modulating the length of the inflationary period, hence produce additional density perturbations in the post-inflationary era. Such perturbations can cancel those generated during inflation as both originate from the same inflaton fluctuations. This allows production of large gravitational waves from small-field inflation, which is normally forbidden by the Lyth bound on the inflaton field excursion. We also find that the field bound is taken over by the light scalar when the inflaton-induced perturbations are suppressed, thus present a generalized form of the Lyth bound that applies to the total field space.

We explore the possibility of testing modified gravity exhibiting the Vainshtein mechanism against observations of cluster lensing. We work in the most general scalar-tensor theory with second-order field equations (Horndeski's theory), and derive static and spherically symmetric solutions, for which the scalar field is screened below a certain radius. It is found that the essential structure of the problem in the most general case can be captured by the program of classifying Vainshtein solutions out of different solutions to a quintic equation, as has been performed in the context of massive gravity. The key effect on gravitational lensing is that the second derivative of the scalar field can substantially be large at the transition from screened to unscreened regions, leaving a dip in the convergence. This allows us to put observational constraints on parameters characterizing the general scalar-tensor modification of gravity. We demonstrate how this occurs in massive gravity as an example, and discuss its observational signatures in cluster lensing.

We develop a theory of nonlinear cosmological perturbations on superhorizon scales for generic single-field inflation. Our inflaton is described by the Lagrangian of the form $W(X,\phi)-G(X,\phi)\Box\phi$ with $X=-\partial^{\mu}\phi\partial_{\mu}\phi/2$, which is no longer equivalent to a perfect fluid. This model is more general than k-inflation, and is called G-inflation. A general nonlinear solution for the metric and the scalar field is obtained at second order in gradient expansion. We derive a simple master equation governing the large-scale evolution of the nonlinear curvature perturbation. It turns out that the nonlinear evolution equation is deduced as a straightforward extension of the corresponding linear equation for the curvature perturbation on uniform $\phi$ hypersurfaces.

We study general constraints on spontaneous R-symmetry breaking models coming from the cosmological effects of the pseudo Nambu-Goldstone bosons, R-axions. They are substantially produced in the early Universe and may cause several cosmological problems. We focus on relatively long-lived R-axions and find that in a wide range of parameter space, models are severely constrained. In particular, R-axions with mass less than 1 MeV are generally ruled out for relatively high reheating temperature, $T_R>10$ GeV.

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

We have observed the primary electron spectrum from 30 GeV to 3 TeV using emulsion chambers flown by balloons at the top of the atmosphere, for the purpose of exploring the origin of cosmic rays in the Galaxy. The atmospheric gamma rays have been simultaneously observed in the 30 GeV $\sim$ 8 TeV energy range. In this paper, we estimate the atmospheric electron spectrum in the upper atmosphere ($<$ 10 ${\rm g/cm^2}$) from our observed gamma-ray spectrum using the electromagnetic shower theory in order to derive the primary cosmic-ray electron spectrum. The transport equations of the electron and gamma-ray spectrum are analytically solved and the results are compared with those of Monte Carlo simulation (MC). Since we used the observed atmospheric gamma rays as the source of atmospheric electrons, our solutions are free from ambiguities on the primary cosmic-ray nuclear spectra and nuclear interaction models included in MC. In the energy range above several hundred GeV, the Dalitz electrons produced directly from neutral pions contribute around 10 percent to the whole atmospheric electron spectrum at the depth of 4 ${\rm g/cm^2}$, which increases in importance at the higher altitude and cannot be ignored in TeV electron balloon observations.

We have performed a series of cosmic-ray electron observations using the balloon-borne emulsion chambers since 1968. While we previously reported the results from subsets of the exposures, the final results of the total exposures up to 2001 are presented here. Our successive experiments have yielded the total exposure of 8.19 m^2 sr day at the altitudes of 4.0 - 9.4 g/cm^2. The performance of the emulsion chambers was examined by accelerator beam tests and Monte-Carlo simulations, and the on-board calibrations were carried out by using the flight data. In this work we present the cosmic-ray electron spectrum in the energy range from 30 GeV to 3 TeV at the top of the atmosphere, which is well represented by a power-law function with an index of -3.28+-0.10. The observed data can be also interpreted in terms of diffusive propagation models. The evidence of cosmic-ray electrons up to 3 TeV suggests the existence of cosmic-ray electron sources at distances within ~1 kpc and times within ~1x10^5 yr ago.

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

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.

We investigate the scale-dependence, or the runnings, of linear and second order density perturbations generated in various curvaton scenarios. We argue that the second order perturbations, i.e. non-Gaussianity, can strongly depend on the scale, even when the linear perturbations are nearly scale-invariant. We present analytic formulae for the runnings from curvatons with general energy potentials, and clarify the conditions under which fNL becomes strongly scale-dependent. From the point of view of the fNL running, curvaton potentials can be classified into roughly two categories by whether the potential flattens or steepens compared to a quadratic one. As such examples, we study pseudo-Nambu-Goldstone curvatons, and self-interacting curvatons, respectively. The dynamics of non-quadratic curvatons and the behaviors of the resulting density perturbations are clarified by analytical methods. Then we also study models where multiple source can be responsible for density perturbations such as the multi-curvaton, and mixed curvaton and inflaton models where the running of fNL can also be large due to their multi-source nature. We make quantitative analysis for each curvaton scenario and discuss in what cases the scale-dependence, in particular, of fNL can be large enough to be probed with future CMB experiments.

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.

The cross-correlation between the integrated Sachs-Wolfe (ISW) effect and the large scale structure (LSS) is a powerful tool to constrain dark energy and alternative theories of gravity. In this paper, we obtain observational constraints on kinetic gravity braiding from the ISW-LSS cross-correlation. We find that the late-time ISW effect in the kinetic gravity braiding model anti-correlates with large scale structures in a wide range of parameters, which clearly demonstrates how one can distinguish modified gravity theories from the LCDM model using the ISW effect. In addition to the analysis based on a concrete model, we investigate a future prospect of the ISW-LSS cross-correlation by using a phenomenological parameterization of modified gravity models.

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.

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce \nu_{e}, \nu_{\mu} and \bar{\nu}_{\mu}. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(\sigma_0/1 mb) < 6.3 \times 10^{-24} (\beta_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where \beta_M is the monopole velocity in units of the speed of light and \sigma_0 is the catalysis cross section at \beta_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for \beta_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

A generic second-order scalar-tensor theory contains a nonlinear derivative self-interaction of the scalar degree of freedom $\phi$ \`{a} la Galileon models, which allows for the Vainshtein screening mechanism. We investigate this effect on subhorizon scales in a cosmological background, based on the most general second-order scalar-tensor theory. Our analysis takes into account all the relevant nonlinear terms and the effect of metric perturbations consistently. We derive an explicit form of Newton's constant, which in general is time-dependent and hence is constrained from observations, as suggested earlier. It is argued that in the most general case the inverse-square law cannot be reproduced on the smallest scales. Some applications of our results are also presented.

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

In the Horndeski's most general scalar-tensor theories the equations of scalar density perturbations are derived in the presence of non-relativistic matter minimally coupled to gravity. Under a quasi-static approximation on sub-horizon scales we obtain the effective gravitational coupling $G_{eff}$ associated with the growth rate of matter perturbations as well as the effective gravitational potential $\Phi_{eff}$ relevant to the deviation of light rays. We then apply our formulas to a number of modified gravitational models of dark energy--such as those based on f(R) theories, Brans-Dicke theories, kinetic gravity braidings, covariant Galileons, and field derivative couplings with the Einstein tensor. Our results are useful to test the large-distance modification of gravity from the future high-precision observations of large-scale structure, weak lensing, and cosmic microwave background.

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.

We study properties of moduli stabilization in the four dimensional N = 1 supergravity theory with heavy moduli and would-be saxion-axion multiplets including light string-theoretic axions. We give general formulation for the scenario that heavy moduli and saxions are stabilized while axions remain light, assuming that moduli are stabilized near the supersymmetric solution. One can find stable vacuum, i.e. non-tachyonic saxions, in the non-supersymmetric Minkowski vacua. We also discuss the cases, where the moduli are coupled to the supersymmetry breaking sector and/or moduli have contributions to supersymmetry breaking. We also study the models with axions originating from matter-like fields. Our analysis on moduli stabilization is applicable even if there are not light axion multiplets.

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.

We study a vacuum Bianchi IX universe in the context of Ho\v{r}ava-Lifshitz (HL) gravity. In particular, we focus on the classical dynamics of the universe and analyze how anisotropy changes the history of the universe. For small anisotropy, we find an oscillating universe as well as a bounce universe just as the case of the Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime. However, if the initial anisotropy is large, we find the universe which ends up with a big crunch after oscillations if a cosmological constant $\Lambda$ is zero or negative. For $\Lambda>0$, we find a variety of histories of the universe, that is de Sitter expanding universe after oscillations in addition to the oscillating solution and the previous big crunch solution. This fate of the universe shows sensitive dependence of initial conditions, which is one of the typical properties of a chaotic system. If the initial anisotropy is near the upper bound, we find the universe starting from a big bang and ending up with a big crunch for $\Lambda\leq 0$, while de Sitter expanding universe starting from a big bang for $\Lambda>0$.

We study how the dynamical vector degree of freedom in modified gravity affects the CMB B-mode polarization in terms of the Einstein-aether theory. In this theory, vector perturbations can be generated from inflation, which can grow on superhorizon scales in the subsequent epochs and thereby leaves imprints on the CMB B-mode polarization. We derive the linear perturbation equations in a covariant formalism, and compute the CMB B-mode polarization using the CAMB code modified so as to incorporate the effect of the aether vector field. We find that the amplitude of the B-mode signal from the aether field can surpass the contribution from the inflationary gravitational waves for a viable range of model parameters. We also give an analytic argument explaining the shape of the spectrum based on the tight coupling approximation.

We argue that a large negative running spectral index, if confirmed, might suggest that there are abundant structures in the inflaton potential, which result in a fairly large (both positive and negative) running of the spectral index at all scales. It is shown that the center value of the running spectral index suggested by the recent CMB data can be easily explained by an inflaton potential with superimposed periodic oscillations. In contrast to cases with constant running, the perturbation spectrum is enhanced at small scales, due to the repeated modulations. We mention that such features at small scales may be seen by 21 cm observations in the future.

We study the impact of modifying the vector sector of gravity on the CMB polarization. We employ the Einstein-aether theory as a concrete example. The Einstein-aether theory admits dynamical vector perturbations generated during inflation, leaving imprints on the CMB polarization. We derive the perturbation equations of the aether vector field in covariant formalism and compute the CMB B-mode polarization using the modified CAMB code. It is found that the amplitude of the B-mode signal from the aether field can surpass the one from the inflationary gravitational waves. The shape of the spectrum is clearly understood in an analytic way using the tight coupling approximation.

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 study the cosmological inflation from the viewpoint of the moduli stabilization. We study the scenario that the superpotential has a large value during the inflation era enough to stabilize moduli, but it is small in the true vacuum. This scenario is discussed by using a simple model, one type of hybrid models.

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.

A Lifshitz scalar with the dynamical critical exponent z = 3 obtains scale-invariant, super-horizon field fluctuations without the need of an inflationary era. Since this mechanism is due to the special scaling of the Lifshitz scalar and persists in the presence of unsuppressed self-couplings, the resulting fluctuation spectrum can deviate from a Gaussian distribution. We study the non-Gaussian nature of the Lifshitz scalar's intrinsic field fluctuations, and show that primordial curvature perturbations sourced from such field fluctuations can have large non-Gaussianity of order f_NL = O(100), which will be detected by upcoming CMB observations. We compute the bispectrum and trispectrum of the fluctuations, and discuss their configurations in momentum space. In particular, the bispectrum is found to take various shapes, including the local, equilateral, and orthogonal shapes. Intriguingly, all integrals in the in-in formalism can be performed analytically.

We study electromagnetic perturbations around a Kerr black hole surrounded by a thin disk on the equatorial plane. Our main purpose is to reveal the black hole superradiance of electromagnetic waves emitted from the disk surface. The outgoing Kerr-Schild field is used to describe the disk emission, and the superradiant scattering is represented by a vacuum wave field which is added to satisfy the ingoing condition on the horizon. The formula to calculate the energy flux on the disk surface is presented, and the energy transport in the disk-black hole system is investigated. Within the low-frequency approximation we find that the energy extracted from the rotating black hole is mainly transported back to the disk, and the energy spectrum of electromagnetic waves observed at infinity is also discussed.

We study superparticle spectra in the superconformal flavor scenario with non-universal gaugino masses. The non-universality of gaugino masses can lead to the wino-like or higgsino-like neutralino LSP. Furthermore, it is shown that the parameter space for the higgsino-like LSP includes the region where the fine-tuning problem can be improved. The degeneracy of soft scalar masses squared does not drastically change by taking ratios of gaugino masses of order one. The degeneracy of scalar masses for squarks and left-handed sleptons would be good to avoid the FCNC problem but that of right-handed slepton masses is weak. However, the overall size of right-handed slepton masses become larger when the bino becomes heavier. It is also pointed out that such region can be realized, and thus, that would be favorable to avoid the FCNC problem for soft scalar masses as well as A-terms.

We explore the possibility that the dark matter and dark energy are mimicked by a single fluid or by a single $k$-essence-like scalar field. The so called Chaplygin gas unified dark matter models can reproduce the observed matter power spectrum by adding a baryon component. It has been argued that the evolution of the baryon fluctuations is particularly favoured for the "superluminal" case where the sound speed of the Chaplygin gas exceeds the speed of light at late times, as well as for the models with the negligibly small sound speed. In this note we compute the integrated Sachs-Wolfe signal in the Chaplygin gas models, focusing on the superluminal case which has not been investigated before because of the premature understanding of causality. It is shown that the superluminal model leads to large enhancement of the integrated Sachs-Wolfe effect, which is inconsistent with the CMB measurements.

We study the dynamics of isotropic and homogeneous universes in the generalized Ho\v{r}ava-Lifshitz gravity, and classify all possible evolutions of vacuum spacetime. In the case without the detailed balance condition, we find a variety of phase structures of vacuum spacetimes depending on the coupling constants as well as the spatial curvature $K$ and a cosmological constant $\Lambda$. A bounce universe solution is obtained for $\Lambda> 0, K=\pm 1$ or $\Lambda= 0, K=- 1$, while an oscillation spacetime is found for $\Lambda\geq 0, K=1$, or $\Lambda< 0, K=\pm 1$. We also propose a quantum tunneling scenario from an oscillating spacetime to an inflationary universe, resulting in a macroscopic cyclic universe.

We study models of late-time cosmic acceleration in terms of scalar-tensor theories generalized to include a certain class of non-linear derivative interaction of the scalar field. The non-linear effect suppress the scalar-mediated force at short distances to pass solar-system tests of gravity. It is found that the expansion history until today is almost indistinguishable from that of the $\Lambda$CDM model or some (phantom) dark energy models, but the fate of the universe depends clearly on the model parameter. The growth index of matter density perturbations is computed to show that its past asymptotic value is given by 9/16, while the value today is as small as 0.4.

In the inflationary universe, there can be light fields other than the inflaton. We explore a possibility that such light fields source the primordial perturbations, while minimally affecting the inflaton dynamics. We show that during inflation, fluctuations of the light fields can be converted to adiabatic curvature perturbations, which accumulate and become significant by the end of the inflationary era. An additional goal of this work is to distinguish between light fields which can/cannot be ignored during inflation. Such criteria become useful for examining cosmological scenarios with multiple fields. As concrete examples, our results are applied to D-brane inflation models. We consider effects from KK modes (oscillation modes) of wrapped branes in monodromy-driven large-field models, and angular directions of throat geometries in warped D-brane inflation.

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

A scalar-tensor theory of gravity can be made not only to account for the current cosmic acceleration, but also to satisfy solar-system and laboratory constraints, by introducing a non-linear derivative interaction for the scalar field. Such an additional scalar degree of freedom is called "Galileon". The basic idea is inspired by the DGP braneworld, but one can construct a ghost-free model that admits a self-accelerating solution. We perform a fully relativistic analysis of linear perturbations in Galileon cosmology. Although the Galileon model can mimic the background evolution of standard $\Lambda$CDM cosmology, the behavior of perturbation is quite different. It is shown that there exists a super-horizon growing mode in the metric and Galileon perturbations at early times, suggesting that the background is unstable. A fine-tuning of the initial condition for the Galileon fluctuation is thus required in order to promote a desirable evolution of perturbations at early times. Assuming the safe initial condition, we then compute the late-time evolution of perturbations and discuss observational implications in Galileon cosmology. In particular, we find anticorrelations in the cross-correlation of the integrated Sachs-Wolfe effect and large scale structure, similar to the normal branch of the DGP model.

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 obtain cosmological constraints on models of inflation which exhibit rapid roll solutions. Rapid roll attractors exist for potentials with large mass terms and are thus of interest for inflationary model building within string theory. We constrain a general ansatz for the power spectrum arising from rapid roll inflation that, in the small field limit, can be associated with tree level hybrid potentials with variable mass terms and nonminimal gravitational coupling $\xi R\phi^2$. We consider perturbations generated through modulated reheating and/or curvaton mechanisms in place of the observationally unacceptable primary spectra generated by inflaton fluctuations in these models. The lack of a hierarchy amongst higher-order $k$-dependencies of the power spectrum results in models with potentially large runnings, allowing us to impose tight constraints on such models using CMB and LSS data. In particular, we find $n_s <1$ and $|\alpha| < 0.01$. We conclude with a concrete realization of rapid roll inflation within warped throat brane inflation that is in good agreement with current data.

We study the spectral tilt of primordial perturbations in Horava-Lifshitz cosmology. The uniform approximation, which is a generalization of the familiar Wentzel-Kramers-Brillouin (WKB) method, is employed to compute the spectral index both numerically and analytically in a closed form. We clarify how the spectral index depends on the inflation model and parameters in the modified dispersion relation.

We present a curvaton model from type IIB string theory compactified on a warped throat with approximate isometries. Considering an (anti-)D3-brane sitting at the throat tip as a prototype standard model brane, we show that the brane's position in the isometry directions can play the role of curvatons. The basic picture is that the fluctuations of the (anti-)D3-brane in the angular isometry directions during inflation eventually turns into the primordial curvature perturbations, and subsequently the brane's oscillation excites other open string modes on the brane and reheat the universe. We find in the explicit case of the KS throat that a wide range of parameters allows a consistent curvaton scenario. It is also shown that the oscillations of branes at throat tips are capable of producing large non-Gaussianity, either through curvature or isocurvature perturbations. Since such setups naturally arise in warped (multi-)throat compactifications and are constrained by observational data, the model can provide tests for compactification scenarios. This work gives an explicit example of string theory providing light fields for generating curvature perturbations. Such mechanisms free the inflaton from being responsible for the perturbations, thus open up new possibilities for inflation models.

We explore the possibility of generating large-scale magnetic fields from second-order cosmological perturbations during the pre-recombination era. The key process for this is Thomson scattering between the photons and the charged particles within the cosmic plasma. To tame the multi-component interacting fluid system, we employ the tight coupling approximation. It is shown that the source term for the magnetic field is given by the vorticity, which signals the intrinsically second-order quantities, and the product of the first order perturbations. The vorticity itself is sourced by the product of the first-order quantities in the vorticity evolution equation. The magnetic fields generated by this process are estimated to be $\sim 10^{-29} $Gauss on the horizon scale at the recombination epoch. Although our rough estimate suggests that the current generation mechanism can work even on smaller scales, more careful investigation is needed to make clear whether it indeed works in a wide range of spatial scales.

We investigate the effect of modified gravity on cluster abundance and the Sunyaev-Zel'dovich angular power spectrum. Our modified gravity is based on a phenomenological extension of the Dvali-Gabadadze-Porrati model which includes two free parameters characterizing deviation from Lambda-CDM cosmology. Assuming that Birkhoff's theorem gives a reasonable approximation, we study the spherical collapse model of structure formation and show that while the growth function changes to some extent, modified gravity gives rise to no significant change in the linear density contrast at collapse time. The growth function is enhanced in the so called normal branch, while in the "self-accelerating" branch it is suppressed. The Sunyaev-Zel'dovich angular power spectrum is computed in the normal branch, which turns out to be amplified compared to the Lambda-CDM case. This allows us to put observational constraints on the parameters of the modified gravity model using small scale CMB observation data.

Recently it was shown that an inflationary background can be realized by any $p$-form field non-minimally coupled to gravity. In this paper, we study gravitational waves generated during p-form inflation. Even though the background evolution is identical to that in conventional scalar field inflation, the behavior of gravitational waves is different in p-form inflation. In particular, we find that the propagation speed of gravitational waves differs from unity in 2- and 3-form inflationary models. We point out that the squared speed becomes negative in the large field models. The small field models are free from pathologies and the correction to the spectrum of gravitational waves turns out to be very small.

We investigate density perturbations generated through modulated reheating while inflation is driven by a conformally coupled scalar field. A large running of the spectral index is obtained, which reflects the basic nature of conformal inflation that higher-order time derivatives of the Hubble parameter during inflation are not necessarily small. This feature may allow us to distinguish between conformal inflation models and standard minimally coupled ones. We also investigate how the resulting fluctuations are modified when there is a deviation from an exact conformal coupling between the inflaton and gravity. Finally, we apply our results to the warped brane inflation model and see that observational bounds from the WMAP5 data suggest a blue tilted density perturbation spectrum.

Although $f(R)$ modified gravity models can be made to satisfy solar system and cosmological constraints, it has been shown that they have the serious drawback of the nonexistence of stars with strong gravitational fields. In this paper, we discuss whether or not higher curvature corrections can remedy the nonexistence consistently. The following problems are shown to arise as the costs one must pay for the $f(R)$ models that allow for neutrons stars: (i) the leading correction must be fine-tuned to have the typical energy scale $\mu \lesssim 10^{-19}$ GeV, which essentially comes from the free fall time of a relativistic star; (ii) the leading correction must be further fine-tuned so that it is not given by the quadratic curvature term. The second problem is caused because there appears an intermediate curvature scale and laboratory experiments of gravity will be under the influence of higher curvature corrections. Our analysis thus implies that it is a challenge to construct viable $f(R)$ models without very careful and unnatural fine-tuning.

In the context of six-dimensional conical braneworlds we consider a simple and explicit model that incorporates long distance modification of gravity and regularization of codimension-2 singularities. To resolve the conical singularities we replace the codimension-2 branes with ring-like codimension-1 branes, filling in the interiors with regular caps. The six-dimensional Planck scale in the cap is assumed to be much greater than the bulk Planck scale, which gives rise to the effect analogous to brane-induced gravity. Weak gravity on the regularized brane is studied in the case of a sharp conical bulk. We show by a linear analysis that gravity at short distances is effectively described by the four-dimensional Brans-Dicke theory, while the higher dimensional nature of gravity emerges at long distances. The linear analysis breaks down at some intermediate scale, below which four-dimensional Einstein gravity is shown to be recovered thanks to the second-order effects of the brane bending.

We have observed cosmic-ray electrons from 10 GeV to 800 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an average altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates with 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of Iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. From the long duration balloon observations, we derived the energy spectrum of cosmic-ray electrons in the energy range from 100 GeV to 800 GeV. In addition, for the first time we derived the electron arrival directions above 100 GeV, which is consistent with the isotropic distribution.

Several f(R) modified gravity models have been proposed which realize the correct cosmological evolution and satisfy solar system and laboratory tests. Although nonrelativistic stellar configurations can be constructed, we argue that relativistic stars cannot be present in such f(R) theories. This problem appears due to the dynamics of the effective scalar degree of freedom in the strong gravity regime. Our claim thus raises doubts on the viability of f(R) models.

We investigate conformal scalar, electromagnetic, and massless Dirac quasinormal modes of a brane-localized black hole. The background solution is the four-dimensional black hole on a 2-brane that has been constructed by Emparan, Horowitz, and Myers in the context of a lower dimensional version of the Randall-Sundrum model. The conformally transformed metric admits a Killing tensor, allowing us to obtain separable field equations. We find that the radial equations take the same form as in the four-dimensional "braneless" Schwarzschild black hole. The angular equations are, however, different from the standard ones, leading to a different prediction for quasinormal frequencies.

Using the Kerr-Schild formalism to solve the Einstein-Maxwell equations, we study energy transport due to time-dependent electromagnetic perturbations around a Kerr black hole, which may work as a mechanism to illuminate a disk located on the equatorial plane. For such a disk-hole system it is found that the energy extraction from the hole can occur under the well-known superradiance condition for wave frequency, even though the energy absorption into the hole should be rather dominant near the polar region of the horizon. We estimate the efficiency of the superradiant amplification of the disk illumination. Further we calculate the time-averaged energy density distribution to show explicitly the existence of a negative energy region near the horizon and to discuss the possible generation of a hot spot on the disk.

We derive constraints on the tensor to scalar ratio and on the background charge of the warped throat for DBI inflation driven by D5- and D7-branes wrapped over cycles of the throat. It is shown that the background charge well beyond the known maximal value is required in most cases for DBI inflation to generate cosmological observables compatible with the WMAP3 data. Most of the results derived in this paper are insensitive to the details of the inflaton potential, and could be applied to generic warped throats.

We study the effective gravitational theory on a brane in a six-dimensional Einstein-Maxwell model of flux compactification, regularizing a conical defect as a codimension-one brane. We employ the gradient expansion technique valid at low energies. A lowest order analysis showed that standard four-dimensional Einstein gravity is reproduced on the brane. We extend this study to include second order corrections in the effective equations, and show that the correction term is given by a quadratic energy-momentum tensor. Taking the thin-brane limit where the regularized brane shrinks to the pole, we find that the second order metric diverges logarithmically on the brane, giving rise to divergences in the brane effective action. Away from the branes, the effective action is however well-defined.

We derive the low energy effective theory on a brane in six-dimensional chiral supergravity. The conical 3-brane singularities are resolved by introducing cylindrical codimension one 4-branes whose interiors are capped by a regular spacetime. The effective theory is described by the Brans-Dicke (BD) theory with the BD parameter given by $\omega_{\rm BD}=1/2$. The BD field is originated from a modulus which is associated with the scaling symmetry of the system. If the dilaton potentials on the branes preserve the scaling symmetry, the scalar field has an exponential potential in the Einstein frame. We show that the time dependent solutions driven by the modulus in the four-dimensional effective theory can be lifted up to the six-dimensional exact solutions found in the literature. Based on the effective theory, we discuss a possible way to stabilize the modulus to recover standard cosmology and also study the implication for the cosmological constant problem.

Conical brane singularities in six-dimensional flux compactification models can be resolved by introducing cylindrical codimension-one branes with regular caps instead of 3-branes (a la Kaluza-Klein braneworlds with fluxes). In this paper, we consider such a regularized braneworld with axial symmetry in six-dimensional Einstein-Maxwell theory. We derive a low energy effective theory on the regularized brane by employing the gradient expansion approach, and show that standard four-dimensional Einstein gravity is recovered at low energies. Our effective equations extend to the nonlinear gravity regime, implying that conventional cosmology can be reproduced in the regularized braneworld.

We study cosmology on a conical brane in the six-dimensional Einstein-Maxwell-dilaton system, where the extra dimensions are compactified by a magnetic flux. We systematically construct exact cosmological solutions using the fact that the system is equivalently described by (6+n)-dimensional pure Einstein-Maxwell theory via dimensional reduction. In particular, we find a power-law inflationary solution for a general dilatonic coupling. When the dilatonic coupling is given by that of Nishino-Sezgin chiral supergravity, this reduces to the known solution which is not inflating. The power-law solution is shown to be the late-time attractor. We also investigate cosmological tensor perturbations in this model using the (6+n)-dimensional description. We obtain the separable equation of motion and find that there always exist a zero mode, while tachyonic modes are absent in the spectrum. The mass spectrum of Kaluza-Klein modes is obtained numerically.

In the standard cosmological model, magnetic fields and vorticity are generated during the radiation era via second-order density perturbations. In order to clarify the complicated physics of this second-order magnetogenesis, we use a covariant approach and present the electromagneto-dynamical equations in the nonlinear regime. We use the tight-coupling approximation to analyze Thomson and Coulomb scattering. At the zero-order limit of exact tight-coupling, we show that the vorticity is zero and no magnetogenesis takes place at any nonlinear order. We show that magnetogenesis also fails at all orders if either protons or electrons have the same velocity as the radiation, and momentum transfer is neglected. Then we prove a key no-go result: at first-order in the tight-coupling approximation, magnetic fields and vorticity still cannot be generated even via nonlinear effects. The tight-coupling approximation must be broken at first order, for the generation of vorticity and magnetic fields, and we derive a closed set of nonlinear evolution equations that governs this generation. We estimate that the amplitude of the magnetic field at recombination on the horizon scale is $\sim 10^{-27} $G.

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.

We study scalar cosmological perturbations in a braneworld model with a bulk Gauss-Bonnet term. For an anti-de Sitter bulk, the five-dimensional perturbation equations share the same form as in the Randall-Sundrum model, which allows us to obtain metric perturbations in terms of a master variable. We derive the boundary conditions for the master variable from the generalized junction conditions on the brane. We then investigate several limiting cases in which the junction equations are reduced to a feasible level. In the low energy limit, we confirm that the standard result of four-dimensional Einstein gravity is reproduced on large scales, whereas on small scales we find that the perturbation dynamics is described by the four-dimensional Brans-Dicke theory. In the high energy limit, all the non-local contributions drop off from the junction equations, leaving a closed system of equations on the brane. We show that, for inflation models driven by a scalar field on the brane, the Sasaki-Mukhanov equation holds on the high energy brane in its original four-dimensional form.

We have observed atmospheric gamma rays from 30GeV to 8TeV, using emulsion chambers at balloon altitudes, accumulating the largest total exposure in this energy range to date, SOT ~ 6.66m^2.sr.day. At very high altitudes, with residual overburden only a few g/cm^2, atmospheric gamma rays are mainly produced by a single interaction of primary cosmic rays with overlying atmospheric nuclei. Thus, we can use these gamma rays to study the spectrum of primary cosmic rays and their products in the atmosphere. From the observed atmospheric gamma ray spectrum, we deconvolved the primary cosmic-ray proton spectrum, assuming appropriate hadronic interaction models. Our deconvolved proton spectrum covers the energy range from 200GeV to 50TeV, which fills a gap in the currently available primary cosmic-ray proton spectra. We also estimated the atmospheric muon spectrum above 30GeV at high altitude from our gamma-ray spectrum, almost without reference to the primary cosmic rays, and compared the estimated flux with direct muon observations below 10GeV.

We study the holographic representation of the entanglement entropy, recently proposed by Ryu and Takayanagi, in a braneworld context. The holographic entanglement entropy of a de Sitter brane embedded in an anti-de Sitter (AdS) spacetime is evaluated using geometric quantities, and it is compared with two kinds of de Sitter entropy: a quarter of the area of the cosmological horizon on the brane and entropy calculated from the Euclidean path integral. We show that the three entropies coincide with each other in a certain limit. Remarkably, the entropy obtained from the Euclidean path integral is in precise agreement with the holographic entanglement entropy in all dimensions. We also comment on the case of a five-dimensional braneworld model with the Gauss-Bonnet term in the bulk.

The braneworld model of Dvali-Gabadadze-Porrati (DGP) provides an interesting alternative to a positive cosmological constant by modifying gravity at large distances. We investigate the asymptotic behavior of homogeneous and anisotropic cosmologies on the DGP brane. It is shown that all Bianchi models except type IX isotropize, as in general relativity, if the so called $E_{\mu\nu}$ term satisfies some energy condition. Isotropization can proceed slower in DGP gravity than in general relativity.

We study the spectrum of gravitational waves generated from inflation in the Randall-Sundrum braneworld. Since the inflationary gravitational waves are of quantum-mechanical origin, the initial configuration of perturbations in the bulk includes Kaluza-Klein quantum fluctuations as well as fluctuations in the zero mode. We show, however, that the initial fluctuations in Kaluza-Klein modes have no significant effect on the late time spectrum, irrespective of the energy scale of inflation and the equation of state parameter in the post-inflationary stage. This is done numerically, using the Wronskian formulation.

We study the generation and evolution of gravitational waves (tensor perturbations) in the context of Randall-Sundrum braneworld cosmology. We assume that the initial and final stages of the background cosmological model are given by de Sitter and Minkowski phases, respectively, and they are connected smoothly by a radiation-dominated phase. This setup allows us to discuss the quantum-mechanical generation of the perturbations and to see the final amplitude of the well-defined zero mode. Using the Wronskian formulation, we numerically compute the power spectrum of gravitational waves, and find that the effect of initial vacuum fluctuations in the Kaluza-Klein modes is subdominant, contributing not more than 10% of the total power spectrum. Thus it is confirmed that the damping due to the Kaluza-Klein mode generation and the enhancement due to the modification of the background Friedmann equation are the two dominant effects, but they cancel each other, leading to the same spectral tilt as the standard four-dimensional result. Kaluza-Klein gravitons that escape from the brane contribute to the energy density of the dark radiation at late times. We show that a tiny amount of the dark radiation is generated due to this process.

We study the quantum-mechanical generation of gravitational waves during inflation on a brane embedded in a five-dimensional anti-de Sitter bulk. To make the problem well-posed, we consider the setup in which both initial and final phases are given by a de Sitter brane with different values of the Hubble expansion rate. Assuming that the quantum state is in a de Sitter invariant vacuum in the initial de Sitter phase, we numerically evaluate the amplitude of quantum fluctuations of the growing solution of the zero mode in the final de Sitter phase. We find that the vacuum fluctuations of the initial Kaluza-Klein gravitons as well as of the zero mode gravitons contribute to the final amplitude of the zero mode on small scales, and the power spectrum is quite well approximated by what we call the rescaled spectrum, which is obtained by rescaling the standard four-dimensional calculation following a simple mapping rule. Our results confirm the speculation raised in Ref. \cite{Kobayashi:2003cn} before.

We consider a Vaidya-type radiating spacetime in Einstein gravity with the Gauss-Bonnet combination of quadratic curvature terms. Simply generalizing the known static black hole solutions in Einstein-Gauss-Bonnet gravity, we present an exact solution in arbitrary dimensions with the energy-momentum tensor given by a null fluid form. As an application, we derive an evolution equation for the ``dark radiation'' in the Gauss-Bonnet braneworld.

Tensor type perturbations in the expanding brane world of the Randall Sundrum type are investigated. We consider a model composed of slow-roll inflation phase and the succeeding radiation phase. The effect of the presence of an extra dimension through the transition to the radiation phase is studied, giving an analytic formula for leading order corrections.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

We study the $D$-term inflation scenario with a nonperturbative K\"ahler potential of the dilaton field. Although the FI term which leads an inflationary expansion is given by the derivative of the K\"ahler potential with respect to the dilaton in heterotic string models with anomalous U(1), the too large magnitude is problematic for a viable $D$-term inflation. In this paper, we point out that the K\"ahler potential with a nonperturbative term can reduce the magnitude of FI term to desired values while both the dilaton stabilization and $D$-term domination in the potential are realized by nonperturbative superpotential.

Evidences of non-thermal X-ray emission and TeV gamma-rays from the supernova remnants (SNRs) has strengthened the hypothesis that primary Galactic cosmic-ray electrons are accelerated in SNRs. High energy electrons lose energy via synchrotron and inverse Compton processes during propagation in the Galaxy. Due to these radiative losses, TeV electrons liberated from SNRs at distances larger than ~1 kpc, or times older than ~10^5 yr, cannot reach the solar system. We investigated the cosmic-ray electron spectrum observed in the solar system using an analytical method, and considered several candidate sources among nearby SNRs which may contribute to the high energy electron flux. Especially, we discuss the effects for the release time from SNRs after the explosion, as well as the deviation of a source spectrum from a simple power-law. From this calculation, we found that some nearby sources such as the Vela, Cygnus Loop, or Monogem could leave unique signatures in the form of identifiable structure in the energy spectrum of TeV electrons and show anisotropies towards the sources, depending on when the electrons are liberated from the remnant. This suggests that, in addition to providing information on the mechanisms of acceleration and propagation of cosmic-rays, specific cosmic-ray sources can be identified through the precise electron observation in the TeV region.

We investigate the possibility of $D$-term inflation within the framework of type I string-inspired models. Although $D$-term inflation model has the excellent property that it is free from the so-called $\eta$- problem, two serious problems appear when we embed $D$-term inflation in string theory, the magnitude of FI term and the rolling motion of the dilation. In the present paper, we analyze the potential of $D$-term inflation in type I inspired models and study the behavior of dilaton and twisted moduli fields. Adopting the non-perturbative superpotential induced by gaugino condensation, the twisted moduli can be stabilized. If the dilaton is in a certain range, it evolves very slowly and does not run away to infinity. Thus $D$-term dominated vacuum energy becomes available for driving inflation. By studying the density perturbation generated by the inflation model, we derive the constraints on model parameters and give some implications on $D$-term inflation in type I inspired models.

We compare predictions for the spin-independent contribution to the neutralino-proton scattering cross section $\sig$ and for the anomalous magnetic dipole moment of the muon, $a_\mu = (g_\mu-2)/2$, in models with gravity-mediated supersymmetry breaking. We nearly always find a positive correlation between these two measurables, i.e. scenarios with larger $a_\mu$ also tend to have larger $\sig$, but the detailed prediction differs greatly between models. In particular, we find that for the popular mSUGRA scenario with universal soft breaking masses at the scale of Grand Unification, measurements of $a_\mu$ currently seem more promising. On the other hand, if scalar soft breaking masses at the GUT scale receive sizable contributions from SO(10) D-terms, one often finds scenarios with large $\sig$ but $a_\mu$ below the currently foreseen sensitivity. A string-inspired model with non-universal scalar spectrum at the GUT scale falls between these two cases.

We show that LHC experiments might well be able to determine all the parameters required for a prediction of the present density of thermal LSP relics from the Big Bang era. If the LSP is an almost pure bino we usually only need to determine its mass and the mass of the SU(2) singlet sleptons. This information can be obtained by reconstructing the cascade $\tilde{q}_L \to \tilde{\chi}_2^0 q \to \tilde{\ell}_R \ell q \to \tilde{\chi}_1^0 \ell^+ \ell^- q$. The only requirement is that $m_{\tilde{\ell}_R} < m_{\tilde{\chi}_2^0}$, which is true for most of the cosmologically interesting parameter space. If the LSP has a significant higgsino component, its predicted thermal relic density is smaller than for an equal--mass bino. We show that in this case squark decays also produce significant numbers of $\tilde{\chi}_4^0$ and $\tilde{\chi}_2^\pm$. Reconstructing the corresponding decay cascades then allows to determine the higgsino component of the LSP.

We apply the renormalization group (RG) method to examine the observable scaling properties in Newtonian cosmology. The original scaling properties of the equations of motion in our model are modified for averaged observables on constant time slices. In the RG flow diagram, we find three robust fixed points: Einstein-de Sitter, Milne and Quiescent fixed points. Their stability (or instability) property does not change under the effect of fluctuations. Inspired by the inflationary scenario in the early Universe, we set the Einstein-de Sitter fixed point with small fluctuations, as the boundary condition at the horizon scale. Solving the RG equations under this boundary condition toward the smaller scales, we find generic behavior of observables such that the density parameter $\Omega$ decreases, while the Hubble parameter $H$ increases for smaller averaging volume. The quantitative scaling properties are analyzed by calculating the characteristic exponents around each fixed point. Finally we argue the possible fractal structure of the Universe beyond the horizon scale.

The $\beta$-delayed neutron decay of $^{17}$B was studied using a radioactive ion beam. The neutron energies, measured via time-of-flight, give information on states in $^{17}$C above the $^{16}$C + neutron threshold. States in $^{17}$C were found at excitation energies of 2.25(2) MeV, 2.64(2) MeV and 3.82(5) MeV, and possibly at 1.18(1) MeV. These low lying states are of possible interest for nuclear physics as well as for astrophysics.