Showing 1–50 of 109 results for author: Okamoto, T

The evolution of supermassive blackhole mass--bulge mass relation by a semi-analytic model, $ν^2$GC

Authors: Tatsuki Shimizu, Taira Oogi, Takashi Okamoto, Masahiro Nagashima, Motohiro Enoki

Abstract: We have investigated the redshift evolution of the relationship between supermassive black hole (SMBH) mass and host bulge mass using a semi-analytical galaxy formation model $ν^2$GC. Our model reproduces the relation in the local universe well. We find that, at high redshift ($z \gtrsim 3$), two sequences appear in the SMBH mass--bulge mass plane. The emergence of these two sequences can be attri… ▽ More We have investigated the redshift evolution of the relationship between supermassive black hole (SMBH) mass and host bulge mass using a semi-analytical galaxy formation model $ν^2$GC. Our model reproduces the relation in the local universe well. We find that, at high redshift ($z \gtrsim 3$), two sequences appear in the SMBH mass--bulge mass plane. The emergence of these two sequences can be attributed to the primary triggers of the growth of the SMBHs and bulges: galaxy mergers and disc instabilities. The growth of SMBHs and bulges as a result of galaxy mergers is responsible for giving rise to the high-mass sequence, in which SMBHs are more massive for a given host bulge mass than in the low-mas sequence. Conversely, disc instabilities are accountable for the emergence of the low-mass sequence. At lower redshifts, galaxy mergers tend to become increasingly deficient in gas, resulting in a preferential increase of bulge mass without a corresponding growth in SMBH mass. This has the effect of causing galaxies in the upper sequence to shift towards the lower one on the SMBH mass-bulge mass plane. The galaxies that undergo dry mergers serve to bridge the gap between the two sequences, eventually leading to convergence into a single relation known in the local universe. Our results suggest that the observations of the SMBH mass-bulge mass relation in high redshifts can provide insight into their growth mechanisms. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 7 pages, 7 figure, accepted for publication in MNRAS

Journal ref: MNRAS, 531, 851 (2024)

Cloud-cloud collisions triggering star formation in galaxy simulations

Authors: Shu Horie, Takashi Okamoto, Asao Habe

Abstract: Cloud-cloud collisions (CCCs) are expected to compress gas and trigger star formation. However, it is not well understood how the collisions and the induced star formation affect galactic-scale properties. By developing an on-the-fly algorithm to identify CCCs at each timestep in a galaxy simulation and a model that relates CCC-triggered star formation to collision speeds, we perform simulations o… ▽ More Cloud-cloud collisions (CCCs) are expected to compress gas and trigger star formation. However, it is not well understood how the collisions and the induced star formation affect galactic-scale properties. By developing an on-the-fly algorithm to identify CCCs at each timestep in a galaxy simulation and a model that relates CCC-triggered star formation to collision speeds, we perform simulations of isolated galaxies to study the evolution of galaxies and giant molecular clouds (GMCs) with prescriptions of self-consistent CCC-driven star formation and stellar feedback. We find that the simulation with the CCC-triggered star formation produces slightly higher star formation rates and a steeper Kennicutt-Schmidt relation than that with a more standard star formation recipe, although collision speeds and frequencies are insensitive to the star formation models. In the simulation with the CCC model, about 70 per cent of the stars are born via CCCs, and colliding GMCs with masses of $\approx 10^{5.5}\,M_{\odot}$ are the main drivers of CCC-driven star formation. In the simulation with the standard star formation recipe, about 50 per cent of stars are born in colliding GMCs even without the CCC-triggered star formation model. These results suggest that CCCs may be one of the most important star formation processes in galaxy evolution. Furthermore, we find that a post-processing analysis of CCCs, as used in previous studies in galaxy simulations, may lead to slightly greater collision speeds and significantly lower collision frequencies than the on-the-fly analysis. △ Less

Submitted 8 March, 2024; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: 18 pages, 16 figures, published in MNRAS

Journal ref: MNRAS, 527, 10077 (2024)

The impact of bar origin and morphology on stellar migration

Authors: Elizabeth J. Iles, Alex R. Pettitt, Takashi Okamoto, Daisuke Kawata

Abstract: Different mechanisms driving bar structure formation indicate that bar origins should be distinguishable in the stellar populations of galaxies. To study how these origins affect different bar morphologies and impact stellar orbits and migration, we analyse three simulated discs which are representative of bar formation under isolated evolution motivated by disc instability, and interaction driven… ▽ More Different mechanisms driving bar structure formation indicate that bar origins should be distinguishable in the stellar populations of galaxies. To study how these origins affect different bar morphologies and impact stellar orbits and migration, we analyse three simulated discs which are representative of bar formation under isolated evolution motivated by disc instability, and interaction driven tidal development. The first isolated disc and the tidally driven disc produce similar bar structure, while the second isolated disc, generated by the tidal initial condition without the companion, is visibly dissimilar. Changes to radial and vertical positions, angular momentum in the disc-plane, orbital eccentricity and the subsequent disc metallicities are assessed, as is the dependence on stellar age and formation radii. Bar origin is distinguishable, with the tidal disc displaying larger migration overall, higher metallicity difference between the inner and outer disc, as well as a population of inner disc stars displaced to large radii and below the disc-plane. The affect of closest approach on populations of stars formed before, after and during this period is evident. However, bar morphology is also found to be a significant factor in the evolution of disc stellar properties, with similar bars producing similar traits in migration tendency with radius, particularly in vertical stellar motion and in the evolution of central metallicity features. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 18 pages, 12 figures, submitted to the Monthly Notices of the Royal Astronomical Society

Effective reaction temperatures of irreversible dust chemical reactions in a protoplanetary disk

Authors: Lily Ishizaki, Shogo Tachibana, Tamami Okamoto, Daiki Yamamoto, Shigeru Ida

Abstract: Dust particles in protoplanetary disks experience various chemical reactions under different physicochemical conditions through their accretion and diffusion, which results in the radial chemical gradient of dust. We performed three-dimensional Monte Carlo simulations to evaluate the dust trajectories and the progress of fictitious irreversible reactions, of which kinetics is expressed by the John… ▽ More Dust particles in protoplanetary disks experience various chemical reactions under different physicochemical conditions through their accretion and diffusion, which results in the radial chemical gradient of dust. We performed three-dimensional Monte Carlo simulations to evaluate the dust trajectories and the progress of fictitious irreversible reactions, of which kinetics is expressed by the Johnson-Mehl-Avrami equation. The distribution of the highest temperature that each particle experiences before the degree of reaction exceeds a certain level shows the log-normal distribution, and its mode temperature was used as the effective reaction temperature. Semi-analytical prediction formulas of the effective reaction temperature and its dispersion were derived by comparing a reaction timescale with a diffusive transport timescale of dust as a function of the reaction parameters and the disk parameters. The formulas reproduce the numerical results of the effective reaction temperatures and their dispersions within 5.5 and 24 %, respectively, in a wide temperature range (200-1400 K). We applied the formulas for the crystallization of amorphous silicate dust and its oxygen isotope exchange with the H2O vapor based on the experimentally determined kinetics. For sub-micron sized amorphous forsterite dust, the predicted effective reaction temperature for the oxygen isotope exchange was lower than that of crystallization without overlap even considering their dispersions. This suggests that the amorphous silicate dust in the protosolar disk could exchange their oxygen isotopes efficiently with the 16O-poor H2O vapor, resulting in the distinct oxygen isotope compositions from the Sun. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 28 pages, 7 figures, accepted for publication in ApJ

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

Evidence for the Operation of the Hanle and Magneto-Optical Effects in the Scattering Polarization Signals Observed by CLASP2 Across the Mg II h and k Lines

Authors: Ryohko Ishikawa, Javier Trujillo Bueno, Ernest Alsina Ballester, Luca Belluzzi, Tanausú del Pino Alemán, David E. McKenzie, Frédéric Auchère, Ken Kobayashi, Takenori J. Okamoto, Laurel A. Rachmeler, Donguk Song

Abstract: Radiative transfer investigations of the solar Mg II h and k resonance lines around 280~nm showed that, while their circular polarization (Stokes V) signals arise from the Zeeman effect, the linear polarization profiles (Stokes Q and U) are dominated by the scattering of anisotropic radiation and the Hanle and magneto-optical (MO) effects. Using the unprecedented observations of the Mg II and Mn I… ▽ More Radiative transfer investigations of the solar Mg II h and k resonance lines around 280~nm showed that, while their circular polarization (Stokes V) signals arise from the Zeeman effect, the linear polarization profiles (Stokes Q and U) are dominated by the scattering of anisotropic radiation and the Hanle and magneto-optical (MO) effects. Using the unprecedented observations of the Mg II and Mn I resonance lines obtained by the Chromospheric LAyer Spectro-Polarimeter (CLASP2), here we investigate how the linear polarization signals at different wavelengths (i.e., at the center, and at the near and far wings of the k line) vary with the longitudinal component of the magnetic field ($B_{L}$) at their approximate height of formation. The $B_{L}$ is estimated from the V signals in the aforementioned spectral lines. Particular attention is given to the following quantities that are expected to be influenced by the presence of magnetic fields through the Hanle and MO effects: the sign of the U signals, the total linear polarization amplitude ($LP$) and its direction ($χ$) with respect to a reference direction. We find that at the center and near wings of the $k$ line, the behavior of these quantities is significantly different in the observed quiet and plage regions, and that both $LP$ and $χ$ seem to depend on $B_{L}$. These observational results are indicative of the operation of the Hanle effect △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 26 pages, 18 figures, accepted for publication in the Astrophysical Journal

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

Tomography of a solar plage with the Tenerife Inversion Code

Authors: Hao Li, T. del Pino Alemán, J. Trujillo Bueno, R. Ishikawa, E. Alsina Ballester, David E. McKenzie, Frédéric Auchère, Ken Kobayashi, Takenori J. Okamoto, Laurel A. Rachmeler, Donguk Song

