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.

In this paper we focus on understanding the physical processes that lead to stable or unstable ionization fronts (I-fronts) observed in simulations of moving black holes (BHs). The front instability may trigger bursts of gas accretion, rendering the BH significantly more luminous than at steady-state. We perform a series of idealized three dimensional radiation hydrodynamics simulations resolving the I-fronts around BHs of mass $M_\mathrm{BH}$ and velocity $v_\infty$ accreting from a medium of density $n_\mathrm{H}$. The I-front, with radius $R_\mathrm{I}$, transitions from D-type to R-type as the BH velocity becomes larger than a critical value $v_\mathrm{R}\sim 40\,\mathrm{km/s}$. The D-type front is preceded by a bow-shock of thickness $\Delta R_\mathrm{I}$ that decreases as $v_\infty$ approaches $v_\mathrm{R}$. We find that both D-type and R-type fronts can be unstable given the following two conditions: i) for D-type fronts the shell thickness must be $\Delta R_\mathrm{I}/R_\mathrm{I}<0.05$ (i.e., $v_\infty \gtrsim 20\,\mathrm{km/s}$.), while no similar restriction holds for R-type fronts; ii) the temperature jump across the I-front must be $T_\mathrm{II}/T_\mathrm{I}>3$. This second condition is satisfied if $T_\mathrm{I}<5000\,\mathrm{K}$ or if $n_\mathrm{H}\,M_\mathrm{BH} \gtrsim 10^6\,M_\odot\,\mathrm{cm^{-3}}$. Due to X-ray pre-heating typically $T_\mathrm{I} \sim 10^4\,\mathrm{K}$, unless the D-type shell is optically thick to X-rays, which also happens when $n_\mathrm{H}\,M_\mathrm{BH}$ is greater than a metallicity-dependent critical value. We thus conclude that I-fronts around BHs are unstable only for relatively massive BHs moving trough very dense molecular clouds. We briefly discuss the observational consequences of the X-ray luminosity bursts likely associated with this instability.

We study star cluster formation in various environments with different metallicities and column densities by performing a suite of three-dimensional radiation hydrodynamics simulations. We find that the photoionization feedback from massive stars controls the star formation efficiency (SFE) in a star-forming cloud, and its impact sensitively depends on the gas metallicity $Z$ and initial cloud surface density $\Sigma$. At $Z=1~Z_{\odot}$, for instance, the SFE increases from 0.02 to 0.3 with $\Sigma$ from $10$ to $300~{\rm M_{\odot} {\rm pc^{-2}}}$. In low-metallicity cases $10^{-2}- 10^{-1} Z_{\odot}$, star clusters form from atomic warm gases because the molecule formation time is not short enough with respect to the cooling or dynamical time. In addition, the whole cloud is disrupted more easily by expanding HII bubbles which have higher temperature owing to less efficient cooling. With smaller dust attenuation, the ionizing radiation feedback from nearby massive stars is stronger and terminate star formation in dense clumps. These effects result in inefficient star formation in low-metallicity environments: the SFE is reduced by a factor of $\sim 0.3$ for the cases with $Z=10^{-2}~Z_{\odot}$ regardless of $\Sigma$. Newborn star clusters are also gravitationally less bound. We further develop a new semi-analytical model that can reproduce the simulation results well, particularly the observed dependencies of the SFEs on the cloud surface densities and metallicities.

