Showing 1–50 of 69 results for author: Shimizu, I

Combined Pre-Supernova Alert System with Kamland and Super-Kamiokande

Authors: KamLAND, Super-Kamiokande Collaborations, :, Seisho Abe, Minori Eizuka, Sawako Futagi, Azusa Gando, Yoshihito Gando, Shun Goto, Takahiko Hachiya, Kazumi Hata, Koichi Ichimura, Sei Ieki, Haruo Ikeda, Kunio Inoue, Koji Ishidoshiro, Yuto Kamei, Nanami Kawada, Yasuhiro Kishimoto, Masayuki Koga, Maho Kurasawa, Tadao Mitsui, Haruhiko Miyake, Daisuke Morita, Takeshi Nakahata , et al. (290 additional authors not shown)

Abstract: Preceding a core-collapse supernova, various processes produce an increasing amount of neutrinos of all flavors characterized by mounting energies from the interior of massive stars. Among them, the electron antineutrinos are potentially detectable by terrestrial neutrino experiments such as KamLAND and Super-Kamiokande via inverse beta decay interactions. Once these pre-supernova neutrinos are ob… ▽ More Preceding a core-collapse supernova, various processes produce an increasing amount of neutrinos of all flavors characterized by mounting energies from the interior of massive stars. Among them, the electron antineutrinos are potentially detectable by terrestrial neutrino experiments such as KamLAND and Super-Kamiokande via inverse beta decay interactions. Once these pre-supernova neutrinos are observed, an early warning of the upcoming core-collapse supernova can be provided. In light of this, KamLAND and Super-Kamiokande, both located in the Kamioka mine in Japan, have been monitoring pre-supernova neutrinos since 2015 and 2021, respectively. Recently, we performed a joint study between KamLAND and Super-Kamiokande on pre-supernova neutrino detection. A pre-supernova alert system combining the KamLAND detector and the Super-Kamiokande detector was developed and put into operation, which can provide a supernova alert to the astrophysics community. Fully leveraging the complementary properties of these two detectors, the combined alert is expected to resolve a pre-supernova neutrino signal from a 15 M$_{\odot}$ star within 510 pc of the Earth, at a significance level corresponding to a false alarm rate of no more than 1 per century. For a Betelgeuse-like model with optimistic parameters, it can provide early warnings up to 12 hours in advance. △ Less

Submitted 1 July, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: Resubmitted to ApJ. 22 pages, 16 figures, for more information about the combined pre-supernova alert system, see https://www.lowbg.org/presnalarm/

Constraining dark matter model using 21cm line intensity mapping

Authors: Koya Murakami, Kenji Kadota, Atsushi J. Nishizawa, Kentaro Nagamine, Ikkoh Shimizu

Abstract: We apply the Convolutional Neural Networks (CNNs) to the mock 21cm maps from the post-reionization epoch to show that the $Λ$CDM and warm dark matter (WDM) model can be distinguished for WDM particle masses $m_{FD}<3$ keV, under the assumption of thermal production of WDM following the Fermi-Dirac (FD) distribution. We demonstrate that the CNN is a potent tool in distinguishing the dark matter mas… ▽ More We apply the Convolutional Neural Networks (CNNs) to the mock 21cm maps from the post-reionization epoch to show that the $Λ$CDM and warm dark matter (WDM) model can be distinguished for WDM particle masses $m_{FD}<3$ keV, under the assumption of thermal production of WDM following the Fermi-Dirac (FD) distribution. We demonstrate that the CNN is a potent tool in distinguishing the dark matter masses, highlighting its sensitivity to the subtle differences in the 21cm maps produced by varying dark matter masses. Furthermore, we extend our analysis to encompass different WDM production mechanisms, recognizing that the dark matter production mechanism in the early Universe is among the sources of the most significant uncertainty for the dark matter model building. In this work, given the mass of the dark matter, we discuss the feasibility of discriminating four different WDM models: Fermi-Dirac (FD) distribution model, neutrino Minimal Standard Model ($ν$MSM), Dodelson-Widrow (DW), and Shi-Fuller (SF) model. For instance, when the WDM mass is 2 keV, we show that one can differentiate between CDM, FD, $ν$MSM, and DW models while discerning between the DW and SF models turns out to be challenging. Our results reinforce the viability of the CNN as a robust analysis for 21cm maps and shed light on its potential to unravel the features associated with different dark matter production mechanisms. △ Less

Submitted 10 March, 2024; originally announced March 2024.

Comments: 16 pages, 13 figures

The AGORA High-resolution Galaxy Simulations Comparison Project IV: Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at $z\le2$

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Joel R. Primack, Minyong Jung, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Clayton Strawn, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Thomas R. Quinn, Eun-jin Shin, Alvaro Segovia-Otero, Oscar Agertz, Kirk S. S. Barrow, Corentin Cadiou, Avishai Dekel, Cameron Hummels , et al. (3 additional authors not shown)

Abstract: In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at… ▽ More In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org <http://www.agorasimulations.org/>) for more information

The AGORA High-resolution Galaxy Simulations Comparison Project. V: Satellite Galaxy Populations In A Cosmological Zoom-in Simulation of A Milky Way-mass Halo

Authors: Minyong Jung, Santi Roca-Fàbrega, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Daniel Ceverino, Joel R. Primack, Thomas R. Quinn, Clayton Strawn, Tom Abel, Avishai Dekel, Bili Dong, Boon Kiat Oh, Romain Teyssier

Abstract: We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy… ▽ More We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it AGORA} simulations, and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called ``missing satellite problem' is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ proper pc at $z=2$). We also compare other properties such as the stellar mass$-$halo mass relation and the mass$-$metallicity relation. Our work highlights the value of comparison studies such as {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

Journal ref: ApJ 964 123 (2024)

The AGORA High-resolution Galaxy Simulations Comparison Project. VI. Similarities and Differences in the Circumgalactic Medium

Authors: Clayton Strawn, Santi Roca-Fàbrega, Joel R. Primack, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Minyong Jung, Thomas R. Quinn, Eun-jin Shin, Kirk S. S. Barrow, Avishai Dekel, Boon Kiat Oh, Nir Mandelker, Romain Teyssier, Cameron Hummels , et al. (4 additional authors not shown)

Abstract: We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the str… ▽ More We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 32 pages, 17 figures. Published in ApJ

Journal ref: ApJ 962 29 (2024)

Impact of astrophysical effects on the dark matter mass constraint with 21cm intensity mapping

Authors: Koya Murakami, Atsushi J. Nishizawa, Kentaro Nagamine, Ikko Shimizu

Abstract: We present an innovative approach to constraining the non-cold dark matter model using a convolutional neural network (CNN). We perform a suite of hydrodynamic simulations with varying dark matter particle masses and generate mock 21cm radio intensity maps to trace the dark matter distribution. Our proposed method complements the traditional power spectrum analysis. We compare our CNN classificati… ▽ More We present an innovative approach to constraining the non-cold dark matter model using a convolutional neural network (CNN). We perform a suite of hydrodynamic simulations with varying dark matter particle masses and generate mock 21cm radio intensity maps to trace the dark matter distribution. Our proposed method complements the traditional power spectrum analysis. We compare our CNN classification results with those from the power spectrum of the differential brightness temperature map of 21cm radiation, and find that the CNN outperforms the latter. Moreover, we investigate the impact of baryonic physics on the dark matter model constraint, including star formation, self-shielding of HI gas, and UV background model. We find that these effects may introduce some contamination in the dark matter constraint, but they are insignificant when compared to the realistic system noise of the SKA instruments. △ Less

Submitted 21 April, 2024; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: 16 pages, 15 figures

The 300 pc resolution imaging of a z = 8.31 galaxy: Turbulent ionized gas and potential stellar feedback 600 million years after the Big Bang

Authors: Yoichi Tamura, Tom J. L. C. Bakx, Akio K. Inoue, Takuya Hashimoto, Tsuyoshi Tokuoka, Chihiro Imamura, Bunyo Hatsukade, Minju M. Lee, Kana Moriwaki, Takashi Okamoto, Kazuaki Ota, Hideki Umehata, Naoki Yoshida, Erik Zackrisson, Masato Hagimoto, Hiroshi Matsuo, Ikkoh Shimizu, Yuma Sugahara, Tsutomu T. Takeuchi

Abstract: We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $μ$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity… ▽ More We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $μ$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity gradient unlike what was found in [CII] 158 $μ$m at a larger scale, suggesting random bulk motion of ionized gas clouds inside the galaxy. In contrast, dust emission appears as two individual clumps apparently separating or bridging the [OIII]/stellar clumps. The cross correlation coefficient between dust and ultraviolet-related emission (i.e., [OIII] and ultraviolet continuum) is unity on a galactic scale, while it drops at < 1 kpc, suggesting well mixed geometry of multi-phase interstellar media on sub-kpc scales. If the cutoff scale characterizes different stages of star formation, the cutoff scale can be explained by gravitational instability of turbulent gas. We also report on a kpc-scale off-center cavity embedded in the dust continuum image. This could be a superbubble producing galactic-scale outflows, since the energy injection from the 4 Myr starburst suggested by a spectral energy distribution analysis is large enough to push the surrounding media creating a kpc-scale cavity. △ Less

Submitted 16 July, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 18 pages, 6 figures, 1 table, published in ApJ

Updated measurements of [O III] 88 $μ$m, [C II] 158 $μ$m, and Dust Continuum Emission from a z=7.2 Galaxy

Authors: Yi W. Ren, Yoshinobu Fudamoto, Akio K. Inoue, Yuma Sugahara, Tsuyoshi Tokuoka, Yoichi Tamura, Hiroshi Matsuo, Kotaro Kohno, Hideki Umehata, Takuya Hashimoto, Rychard J. Bouwens, Renske Smit, Nobunari Kashikawa, Takashi Okamoto, Takatoshi Shibuya, Ikkoh Shimizu

