The impact of the streaming between baryons and dark matter on the first structures has been actively explored by recent studies. We investigate how the key results are affected by two popular approximations. One is to implement the streaming by accounting for only the relative motion while assuming "baryons trace dark matter" spatially at the initialization of simulation. This neglects the smoothing on the gas density taking place before the initialization. In our simulation initialized at $z_i=200$, it overestimates the gas density power spectrum up to 40\% at $k\approx10^2~h~\mbox{Mpc}^{-1}$ at $z=20$. Halo mass ($M_h$) and baryonic fraction in halos ($f_{b,h}$) are also overestimated, but the relation between the two remains unchanged. The other approximation tested is to artificially amplify the density/velocity fluctuations in the cosmic mean density to simulate the first minihalos that form in overdense regions. This gives a head start to the halo growth while the subsequent growth is similar to that in the mean density. The growth in a true overdense region, on the other hand, is accelerated gradually in time. For example, raising $\sigma_8$ by 50\% effectively transforms $z\rightarrow\sqrt{1.5}z$ in the halo mass growth history while in 2-$\sigma$ overdensity, the growth is accelerated by a constant in redshift: $z\rightarrow{z+4.8}$. As a result, halos are more grown in the former than in the latter before $z\approx27$ and vice versa after. The $f_{b,h}$-$M_h$ relation is unchanged in those cases as well, suggesting that the Pop III formation rate for a given $M_h$ is insensitive to the tested approximations.

The formation of circumstellar disks is investigated using three-dimensional resistive magnetohydrodynamic simulations, in which the initial prestellar cloud has a misaligned rotation axis with respect to the magnetic field. We examine the effects of (i) the initial angle difference between the global magnetic field and the cloud rotation axis ($\theta_0$) and (ii) the ratio of the thermal to gravitational energy ($\alpha_0$). We study $16$ models in total and calculate the cloud evolution until $\sim \! 5000$ yr after protostar formation. Our simulation results indicate that an initial non-zero $\theta_0$ ($> 0$) promotes the disk formation but tends to suppress the outflow driving, for models that are moderately gravitationally unstable, $\alpha_0 \lesssim 1$. In these models, a large-sized rotationally-supported disk forms and a weak outflow appears, in contrast to a smaller disk and strong outflow in the aligned case ($\theta_0 = 0$). Furthermore, we find that when the initial cloud is highly unstable with small $\alpha_0$, the initial angle difference $\theta_0$ does not significantly affect the disk formation and outflow driving.

We study the formation of massive Population III binary stars using a newly developed radiation hydrodynamics code with the adaptive mesh refinement and adaptive ray-tracing methods. We follow the evolution of a typical primordial star-forming cloud obtained from a cosmological hydrodynamics simulation. Several protostars form as a result of disk fragmentation and grow in mass by the gas accretion, which is finally quenched by the radiation feedback from the protostars. Our code enables us, for the first time, to consider the feedback by both the ionizing and dissociating radiation from the multiple protostars, which is essential for self-consistently determining their final masses. At the final step of the simulation, we observe a very wide ($\gtrsim 10^4\,\mathrm{au}$) binary stellar system consisting of $60$ and $70\,M_\odot$ stars. One of the member stars also has two smaller mass ($10\,M_\odot$) companion stars orbiting at $200$ and $800\,\mathrm{au}$, making up a mini-triplet system. Our results suggest that massive binary or multiple systems are common among Population III stars.

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.

Using resistive magnetohydrodynamics simulations, the propagation of protostellar jets, the formation of circumstellar discs and the configuration of magnetic fields are investigated from the prestellar cloud phase until $\sim$500 yr after protostar formation. As the initial state, we prepare magnetized rotating clouds, in which the rotation axis is misaligned with the global magnetic field by an angle $\theta_0$. We calculate the cloud evolution for nine models with different $\theta_0 (=$ 0, 5, 10, 30, 45, 60, 80, 85, 90$^\circ$). Our simulations show that there is no significant difference in the physical quantities of the protostellar jet, such as the mass and momentum, among the models except for the model with $\theta_0=90^\circ$. On the other hand, the directions of the jet, disc normal and magnetic field are never aligned with each other during the early phase of star formation except for the model with $\theta_0=0^\circ$. Even when the rotation axis of the prestellar cloud is slightly inclined to the global magnetic field, the directions of the jet, disc normal and local magnetic field differ considerably, and they randomly change over time. Our results indicate that it is very difficult to extract any information from the observations of the directions of the outflow, disc and magnetic field at the scale of $\lesssim 1000$ au. Thus, we cannot use such observations to derive any restrictions on the star formation process.

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.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

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.