Abstract: We apply the Tenerife Inversion Code (TIC) to the plage spectropolarimetric observations obtained by the Chromospheric LAyer SpectroPolarimeter (CLASP2). These unprecedented data consist of full Stokes profiles in the spectral region around the Mg II h and k lines for a single slit position, with around two thirds of the 200 arcsec slit crossing a plage region and the rest crossing an enhanced net… ▽ More We apply the Tenerife Inversion Code (TIC) to the plage spectropolarimetric observations obtained by the Chromospheric LAyer SpectroPolarimeter (CLASP2). These unprecedented data consist of full Stokes profiles in the spectral region around the Mg II h and k lines for a single slit position, with around two thirds of the 200 arcsec slit crossing a plage region and the rest crossing an enhanced network. A former analysis of these data had allowed us to infer the longitudinal component of the magnetic field by applying the weak field approximation (WFA) to the circular polarization profiles, and to assign the inferred magnetic fields to different layers of the solar atmosphere based on the results of previous theoretical radiative transfer investigations. In this work, we apply the recently developed TIC to the same data. We obtain the stratified model atmosphere that fits the intensity and circular polarization profiles at each position along the spectrograph slit and we compare our results for the longitudinal component of the magnetic field with the previously obtained WFA results, highlighting the generally good agreement in spite of the fact that the WFA is known to produce an underestimation when applied to the outer lobes of the Mg II h and k circular polarization profiles. Finally, we use the inverted model atmospheres to give a rough estimation of the energy that could be carried by Alfvèn waves propagating along the chromosphere in the plage and network regions, showing that it is sufficient to compensate the estimated energy losses in the chromosphere of solar active regions. △ Less

Submitted 31 January, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: Accepted for publication in The Astrophysical Journal

CO($J$=1-0) mapping survey of 64 galaxies in the Fornax cluster with the ALMA Morita array

Authors: Kana Morokuma-Matsui, Kenji Bekki, Jing Wang, Paolo Serra, Yusei Koyama, Tomoki Morokuma, Fumi Egusa, Bi-Qing For, Kouichiro Nakanishi, Bäbel S. Koribalski, Takashi Okamoto, Tadayuki Kodama, Bumhyun Lee, Filippo M. Maccagni, Rie E. Miura, Daniel Espada, Tsutomu T. Takeuchi, Dong Yang, Minju M. Lee, Masaki Ueda, Kyoko Matsushita

Abstract: We conduct a $^{12}$C$^{16}$O($J$=1-0) (hereafter CO) mapping survey of 64 galaxies in the Fornax cluster using the ALMA Morita array in cycle 5. CO emission is detected from 23 out of the 64 galaxies. Our sample includes dwarf, spiral and elliptical galaxies with stellar masses of $M_{\rm star}\sim10^{6.3-11.6}$~M$_\odot$. The achieved beam size and sensitivity are $15''\times8''$ and $\sim12$~mJ… ▽ More We conduct a $^{12}$C$^{16}$O($J$=1-0) (hereafter CO) mapping survey of 64 galaxies in the Fornax cluster using the ALMA Morita array in cycle 5. CO emission is detected from 23 out of the 64 galaxies. Our sample includes dwarf, spiral and elliptical galaxies with stellar masses of $M_{\rm star}\sim10^{6.3-11.6}$~M$_\odot$. The achieved beam size and sensitivity are $15''\times8''$ and $\sim12$~mJy~beam$^{-1}$ at the velocity resolution of $\sim10$~km~s$^{-1}$, respectively. We study the cold-gas (molecular- and atomic-gas) properties of 38 subsamples with $M_{\rm star}>10^9$~M$_\odot$ combined with literature HI data. We find that: (1) the low star-formation (SF) activity in the Fornax galaxies is caused by the decrease in the cold-gas mass fraction with respect to stellar mass (hereafter, gas fraction) rather than the decrease of the SF efficiency from the cold gas; (2) the atomic-gas fraction is more heavily reduced than the molecular-gas fraction of such galaxies with low SF activity. A comparison between the cold-gas properties of the Fornax galaxies and their environmental properties suggests that the atomic gas is stripped tidally and by the ram pressure, which leads to the molecular gas depletion with an aid of the strangulation and consequently SF quenching. Pre-processes in the group environment would also play a role in reducing cold-gas reservoirs in some Fornax galaxies. △ Less

Submitted 16 October, 2022; originally announced October 2022.

Comments: 53 pages, 41 figures, accepted for publication in ApJS

The Missing Satellite Problem Outside of the Local Group. II. Statistical Properties of Satellites of Milky Way-like Galaxies

Authors: Masashi Nashimoto, Masayuki Tanaka, Masashi Chiba, Kohei Hayashi, Yutaka Komiyama, Takashi Okamoto

Abstract: We present a new observation of satellite galaxies around seven Milky Way (MW)-like galaxies located outside of the Local Group (LG) using Subaru/Hyper Suprime-Cam imaging data to statistically address the missing satellite problem. We select satellite galaxy candidates using magnitude, surface brightness, Sérsic index, axial ratio, full width half maximum, and surface brightness fluctuation cuts,… ▽ More We present a new observation of satellite galaxies around seven Milky Way (MW)-like galaxies located outside of the Local Group (LG) using Subaru/Hyper Suprime-Cam imaging data to statistically address the missing satellite problem. We select satellite galaxy candidates using magnitude, surface brightness, Sérsic index, axial ratio, full width half maximum, and surface brightness fluctuation cuts, followed by visual screening of false-positives such as optical ghosts of bright stars. We identify 51 secure dwarf satellite galaxies within the virial radius of nine host galaxies, two of which are drawn from the pilot observation presented in Paper I. We find that the average luminosity function of the satellite galaxies is consistent with that of the MW satellites, although the luminosity function of each host galaxy varies significantly. We observe an indication that more massive hosts tend to have a larger number of satellites. Physical properties of the satellites such as the size-luminosity relation is also consistent with the MW satellites. However, the spatial distribution is different; we find that the satellite galaxies outside of LG shows no sign of concentration or alignment, while that of the MW satellites is more concentrated around the host and exhibits a significant alignment. As we focus on relatively massive satellites with $M_V<-10$, we do not expect that the observational incompleteness can be responsible here. This trend might represent a peculiarity of the MW satellites, and further work is needed to understand its origin. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 20 pages, 11 figures, accepted for publication in The Astrophysical Journal

Classification of cosmic structures for galaxies with deep learning: connecting cosmological simulations with observations

Authors: Shigeki Inoue, Xiaotian Si, Takashi Okamoto, Moka Nishigaki

Abstract: We explore the capability of deep learning to classify cosmic structures. In cosmological simulations, cosmic volumes are segmented into voids, sheets, filaments and knots, according to the distribution and kinematics of dark matter (DM), and galaxies are also classified according to the segmentation. However, observational studies cannot adopt this classification method using DM. In this study, w… ▽ More We explore the capability of deep learning to classify cosmic structures. In cosmological simulations, cosmic volumes are segmented into voids, sheets, filaments and knots, according to the distribution and kinematics of dark matter (DM), and galaxies are also classified according to the segmentation. However, observational studies cannot adopt this classification method using DM. In this study, we demonstrate that deep learning can bridge the gap between simulations and observations. Our models are based on three-dimensional convolutional neural networks and trained with data of the distribution of galaxies in a simulation to deduce the structure classes from the galaxies rather than DM. Our model can predict the class labels as accurate as a previous study using DM distribution for the training and prediction. This means that galaxy distribution can be a substitution for DM for the cosmic-structure classification, and our models using galaxies can be directly applied to wide-field survey observations. When observational restrictions are ignored, our model can classify simulated galaxies into the four classes with an accuracy (macro-averaged $F_{\rm 1}$-score) of 64 per cent. If restrictions such as limiting magnitude are considered, our model can classify SDSS galaxies at $\sim100~{\rm Mpc}$ with an accuracy of 60 per cent. In the binary classification distinguishing void galaxies from the others, our model can achieve an accuracy of 88 per cent. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: 15 pages (+Appendices), 21 figures, accepted for publication in MNRAS

Quiet Sun Center to Limb Variation of the Linear Polarization Observed by CLASP2 Across the Mg II h & k Lines

Authors: L. A. Rachmeler, J. Trujillo Bueno, D. E. McKenzie, R. Ishikawa, F. Auchere, K. Kobayashi, R. Kano, T. J. Okamoto, C. W. Bethge, D. Song, E. Alsina Ballester, L. Belluzzi, T. del Pino Aleman, M. Yoshida, T. Shimizu, A. Winebarger, A. R. Kobelski, G. D. Vigil, B. De Pontieu, Y. Suematsu, N. Narukage, M. Kubo, T. Sakao, H. Hara, J. Stepan , et al. (3 additional authors not shown)

Abstract: The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region.… ▽ More The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region. The focus of the work presented here is the center-to-limb variation of the linear polarization across these resonance lines, which is produced by the scattering of anisotropic radiation in the solar atmosphere. The linear polarization signals of the Mg II h & k lines are sensitive to the magnetic field from the low to the upper chromosphere through the Hanle and magneto-optical effects. We compare the observations to theoretical predictions from radiative transfer calculations in unmagnetized semi-empirical models, arguing that magnetic fields and horizontal inhomogeneities are needed to explain the observed polarization signals and spatial variations. This comparison is an important step in both validating and refining our understanding of the physical origin of these polarization signatures, and also in paving the way toward future space telescopes for probing the magnetic fields of the solar upper atmosphere via ultraviolet spectropolarimetry. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 14 pages, 5 figures, to be published in the Astrophysical Journal (ApJ), initially submitted May 2022, revised submission July 2022

Journal ref: ApJ 936 67 (2022)

Origin of highly $r$-process-enhanced stars in a cosmological zoom-in simulation of a Milky Way-like galaxy

Authors: Yutaka Hirai, Timothy C. Beers, Masashi Chiba, Wako Aoki, Derek Shank, Takayuki R. Saitoh, Takashi Okamoto, Junichiro Makino

