Source author record

Shingo Hirano

Shingo Hirano appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2026arXiv

Rapid emergence of overmassive black holes in the early Universe

The origin of supermassive black holes (SMBHs) remains a long-standing problem in astrophysics. Recent JWST observations reveal an unexpectedly abundant population of overmassive black holes at z>4-6, where the BH masses lie far above local scaling relations and not reproduced by current cosmological models. How such overmassive black holes form and rapidly grow within young galaxies has remained unclear. Here we present fully cosmological radiation-hydrodynamic simulations that, for the first time, self-consistently follow the birth, early growth, and emergent observable signatures of SMBHs in proto-cluster environments. We find that heavy seeds of order $10^6 M_\text{sun}$ naturally form, exceeding typical theoretical expectations by an order of magnitude. These seeds rapidly develop dense, optically thick disks whose strong electron scattering produces broad H$α$ emission comparable to that seen in little red dots (LRDs). Sustained super-Eddington accretion then drives fast growth to $\sim 3 \times 10^7 ~M_\text{sun}$ by $z \sim 8$. These results provide a unified physical scenario in which LRDs correspond to a short-lived, enshrouded phase of heavy-seed formation, naturally evolving into the overmassive quasars detected by JWST and ultimately the progenitors of today's SMBHs.

preprint2022arXiv

Exponentially amplified magnetic field eliminates disk fragmentation around the Population III protostar

One critical remaining issue to unclear the initial mass function of the first (Population III) stars is the final fate of secondary protostars formed in the accretion disk, specifically whether they merge or survive. We focus on the magnetic effects on the first star formation under the cosmological magnetic field. We perform a suite of ideal magnetohydrodynamic simulations until 1000 years after the first protostar formation. Instead of the sink particle technique, we employ a stiff equation of state approach to represent the magnetic field structure connecting to protostars. Ten years after the first protostar formation in the cloud initialized with $B_0 = 10^{-20}$ G at $n_0 = 10^4\,{\rm cm^{-3}}$, the magnetic field strength around protostars amplifies from pico- to kilo-gauss, which is the same strength as the present-day star. The magnetic field rapidly winds up since the gas in the vicinity of the protostar ($\leq\!10$ au) has undergone several tens orbital rotations in the first decade after protostar formation. As the mass accretion progresses, the vital magnetic field region extends outward, and the magnetic braking eliminates fragmentation of the disk that would form in the unmagnetized model. On the other hand, assuming a gas cloud with small angular momentum, this amplification might not work because the rotation would be slower. However, disk fragmentation would not occur in that case. We conclude that the exponential amplification of the cosmological magnetic field strength, about $10^{-18}$ G, eliminates disk fragmentation around the Population III protostars.

preprint2022arXiv

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.

preprint2021arXiv

Large-scale variation in reionization history caused by Baryon-dark matter streaming velocity

At cosmic recombination, there was supersonic relative motion between baryons and dark matter, which originated from the baryonic acoustic oscillations in the early universe. This motion has been considered to have a negligible impact on the late stage of cosmic reionization because the relative velocity quickly decreases. However, recent studies have suggested that the recombination in gas clouds smaller than the local Jeans mass ($\lesssim$ $10^8~M_\odot$) can affect the reionization history by boosting the number of ultraviolet photons required for ionizing the intergalactic medium. Motivated by this, we performed a series of radiation-hydrodynamic simulations to investigate whether the streaming motion can generate variation in the local reionization history by smoothing out clumpy small-scale structures and lowering the ionizing photon budget. We found that the streaming velocity can add a variation of $Δz_e$ $\sim$ $0.05$ $-$ $0.5$ in the end-of-reionization redshift, depending on the level of X-ray preheating and the time evolution of ionizing sources. The variation tends to be larger when the ionizing efficiency of galaxies decreases toward later times. Given the long spatial fluctuation scales of the streaming motion ($\gtrsim 100$ Mpc), it can help to explain the Ly$α$ opacity variation observed from quasars and leave large-scale imprints on the ionization field of the intergalactic medium during the reionization. The pre-reionization heating by X-ray sources is another critical factor that can suppress small-scale gas clumping and can diminish the variation in $z_e$ introduced by the streaming motion.

preprint2020arXiv

First Structure Formation under the Influence of Gas-Dark Matter Streaming Velocity and Density: Impact of the Baryons-trace-dark matter Approximation

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 by 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 $σ_8$ by 50\% effectively transforms $z\rightarrow\sqrt{1.5}z$ in the halo mass growth history while in 2-$σ$ overdensity, the growth is accelerated by a constant in redshift: $z\rightarrow{z+4.8}$. As a result, halos have grown more 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.

preprint2020arXiv

The Birth of a Massive First-Star Binary

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.

preprint2020arXiv

The Effect of Misalignment between Rotation Axis and Magnetic Field on Circumstellar Disk

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 ($θ_0$) and (ii) the ratio of the thermal to gravitational energy ($α_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 $θ_0$ ($> 0$) promotes the disk formation but tends to suppress the outflow driving, for models that are moderately gravitationally unstable, $α_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 ($θ_0 = 0$). Furthermore, we find that when the initial cloud is highly unstable with small $α_0$, the initial angle difference $θ_0$ does not significantly affect the disk formation and outflow driving.

preprint2019arXiv

Misalignment of Magnetic Fields, Outflows and Discs in Star-forming Clouds

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 $θ_0$. We calculate the cloud evolution for nine models with different $θ_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 $θ_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 $θ_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.

preprint2016arXiv

Cosmological Simulations of Early Blackhole Formation: Halo Mergers, Tidal Disruption, and the Conditions for Direct Collapse

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.

preprint2016arXiv

Formation of Massive Primordial Stars: Intermittent UV Feedback with Episodic Mass Accretion

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.

preprint2016arXiv

The Hydrodynamic Feedback of Cosmic Reionization on Small-Scale Structures and Its Impact on Photon Consumption during the Epoch of Reionization

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δ$, 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δ]^{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.

preprint2015arXiv

Early structure formation from primordial density fluctuations with a blue-tilted power spectrum

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.

preprint2015arXiv

Primordial Star Formation under the Influence of Far Ultraviolet Radiation: 1540 Cosmological Halos and the Stellar Mass Distribution

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.

preprint2014arXiv

One Hundred First Stars : Protostellar Evolution and the Final Masses

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.

preprint2012arXiv

Radiative cooling implementations in simulations of primordial star formation

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.

preprint2012arXiv

WIMP DM and first stars: suppression of fragmentation in primordial star formation

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.

preprint2011arXiv

Evolution of Primordial Stars Powered by Dark Matter Annihilation up to the Main-Sequence Stage

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.