Abstract: We present updated measurements of the [O III] 88 $μ$m, [C II] 158 $μ$m, and dust continuum emission from a star-forming galaxy at $z=7.212$, SXDF-NB1006-2, by utilizing Atacama Large Millimeter/submillimeter Array (ALMA) archival data sets analysed in previous studies and data sets that have not been analysed before. The follow-up ALMA observations with higher angular resolution and sensitivity r… ▽ More We present updated measurements of the [O III] 88 $μ$m, [C II] 158 $μ$m, and dust continuum emission from a star-forming galaxy at $z=7.212$, SXDF-NB1006-2, by utilizing Atacama Large Millimeter/submillimeter Array (ALMA) archival data sets analysed in previous studies and data sets that have not been analysed before. The follow-up ALMA observations with higher angular resolution and sensitivity reveal a clumpy structure of the [O III] emission on a scale of $0.32-0.85\,\rm{kpc}$. We also combined all the ALMA [O III] ([C II]) data sets and updated the [O III] ([C II]) detection to $5.9σ$ ($3.6σ-4.5σ$). The non-detection of [C II] with data from the REBELS large program implies the incompleteness of spectral-scan surveys using [C II] to detect galaxies with high star formation rates (SFRs) but marginal [C II] emission at high-$z$. The dust continuum at 90 $μ$m and 160 $μ$m remains undetected, indicating little dust content of $<3.9\times10^{6}\,M_\odot\,(3σ)$, and we obtained a more stringent constraint on the total infrared luminosity. We updated the [O III]/[C II] luminosity ratios to $10.2\pm4.7~(6.1\pm3.5$) and $20\pm12~(9.6\pm6.1$) for $4.5σ$ and $3.6σ$ [C II] detections, respectively, where the ratios in the parentheses are corrected for the surface brightness dimming effect on the extended [C II] emission. We also found a strong [C II] deficit ($0.6-1.3$ dex) between SXDF-NB1006-2 and the mean $L_{\rm{[CII]}}-\rm{SFR}$ relation of galaxies at $0<z<9$. △ Less

Submitted 9 March, 2023; v1 submitted 5 February, 2023; originally announced February 2023.

Comments: 13 pages, 10 figures, 8 tables, published in ApJ

Star Formation and Chemical Enrichment in Protoclusters

Authors: Keita Fukushima, Kentaro Nagamine, Ikkoh Shimizu

Abstract: We examine star formation and chemical enrichment in protoclusters (PCs) using cosmological zoom-in hydrodynamic simulations. We find that the total star formation rate (SFR) in all PC ($>10^{14.4}\,h^{-1}$M$_\odot$) reaches $>10^4\,\mathrm{M}_\odot \mathrm{yr}^{-1}$ at $z=3$, equivalent to the observed PCs. The SFR in the Core region accounts for about $30\%$ of the total star formation in the PC… ▽ More We examine star formation and chemical enrichment in protoclusters (PCs) using cosmological zoom-in hydrodynamic simulations. We find that the total star formation rate (SFR) in all PC ($>10^{14.4}\,h^{-1}$M$_\odot$) reaches $>10^4\,\mathrm{M}_\odot \mathrm{yr}^{-1}$ at $z=3$, equivalent to the observed PCs. The SFR in the Core region accounts for about $30\%$ of the total star formation in the PC at $z\gtrsim1$, suggesting the importance of the outer regions to reveal the evolution of galaxy clusters. We find that the total SFR of PC is dominated by galaxies with $10^{10}\,\le\,(\mathrm{M}_\star/M_\odot)\,\le\,10^{11}$, while more massive galaxies dominate the SFR in the Core. For the chemical abundance evolution, we find that the higher-density region has a higher metallicity and faster evolution. We show that the [O/Fe] vs. [Fe/H] relation turns down in the Core at $z=3.4$ due to the enrichment of Fe by Type Ia supernovae. We find no environmental effects for the mass--metallicity relations (MZR) or $\log$(N/O) vs. $12+\log$(O/H) for galaxies. We find that the chemical enrichment in galaxy clusters proceeds faster in the high redshift Universe ($z>1$). Our work will benefit future tomographic observations, particularly using PCs as unique probes of accelerated structure formation and evolution in high-density regions of the universe. △ Less

Submitted 13 September, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: 19 pages, 18 figures, 2 tables, 2 appendices, Accepted for publication in MNRAS

JWST reveals a possible $z \sim 11$ galaxy merger in triply-lensed MACS0647$-$JD

Authors: Tiger Yu-Yang Hsiao, Dan Coe, Abdurro'uf, Lily Whitler, Intae Jung, Gourav Khullar, Ashish Kumar Meena, Pratika Dayal, Kirk S. S. Barrow, Lillian Santos-Olmsted, Adam Casselman, Eros Vanzella, Mario Nonino, Yolanda Jimenez-Teja, Masamune Oguri, Daniel P. Stark, Lukas J. Furtak, Adi Zitrin, Angela Adamo, Gabriel Brammer, Larry Bradley, Jose M. Diego, Erik Zackrisson, Steven L. Finkelstein, Rogier A. Windhorst , et al. (41 additional authors not shown)

Abstract: MACS0647$-$JD is a triply-lensed $z\sim11$ galaxy originally discovered with the Hubble Space Telescope. Here we report new JWST imaging, which clearly resolves MACS0647$-$JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. Both are very small, with stellar masses $\sim10^8\,M_\odot$ and radii $r<100\,\rm pc$. The brighter larger component "A"… ▽ More MACS0647$-$JD is a triply-lensed $z\sim11$ galaxy originally discovered with the Hubble Space Telescope. Here we report new JWST imaging, which clearly resolves MACS0647$-$JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. Both are very small, with stellar masses $\sim10^8\,M_\odot$ and radii $r<100\,\rm pc$. The brighter larger component "A" is intrinsically very blue ($β\sim-2.6$), likely due to very recent star formation and no dust, and is spatially extended with an effective radius $\sim70\,\rm pc$. The smaller component "B" appears redder ($β\sim-2$), likely because it is older ($100-200\,\rm Myr$) with mild dust extinction ($A_V\sim0.1\,\rm mag$), and a smaller radius $\sim20\,\rm pc$. We identify galaxies with similar colors in a high-redshift simulation, finding their star formation histories to be out of phase. With an estimated stellar mass ratio of roughly 2:1 and physical projected separation $\sim400\,\rm pc$, we may be witnessing a galaxy merger 400 million years after the Big Bang. We also identify a candidate companion galaxy C $\sim3\,{\rm kpc}$ away, likely destined to merge with galaxies A and B. The combined light from galaxies A+B is magnified by factors of $\sim$8, 5, and 2 in three lensed images JD1, 2, and 3 with F356W fluxes $\sim322$, $203$, $86\,\rm nJy$ (AB mag 25.1, 25.6, 26.6). MACS0647$-$JD is significantly brighter than other galaxies recently discovered at similar redshifts with JWST. Without magnification, it would have AB mag 27.3 ($M_{UV}=-20.4$). With a high confidence level, we obtain a photometric redshift of $z=10.6\pm0.3$ based on photometry measured in 6 NIRCam filters spanning $1-5\rmμm$, out to $4300\,Å$ rest-frame. JWST NIRSpec observations planned for January 2023 will deliver a spectroscopic redshift and a more detailed study of the physical properties of MACS0647$-$JD. △ Less

Submitted 31 May, 2023; v1 submitted 25 October, 2022; originally announced October 2022.

Comments: 27 pages, 14 figures, accepted for publication in ApJL

Possible Systematic Rotation in the Mature Stellar Population of a $z=9.1$ Galaxy

Authors: Tsuyoshi Tokuoka, Akio K. Inoue, Takuya Hashimoto, Richard S. Ellis, Nicolas Laporte, Yuma Sugahara, Hiroshi Matsuo, Yoichi Tamura, Yoshinobu Fudamoto, Kana Moriwaki, Guido Roberts-Borsani, Ikkoh Shimizu, Satoshi Yamanaka, Naoki Yoshida, Erik Zackrisson, Wei Zheng

Abstract: We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$μ$m emission is detected at 10$σ$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a… ▽ More We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$μ$m emission is detected at 10$σ$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $ΔV_{\rm obs}/2σ_{\rm tot}=0.84\pm0.23$, where $ΔV_{\rm obs}$ is the observed maximum velocity difference and $σ_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/σ_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $σ_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$. △ Less

Submitted 28 May, 2022; originally announced May 2022.

Comments: ApJL accepted

Search for supernova neutrinos and constraint on the galactic star formation rate with the KamLAND data

Authors: S. Abe, S. Asami, M. Eizuka, S. Futagi, A. Gando, Y. Gando, T. Gima, A. Goto, T. Hachiya, K. Hata, K. Hosokawa, K. Ichimura, S. Ieki, H. Ikeda, K. Inoue, K. Ishidoshiro, Y. Kamei, N. Kawada, Y. Kishimoto, M. Koga, M. Kurasawa, N. Maemura, T. Mitsui, H. Miyake, T. Nakahata , et al. (42 additional authors not shown)

Abstract: We present the results of a search for core-collapse supernova neutrinos, using long-term KamLAND data from 2002 March 9 to 2020 April 25. We focus on the electron antineutrinos emitted from supernovae in the energy range of 1.8--111 MeV. Supernovae will make a neutrino event cluster with the duration of $\sim$10 s in the KamLAND data. We find no neutrino clusters and give the upper limit on the s… ▽ More We present the results of a search for core-collapse supernova neutrinos, using long-term KamLAND data from 2002 March 9 to 2020 April 25. We focus on the electron antineutrinos emitted from supernovae in the energy range of 1.8--111 MeV. Supernovae will make a neutrino event cluster with the duration of $\sim$10 s in the KamLAND data. We find no neutrino clusters and give the upper limit on the supernova rate as to be 0.15 yr$^{-1}$ with a 90% confidence level. The detectable range, which corresponds to a >95% detection probability, is 40--59 kpc and 65--81 kpc for core-collapse supernovae and failed core-collapse supernovae, respectively. This paper proposes to convert the supernova rate obtained by the neutrino observation to the Galactic star formation rate. Assuming a modified Salpeter-type initial mass function, the upper limit on the Galactic star formation rate is <(17.5--22.7) $M_{\odot} \mathrm{yr}^{-1}$ with a 90% confidence level. △ Less

Submitted 29 July, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 10 pages, 4 figures, accepted for publication in Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 934, Number 1, Page 85 (2022)

KamLAND's search for correlated low-energy electron antineutrinos with astrophysical neutrinos from IceCube

Authors: S. Abe, S. Asami, M. Eizuka, S. Futagi, A. Gando, Y. Gando, T. Gima, A. Goto, T. Hachiya, K. Hata, K. Hosokawa, K. Ichimura, S. Ieki, H. Ikeda, K. Inoue, K. Ishidoshiro, Y. Kamei, N. Kawada, Y. Kishimoto, T. Kinoshita, M. Koga, M. Kurasawa, N. Maemura, T. Mitsui, H. Miyake , et al. (45 additional authors not shown)