Recent observations uncover various phenomena around the protostar such as misalignment between the outflow and magnetic field, precession of the jet, and time variability of the ejected clumps, whose origins are under debate. We perform a three-dimensional resistive magnetohydrodynamics simulation of the protostar formation in a star-forming core whose rotation axis is tilted at an angle $45^\circ$ with respect to the initial magnetic field, in which the protostar is resolved with a spatial resolution of $0.01\,$au. In low-dense outer region, the prestellar core contracts along the magnetic field lines due to the flux freezing. In high-dense inner region, on the other hand, the magnetic dissipation becomes efficient and weakens the magnetic effects when the gas number density exceeds about $10^{11}\,{\rm cm^{-3}}$. Then, the normal direction of the flattened disc is aligned with the angular momentum vector. The outflow, jet, and protostellar ejection are driven from different scales of the circumstellar disc and spout in different directions normal to the warped disc. These axes do not coincide with the global magnetic field direction and vary with time. This study demonstrates that a couple of misalignment natures reported by observations can be simultaneously reproduced only by assuming the star-forming core rotating around a different direction from the magnetic field.

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.

Elucidating the particle physics nature of dark matter (DM) is one of the great challenges in modern science. The current lack of any direct DM detections in the laboratory heightens the need for astrophysical constraints, extending the search to DM models beyond the popular weakly interacting massive particle (WIMP) scenario. We here apply the classical Rees-Ostriker-Silk cooling criterion for galaxy formation to models with ultralight axion DM, also known as fuzzy dark matter (FDM). The resulting constraints provide a heuristic framework for upcoming observations, and our approximate analysis motivates the need for future, self-consistent simulations of FDM structure formation. We use observational constraints for the DM hosts of ultra faint dwarf (UFD) galaxies in the Local Group, together with the redshift constraints for the onset of primordial star formation from the recent EDGES 21-cm cosmology measurement, to illustrate this approach. We find that the existing constraints are straightforward to reconcile with standard $\Lambda$CDM, but disfavour FDM axion masses below $\sim 10^{-21}\,{\rm eV}/c^2$. The future potential for harnessing astrophysical probes of DM particle physics is compelling.

We study gravitational collapse of low-metallicity gas clouds and the formation of protostars by three-dimensional hydrodynamic simulations. Grain growth, non-equilibrium chemistry, molecular cooling, and chemical heating are solved in a self-consistent manner for the first time. We employ the realistic initial conditions for the abundances of metal and dust, and the dust size distribution obtained from recent Population III supernova calculations. We also introduce the state-of-the-art particle splitting method based on the Voronoi tessellation and achieve an extremely high mass resolution of ~10^{-5} Msun (10 earth masses) in the central region. We follow the thermal evolution of several early clouds with various metallicities. We show that the condition for cloud fragmentation depends not only on the gas metallicity but also on the collapse timescale. In many cases, the cloud fragmentation is prevented by the chemical heating owing to molecular hydrogen formation even though dust cooling becomes effective. Meanwhile, in several cases, efficient OH and H2O cooling promotes the cloud elongation, and then cloud "filamentation" is driven by dust thermal emission as a precursor of eventual fragmentation. While the filament fragmentation is driven by rapid gas cooling with >10^{-5} Zsun, fragmentation occurs in a different manner by the self-gravity of a circumstellar disk with <10^{-5} Zsun. We use a semi-analytic model to estimate the number fraction of the clouds which undergo the filament fragmentation to be a several percents with 10^{-5}--10^{-4} Zsun. Overall, our simulations show a viable formation path of the recently discovered Galactic low-mass stars with extremely small metallicities.