Coalescence of intermediate-mass black holes (IMBHs) as a result of the migration toward galactic centers via dynamical friction may contribute to the formation of supermassive BHs. Here we reinvestigate the gaseous dynamical friction, which was claimed to be inefficient with radiative feedback from BHs in literature, by performing 3D radiation-hydrodynamics simulations that solve the flow structure in the vicinity of BHs. We consider a $10^4~M_\odot$ BH moving at the velocity $V_{\rm flow}$ through the homogeneous medium with metallicity $Z$ in the range of $0-0.1~Z_\odot$ and density $n_{\infty}$. We show that, if $n_{\infty} \lesssim 10^{6}~{\rm cm^{-3}}$ and $V_{\rm flow} \lesssim 60~{\rm km~s^{-1}}$, the BH is accelerated forward because of the gravitational pull from a dense shell ahead of an ionized bubble around the BH, regardless of the value of $Z$. If $n_{\infty} \gtrsim 10^{6}~{\rm cm^{-3}}$, however, our simulation shows the opposite result. The ionized bubble and associating shell temporarily appear, but immediately go downstream with significant ram pressure of the flow. They eventually converge into a massive downstream wake, which gravitationally drags the BH backward. The BH decelerates over the timescale of $\sim 0.01$~Myr, much shorter than the dynamical timescale in galactic disks. Our results suggest that IMBHs that encounter the dense clouds rapidly migrate toward galactic centers, where they possibly coalescence with others.

The super-Eddington accretion onto intermediate seed BHs is a potential formation mode of supermassive black holes exceeding $10^9~M_\odot$ in the early universe. We here investigate how such rapid accretion may occur with finite amounts of heavy elements contained in the gas and dust. In our 1D radiation-hydrodynamics simulations, the radiative transfer is solved for both the direct UV lights emitted by an accretion disk and the diffuse IR lights thermally emitted by dust grains. Our results show that the radiative force by the IR lights causes a strong feedback to regulate the mass accretion. The resulting mean accretion rate is lower with the higher metallicity, and there is the critical metallicity $Z \sim 10^{-2}~Z_\odot$, above which the super-Eddington accretion is prevented by the radiation pressure of the IR lights. With this taken into account, we examine if the dusty super-Eddington accretion occurs in young galaxies using a simple model. We show that a sufficient number of galaxies at $z \gtrsim 10$ can be such potential sites if BHs accrete the cold dense gas with $T \sim 10^2$ K, approximately the thermal equilibrium value at $Z = 10^{-2}~Z_\odot$. We argue that the efficiency of the BH growth via the rapid accretion depends on the metallicity, and that the metallicity slightly lower than $10^{-2}~Z_\odot$ provides a chance for the most efficient growth.

Supermassive stars (SMSs) with mass $\sim10^{5}~\rm{M}_{\odot}$ are promising candidates for the origin of supermassive black holes observed at redshift $\gtrsim6$. They are supposed to form as a result of rapid accretion of primordial gas, although it can be obstructed by the time variation caused by circum-stellar disc fragmentation due to gravitational instability. To assess the occurrence of fragmentation, we study the structure of marginally gravitationally unstable accretion discs, by using a steady one-dimensional thin disc model with detailed treatment of chemical and thermal processes. Motivated by two SMS formation scenarios, i.e., those with strong ultraviolet radiation background or with large velocity difference between the baryon and the dark matter, we consider two types of flows, i.e., atomic and molecular flows, respectively, for a wide range of the central stellar mass $10-10^5~\rm{M}_{\odot}$ and the accretion rate $10^{-3}-1~\rm{M}_{\odot}~\rm{yr}^{-1}$. In the case of a mostly atomic gas flowing to the disc outer boundary, the fragmentation condition is expressed as the accretion rate being higher than the critical value of $10^{-1}~\rm{M}_{\odot}~\rm{yr}^{-1}$ regardless of the central stellar mass. On the other hand, in the case of molecular flows, there is a critical disc radius outside of which the disc becomes unstable. Those conditions appears to be marginally satisfied according to numerical simulations, suggesting that disc fragmentation can be common during SMS formation.