Abstract: We report the results of a search for MeV-scale astrophysical neutrinos in KamLAND presented as an excess in the number of coincident neutrino interactions associated with the publicly available high-energy neutrino datasets from the IceCube Neutrino Observatory. We find no statistically significant excess in the number of observed low-energy electron antineutrinos in KamLAND, given a coincidence… ▽ More We report the results of a search for MeV-scale astrophysical neutrinos in KamLAND presented as an excess in the number of coincident neutrino interactions associated with the publicly available high-energy neutrino datasets from the IceCube Neutrino Observatory. We find no statistically significant excess in the number of observed low-energy electron antineutrinos in KamLAND, given a coincidence time window of $\pm$500s, $\pm$1,000s, $\pm$3,600s, and $\pm$10,000s around each of the IceCube neutrinos. We use this observation to present limits from 1.8 MeV to 100 MeV on the electron antineutrino fluence, assuming a mono-energetic flux. We then compare the results to several astrophysical measurements performed by IceCube and place a limit at the 90% confidence level on the electron antineutrino isotropic thermal luminosity from the TXS 0506+056 blazar. △ Less

Submitted 21 July, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 12 pages, 5 figures

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

Authors: Yuri Oku, Kengo Tomida, Kentaro Nagamine, Ikkoh Shimizu, Renyue Cen

Abstract: 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 scale… ▽ More 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. △ Less

Submitted 26 July, 2022; v1 submitted 3 January, 2022; originally announced January 2022.

Comments: 36 pages, 22 figures, accepted for publication in The Astrophysical Journal Supplement Series. Movies available at https://www.yurioku.com/research/

A search for correlated low-energy electron antineutrinos in KamLAND with gamma-ray bursts

Authors: S. Abe, S. Asami, A. Gando, Y. Gando, T. Gima, A. Goto, T. Hachiya, K. Hata, K. Hosokawa, K. Ichimura, S. Ieki, H. Ikeda, K. Inoue, K. Ishidoshiro, Y. Kamei, N. Kawada, Y. Kishimoto, T. Kinoshita, M. Koga, N. Maemura, T. Mitsui, H. Miyake, K. Nakamura, K. Nakamura, R. Nakamura , et al. (43 additional authors not shown)

Abstract: We present the results of a time-coincident event search for low-energy electron antineutrinos in the KamLAND detector with gamma-ray bursts from the Gamma-ray Coordinates Network and Fermi Gamma-ray Burst Monitor. Using a variable coincidence time window of $\pm$500s plus the duration of each gamma-ray burst, no statistically significant excess above background is observed. We place the world's m… ▽ More We present the results of a time-coincident event search for low-energy electron antineutrinos in the KamLAND detector with gamma-ray bursts from the Gamma-ray Coordinates Network and Fermi Gamma-ray Burst Monitor. Using a variable coincidence time window of $\pm$500s plus the duration of each gamma-ray burst, no statistically significant excess above background is observed. We place the world's most stringent 90% confidence level upper limit on the electron antineutrino fluence below 17.5 MeV. Assuming a Fermi-Dirac neutrino energy spectrum from the gamma-ray burst source, we use the available redshift data to constrain the electron antineutrino luminosity and effective temperature. △ Less

Submitted 24 January, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 9 pages, 3 figures

Limits on astrophysical antineutrinos with the KamLAND experiment

Authors: S. Abe, S. Asami, A. Gando, Y. Gando, T. Gima, A. Goto, T. Hachiya, K. Hata, S. Hayashida, K. Hosokawa, K. Ichimura, S. Ieki, H. Ikeda, K. Inoue, K. Ishidoshiro, Y. Kamei, N. Kawada, T. Kinoshita, Y. Kishimoto, M. Koga, N. Maemura, T. Mitsui, H. Miyake, K. Nakamura, K. Nakamura , et al. (45 additional authors not shown)

Abstract: We report on a search for electron antineutrinos ($\barν_e$) from astrophysical sources in the neutrino energy range 8.3 to 30.8 MeV with the KamLAND detector. In an exposure of 6.72 kton-year of the liquid scintillator, we observe 18 candidate events via the inverse beta decay reaction. Although there is a large background uncertainty from neutral current atmospheric neutrino interactions, we fin… ▽ More We report on a search for electron antineutrinos ($\barν_e$) from astrophysical sources in the neutrino energy range 8.3 to 30.8 MeV with the KamLAND detector. In an exposure of 6.72 kton-year of the liquid scintillator, we observe 18 candidate events via the inverse beta decay reaction. Although there is a large background uncertainty from neutral current atmospheric neutrino interactions, we find no significant excess over background model predictions. Assuming several supernova relic neutrino spectra, we give upper flux limits of 60--110 cm$^{-2}$ s$^{-1}$ (90% CL) in the analysis range and present a model-independent flux. We also set limits on the annihilation rates for light dark matter pairs to neutrino pairs. These data improves on the upper probability limit of $^{8}$B solar neutrinos converting into $\barν_e$'s, $P_{ν_e \rightarrow \barν_e} < 3.5\times10^{-5}$ (90% CL) assuming an undistorted $\barν_e$ shape. This corresponds to a solar $\barν_e$ flux of 60 cm$^{-2}$ s$^{-1}$ (90% CL) in the analysis energy range. △ Less

Submitted 22 October, 2021; v1 submitted 19 August, 2021; originally announced August 2021.

Comments: 21 pages, 9 figures, 4 tables, accepted for publication in Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 925, Number 1, Page 14 (2022)

Mapping Lyman-alpha forest three-dimensional large scale structure in real and redshift space

Authors: Francesco Sinigaglia, Francisco-Shu Kitaura, Andrés Balaguera-Antolínez, Ikkoh Shimizu, Kentaro Nagamine, Manuel Sánchez-Benavente, Metin Ata