The recent detection of the sky-averaged 21-cm cosmological signal indicates a stronger absorption than the maximum allowed value based on the standard model. One explanation for the required colder primordial gas is the energy transfer between the baryon and dark matter fluids due to non-gravitational scattering. Here, we explore the thermal evolution of primordial gas, collapsing to form Population III (Pop III) stars, when this energy transfer is included. Performing a series of one-zone calculations, we find that the evolution results in stars more massive than in the standard model, provided that the dark matter is described by the best-fit parameters inferred from the 21-cm observation. On the other hand, a significant part of the dark matter parameter space can be excluded by the requirement to form massive Pop III stars sufficiently early in cosmic history. Otherwise, the radiation background needed to bring about the strong Wouthuysen-Field coupling at z >~ 17, inferred to explain the 21-cm absorption feature, could not be builtup. Intriguingly, the independent constraint from the physics of first star formation at high densities points to a similarly narrow range in dark matter properties. This exploratory study has to be followed-up with self-consistent three-dimensional simulations for a more rigorous derivation.

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 perform a set of cosmological simulations of early structure formation with incorporating baryonic streaming motions. We present a case where a significantly elongated gas cloud with $\sim\!10^4\,$solar masses (${\rm M_\odot}$) is formed in a pre-galactic ($\sim\!10^7\,{\rm M_\odot}$) dark halo. The gas streaming into the halo compresses and heats the massive filamentary cloud to a temperature of $\sim\!10,000\,$Kelvin. The gas cloud cools rapidly by atomic hydrogen cooling, and then by molecular hydrogen cooling down to $\sim\!400\,$Kelvin. The rapid decrease of the temperature and hence of the Jeans mass triggers fragmentation of the filament to yield multiple gas clumps with a few hundred solar masses. We estimate the mass of the primordial star formed in each fragment by adopting an analytic model based on a large set of radiation hydrodynamics simulations of protostellar evolution. The resulting stellar masses are in the range of $\sim\!50$-$120\,{\rm M_\odot}$. The massive stars gravitationally attract each other and form a compact star cluster. We follow the dynamics of the star cluster using a hybrid $N$-body simulation. We show that massive star binaries are formed in a few million years through multi-body interactions at the cluster center. The eventual formation of the remnant black holes will leave a massive black hole binary, which can be a progenitor of strong gravitational wave sources similar to those recently detected by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO).

We investigate the rotation velocity of the first stars by modelling the angular momentum transfer in the primordial accretion disc.Assessing the impact of magnetic braking, we consider the transition in angular momentum transport mode at the Alfv$\acute{\rm e}$n radius, from the dynamically dominated free-fall accretion to the magnetically dominated solid-body one.The accreting protostar at the centre of the primordial star-forming cloud rotates with close to breakup speed in the case without magnetic fields.Considering a physically-motivated model for small-scale turbulent dynamo amplification, we find that stellar rotation speed quickly declines if a large fraction of the initial turbulent energy is converted to magnetic energy ($\gtrsim 0.14$). Alternatively, if the dynamo process were inefficient, for amplification due to flux-freezing, stars would become slow rotators if the pre-galactic magnetic field strength is above a critical value, $\simeq 10^{-8.2}$G, evaluated at a scale of $n_{\rm H} = 1 {\rm cm}^{-3}$, which is significantly higher than plausible cosmological seed values ($\sim 10^{-15}$G). Because of the rapid decline of the stellar rotational speed over a narrow range in model parameters, the first stars encounter a bimodal fate: rapid rotation at almost the breakup level, or the near absence of any rotation.

Metal enrichment by the first-generation (Pop III) stars is the very first step of the matter cycle in the structure formation and it is followed by the formation of extremely metal-poor (EMP) stars. To investigate the enrichment process by the Pop III stars, we carry out a series of numerical simulations including the feedback effects of photoionization and supernovae (SNe) of Pop III stars with a range of masses of minihaloes (MHs), M_halo , and Pop III stars, M_PopIII . We find that the metal-rich ejecta reaches neighbouring haloes and external enrichment (EE) occurs when the halo binding energy is sufficiently below the SN explosion energy, E_SN . The neighbouring haloes are only superficially enriched, and the metallicity of the clouds is [Fe/H] < -5. Otherwise, the SN ejecta falls back and recollapses to form enriched cloud, i.e. internal enrichment (IE) process takes place. In case that a Pop III star explodes as a core-collapse SNe (CCSNe), MHs undergo IE, and the metallicity in the recollapsing region is -5 < [Fe/H] < -3 in most cases. We conclude that IE from a single CCSN can explain the formation of EMP stars. For pair-instability SNe (PISNe), EE takes place for all relevant mass range of MHs, consistent with no observational sign of PISNe among EMP stars.