Discovery of high-redshift (z > 6) supermassive black holes (BHs) may indicate that the rapid (or super-Eddington) gas accretion has aided their quick growth. Here, we study such rapid accretion of the primordial gas on to intermediate-mass (10^2 - 10^5 M_sun) BHs under anisotropic radiation feedback. We perform two-dimensional radiation hydrodynamics simulations that solve the flow structure across the Bondi radius, from far outside of the Bondi radius down to a central part which is larger than a circum-BH accretion disc. The radiation from the unresolved circum-BH disc is analytically modeled considering self-shadowing effect. We show that the flow settles into a steady state, where the flow structure consists of two distinct parts: (1) bipolar ionized outflowing regions, where the gas is pushed outward by thermal gas pressure and super-Eddington radiation pressure, and (2) an equatorial neutral inflowing region, where the gas falls toward the central BH without affected by radiation feedback. The resulting accretion rate is much higher than that in the case of isotropic radiation, far exceeding the Eddington-limited rate to reach a value slightly lower than the Bondi one. The opening angle of the equatorial inflowing region is determined by the luminosity and directional dependence of the central radiation. We find that photoevaporation from its surfaces set the critical opening angle of about ten degrees below which the accretion to the BH is quenched. We suggest that the shadowing effect allows even stellar-remnant BHs to grow rapidly enough to become high-redshift supermassive BHs.

Accretion on to seed black holes (BHs) is believed to play a crucial role in formation of supermassive BHs observed at high-redshift (z>6). Here, we investigate the combined effect of gas angular momentum and radiation feedback on the accretion flow, by performing 2D axially symmetric radiation hydrodynamics simulations that solve the flow structure across the Bondi radius and the outer part of the accretion disc simultaneously. The accreting gas with finite angular momentum forms a rotationally-supported disc inside the Bondi radius, where the accretion proceeds by the angular momentum transport due to assumed alpha-type viscosity. We find that the interplay of radiation and angular momentum significantly suppresses accretion even if the radiative feedback is weakened in an equatorial shadowing region. The accretion rate is O(alpha)\sim O(0.01-0.1) times the Bondi value, where alpha is the viscosity parameter. By developing an analytical model, we show that such a great reduction of the accretion rate persists unless the angular momentum is so small that the corresponding centrifugal radius is \lesssim 0.04 times the Bondi radius. We argue that BHs are hard to grow quickly via rapid mass accretion considering the angular momentum barrier presented in this paper.

Understanding Ly$\alpha$ emitting galaxies (LAEs) can be a key to reveal cosmic reionization and galaxy formation in the early Universe. Based on halo merger trees and Ly$\alpha$ radiation transfer calculations, we model redshift evolution of LAEs and their observational properties at $z \ge 6$. We consider ionized bubbles associated with individual LAEs and IGM transmission of Ly$\alpha$ photons. We find that Ly$\alpha$ luminosity tightly correlates with halo mass and stellar mass, while the relation with star formation rate has a large dispersion. Comparing our models with the observed luminosity function by Konno et al., we suggest that LAEs at $z \sim 7$ have galactic wind of $V_{\rm out} \gtrsim 100~\rm km\, s^{-1}$ and HI column density of $N_{\rm HI} \gtrsim 10^{20}~\rm cm^{-2}$. Number density of bright LAEs rapidly decreases as redshift increases, due to both lower star formation rate and smaller HII bubbles. Our model predicts future wide deep surveys with next generation telescopes, such as JWST, E-ELT and TMT, can detect LAEs at $z \sim 10$ with a number density of $n_{\rm LAE} \sim {\rm a~few~} \times10^{-6} ~\rm Mpc^{-3}$ for the flux sensitivity of $10^{-18} ~\rm erg\, cm^{-2}\, s^{-1}$. By combining these surveys with future 21-cm observations, it could be possible to detect both LAEs with $L_{Ly\alpha} \gtrsim 10^{42}~\rm erg~s^{-1}$ and their associated giant HII bubbles with the size $\gtrsim 250 ~\rm kpc$ at $z \sim 10$.