Abstract: This work presents a new physically-motivated supervised machine learning method, Hydro-BAM, to reproduce the three-dimensional Lyman-$α$ forest field in real and in redshift space learning from a reference hydrodynamic simulation, thereby saving about 7 orders of magnitude in computing time. We show that our method is accurate up to $k\sim1\,h\,\rm{Mpc}^{-1}$ in the one- (PDF), two- (power-spectr… ▽ More This work presents a new physically-motivated supervised machine learning method, Hydro-BAM, to reproduce the three-dimensional Lyman-$α$ forest field in real and in redshift space learning from a reference hydrodynamic simulation, thereby saving about 7 orders of magnitude in computing time. We show that our method is accurate up to $k\sim1\,h\,\rm{Mpc}^{-1}$ in the one- (PDF), two- (power-spectra) and three-point (bi-spectra) statistics of the reconstructed fields. When compared to the reference simulation including redshift space distortions, our method achieves deviations of $\lesssim2\%$ up to $k=0.6\,h\,\rm{Mpc}^{-1}$ in the monopole, $\lesssim5\%$ up to $k=0.9\,h\,\rm{Mpc}^{-1}$ in the quadrupole. The bi-spectrum is well reproduced for triangle configurations with sides up to $k=0.8\,h\,\rm{Mpc}^{-1}$. In contrast, the commonly-adopted Fluctuating Gunn-Peterson approximation shows significant deviations already neglecting peculiar motions at configurations with sides of $k=0.2-0.4\,h\,\rm{Mpc}^{-1}$ in the bi-spectrum, being also significantly less accurate in the power-spectrum (within 5$\%$ up to $k=0.7\,h\,\rm{Mpc}^{-1}$). We conclude that an accurate analysis of the Lyman-$α$ forest requires considering the complex baryonic thermodynamical large-scale structure relations. Our hierarchical domain specific machine learning method can efficiently exploit this and is ready to generate accurate Lyman-$α$ forest mock catalogues covering large volumes required by surveys such as DESI and WEAVE. △ Less

Submitted 3 February, 2022; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: Accepted for publication by ApJ

The AGORA High-resolution Galaxy Simulations Comparison Project. III: Cosmological zoom-in simulation of a Milky Way-mass halo

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Loic Hausammann, Kentaro Nagamine, Johnny W. Powell, Ikkoh Shimizu, Daniel Ceverino, Alessandro Lupi, Joel R. Primack, Thomas Quinn, Yves Revaz, Héctor Velázquez, Tom Abel, Michael Buehlmann, Avishai Dekel, Bili Dong, Oliver Hahn, Cameron B. Hummels, Ki-won Kim, Britton D. Smith, Clayton J. Strawn, Romain Teyssier, Matthew Turk

Abstract: We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now ac… ▽ More We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now account for the effects of cosmological processes. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. After all the participating code groups successfully completed the calibration steps, we reach a suite of cosmological simulations with similar mass assembly histories down to $z=4$. With numerical accuracy that resolves the internal structure of a target halo, we find that the codes overall agree well with one another in e.g., gas and stellar properties, but also show differences in e.g., circumgalactic medium properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, the results of high-resolution cosmological zoom-in simulations can be robust and reproducible. New code groups are invited to join this comparison by generating equivalent models by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the simulations presented here will be in forthcoming reports from our Collaboration. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: Accepted for publication to the ApJ. The simulation snapshots used in this paper will be made publicly available in a user-friendly interface. The announcement of the data release will be posted on the AGORA Project webpage (www.agorasimulations.org). Early access to the data can be provided upon request to the authors, or to the project coordinator (sroca01@ucm.es, santacruzgalaxy@gmail.com)

Search for Solar Flare Neutrinos with the KamLAND detector

Authors: S. Abe, S. Asami, A. Gando, Y. Gando, T. Gima, A. Goto, T. Hachiya, K. Hata, S. Hayashida, K. Hosokawa, K. Ichimura, S. Ieki, H. Ikeda, K. Inoue, K. Ishidoshiro, Y. Kamei, N. Kawada, Y. Kishimoto, T. Kinoshita, M. Koga, N. Maemura, T. Mitsui, H. Miyake, K. Nakamura, K. Nakamura , et al. (44 additional authors not shown)

Abstract: We report the result of a search for neutrinos in coincidence with solar flares from the GOES flare database. The search was performed on a 10.8 kton-year exposure of KamLAND collected from 2002 to 2019. This large exposure allows us to explore previously unconstrained parameter space for solar flare neutrinos. We found no statistical excess of neutrinos and established 90% confidence level upper… ▽ More We report the result of a search for neutrinos in coincidence with solar flares from the GOES flare database. The search was performed on a 10.8 kton-year exposure of KamLAND collected from 2002 to 2019. This large exposure allows us to explore previously unconstrained parameter space for solar flare neutrinos. We found no statistical excess of neutrinos and established 90% confidence level upper limits of $8.4 \times 10^7$ cm$^{-2}$ ($3.0 \times 10^{9}$ cm$^{-2}$) on electron anti-neutrino (electron neutrino) fluence at 20 MeV normalized to the X12 flare, assuming that the neutrino fluence is proportional to the X-ray intensity. △ Less

Submitted 26 October, 2021; v1 submitted 6 May, 2021; originally announced May 2021.

Comments: 13 pages, 9 figures, accepted October 27, 2021

Journal ref: The Astrophysical Journal, Volume 924, Number 2, Page 103 (2022)

Catch Me if You Can: Biased Distribution of Ly$α$-emitting Galaxies according to the Viewing Direction

Authors: Rieko Momose, Kazuhiro Shimasaku, Kentaro Nagamine, Ikkoh Shimizu, Nobunari Kashikawa, Makoto Ando, Haruka Kusakabe

Abstract: We report that Ly$α$-emitting galaxies (LAEs) may not faithfully trace the cosmic web of neutral hydrogen (HI), but their distribution is likely biased depending on the viewing direction. We calculate the cross-correlation (CCF) between galaxies and Ly$α$ forest transmission fluctuations on the near and far sides of the galaxies separately, for three galaxy samples at $z\sim2$: LAEs, [OIII] emitte… ▽ More We report that Ly$α$-emitting galaxies (LAEs) may not faithfully trace the cosmic web of neutral hydrogen (HI), but their distribution is likely biased depending on the viewing direction. We calculate the cross-correlation (CCF) between galaxies and Ly$α$ forest transmission fluctuations on the near and far sides of the galaxies separately, for three galaxy samples at $z\sim2$: LAEs, [OIII] emitters (O3Es), and continuum-selected galaxies. We find that only LAEs have anisotropic CCFs, with the near side one showing lower signals up to $r=3-4~h^{-1}$ comoving Mpc. This means that the average HI density on the near side of LAEs is lower than that on the far-side by a factor of $2.1$ under the Fluctuating Gunn-Peterson Approximation. Mock LAEs created by assigning Ly$α$ equivalent width ($EW_\text{Ly$α$}^\text{obs}$) values to O3Es with an empirical relation also show similar, anisotropic CCFs if we use only objects with higher $EW_\text{Ly$α$}^\text{obs}$ than a certain threshold. These results indicate that galaxies on the far side of a dense region are more difficult to be detected ("hidden") in Ly$α$ because Ly$α$ emission toward us is absorbed by dense neutral hydrogen. If the same region is viewed from a different direction, a different set of LAEs will be selected as if galaxies are playing hide-and-seek using HI gas. Care is needed when using LAEs to search for overdensities. △ Less

Submitted 21 April, 2021; originally announced April 2021.

Comments: 10 pages, 5 figures, Accepted for publication on ApJL

Supernova Model Discrimination with Hyper-Kamiokande

Authors: Hyper-Kamiokande Collaboration, :, K. Abe, P. Adrich, H. Aihara, R. Akutsu, I. Alekseev, A. Ali, F. Ameli, I. Anghel, L. H. V. Anthony, M. Antonova, A. Araya, Y. Asaoka, Y. Ashida, V. Aushev, F. Ballester, I. Bandac, M. Barbi, G. J. Barker, G. Barr, M. Batkiewicz-Kwasniak, M. Bellato, V. Berardi, M. Bergevin , et al. (478 additional authors not shown)

Abstract: Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-colla… ▽ More Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-collapse supernovae is not yet well understood. Hyper-Kamiokande is a next-generation neutrino detector that will be able to observe the neutrino flux from the next galactic core-collapse supernova in unprecedented detail. We focus on the first 500 ms of the neutrino burst, corresponding to the accretion phase, and use a newly-developed, high-precision supernova event generator to simulate Hyper-Kamiokande's response to five different supernova models. We show that Hyper-Kamiokande will be able to distinguish between these models with high accuracy for a supernova at a distance of up to 100 kpc. Once the next galactic supernova happens, this ability will be a powerful tool for guiding simulations towards a precise reproduction of the explosion mechanism observed in nature. △ Less

Submitted 20 July, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 21 pages, 7 figures. Article based on thesis published as arXiv:2002.01649. v2: added references and some explanations in response to reviewer comments

Journal ref: Astrophys.J. 916 (2021) 15

Search for Low-energy Electron Antineutrinos in KamLAND Associated with Gravitational Wave Events

Authors: S. Abe, S. Asami, A. Gando, Y. Gando, T. Gima, A. Goto, T. Hachiya, K. Hata, S. Hayashida, K. Hosokawa, K. Ichimura, S. Ieki, H. Ikeda, K. Inoue, K. Ishidoshiro, Y. Kamei, N. Kawada, Y. Kishimoto, T. Kinoshita, M. Koga, N. Maemura, T. Mitsui, H. Miyake, K. Nakamura, K. Nakamura , et al. (44 additional authors not shown)

Abstract: We present the results of a search for MeV-scale electron antineutrino events in KamLAND in coincident with the 60 gravitational wave events/candidates reported by the LIGO/Virgo collaboration during their second and third observing runs. We find no significant coincident signals within a $\pm$ 500 s timing window from each gravitational wave and present 90% C.L. upper limits on the electron antin… ▽ More We present the results of a search for MeV-scale electron antineutrino events in KamLAND in coincident with the 60 gravitational wave events/candidates reported by the LIGO/Virgo collaboration during their second and third observing runs. We find no significant coincident signals within a $\pm$ 500 s timing window from each gravitational wave and present 90% C.L. upper limits on the electron antineutrino fluence between $10^{8}$-$10^{13}\,{\mathrm cm^2}$ for neutrino energies in the energy range of 1.8-111 MeV. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: 11 pages, 2 figures, 2 tables, Accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 909, Number 2 (2021)

The bias from hydrodynamic simulations: mapping baryon physics onto dark matter fields

Authors: Francesco Sinigaglia, Francisco-Shu Kitaura, Andrés Balaguera-Antolínez, Kentaro Nagamine, Metin Ata, Ikkoh Shimizu, Manuel Sánchez-Benavente

Abstract: This paper investigates the hierarchy of baryon physics assembly bias relations obtained from state-of-the-art hydrodynamic simulations with respect to the underlying cosmic web spanned by the dark matter field. Using the Bias Assignment Method (BAM) we find that non-local bias plays a central role. We classify the cosmic web based on the invariants of the curvature tensor defined not only by the… ▽ More This paper investigates the hierarchy of baryon physics assembly bias relations obtained from state-of-the-art hydrodynamic simulations with respect to the underlying cosmic web spanned by the dark matter field. Using the Bias Assignment Method (BAM) we find that non-local bias plays a central role. We classify the cosmic web based on the invariants of the curvature tensor defined not only by the gravitational potential, but especially by the over-density, as small scale clustering becomes important in this context. First, the gas density bias relation can be directly mapped onto the dark matter density field to high precision exploiting the strong correlation between them. In a second step, the neutral hydrogen is mapped based on the dark matter and the gas density fields. Finally, the temperature is mapped based on the previous quantities. This permits us to statistically reconstruct the baryon properties within the same simulated volume finding percent-precision in the two-point statistics and compatible results in the three-point statistics, in general within 1-$σ$, with respect to the reference simulation (with 5 to 6 orders of magnitude less computing time). This paves the path to establish the best set-up for the construction of mocks probing the intergalactic medium for the generation of such key ingredients in the statistical analysis of large forthcoming missions such as DESI, Euclid, J-PAS and WEAVE. △ Less

Submitted 19 July, 2021; v1 submitted 12 December, 2020; originally announced December 2020.

Comments: Accepted for publication by ApJ

A puzzling non-detection of [O III] and [C II] from a $z \approx 7.7$ galaxy observed with ALMA

Authors: C. Binggeli, A. K. Inoue, T. Hashimoto, M. C. Toribio, E. Zackrisson, S. Ramstedt, K. Mawatari, Y. Harikane, H. Matsuo, T. Okamoto, K. Ota, I. Shimizu, Y. Tamura, Y. Taniguchi, H. Umehata

Abstract: Characterizing the galaxy population in the early Universe holds the key to understanding the evolution of these objects and the role they played in cosmic reionization. However, the number of observations at the very highest redshifts are to date, few. In order to shed light on the properties of galaxies in the high-redshift Universe and their interstellar media, we observe the Lyman-$α$ emitting… ▽ More Characterizing the galaxy population in the early Universe holds the key to understanding the evolution of these objects and the role they played in cosmic reionization. However, the number of observations at the very highest redshifts are to date, few. In order to shed light on the properties of galaxies in the high-redshift Universe and their interstellar media, we observe the Lyman-$α$ emitting galaxy z7\_GSD\_3811 at $z=7.664$ with band 6 and 8 at the Atacama Large Millimeter/submillimeter Array (ALMA). We target the far-infrared [O III] 88 $μm$, [C II] 158 $μm$ emission lines and dust continuum in the star-forming galaxy z7\_GSD\_3811 with ALMA. We combine these measurements with earlier observations in the rest-frame ultraviolet (UV) in order to characterize the object, and compare results to those of earlier studies observing [O III] and [C II] emission in high-redshift galaxies. The [O III] 88 $μm$ and [C II] 158 $μm$ emission lines are undetected at the position of z7\_GSD\_3811, with $3σ$ upper limits of $1.6 \ \times \ 10^{8} \ L_{\odot}$ and $4.0 \ \times \ 10^{7} \ L_{\odot}$, respectively. We do not detect any dust continuum in band 6 nor band 8. The measured rms in the band 8 and band 6 continuum is 26 and 9.9 $μJy \ beam^{-1}$, respectively. Similar to several other high-redshift galaxies, z7\_GSD\_3811 exhibits low [C II] emission for its star formation rate compared to local galaxies. Furthermore, our upper limit on the [O III] line luminosity is lower than all the previously observed [O III] lines in high-redshift galaxies with similar ultraviolet luminosities. Our ALMA band 6 and 8 dust continuum observations imply that z7\_GSD\_3811 likely has a low dust content, and our non-detections of the [O III] and [C II] lines could indicate that z7\_GSD\_3811 has a low metallicity ($Z \lesssim 0.1 \ Z_{\odot}$). △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: 13 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 646, A26 (2021)

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

Authors: Yongming Liang, Nobunari Kashikawa, Zheng Cai, Xiaohui Fan, J. Xavier Prochaska, Kazuhiro Shimasaku, Masayuki Tanaka, Hisakazu Uchiyama, Kei Ito, Rhythm Shimakawa, Kentaro Nagamine, Ikkoh Shimizu, Masafusa Onoue, Jun Toshikawa

Abstract: 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 Telesco… ▽ More 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$. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Comments: 27 pages, 15 figures; Resubmitted to ApJ after the first referee's report. Comments are welcome

Probing Feedback via IGM tomography and Ly$α$ forest with Subaru PFS, TMT/ELT, and JWST

Authors: Kentaro Nagamine, Ikkoh Shimizu, Katsumi Fujita, Nao Suzuki, Khee-Gan Lee, Rieko Momose, Shiro Mukae, Yongming Liang, Nobunari Kashikawa, Masami Ouchi, John Silverman

Abstract: In preparation for the IGM tomography study by Subaru Prime Focus Spectrograph (PFS) survey and other large future telescopes such as TMT/ELT/GMT, we present the results of our pilot study on Ly$α$ forest and IGM tomography statistics using the GADGET3-OSAKA cosmological smoothed particle hydrodynamical simulation. Our simulation includes models for star formation and supernova feedback, which ena… ▽ More In preparation for the IGM tomography study by Subaru Prime Focus Spectrograph (PFS) survey and other large future telescopes such as TMT/ELT/GMT, we present the results of our pilot study on Ly$α$ forest and IGM tomography statistics using the GADGET3-OSAKA cosmological smoothed particle hydrodynamical simulation. Our simulation includes models for star formation and supernova feedback, which enables more realistic cross-correlation studies between galaxies, neutral hydrogen (HI) and metals in circumgalactic and intergalactic medium. We create a light-cone data set at $z=2-3$ from our simulations and generate mock Ly$α$ forest data. As a first step, in this paper, we focus on the distribution of HI and galaxies, and present statistical results on 1D flux PDF, 1D power spectrum, flux contrast vs. impact parameter, HI$-$galaxy cross-correlations. Our results show overall agreement with current observational data, with some interesting discrepancies on small scales that are due to either feedback effects or varying observational conditions. Our simulation shows stronger HI absorption with decreasing transverse distance from galaxies. We find that the massive galaxies with $M_\star \ge 10^{10}\,M_\odot$ contribute strongly to the flux contrast signal, and that the lower-mass galaxies with $M_\star \sim 10^8-10^{10}\,M_\odot$ tend to dilute the flux contrast signal from massive galaxies. On large scales, the average flux contrast smoothly connects to the IGM level, supporting the concordance $Λ$ cold dark matter model. We also find an increase in the HI absorption toward the center of a protocluster. △ Less

Submitted 7 July, 2021; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 26 pages, 16 figures. Published in ApJ (matched to the published version)

Journal ref: ApJ, 914, 66 (2021)

Environmental Dependence of Galactic Properties Traced by Ly$α$ Forest Absorption: Diversity among Galaxy Populations

Authors: Rieko Momose, Kazuhiro Shimasaku, Nobunari Kashikawa, Kentaro Nagamine, Ikkoh Shimizu, Kimihiko Nakajima, Yasunori Terao, Haruka Kusakabe, Makoto Ando, Kentaro Motohara, Lee Spitler

Abstract: In order to shed light on how galactic properties depend on the intergalactic medium (IGM) environment traced by the Ly$α$ forest, we observationally investigate the IGM-galaxy connection using the publicly available 3D IGM tomography data (CLAMATO) and several galaxy catalogs in the COSMOS field. We measure the cross-correlation function (CCF) for $570$ galaxies with spec-$z$ measurements and det… ▽ More In order to shed light on how galactic properties depend on the intergalactic medium (IGM) environment traced by the Ly$α$ forest, we observationally investigate the IGM-galaxy connection using the publicly available 3D IGM tomography data (CLAMATO) and several galaxy catalogs in the COSMOS field. We measure the cross-correlation function (CCF) for $570$ galaxies with spec-$z$ measurements and detect a correlation with the IGM up to $50$ $h^{-1}$ comoving Mpc. We show that galaxies with stellar masses of $10^9-10^{10}$ M$_\odot$ are the dominant contributor to the total CCF signal. We also investigate CCFs for several galaxy populations: Ly$α$ emitters (LAEs), H$α$ emitters (HAEs), [OIII] emitters (O3Es), active galactic nuclei (AGNs), and submillimeter galaxies (SMGs), and we detect the highest signal in AGNs and SMGs at large scales ($r\geq5$ $h^{-1}$ Mpc), but in LAEs at small scales ($r<5$ $h^{-1}$ Mpc). We find that they live in various IGM environments -- HAEs trace the IGM in a similar manner to the continuum-selected galaxies, but LAEs and O3Es tend to reside in higher-density regions. Additionally, LAEs' CCF is flat up to $r\sim3$ $h^{-1}$ Mpc, indicating that they tend to avoid the highest-density regions. For AGNs and SMGs, the CCF peak at $r=5-6$ $h^{-1}$ Mpc implies that they tend to be in locally lower-density regions. We suspect that it is due to the photoionization of IGM HI by AGNs, i.e., the proximity effect. △ Less

Submitted 8 March, 2021; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 27 pages, 16 figures, Accepted for publication on ApJ

Connection Between Galaxies and HI in the Circumgalactic and Intergalactic Media: Variation According to Galaxy Stellar Mass and Star-formation Activity

Authors: Rieko Momose, Ikkoh Shimizu, Kentaro Nagamine, Kazuhiro Shimasaku, Nobunari Kashikawa, Haruka Kusakabe

Abstract: This paper systematically investigates comoving Mpc scale intergalactic medium (IGM) environment around galaxies traced by the Ly$α$ forest. Using our cosmological hydrodynamic simulations, we investigate the IGM-galaxy connection at $z=2$ by two methods: (I) cross-correlation analysis between galaxies and the fluctuation of Ly$α$ forest transmission ($δ_\text{F}$); and (II) comparing the overdens… ▽ More This paper systematically investigates comoving Mpc scale intergalactic medium (IGM) environment around galaxies traced by the Ly$α$ forest. Using our cosmological hydrodynamic simulations, we investigate the IGM-galaxy connection at $z=2$ by two methods: (I) cross-correlation analysis between galaxies and the fluctuation of Ly$α$ forest transmission ($δ_\text{F}$); and (II) comparing the overdensity of neutral hydrogen (HI) and galaxies. Our simulations reproduce observed cross-correlation functions (CCF) between Ly$α$ forest and Lyman-break galaxies. We further investigate the variation of the CCF using subsamples divided by dark matter halo mass ($M_\text{DH}$), galaxy stellar mass ($M_\star$), and star-formation rate (SFR), and find that the CCF signal becomes stronger with increasing $M_\text{DH}$, $M_\star$, and SFR. The CCFs between galaxies and gas-density fluctuation are also found to have similar trends. Therefore, the variation of the $δ_\text{F}$-CCF depending on $M_\text{DH}$, $M_\star$, and SFR is due to varying gas density around galaxies. We find that the correlation between galaxies and the IGM HI distribution strongly depends on $M_\text{DH}$ as expected from the linear theory. Our results support the $Λ$CDM paradigm, finding a spatial correlation between galaxies and IGM HI, with more massive galaxies being clustered in higher-density regions. △ Less

Submitted 8 March, 2021; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 21 pages, 9 figures, Accepted for publication on ApJ

ALMA uncovers the [CII] emission and warm dust continuum in a z = 8.31 Lyman break galaxy

Authors: Tom J. L. C. Bakx, Yoichi Tamura, Takuya Hashimoto, Akio K. Inoue, Minju M. Lee, Ken Mawatari, Kazuaki Ota, Hideki Umehata, Erik Zackrisson, Bunyo Hatsukade, Kotaro Kohno, Yuichi Matsuda, Hiroshi Matsuo, Takashi Okamoto, Takatoshi Shibuya, Ikkoh Shimizu, Yoshiaki Taniguchi, Naoki Yoshida

Abstract: We report on the detection of the [CII] 157.7 $μ$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $μ$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emiss… ▽ More We report on the detection of the [CII] 157.7 $μ$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $μ$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $μ$m dust emission (90 $μ$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $μ$m rest-frame) down to 18 $μ$Jy (3$σ$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $μ$Jy continuum emission at 850 $μ$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($β_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $β_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization. △ Less

Submitted 18 February, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 15 pages, 8 figures, Accepted by MNRAS

Balmer breaks in simulated galaxies at z>6

Authors: C. Binggeli, E. Zackrisson, X. Ma, A. K. Inoue, A. Vikaeus, T. Hashimoto, K. Mawatari, I. Shimizu, D. Ceverino

Abstract: Photometric observations of the spectroscopically confirmed $z\approx 9.1$ galaxy MACS1149-JD1 have indicated the presence of a prominent Balmer break in its spectral energy distribution, which may be interpreted as due to very large fluctuations in its past star formation activity. In this paper, we investigate to what extent contemporary simulations of high-redshift galaxies produce star formati… ▽ More Photometric observations of the spectroscopically confirmed $z\approx 9.1$ galaxy MACS1149-JD1 have indicated the presence of a prominent Balmer break in its spectral energy distribution, which may be interpreted as due to very large fluctuations in its past star formation activity. In this paper, we investigate to what extent contemporary simulations of high-redshift galaxies produce star formation rate variations sufficiently large to reproduce the observed Balmer break of MACS1149-JD1. We find that several independent galaxy simulations are unable to account for Balmer breaks of the inferred size, suggesting that MACS1149-JD1 either must be a very rare type of object or that our simulations are missing some key ingredient. We present predictions of spectroscopic Balmer break strength distributions for $z\approx 7-9$ galaxies that may be tested through observations with the upcoming James Webb Space Telescope and also discuss the impact that various assumptions on dust reddening, Lyman continuum leakage and deviations from a standard stellar initial mass function would have on the results. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: 11 pages, 8 figures. Accepted for publication in MNRAS

Bubble mapping with the Square Kilometer Array -- I. Detecting galaxies with Euclid, JWST, WFIRST and ELT within ionized bubbles in the intergalactic medium at z>6

Authors: Erik Zackrisson, Suman Majumdar, Rajesh Mondal, Christian Binggeli, Martin Sahlén, Tirthankar Roy Choudhury, Benedetta Ciardi, Abhirup Datta, Kanan K. Datta, Pratika Dayal, Andrea Ferrara, Sambit K. Giri, Umberto Maio, Sangeeta Malhotra, Garrelt Mellema, Andrei Mesinger, James Rhoads, Claes-Erik Rydberg, Ikkoh Shimizu

Abstract: The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be cente… ▽ More The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts. △ Less

Submitted 16 February, 2020; v1 submitted 1 May, 2019; originally announced May 2019.

Comments: 18 pages, 7 figures, v.2: Accepted for publications in MNRAS

Osaka Feedback Model: Isolated Disk Galaxy Simulations

Authors: Ikkoh Shimizu, Keita Todoroki, Hidenobu Yajima, Kentaro Nagamine

Abstract: We study various implementations of supernova feedback model and present the results of our `Osaka feedback model' using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {\small GADGET-3}. Our model is a modified version of Stinson et al.'s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical… ▽ More We study various implementations of supernova feedback model and present the results of our `Osaka feedback model' using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {\small GADGET-3}. Our model is a modified version of Stinson et al.'s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical properties of SN remnants from the results of Chevalier and McKee \& Ostriker, such as the effective radius of SN bubble and the remnant life-time, in the form of Sedov-Taylor (ST)-like solutions with the effect of radiative cooling. Our model utilizes the local, physical parameters such as density and temperature of the ISM rather than galactic or halo properties to determine the galactic wind velocity or mass-loading factor. The Osaka model succeeds in self-regulating star formation, and naturally produces galactic outflow with variable velocities depending on the local environment and available SN energy as a function of time.An important addition to our previous work by Aoyama et al. is the implementation of the {\small CELib} chemistry library which allows us to deal with the time-dependent input of energy and metal yields for type Ia \& II supernovae (SNe) and asymptotic giant branch (AGB) stars. As initial tests of our model, we apply it to isolated galaxy simulations, and examine various galactic properties and compare with observational data including metal abundances. △ Less

Submitted 12 January, 2019; originally announced January 2019.

Comments: 25 pages, 34 figures, accepeted for publication in MNRAS

Dust scaling relations in a cosmological simulation

Authors: Kuan-Chou Hou, Shohei Aoyama, Hiroyuki Hirashita, Kentaro Nagamine, Ikkoh Shimizu

Abstract: To study the dust evolution in the cosmological structure formation history, we perform a smoothed particle hydrodynamic simulation with a dust enrichment model in a cosmological volume. We adopt the dust evolution model that represents the grain size distribution by two sizes and takes into account stellar dust production and interstellar dust processing. We examine the dust mass function and the… ▽ More To study the dust evolution in the cosmological structure formation history, we perform a smoothed particle hydrodynamic simulation with a dust enrichment model in a cosmological volume. We adopt the dust evolution model that represents the grain size distribution by two sizes and takes into account stellar dust production and interstellar dust processing. We examine the dust mass function and the scaling properties of dust in terms of the characteristics of galaxies. The simulation broadly reproduces the observed dust mass functions at redshift $z = 0$, except that it overproduces the massive end at dust mass $M_\mathrm{d} \gtrsim 10^{8}$ ${\rm M}_\odot$. This overabundance is due to overproducing massive gas/metal-rich systems, but we also note that the relation between stellar mass and gas-phase metallicity is reproduced fairly well by our recipe. The relation between dust-to-gas ratio and metallicity shows a good agreement with the observed one at $z=0$, which indicates successful implementation of dust evolution in our cosmological simulation. Star formation consumes not only gas but also dust, causing a decreasing trend of the dust-to-stellar mass ratio at the high-mass end of galaxies. We also examine the redshift evolution up to $z \sim~ 5$, and find that the galaxies have on average the highest dust mass at $z = 1-2$. For the grain size distribution, we find that galaxies with metallicity $\sim 0.3~ Z_\odot$ tend to have the highest small-to-large grain abundance ratio; consequently, the extinction curves in those galaxies have the steepest ultraviolet slopes. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: 18 pages, 9 figures, accepted for publication in MNRAS

Comparison of cosmological simulations and deep submillimetre galaxy surveys

Authors: Shohei Aoyama, Hiroyuki Hirashita, Chen-Fatt Lim, Yu-Yen Chang, Wei-Hao Wang, Kentaro Nagamine, Kuan-Chou Hou, Ikkoh Shimizu, Hui-Hsuan Chung, Chien-Hsiu Lee, Xian-Zhong Zheng

Abstract: Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of… ▽ More Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of infrared (IR)-emitting region, IR luminosity function and IRX--$β$ relation. After confirming that our models reproduce the local galaxy properties, we examine the STUDIES sample at $z\approx 1-4$, finding that the simulation reproduces the aforementioned quantities except for the $z\gtrsim 2$ IR luminosity function at the massive end ($\sim 10^{13}$ L$_{\odot}$). This means that the current simulation correctly reproduces the overall scaling between the size and luminosity (or star formation rate) of dusty region, but lacks extreme starburst phenomena at $z\gtrsim 2$. We also discuss extinction curves and possible AGN contribution. △ Less

Submitted 21 January, 2019; v1 submitted 27 September, 2018; originally announced September 2018.

Comments: 13 pages, 10 figures, published in MNRAS

Redshift Space Distortion of 21cm line at 1<z<5 with Cosmological Hydrodynamic Simulations

Authors: Rika Ando, Atsushi J. Nishizawa, Kenji Hasegawa, Ikko Shimizu, Kentaro Nagamine

Abstract: We measure the scale dependence and redshift dependence of 21 cm line emitted from the neutral hydrogen gas at redshift 1<z<5 using full cosmological hydrodynamic simulations by taking the ratios between the power spectra of HI-dark matter cross correlation and dark matter auto-correlation. The neutral hydrogen distribution is computed in full cosmological hydrodynamic simulations including star f… ▽ More We measure the scale dependence and redshift dependence of 21 cm line emitted from the neutral hydrogen gas at redshift 1<z<5 using full cosmological hydrodynamic simulations by taking the ratios between the power spectra of HI-dark matter cross correlation and dark matter auto-correlation. The neutral hydrogen distribution is computed in full cosmological hydrodynamic simulations including star formation and supernova feedback under a uniform ultra-violet background radiation. We find a significant scale dependence of HI bias at z>3 on scales of k>1h/Mpc, but it is roughly constant at lower redshift z<3. The redshift evolution of HI bias is relatively slow compared to that of QSOs at similar redshift range. We also measure a redshift space distortion (RSD) of HI gas to explore the properties of HI clustering. Fitting to a widely applied theoretical prediction, we find that the constant bias is consistent with that measured directly from the real-space power spectra, and the velocity dispersion is marginally consistent with the linear perturbation prediction. Finally we compare the results obtained from our simulation and the Illustris simulation, and conclude that the detailed astrophysical effects do not affect the scale dependence of HI bias very much, which implies that the cosmological analysis using 21 cm line of HI will be robust against the uncertainties arising from small-scale astrophysical processes such as star formation and supernova feedback. △ Less

Submitted 28 January, 2019; v1 submitted 3 August, 2018; originally announced August 2018.

Comments: 10 pages, 6 figures, submitted to MNRAS

Detection of the Far-infrared [O III] and Dust Emission in a Galaxy at Redshift 8.312: Early Metal Enrichment in the Heart of the Reionization Era

Authors: Yoichi Tamura, Ken Mawatari, Takuya Hashimoto, Akio K. Inoue, Erik Zackrisson, Lise Christensen, Christian Binggeli, Yuichi Matsuda, Hiroshi Matsuo, Tsutomu T. Takeuchi, Ryosuke S. Asano, Kaho Sunaga, Ikkoh Shimizu, Takashi Okamoto, Naoki Yoshida, Minju Lee, Takatoshi Shibuya, Yoshiaki Taniguchi, Hideki Umehata, Bunyo Hatsukade, Kotaro Kohno, Kazuaki Ota

Abstract: We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $μ$m line and rest-frame 90 $μ$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthes… ▽ More We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $μ$m line and rest-frame 90 $μ$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $μ$m flux density of $137 \pm 26$ $μ$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $β= 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($τ_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $τ_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $τ_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass. △ Less

Submitted 11 February, 2019; v1 submitted 11 June, 2018; originally announced June 2018.

Comments: Accepted by ApJ. 18 pages, 10 figures, 5 tables

"Big Three Dragons": a z = 7.15 Lyman BreakGalaxy Detected in [OIII] 88 $μ$m, [CII] 158 $μ$m, and Dust Continuum with ALMA

Authors: Takuya Hashimoto, Akio K. Inoue, Ken Mawatari, Yoichi Tamura, Hiroshi Matsuo, Hisanori Furusawa, Yuichi Harikane, Takatoshi Shibuya, Kirsten K. Knudsen, Kotaro Kohno, Yoshiaki Ono, Erik Zackrisson, Takashi Okamoto, Nobunari Kashikawa, Pascal A. Oesch, Masami Ouchi, Kazuaki Ota, Ikkoh Shimizu, Yoshiaki Taniguchi, Hideki Umehata, Darach Watson

Abstract: We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$α$ with a small rest-frame equivalent width of $\approx4$ Å. Previous HST image has shown that the target is comprised of two spatially separated clumps in the… ▽ More We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$α$ with a small rest-frame equivalent width of $\approx4$ Å. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $μ$m, [CII] 158 $μ$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $β_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio. △ Less

Submitted 8 April, 2019; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: 24 pages, 13 figures, 9 tables. Accepted for publication in PASJ. In addition to the dust continuum detection in ALMA Band 6 reported in the original manuscript, we newly report the dust continuum detection in ALMA Band 8 based on our new ALMA Cycle 5 observations

The distribution and physical properties of high-redshift [OIII] emitters in a cosmological hydrodynamics simulation

Authors: Kana Moriwaki, Naoki Yoshida, Ikkoh Shimizu, Yuichi Harikane, Yuichi Matsuda, Hiroshi Matsuo, Takuya Hashimoto, Akio K. Inoue, Yoichi Tamura, Tohru Nagao

Abstract: Recent observations with the Atacama Large Millimeter/submillimeter Array (ALMA) detected far-infrared emission lines such as the [OIII] 88 μm line from galaxies at $z \sim 7 - 9$. Far-infrared lines can be used to probe the structure and kinematics of such high-redshift galaxies as well as to accurately determine their spectroscopic redshifts. We use a cosmological simulation of galaxy formation… ▽ More Recent observations with the Atacama Large Millimeter/submillimeter Array (ALMA) detected far-infrared emission lines such as the [OIII] 88 μm line from galaxies at $z \sim 7 - 9$. Far-infrared lines can be used to probe the structure and kinematics of such high-redshift galaxies as well as to accurately determine their spectroscopic redshifts. We use a cosmological simulation of galaxy formation to study the physical properties of [OIII] 88 μm emitters. In a comoving volume of 50 $h^{-1}$ Mpc on a side, we locate 34 galaxies with stellar masses greater than $10^8\ {\rm M_\odot}$ at $z = 9$, and more than 270 such galaxies at $z = 7$. We calculate the [OIII] 88 μm luminosities ($L_{\rm OIII}$) by combining a physical model of HII regions with emission line calculations using the photoionization code CLOUDY. We show that the resulting $L_{\rm OIII}$, for a given star formation rate, is slightly higher than predicted from the empirical relation for local galaxies, and is consistent with recent observations of galaxies at redshifts 7 - 9. Bright [OIII] emitters with $L_{\rm OIII} > 10^8 {\rm L_\odot}$ have stellar masses greater than $10^9\ {\rm M_\odot}$, star formation rates higher than $3\ {\rm M_\odot\ yr}^{-1}$, and the typical metallicity is $\sim 0.1\ {\rm Z_\odot}$. The galaxies are hosted by dark matter halos with masses greater than $10^{10.5}\ {\rm M_\odot}$. Massive galaxies show characteristic structure where the [OIII] emitting gas largely overlaps with young stars, but the emission peak is separated from the main stellar population, suggesting the stochastic and localized nature of star formation in the first galaxies. We propose to use the [OIII] 5007 Å line, to be detected by James Webb Space Telescope (JWST), to study the properties of galaxies whose [OIII] 88 μm line emission has been already detected with ALMA. △ Less

Submitted 4 September, 2018; v1 submitted 18 May, 2018; originally announced May 2018.

Comments: 5 pages, 4 figures. Accepted for publication in MN Letters

The onset of star formation 250 million years after the Big Bang

Authors: Takuya Hashimoto, Nicolas Laporte, Ken Mawatari, Richard S. Ellis, Akio. K. Inoue, Erik Zackrisson, Guido Roberts-Borsani, Wei Zheng, Yoichi Tamura, Franz E. Bauer, Thomas Fletcher, Yuichi Harikane, Bunyo Hatsukade, Natsuki H. Hayatsu, Yuichi Matsuda, Hiroshi Matsuo, Takashi Okamoto, Masami Ouchi, Roser Pello, Claes-Erik Rydberg, Ikkoh Shimizu, Yoshiaki Taniguchi, Hideki Umehata, Naoki Yoshida

Abstract: A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies migh… ▽ More A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: To appear in the 17 May issue of Nature

Journal ref: Proc. IAU 15 (2019) 221-225

Hyper-Kamiokande Design Report

Authors: Hyper-Kamiokande Proto-Collaboration, :, K. Abe, Ke. Abe, H. Aihara, A. Aimi, R. Akutsu, C. Andreopoulos, I. Anghel, L. H. V. Anthony, M. Antonova, Y. Ashida, V. Aushev, M. Barbi, G. J. Barker, G. Barr, P. Beltrame, V. Berardi, M. Bergevin, S. Berkman, L. Berns, T. Berry, S. Bhadra, D. Bravo-Berguño, F. d. M. Blaszczyk , et al. (291 additional authors not shown)

Abstract: On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from th… ▽ More On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation. △ Less

Submitted 28 November, 2018; v1 submitted 9 May, 2018; originally announced May 2018.

Comments: 325 pages

Lyman continuum leakage versus quenching with the James Webb Space Telescope: The spectral signatures of quenched star formation activity in reionization-epoch galaxies

Authors: C. Binggeli, E. Zackrisson, K. Pelckmans, R. Cubo, H. Jensen, I. Shimizu

Abstract: In this paper, we study the effects of a recent drop in star formation rate (SFR) on the spectra of epoch of reionization (EoR) galaxies, and the resulting degeneracy with the spectral features produced by extreme Lyman continuum leakage. In order to study these effects in the wavelength range relevant for the upcoming James Webb Space Telescope (JWST), we utilize synthetic spectra of simulated Eo… ▽ More In this paper, we study the effects of a recent drop in star formation rate (SFR) on the spectra of epoch of reionization (EoR) galaxies, and the resulting degeneracy with the spectral features produced by extreme Lyman continuum leakage. In order to study these effects in the wavelength range relevant for the upcoming James Webb Space Telescope (JWST), we utilize synthetic spectra of simulated EoR galaxies from cosmological simulations together with synthetic spectra of partially quenched mock galaxies. We find that rapid declines in the SFR of EoR galaxies could seriously affect the applicability of methods that utilize the equivalent width of Balmer lines and the ultraviolet spectral slope to assess the escape fraction of EoR galaxies. In order to determine if the aforementioned degeneracy can be avoided by using the overall shape of the spectrum, we generate mock NIRCam observations and utilize a classification algorithm to identify galaxies that have undergone quenching. We find that while there are problematic cases, JWST/NIRCam or NIRSpec should be able to reliably identify galaxies with redshifts $z\sim 7$ that have experienced a significant decrease in the SFR (by a factor 10-100) in the past 50-100 Myr with a success rate $\gtrsim 85\%$. We also find that uncertainties in the dust-reddening effects on EoR galaxies significantly affect the performance of the results of the classification algorithm. We argue that studies that aim to characterize the dust extinction law most representative in the EoR would be extremely useful. △ Less

Submitted 24 July, 2018; v1 submitted 24 April, 2018; originally announced April 2018.

Comments: 10 pages, 7 figures. Accepted for publication in MNRAS

Cosmological simulation with dust formation and destruction

Authors: Shohei Aoyama, Kuan-Chou Hou, Hiroyuki Hirashita, Kentaro Nagamine, Ikkoh Shimizu

Abstract: To investigate the evolution of dust in a cosmological volume, we perform hydrodynamic simulations, in which the enrichment of metals and dust is treated self-consistently with star formation and stellar feedback. We consider dust evolution driven by dust production in stellar ejecta, dust destruction by sputtering, grain growth by accretion and coagulation, and grain disruption by shattering, and… ▽ More To investigate the evolution of dust in a cosmological volume, we perform hydrodynamic simulations, in which the enrichment of metals and dust is treated self-consistently with star formation and stellar feedback. We consider dust evolution driven by dust production in stellar ejecta, dust destruction by sputtering, grain growth by accretion and coagulation, and grain disruption by shattering, and treat small and large grains separately to trace the grain size distribution. After confirming that our model nicely reproduces the observed relation between dust-to-gas ratio and metallicity for nearby galaxies, we concentrate on the dust abundance over the cosmological volume in this paper. The comoving dust mass density has a peak at redshift $z\sim 1$--2, coincident with the observationally suggested dustiest epoch in the Universe. {In the local Universe}, roughly 10 per cent of the dust is contained in the intergalactic medium (IGM), where only 1/3--1/4 of the dust survives against dust destruction by sputtering. We also show that the dust mass function is roughly reproduced at $\lesssim 10^8$ M$_\odot$, while the massive end still has a discrepancy, which indicates {the necessity of stronger feedback in massive galaxies}. %%The relation showed that accretion is essential for dusty galaxies. In addition, our model broadly reproduces the observed radial profile of dust surface density in the circum-galactic medium (CGM). While our model satisfies the observational constraints for the dust extinction {on cosmological scales}, it predicts that the dust in the CGM and IGM is dominated by large ($> 0.03~μ$m) grains, which is in tension with the steep reddening curves {observed} in the CGM. △ Less

Submitted 22 July, 2018; v1 submitted 12 February, 2018; originally announced February 2018.

Comments: 17 pages, 9 figures. published on MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2018, Volume 478, Issue 4, Pages 4905-4921

SILVERRUSH. VI. A simulation of Ly$α$ emitters in the reionization epoch and a comparison with Subaru Hyper Suprime-Cam survey early data

Authors: Akio K. Inoue, Kenji Hasegawa, Tomoaki Ishiyama, Hidenobu Yajima, Ikkoh Shimizu, Masayuki Umemura, Akira Konno, Yuichi Harikane, Takatoshi Shibuya, Masami Ouchi, Kazuhiro Shimasaku, Yoshiaki Ono, Haruka Kusakabe, Ryo Higuchi, Chien-Hsiu Lee

Abstract: The survey of Lyman $α$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes… ▽ More The survey of Lyman $α$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined $2^3$ models depending on the presence or absence of dispersion of halo Ly$α$ emissivity, dispersion of the halo Ly$α$ optical depth, $τ_α$, and halo mass dependence of $τ_α$. The unique free parameter in our model, a pivot value of $τ_α$, is calibrated so as to reproduce the $z=5.7$ Ly$α$ luminosity function (LF). We compare our model predictions with Ly$α$ LFs at $z=6.6$ and $7.3$, LAE angular auto-correlation functions (ACFs) at $z=5.7$ and $6.6$, and LAE fractions in Lyman break galaxies at $5<z<7$. The Ly$α$ LFs and ACFs are reproduced by multiple models, but the LAE fraction turns out to be the most critical test. The dispersion of $τ_α$ and the halo mass dependence of $τ_α$ are essential to explain all observations reasonably. Therefore, a simple model of one-to-one correspondence between halo mass and Ly$α$ luminosity with a constant Ly$α$ escape fraction has been ruled out. Based on our best model, we present a formula to estimate the intergalactic neutral hydrogen fraction, $x_{\rm HI}$, from the observed Ly$α$ luminosity density at $z\gtrsim6$. We finally obtain $x_{\rm HI}=0.5_{-0.3}^{+0.1}$ as a volume-average at $z=7.3$. △ Less

Submitted 3 April, 2018; v1 submitted 29 December, 2017; originally announced January 2018.

Comments: PASJ in press, 30 pages, 20 figures

Populating H$_2$ and CO in galaxy simulation with dust evolution

Authors: Li-Hsin Chen, Hiroyuki Hirashita, Kuan-Chou Hou, Shohei Aoyama, Ikkoh Shimizu, Kentaro Nagamine

Abstract: There are two major theoretical issues for the star formation law (the relation between the surface densities of molecular gas and star formation rate on a galaxy scale): (i) At low metallicity, it is not obvious that star-forming regions are rich in H$_2$ because the H$_2$ formation rate depends on the dust abundance; and (ii) whether or not CO really traces H$_2$ is uncertain, especially at low… ▽ More There are two major theoretical issues for the star formation law (the relation between the surface densities of molecular gas and star formation rate on a galaxy scale): (i) At low metallicity, it is not obvious that star-forming regions are rich in H$_2$ because the H$_2$ formation rate depends on the dust abundance; and (ii) whether or not CO really traces H$_2$ is uncertain, especially at low metallicity. To clarify these issues, we use a hydrodynamic simulation of an isolated disc galaxy with a spatial resolution of a few tens parsecs. The evolution of dust abundance and grain size distribution is treated consistently with the metal enrichment and the physical state of the interstellar medium. We compute the H$_2$ and CO abundances using a subgrid post-processing model based on the dust abundance and the dissociating radiation field calculated in the simulation. We find that when the metallicity is $\lesssim 0.4$ Z$_\odot$ ($t<1$ Gyr), H$_2$ is not a good tracer of star formation rate because H$_2$-rich regions are limited to dense compact regions. At $Z\gtrsim 0.8$ Z$_\odot$, a tight star formation law is established for both H$_2$ and CO. At old ($t \sim 10$ Gyr) ages, we also find that adopting the so-called MRN grain size distribution with an appropriate dust-to-metal ratio over the entire disc gives reasonable estimates for the H$_2$ and CO abundances. For CO, improving the spatial resolution of the simulation is important while the H$_2$ abundance is not sensitive to sub-resolution structures at $Z\gtrsim 0.4$ Z$_\odot$. △ Less

Submitted 1 November, 2017; originally announced November 2017.

Comments: 21 pages, 14 figures, accepted for publication in MNRAS

Evolution of dust extinction curves in galaxy simulation

Authors: Kuan-Chou Hou, Hiroyuki Hirashita, Kentaro Nagamine, Shohei Aoyama, Ikkoh Shimizu

Abstract: To understand the evolution of extinction curve, we calculate the dust evolution in a galaxy using smoothed particle hydrodynamics simulations incorporating stellar dust production, dust destruction in supernova shocks, grain growth by accretion and coagulation, and grain disruption by shattering. The dust species are separated into carbonaceous dust and silicate. The evolution of grain size distr… ▽ More To understand the evolution of extinction curve, we calculate the dust evolution in a galaxy using smoothed particle hydrodynamics simulations incorporating stellar dust production, dust destruction in supernova shocks, grain growth by accretion and coagulation, and grain disruption by shattering. The dust species are separated into carbonaceous dust and silicate. The evolution of grain size distribution is considered by dividing grain population into large and small gains, which allows us to estimate extinction curves. We examine the dependence of extinction curves on the position, gas density, and metallicity in the galaxy, and find that extinction curves are flat at $t \lesssim 0.3$ Gyr because stellar dust production dominates the total dust abundance. The 2175 Å bump and far-ultraviolet (FUV) rise become prominent after dust growth by accretion. At $t \gtrsim 3$ Gyr, shattering works efficiently in the outer disc and low density regions, so extinction curves show a very strong 2175 Å bump and steep FUV rise. The extinction curves at $t\gtrsim 3$ Gyr are consistent with the Milky Way extinction curve, which implies that we successfully included the necessary dust processes in the model. The outer disc component caused by stellar feedback has an extinction curves with a weaker 2175 Å bump and flatter FUV slope. The strong contribution of carbonaceous dust tends to underproduce the FUV rise in the Small Magellanic Cloud extinction curve, which supports selective loss of small carbonaceous dust in the galaxy. The snapshot at young ages also explain the extinction curves in high-redshift quasars. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test

Authors: Ji-hoon Kim, Oscar Agertz, Romain Teyssier, Michael J. Butler, Daniel Ceverino, Jun-Hwan Choi, Robert Feldmann, Ben W. Keller, Alessandro Lupi, Thomas Quinn, Yves Revaz, Spencer Wallace, Nickolay Y. Gnedin, Samuel N. Leitner, Sijing Shen, Britton D. Smith, Robert Thompson, Matthew J. Turk, Tom Abel, Kenza S. Arraki, Samantha M. Benincasa, Sukanya Chakrabarti, Colin DeGraf, Avishai Dekel, Nathan J. Goldbaum , et al. (18 additional authors not shown)

Abstract: Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the st… ▽ More Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes. △ Less

Submitted 8 March, 2018; v1 submitted 10 October, 2016; originally announced October 2016.

Comments: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at http://www.jihoonkim.org/agora/AGORA_Paper4_draft.pdf, The first paper of the AGORA Initiative is at http://adsabs.harvard.edu/abs/2014ApJS..210...14K, More information on AGORA is at http://www.AGORAsimulations.org/

Journal ref: Astrophysical Journal, 833, 202, 2016

Galaxy Simulation with Dust Formation and Destruction

Authors: Shohei Aoyama, Kuan-Chou Hou, Ikkoh Shimizu, Hiroyuki Hirashita, Keita Todoroki, Jun-Hwan Choi, Kentaro Nagamine

Abstract: We perform smoothed particle hydrodynamics (SPH) simulations of an isolated galaxy with a new treatment for dust formation and destruction. To this aim, we treat dust and metal production self-consistently with star formation and supernova feedback. For dust, we consider a simplified model of grain size distribution by representing the entire range of grain sizes with large and small grains. We in… ▽ More We perform smoothed particle hydrodynamics (SPH) simulations of an isolated galaxy with a new treatment for dust formation and destruction. To this aim, we treat dust and metal production self-consistently with star formation and supernova feedback. For dust, we consider a simplified model of grain size distribution by representing the entire range of grain sizes with large and small grains. We include dust production in stellar ejecta, dust destruction by supernova (SN) shocks, grain growth by accretion and coagulation, and grain disruption by shattering. We find that the assumption of fixed dust-to-metal mass ratio becomes no longer valid when the galaxy is older than 0.2 Gyr, at which point the grain growth by accretion starts to contribute to the nonlinear rise of dust-to-gas ratio. As expected in our previous one-zone model, shattering triggers grain growth by accretion since it increases the total surface area of grains. Coagulation becomes significant when the galaxy age is greater than $\sim$ 1 Gyr: at this epoch the abundance of small grains becomes high enough to raise the coagulation rate of small grains. We further compare the radial profiles of dust-to-gas ratio $(\mathcal{D})$ and dust-to-metal ratio $(\mathcal{D}/Z)$ (i.e., depletion) at various ages with observational data. We find that our simulations broadly reproduce the radial gradients of dust-to-gas ratio and depletion. In the early epoch ($\lesssim 0.3$ Gyr), the radial gradient of $\mathcal{D}$ follows the metallicity gradient with $\mathcal{D}/Z$ determined by the dust condensation efficiency in stellar ejecta, while the $\mathcal{D}$ gradient is steeper than the $Z$ gradient at the later epochs because of grain growth by accretion. The framework developed in this paper is applicable to any SPH-based galaxy evolution simulations including cosmological ones. △ Less

Submitted 18 January, 2017; v1 submitted 23 September, 2016; originally announced September 2016.

Comments: 18 pages, 16 figures, published on MNRAS

Journal ref: MNRAS (2017) 466 (1): 105-121

The spectral evolution of the first Galaxies. III. Simulated James Webb Space Telescope spectra of reionization-epoch galaxies with Lyman continuum leakage

Authors: E. Zackrisson, C. Binggeli, K. Finlator, N. Y. Gnedin, J. -P. Paardekooper, I. Shimizu, A. K. Inoue, H. Jensen, G. Micheva, S. Khochfar, C. Dalla Vecchia

Abstract: Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta dia… ▽ More Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 Å to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement. △ Less

Submitted 12 January, 2017; v1 submitted 29 August, 2016; originally announced August 2016.

Comments: 12 pages, 8 figures, v.2: ApJ, accepted. Model grids are available from http://www.astro.uu.se/~ez/lycan/lycan.html

Search for electron antineutrinos associated with gravitational wave events GW150914 and GW151226 using KamLAND

Authors: KamLAND Collaboration, A. Gando, Y. Gando, T. Hachiya, A. Hayashi, S. Hayashida, H. Ikeda, K. Inoue, K. Ishidoshiro, Y. Karino, M. Koga, S. Matsuda, T. Mitsui, K. Nakamura, S. Obara, T. Oura, H. Ozaki, I. Shimizu, Y. Shirahata, J. Shirai, A. Suzuki, T. Takai, K. Tamae, Y. Teraoka, K. Ueshima , et al. (23 additional authors not shown)

Abstract: We present a search for low energy antineutrino events coincident with the gravitational wave events GW150914 and GW151226, and the candidate event LVT151012 using KamLAND, a kiloton-scale antineutrino detector. We find no inverse beta-decay neutrino events within $\pm 500$ seconds of either gravitational wave signal. This non-detection is used to constrain the electron antineutrino fluence and th… ▽ More We present a search for low energy antineutrino events coincident with the gravitational wave events GW150914 and GW151226, and the candidate event LVT151012 using KamLAND, a kiloton-scale antineutrino detector. We find no inverse beta-decay neutrino events within $\pm 500$ seconds of either gravitational wave signal. This non-detection is used to constrain the electron antineutrino fluence and the luminosity of the astrophysical sources. △ Less

Submitted 3 October, 2016; v1 submitted 22 June, 2016; originally announced June 2016.

Comments: 11 pages, 4 figures, as published in ApJL. Updated to replace power law spectrum with Fermi Dirac spectrum

Journal ref: Astrophys.J.829, L34, 2016

Detection of an oxygen emission line from a high redshift galaxy in the reionization epoch

Authors: Akio K. Inoue, Yoichi Tamura, Hiroshi Matsuo, Ken Mawatari, Ikkoh Shimizu, Takatoshi Shibuya, Kazuaki Ota, Naoki Yoshida, Erik Zackrisson, Nobunari Kashikawa, Kotaro Kohno, Hideki Umehata, Bunyo Hatsukade, Masanori Iye, Yuichi Matsuda, Takashi Okamoto, Yuki Yamaguchi

Abstract: The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abun… ▽ More The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: First Release on June 16, 2016 in Science