The origin of super-massive black holes in the early universe remains poorly understood.Gravitational collapse of a massive primordial gas cloud is a promising initial process,but theoretical studies have difficulty growing the black hole fast enough.We report numerical simulations of early black hole formation starting from realistic cosmological conditions.Supersonic gas motions left over from the Big Bang prevent early gas cloud formation until rapid gas condensation is triggered in a proto-galactic halo. A protostar is formed in the dense, turbulent gas cloud, and it grows by sporadic mass accretion until it acquires 34,000 solar masses.The massive star ends its life with a catastrophic collapse to leave a black hole -- a promising seed for the formation of a monstrous black hole.

The formation of the first stars in the high-redshift Universe is a sensitive probe of the small-scale, particle physics nature of dark matter (DM). We carry out cosmological simulations of primordial star formation in ultra-light, axion-like particle DM cosmology, with masses of $10^{-22}$ and $10^{-21}\,{\rm eV}$, with de Broglie wavelengths approaching galactic scales ($\sim$kpc). The onset of star formation is delayed, and shifted to more massive host structures. For the lightest DM particle mass explored here, first stars form at $z \sim 7$ in structures with $\sim 10^9\,{\rm M}_\odot$, compared to the standard minihalo environment within the $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology, where $z \sim 20 - 30$ and $\sim 10^5 - 10^6\,{\rm M}_\odot$. Despite this greatly altered DM host environment, the thermodynamic behaviour of the metal-free gas as it collapses into the DM potential well asymptotically approaches a very similar evolutionary track. Thus, the fragmentation properties are predicted to remain the same as in $\Lambda$CDM cosmology, implying a similar mass scale for the first stars. These results predict intense starbursts in the axion cosmologies, which may be amenable to observations with the {\it James Webb Space Telescope}.

We study the formation of massive black holes in the first star clusters. We first locate star-forming gas clouds in proto-galactic haloes of $\gtrsim \!10^7\,{\rm M}_{\odot}$ in cosmological hydrodynamics simulations and use them to generate the initial conditions for star clusters with masses of $\sim \!10^5\,{\rm M}_{\odot}$. We then perform a series of direct-tree hybrid $N$-body simulations to follow runaway stellar collisions in the dense star clusters. In all the cluster models except one, runaway collisions occur within a few million years, and the mass of the central, most massive star reaches $\sim \!400-1900\,{\rm M}_{\odot}$. Such very massive stars collapse to leave intermediate-mass black holes (IMBHs). The diversity of the final masses may be attributed to the differences in a few basic properties of the host haloes such as mass, central gas velocity dispersion, and mean gas density of the central core. Finally, we derive the IMBH mass to cluster mass ratios, and compare them with the observed black hole to bulge mass ratios in the present-day Universe.

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.

The initial mass function of the first, Population III (Pop III), stars plays a vital role in shaping galaxy formation and evolution in the early Universe. One key remaining issue is the final fate of secondary protostars formed in the accretion disc, specifically whether they merge or survive. We perform a suite of hydrodynamic simulations of the complex interplay between fragmentation, protostellar accretion, and merging inside dark matter minihaloes. Instead of the traditional sink particle method, we employ a stiff equation of state approach, so that we can more robustly ascertain the viscous transport inside the disc. The simulations show inside-out fragmentation because the gas collapses faster in the central region. Fragments migrate on the viscous timescale, over which angular momentum is lost, enabling them to move towards the disc centre, where merging with the primary protostar can occur. This process depends on the fragmentation scale, such that there is a maximum scale of $(1 - 5) \times 10^4$ au, inside which fragments can migrate to the primary protostar. Viscous transport is active until radiative feedback from the primary protostar destroys the accretion disc. The final mass spectrum and multiplicity thus crucially depends on the effect of viscosity in the disc. The entire disc is subjected to efficient viscous transport in the primordial case with viscous parameter $\alpha \le 1$. An important aspect of this question is the survival probability of Pop III binary systems, possible gravitational wave sources to be probed with the Advanced LIGO detectors.