Fragmentation of filaments into dense cores is thought to be an important step in forming stars. The bar-mode instability of spherically collapsing cores found in previous linear analysis invokes a possibility of re-fragmentation of the cores due to their ellipsoidal (prolate or oblate) deformation. To investigate this possibility, here we perform three-dimensional self-gravitational hydrodynamics simulations that follow all the way from filament fragmentation to subsequent core collapse. We assume the gas is polytropic with index \gamma, which determines the stability of the bar-mode. For the case that the fragmentation of isolated hydrostatic filaments is triggered by the most unstable fragmentation mode, we find the bar mode grows as collapse proceeds if \gamma < 1.1, in agreement with the linear analysis. However, it takes more than ten orders-of-magnitude increase in the central density for the distortion to become non-linear. In addition to this fiducial case, we also study non-fiducial ones such as the fragmentation is triggered by a fragmentation mode with a longer wavelength and it occurs during radial collapse of filaments and find the distortion rapidly grows. In most of astrophysical applications, the effective polytropic index of collapsing gas exceeds 1.1 before ten orders-of-magnitude increase in the central density. Thus, supposing the fiducial case of filament fragmentation, re-fragmentation of dense cores would not be likely and their final mass would be determined when the filaments fragment.

Evidence for dust around supermassive black holes (SMBHs) in the early Universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black-holes (IMBHs) of $\sim 10^{5}~M_{\odot}$ in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the accretion of dusty gas onto IMBHs proceeds gently with small fluctuations of the accretion rate, whereas that of pristine gas causes more violent periodic bursts. At dust-to-gas mass ratios similar to the solar neighborhood, the time averaged luminosity becomes smaller than that for primordial gas by one order of magnitude and the time-averaged Eddington ratio ranges from $\sim 10^{-4}$ to $\sim 10^{-2}$ in clouds with initial gas densities of $n_{\rm H} = 10 - 1000~\rm cm^{-3}$. Our calculations show that the effect of dust opacity alone is secondary compared to the radiation pressure on dust in regulating the BH growth. We also derive spectral energy distributions at IR bands by calculating dust thermal emission and show that the flux ratio between $\lambda \lesssim 20~\rm \mu m$ and $\gtrsim 100~\rm \mu m$ is closely related to the Eddington ratio. Thermal emission from hot dust near the BH dominates only during the high accretion phase, producing higher flux density at $\lesssim 20~\rm \mu m$. Therefore, we suggest that the combinations of MIR observations by JWST and FIR observation by ALMA or Spitzer can be used to estimate the Eddington ratio of massive BHs.

We investigate the role of the H_2^+ channel on H_2 molecule formation during the collapse of primordial gas clouds immersed in strong radiation fields which are assumed to have the shape of a diluted black-body spectra with temperature T_rad. Since the photodissociation rate of H_2^+ depends on its level population, we take full account of the vibrationally-resolved H_2^+ kinetics. We find that in clouds under soft but intense radiation fields with spectral temperature T_rad < 7000 K, the H_2^+ channel is the dominant H_2 formation process. On the other hand, for harder spectra with T_rad > 7000 K, the H^- channel takes over H_2^+ in the production of molecular hydrogen. We calculate the critical radiation intensity needed for supermassive star formation by direct collapse and examine its dependence on the H_2^+ level population. Under the assumption of local thermodynamic equilibrium (LTE) level population, the critical intensity is underestimated by a factor of a few for soft spectra with T_rad < 7000 K. For harder spectra, the value of the critical intensity is not affected by the level population of H_2^+. This result justifies previous estimates of the critical intensity assuming LTE populations since radiation sources like young and/or metal-poor galaxies are predicted to have rather hard spectra.

It has been proposed that supermassive black holes (SMBHs) are originated from direct-collapse black holes (DCBHs) that are formed at z gtrsim 10 in the primordial gas in the case that H2 cooling is suppressed by strong external radiation. In this work, we study the critical specific intensity J^crit required for DCBH formation for various radiation spectral shapes by a series of one-zone calculations of a collapsing primordial- gas cloud. We calculate the critical specific intensity at the Lyman-Werner (LW) bands J^crit_LW,21 (in units of 10^-21 erg s^-1 Hz^-1 sr^-1 cm^-2) for realistic spectra of metal-poor galaxies. We find J^crit is not sensitive to the age or metallicity for the constant star formation galaxies with J^crit_LW,21 = 1300-1400, while J^crit decreases as galaxies become older or more metal-enriched for the instantaneous starburst galaxies. However, such dependence for the instantaneous starburst galaxies is weak for the young or extremely metal-poor galaxies: J^crit_LW,21 = 1000-1400 for the young galaxies and J^crit_LW,21 approx 1400 for the extremely metal-poor galaxies. The typical value of J^crit for the realistic spectra is higher than those expected in the literature, which affects the estimated DCBH number density n_DCBH. By extrapolating the result of Dijkstra, Ferrara and Mesinger, we obtain n_DCBH sim 10^-10 cMpc^-3 at z = 10, although there is still large uncertainty in this estimation.

