Researcher profile

Kentaro Nagamine

Kentaro Nagamine contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2022arXiv

Dust diffusion in SPH simulations of an isolated galaxy

We compute the evolution of the grain size distribution (GSD) in a suite of numerical simulations of an isolated Milky-Way-like galaxy using the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4}. The full GSD is sampled on a logarithmically spaced grid with 30 bins, and its evolution is calculated self-consistently with the hydrodynamical and chemical evolution of the galaxy using a state-of-the-art star formation and feedback model. In previous versions of this model, the GSD tended to be slightly biased towards larger grains and the extinction curve had a tendency to be flatter than the observations. This work addresses these issues by considering the diffusion of dust and metals through turbulence on subgrid scales and introducing a multi-phase subgrid model that enables a smoother transition from diffuse to dense gas. We show that diffusion can significantly enhance the production of small grains and improve the agreement with the observed dust extinction curve in the Milky Way.

preprint2022arXiv

FOREVER22: galaxy formation in protocluster regions

We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization.

preprint2022arXiv

Osaka Feedback Model II: Modeling Supernova Feedback Based on High-Resolution Simulations

Feedback from supernovae (SNe) is an essential mechanism that self-regulates the growth of galaxies, and a better model of SN feedback is still needed in galaxy formation simulations. In the first part of this paper, using an Eulerian hydrodynamic code Athena++, we find universal scaling relations for the time evolution of momentum and radius for a superbubble, when the momentum and time are scaled by those at the shell-formation time. In the second part of this paper, we develop an SN feedback model based on the Athena++ simulation results utilizing Voronoi tessellation around each star particle, and implement it into the GADGET3-Osaka smoothed particle hydrodynamic code. Our feedback model was demonstrated to be isotropic and conservative in terms of energy and momentum. We examined the mass/energy/metal loading factors and find that our stochastic thermal feedback model produced galactic outflow that carries metals high above the galactic plane but with weak suppression of star formation. Additional mechanical feedback further suppressed star formation and brought the simulation results in better agreement with the observations of the Kennicutt--Schmidt relation, with all the results being within the uncertainties of observed data. We argue that both thermal and mechanical feedback are necessary for the SN feedback model of galaxy evolution when an individual SN bubble is unresolved.

preprint2022arXiv

The co-evolution of molecular hydrogen and the grain size distribution in an isolated galaxy

Understanding the evolution of dust and molecular hydrogen (H$_2$) is a critical aspect of galaxy evolution, as they affect star formation and the spectral energy distribution of galaxies. We use the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4} to compute the evolution of dust and H$_2$ in a suite of numerical simulations of an isolated Milky-Way-like galaxy. The evolution of the full grain size distribution (GSD) is solved by sampling the grain size on a logarithmically spaced grid with 30 bins. The evolution of a primordial chemistry network with twelve species is solved consistently with the hydrodynamic evolution of the system, including star formation, metal and energy ejections from stars into the interstellar medium through supernova feedback and stellar winds. The formation model for H$_2$ considers the GSD and photo-dissociation through the UV radiation of young stars. We identify the processes needed for producing a sizeable amount of H$_2$, verify that the resulting star formation law in the later stages of galaxy evolution is consistent with observations of local spirals, and show that our model manages to produce a galactic molecular gas fraction in line with observations of Milky-Way-like galaxies. We stress the importance of the co-evolution of the GSD and H$_2$, as models assuming a fixed MRN shape for the GSD overestimate the production of H$_2$ in regimes where the dust abundance is dominated by large grains and underestimate it in the regime where the dust is dominated by small grains, both of which are realized in simulations of dust evolution.

preprint2020arXiv

3D Distribution Map of HI Gas and Galaxies Around an Enormous Ly$α$ Nebula and Three QSOs at $z=2.3$ Revealed by the HI Tomographic Mapping Technique

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$α$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$α$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

preprint2020arXiv

Direct collapse to supermassive black hole seeds: the critical conditions for suppression of $\rm H_2$ cooling

