We study the possibility to directly detect the boosted dark matter generated from the scatterings with high energetic cosmic particles such as protons and electrons. As a concrete example, we consider the sub-GeV dark matter mediated by a $U(1)_D$ gauge boson which has mixing with $U(1)_Y$ gauge boson in the standard model. The enhanced kinetic energy of the light dark matter from the collision with the cosmic rays can recoil the target nucleus and electron in the underground direct detection experiments transferring enough energy to them to be detectable. We show the impact of BDM with existing direct detection experiments as well as collider and beam-dump experiments.

Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

Recently it has been noticed that the Fermi-LAT data of gamma-rays from some galactic pulsars and supernova remnants reveal spectral modulations that might be explained by the conversion of photons to ALPs (axion-like particles) induced by the conventional ALP coupling to photon in the presence of galactic magnetic fields. However the corresponding ALP mass and coupling are in a severe tension with the observational constraints from CAST, SN1987A, and other gamma-ray observations. Motivated by this, we examine an alternative possibility that those spectral modulations are explained by other type of ALP coupling involving both the ordinary photon and a massless dark photon, when nonzero background dark photon gauge fields are assumed. We find that our scheme results in oscillations among the photon, ALP, and dark photon, which can explain the gamma-ray spectral modulations of galactic pulsars or supernova remnants, while satisfying the known observational constraints.

We propose a new method to detect gravitational waves, based on spatial coherence interferometry with stellar light, as opposed to the conventional temporal coherence interferometry with laser sources. The proposed method detects gravitational waves by combining two coherent beams of light from a single distant star measured at two separate detectors, each on-board a satellite. For this experiment, no complex detector system is required, hence any associated acceleration noise can be minimized, allowing for the search of gravitational waves in the lower frequency range of $10^{-7} - 10^{-4}$ Hz, not covered by any experiment at the time of this paper. In this work, we present the detection sensitivity of the proposed stellar interferometer for gravitational waves by taking the detection response and potential noises into account. Furthermore, the proposed experimental setup is capable of searching for primordial black holes and studying the size of the target neutron star, which are also discussed in the paper.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