Gravitational collapse of a massive primordial gas cloud is thought to be a promising path for the formation of supermassive blackholes in the early universe. We study conditions for the so-called direct collapse (DC) blackhole formation in a fully cosmological context. We combine a semianalytic model of early galaxy formation with halo merger trees constructed from dark matter $N$-body simulations. We locate a total of 68 possible DC sites in a volume of $20\;h^{-1}\;\mathrm{Mpc}$ on a side. We then perform hydrodynamics simulations for 42 selected halos to study in detail the evolution of the massive clouds within them. We find only two successful cases where the gas clouds rapidly collapse to form stars. In the other cases, gravitational collapse is prevented by the tidal force exerted by a nearby massive halo, which otherwise should serve as a radiation source necessary for DC. Ram pressure stripping disturbs the cloud approaching the source. In many cases, a DC halo and its nearby light source halo merge before the onset of cloud collapse. Only when the DC halo is assembled through major mergers, the gas density increases rapidly to trigger gravitational instability. Based on our cosmological simulations, we conclude that the event rate of DC is an order of magnitude smaller than reported in previous studies, although the absolute rate is still poorly constrained. It is necessary to follow the dynamical evolution of a DC cloud and its nearby halo(s) in order to determine the critical radiation flux for DC.

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.

Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reionization on it. Towards that end, we perform a set of fully-coupled, radiation-hydrodynamics simulations from cosmological initial conditions, extending the mass resolution of previous work to the scale of minihalos. Pre-reionization structure is evolved until a redshift $z_i$ at which the ionizing radiation from external sources arrives to sweep an R-type ionization front supersonically across the volume in a few Myr, until it is trapped on the surfaces of minihalos and converted to D-type, after which the minihalo gas is removed by photoevaporative winds. Small-scale density structures during this time lead to a high ($>$10) clumping factor for ionized gas, which hugely boosts the recombination rate until the structures are disrupted by the hydrodynamic feedback after $\sim 10-100~\rm{Myr}$. For incoming stellar radiation with intensity $J_{21}$ in a $200~h^{-1}~\rm{kpc}$ box with the mean density contrast $\bar\delta$, the number of extra recombinations per H atom, on top of what is expected from homogeneously distributed gas, is given by $0.32[J_{21}]^{0.12}[(1+z_i)/11]^{-1.7}[1+\bar\delta]^{2.5}$. In models in which most of the volume is ionized toward the end of reionization, this can add more than one recombination per H atom to the ionizing photon budget to achieve reionization.

We present coupled stellar evolution (SE) and 3D radiation-hydrodynamic (RHD) simulations of the evolution of primordial protostars, their immediate environment, and the dynamic accretion history under the influence of stellar ionizing and dissociating UV feedback. Our coupled SE-RHD calculations result in a wide diversity of final stellar masses covering 10 Msun $\lesssim M_* \lesssim$ 1000 Msun. The formation of very massive ($\gtrsim$ 250 Msun) stars is possible under weak UV feedback, whereas ordinary massive (a few x 10 Msun) stars form when UV feedback can efficiently halt the accretion. This may explain the peculiar abundance pattern of a Galactic metal-poor star recently reported by Aoki et al. (2014), possibly the observational signature of very massive precursor primordial stars. Weak UV feedback occurs in cases of variable accretion, in particular when repeated short accretion bursts temporarily exceed 0.01 Msun/yr, causing the protostar to inflate. In the bloated state, the protostar has low surface temperature and UV feedback is suppressed until the star eventually contracts, on a thermal adjustment timescale, to create an HII region. If the delay time between successive accretion bursts is sufficiently short, the protostar remains bloated for extended periods, initiating at most only short periods of UV feedback. Disk fragmentation does not necessarily reduce the final stellar mass. Quite the contrary, we find that disk fragmentation enhances episodic accretion as many fragments migrate inward and are accreted onto the star, thus allowing continued stellar mass growth under conditions of intermittent UV feedback. This trend becomes more prominent as we improve the resolution of our simulations. We argue that simulations with significantly higher resolution than reported previously are needed to derive accurate gas mass accretion rates onto primordial protostars.