Abstract: The $r$-process-enhanced (RPE) stars provide fossil records of the assembly history of the Milky Way and the nucleosynthesis of the heaviest elements. Observations by the $R$-Process Alliance (RPA) and others have confirmed that many RPE stars are associated with chemo-dynamically tagged groups, which likely came from accreted dwarf galaxies of the Milky Way (MW). However, we do not know how RPE s… ▽ More The $r$-process-enhanced (RPE) stars provide fossil records of the assembly history of the Milky Way and the nucleosynthesis of the heaviest elements. Observations by the $R$-Process Alliance (RPA) and others have confirmed that many RPE stars are associated with chemo-dynamically tagged groups, which likely came from accreted dwarf galaxies of the Milky Way (MW). However, we do not know how RPE stars are formed. Here, we present the result of a cosmological zoom-in simulation of an MW-like galaxy with $r$-process enrichment, performed with the highest resolution in both time and mass. Thanks to this advancement, unlike previous simulations, we find that most highly RPE ($r$-II; [Eu/Fe] $> +0.7$) stars are formed in low-mass dwarf galaxies that have been enriched in $r$-process elements for [Fe/H] $\,<-2.5$, while those with higher metallicity are formed in situ, in locally enhanced gas clumps that were not necessarily members of dwarf galaxies. This result suggests that low-mass accreted dwarf galaxies are the main formation site of $r$-II stars with [Fe/H] $\,<-2.5$. We also find that most low-metallicity $r$-II stars exhibit halo-like kinematics. Some $r$-II stars formed in the same halo show low dispersions in [Fe/H] and somewhat larger dispersions of [Eu/Fe], similar to the observations. The fraction of simulated $r$-II stars is commensurate with observations from the RPA, and the distribution of the predicted [Eu/Fe] for halo $r$-II stars matches that observed. These results demonstrate that RPE stars can be valuable probes of the accretion of dwarf galaxies in the early stages of their formation. △ Less

Submitted 13 December, 2022; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: 20 pages, 15 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 517, Issue 4, December 2022, Pages 4856-4874

Modelling self-interacting dark matter substructures I: Calibration with N-body simulations of a Milky-Way-sized halo and its satellite

Authors: Masato Shirasaki, Takashi Okamoto, Shin'ichiro Ando

Abstract: We study evolution of single subhaloes with their masses of $\sim10^9 M_\odot$ in a Milky-Way-sized host halo for self-interacting dark matter (SIDM) models. We perform dark-matter-only N-body simulations of dynamical evolution of individual subhaloes orbiting its host by varying self-scattering cross sections (including a velocity-dependent scenario), subhalo orbits, and internal properties of th… ▽ More We study evolution of single subhaloes with their masses of $\sim10^9 M_\odot$ in a Milky-Way-sized host halo for self-interacting dark matter (SIDM) models. We perform dark-matter-only N-body simulations of dynamical evolution of individual subhaloes orbiting its host by varying self-scattering cross sections (including a velocity-dependent scenario), subhalo orbits, and internal properties of the subhalo. We calibrate a gravothermal fluid model to predict time evolution in spherical mass density profiles of isolated SIDM haloes with the simulations. We find that tidal effects of SIDM subhaloes can be described with a framework developed for the case of collision-less cold dark matter (CDM), but a shorter typical time scale for the mass loss due to tidal stripping is required to explain our SIDM simulation results. As long as the cross section is less than $\sim10\, \mathrm{cm}^2/\mathrm{g}$ and initial states of subhaloes are set within a $2σ$-level scatter at redshifts of $\sim2$ predicted by the standard $Λ$CDM cosmology, our simulations do not exhibit a prominent feature of gravothermal collapse in the subhalo central density for 10 Gyr. We develop a semi-analytic model of SIDM subhaloes in a time-evolving density core of the host with tidal stripping and self-scattering ram pressure effects. Our semi-analytic approach provides a simple, efficient and physically-intuitive prediction of SIDM subhaloes, but further improvements are needed to account for baryonic effects in the host and the gravothermal instability accelerated by tidal stripping effects. △ Less

Submitted 5 September, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: 19 pages, 12 figures, 1 table. Accepted for publication in MNRAS. Note that results in the previous manuscript (v1) were affected by a bug in our SIDM implementation. We fixed the bug and updated all results accordingly

Massive core/star formation triggered by cloud-cloud collision: II High-speed collisions of magnetized clouds

Authors: Nirmit Sakre, Asao Habe, Alex R. Pettitt, Takashi Okamoto, Rei Enokiya, Yasuo Fukui, Takashi Hosokawa

Abstract: We study the effects of the magnetic fields on the formation of massive, self-gravitationally bound cores (MBCs) in high-speed cloud-cloud collisions (CCCs). Extending our previous work (Sakre et al. 2021), we perform magnetohydrodynamic simulations following the high-speed (20 - 40 km s$^{-1}$) collisions between two magnetized (4 $μ$G initially), turbulent clouds of different sizes in the range… ▽ More We study the effects of the magnetic fields on the formation of massive, self-gravitationally bound cores (MBCs) in high-speed cloud-cloud collisions (CCCs). Extending our previous work (Sakre et al. 2021), we perform magnetohydrodynamic simulations following the high-speed (20 - 40 km s$^{-1}$) collisions between two magnetized (4 $μ$G initially), turbulent clouds of different sizes in the range of 7 - 20 pc. We show that a magnetic field effect hinders the core growth, particularly after a short-duration collision during which cores cannot get highly bound. In such a case, a shocked region created by the collision rapidly expands to the ambient medium owing to the enhanced magnetic pressure, resulting in the destruction of the highly unbound cores and suppression of gas accretion to massive cores. This negative effect on the MBC formation is a phenomenon not seen in the past hydrodynamic simulations of similar CCC models. Together with our previous work, we conclude that the magnetic fields provide the two competing effects on the MBC formation in CCC; while they promote the mass accumulation into cores during a collision, they operate to destroy cores or hinder the core growth after the collision. The duration of collision determines which effect prevails, providing the maximum collision speed for the MBC formation with given colliding clouds. Our results agree with the observed trend among CCC samples in the corresponding column density range; clouds with higher relative velocity require higher column density for the formation of massive stars (Enokiya et al. 2021). △ Less

Submitted 14 May, 2022; originally announced May 2022.

Comments: Submitted to MNRAS journal

Electrostatic Dust Ejection From Asteroid (3200) Phaethon With the Aid of Mobile Alkali Ions at Perihelion

Authors: Hiroshi Kimura, Katsuhito Ohtsuka, Shota Kikuchi, Keiji Ohtsuki, Tomoko Arai, Fumi Yoshida, Naoyuki Hirata, Hiroki Senshu, Koji Wada, Takayuki Hirai, Peng K. Hong, Masanori Kobayashi, Ko Ishibashi, Manabu Yamada, Takaya Okamoto

Abstract: The asteroid (3200) Phaethon is known to be the parent body of the Geminids, although meteor showers are commonly associated with the activity of periodic comets. What is most peculiar to the asteroid is its comet-like activity in the ejection of micrometer-sized dust particles at every perihelion passage, while the activity of the asteroid has never been identified outside the near-perihelion zon… ▽ More The asteroid (3200) Phaethon is known to be the parent body of the Geminids, although meteor showers are commonly associated with the activity of periodic comets. What is most peculiar to the asteroid is its comet-like activity in the ejection of micrometer-sized dust particles at every perihelion passage, while the activity of the asteroid has never been identified outside the near-perihelion zone at $0.14~\mathrm{au}$ from the Sun. From the theoretical point of view, we argue that the activity of the asteroid is well explained by the electrostatic lofting of micrometer-sized dust particles with the aid of mobile alkali ions at high temperatures. The mass-loss rates of micrometer-sized particles from the asteroid in our model is entirely consistent with the values inferred from visible observations of Phaethon's dust tail. For millimeter-sized particles, we predict three orders of magnitudes higher mass-loss rates, which could also account for the total mass of the Geminid meteoroid stream by the electrostatic lofting mechanism. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: 7 pages. 2 figures, to appear in Icarus Letters

Monte Carlo Simulation of Dust Particles in a Protoplanetary Disk: Crystalline to Amorphous Silicate Ratio in Comets

Authors: Tamami Okamoto, Shigeru Ida

Abstract: Observationally inferred crystalline abundance in silicates in comets, which should have been formed in the outer region of a protoplanetary disk, is relatively high (~ 10-60%), although crystalline silicates would be formed by annealing of amorphous precursors in the disk inner region. In order to quantitatively address this puzzle, we have performed Monte Carlo simulation of advection/diffusion… ▽ More Observationally inferred crystalline abundance in silicates in comets, which should have been formed in the outer region of a protoplanetary disk, is relatively high (~ 10-60%), although crystalline silicates would be formed by annealing of amorphous precursors in the disk inner region. In order to quantitatively address this puzzle, we have performed Monte Carlo simulation of advection/diffusion of silicate particles in a turbulent disk, in the setting based on pebble accretion model: pebbles consisting of many small amorphous silicates embedded in icy mantle are formed in the disk outer region, silicate particles are released at the snow line, crystalline silicate particles are produced at the annealing line, the silicate particles diffused beyond the snow line, and they eventually stick to drifting pebbles to come back to the snow line. In a simple case without the sticking and with a steady pebble flux, we show through the simulations and analytical arguments that crystalline components in silicate materials beyond the snow line is robustly and uniformly ~ 5%. On the other hand, in a more realistic case with the sticking and with a decaying pebble flux, the crystalline abundance is raised up to ~ 20-25%, depending on the ratio of decay and diffusion timescales. This abundance is consistent with the observations. In this investigation, we assume a simple steady accretion disk. The simulations coupled with the disk evolution is needed for more detailed comparison with observed data. △ Less

Submitted 20 January, 2022; v1 submitted 14 January, 2022; originally announced January 2022.