The Borexino detector measures solar neutrino fluxes via neutrino-electron elastic scattering. Observed spectra are determined by the solar-$\nu_{e}$ survival probability $P_{ee}(E)$, and the chiral couplings of the neutrino and electron. Some theories of physics beyond the Standard Model postulate the existence of Non-Standard Interactions (NSI's) which modify the chiral couplings and $P_{ee}(E)$. In this paper, we search for such NSI's, in particular, flavor-diagonal neutral current interactions that modify the $\nu_e e$ and $\nu_\tau e$ couplings using Borexino Phase II data. Standard Solar Model predictions of the solar neutrino fluxes for both high- and low-metallicity assumptions are considered. No indication of new physics is found at the level of sensitivity of the detector and constraints on the parameters of the NSI's are placed. In addition, with the same dataset the value of $\sin^2\theta_W$ is obtained with a precision comparable to that achieved in reactor antineutrino experiments.

We investigate the detectability of axino-like particle, which is defined as a supersymmetric partner of axion-like particle and can be a good candidate for dark matter in our Universe. Especially, we consider the fixed target experiments to search for the light axino-like particle with a neutralino as the next-to-lightest supersymmetric particle. We calculate the production and decay rate of neutralinos and the consequent number of events (such as photons and charged leptons) that are produced when the neutralinos decay to the axino-like particles.

We present the first simultaneous measurement of the interaction rates of $pp$, $^7$Be, and $pep$ solar neutrinos performed with a global fit to the Borexino data in an extended energy range (0.19-2.93)$\,$MeV. This result was obtained by analyzing 1291.51$\,$days of Borexino Phase-II data, collected between December 2011 and May 2016 after an extensive scintillator purification campaign. We find: rate($pp$)$\,$=$\,$$134$$\,$$\pm$$\,$$10$$\,$($stat$)$\,$$^{\rm +6}_{\rm -10}$$\,$($sys$)$\,$cpd/100$\,$t, rate($^7$Be)$\,$=$\,$$48.3$$\,$$\pm$$\,$$1.1$$\,$($stat$)$\,$$^{\rm +0.4}_{\rm -0.7}$$\,$($sys$)$\,$cpd/100$\,$t, and rate($pep$)$\,$=$\,$$2.43$$\pm$$\,$$0.36$$\,$($stat$)$^{+0.15}_{-0.22}$$\,$($sys$)$\,$cpd/100$\,$t. These numbers are in agreement with and improve the precision of our previous measurements. In particular, the interaction rate of $^7$Be $\nu$'s is measured with an unprecedented precision of 2.7%, showing that discriminating between the high and low metallicity solar models is now largely dominated by theoretical uncertainties. The absence of $pep$ neutrinos is rejected for the first time at more than 5$\,$$\sigma$. An upper limit of $8.1$$\,$cpd/100$\,$t (95%$\,$C.L.) on the CNO neutrino rate is obtained by setting an additional constraint on the ratio between the $pp$ and $pep$ neutrino rates in the fit. This limit has the same significance as that obtained by the Borexino Phase-I (currently providing the tightest bound on this component), but is obtained by applying a less stringent constraint on the $pep$ $\nu$ flux.

With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos.

We examine the resonant conversion of axion-like particle (ALP) or dark photon to the electromagnetic photon in the early Universe, which takes place due to the ALP-photon-dark photon oscillations in background dark photon gauge fields. It is noted that the corresponding conversion probability can have an unusual spectral feature which allows strong conversion at low frequency domain, but has negligible conversion at high frequencies above certain critical frequency which is determined by the ALP coupling to dark photon and the strength of background dark photon gauge field. We apply this scheme to heat up the 21 cm photons without affecting the Cosmic Microwave Background, which can explain the tentative absorption signal of 21 cm photons detected recently by the EDGES experiment.

Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $\tau$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_\nu^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_\nu/{\rm d}E_\nu < 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

We study neutrino oscillations in a medium of dark matter which generalizes the standard matter effect. A general formula is derived to describe the effect of various mediums and their mediators to neutrinos. Neutrinos and anti-neutrinos receive opposite contributions from asymmetric distribution of (dark) matter and anti-matter, and thus it could appear in precision measurement of neutrino or anti-neutrino oscillations. Furthermore, the standard neutrino oscillation can occur from the symmetric dark matter effect even for massless neutrinos.

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.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

Astrophysical neutrinos travel long distances from their sources to the Earth traversing dark matter halos of clusters of galaxies and that of our own Milky Way. The interaction of neutrinos with dark matter may affect the flux of neutrinos. The recent multi-messenger observation of a high energy neutrino, IceCube-170922A, can give a robust upper bound $\sigma /M_{dm} \lesssim 5.1\times 10^{-23} {\rm cm}^2 /$GeV on the interaction between neutrino and dark matter at a neutrino energy of 290 TeV allowing 90\% suppression. Combining the constraints from CMB and LSS at different neutrino energies, we can constrain models of dark matter-neutrino interactions.

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.

We consider the thermal production of axino dark matter in high-scale supersymmetry where all the superpartners except the axino are heavier than the maximum and reheating temperatures. In this case, the axinos are produced dominantly in pairs from the scattering of SM particles in thermal plasma in the early Universe after inflation. We find that the thermal averaged scattering cross section for the axino pair production is given by $\langle\sigma v\rangle \propto T^4$ in Kim-Shifman-Vainstein-Zakharov (KSVZ) axion model, while it does not depend on the temperature in Dine-Fischler-Srednicki-Zhitnitski (DFSZ) axion model. As a result, the axinos produced during the early matter domination is diluted by the entropy production, so the axino abundance is determined mainly by the reheating temperature, unlike the case with gravitino dark matter. We show that the axino pair production in DFSZ model opens up new parameter space for axino dark matter, due to non-decoupled Higgsino interactions at tree level.

We propose a mechanism generating primordial magnetic fields after the $e^+e^-$ annihilations. Our mechanism involves an ultra-light axion-like particle (ALP) which constitutes the dark matter, and a dark $U(1)_X$ gauge boson introduced to bypass the obstacle placed by the conductivity of cosmic plasma. In our scheme, a coherently oscillating ALP amplifies the dark photon field, and part of the amplified dark photon field is concurrently converted to the ordinary magnetic field through the ALP-induced magnetic mixing. For the relevant ALP mass range $10^{-21} {\rm eV}\lesssim m_\phi\lesssim 10^{-17}{\rm eV}$, our mechanism can generate $B\sim 10^{-24} \,{\rm G} \,(m_\phi/10^{-17} {\rm eV})^{5/4}$ with a coherent length $\lambda \sim (m_\phi/10^{-17} {\rm eV})^{-1/2}$ kpc, which is large enough to provide a seed of the galactic magnetic fields. The mechanism also predicts a dark $U(1)_X$ electromagnetic field $E_X \sim B_X\sim 80\,{\rm nG}\, (m_\phi/10^{-17}{\rm eV})^{-1/4}$, which can result in interesting astrophysical/cosmological phenomena by inducing the mixings between the ALP, ordinary photon, and dark photon states.

We propose a WIMP baryogensis achieved by the annihilation of non-thermally produced WIMPs from decay of heavy particles, which can result in low reheating temerature. Dark matter (DM) can be produced non-thermally during a reheating period created by the decay of long-lived heavy particle, and subsequently re-annihilate to lighter particles even after the thermal freeze-out. The re-annihilation of DM provides the observed baryon asymmetry as well as the correct relic density of DM. We investigate how wahout effects can affect the generation of the baryon asymmetry and study a model suppressing them. In this scenario, we find that DM can be heavy enough and its annihilation cross section can also be larger than that adopted in the usual thermal WIMP baryogenesis.

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 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 present the results of a low-energy neutrino search using the Borexino detector in coincidence with the gravitational wave (GW) events GW150914, GW151226 and GW170104. We searched for correlated neutrino events with energies greater than 250 keV within a time window of $\pm500$ s centered around the GW detection time. A total of five candidates were found for all three GW150914, GW151226 and GW170104. This is consistent with the number of expected solar neutrino and background events. As a result, we have obtained the best current upper limits on the GW event neutrino fluence of all flavors ($\nu_e, \nu_{\mu}, \nu_{\tau}$) in the energy range $(0.5 - 5.0)$ MeV.

We investigate a new bound on the low reheating temperature in a scenario where the Universe experiences early matter-domination before reheating after which the standard big bang cosmology begins. In many models of dark matter (DM), the small scale fluctuations of DM grow during the early matter-domination era and seed the formation of the ultracompact minihalos (UCMHs). Using the constraints on the number of UCMHs from gamma-ray observations, we find a lower bound on the reheating temperature between ${\cal O}(10)~{\rm MeV} - {\cal O}(100)~{\rm MeV}$ for WIMP dark matter depending on the nature of DM. A similar bound could be obtained for non-WIMP dark matter by observing UCMHs gravitationally such as pulsar timing, microlensing and so on in some future observations.

Many extensions of Standard Model (SM) include a dark sector which can interact with the SM sector via a light mediator. We explore the possibilities to probe such a dark sector by studying the distortion of the CMB spectrum from the blackbody shape due to the elastic scatterings between the dark matter and baryons through a hidden light mediator. We in particular focus on the model where the dark sector gauge boson kinetically mixes with the SM and present the future experimental prospect for a PIXIE-like experiment along with its comparison to the existing bounds from complementary terrestrial experiments.

We examine the low energy phenomenology of the relaxion solution to the weak scale hierarchy problem. Assuming that the Hubble friction is responsible for a dissipation of the relaxion energy, we identify the cosmological relaxion window which corresponds to the parameter region compatible with a given value of the acceptable number of inflationary $e$-foldings. We then discuss a variety of observational constraints on the relaxion window, including those from astrophysical and cosmological considerations. We find that majority of the parameter space with a relaxion mass $m_\phi\gtrsim 100$ eV or a relaxion decay constant $f\lesssim 10^7$ GeV is excluded by existing constraints. There is an interesting parameter region with $m_\phi\sim \,0.2-10$ GeV and $f\sim\, {\rm few}-200$ TeV, which is allowed by existing constraints, but can be probed soon by future beam dump experiments such as the SHiP experiment, or by improved EDM experiments.

A search for neutrino and antineutrino events correlated with 2,350 gamma-ray bursts (GRBs) is performed with Borexino data collected between December 2007 and November 2015. No statistically significant excess over background is observed. We look for electron antineutrinos ($\bar{\nu}_e$) that inverse beta decay on protons with energies from 1.8\,MeV to 15\,MeV and set the best limit on the neutrino fluence from GRBs below 8\,MeV. The signals from neutrinos and antineutrinos from GRBs that scatter on electrons are also searched for, a detection channel made possible by the particularly radio-pure scintillator of Borexino. We obtain currently the best limits on the neutrino fluence of all flavors and species below 7\,MeV. Finally, time correlations between GRBs and bursts of events are investigated. Our analysis combines two semi-independent data acquisition systems for the first time: the primary Borexino readout optimized for solar neutrino physics up to a few MeV, and a fast waveform digitizer system tuned for events above 1\,MeV.

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

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.

The weak gravity conjecture applied for the aligned natural inflation indicates that generically there can be a modulation of the inflaton potential, with a period determined by sub-Planckian axion scale. We study the oscillations in the primordial power spectrum induced by such modulation, and discuss the resulting observational constraints on the model.

We study the primordial scalar and tensor perturbations in inflation scenario involving a spectator dilaton field. In our setup, the rolling spectator dilaton causes a tachyonic instability of gauge fields, leading to a copious production of gauge fields in the superhorizon regime, which generates additional scalar and tensor perturbations through gravitational interactions. Our prime concern is the possibility to enhance the tensor-to-scalar ratio $r$ relative to the standard result, while satisfying the observational constraints. To this end, we allow the dilaton field to be stabilized before the end of inflation, but after the CMB scales exit the horizon. We show that for the inflaton slow roll parameter $\epsilon \gtrsim 10^{-3}$, the tensor-to-scalar ratio in our setup can be enhanced only by a factor of ${\cal O}(1)$ compared to the standard result. On the other hand, for smaller $\epsilon$ corresponding to a lower inflation energy scale, a much bigger enhancement can be achieved, so that our setup can give rise to an observably large $r\gtrsim 10^{-2}$ even when $\epsilon\ll 10^{-3}$. The tensor perturbation sourced by the spectator dilaton can have a strong scale dependence, and is generically red-tilted. We also discuss a specific model to realize our scenario, and identify the parameter region giving an observably large $r$ for relatively low inflation energy scales.

We study the primordial scalar and tensor perturbations in inflation scenario involving a spectator dilaton field. In our setup, the rolling spectator dilaton causes a tachyonic instability of gauge fields, leading to a copious production of gauge fields in the superhorizon regime, which generates additional scalar and tensor perturbations through gravitational interactions. Our prime concern is the possibility to enhance the tensor-to-scalar ratio $r$ relative to the standard result, while satisfying the observational constraints. To this end, we allow the dilaton field to be stabilized before the end of inflation, but after the CMB scales exit the horizon. We show that for the inflaton slow roll parameter $\epsilon \gtrsim 10^{-3}$, the tensor-to-scalar ratio in our setup can be enhanced only by a factor of ${\cal O}(1)$ compared to the standard result. On the other hand, for smaller $\epsilon$ corresponding to a lower inflation energy scale, a much bigger enhancement can be achieved, so that our setup can give rise to an observably large $r\gtrsim 10^{-2}$ even when $\epsilon\ll 10^{-3}$. The tensor perturbation sourced by the spectator dilaton can have a strong scale dependence, and is generically red-tilted. We also discuss a specific model to realize our scenario, and identify the parameter region giving an observably large $r$ for relatively low inflation energy scales.

The adiabatic perturbation of dark matter is damped during the kinetic decoupling due to the collision with relativistic component on sub-horizon scales. However the isocurvature part is free from damping and could be large enough to make a substantial contribution to the formation of small scale structure. We explicitly study the weakly interacting massive particles as dark matter with an early mater dominated period before radiation domination and show that the isocurvature perturbation is generated during the phase transition and leaves imprint in the observable signatures for small scale structure.

We propose a new mechanism to suppress the axion isocurvature perturbation, while producing the right amount of axion dark matter, within the framework of supersymmetric axion models with the axion scale induced by supersymmetry breaking. The mechanism involves an intermediate phase transition to generate the Higgs \mu-parameter, before which the weak scale is comparable to the axion scale and the resulting stronger QCD yields an axion mass heavier than the Hubble scale over a certain period. Combined with that the Hubble-induced axion scale during the primordial inflation is well above the intermediate axion scale at present, the stronger QCD in the early Universe suppresses the axion fluctuation to be small enough even when the inflationary Hubble scale saturates the current upper bound, while generating an axion misalignment angle of order unity.

We present a perspective on the inflation paths in 2-, 3-,,, N-flation models based on the ultraviolet completion in heterotic string theory, where a number of grand unification scale axions are used. The number of non-Abelian gauge groups for a natural inflation is restricted in string compactification, and we argue that the most plausible completion of natural inflation from a theory perspective is the 2-flation.

Increasingly stringent limits from LHC searches for new physics, coupled with lack of convincing signals of weakly interacting massive particle (WIMP) in dark matter searches, have tightly constrained many realizations of the standard paradigm of thermally produced WIMPs as cold dark matter. In this article, we review more generally both thermally and non-thermally produced dark matter (DM). One may classify DM models into two broad categories: one involving bosonic coherent motion (BCM) and the other involving WIMPs. BCM and WIMP candidates need, respectively, some approximate global symmetries and almost exact discrete symmetries. Supersymmetric axion models are highly motivated since they emerge from compelling and elegant solutions to the two fine-tuning problems of the Standard Model: the strong CP problem and the gauge hierarchy problem. We review here non-thermal relics in a general setup, but we also pay particular attention to the rich cosmological properties of various aspects of mixed SUSY/axion dark matter candidates which can involve both WIMPs and BCM in an interwoven manner. We also review briefly a panoply of alternative thermal and non-thermal DM candidates.

We study a supersymmetric neutrinophilic Higgs model with large neutrino Yukawa couplings where neutrinos are Dirac particles and the lightest right-handed (RH) sneutrino is the lightest supersymmetric particle (LSP) as a dark matter candidate. Neutrinophilic Higgs bosons need to be rather heavy by the precise determination of the muon decay width and dark radiation constraints for large Yukawa couplings. From the Large Hadron Collider constraints, neutrinophilic Higgsino mass need to be heavier than several hundred GeV or close to the RH sneutrino LSP mass. The latter case is interesting because the muon anomalous magnetic dipole moment can be explained with a relatively large lightest neutrino mass, if RH sneutrino mixings are appropriately fine tuned in order to avoid stringent lepton flavor violation constraints. Dark matter is explained by asymmetric RH sneutrino dark matter in the favoured region by the muon anomalous magnetic dipole moment. In other regions, RH sneutrino could be an usual WIMP dark matter.

We extend the Kim-Nilles-Peloso (KNP) alignment mechanism for natural inflation to models with $N>2$ axions, which obtains a super-Planckian effective axion decay constant $f_{\textrm{eff}}\gg M_{Pl}$ through an alignment of the anomaly coefficients of multiple axions having sub-Planckian fundamental decay constants $f_0\ll M_{Pl}$. The original version of the KNP mechanism realized with two axions requires that some of the anomaly coefficients should be of the order of $f_{\textrm{eff}}/f_0$, which would be uncomfortably large if $f_{\rm eff}/f_0 \gtrsim {\cal O}(100)$ as suggested by the recent BICEP2 results. We note that the KNP mechanism can be realized with the anomaly coefficients of $\mathcal{O}(1)$ if the number of axions $N$ is large as $N\ln N\gtrsim 2\ln (f_{\textrm{eff}}/f_0)$, in which case the effective decay constant can be enhanced as $f_{\rm eff}/f_0 \sim \sqrt{N !}\,n^{N-1}$ for $n$ denoting the typical size of the integer-valued anomaly coefficients. Comparing to the other multiple axion scenario, the N-flation scenario which requires $N \sim f_{\textrm{eff}}^2/f_0^2$, the KNP mechanism has a virtue of not invoking to a too large number of axions, although it requires a specific alignment of the anomaly coefficients, which can be achieved with a probability of ${\cal O}(f_0/f_{\rm eff})$ under a random choice of the anomaly coefficients. We also present a simple model realizing a multiple axion monodromy along the inflaton direction.

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

The large tensor spectrum recently observed by the BICEP2 Collaboration requires a super-Planckian field variation of the inflaton in the single-field inflationary scenario. The required slow-roll parameter \epsilon \approx 0.01 would restrict the e-folding number to around 7 in (sub-)Planckian inflationary models. To overcome such problems, we consider a two-field scenario based on the natural assisted supersymmetric (SUSY) hybrid model ("natural SUSY hybrid inflation" [1]), which combines the SUSY hybrid and the natural inflation models. The axionic inflaton field from the natural inflation sector can admit the right values for the tensor spectrum as well as a spectral index of 0.96 with a decay constant smaller than the Planck scale, f \lesssim M_P. On the other hand, the vacuum energy of 2 x 10^{16} GeV with 50 e-folds is provided by the inflaton coming from the SUSY hybrid sector, avoiding the eta problem. These are achieved by introducing both the U(1)_R and a shift symmetry, and employing the minimal Kahler potential.

We consider axino warm dark matter in a supersymmetric axion model with R-parity violation. In this scenario, axino with the mass $m_\axino\simeq 7$ keV can decay into photon and neutrino resulting in the X-ray line signal at $3.5$ keV, which might be the origin of unidentified X-ray emissions from galaxy clusters and Andromeda galaxy detected by the XMM-Newton X-ray observatory.

Dark matter could be captured in the Sun and self-annihilate, giving rise to an observable neutrino flux. Indirect searches for dark matter looking for this signal with neutrino telescopes have resulted in tight constraints on the interaction cross-section of dark matter with ordinary matter. We investigate how robust limits are against astro-physical uncertainties. We study the effect of the velocity distribution of dark matter in our Galaxy on capture rates in the Sun. We investigate four sources of uncertainties: orbital speed of the Sun, escape velocity of dark matter from the halo, dark matter velocity distribution functions and existence of a dark disc. We find that even extreme cases currently discussed do not decrease the sensitivity of indirect detection significantly because the capture is achieved over a broad range of the velocity distribution by integration over the velocity distribution. The effect of the uncertainty in the high-velocity tail of dark matter halo is very marginal as the capture process is rather inefficient at this region. The difference in capture rate in the Sun for various scenarios is compared to the expected change in event rates for direct detection. The possibility of co-rotating structure with the Sun can largely boost the signal and hence makes the interpretation of indirect detection conservative compared to direct detection.

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.

Recently, it was shown, that the absence of circular polarization of visible light from quasars, severely constrains the interpretation of axion like particles (ALPs) as a solution for the generation of the linear polarization. Furthermore, the new observation of the linear polarization in the radio wavelength from quasars, similar to the earlier observation performed in the optical band, makes the ALPs scenario inconsistent with at least one of the two observations. In this study we extend this scenario, to two axion like particles, one scalar and another pseudoscalar. We find, that the effects from scalar and pseudoscalar cancel out each other, thereby suppressing the circular polarization, while preserving consistent linear polarization, as observed in both visible and radio waves bands.

We review the status of axino dark matter. Two hierarchy problems, the strong CP problem and the gauge hierarchy problem, have led to introducing into particle physics a spontaneously broken global Peccei-Quinn symmetry and a softly broken supersymmetry, respectively. Combining them implies the presence of not only an axion, but also of its scalar component, saxion, and their fermionic partner, axino. Among these, the axion and the axino are attractive dark matter candidates. Various possibilities for the axino as dark matter are discussed.

We propose a renormalizable model for vector dark matter with extra U(1) gauge symmetry, which is broken by the VEV of a complex singlet scalar. When the singlet scalar has a quartic coupling to a heavy charged scalar, the resonance effect enhances the annihilation cross section of vector dark matter into two photons such that Fermi-gamma ray line at about 130 GeV is obtained. In the presence of a tiny mixing between the singlet scalar and the Standard Model Higgs doublet, the relic density is determined dominantly by WW/ZZ and two-photon channels near the resonance pole of the singlet scalar. We also show that various phenomenological bounds coming from the Higgs to diphoton decay rate, precision data and collider searches for the charged scalar and vacuum stability are satisfied in the model.

The nucleon form factors in free space are usually thought to be modified when a nucleon is bound in a nucleus or immersed in a nuclear medium. We investigate effects of the density-dependent axial and weak-vector form factors on the electro-neutrino ($\nu_e$) and anti-electro-neutrino $({\bar \nu_e})$ reactions via neutral current (NC) for a nucleon in nuclear medium or $^{12}$C. For the density-dependent form factors, we exploit the quark-meson-coupling (QMC) model, and apply them to the $\nu_e$ and ${\bar \nu_e}$ induced reactions by NC. About 12% decrease of the total cross section by $\nu_e$ reaction on the nucleon is obtained at normal density, $\rho = \rho_0 \sim 0.15 {fm}^{-3} $, as well as about 18% reduction of total ${\nu}_e$ cross section on $^{12}$C, by the modification of the weak form factors of the bound nucleon. However, similarly to the charged current reaction, effects of the nucleon property change in the ${\bar \nu}_e$ reaction reduce significantly the cross sections about 30% for the nucleon in matter and $^{12}$C cases. Such a large asymmetry in the ${\bar \nu}_e$ cross sections is addressed to originate from the different helicities of ${\bar \nu}_e$ and ${\nu}_e$.

We study the effect of the density-dependent axial and vector form factors on the electro-neutrino ($\nu_e$) and anti-neutrino $({\bar \nu}_e)$ reactions for a nucleon in nuclear matter or in $^{12}$C. The nucleon form factors in free space are presumed to be modified for a bound nucleon in a nuclear medium. We adopt the density-dependent form factors calculated by the quark-meson coupling (QMC) model, and apply them to the $\nu_e$ and ${\bar \nu}_e$ induced reactions with the initial energy $E = $ 8 $\sim$ 80 MeV. We find that the total ${\nu}_e$ cross sections on $^{12}$C as well as a nucleon in nuclear matter are reduced by about 5% at the nuclear saturation density, $\rho_0$. This reduction is caused by the modification of the nucleon structure in matter. Although the density effect for both cases is relatively small, it is comparable with the effect of Coulomb distortion on the outgoing lepton in the $\nu$-reaction. In contrast, the density effect on the ${\bar \nu}_e$ reaction reduces the cross section significantly in both nuclear matter and $^{12}$C cases, and the amount maximally becomes of about 35% around $\rho_0$. Such large asymmetry in the $\nu_e$ and ${\bar \nu}_e$ cross sections, which seems to be nearly independent of the target, is originated from the difference in the helicities of ${\bar \nu}_e$ and ${\nu}_e$. It is expected that the asymmetry influences the r-process and also the neutrino-process nucleosynthesis in core-collapse supernovae.

We study the Standard Model Higgs field as a source for the primordial curvature perturbation, particularly in the curvaton and modulated reheating scenario. We conclude that the Higgs cannot play as a curvaton due to the small energy density when it decays, however the modulated reheating by Higgs can be a viable scenario. In the latter case, the non-Gaussianity is inevitably generated and strongly constrains the type of potential of inflaton field and Higgs-dependent interaction term. For the quadratic potential of the inflaton field with decay rate which non-linearly depends on the Higgs vacuum expectation value, the contribution of Higgs field to the primordial curvature perturbation must be less than 8%.

We discuss asymmetric or symmetric dark matter candidate in the supersymmetric Dirac leptogenesis scenario. By introducing a singlet superfield coupling to right-handed neutrinos, the overabundance problem of dark matter can be evaded and various possibilities for dark matter candidate arise. If the singlino is the lightest supersymmetric particle (LSP), it becomes naturally asymmetric dark matter. On the other hand, the right-handed sneutrino is a symmetric dark matter candidate whose relic density can be determined by the usual thermal freeze-out process. The conventional neutralino or gravitino LSP can be also a dark matter candidate as its non-thermal production from the right-handed sneutrino can be controlled appropriately. In our scenario, the late-decay of heavy supersymmetric particles mainly produce the right-handed sneutrino and neutrino which is harmless to the standard prediction of the big-bang nucleosynthesis.

Although the standard {\Lambda}CDM model describes the cosmic microwave background radiation and the large scale structure of the Universe with great success, it has some tensions with observations in the effective number of neutrino species (dark radiation) and the number of small scale structures (overabundance problem). Here we propose a scenario which can relax these tensions by producing both dark matter and dark radiation by late decays of heavy particle. Thanks to the generation mechanism, dark matters are rather warm so that the small-scale structure problem is resolved. This scenario can be naturally realized in supersymmetric axion model, in which axions produced by saxion decays provide dark radiation, while axinos from saxion decays form warm dark matter. We identify a parameter region of supersymmetric axion model satisfying all known cosmological constraints.

Although the standard {\Lambda}CDM model describes the cosmic microwave background radiation and the large scale structure of the Universe with great success, it has some tensions with observations in the effective number of neutrino species (dark radiation) and the number of small scale structures (overabundance problem). Here we propose a scenario which can relax these tensions by producing both dark matter and dark radiation by late decays of heavy particle. Thanks to the generation mechanism, dark matters are rather warm so that the small-scale structure problem is resolved. This scenario can be naturally realized in supersymmetric axion model, in which axions produced by saxion decays provide dark radiation, while axinos from saxion decays form warm dark matter. We identify a parameter region of supersymmetric axion model satisfying all known cosmological constraints.

The almost conformal dynamics of walking technicolor (TC) implies the existence of the approximate scale invariance, which breaks down spontaneously by the condensation of anti-techni and techni-fermions. According to the Goldstone theorem, a spinless, parity-even particle, called techni-dilaton (TD), then emerges at low energy. If TC exhibits an extreme walking, TD mass is parametrically much smaller than that of techni-fermions (around 1 TeV), while its decay constant is comparable to the cutoff scale of walking TC. We analyze the light, decoupled TD as a dark matter candidate and study cosmological productions of TD, both thermal and non-thermal, in the early Universe. The thermal population is governed dominantly by single TD production processes involving vertices breaking the scale symmetry, while the non-thermal population is by the vacuum misalignment and is accumulated via harmonic and coherent oscillations of misaligned classical TD fields. The non-thermal population turns out to be dominant and large enough to explain the abundance of presently observed dark matter, while the thermal population is highly suppressed due to the large TD decay constant. Several cosmological and astrophysical limits on the light, decoupled TD are examined to find that the TD mass is constrained to be in a range between 0.01 eV and 500 eV. From the combined constraints on cosmological productions and astrophysical observations, we find that the light, decoupled TD can be a good dark matter candidate with the mass around a few hundreds of eV for typical models of (extreme) walking TC. We finally mention possible designated experiments to detect the TD dark matter.

We show that a Dirac right-handed scalar neutrino can be weakly interacting massive particle in the neutrinophilic Higgs model. When the additional Higgs fields couple only to the leptonic sector through neutrino Yukawa couplings, the right number of relic density of dark matter can be obtained from thermal freeze-out of the dark matter annihilation into charged leptons and neutrinos. At present, this annihilation is suppressed by the velocity of dark matter. However one-loop annihilation cross section into $\gamma\gamma$ can be larger than that of the helicity suppressed annihilation into fermions, because relevant coupling constants are different. Hence, gamma-ray line signal which might have been observed in the Fermi-LAT is also able to be explained by its annihilation.

There might be a light scalar field during inflation which is not responsible for the accelerating inflationary expansion. Then, its quantum fluctuation is stretched during inflation. This scalar field could be a curvaton, if it decays at a late time. In addition, if the inflaton decay rate depends on the light scalar field expectation value by interactions between them, density perturbations could be generated by the quantum fluctuation of the light field when the inflaton decays. This is modulated reheating mechanism. We study curvature perturbation in models where a light scalar field does not only play a role of curvaton but also induce modulated reheating at the inflaton decay. We calculate the non-linearity parameters as well as the scalar spectral index and the tensor-to-scalar ratio. We find that there is a parameter region where non-linearity parameters are also significantly enhanced by the cancellation between the modulated effect and the curvaton contribution. For the simple quadratic potential model of both inflaton and curvaton, both tensor-to-scalar ratio and nonlinearity parameters could be simultaneously large.

Axino arises in supersymmetric versions of axion models and is a natural candidate for cold or warm dark matter. Here we revisit axino dark matter produced thermally and non-thermally in light of recent developments. First we discuss the definition of axino relative to low energy axion one for several KSVZ and DFSZ models of the axion. Then we review and refine the computation of the dominant QCD production in order to avoid unphysical cross-sections and, depending on the model, to include production via SU(2) and U(1) interactions and Yukawa couplings.

We study the generation of the primordial curvature perturbation in multi-field inflation. Considering both the evolution of the perturbation during inflation and the effects generated at the end of inflation, we present a general formula for the curvature perturbation. We provide the analytic expressions of the power spectrum, spectral tilt and non-Gaussianity for the separable potentials of two inflaton scalars, and apply them to some specific models.

We propose an inflationary model ("natural hybrid model"), which combines the supersymmetric hybrid model and the natural inflation model to achieve the spectral index of 0.96, and the axion decay constant smaller than the Planck scale, f<< M_P. By introducing both U(1)_R and a shift symmetry and employing the minimal Kahler potential, the eta-problem can be still avoided. The two inflaton fields in this model can admit large non-Gaussianity.

We discuss the effective interactions of axion supermultiplet, which might be important for analyzing the cosmological aspect of supersymmetric axion model. Related to axino cosmology, it is stressed that three seemingly similar but basically different quantities, the Wilsonian axino-gluino-gluon coupling, the 1PI axino-gluino-gluon amplitude, and the PQ anomaly coefficient, should be carefully distinguished from each other for correct analysis of the thermal production of axinos in the early Universe. It is then noticed that the 1PI axino-gluino-gluon amplitude at energy scale p in the range M_\Phi < p < v_{PQ} is suppressed by (M_\Phi ^2)/(p^2) in addition to the well-known suppression by p/(16\pi^2 v_{PQ}), where M_\Phi is the mass of the heaviest PQ-charged and gauge-charged matter supermultiplet in the model, which can be well below the PQ scale v_{PQ}. As a result, axino production at temperature T>M_\Phi is dominated by the production by matter supermultiplet, not by the production by gauge supermultiplet. Still the axino production rate is greatly reduced if M_\Phi << v_{PQ}, which would make the subsequent cosmology significantly altered. This would be most notable in the supersymmetric DFSZ model in which M_\Phi corresponds to the Higgsino mass which is around the weak scale, however a similar reduction is possible in the KSVZ model also. We evaluate the relic axino density for both the DFSZ and KSVZ models while including the axino production in the processes involving the heaviest PQ-charged and gauge-charged matter supermultiplet.

In this paper we analyse three models of the early universe, for which the respective mechanisms for generating the curvature perturbation are considered disparate. We find that in fact the mechanisms are very similar, and hence explain why they give rise to a large non-gaussianity. We show that the mechanism for generating the primordial curvature perturbation, and hence the observable non-gaussianity, is similar in both the Curvaton and Modulated Reheating models. In both cases the model can be written in terms of an energy transfer between the constituting fluids. We then show that this is also true for the mechanism of generating the curvature perturbation by symmetry breaking the end of inflation. We then relate this to the non-gaussian contribution to the curvature perturbation and find that it is inversely proportional to the efficiency with which the curvature perturbation is transferred between the fluids. For the first time, we generalise models of modulated reheating to allow for a non-linear energy transfer rate.

We study the density perturbation by a curvaton with a double well potential and estimate the nonlinear parameters for non-Gaussianity and the amplitude of gravitational wave background generated during inflation. The predicted nonlinear parameters strongly depend on the size of a curvaton self-coupling constant as well as the reheating temperature after inflation for a given initial amplitude of the curvaton. The difference from usual massive self-interacting curvaton is also emphasized.

We review models which generate a large non-Gaussianity of the local form. We first briefly consider three models which generate the non-Gaussianity either at or after the end of inflation; the curvaton scenario, modulated (p)reheating and an inhomogeneous end of inflation. We then focus on ways of generating the non-Gaussianity during inflation. We derive general conditions which a product or sum separable potential must satisfy in order to generate a large local bispectrum during slow-roll inflation. As an application we consider two-field hybrid inflation. We then derive a formalism not based on slow roll which can be applied to models in which the slow-roll parameters become large before inflation ends. An exactly soluble two-field model is given in which this happens. Finally we also consider further non-Gaussian observables; a scale dependence of f_NL and the trispectrum.

In this work, we explore the "degenerate gravitino" scenario where the mass difference between the gravitino and the lightest MSSM particle is much smaller than the gravitino mass itself. In this case, the energy released in the decay of the next to lightest sypersymmetric particle (NLSP) is reduced. Consequently the cosmological and astrophysical constraints on the gravitino abundance, and hence on the reheating temperature, become softer than in the usual case. On the other hand, such small mass splittings generically imply a much longer lifetime for the NLSP. We find that, in the constrained MSSM (CMSSM), for neutralino LSP or NLSP, reheating temperatures compatible with thermal leptogenesis are reached for small splittings of order 10^{-2} GeV. While for stau NLSP, temperatures of 4x10^9 GeV can be obtained even for splittings of order of tens of GeVs. This "degenerate gravitino" scenario offers a possible way out to the gravitino problem for thermal leptogenesis in supersymmetric theories.

We consider the three-body decays of gravitino dark matter in supersymmetric scenarios with bilinear R-parity violation. In particular, gravitino decays into lepton+W^* (lepton f\bar f') and \nu+Z^* (\nu f\bar f) are examined for gravitino masses below Mw. After computing the gravitino decay rates into these three-body final states and studying their dependence on supersymmetric parameters, we find that these new decay modes are often more important than the two-body decay, into a photon and a neutrino, considered in previous works. Consequently, the gravitino lifetime and its branching ratios are substantially modified, with important implications for the indirect detection of gravitino dark matter.

The $\mu\nu$SSM provides a solution to the $\mu$-problem of the MSSM and explains the origin of neutrino masses by simply using right-handed neutrino superfields. Given that R-parity is broken in this model, the gravitino is a natural candidate for dark matter since its lifetime becomes much longer than the age of the Universe. We consider the implications of gravitino dark matter in the $\mu\nu$SSM, analyzing in particular the prospects for detecting gamma rays from decaying gravitinos. If the gravitino explains the whole dark matter component, a gravitino mass larger than 20 GeV is disfavored by the isotropic diffuse photon background measurements. On the other hand, a gravitino with a mass range between 0.1-20 GeV gives rise to a signal that might be observed by the FERMI satellite. In this way important regions of the parameter space of the $\mu\nu$SSM can be checked.

We re-consider the gravitino as dark matter in the framework of the Constrained MSSM. We include several recently suggested improvements on: (i) the thermal production of gravitinos, (ii) the calculation of the hadronic spectrum from NLSP decay and (iii) the BBN calculation including stau bound-state effects. In most cases we find an upper bound on the reheating temperature $\treh\lsim {\rm a few}\times 10^7\gev$ from over-production of $^6{\rm Li}$ from bound state effects. We also find an upper limit on the stau lifetime of $3\times10^4 \sec$, which is nearly an order of magnitude larger than the simple limit $5\times 10^3 \sec$ often used to avoid the effect of bound-state catalysis. The bound on $\treh$ is relaxed to $\lsim {\rm a few}\times 10^8\gev$ when we use a more conservative bound on ${^6}{\rm Li}/{^7}{\rm Li}$, in which case a new region at small stau mass at $\sim 100\gev$ and much longer lifetimes opens up. Such a low stau mass region can be easily tested at the LHC.

We study the generation of non-Gaussianity in models of hybrid inflation with two inflaton fields, (2-brid inflation). We analyse the region in the parameter and the initial condition space where a large non-Gaussianity may be generated during slow-roll inflation which is generally characterised by a large f_NL, tau_NL and a small g_NL. For certain parameter values we can satisfy tau_NL>>f_NL^2. The bispectrum is of the local type but may have a significant scale dependence. We show that the loop corrections to the power spectrum and bispectrum are suppressed during inflation, if one assume that the fields follow a classical background trajectory. We also include the effect of the waterfall field, which can lead to a significant change in the observables after the waterfall field is destabilised, depending on the couplings between the waterfall and inflaton fields.

We study the evolution of the curvature perturbation on the super-horizon scales starting from the inflationary epoch until there remains only a single dynamical degree of freedom, presureless matter, in the universe. We consider the cosmic inflation driven by a multiple number of the inflaton fields, which decay into both radiation and pressureless matter components. We present a complete set of the exact background and perturbation equations which describe the evolution of the universe throughout its history. By applying these equations to the simple but reasonable model of multi-field chaotic inflation, we explicitly show that the total curvature perturbation is continuously varying because of the non-adiabatic components of the curvature perturbation generated by the multiple inflaton fields throughout the whole evolution of the universe. We also provide an useful analytic estimation of the total as well as matter and radiation curvature perturbations, assuming that matter is completely decoupled from radiation from the beginning. The resulting isocurvature perturbation between matter and radiation is at most sub-percent level when the masses of the inflaton fields are distributed between $10^{-6}m_{pl}$ and $10^{-5}m_{pl}$. We find that this result is robust unless we use non-trivial decay rates, and that thus, in general, it is hard to obtain large matter-radiation isocurvature perturbation. Also, by using the $\delta{N}$ formalism, we point out that the inflationary calculation, especially when involving multiple inflaton fields, is likely to lose the potentially important post-inflationary evolution which can modify the resulting curvature perturbation.

Standard Big Bang nucleosynthesis at baryonic density as inferred by WMAP implies a primordial Li7 abundance factor of two to three larger than that inferred by observations of low--metallicity halo stars. Recent observations of Li6 in halo stars suggest a pre--galactic origin of this isotope, which is exceedingly difficult to explain by putative high redshift cosmic ray populations. We investigate if one or both of these lithium problems may be solved by late decaying relic particles in the Constrained Minimal Supersymmetric Standard Model (CMSSM) coupled to gravity. Assuming that the gravitino is the lightest supersymmetric particle (LSP) we find that in large parts of the CMSSM parameter space either of these problems can be potentially solved. In particular, we also find solutions where both lithium problems may be solved simultaneously. These entail the hadronic decay of relic 1 TeV staus into 100 GeV gravitinos at 1000 sec after the Big Bang, as proposed by one of us before. Moreover, the produced gravitinos naturally contribute a large fraction, or all, to the by WMAP required dark matter density. A prediction of this model is the dark matter to be lukewarm with free--streaming lengths of the order of a 4 keV early freezing--out relic particle. Such lukewarm gravitino dark matter may have beneficial effects to the formation of galactic structures and may be constrained by the reionisation history of the Universe. The current version of the paper presents results for both cases (a) when catalytic nuclear reations are included and (b) when they are neglected.

We study the level of primordial non-Gaussianity in slow-roll two-field inflation. Using an analytic formula for the nonlinear parameter f_nl in the case of a sum or product separable potential, we find that it is possible to generate significant non-Gaussianity even during slow-roll inflation with Gaussian perturbations at Hubble exit. In this paper we give the general conditions to obtain large non-Gaussianity and calculate the level of fine-tuning required to obtain this. We present explicit models in which the non-Gaussianity at the end of inflation can exceed the current observational bound of |f_nl|<100.

We consider cosmological consequences of a heavy axino, decaying to the neutralino in R-parity conserving models. The importance and influence of the axino decay on the resultant abundance of neutralino dark matter depends on the lifetime and the energy density of axino. For a high reheating temperature after inflation, copiously produced axinos dominate the energy density of the universe and its decay produces a large amount of entropy. As a bonus, we obtain that the upper bound on the reheating temperature after inflation via gravitino decay can be moderated, because the entropy production by the axino decay more or less dilutes the gravitinos.

After a period of inflationary expansion, the Universe reheated and reached full thermal equilibrium at the reheating temperature T_R. In this work we point out that, in the context of effective low-energy supersymmetric models, LHC measurements may allow one to determine T_R as a function of the mass of the dark matter particle assumed to be either an axino or a gravitino. An upper bound on their mass may also be derived.

In the framework of the Constrained MSSM we re--examine the gravitino as the lightest superpartner and a candidate for cold dark matter in the Universe. Unlike in other recent studies, we include both a thermal contribution to its relic population from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. Relative to a previous analysis [1] we update, extend and considerably improve our treatment of constraints from observed light element abundances on additional energy released during BBN in association with late gravitino production. Assuming the gravitino mass in the GeV to TeV range, and for natural ranges of other supersymmetric parameters, the neutralino region is excluded, while for smaller values of the gravitino mass it becomes allowed again. The gravitino relic abundance is consistent with observational constraints on cold dark matter from BBN and CMB in some well defined domains of the stau region but, in most cases, only due to a dominant contribution of the thermal population. This implies, depending on the gravitino mass, a large enough reheating temperature. If $\mgravitino>1$ GeV then $T_R>10^7$ GeV, if allowed by BBN and other constraints but, for light gravitinos, if $\mgravitino>100$ keV then $T_R>3\times 10^3$ GeV. On the other hand, constraints mostly from BBN imply an upper bound $T_R \lsim {a few}x 10^8\times10^9$ GeV which appears inconsistent with thermal leptogenesis. Finally, most of the preferred stau region corresponds to the physical vacuum being a false vacuum. The scenario can be partially probed at the LHC.

We study the evolution of the universe which contains a multiple number of non-relativistic scalar fields decaying into both radiation and pressureless matter. We present a powerful analytic formalism to calculate the matter and radiation curvature perturbations, and find that our analytic estimates agree with full numerical results within an error of less than one percent. Also we discuss the isocurvature perturbation between matter and radiation components, which may be detected by near future cosmological observations, and point out that it crucially depends on the branching ratio of the decay rate of the scalar fields and that it is hard to make any model independent predictions.

Theories beyond the standard model such as string theory motivate low energy effective field theories with several scalar fields which are not only coupled through a potential but also through their kinetic terms. For such theories we derive the general formulae for the running of the spectral indices for the adiabatic, isocurvature and correlation spectra in the case of two field inflation. We also compute the expected non-Gaussianity in such models for specific forms of the potentials. We find that the coupling has little impact on the level of non-Gaussianity during inflation.

In the framework of the CMSSM we study the gravitino as the lightest supersymmetric particle and the dominant component of cold dark matter in the Universe. We include both a thermal contribution to its relic abundance from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. In general both contributions can be important, although in different regions of the parameter space. We further include constraints from BBN on electromagnetic and hadronic showers, from the CMB blackbody spectrum and from collider and non--collider SUSY searches. The region where the neutralino is the next--to--lightest superpartner is severely constrained by a conservative bound from excessive electromagnetic showers and probably basically excluded by the bound from hadronic showers, while the stau case remains mostly allowed. In both regions the constraint from CMB is often important or even dominant. In the stau case, for the assumed reasonable ranges of soft SUSY breaking parameters, we find regions where the gravitino abundance is in agreement with the range inferred from CMB studies, provided that, in many cases, a reheating temperature $\treh$ is large, $\treh\sim10^{9}\gev$. On the other side, we find an upper bound $\treh\lsim 5\times 10^{9}\gev$. Less conservative bounds from BBN or an improvement in measuring the CMB spectrum would provide a dramatic squeeze on the whole scenario, in particular it would strongly disfavor the largest values of $\treh\sim 10^{9}\gev$. The regions favored by the gravitino dark matter scenario are very different from standard regions corresponding to the neutralino dark matter, and will be partly probed at the LHC.

A slow-rolling scalar field ($Q\equiv$ Quintessence) with potential energy $V_Q\sim (3\times 10^{-3} {\rm eV})^4$ has been proposed as the origin of accelerating universe at present. We investigate the effective potential of $Q$ in the framework of supergravity model including the quantum corrections induced by generic (nonrenormalizable) couplings of $Q$ to the gauge and charged matter multiplets. It is argued that the K\"ahler potential, superpotential and gauge kinetic functions of the underlying supergravity model are required to be invariant under the variation of $Q$ with an extremely fine accuracy in order to provide a working quintessence potential. Applying these results for string or $M$-theory, we point out that the heterotic $M$-theory or Type I string axion can be a plausible candidate for quintessence if (i) it does not couple to the instanton number of gauge interactions not weaker than those of the standard model and (ii) the modulus partner ${\rm Re}(Z)$ of the periodic quintessence axion ${\rm Im}(Z)\equiv {\rm Im}(Z)+1$ has a large VEV: ${\rm Re}(Z)\sim \frac{1}{2\pi}\ln(m_{3/2}^2 M_{Planck}^2/V_Q)$. It is stressed that such a large ${\rm Re}(Z)$ gives the gauge unification scale at around the phenomenologically favored value $3\times 10^{16}$ GeV. To provide an accelerating universe, the quintessence axion should be at near the top of its effective potential at present, which requires a severe fine tuning of the initial condition of $Q$ and $\dot{Q}$ in the early universe. We discuss a late time inflation scenario based on the modular and CP invariance of the moduli effective potential, yielding the required initial condition in a natural manner if the K\"ahler metric of the quintessence axion superfield receives a sizable nonperturbative contribution.

A slow-rolling scalar field (Quintessence) has been proposed as the origin of accelerating universe at present. We discuss some features of quintessence, particularly those related with the {\it flatness} of its potential. We distinguish two types of quintessences, a metrical quintessence which couples to the particle physics sector {\it only} through the mixing with the spacetime metric and a non-metrical one without such property. It is stressed that these two types of quintessences have quite different features in regard to the fine tuning of parameters required for the flat potential, as well as different implications for the quintessence-mediated long range force. The fine tuning argument strongly suggests that non-metrical quintessence should correspond to the Goldstone boson of an almost exact nonlinear global symmetry whose explicit breaking by quantum gravity effects is highly suppressed. It is briefly discussed whether the radion field in higher dimensional theories can be a candidate for metrical quintessence. We finally discuss the possibility that string/$M$ theory axion corresponds to a Goldstone-type quintessence, and argue that a certain combination of the heterotic $M$ theory or Type I string theory axions can be a good candidate for quintessence if some conditions on the moduli dynamics are satisfied.

We examine the effects of the axion coupling to $R\tilde{R}$ on the evolution of cosmic structures. It is shown that the evolutions of the scalar- and vector-type perturbations are not affected by this axion coupling. However the axion coupling causes an asymmetric evolution of the two polarization states of the tensor-type perturbation, which may lead to a sizable polarization asymmetry in the cosmological gravitational wave if inflation involves a period in which the axion coupling is important. The polarization asymmetry produced during inflation are conserved over the subsequent evolution as long as the scales remain in the large-scale limit, and thus this may lead to an observable trace in the cosmic microwave background radiation.

In string/M-theory with a large compactification radius, some axion-like moduli can be much lighter than the gravitino. Generic moduli in gauge-mediated supersymmetry breaking models also have a mass far below the weak scale. Motivated by these, we examine the cosmological implications of light moduli for the mass range from the weak scale to an extremely small scale of order 10^{-26} eV, and obtain an upper bound on the initial moduli misalignment for both cases with and without a late entropy production.

We examine in the context of supersymmetric models whether the usual cosmological upper bound on the axion decay constant can be relaxed by assuming a period of stronger QCD in the early universe. By evaluating the axion potential in the early universe and also taking into account the dilaton potential energy, it is argued that a stronger QCD is not useful for raising up the bound.

Hadronic axions with the decay constant $f_a\simeq 10^{6}$ GeV may fulfill all astrophysical and laboratory constraints discussed so far. In this paper, we reexamine the possibility of the hadronic axion window while taking into account the uncertainties of some parameters describing low energy axion dynamics. It is found that $f_a$ in the range from $3\times 10^5$ GeV to $3\times 10^6$ can not be excluded by existing arguments. We then examine the implication of this hadronic axion window for the big-bang nucleosynthesis (NS) by evaluating the energy density of thermal axions at the nucleosynthesis epoch. Our analysis yields $(\rho_a/\rho_{\nu})_{NS}=0.4\sim 0.5$ which exceeds slightly the current best bound $(\rho_a/\rho_{\nu})_{NS}\leq 0.3$.

We examine whether observable majoron emission in double beta decay can be compatible with the big-bang nucleosynthesis (NS) and the observed neutrino flux from SN1987A. It is found that the NS upper bound on $^4$He abundance implies that the majoron-neutrino Yukawa coupling constant $g\leq 9\times 10^{-6}$ and its maximal value is allowed only when the scalar quartic coupling constant $\lambda$ is extremely small, $\lambda\leq 100g^2$. It is also observed that, although quite less restrictive, SN1987A also provides independent constraints on coupling constants.