While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10~Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the standard scale-invariant PPS is modified at small length scales as $P(k) \sim k^{m_{\rm s}}$ with $m_{\rm s} > 1$. We run a series of cosmological hydrodynamic simulations to examine the dependence of the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models with $m_{\rm s} > 1$, star-forming gas clouds are formed at $z > 100$, when formation of hydrogen molecules is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars formed in such "hot" clouds grow very rapidly by accretion to become extremely massive stars that may leave massive black holes with a few hundred solar-masses at $z > 100$. The shape of the PPS critically affects the properties and the formation epoch of the first generation of stars. Future experiments of the CMB polarization and the spectrum distortion may provide important information on the nature of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum.

We perform a large set of cosmological simulations of early structure formation and follow the formation and evolution of 1540 star-forming gas clouds to derive the mass distribution of primordial stars. The star formation in our cosmological simulations is characterized by two distinct populations, the so-called Population III.1 stars and primordial stars formed under the influence of far ultraviolet (FUV) radiation (Population III.2D stars). In this work, we determine the stellar masses by using the dependences on the physical properties of star-forming cloud and/or the external photodissociating intensity from nearby primordial stars, which are derived from the results of two-dimensional radiation hydrodynamic simulations of protostellar feedback. The characteristic mass of the Pop III stars is found to be a few hundred solar masses at z ~ 25, and it gradually shifts to lower masses with decreasing redshift. At high redshifts z > 20, about half of the star-forming gas clouds are exposed to intense FUV radiation and thus give birth to massive Pop III.2D stars. However, the local FUV radiation by nearby Pop III stars becomes weaker at lower redshifts, when typical Pop III stars have smaller masses and the mean physical separation between the stars becomes large owing to cosmic expansion. Therefore, at z < 20, a large fraction of the primordial gas clouds host Pop III.1 stars. At z =< 15, the Pop III.1 stars are formed in relatively cool gas clouds due to efficient radiative cooling by H_2 and HD molecules; such stars have masses of a few x 10 Msun. Since the stellar evolution and the final fate are determined by the stellar mass, Pop III stars formed at different epochs play different roles in the early universe.

n-DBI gravity explicitly breaks Lorentz invariance by the introduction of a unit time-like vector field, thereby giving rise to an extra (scalar) degree of freedom. We look for observational consequences of this mode in two setups. Firstly, we compute the parametrized post-Newtonian (PPN) expansion of the metric to first post-Newtonian order. Surprisingly, we find that the PPN parameters are exactly the same as in General Relativity (GR), and no preferred-frame effects are produced. In particular this means that n-DBI gravity is consistent with all GR solar system experimental tests. We discuss the origin of such degeneracy between n-DBI gravity and GR, and suggest it may also hold in higher post-Newtonian order. Secondly, we study gravitational scalar perturbations of a Friedmann-Robertson-Walker space-time with a cosmological constant $\Lambda \geq 0$. In the case of de Sitter space, we show that the scalar mode grows as the universe expands and, in contrast with a canonical scalar field coupled to GR, it does not freeze on super horizon scales.

Bosonic superweakly interacting massive particles (super-WIMPs) are a candidate for warm dark matter. With the absorption of such a boson by a xenon atom these dark matter candidates would deposit an energy equivalent to their rest mass in the detector. This is the first direct detection experiment exploring the vector super-WIMPs in the mass range between 40 and 120 keV. Using 165.9 days of data no significant excess above background was observed in the fiducial mass of 41 kg. The present limit for the vector super-WIMPs excludes the possibility that such particles constitute all of dark matter. The absence of a signal also provides the most stringent direct constraint on the coupling constant of pseudoscalar super-WIMPs to electrons. The unprecedented sensitivity was achieved exploiting the low background at a level $10^{-4}$ kg$^{-1}$keV$_{ee}^{-1}$day$^{-1}$ in the detector.

A search for inelastic scattering of Weakly Interacting Massive Particles (WIMPs) on the isotope $^{129}$Xe was done in data taken with the single phase liquid xenon detector XMASS at the Kamioka Observatory. Using a restricted volume containing 41 kg of LXe at the very center of our detector we observed no significant excess of events in 165.9 days of data. Our background reduction allowed us to derive our limits without explicitly subtracting the remaining events which are compatible with background expectations and derive for e.g. a 50 GeV WIMP an upper limit for its inelastic cross section on $^{129}$Xe nuclei of 3.2 pb at the 90% confidence level.