Comments: Accepted for publication in the Astrophysical Journal

Differences in star formation activity between tidally triggered and isolated bars: a case study of NGC 4303 and NGC 3627

Authors: Elizabeth J. Iles, Alex R. Pettitt, Takashi Okamoto

Abstract: Galactic bars are important drivers of galactic evolution, and yet how they impact the interstellar medium and correspondingly star formation, remains unclear. We present simulation results for two barred galaxies with different formation mechanisms, bars formed in isolation or via a tidal interaction, to consider the spatially and temporally varying trends of star formation. We focus on the early… ▽ More Galactic bars are important drivers of galactic evolution, and yet how they impact the interstellar medium and correspondingly star formation, remains unclear. We present simulation results for two barred galaxies with different formation mechanisms, bars formed in isolation or via a tidal interaction, to consider the spatially and temporally varying trends of star formation. We focus on the early (< 1Gyr) epoch of bar formation so that the interaction is clearly identifiable. The nearby NGC 4303 (isolated) and NGC 3627 (interaction history) are selected as observational analogues to tailor these simulations. Regardless of formation mechanism, both models show similar internal dynamical features, although the interaction appears to promote bar-arm disconnection in the outer disc velocity structure. Both bars trigger similar boosts in star formation (79%; 66%), while the interaction also triggers an earlier 31% burst. Significant morphological dependence is observed in the relation between surface gas and star formation rate. In both cases, the bar component is notably steepest; the arm is similar to the overall disc average; and the inter-arm clearly the shallowest. A distinguishable feature of the tidal disc is the presence of moderately dense, inefficiently star forming gas mostly confined to tidal debris outside the optical disc. The tidal disc also exhibits a unique trend of radially increasing star formation efficiency and a clear dearth of star formation which persists along the bar between the centre and bar ends. These are potential signatures for identifying a barred system post-interaction. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: 19 pages, 15 figures, accepted for publication in MNRAS

Testing the effect of resolution on gravitational fragmentation with Lagrangian hydrodynamic schemes

Authors: Yasuyoshi Yamamoto, Takashi Okamoto, Takayuki R. Saitoh

Abstract: To study the resolution required for simulating gravitational fragmentation with newly developed Lagrangian hydrodynamic schemes, Meshless Finite Volume method (MFV) and Meshless Finite Mass method (MFM), we have performed a number of simulations of the Jeans test and compared the results with both the expected analytic solution and results from the more standard Lagrangian approach: Smoothed Part… ▽ More To study the resolution required for simulating gravitational fragmentation with newly developed Lagrangian hydrodynamic schemes, Meshless Finite Volume method (MFV) and Meshless Finite Mass method (MFM), we have performed a number of simulations of the Jeans test and compared the results with both the expected analytic solution and results from the more standard Lagrangian approach: Smoothed Particle Hydrodynamics (SPH). We find that the different schemes converge to the analytic solution when the diameter of a fluid element is smaller than a quarter of the Jeans wavelength, $λ_\mathrm{J}$. Among the three schemes, SPH/MFV shows the fastest/slowest convergence to the analytic solution. Unlike the well-known behaviour of Eulerian schemes, none of the Lagrangian schemes investigated displays artificial fragmentation when the perturbation wavelength, $λ$, is shorter than $λ_\mathrm{J}$, even at low numerical resolution. For larger wavelengths ($λ> λ_\mathrm{J}$) the growth of the perturbation is delayed when it is not well resolved. Furthermore, with poor resolution, the fragmentation seen with the MFV scheme proceeds very differently compared to the converged solution. All these results suggest that, when unresolved, the ratio of the magnitude of hydrodynamic force to that of self-gravity at the sub-resolution scale is the largest/smallest in MFV/SPH, the reasons for which we discussed in detail. These tests are repeated to investigate the effect of kernels of higher-order than the fiducial cubic spline. Our results indicate that the standard deviation of the kernel is a more appropriate definition of the size of a fluid element than its compact support radius. △ Less

Submitted 20 April, 2021; v1 submitted 18 April, 2021; originally announced April 2021.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS. The typos in the authors' names have been corrected

A phase-space view of cold-gas properties of Virgo-cluster galaxies: multiple quenching processes at work?

Authors: Kana Morokuma-Matsui, Tadayuki Kodama, Tomoki Morokuma, Kouichiro Nakanishi, Yusei Koyama, Takuji Yamashita, Shuhei Koyama, Takashi Okamoto

Abstract: We investigate the cold-gas properties of massive Virgo galaxies ($>10^9$ M$_\odot$) at $<3R_{200}$ ($R_{200}$ is the radius where the mean interior density is 200 times the critical density) on the projected phase-space diagram (PSD) with the largest archival dataset to date to understand the environmental effect on galaxy evolution in the Virgo cluster. We find: lower HI and H$_2$ mass fractions… ▽ More We investigate the cold-gas properties of massive Virgo galaxies ($>10^9$ M$_\odot$) at $<3R_{200}$ ($R_{200}$ is the radius where the mean interior density is 200 times the critical density) on the projected phase-space diagram (PSD) with the largest archival dataset to date to understand the environmental effect on galaxy evolution in the Virgo cluster. We find: lower HI and H$_2$ mass fractions and higher star-formation efficiencies (SFEs) from HI and H$_2$ in the Virgo galaxies than the field galaxies for matched stellar masses; the Virgo galaxies generally follow the field relationships between the offset from the main sequence of the star-forming galaxies [$Δ$(MS)] with gas fractions and SFEs but slightly offset to lower gas fractions or higher SFEs than field galaxies at $Δ({\rm MS})< 0$; lower gas fractions in galaxies with smaller clustocentric distance and velocity; lower gas fractions in the galaxies in the W cloud, a substructure of the Virgo cluster. Our results suggest the cold-gas properties of some Virgo galaxies are affected by their environment at least at $3 R_{200}$ maybe via strangulation and/or pre-processes and HI and H$_2$ in some galaxies are removed by ram pressure at $<1.5 R_{200}$. Our data cannot rule the possibility of the other processes such as strangulation and galaxy harassment accounting for the gas reduction in some galaxies at $<1.5 R_{200}$. Future dedicated observations of a mass-limited complete sample are required for definitive conclusions. △ Less

Submitted 24 March, 2021; v1 submitted 9 March, 2021; originally announced March 2021.

Comments: 32 pages, 16 figures, accepted for publication in ApJ; Updated stellar masses of galaxies after the first acceptance, relevant texts, and figures

Mapping Solar Magnetic Fields from the Photosphere to the Base of the Corona

Authors: Ryohko Ishikawa, Javier Trujillo Bueno, Tanausu del Pino Aleman, Takenori J. Okamoto, David E. McKenzie, Frederic Auchere, Ryouhei Kano, Donguk Song, Masaki Yoshida, Laurel A. Rachmeler, Ken Kobayashi, Hirohisa Hara, Masahito Kubo, Noriyuki Narukage, Taro Sakao, Toshifumi Shimizu, Yoshinori Suematsu, Christian Bethge, Bart De Pontieu, Alberto Sainz Dalda, Genevieve D. Vigil, Amy Winebarger, Ernest Alsina Ballester, Luca Belluzzi, Jiri Stepan , et al. (3 additional authors not shown)

Abstract: Routine ultraviolet imaging of the Sun's upper atmosphere shows the spectacular manifestation of solar activity; yet we remain blind to its main driver, the magnetic field. Here we report unprecedented spectropolarimetric observations of an active region plage and its surrounding enhanced network, showing circular polarization in ultraviolet (Mg II $h$ & $k$ and Mn I) and visible (Fe I) lines. We… ▽ More Routine ultraviolet imaging of the Sun's upper atmosphere shows the spectacular manifestation of solar activity; yet we remain blind to its main driver, the magnetic field. Here we report unprecedented spectropolarimetric observations of an active region plage and its surrounding enhanced network, showing circular polarization in ultraviolet (Mg II $h$ & $k$ and Mn I) and visible (Fe I) lines. We infer the longitudinal magnetic field from the photosphere to the very upper chromosphere. At the top of the plage chromosphere the field strengths reach more than 300 gauss, strongly correlated with the Mg II $k$ line core intensity and the electron pressure. This unique mapping shows how the magnetic field couples the different atmospheric layers and reveals the magnetic origin of the heating in the plage chromosphere. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: 50 pages, 11 figures, 1 table, published in Science Advances

Journal ref: Science Advances, Vol. 7, no. 8, eabe8406, 2021

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)

Revisiting Soltan's argument based on a semi-analytical model for galaxy and black hole evolution

Authors: Hikari Shirakata, Toshihiro Kawaguchi, Takashi Okamoto, Masahiro Nagashima, Taira Oogi

Abstract: We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired dur… ▽ More We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired during the super-Eddington accretion phases to the total SMBH mass becomes larger for less massive black holes and at higher redshift. Even at z = 0, SMBHs with > 1e+9 Msun have acquired more than 50% of their mass by super-Eddington accretions, which is apparently inconsistent with classical Soltan's argument. However, the mass-weighted radiation efficiency of SMBHs with > 1e+8 Msun obtained with our model, is about 0.08 at z = 0, which is consistent with Soltan's argument within the observational uncertainties. We, therefore, conclude that Soltan's argument cannot reject the possibility that SMBHs are grown mainly by super-Eddington accretions. △ Less

Submitted 3 June, 2020; originally announced June 2020.

Comments: 10 pages, 6 figures. Accepted for publication in ApJ

Impact Ejecta near the Impact Point Observed using Ultra-high-speed Imaging and SPH Simulations, and a Comparison of the Two Methods

Authors: Takaya Okamoto, Kosuke Kurosawa, Hidenori Genda, Takafumi Matsui