Observations of high-redshift quasars imply the presence of supermassive black holes already at z~ 7.5. An appealing and promising pathway to their formation is the direct collapse scenario of a primordial gas in atomic-cooling haloes at z ~ 10 - 20, when the $\rm H_2$ formation is inhibited by a strong background radiation field, whose intensity exceeds a critical value, $J_{\rm crit}$. To estimate $J_{\rm crit}$, typically, studies have assumed idealized spectra, with a fixed ratio of $\rm H_{2}$ photo-dissociation rate $k_{\rm H_2}$ to the $\rm H^-$ photo-detachment rate $k_{\rm H^-}$. This assumption, however, could be too narrow in scope as the nature of the background radiation field is not known precisely. In this work we argue that the critical condition for suppressing the $\rm H_2$ cooling in the collapsing gas could be described in a more general way by a combination of $k_{\rm H_2}$ and $k_{\rm H^-}$ parameters. By performing a series of cosmological zoom-in simulations with an encompassing set of $k_{\rm H_2}$ and $k_{\rm H^-}$, we examine the gas flow by following evolution of basic parameters of the accretion flow. We test under what conditions the gas evolution is dominated by $\rm H_{2}$ and/or atomic cooling. We confirm the existence of a critical curve in the $k_{\rm H_2}-k_{\rm H^-}$ plane, and provide an analytical fit to it. This curve depends on the conditions in the direct collapse, and reveals domains where the atomic cooling dominates over the molecular cooling. Furthermore, we have considered the effect of $\rm H_{2}$ self-shielding on the critical curve, by adopting three methods for the effective column density approximation in $\rm H_{2}$. We find that the estimate of the characteristic length-scale for shielding can be improved by using $λ_{\rm Jeans25}$, which is 0.25 times that of the local Jeans length.

preprint2020arXiv

Statistical correlation between the distribution of Ly$α$ emitters and IGM HI at $z\sim2.2$ mapped by Subaru/Hyper Suprime-Cam

The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telescope. We select out 2,642 Ly$α$ emitter (LAE) candidates at $z=2.177\pm0.023$ down to the Ly$α$ luminosity of $L_{\text{Ly}α}\approx 2 \times 10^{42} {\rm erg~s}^{-1}$ to construct the galaxy overdensity maps, covering an effective area of 5.39 deg$^2$. Combining the sample with the Ly$α$ absorption estimated from 64 (e)BOSS quasar spectra, we find a moderate to strong correlation between the LAE overdensity $δ_{\rm LAE}$ and the effective optical depth $τ_{\rm LoS}$ in line-of-sights, with $P$-value$=0.09\%$ ($<0.01\%$) when the field that contains a significant quasar overdensity is in(ex)cluded. The cross-correlation analysis also clearly suggests that up to $4\pm1$ pMpc, LAEs tend to cluster in the regions rich in HI gas, indicated by the high $τ_{\rm LoS}$, and avoid the low $τ_{\rm LoS}$ region where the HI gas is deficient. By averaging the $τ_{\rm LoS}$ as a function of the projected distance ($d$) to LAEs, we find a $30\%$ excess signal at $2σ$ level at $d<200$ pkpc, indicating the dense HI in circumgalactic medium, and a tentative excess at $400<d<600$ pkpc in IGM regime, corroborating the cross-correlation signal detected at about $0.5$ pMpc. These statistical analyses indicate that galaxy$-$IGM HI correlations exist on scales ranging from several hundred pkpc to several pMpc at $z\sim2.2$.

preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

preprint2019arXiv

Galaxy evolution and radiative properties in the early Universe: multi-wavelength analysis in cosmological simulations

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

preprint2019arXiv

Galaxy simulation with the evolution of grain size distribution

We compute the evolution of interstellar dust in a hydrodynamic simulation of an isolated disc galaxy. We newly implement the evolution of full grain size distribution by sampling 32 grid points on the axis of the grain radius. We solve it consistently with the chemical enrichment and hydrodynamic evolution of the galaxy. This enables us to theoretically investigate spatially resolved evolution of grain size distribution in a galaxy. The grain size distribution evolves from a large-grain-dominated ($\gtrsim 0.1~μ$m) phase to a small-grain production phase, eventually converging to a power-law-like grain size distribution similar to the so-called MRN distribution. We find that the small-grain abundance is higher in the dense ISM in the early epoch ($t\lesssim 1$ Gyr) because of efficient dust growth by accretion, while coagulation makes the small-grain abundance less enhanced in the dense ISM later. This leads to steeper extinction curves in the dense ISM than in the diffuse ISM in the early phase, while they show the opposite trend later. The radial trend is also described by faster evolution in the inner part. We also confirm that the simulation reproduces the observed trend in the relation between dust-to-gas ratio and metallicity, and in the radial gradients of dust-to-gas ratio and dust-to-metal ratio. Since the above change in the grain size distribution occurs at $t\sim 1$ Gyr, the age and density dependence of grain size distribution has a significant impact on the extinction curves even at high redshift.