We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in each gas cloud and the subsequent protostellar evolution until the gas mass accretion onto the protostar is halted by stellar radiative feedback. The accretion rates differ significantly among the primordial gas clouds which largely determine the final stellar masses. For low accretion rates the growth of a protostar is self-regulated by radiative feedback effects and the final mass is limited to several tens of solar masses. At high accretion rates the protostar's outer envelope continues to expand and the effective surface temperature remains low; such protostars do not exert strong radiative feedback and can grow in excess to one hundred solar masses. The obtained wide mass range suggests that the first stars play a variety of roles in the early universe, by triggering both core-collapse supernovae and pair-instability supernovae as well as by leaving stellar mass black holes. We find certain correlations between the final stellar mass and the physical properties of the star-forming cloud. These correlations can be used to estimate the mass of the first star from the properties of the parent cloud or of the host halo, without following the detailed protostellar evolution.

XMASS, a low-background, large liquid-xenon detector, was used to search for solar axions that would be produced by bremsstrahlung and Compton effects in the Sun. With an exposure of 5.6ton days of liquid xenon, the model-independent limit on the coupling for mass $\ll$ 1keV is $|g_{aee}|< 5.4\times 10^{-11}$ (90% C.L.), which is a factor of two stronger than the existing experimental limit. The bounds on the axion masses for the DFSZ and KSVZ axion models are 1.9 and 250eV, respectively. In the mass range of 10-40keV, this study produced the most stringent limit, which is better than that previously derived from astrophysical arguments regarding the Sun to date.

A search for light dark matter using low-threshold data from the single phase liquid xenon scintillation detector XMASS, has been conducted. Using the entire 835 kg inner volume as target, the analysis threshold can be lowered to 0.3 keVee (electron-equivalent) to search for light dark matter. With low-threshold data corresponding to a 5591.4 kg$\cdot$day exposure of the detector and without discriminating between nuclear-recoil and electronic events, XMASS excludes part of the parameter space favored by other experiments.

The XMASS project aims to detect dark matter, pp and $^{7}$Be solar neutrinos, and neutrinoless double beta decay using ultra pure liquid xenon. The first phase of the XMASS experiment searches for dark matter. In this paper, we describe the XMASS detector in detail, including its configuration, data acquisition equipment and calibration system.

We present the first 3D simulations to include the effects of dark matter annihilation feedback during the collapse of primordial mini-halos. We begin our simulations from cosmological initial conditions and account for dark matter annihilation in our treatment of the chemical and thermal evolution of the gas. The dark matter is modelled using an analytical density profile that responds to changes in the peak gas density. We find that the gas can collapse to high densities despite the additional energy input from the dark matter. No objects supported purely by dark matter annihilation heating are formed in our simulations. However, we find that the dark matter annihilation heating has a large effect on the evolution of the gas following the formation of the first protostar. Previous simulations without dark matter annihilation found that protostellar discs around Population III stars rapidly fragmented, forming multiple protostars that underwent mergers or ejections. When dark matter annihilation is included, however, these discs become stable to radii of 1000 AU or more. In the cases where fragmentation does occur, it is a wide binary that is formed.

We study the thermal evolution of primordial star-forming gas clouds using three-dimensional cosmological simulations. We critically examine how assumptions and approximations made in calculating radiative cooling rates affect the dynamics of the collapsing gas clouds. We consider two important molecular hydrogen cooling processes that operate in a dense primordial gas; H_2 line cooling and continuum cooling by H_2 collision-induced emission. To calculate the optically thick cooling rates, we follow the Sobolev method for the former, whereas we perform ray-tracing for the latter. We also run the same set of simulations using simplified fitting functions for the net cooling rates. We compare the simulation results in detail. We show that the time- and direction-dependence of hydrodynamic quantities such as gas temperature and local velocity gradients significantly affects the optically thick cooling rates. Gravitational collapse of the cloud core is accelerated when the cooling rates are calculated by using the fitting functions. The structure and evolution of the central pre-stellar disk are also affected. We conclude that physically motivated implementations of radiative transfer are necessary to follow accurately the thermal and chemical evolution of a primordial gas to high densities.