Abstract: High-speed impact ejecta at velocities comparable to the impact velocity are expected to contribute to material transport between planetary bodies and deposition of ejecta far from the impact crater. We investigated the behavior of high-speed ejecta produced at angles of 45 and 90 degrees, using both experimental and numerical methods. The experimental system developed at the Planetary Exploration… ▽ More High-speed impact ejecta at velocities comparable to the impact velocity are expected to contribute to material transport between planetary bodies and deposition of ejecta far from the impact crater. We investigated the behavior of high-speed ejecta produced at angles of 45 and 90 degrees, using both experimental and numerical methods. The experimental system developed at the Planetary Exploration Research Center of Chiba Institute of Technology (Japan) allowed us to observe the initial growth of the ejecta. We succeeded in imaging high-speed ejecta at 0.2 $\mathrm{μ}$s intervals for impacts of polycarbonate projectiles of 4.8 mm diameter onto a polycarbonate plate at an impact velocity of ~4 km s$^{-1}$. Smoothed particle hydrodynamics (SPH) simulations of various numerical resolutions were conducted for the same impact conditions as pertaining to the experiments. We compared the morphology and velocities of the ejecta for the experiments and simulations, and we confirmed a close match for high-resolution simulations (with $\geq10^6$ SPH particles representing the projectile). According to the ejecta velocity distributions obtained from our high-resolution simulations, the ejection velocities of the high-speed ejecta for oblique impacts are much greater than those for vertical impacts. The translational motion of penetrating projectiles parallel to the target surface in oblique impacts could cause long-term, sustained acceleration at the root of the ejecta. △ Less

Submitted 18 March, 2020; originally announced March 2020.

Comments: 33 pages, 16 figures, 1 supporting information, accepted for publication in Journal of Geophysical Research: Planets

Massive core/star formation triggered by cloud-cloud collision: Effect of magnetic field

Authors: Nirmit Sakre, Asao Habe, Alex R. Pettitt, Takashi Okamoto

Abstract: We study effect of magnetic field on massive dense core formation in colliding unequal molecular clouds by performing magnetohydrodynamic simulations with sub-parsec resolution (0.015 pc) that can resolve the molecular cores. Initial clouds with the typical gas density of the molecular clouds are immersed in various uniform magnetic fields. The turbulent magnetic fields in the clouds consistent wi… ▽ More We study effect of magnetic field on massive dense core formation in colliding unequal molecular clouds by performing magnetohydrodynamic simulations with sub-parsec resolution (0.015 pc) that can resolve the molecular cores. Initial clouds with the typical gas density of the molecular clouds are immersed in various uniform magnetic fields. The turbulent magnetic fields in the clouds consistent with the observation by Crutcher et al. (2010) are generated by the internal turbulent gas motion before the collision, if the uniform magnetic field strength is 4.0 $μ$G. The collision speed of 10 km s$^{-1}$ is adopted, which is much larger than the sound speeds and the Alfvén speeds of the clouds. We identify gas clumps with gas densities greater than 5 $\times$ 10$^{-20}$ g cm$^{-3}$ as the dense cores and trace them throughout the simulations to investigate their mass evolution and gravitational boundness. We show that a greater number of massive, gravitationally bound cores are formed in the strong magnetic field (4.0 $μ$G) models than the weak magnetic field (0.1 $μ$G) models. This is partly because the strong magnetic field suppresses the spatial shifts of the shocked layer that should be caused by the nonlinear thin shell instability. The spatial shifts promote formation of low-mass dense cores in the weak magnetic field models. The strong magnetic fields also support low-mass dense cores against gravitational collapse. We show that the numbers of massive, gravitationally bound cores formed in the strong magnetic field models are much larger than the isolated, non-colliding cloud models, which are simulated for comparison. We discuss the implications of our numerical results on massive star formation. △ Less

Submitted 22 July, 2020; v1 submitted 5 March, 2020; originally announced March 2020.