We calculate the bispectrum of primordial curvature perturbations, \zeta, generated during "open inflation." Inflation occurs inside a bubble nucleated via quantum tunneling from the background false vacuum state. Our universe lives inside the bubble, which can be described as a Friedman-Lema\^itre-Robertson-Walker (FLRW) universe with negative spatial curvature, undergoing slow-roll inflation. We pay special attention to the issue of an initial state for quantum fluctuations. A "vacuum state" defined by a positive-frequency mode in de Sitter space charted by open coordinates is different from the Euclidean vacuum (which is equivalent to the so-called "Bunch-Davies vacuum"). Quantum tunneling then modifies the initial state away from the original Euclidean vacuum. While most of the previous study on modifications of the initial quantum state introduces, by hand, an initial time at which the quantum state is modified as well as the form of the modification, an effective initial time naturally emerges and the form is fixed by quantum tunneling in open inflation models. Therefore, open inflation enables a self-consistent computation of the effect of a modified initial state on the bispectrum. We find a term which goes as <\zeta_{k_1} \zeta_{k_2} \zeta_{k_3}> \propto 1/k_1^2k_3^4 in the so-called squeezed configurations, k_3\ll k_1\approx k_2, in agreement with the previous study. The bispectrum in the exact folded limit, e.g., k_1=k_2+k_3, is also enhanced and remains finite. However, these terms are exponentially suppressed when the wavelength of \zeta is smaller than the curvature radius of the universe. The leading-order bispectrum is equal to the usual one from single-field slow-roll inflation; the terms specific for open inflation arise only in the sub-leading order when the wavelength of \zeta is smaller than the curvature radius.

We reformulate quantum tunneling in a multi-dimensional system where the tunneling sector is non-linearly coupled to oscillators. The WKB wave function is explicitly constructed under the assumption that the system was in the ground state before tunneling. We find that the quantum state after tunneling can be expressed in the language of the conventional in-in formalism. Some implications of the result to cosmology are discussed.

We point out a possible generation mechanism of non-Gaussian bubbles in the sky due to bubble nucleation in the early universe. We consider a curvaton scenario for inflation and assume that the curvaton field phi, whose energy density is subdominant during inflation but which is responsible for the curvature perturbation of the universe, is coupled to another field sigma which undergoes false vacuum decay through quantum tunneling. For this model, we compute the skewness of the curvaton fluctuations due to its interaction with sigma during tunneling, that is, on the background of an instanton solution that describes false vacuum decay. We find that the resulting skewness of the curvaton can become large in the spacetime region inside the bubble. We then compute the corresponding skewness in the statistical distribution of the cosmic microwave background (CMB) temperature fluctuations. We find a non-vanishing skewness in a bubble-shaped region in the sky. It can be large enough to be detected in the near future, and if detected it will bring us invaluable information about the physics in the early universe.

We consider a multi-field open inflation model, in which one of the fields dominates quantum tunneling from a false vacuum while the other field governs slow-roll inflation within the bubble nucleated from false vacuum decay. We call the former the tunneling field and the latter the inflaton field. In the limit of a negligible interaction between the two fields, the false vacuum decay is described by a Coleman-De Luccia instanton. Here we take into account the coupling between the two fields and construct explicitly a multi-field instanton for a simple quartic potential model. We also solve the evolution of the scalar fields within the bubble. We find our model realizes open inflation successfully. This is the first concrete, viable model of open inflation realized with a simple potential. We then study the effect of the multi-field dynamics on the false vacuum decay, specifically on the tunneling rate. We find the tunneling rate increases in general in comparison with the single field case, though the increase is small unless the inflaton affects the instanton solution substantially.