n-DBI gravity is a gravitational theory which yields near de Sitter inflation spontaneously at the cost of breaking Lorentz invariance by a preferred choice of foliation. We show that this breakdown endows n-DBI gravity with one extra physical gravitational degree of freedom: a scalar graviton. Its existence is established by Dirac's theory of constrained systems. Firstly, studying scalar perturbations around Minkowski space-time, we show that there exists one scalar degree of freedom and identify it in terms of the metric perturbations. Then, a general analysis is made in the canonical formalism, using ADM variables. It is useful to introduce an auxiliary scalar field, which allows recasting n-DBI gravity in an Einstein-Hilbert form but in a Jordan frame. Identifying the constraints and their classes we confirm the existence of an extra degree of freedom in the full theory, besides the two usual tensorial modes of the graviton. We then argue that, unlike the case of (the original proposal for) Horava-Lifschitz gravity, there is no evidence that the extra degree of freedom originates pathologies, such as vanishing lapse, instabilities and strong self-coupling at low energy scales.

We propose a 3 + 1 dimensional model of gravity which results in inflation at early times, followed by radiation- and matter-dominated epochs and a subsequent acceleration at late times. Both the inflation and late time acceleration are nearly de Sitter with a large hierarchy between the effective cosmological constants. There is no scalar field agent of inflation, and the transition from the inflation to the radiation-dominated period is smooth. This model is designed so that it yields, at the cost of giving up on Lorentz invariance in the gravitational sector, the Dirac-Born-Infeld type conformal scalar theory when the universe is conformally flat. It, however, resembles Einstein's gravity with the Gibbons-Hawking-York boundary term in weakly curved space-times.

Primordial stars formed in the early universe are thought to be hosted by compact dark matter (DM) halos. If DM consists of Weakly Interacting Massive Particles (WIMPs), such stars may be powered by DM annihilation during the early phases of their evolutions. We study the pre-main sequence evolutions of the primordial star using a detailed stellar evolution code under the assumption that the annihilation of adiabatically contracted WIMPs DM within the star provides a sufficient energy to sustain the stellar equilibrium. We follow the evolution of accreting stars using several gas mass accretion rates derived from cosmological simulations. We show that the stellar mass becomes very large, up to 900 - 1000 M_sun when the star reaches the main-sequence phase for a reasonable set of model parameters such as DM particle mass and the annihilation cross section. During the dark star phase, the star expands over a thousand solar-radii, while the surface temperature remains below 10^4 K. The energy generated by nuclear reactions is not dominant during this phase. We also study models with different gas mass accretion rates and the DM particle masses. All our models for different DM particle masses pass the dark star phase. The final mass of the dark stars is essentially unchanged for DM mass of m_DM <= 10 GeV. Gravitational collapse of the massive dark stars will leave massive black holes with mass as large as 1000 M_sun in the early universe.

The origin of Big Bang singularity in 3+1 dimensions can be understood in an exact string theory background obtained by an analytic continuation of a cigar like geometry with a nontrivial dilaton. In a T-dual conformal field theory picture there exists a closed string tachyon potential which excises the singular space-time of a strongly coupled regime to ensure that a higher dimensional universe has no curvature singularity. However in 3+1 dimensions the universe exhibits all the pathology of a standard Big Bang cosmology. The emergence of a singularity now owes to a higher dimensional orbifold singularity which does not have a curvature singularity in higher dimensions, suggesting that close to the compactification scale an effective description of 3+1 dimensions breaks down and bouncing universe emerges in 5 and higher dimensions.

A proposal is made for a cosmological D3/D7 model with a constant magnetic flux along the D7 world-volume. It describes an N=2 gauge model with Fayet-Iliopoulos terms and the potential of the hybrid P-term inflation. The motion of the D3-brane towards D7 in a phase with spontaneously broken supersymmetry provides a period of slow-roll inflation in the de Sitter valley, the role of the inflaton being played by the distance between D3 and D7-branes. After tachyon condensation a supersymmetric ground state is formed: a D3/D7 bound state corresponding to an Abelian non-linear (non-commutative) instanton. In this model the existence of a non-vanishing cosmological constant is associated with the resolution of the instanton singularity. We discuss a possible embedding of this model into a compactified M-theory setup.