Comments: Published in PASJ (for special issue on cloud-cloud collision). This arXiv version is similar to the published paper ( https://academic.oup.com/pasj/article/doi/10.1093/pasj/psaa059/5867032 )

Journal ref: PASJ psaa059 (2020)

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

Semi-analytic modelling of AGNs: auto-correlation function and halo occupation

Authors: Taira Oogi, Hikari Shirakata, Masahiro Nagashima, Takahiro Nishimichi, Toshihiro Kawaguchi, Takashi Okamoto, Tomoaki Ishiyama, Motohiro Enoki

Abstract: The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving… ▽ More The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving volume. We have introduced a new time-scale of gas accretion on to the supermassive black holes to account for the loss of the angular momentum on small scales, which is required to match the faint end of the observed X-ray LF. The large simulation box allows us accurate determination of the auto-correlation function of the AGNs. The model prediction indicates that this time-scale plays a significant role in allowing massive haloes to host relatively faint population of AGNs, leading to a higher bias factor for those AGNs. The model predictions are in agreement with observations of X-ray selected AGNs in the luminosity range $10^{41.5}~\mathrm{erg} \ \mathrm{s}^{-1} \leq L_{2-10\mathrm{keV}} \leq 10^{44.5}~\mathrm{erg} \ \mathrm{s}^{-1}$, with the typical host halo mass of $10^{12.5-13.5} h^{-1}\,{\rm M}_{\odot}$ at $z \lesssim 1$. This result shows that the observational clustering measurements impose an independent constraint on the accretion time-scale complementary to the LF measurements. Moreover, we find that not only the effective halo mass corresponding to the overall bias factor, but the extended shape of the predicted AGN correlation function shows remarkable agreement with those from observations. Further observational efforts towards the low luminosity end at $z \sim 1$ would give us stronger constraints on the triggering mechanisms of AGN activities through their clustering. △ Less

Submitted 4 July, 2020; v1 submitted 25 December, 2019; originally announced December 2019.

Comments: 20 pages, 13 figures, accepted in MNRAS. Error estimates of the 2PCFs and the fitting to them have been improved. Further discussions and references have been added. Main results remain unchanged

Report number: YITP-19-125

Estimating the temperature and density of a spicule from 100 GHz data obtained with ALMA

Authors: Masumi Shimojo, Tomoko Kawate, Takenori J. Okamoto, Takaaki Yokoyama, Noriyuki Narukage, Taro Sakao, Kazumasa Iwai, Gregory D. Fleishman, Kazunari Shibata

Abstract: We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304Å images (MgII/304A spicule). The other one is a spicule seen in the 100GHz im… ▽ More We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304Å images (MgII/304A spicule). The other one is a spicule seen in the 100GHz images obtained with ALMA (100GHz spicule). Although the 100GHz spicule overlapped with the MgII/304A spicule in the early phase, it did not show any corresponding structures in the IRIS Mg II and AIA 304A images after the early phase. It suggests that the spicules are individual events and do not have a physical relationship. To obtain the physical parameters of the 100GHz spicule, we estimate the optical depths as a function of temperature and density using two different methods. One is using the observed brightness temperature by assuming a filling factor, and the other is using an emission model for the optical depth. As a result of comparing them, the kinetic temperature of the plasma and the number density of ionized hydrogens in the 100GHz spicule are ~6800 K and 2.2 x 10^10 cm^-3. The estimated values can explain the absorbing structure in the 193A image, which appear as a counterpart of the 100GHz spicule. These results suggest that the 100GHz spicule presented in this paper is classified to a macrospicule without a hot sheath in former terminology. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 8 pages, 4 figures, Accepted to ApJL

Journal ref: Astrophysical Journal Letters, Vol. 888, L28 (2020)

IRIS Mg II Observations and Non-LTE Modeling of Off-limb Spicules in a Solar Polar Coronal Hole

Authors: Akiko Tei, Stanislav Gunar, Petr Heinzel, Takenori J. Okamoto, Jiri Stepan, Sonja Jejcic, Kazunari Shibata

Abstract: We investigated the off-limb spicules observed in the Mg II h and k lines by IRIS in a solar polar coronal hole. We analyzed the large dataset of obtained spectra to extract quantitative information about the line intensities, shifts, and widths. The observed Mg II line profiles are broad and double-peaked at lower altitudes, broad but flat-topped at middle altitudes, and narrow and single-peaked… ▽ More We investigated the off-limb spicules observed in the Mg II h and k lines by IRIS in a solar polar coronal hole. We analyzed the large dataset of obtained spectra to extract quantitative information about the line intensities, shifts, and widths. The observed Mg II line profiles are broad and double-peaked at lower altitudes, broad but flat-topped at middle altitudes, and narrow and single-peaked with the largest Doppler shifts at higher altitudes. We use 1D non-LTE vertical slab models (i.e. models which consider departures from Local Thermodynamic Equilibrium) in single-slab and multi-slab configurations to interpret the observations and to investigate how a superposition of spicules along the line of sight (LOS) affects the synthetic Mg II line profiles. The used multi-slab models are either static, i.e. without any LOS velocities, or assume randomly assigned LOS velocities of individual slabs, representing the spicule dynamics. We conducted such single-slab and multi-slab modeling for a broad set of model input parameters and showed the dependence of the Mg II line profiles on these parameters. We demonstrated that the observed line widths of the h and k line profiles are strongly affected by the presence of multiple spicules along the LOS. We later showed that the profiles obtained at higher altitudes can be reproduced by single-slab models representing individual spicules. We found that the multi-slab model with a random distribution of the LOS velocities ranging from -25 to 25 km s$^{-1}$ can well reproduce the width and the shape of Mg II profiles observed at middle altitudes. △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: 19 pages, 11 figures, accepted in the Astrophysical Journal

Shock vaporization/devolatilization of evaporitic minerals, halite and gypsum, in an open system investigated by a two-stage light gas gun

Authors: Kosuke Kurosawa, Ryota Moriwaki, Goro Komatsu, Takaya Okamoto, Hiroshi Sakuma, Hikaru Yabuta, Takafumi Matsui

Abstract: Dry lakebeds might constitute large volatile reservoirs on Mars. Hypervelocity impacts onto ancient dry lakebeds would have affected the volatile distribution on Mars. We developed a new experimental method to investigate the response of evaporitic minerals (halite and gypsum) to impact shocks in an open system. This technique does not result in chemical contamination from the operation of the gas… ▽ More Dry lakebeds might constitute large volatile reservoirs on Mars. Hypervelocity impacts onto ancient dry lakebeds would have affected the volatile distribution on Mars. We developed a new experimental method to investigate the response of evaporitic minerals (halite and gypsum) to impact shocks in an open system. This technique does not result in chemical contamination from the operation of the gas gun. The technique is termed the two-valve method and the gun system is located in the Planetary Exploration Research Center, Chiba Institute of Technology, Japan. We detected the vaporization of halite at 31 GPa and devolatilization from gypsum at 11 GPa, suggesting that impact-induced volatile release from dry lakebeds has periodically occurred throughout Martian history. The vaporization of halite deposits might have enhanced the production of perchlorates, which are found globally on Mars. The water loss from gypsum possibly explains the coexisting types of Ca-sulfates found in Gale Crater. △ Less

Submitted 10 June, 2019; originally announced June 2019.

Comments: 17 pages, 4 figures, 1 supporting information, accepted for publication in Geophysical Research Letters

Slowing down of cosmic growth of supermassive black holes: Theoretical prediction of the Eddington ratio distribution

Authors: Hikari Shirakata, Toshihiro Kawaguchi, Taira Oogi, Takashi Okamoto, Masahiro Nagashima

Abstract: We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the hos… ▽ More We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the host galaxies is higher at higher redshift. We also find that the super- Eddington growth is more common for less massive supermassive black holes, supporting an idea that supermassive black holes have been formed via super-Eddington accretion. These results indicate the "slowing down" of cosmic growth of supermassive black holes: the growth of supermassive black holes with a higher Eddington ratio peaks at higher redshift. We also show the effect of the sample selection on the shape of the Eddington ratio distribution functions and find that shallower observations will miss active galactic nuclei with not only the smaller but also higher Eddington ratios. △ Less

Submitted 7 May, 2019; originally announced May 2019.

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

The Missing Satellite Problem Outside of the Local Group: I -- Pilot Observation

Authors: Masayuki Tanaka, Masashi Chiba, Kohei Hayashi, Yutaka Komiyama, Takashi Okamoto, Andrew Cooper, Sakurako Okamoto, Lee Spitler

Abstract: We present results from a pilot observation of nearby (~20 Mpc) galaxies with mass similar to that of the Milky Way (MW) to address the missing satellite problem. This is the first paper from an on-going project to address the problem with a statistical sample of galaxies outside of the Local Group (LG) without employing an assumption that the LG is a typical halo in the Universe. Thanks to the cl… ▽ More We present results from a pilot observation of nearby (~20 Mpc) galaxies with mass similar to that of the Milky Way (MW) to address the missing satellite problem. This is the first paper from an on-going project to address the problem with a statistical sample of galaxies outside of the Local Group (LG) without employing an assumption that the LG is a typical halo in the Universe. Thanks to the close distances of our targets, dwarf galaxies around them can be identified as extended, diffuse galaxies. By applying a surface brightness cut together with a careful visual screening to remove artifacts and background contamination, we construct a sample of dwarf galaxies. The luminosity function (LF) of one of the targets is broadly consistent with that of the MW, but the other has a more abundant dwarf population. Numerical simulations by Okamoto (2013) seem to overpredict the number of dwarfs on average, while more recent predictions from Copernicus Complexio are in a better agreement. In both observations and simulations, there is a large diversity in the LFs, demonstrating the importance of addressing the missing satellite problem with a statistically representative sample. We also characterize the projected spatial distributions of the satellites and do not observe strong evidence for alignments around the central galaxies. Based on this successful pilot observation, we are carrying out further observations to increase the sample of nearby galaxies, which we plan to report in our future paper. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: 13 pages, 9 figures, accepted for publication in ApJ

ALMA Observations of the Solar Chromosphere on the Polar Limb

Authors: Takaaki Yokoyama, Masumi Shimojo, Takenori J. Okamoto, Haruhisa Iijima

Abstract: We report the results of the Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the solar chromosphere on the southern polar limb. Coordinated observations with the Interface Region Imaging Spectrograph (IRIS) are also conducted. ALMA provided unprecedented high spatial resolution in the millimeter band ($\approx$ 2.0 arcsec) at 100 GHz frequency with a moderate cadence (20 s). T… ▽ More We report the results of the Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the solar chromosphere on the southern polar limb. Coordinated observations with the Interface Region Imaging Spectrograph (IRIS) are also conducted. ALMA provided unprecedented high spatial resolution in the millimeter band ($\approx$ 2.0 arcsec) at 100 GHz frequency with a moderate cadence (20 s). The results are as follows: (1) The ALMA 100 GHz images show saw-tooth patterns on the limb, and a comparison with SDO/AIA 171Å images shows a good correspondence of the limbs with each other. (2) The ALMA 100 GHz movie shows a dynamic thorn-like structure elongating from the saw-tooth patterns on the limb, with lengths reaching at least 8 arcsec, thus suggesting jet-like activity in the ALMA microwave range. These ALMA jets are in good correspondence with IRIS jet clusters. (3) A blob ejection event is observed. By comparing with the IRIS Mg II slit-jaw images, the trajectory of the blob is located along the spicular patterns. △ Less

Submitted 3 July, 2018; originally announced July 2018.

Comments: Accepted for publication in The Astrophysical Journal

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 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

Blue wing enhancement of the chromospheric Mg II h and k lines in a solar flare

Authors: Akiko Tei, Takahito Sakaue, Takenori J. Okamoto, Tomoko Kawate, Petr Heinzel, Satoru Ueno, Ayumi Asai, Kiyoshi Ichimoto, Kazunari Shibata

Abstract: We performed coordinated observations of AR 12205, which produced a C-class flare on 2014 November 11, with the Interface Region Imaging Spectrograph (IRIS) and the Domeless Solar Telescope (DST) at Hida Observatory. Using spectral data in the Si IV 1403 Å, C II 1335 Å, and Mg II h and k lines from IRIS and the Ca II K, Ca II 8542 Å, and H$α$ lines from DST, we investigated a moving flare kernel d… ▽ More We performed coordinated observations of AR 12205, which produced a C-class flare on 2014 November 11, with the Interface Region Imaging Spectrograph (IRIS) and the Domeless Solar Telescope (DST) at Hida Observatory. Using spectral data in the Si IV 1403 Å, C II 1335 Å, and Mg II h and k lines from IRIS and the Ca II K, Ca II 8542 Å, and H$α$ lines from DST, we investigated a moving flare kernel during the flare. In the Mg II h line, the leading edge of the flare kernel showed the intensity enhancement in the blue wing, and the smaller intensity of the blue-side peak (h2v) than that of the red-side one (h2r). The blueshift lasted for 9-48 s with a typical speed of 10.1 $\pm$ 2.6 km s$^{-1}$ and it was followed by the high intensity and the large redshift with a speed of up to 51 km s$^{-1}$ detected in the Mg II h line. The large redshift was a common property for all six lines but the blueshift prior to it was found only in the Mg II lines. A cloud modeling of the Mg II h line suggests that the blue wing enhancement with such peak difference can be caused by a chromospheric-temperature (cool) upflow. We discuss a scenario in which an upflow of cool plasma is lifted up by expanding hot plasma owing to the deep penetration of non-thermal electrons into the chromosphere. Furthermore, we found that the blueshift persisted without any subsequent redshift in the leading edge of the flare kernel during its decaying phase. The cause of such long-lasting blueshift is also discussed. △ Less

Submitted 15 March, 2018; v1 submitted 14 March, 2018; originally announced March 2018.

Comments: 19 pages, 13 figures, accepted for publication in PASJ, replaced for the correction of the name of one author

The New Numerical Galaxy Catalogue (ν^2 GC): Properties of Active Galactic Nuclei and Their Host Galaxies

Authors: Hikari Shirakata, Takashi Okamoto, Toshihiro Kawaguchi, Masahiro Nagashima, Tomoaki Ishiyama, Ryu Makiya, Masakazu A. R. Kobayashi, Motohiro Enoki, Taira Oogi, Katsuya Okoshi

Abstract: We present the latest results of a semi-analytic model of galaxy formation, "New Numerical Galaxy Catalogue", which is combined with large cosmological N-body simulations. This model can reproduce statistical properties of galaxies at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and supermassive black holes, especially on the accretion timescale onto black holes. We find th… ▽ More We present the latest results of a semi-analytic model of galaxy formation, "New Numerical Galaxy Catalogue", which is combined with large cosmological N-body simulations. This model can reproduce statistical properties of galaxies at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and supermassive black holes, especially on the accretion timescale onto black holes. We find that the number density of AGNs at z < 1.5 and at hard X-ray luminosity 10^{ 44 }< erg/s is underestimated compared with recent observational estimates when we assume the exponentially decreasing accretion rate and the accretion timescale which is proportional to the dynamical time of the host halo or the bulge, as is often assumed in semi-analytic models. We show that to solve this discrepancy, the accretion timescale of such less luminous AGNs instead should be a function of the black hole mass and the accreted gas mass. This timescale can be obtained from a phenomenological modelling of the gas angular momentum loss in the circumnuclear torus and/or the accretion disc. Such models predict a longer accretion timescale for less luminous AGNs at z < 1.0 than bright QSOs whose accretion timescale would be 10^{ 7-8 } yr. With this newly introduced accretion timescale, our model can explain the observed luminosity functions of AGNs at z < 6.0. △ Less

Submitted 31 October, 2018; v1 submitted 6 February, 2018; originally announced February 2018.

Comments: 31 pages and 29 figures. Accepted for publication in MNRAS

Super-strong Magnetic Field in Sunspots

Authors: Takenori J. Okamoto, Takashi Sakurai

Abstract: Sunspots are the most notable structure on the solar surface with strong magnetic fields. The field is generally strongest in a dark area (umbra), but sometimes stronger fields are found in non-dark regions such as a penumbra and a light bridge. The formation mechanism of such strong fields outside umbrae is still puzzling. Here we report clear evidence of the magnetic field of 6,250 G, which is t… ▽ More Sunspots are the most notable structure on the solar surface with strong magnetic fields. The field is generally strongest in a dark area (umbra), but sometimes stronger fields are found in non-dark regions such as a penumbra and a light bridge. The formation mechanism of such strong fields outside umbrae is still puzzling. Here we report clear evidence of the magnetic field of 6,250 G, which is the strongest field among Stokes I profiles with clear Zeeman splitting ever observed on the Sun. The field was almost parallel to the solar surface and located in a bright region sandwiched by two opposite-polarity umbrae. Using a time series of spectral datasets, we discussed the formation process of the super-strong field and suggested that this strong field region was generated as a result of compression of one umbra pushed by the horizontal flow from the other umbra, like the subduction of the Earth's crust in plate tectonics. △ Less

Submitted 22 December, 2017; originally announced December 2017.

Comments: 11 pages, 5 figures. Accepted for publication in ApJL

The Impact of Galactic Disc Environment on Star-Forming Clouds

Authors: Ngan K. Nguyen, Alex R. Pettitt, Elizabeth J. Tasker, Takashi Okamoto

Abstract: We explore the effect of different galactic disc environments on the properties of star-forming clouds through variations in the background potential in a set of isolated galaxy simulations. Rising, falling and flat rotation curves expected in halo dominated, disc dominated and Milky Way-like galaxies were considered, with and without an additional two-arm spiral potential. The evolution of each d… ▽ More We explore the effect of different galactic disc environments on the properties of star-forming clouds through variations in the background potential in a set of isolated galaxy simulations. Rising, falling and flat rotation curves expected in halo dominated, disc dominated and Milky Way-like galaxies were considered, with and without an additional two-arm spiral potential. The evolution of each disc displayed notable variations that are attributed to different regimes of stability, determined by shear and gravitational collapse. The properties of a typical cloud were largely unaffected by the changes in rotation curve, but the production of small and large cloud associations was strongly dependent on this environment. This suggests that while differing rotation curves can influence where clouds are initially formed, the average bulk properties are effectively independent of the global environment. The addition of a spiral perturbation made the greatest difference to cloud properties, successfully sweeping the gas into larger, seemingly unbound, extended structures and creating large arm-interarm contrasts. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: Accepted to MNRAS on 3rd December, 2017

Hydrocode modeling of the spallation process during hypervelocity impacts: Implications for the ejection of Martian meteorites

Authors: Kosuke Kurosawa, Takaya Okamoto, Hidenori Genda

Abstract: Hypervelocity ejection of material by impact spallation is considered a plausible mechanism for material exchange between two planetary bodies. We have modeled the spallation process during vertical impacts over a range of impact velocities from 6 to 21 km/s using both grid- and particle-based hydrocode models. The Tillotson equations of state, which are able to treat the nonlinear dependence of d… ▽ More Hypervelocity ejection of material by impact spallation is considered a plausible mechanism for material exchange between two planetary bodies. We have modeled the spallation process during vertical impacts over a range of impact velocities from 6 to 21 km/s using both grid- and particle-based hydrocode models. The Tillotson equations of state, which are able to treat the nonlinear dependence of density on pressure and thermal pressure in the strongly shocked matter, were used to study the hydrodynamic and thermodynamic response after impacts. The effects of material strength and gravitational acceleration were not considered. A two-dimensional time-dependent pressure field within a 1.5-fold projectile radius from the impact point was investigated in cylindrical coordinates to address the generation of spalled material. A resolution test was also performed to reject ejected materials with peak pressures that were too low due to artificial viscosity. The relationship between ejection velocity veject and peak pressure Ppeak was also derived. Our approach shows that late stage acceleration in an ejecta curtain occurs due to the compressible nature of the ejecta, resulting in an ejection velocity that can be higher than the ideal maximum of the resultant particle velocity after passage of a shock wave. We also calculate the ejecta mass that can escape from a planet like Mars (i.e., veject higher than 5 km/s) that matches the petrographic constraints from Martian meteorites, and which occurs when Ppeak from 30-50 GPa. Although the mass of such ejecta is limited to from 0.1-1 percent of the projectile mass in vertical impacts, this is sufficient for spallation to have been a plausible mechanism for the ejection of Martian meteorites. Finally, we propose that impact spallation is a plausible mechanism for the generation of tektites. △ Less

Submitted 2 September, 2017; originally announced September 2017.

Comments: 67 pages, 28 figures, accepted for publication in Icarus

The effect of radiation pressure on spatial distribution of dust inside HII regions

Authors: Shohei Ishiki, Takashi Okamoto, Akio K. Inoue

Abstract: We investigate the impact of radiation pressure on spatial dust distribution inside H$_\mathrm{II}$ regions using one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. In order to investigate grain size effects as well, we introduce two additional fluid components describing large and small dust grains in the simulations. Relative velocity… ▽ More We investigate the impact of radiation pressure on spatial dust distribution inside H$_\mathrm{II}$ regions using one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. In order to investigate grain size effects as well, we introduce two additional fluid components describing large and small dust grains in the simulations. Relative velocity between dust and gas strongly depends on the drag force. We include collisional drag force and coulomb drag force. We find that, in a compact H$_\mathrm{II}$ region, a dust cavity region is formed by radiation pressure. Resulting dust cavity sizes (~0.2 pc) agree with observational estimates reasonably well. Since dust inside an H$_\mathrm{II}$ region is strongly charged, relative velocity between dust and gas is mainly determined by the coulomb drag force. Strength of the coulomb drag force is about 2-order of magnitude larger than that of the collisional drag force. In addition, in a cloud of mass $10^5$ $M_{\odot}$, we find that the radiation pressure changes the grain size distribution inside H$_\mathrm{II}$ regions. Since large (0.1 $\mathrm{μm}$) dust grains are accelerated more efficiently than small (0.01 $\mathrm{μm}$) grains, the large to small grain mass ratio becomes smaller by an order of magnitude compared with the initial one. Resulting dust size distributions depend on the luminosity of the radiation source. The large and small grain segregation becomes weaker when we assume stronger radiation source, since dust grain charges become larger under stronger radiation and hence coulomb drag force becomes stronger. △ Less

Submitted 20 November, 2017; v1 submitted 23 August, 2017; originally announced August 2017.

Comments: 10 pages, 4 figures, 3 tables, accepted for publication in MNRAS

Imprints of the super-Eddington accretion on the quasar clustering

Authors: Taira Oogi, Motohiro Enoki, Tomoaki Ishiyama, Masakazu A. R. Kobayashi, Ryu Makiya, Masahiro Nagashima, Takashi Okamoto, Hikari Shirakata

Abstract: Super-Eddington mass accretion has been suggested as an efficient mechanism to grow supermassive black holes (SMBHs). We investigate the imprint left by the radiative efficiency of the super-Eddington accretion process on the clustering of quasars using a new semi-analytic model of galaxy and quasar formation based on large-volume cosmological $N$-body simulations. Our model includes a simple mode… ▽ More Super-Eddington mass accretion has been suggested as an efficient mechanism to grow supermassive black holes (SMBHs). We investigate the imprint left by the radiative efficiency of the super-Eddington accretion process on the clustering of quasars using a new semi-analytic model of galaxy and quasar formation based on large-volume cosmological $N$-body simulations. Our model includes a simple model for the radiative efficiency of a quasar, which imitates the effect of photon trapping for a high mass accretion rate. We find that the model of radiative efficiency affects the relation between the quasar luminosity and the quasar host halo mass. The quasar host halo mass has only weak dependence on quasar luminosity when there is no upper limit for quasar luminosity. On the other hand, it has significant dependence on quasar luminosity when the quasar luminosity is limited by its Eddington luminosity. In the latter case, the quasar bias also depends on the quasar luminosity, and the quasar bias of bright quasars is in agreement with observations. Our results suggest that the quasar clustering studies can provide a constraint on the accretion disc model. △ Less

Submitted 21 June, 2017; originally announced June 2017.

Comments: 6 pages, 4 figures, to appear in MNRAS Letters

The metal enrichment of passive galaxies in cosmological simulations of galaxy formation

Authors: Takashi Okamoto, Masahiro Nagashima, Cedric G. Lacey, Carlos S. Frenk

Abstract: Massive early-type galaxies have higher metallicities and higher ratios of $α$ elements to iron than their less massive counterparts. Reproducing these correlations has long been a problem for hierarchical galaxy formation theory, both in semi-analytic models and cosmological hydrodynamic simulations. We show that a simulation in which gas cooling in massive dark haloes is quenched by radio-mode a… ▽ More Massive early-type galaxies have higher metallicities and higher ratios of $α$ elements to iron than their less massive counterparts. Reproducing these correlations has long been a problem for hierarchical galaxy formation theory, both in semi-analytic models and cosmological hydrodynamic simulations. We show that a simulation in which gas cooling in massive dark haloes is quenched by radio-mode active galactic nuclei (AGNs) feedback naturally reproduces the observed trend between $α$/Fe and the velocity dispersion of galaxies, $σ$. The quenching occurs earlier for more massive galaxies. Consequently, these galaxies complete their star formation before $α$/Fe is diluted by the contribution from type Ia supernovae. For galaxies more massive than $\sim 10^{11}~M_\odot$ whose $α$/Fe correlates positively with stellar mass, we find an inversely correlated mass-metallicity relation. This is a common problem in simulations in which star formation in massive galaxies is quenched either by quasar- or radio-mode AGN feedback. The early suppression of gas cooling in progenitors of massive galaxies prevents them from recapturing enriched gas ejected as winds. Simultaneously reproducing the [$α$/Fe]-$σ$ relation and the mass-metallicity relation is, thus, difficult in the current framework of galaxy formation. △ Less

Submitted 27 October, 2016; v1 submitted 20 October, 2016; originally announced October 2016.

Comments: 10 pages, 5 figures, accepted for publication in MNRAS, added references

Journal ref: Mon.Not.Roy.Astron.Soc. 464, 4866-4874, 2017

Helical motions of fine-structure prominence threads observed by Hinode and IRIS

Authors: Takenori J. Okamoto, Wei Liu, Saku Tsuneta

Abstract: Fine-structure dynamics in solar prominences holds critical clues to understanding their physical nature of significant space-weather implications. We report evidence of rotational motions of horizontal helical threads in two active-region prominences observed by the \emph{Hinode} and/or \emph{IRIS} satellites at high resolution. In the first event, we found transverse motions of brightening threa… ▽ More Fine-structure dynamics in solar prominences holds critical clues to understanding their physical nature of significant space-weather implications. We report evidence of rotational motions of horizontal helical threads in two active-region prominences observed by the \emph{Hinode} and/or \emph{IRIS} satellites at high resolution. In the first event, we found transverse motions of brightening threads at speeds up to 55~km~s$^{-1}$ seen in the plane of the sky. Such motions appeared as sinusoidal space--time trajectories with a typical period of $\sim$390~s, which is consistent with plane-of-sky projections of rotational motions. Phase delays at different locations suggest propagation of twists along the threads at phase speeds of 90--270~km~s$^{-1}$. At least 15 episodes of such motions occurred in two days, none associated with any eruption. For these episodes, the plane-of-sky speed is linearly correlated with the vertical travel distance, suggestive of a constant angular speed. In the second event, we found Doppler velocities of 30--40~km~s$^{-1}$ in opposite directions in the top and bottom portions of the prominence, comparable to the plane-of-sky speed. The moving threads have about twice broader line widths than stationary threads. These observations, when taken together, provide strong evidence for rotations of helical prominence threads, which were likely driven by unwinding twists triggered by magnetic reconnection between twisted prominence magnetic fields and ambient coronal fields. △ Less

Submitted 30 July, 2016; originally announced August 2016.

Comments: 25 pages, 11 figures, accepted for publication in ApJ

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

Radiation feedback in dusty clouds

Authors: Shohei Ishiki, Takashi Okamoto

Abstract: We have investigated the impact of photoionization and radiation pressure on a dusty star-forming cloud using one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. We find that, in a cloud of mass 10^5 Msun and radius 17 pc, the effect of radiation pressure is negligible when star formation efficiency is 2%. The importance of radiation pre… ▽ More We have investigated the impact of photoionization and radiation pressure on a dusty star-forming cloud using one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. We find that, in a cloud of mass 10^5 Msun and radius 17 pc, the effect of radiation pressure is negligible when star formation efficiency is 2%. The importance of radiation pressure increases with increasing star formation efficiency or an increasing dust-to-gas mass ratio. The net effect of radiation feedback, however, becomes smaller with the increasing dust-to-gas mass ratio, since the absorption of ultra-violet photons by dust grains suppresses photoionization and hence photoheating. △ Less

Submitted 19 December, 2016; v1 submitted 9 May, 2016; originally announced May 2016.

Comments: 6 pages, 1 figure, 1 table; largely updated from ver 1., accepted for publication in MNRAS

Theoretical reevaluations of the black hole mass - bulge mass relation - I. Effect of the seed black hole mass

Authors: Hikari Shirakata, Toshihiro Kawaguchi, Takashi Okamoto, Ryu Makiya, Tomoaki Ishiyama, Yoshiki Matsuoka, Masahiro Nagashima, Motohiro Enoki, Taira Oogi, Masakazu A. R. Kobayashi

Abstract: We explore the effect of varying the mass of the seed black hole on the resulting black hole mass - bulge mass relation at z ~ 0, using a semi-analytic model of galaxy formation combined with large cosmological N-body simulations. We constrain our model by requiring the observed properties of galaxies at z ~ 0 are reproduced. In keeping with previous semi-analytic models, we place a seed black hol… ▽ More We explore the effect of varying the mass of the seed black hole on the resulting black hole mass - bulge mass relation at z ~ 0, using a semi-analytic model of galaxy formation combined with large cosmological N-body simulations. We constrain our model by requiring the observed properties of galaxies at z ~ 0 are reproduced. In keeping with previous semi-analytic models, we place a seed black hole immediately after a galaxy forms. When the mass of the seed is set at 10^5 M_sun, we find that the model results become inconsistent with recent observational results of the black hole mass - bulge mass relation for dwarf galaxies. In particular, the model predicts that bulges with ~ 10^9 M_sun harbour larger black holes than observed. On the other hand, when we employ seed black holes with 10^3 M_sun, or randomly select their mass within a 10^(3-5) M_sun range, the resulting relation is consistent with observation estimates, including the observed dispersion. We find that to obtain stronger constraints on the mass of seed black holes, observations of less massive bulges at z ~ 0 are a more powerful comparison than the relations at higher redshifts. △ Less

Submitted 22 June, 2016; v1 submitted 18 April, 2016; originally announced April 2016.

Comments: 6 pages, 2 tables, 4 figures. Accepted for publication by MNRAS Letter

Nebular line emission from z > 7 galaxies in cosmological simulations: rest-frame UV to Optical lines

Authors: Ikkoh Shimizu, Akio K. Inoue, Naoki Yoshida, Takashi Okamoto

Abstract: We have performed very large and high resolution cosmological hydrodynamic simulations in order to investigate detectability of nebular lines in the rest-frame ultraviolet (UV) to optical wavelength range from galaxies at $z > 7$. We find that the expected line fluxes are very well correlated with apparent UV magnitudes. The C IV $1549 \rm Å$ and the C III] $1909 \rm Å$ lines of galaxies brighter… ▽ More We have performed very large and high resolution cosmological hydrodynamic simulations in order to investigate detectability of nebular lines in the rest-frame ultraviolet (UV) to optical wavelength range from galaxies at $z > 7$. We find that the expected line fluxes are very well correlated with apparent UV magnitudes. The C IV $1549 \rm Å$ and the C III] $1909 \rm Å$ lines of galaxies brighter than $26~ \rm AB$ are detectable with current facilities such as the VLT/XShooter and the Keck/MOSFIRE. Metal lines such as C IV $1549 \rm Å$, C III] $1909 \rm Å$, [O II] $3727 \rm Å$ and [O III] $4959/5007 \rm Å$ can be good targets for the spectroscopic observation with Thirty Meter Telescope (TMT), European Extremely Large Telescope (E-ELT), Giant Magellan Telescope (GMT) and James Webb Space Telescope (JWST). H$α$ and H$β$ lines are also expected to be detectable with these telescopes. We also predict detectability of the nebular lines for $z > 10$ galaxies that will be found with JWST, Wide-Field Infrared Survey Telescope (WFIRST) and First Light And Reionization Explorer (FLARE) $(11 \leq z \leq 15$). We conclude that the C IV $1549 \rm Å$, C III] $1909 \rm Å$, [O III] $4959/5007 \rm Å$ and H$β$ lines even from $z \sim 15$ galaxies can be strong targets for TMT, ELT and JWST. We also find that the magnification by gravitational lensing is of great help to detect such high-$z$ galaxies. According to our model, the C III] $1909 \rm Å$ line in $z > 9$ galaxy candidates is detectable even using the current facilities. △ Less

Submitted 11 June, 2016; v1 submitted 2 September, 2015; originally announced September 2015.

Comments: 15 pages, 15 figures, 1 table. Accepted for publication to MNRAS

The New Numerical Galaxy Catalog ($ν^2$GC): An Updated Semi-analytic Model of Galaxy and AGN with Large Cosmological N-body Simulation

Authors: Ryu Makiya, Motohiro Enoki, Tomoaki Ishiyama, Masakazu A. R. Kobayashi, Masahiro Nagashima, Takashi Okamoto, Katsuya Okoshi, Taira Oogi, Hikari Shirakata

Abstract: We present a new cosmological galaxy formation model, $ν^2$GC, as an updated version of our previous model $ν$GC. We adopt the so-called "semi-analytic" approach, in which the formation history of dark matter halos is computed by ${\it N}$-body simulations, while the baryon physics such as gas cooling, star formation and supernova feedback are simply modeled by phenomenological equations. Major up… ▽ More We present a new cosmological galaxy formation model, $ν^2$GC, as an updated version of our previous model $ν$GC. We adopt the so-called "semi-analytic" approach, in which the formation history of dark matter halos is computed by ${\it N}$-body simulations, while the baryon physics such as gas cooling, star formation and supernova feedback are simply modeled by phenomenological equations. Major updates of the model are as follows: (1) the merger trees of dark matter halos are constructed in state-of-the-art ${\it N}$-body simulations, (2) we introduce the formation and evolution process of supermassive black holes and the suppression of gas cooling due to active galactic nucleus (AGN) activity, (3) we include heating of the intergalactic gas by the cosmic UV background, and (4) we tune some free parameters related to the astrophysical processes using a Markov chain Monte Carlo method. Our ${\it N}$-body simulations of dark matter halos have unprecedented box size and mass resolution (the largest simulation contains 550 billion particles in a 1.12 Gpc/h box), enabling the study of much smaller and rarer objects. The model was tuned to fit the luminosity functions of local galaxies and mass function of neutral hydrogen. Local observations, such as the Tully-Fisher relation, size-magnitude relation of spiral galaxies and scaling relation between the bulge mass and black hole mass were well reproduced by the model. Moreover, the model also well reproduced the cosmic star formation history and the redshift evolution of rest-frame ${\it K}$-band luminosity functions. The numerical catalog of the simulated galaxies and AGNs is publicly available on the web. △ Less

Submitted 22 October, 2016; v1 submitted 28 August, 2015; originally announced August 2015.

Comments: 25 pages, 21 figures, 2 tables, accepted to PASJ. The numerical catalog of galaxies and AGNs are available at http://cdsarc.u-strasbg.fr/cgi-bin/VizieR?-source=J/PASJ/68/25