Showing 1–50 of 70 results for author: Nagashima, M

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)

The History of The Milky Way: The Evolution of Star Formation, Cosmic Rays, Metallicity, and Stellar Dynamics over Cosmic Time

Authors: Jiro Shimoda, Shu-ichiro Inutsuka, Masahiro Nagashima

Abstract: We study the long-term evolution of the Milky Way (MW) over cosmic time by modeling the star formation, cosmic rays, metallicity, stellar dynamics, outflows and inflows of the galactic system to obtain various insights into the galactic evolution. The mass accretion is modeled by the results of cosmological N-body simulations for the cold dark matter. We find that the star formation rate is about… ▽ More We study the long-term evolution of the Milky Way (MW) over cosmic time by modeling the star formation, cosmic rays, metallicity, stellar dynamics, outflows and inflows of the galactic system to obtain various insights into the galactic evolution. The mass accretion is modeled by the results of cosmological N-body simulations for the cold dark matter. We find that the star formation rate is about half the mass accretion rate of the disk, given the consistency between observed Galactic Diffuse X-ray Emissions (GDXEs) and possible conditions driving the Galactic wind. Our model simultaneously reproduces the quantities of star formation rate, cosmic rays, metals, and the rotation curve of the current MW. The most important predictions of the model are that there is an unidentified accretion flow with a possible number density of $\sim10^{-2}$ cm$^{-3}$ and the part of the GDXEs originates from a hot, diffuse plasma which is formed by consuming about 10 % of supernova explosion energy. The latter is the science case for future X-ray missions; XRISM, Athena, and so on. We also discuss further implications of our results for the planet formation and observations of externalgalaxies in terms of the multimessenger astronomy. △ Less

Submitted 8 December, 2023; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 18 pages, 11 figures, accepted for publication in PASJ

Uchuu-$ν^2$GC galaxies and AGN: Cosmic variance forecasts of high-redshift AGN for JWST, Euclid, and LSST

Authors: Taira Oogi, Tomoaki Ishiyama, Francisco Prada, Manodeep Sinha, Darren Croton, Sofía A. Cora, Eric Jullo, Anatoly A. Klypin, Masahiro Nagashima, J. López Cacheiro, José Ruedas, Masakazu A. R. Kobayashi, Ryu Makiya

Abstract: Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-bo… ▽ More Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$ν^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$ν^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST). △ Less

Submitted 5 August, 2023; v1 submitted 29 July, 2022; originally announced July 2022.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS, Figs. 3, 4, 11, 13, and 14 revised, Figs. 14 and 16 in the original manuscript moved to the Appendix, All Uchuu-$ν^2$GC data used in this paper are publicly available on our Skies & Universes website https://www.skiesanduniverses.org/Simulations/Uchuu/GalaxyCatalogues/

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

Quantifying the effect of field variance on the H$α$ luminosity function with the New Numerical Galaxy Catalogue ($ν^2$GC)

Authors: Kazuyuki Ogura, Masahiro Nagashima, Rhythm Shimakawa, Masao Hayashi, Masakazu A. R. Kobayashi, Taira Oogi, Tomoaki Ishiyama, Yusei Koyama, Ryu Makiya, Katsuya Okoshi, Masato Onodera, Hikari Shirakata

Abstract: We construct a model of H$α$ emitters (HAEs) based on a semi-analytic galaxy formation model, the New Numerical Galaxy Catalog ($ν^2$GC). In this paper, we report our estimate for the field variance of the HAE distribution. By calculating the H$α$ luminosity from the star-formation rate of galaxies, our model well reproduces the observed H$α$ luminosity function (LF) at $z=0.4$. The large volume o… ▽ More We construct a model of H$α$ emitters (HAEs) based on a semi-analytic galaxy formation model, the New Numerical Galaxy Catalog ($ν^2$GC). In this paper, we report our estimate for the field variance of the HAE distribution. By calculating the H$α$ luminosity from the star-formation rate of galaxies, our model well reproduces the observed H$α$ luminosity function (LF) at $z=0.4$. The large volume of the $ν^2$GC makes it possible to examine the spatial distribution of HAEs over a region of (411.8 Mpc)$^3$ in the comoving scale. The surface number density of $z=0.4$ HAEs with $L_{\rm Hα} \geq 10^{40}$ erg s$^{-1}$ is 308.9 deg$^{-2}$. We have confirmed that the HAE is a useful tracer for the large-scale structure of the Universe because of their significant overdensity ($>$ 5$σ$) at clusters and the filamentary structures. The H$α$ LFs within a survey area of $\sim$2 deg$^2$ (typical for previous observational studies) show a significant field variance up to $\sim$1 dex. Based on our model, one can estimate the variance on the H$α$ LFs within given survey areas. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 13 pages, 11 figures, accepted for publication in the ApJ

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

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

New constraints on red-spiral galaxies from their kinematics in clusters of galaxies

Authors: Akinari Hamabata, Taira Oogi, Masamune Oguri, Takahiro Nishimichi, Masahiro Nagashima

Abstract: The distributions of the pairwise line-of-sight velocity between galaxies and their host clusters are segregated according to the galaxy's colour and morphology. We investigate the velocity distribution of red-spiral galaxies, which represents a rare population within galaxy clusters. We find that the probability distribution function of the pairwise line-of-sight velocity $v_{\rm{los}}$ between r… ▽ More The distributions of the pairwise line-of-sight velocity between galaxies and their host clusters are segregated according to the galaxy's colour and morphology. We investigate the velocity distribution of red-spiral galaxies, which represents a rare population within galaxy clusters. We find that the probability distribution function of the pairwise line-of-sight velocity $v_{\rm{los}}$ between red-spiral galaxies and galaxy clusters has a dip at $v_{\rm{los}} = 0$, which is a very odd feature, at 93\% confidence level. To understand its origin, we construct a model of the phase space distribution of galaxies surrounding galaxy clusters in three-dimensional space by using cosmological $N$-body simulations. We adopt a two component model that consists of the infall component, which corresponds to galaxies that are now falling into galaxy clusters, and the splashback component, which corresponds to galaxies that are on their first (or more) orbit after falling into galaxy clusters. We find that we can reproduce the distribution of the line-of-sight velocity of red-spiral galaxies with the dip with a very simple assumption that red-spiral galaxies reside predominantly in the infall component, regardless of the choice of the functional form of their spatial distribution. Our results constrain the quenching timescale of red-spiral galaxies to a few Gyrs, and the radius where the morphological transformation is effective as $r \sim 0.2 h^{-1} \rm{Mpc}$. △ Less

Submitted 17 July, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

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

Report number: YITP-19-07

Lethal radiation from nearby supernovae helps to explain the small cosmological constant

Authors: Tomonori Totani, Hidetoshi Omiya, Takahiro Sudoh, Masakazu A. R. Kobayashi, Masahiro Nagashima

Abstract: The observed value $Λ_{\rm obs}$ of the cosmological constant $Λ$ is extremely smaller than theoretical expectations, and the anthropic argument has been proposed as a solution to this problem because galaxies do not form when $Λ\gg Λ_{\rm obs}$. However, the contemporary galaxy formation theory predicts that stars form even with a high value of $Λ/ Λ_{\rm obs} \sim$ 50, which makes the anthropic… ▽ More The observed value $Λ_{\rm obs}$ of the cosmological constant $Λ$ is extremely smaller than theoretical expectations, and the anthropic argument has been proposed as a solution to this problem because galaxies do not form when $Λ\gg Λ_{\rm obs}$. However, the contemporary galaxy formation theory predicts that stars form even with a high value of $Λ/ Λ_{\rm obs} \sim$ 50, which makes the anthropic argument less persuasive. Here we calculate the probability distribution of $Λ$ using a model of cosmological galaxy formation, considering extinction of observers caused by radiation from nearby supernovae. The life survival probability decreases in a large $Λ$ universe because of higher stellar density. Using a reasonable rate of lethal supernovae, we find that the mean expectation value of $Λ$ can be close to $Λ_{\rm obs}$, and hence this effect may be essential to understand the small but nonzero value of $Λ$. It is predicted that we are located on the edge of habitable regions about stellar density in the Galaxy, which may be tested by future exoplanet studies. △ Less

Submitted 18 July, 2018; v1 submitted 27 April, 2018; originally announced April 2018.

Comments: 13 pages, 2 figures. the version accepted to Astrobiology

Journal ref: Astrobiology, 19, 126 (2019)

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

First results on the cluster galaxy population from the Subaru Hyper Suprime-Cam survey. II. Faint end color-magnitude diagrams and radial profiles of red and blue galaxies at $0.1<z<1.1$

Authors: Atsushi J. Nishizawa, Masamue Oguri, Taira Oogi, Surhud More, Takahiro Nishimichi, Masahiro Nagashima, Yen-Ting Lin, Rachel Mandelbaum, Masahiro Takada, Neta Bahcall, Jean Coupon, Song Huang, Hung-Yu Jian, Yutaka Komiyama, Alexie Leauthaud, Lihwai Lin, Hironao Miyatake, Satoshi Miyazaki, Masayuki Tanaka

Abstract: We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies do… ▽ More We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $σ_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Comments: 18pages, 10 figures, accepted as PASJ special issue

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

Clustering of quasars in a wide luminosity range at redshift 4 with Subaru Hyper Suprime-Cam wide field imaging

Authors: Wanqiu He, Masayuki Akiyama, James Bosch, Motohiro Enoki, Yuichi Harikane, Hiroyuki Ikeda, Nobunari Kashikawa, Toshihiro Kawaguchi, Yutaka Komiyama, Chien-Hsiu Lee, Yoshiki Matsuoka, Satoshi Miyazaki, Tohru Nagao, Masahiro Nagashima, Mana Niida, Atsushi J Nishizawa, Masamune Oguri, Masafusa Onoue, Taira Oogi, Masami Ouchi, Andreas Schulze, Yuji Shirasaki, John D. Silverman, Manobu M. Tanaka, Masayuki Tanaka , et al. (3 additional authors not shown)

Abstract: We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (S… ▽ More We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 20 pages, 10 figures, submitted to PASJ

Mid-infrared Spectroscopic Observations of the Dust-forming Classical Nova V2676 Oph

Authors: Hideyo Kawakita, Takafumi Ootsubo, Akira Arai, Yoshiharu Shinnaka, Masayoshi Nagashima

Abstract: The dust-forming nova V2676 Oph is unique in that it was the first nova to provide evidence of C_2 and CN molecules during its near-maximum phase and evidence of CO molecules during its early decline phase. Observations of this nova have revealed the slow evolution of its lightcurves and have also shown low isotopic ratios of carbon (12C/13C) and nitrogen (14N/15N) in its nova envelope. These beha… ▽ More The dust-forming nova V2676 Oph is unique in that it was the first nova to provide evidence of C_2 and CN molecules during its near-maximum phase and evidence of CO molecules during its early decline phase. Observations of this nova have revealed the slow evolution of its lightcurves and have also shown low isotopic ratios of carbon (12C/13C) and nitrogen (14N/15N) in its nova envelope. These behaviors indicate that the white dwarf (WD) star hosting V2676 Oph is a CO-rich WD rather than an ONe-rich WD (typically larger in mass than the former). We performed mid-infrared spectroscopic and photometric observations of V2676 Oph in 2013 and 2014 (respectively 452 and 782 days after its discovery). No significant [Ne II] emission at 12.8 micron was detected at either epoch. These provided evidence for a CO-rich WD star hosting V2676 Oph. Both carbon-rich and oxygen-rich grains were detected in addition to an unidentified infrared feature at 11.4 micron originating from polycyclic aromatic hydrocarbon molecules or hydrogenated amorphous carbon grains in the envelope of V2676 Oph. △ Less

Submitted 22 January, 2017; originally announced January 2017.

Comments: 14 pages including 2 figures and 4 tables, to appear in the Astronomical Journal, Volume 153, Number 2

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

Testing anthropic reasoning for the cosmological constant with a realistic galaxy formation model

Authors: Takahiro Sudoh, Tomonori Totani, Ryu Makiya, Masahiro Nagashima

Abstract: The anthropic principle is one of the possible explanations for the cosmological constant ($Λ$) problem. In previous studies, a dark halo mass threshold comparable with our Galaxy must be assumed in galaxy formation to get a reasonably large probability of finding the observed small value, $P(<$$Λ_{\rm obs})$, though stars are found in much smaller galaxies as well. Here we examine the anthropic a… ▽ More The anthropic principle is one of the possible explanations for the cosmological constant ($Λ$) problem. In previous studies, a dark halo mass threshold comparable with our Galaxy must be assumed in galaxy formation to get a reasonably large probability of finding the observed small value, $P(<$$Λ_{\rm obs})$, though stars are found in much smaller galaxies as well. Here we examine the anthropic argument by using a semi-analytic model of cosmological galaxy formation, which can reproduce many observations such as galaxy luminosity functions. We calculate the probability distribution of $Λ$ by running the model code for a wide range of $Λ$, while other cosmological parameters and model parameters for baryonic processes of galaxy formation are kept constant. Assuming that the prior probability distribution is flat per unit $Λ$, and that the number of observers is proportional to stellar mass, we find $P(<$$Λ_{\rm obs}) = 6.7 \%$ without introducing any galaxy mass threshold. We also investigate the effect of metallicity; we find $P(<$$Λ_{\rm obs}) = 9.0 \%$ if observers exist only in galaxies whose metallicity is higher than the solar abundance. If the number of observers is proportional to metallicity, we find $P(<$$Λ_{\rm obs}) = 9.7 \%$. Since these probabilities are not extremely small, we conclude that the anthropic argument is a viable explanation, if the value of $Λ$ observed in our universe is determined by a probability distribution. △ Less

Submitted 21 September, 2016; v1 submitted 1 July, 2016; originally announced July 2016.

Comments: 6 pages, 4 figures, 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

Quasar clustering in a galaxy and quasar formation model based on ultra high-resolution N-body simulations

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

Abstract: We investigate clustering properties of quasars using a new version of our semi-analytic model of galaxy and quasar formation with state-of-the-art cosmological N-body simulations. In this study, we assume that a major merger of galaxies triggers cold gas accretion on to a supermassive black hole and quasar activity. Our model can reproduce the downsizing trend of the evolution of quasars. We find… ▽ More We investigate clustering properties of quasars using a new version of our semi-analytic model of galaxy and quasar formation with state-of-the-art cosmological N-body simulations. In this study, we assume that a major merger of galaxies triggers cold gas accretion on to a supermassive black hole and quasar activity. Our model can reproduce the downsizing trend of the evolution of quasars. We find that the median mass of quasar host dark matter haloes increases with cosmic time by an order of magnitude from z=4 (a few 1e+11 Msun) to z=1 (a few 1e+12 Msun), and depends only weakly on the quasar luminosity. Deriving the quasar bias through the quasar--galaxy cross-correlation function in the model, we find that the quasar bias does not depend on the quasar luminosity, similar to observed trends. This result reflects the fact that quasars with a fixed luminosity have various Eddington ratios and thus have various host halo masses that primarily determine the quasar bias. We also show that the quasar bias increases with redshift, which is in qualitative agreement with observations. Our bias value is lower than the observed values at high redshifts, implying that we need some mechanisms that make quasars inactive in low-mass haloes and/or that make them more active in high-mass haloes. △ Less

Submitted 1 December, 2015; originally announced December 2015.

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

Journal ref: Monthly Notices of the Royal Astronomical Society: Letters, Volume 456, Issue 1, p.L30-L34 (2016)

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

The $ν^2$GC Simulations : Quantifying the Dark Side of the Universe in the Planck Cosmology

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

Abstract: We present the evolution of dark matter halos in six large cosmological N-body simulations, called the $ν^2$GC (New Numerical Galaxy Catalog) simulations on the basis of the LCDM cosmology consistent with observational results obtained by the Planck satellite. The largest simulation consists of $8192^3$ (550 billion) dark matter particles in a box of $1.12 \, h^{-1} \rm Gpc$ (a mass resolution of… ▽ More We present the evolution of dark matter halos in six large cosmological N-body simulations, called the $ν^2$GC (New Numerical Galaxy Catalog) simulations on the basis of the LCDM cosmology consistent with observational results obtained by the Planck satellite. The largest simulation consists of $8192^3$ (550 billion) dark matter particles in a box of $1.12 \, h^{-1} \rm Gpc$ (a mass resolution of $2.20 \times 10^{8} \, h^{-1} M_{\odot}$). Among simulations utilizing boxes larger than $1 \, h^{-1} \rm Gpc$, our simulation yields the highest resolution simulation that has ever been achieved. A $ν^2$GC simulation with the smallest box consists of eight billions particles in a box of $70 \, h^{-1} \rm Mpc$ (a mass resolution of $3.44 \times 10^{6} \, h^{-1} M_{\odot}$). These simulations can follow the evolution of halos over masses of eight orders of magnitude, from small dwarf galaxies to massive clusters. Using the unprecedentedly high resolution and powerful statistics of the $ν^2$GC simulations, we provide statistical results of the halo mass function, mass accretion rate, formation redshift, and merger statistics, and present accurate fitting functions for the Planck cosmology. By combining the $ν^2$GC simulations with our new semi-analytic galaxy formation model, we are able to prepare mock catalogs of galaxies and active galactic nuclei, which will be made publicly available in the near future. △ Less

Submitted 5 April, 2016; v1 submitted 9 December, 2014; originally announced December 2014.

Comments: 18 pages, 16 figures, accepted, simulation data can be downloaded on http://hpc.imit.chiba-u.jp/~nngc/

The Impact of Dust in Host Galaxies on Quasar Luminosity Functions

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

Abstract: We have investigated effects of dust attenuation on quasar luminosity functions using a semi-analytic galaxy formation model combined with a large cosmological N-body simulation. We estimate the dust attenuation of quasars self-consistently with that of galaxies by considering the dust in their host bulges. We find that the luminosity of the bright quasars is strongly dimmed by the dust attenuatio… ▽ More We have investigated effects of dust attenuation on quasar luminosity functions using a semi-analytic galaxy formation model combined with a large cosmological N-body simulation. We estimate the dust attenuation of quasars self-consistently with that of galaxies by considering the dust in their host bulges. We find that the luminosity of the bright quasars is strongly dimmed by the dust attenuation, about 2 mag in the B-band. Assuming the empirical bolometric corrections for active galactic nuclei (AGNs) by Marconi et al., we find that this dust attenuation is too strong to explain the B-band and X-ray quasar luminosity functions simultaneously. We consider two possible mechanisms that weaken the dust attenuation. As such a mechanism, we introduce a time delay for AGN activity, that is, gas fueling to a central black hole starts some time after the beginning of the starburst induced by a major merger. The other is the anisotropy in the dust distribution. We find that in order to make the dust attenuation of the quasars negligible, either the gas accretion into the black holes has to be delayed at least three times the dynamical timescale of their host bulges or the dust covering factor is as small as 0.1. △ Less

Submitted 24 February, 2015; v1 submitted 3 December, 2014; originally announced December 2014.

Comments: 5 pages, 4 figures, 1 tables. Accepted for publication by MNRAS Letters

Anti-hierarchical evolution of the Active Galactic Nucleus space density in a hierarchical universe

Authors: Motohiro Enoki, Tomoaki Ishiyama, Masakazu A. R. Kobayashi, Masahiro Nagashima

Abstract: Recent observations show that the space density of luminous active galactic nuclei (AGNs) peaks at higher redshifts than that of faint AGNs. This downsizing trend in the AGN evolution seems to be contradictory to the hierarchical structure formation scenario. In this study, we present the AGN space density evolution predicted by a semi-analytic model of galaxy and AGN formation based on the hierar… ▽ More Recent observations show that the space density of luminous active galactic nuclei (AGNs) peaks at higher redshifts than that of faint AGNs. This downsizing trend in the AGN evolution seems to be contradictory to the hierarchical structure formation scenario. In this study, we present the AGN space density evolution predicted by a semi-analytic model of galaxy and AGN formation based on the hierarchical structure formation scenario. We demonstrate that our model can reproduce the downsizing trend of the AGN space density evolution. The reason for the downsizing trend in our model is a combination of the cold gas depletion as a consequence of star formation, the gas cooling suppression in massive halos and the AGN lifetime scaling with the dynamical timescale. We assume that a major merger of galaxies causes a starburst, spheroid formation, and cold gas accretion onto a supermassive black hole (SMBH). We also assume that this cold gas accretion triggers AGN activity. Since the cold gas is mainly depleted by star formation and gas cooling is suppressed in massive dark halos, the amount of cold gas accreted onto SMBHs decreases with cosmic time. Moreover, AGN lifetime increases with cosmic time. Thus, at low redshifts, major mergers do not always lead to luminous AGNs. Because the luminosity of AGNs is correlated with the mass of accreted gas onto SMBHs, the space density of luminous AGNs decreases more quickly than that of faint AGNs. We conclude that the anti-hierarchical evolution of the AGN space density is not contradictory to the hierarchical structure formation scenario. △ Less

Submitted 1 October, 2014; v1 submitted 16 August, 2014; originally announced August 2014.

Comments: 10 pages, 5 figures. Accepted for publication in ApJ. Minor editorial changes and typos corrected

Journal ref: 2014 ApJ, 794, 69

Galaxy luminosity function and its cosmological evolution: Testing a new feedback model depending on galaxy-scale dust opacity

Authors: Ryu Makiya, Tomonori Totani, Masakazu A. R. Kobayashi, Masahiro Nagashima, Tsutomu T. Takeuchi

Abstract: We present a new version of a semi-analytic model of cosmological galaxy formation, incorporating a star formation law with a feedback depending on the galaxy-scale mean dust opacity and metallicity, motivated by recent observations of star formation in nearby galaxies and theoretical considerations. This new model is used to investigate the effect of such a feedback on shaping the galaxy luminosi… ▽ More We present a new version of a semi-analytic model of cosmological galaxy formation, incorporating a star formation law with a feedback depending on the galaxy-scale mean dust opacity and metallicity, motivated by recent observations of star formation in nearby galaxies and theoretical considerations. This new model is used to investigate the effect of such a feedback on shaping the galaxy luminosity function and its evolution. Star formation activity is significantly suppressed in dwarf galaxies by the new feedback effect, and the faint-end slope of local luminosity functions can be reproduced with a reasonable strength of supernova feedback, which is in contrast to the previous models that require a rather extreme strength of supernova feedback. Our model can also reproduce the early appearance of massive galaxies manifested in the bright-end of high redshift K-band luminosity functions. Though some of the previous models also succeeded in reproducing this, they assumed a star formation law depending on the galaxy-scale dynamical time, which is not supported by observations. We argue that the feedback depending on dust opacity (or metal column density) is essential, rather than that simply depending on gas column density, to get these results. △ Less

Submitted 19 March, 2014; v1 submitted 1 February, 2014; originally announced February 2014.

Comments: 10 pages, 7 figures, Accepted for publication in MNRAS. v2: minor revision, new references added

Infrared Spectral Energy Distribution of Galaxies in the AKARI All Sky Survey: Correlations with Galaxy Properties, and Their Physical Origin

Authors: Tomonori Totani, Tsutomu T. Takeuchi, Masahiro Nagashima, Masakazu A. R. Kobayashi, Ryu Makiya

Abstract: We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm observations. Our goal is to understand the physics determining the infrared spectral energy distribution (SED). The ratio of the total infrared luminosity L_TIR, to the star-formation rate (SFR) i… ▽ More We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm observations. Our goal is to understand the physics determining the infrared spectral energy distribution (SED). The ratio of the total infrared luminosity L_TIR, to the star-formation rate (SFR) is tightly correlated by a power-law to specific SFR (SSFR), and L_TIR is a good SFR indicator only for galaxies with the largest SSFR. We discovered a tight linear correlation for normal galaxies between the radiation field strength of dust heating, estimated by infrared SED fits (U_h), and that of galactic-scale infrared emission (U_TIR ~ L_TIR/R^2), where R is the optical size of a galaxy. The dispersion of U_h along this relation is 0.3 dex, corresponding to 13% dispersion in the dust temperature. This scaling and the U_h/U_TIR ratio can be explained physically by a thin layer of heating sources embedded in a thicker, optically-thick dust screen. The data also indicate that the heated fraction of the total dust mass is anti-correlated to the dust column density, supporting this interpretation. In the large U_TIR limit, the data of circumnuclear starbursts indicate the existence of an upper limit on U_h, corresponding to the maximum SFR per gas mass of ~ 10 Gyr^{-1}. We find that the number of galaxies sharply drops when they become optically thin against dust-heating radiation, suggesting that a feedback process to galaxy formation (likely by the photoelectric heating) is working when dust-heating radiation is not self-shielded on a galactic scale. Implications are discussed for the M_HI-size relation, the Kennicutt-Schmidt relation, and galaxy formation in the cosmological context. △ Less

Submitted 26 January, 2012; v1 submitted 28 March, 2011; originally announced March 2011.

Comments: 29 pages including 28 figures. matches the published version (PASJ 2011 Dec. 25 issue). The E-open option was chosen for this article, i.e., the official version available from PASJ site (http://pasj.asj.or.jp/v63/n6/630613/630613-frame.html) without restriction

Journal ref: Pub. Astron. Soc. Jpn, 63, 1181-1206 (2011)

Lyman Alpha Emitters at z=7 in the Subaru/XMM-Newton Deep Survey Field: Photometric Candidates and Luminosity Function

Authors: Kazuaki Ota, Masanori Iye, Nobunari Kashikawa, Kazuhiro Shimasaku, Masami Ouchi, Tomonori Totani, Masakazu A. R. Kobayashi, Masahiro Nagashima, Atsushi Harayama, Natsuki Kodaka, Tomoki Morokuma, Hisanori Furusawa, Akito Tajitsu, Takashi Hattori

Abstract: We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the S… ▽ More We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the Subaru Deep Field that led to the discovery of a z=6.96 LAE, we detected three probable z=7 LAE candidates. Even if all the candidates are real, the Lyman alpha luminosity function (LF) at z=7 shows a significant deficit from the LF at z=5.7 determined by previous surveys. The LAE number and Lyman alpha luminosity densities at z=7 is ~ 7.7-54% and ~5.5-39% of those at z=5.7 to the Lyman alpha line luminosity limit of L(Ly-alpha) >~ 9.2 x 10^{42} erg s^{-1}. This could be due to evolution of the LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly evolves from z=5.7 to 7. However, even after correcting for this effect of galaxy evolution on the decrease in LAE number density, the z=7 Lyman alpha LF still shows a deficit from z=5.7 LF. This might reflect the attenuation of Lyman alpha emission by neutral hydrogen remaining at the epoch of reionization and suggests that reionization of the universe might not be complete yet at z=7. If we attribute the density deficit to reionization, the intergalactic medium (IGM) transmission for Lyman alpha photons at z=7 would be 0.4 <= T_{Ly-alpha}^{IGM} <= 1, supporting the possible higher neutral fraction at the earlier epochs at z > 6 suggested by the previous surveys of z=5.7-7 LAEs, z ~ 6 quasars and z > 6 gamma-ray bursts. △ Less

Submitted 28 August, 2010; originally announced August 2010.

Comments: Accepted to ApJ for publication

HI-selected Galaxies as a probe of Quasar Absorption Systems

Authors: Katsuya Okoshi, Masahiro Nagashima, Naoteru Gouda, Yousuke Minowa

Abstract: We investigate the properties of HI-rich galaxies detected in blind radio surveys within the hierarchical structure formation scenario using a semi-analytic model of galaxy formation. By drawing a detailed comparison between the properties of HI-selected galaxies and HI absorption systems, we argue a link between the local galaxy population and quasar absorption systems, particularly for Damped… ▽ More We investigate the properties of HI-rich galaxies detected in blind radio surveys within the hierarchical structure formation scenario using a semi-analytic model of galaxy formation. By drawing a detailed comparison between the properties of HI-selected galaxies and HI absorption systems, we argue a link between the local galaxy population and quasar absorption systems, particularly for Damped Ly-alpha absorption (DLA) systems and sub-DLA systems. First, we evaluate how many HI-selected galaxies exhibit HI column densities as high as those of DLA systems. We find that HI-selected galaxies with HI masses M(HI) > 10^8 solar masses have gaseous disks that produce HI column densities comparable to those of DLA systems. We conclude that DLA galaxies where the HI column densities are as high as those of DLA systems, contribute significantly to the population of HI-selected galaxies at M(HI) > 10^8 solar masses. Second, we find that star formation rates (SFRs) correlate tightly with HI masses rather than B- (and J-) band luminosities. In the low-mass range M(HI) < 10^8 solar masses, sub-DLA galaxies replace DLA galaxies as the dominant population. The number fraction of sub-DLA galaxies relative to galaxies reaches 40%-60% at HI masses 10^8 solar masses and 30%-80% at 10^7 solar masses. The HI-selected galaxies at 10^7 solar masses are a strong probe of sub-DLA systems that place stringent constraints on galaxy formation and evolution. △ Less

Submitted 2 February, 2010; originally announced February 2010.

Comments: 25 pages, 13 figures, Accepted for publication in the Astrophysical Journal

Lyman Alpha Emitters in Hierarchical Galaxy Formation II. UV Continuum Luminosity Function and Equivalent Width Distribution

Authors: Masakazu A. R. Kobayashi, Tomonori Totani, Masahiro Nagashima

Abstract: We present theoretical predictions of UV continuum luminosity function (UV LF) and Lya equivalent width (EW) distribution of Lyman alpha emitters (LAEs) in the framework of the hierarchical clustering model of galaxy formation. The model parameters about LAEs were determined by fitting to the observed Lya LF at z=5.7 in our previous study, and the fit indicates that extinction of Lya photons by… ▽ More We present theoretical predictions of UV continuum luminosity function (UV LF) and Lya equivalent width (EW) distribution of Lyman alpha emitters (LAEs) in the framework of the hierarchical clustering model of galaxy formation. The model parameters about LAEs were determined by fitting to the observed Lya LF at z=5.7 in our previous study, and the fit indicates that extinction of Lya photons by dust is significantly less effective than that of UV continuum photons, implying clumpy dust distribution in interstellar medium. We then compare the predictions about UV LFs and EW distributions with a variety of observations at z ~ 3-6, allowing no more free parameters and paying careful attention to the selection conditions of LAEs in each survey. We find that the predicted UV LFs and EW distributions are in nice agreement with observed data, and especially, our model naturally reproduces the existence of large EW LAEs (> 240 A) without introducing Pop III stars or top-heavy initial mass function. We show that both the stellar population (young age and low metallicity) and extinction by clumpy dust are the keys to reproduce large EW LAEs. The evidence of EW enhancement by clumpy dust is further strengthened by the quantitative agreement between our model and recent observations about a positive correlation between EW and extinction. The observed trend that brighter LAEs in UV continuum tend to have smaller mean EW is also reproduced, and the clumpy dust is playing an important role again for this trend. We suggested in our previous study that the transmission of intergalactic medium for Lya emission rapidly decreases from z ~ 6 to 7 by the fitting to Lya LFs, and this evidence is quantitatively strengthened by the comparison with the UV LF and EW distribution at z ~ 6.6. △ Less

Submitted 16 November, 2009; v1 submitted 17 February, 2009; originally announced February 2009.

Comments: 21 pages, 13 figures, 1 table, accepted for publication in ApJ

Journal ref: Astrophys.J.708:1119-1134,2010

Dynamical response to supernova-induced gas removal in spiral galaxies with dark matter halo

Authors: Hiroko Koyama, Masahiro Nagashima, Takayuki Kakehata, Yuzuru Yoshii

Abstract: We investigate the dynamical response, in terms of disc size and rotation velocity, to mass loss by supernovae in the evolution of spiral galaxies. A thin baryonic disc having the Kuzmin density profile embedded in a spherical dark matter halo having a density profile proposed by Navarro, Frenk & White is considered. For a purpose of comparison, we also consider the homogeneous and $r^{-1}$ prof… ▽ More We investigate the dynamical response, in terms of disc size and rotation velocity, to mass loss by supernovae in the evolution of spiral galaxies. A thin baryonic disc having the Kuzmin density profile embedded in a spherical dark matter halo having a density profile proposed by Navarro, Frenk & White is considered. For a purpose of comparison, we also consider the homogeneous and $r^{-1}$ profiles for dark matter in a truncated spherical halo. Assuming for simplicity that the dark matter distribution is not affected by mass loss from discs and the change of baryonic disc matter distribution is homologous, we evaluate the effects of dynamical response in the resulting discs. We found that the dynamical response only for an adiabatic approximation of mass loss can simultaneously account for the rotation velocity and disc size as observed particularly in dwarf spiral galaxies, thus reproducing the Tully-Fisher relation and the size versus magnitude relation over the full range of magnitude. Furthermore, we found that the mean specific angular momentum in discs after the mass loss becomes larger than that before the mass loss, suggesting that the mass loss would occur most likely from the central disc region where the specific angular momentum is low. △ Less

Submitted 31 August, 2008; v1 submitted 10 June, 2008; originally announced June 2008.

Comments: 15 pages, 11 figures, final version for publication in MNRAS

The origin of globular cluster systems from cosmological simulations

Authors: K. Bekki, H. Yahagi, M. Nagashima, D. A. Forbes

Abstract: We investigate the structural, kinematical, and chemical properties of globular cluster systems (GCSs) in galaxies of different Hubble types in a self-consistent manner based on high-resolution cosmological N-body simulations combined with semi-analytic models of galaxy and globular cluster (GC) formation. We focus on correlations between the physical properties of GCSs and those of their host g… ▽ More We investigate the structural, kinematical, and chemical properties of globular cluster systems (GCSs) in galaxies of different Hubble types in a self-consistent manner based on high-resolution cosmological N-body simulations combined with semi-analytic models of galaxy and globular cluster (GC) formation. We focus on correlations between the physical properties of GCSs and those of their host galaxies for about 10^5 simulated galaxies located at the centres of dark matter halos (i.e. we do not consider satellite galaxies in sub-halos). Our principal results, which can be tested against observations, are as follows. The majority (about 90%) of GCs currently in halos are formed in low-mass galaxies at redshifts greater than 3 with mean formation redshifts of z = 5.7 (12.7 Gyrs ago) and 4.3 (12.3 Gyrs ago) for metal-poor GCs (MPC) and metal-rich GCs (MRCs), respectively. About 52 % of galaxies with GCs show clear bimodality in their metallicity distribution functions, though less luminous galaxies with M_B fainter than -17 are much less likely to show bimodality owing to little or no MRCs. The number fraction of MRCs does not depend on Hubble type but is generally smaller for less luminous galaxies. The specific frequencies (S_ N) of GCSs are typically higher in ellipticals (S_ N ~ 4.0) than in spirals (S_ N ~ 1.8), and higher again (S_N ~ 5.0) for galaxies located at the centers of clusters of galaxies. The total number of GCs per unit halo mass does not depend strongly on M_B or Hubble type of the host galaxy. The mean metallicities of MPCs and MRCs depend on M_B such that they are higher in more luminous galaxies, though the dependence is significantly weakerfor MPCs. △ Less

Submitted 11 April, 2008; originally announced April 2008.

Comments: 19 pages, 29 figures, accepted in MNRAS

Star formation history of galaxies from z=0 to z=0.7 A backward approach to the evolution of star-forming galaxies

Authors: V. Buat, S. Boissier, D. Burgarella, T. T. Takeuchi, E. Le Floc'h, D. Marcillac, J. Huang, M. Nagashima, M. Enoki

Abstract: We investigate whether the mean star formation activity of star-forming galaxies from z=0 to z=0.7 in the GOODS-S field can be reproduced by simple evolution models of these systems. In this case, such models might be used as first order references for studies at higher z to decipher when and to what extent a secular evolution is sufficient to explain the star formation history in galaxies. We… ▽ More We investigate whether the mean star formation activity of star-forming galaxies from z=0 to z=0.7 in the GOODS-S field can be reproduced by simple evolution models of these systems. In this case, such models might be used as first order references for studies at higher z to decipher when and to what extent a secular evolution is sufficient to explain the star formation history in galaxies. We selected star-forming galaxies at z=0 and at z=0.7 in IR and in UV to have access to all the recent star formation. We focused on galaxies with a stellar mass ranging between 10^{10} and 10^{11} M_sun for which the results are not biased by the selections. We compared the data to chemical evolution models developed for spiral galaxies and originally built to reproduce the main characteristics of the Milky Way and nearby spirals without fine-tuning them for the present analysis. We find a shallow decrease in the specific star formation rate (SSFR) when the stellar mass increases. The evolution of the SSFR characterizing both UV and IR selected galaxies from z=0 to z=0.7 is consistent with the models built to reproduce the present spiral galaxies. There is no need to strongly modify of the physical conditions in galaxies to explain the average evolution of their star formation from z=0 to z=0.7. We use the models to predict the evolution of the star formation rate and the metallicity on a wider range of redshift and we compare these predictions with the results of semi-analytical models. △ Less

Submitted 4 March, 2008; originally announced March 2008.

Comments: 14 pages, 10 figures. accepted for publication in Astronomy & Astrophysics

Formation of the Galactic globular clusters with He-rich stars in low-mass halos virialized at high redshift

Authors: K. Bekki, H. Yahagi, M. Nagashima, D. A. Forbes

Abstract: Recent observations have reported that the Galactic globular clusters (GCs) with unusually extended horizontal-branch (EHB) morphologies show a significantly lower velocity dispersion compared with that of the entire Galactic GC system. We consider that the observed distinctive kinematics of GCs with EHB has valuable information on the formation epochs of GCs and accordingly discuss this observa… ▽ More Recent observations have reported that the Galactic globular clusters (GCs) with unusually extended horizontal-branch (EHB) morphologies show a significantly lower velocity dispersion compared with that of the entire Galactic GC system. We consider that the observed distinctive kinematics of GCs with EHB has valuable information on the formation epochs of GCs and accordingly discuss this observational result based on cosmological N-body simulations with a model of GC formation. We assume that GCs in galaxies were initially formed in low-mass halos at high redshifts and we investigate final kinematics of GCs in their host halos at $z=0$. We find that GCs formed in halos virialized at z>10 show lower velocity dispersions on average than those formed at z>6 for halos with GCs at z=0. We thus suggest that the origin of the observed lower velocity dispersion for the Galactic GCs with EHBs is closely associated with earlier formation epochs (z>10) of halos initially hosting the GCs in the course of the Galaxy formation. Considering that the origin of EHBs can be due to the presence of helium-enhanced second-generation stars in GCs, we discuss the longstanding second parameter problem of GCs in the context of different degrees of chemical pollution in GC-forming gas clouds within low-mass halos virialized at different redshifts. △ Less

Submitted 22 September, 2007; originally announced September 2007.

Comments: 5 pages, 3 figures, accepted by MNRAS Letters

The Reionization and Galaxy Evolution Probed by z=7 Lyman Alpha Emitters

Authors: Kazuaki Ota, Masanori Iye, Nobunari Kashikawa, Kazuhiro Shimasaku, Masakazu A. R. Kobayashi, Tomonori Totani, Masahiro Nagashima, Tomoki Morokuma, Hisanori Furusawa, Takashi Hattori, Yuichi Matsuda, Tetsuya Hashimoto, Masami Ouchi

Abstract: We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function… ▽ More We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function of LAE is known to decline from z=5.7 to 6.6 in the SDF. L* at z=6.6 is 40-60% of that at z=5.7. We also confirm that the number density of z=7 LAE is only 17% of the density at z=6.6 comparing the latest SDF LAE samples. This series of significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break galaxies, and a LAE evolution model based on the hierarchical clustering, we find that galaxy evolution alone cannot explain all the decrease in density. This extra density deficit can be interpreted as the attenuation of the Lyman alpha photons from LAEs due to a rapid evolution of neutral hydrogen fraction during the ongoing cosmic reionization at z~6.6-7. △ Less

Submitted 22 April, 2008; v1 submitted 11 July, 2007; originally announced July 2007.

Comments: published in Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 677, Issue 1, pp. 12-26 (2008)

Lyman Alpha Emitters in the Hierarchically Clustering Galaxy Formation

Authors: Masakazu A. R. Kobayashi, Tomonori Totani, Masahiro Nagashima

Abstract: We present a new theoretical model for the luminosity functions (LFs) of Lyman alpha (Lya) emitting galaxies in the framework of hierarchical galaxy formation. We extend a semi-analytic model of galaxy formation that reproduces a number of observations for local and high-z galaxies, without changing the original model parameters but introducing a physically-motivated modelling to describe the es… ▽ More We present a new theoretical model for the luminosity functions (LFs) of Lyman alpha (Lya) emitting galaxies in the framework of hierarchical galaxy formation. We extend a semi-analytic model of galaxy formation that reproduces a number of observations for local and high-z galaxies, without changing the original model parameters but introducing a physically-motivated modelling to describe the escape fraction of Lya photons from host galaxies (f_esc). Though a previous study using a hierarchical clustering model simply assumed a constant and universal value of f_esc, we incorporate two new effects on f_esc: extinction by interstellar dust and galaxy-scale outflow induced as a star formation feedback. It is found that the new model nicely reproduces all the observed Lya LFs of the Lya emitters (LAEs) at different redshifts in z ~ 3-6. Especially, the rather surprisingly small evolution of the observed LAE Lya LFs compared with the dark halo mass function is naturally reproduced. Our model predicts that galaxies with strong outflows and f_esc ~ 1 are dominant in the observed LFs. This is also consistent with available observations, while the simple universal f_esc model requires f_esc << 1 not to overproduce the brightest LAEs. On the other hand, we found that our model significantly overpredicts LAEs at z > 6, and absorption of Lya photons by neutral hydrogen in intergalactic medium (IGM) is a reasonable interpretation for the discrepancy. This indicates that the IGM neutral fraction x_HI rapidly evolves from x_HI << 1 at z < 6 to a value of order unity at z ~ 6-7, which is broadly consistent with other observational constraints on the reionization history. △ Less

Submitted 7 August, 2007; v1 submitted 30 May, 2007; originally announced May 2007.

Comments: 14 pages, 7 figures, 1 table; accepted to ApJ; the html abstract is replaced to match the accepted version, the .ps and .pdf files are strictly identical between the 2nd and the 3rd versions

The Effect of Orbital Eccentricity on Gravitational Wave Background Radiation from Supermassive Black Hole Binaries

Authors: Motohiro Enoki, Masahiro Nagashima

Abstract: A compact binary in an eccentric orbit radiates gravitational waves (GWs) at all integer harmonics of its orbital frequency. In this study, we investigate the effect of orbital eccentricity on the expected gravitational background radiation (GWBR) from supermassive black hole (SMBH) binaries in the nuclei of galaxies. For this purpose, we formulate a power spectrum of the GWBR from cosmological… ▽ More A compact binary in an eccentric orbit radiates gravitational waves (GWs) at all integer harmonics of its orbital frequency. In this study, we investigate the effect of orbital eccentricity on the expected gravitational background radiation (GWBR) from supermassive black hole (SMBH) binaries in the nuclei of galaxies. For this purpose, we formulate a power spectrum of the GWBR from cosmological evolving eccentric binaries. Then, we apply this formulation to the case of the GWBR from SMBH binaries. The key to doing this is to correctly estimate the number density of coalescing SMBH binaries. In this study, we use a semi-analytic model of galaxy and SMBH formation. We find that the power spectrum of the GWBR from SMBH binaries on eccentric orbits is suppressed for frequencies $f \lesssim 1~{\rm nHz}$ if the initial eccentricity, $e_0$, satisfies $e_0 > 0.2$ and the initial semi-major axis is 300 times Scwarzschild radius. Our model predicts that while the overall shape and amplitude of the power spectrum depend strongly on the processes of galaxy formation, the eccentricity of binaries can affect the shape of the power spectrum for lower frequencies, i.e., $f \lesssim 1~{\rm nHz}$. Pulsar timing measurements, which can detect GW in this frequency range, could constrain the effect of eccentricity on the power spectrum of the GWBR from SMBH binaries. △ Less

Submitted 23 February, 2007; v1 submitted 14 September, 2006; originally announced September 2006.

Comments: 16 pages, 12 figures, typos correced, matched to the published version

Journal ref: Prog.Theor.Phys. 117 (2007) 241-256

How Long Can Tiny HI Clouds Survive?

Authors: Masahiro Nagashima, Shu-ichiro Inutsuka, Hiroshi Koyama

Abstract: We estimate the evaporation timescale for spherical HI clouds consisting of the cold neutral medium surrounded by the warm neutral medium. We focus on clouds smaller than 1pc, which corresponds to tiny HI clouds recently discovered by Braun & Kanekar and Stanimirovi{ć} & Heiles. By performing one-dimensional spherically symmetric numerical simulations of the two-phase interstellar medium (ISM),… ▽ More We estimate the evaporation timescale for spherical HI clouds consisting of the cold neutral medium surrounded by the warm neutral medium. We focus on clouds smaller than 1pc, which corresponds to tiny HI clouds recently discovered by Braun & Kanekar and Stanimirovi{ć} & Heiles. By performing one-dimensional spherically symmetric numerical simulations of the two-phase interstellar medium (ISM), we derive the timescales as a function of the cloud size and of pressure of the ambient warm medium. We find that the evaporation timescale of the clouds of 0.01 pc is about 1Myr with standard ISM pressure, $p/k_{B}\sim 10^{3.5}$ K cm$^{-3}$, and for clouds larger than about 0.1 pc it depends strongly on the pressure. In high pressure cases, there exists a critical radius for clouds growing as a function of pressure, but the minimum critical size is $\sim$ 0.03 pc for a standard environment. If tiny HI clouds exist ubiquitously, our analysis suggests two implications: tiny HI clouds are formed continuously with the timescale of 1Myr, or the ambient pressure around the clouds is much higher than the standard ISM pressure. We also find that the results agree well with those obtained by assuming quasi-steady state evolution. The cloud-size dependence of the timescale is well explained by an analytic approximate formula derived by Nagashima, Koyama & Inutsuka. We also compare it with the evaporation rate given by McKee & Cowie. △ Less

Submitted 10 March, 2006; originally announced March 2006.

Comments: 4 pages, 3 figures

Journal ref: Astrophys.J.652:L41-L44,2006

Clustering of Lyman Break Galaxies at z=4 and 5 in The Subaru Deep Field: Luminosity Dependence of The Correlation Function Slope

Authors: Nobunari Kashikawa, Makiko Yoshida, Kazuhiro Shimasaku, Masahiro Nagashima, Hideki Yahagi, Masami Ouchi, Yuichi Matsuda, Matthew A. Malkan, Mamoru Doi, Masanori Iye, the SDF team

Abstract: We explored the clustering properties of Lyman Break Galaxies (LBGs) at z=4 and 5 with an angular two-point correlation function on the basis of the very deep and wide Subaru Deep Field data. We found an apparent dependence of the correlation function slope on UV luminosity for LBGs at both z=4 and 5. More luminous LBGs have a steeper correlation function. To compare these observational results,… ▽ More We explored the clustering properties of Lyman Break Galaxies (LBGs) at z=4 and 5 with an angular two-point correlation function on the basis of the very deep and wide Subaru Deep Field data. We found an apparent dependence of the correlation function slope on UV luminosity for LBGs at both z=4 and 5. More luminous LBGs have a steeper correlation function. To compare these observational results, we constructed numerical mock LBG catalogs based on a semianalytic model of hierarchical clustering combined with high-resolution N-body simulation, carefully mimicking the observational selection effects. The luminosity functions for LBGs predicted by this mock catalog were found to be almost consistent with the observation. Moreover, the overall correlation functions of LBGs were reproduced reasonably well. The observed dependence of the clustering on UV luminosity was not reproduced by the model, unless subsamples of distinct halo mass were considered. That is, LBGs belonging to more massive dark haloes had steeper and larger-amplitude correlation functions. With this model, we found that LBG multiplicity in massive dark halos amplifies the clustering strength at small scales, which steepens the slope of the correlation function. The hierarchical clustering model could therefore be reconciled with the observed luminosity-dependence of the angular correlation function, if there is a tight correlation between UV luminosity and halo mass. Our finding that the slope of the correlation function depends on luminosity could be an indication that massive dark halos hosted multiple bright LBGs (abridged). △ Less

Submitted 13 January, 2006; v1 submitted 19 September, 2005; originally announced September 2005.

Comments: 16 pages, 17 figures, Accepted for publication in ApJ, Full resolution version is available at http://zone.mtk.nao.ac.jp/~kashik/sdf/acf/sdf_lbgacf.pdf

Journal ref: Astrophys.J.637:631-647,2006

Numerical Galaxy Catalog -I. A Semi-analytic Model of Galaxy Formation with N-body simulations

Authors: Masahiro Nagashima, Hideki Yahagi, Motohiro Enoki, Yuzuru Yoshii, Naoteru Gouda

Abstract: We construct the Numerical Galaxy Catalog ($ν$GC), based on a semi-analytic model of galaxy formation combined with high-resolution N-body simulations in a $Λ$-dominated flat cold dark matter ($Λ$CDM) cosmological model. The model includes several essential ingredients for galaxy formation, such as merging histories of dark halos directly taken from N-body simulations, radiative gas cooling, sta… ▽ More We construct the Numerical Galaxy Catalog ($ν$GC), based on a semi-analytic model of galaxy formation combined with high-resolution N-body simulations in a $Λ$-dominated flat cold dark matter ($Λ$CDM) cosmological model. The model includes several essential ingredients for galaxy formation, such as merging histories of dark halos directly taken from N-body simulations, radiative gas cooling, star formation, heating by supernova explosions (supernova feedback), mergers of galaxies, population synthesis, and extinction by internal dust and intervening HI clouds. As the first paper in a series using this model, we focus on basic photometric, structural and kinematical properties of galaxies at present and high redshifts. Two sets of model parameters are examined, strong and weak supernova feedback models, which are in good agreement with observational luminosity functions of local galaxies in a range of observational uncertainty. Both models agree well with many observations such as cold gas mass-to-stellar luminosity ratios of spiral galaxies, HI mass functions, galaxy sizes, faint galaxy number counts and photometric redshift distributions in optical pass-bands, isophotal angular sizes, and cosmic star formation rates. In particular, the strong supernova feedback model is in much better agreement with near-infrared (K'-band) faint galaxy number counts and redshift distribution than the weak feedback model and our previous semi-analytic models based on the extended Press-Schechter formalism. (Abridged) △ Less

Submitted 2 August, 2005; originally announced August 2005.

Comments: 26 pages including 27 figures, accepted for publication in ApJ, full-resolution version is available at http://grape.astron.s.u-tokyo.ac.jp/~yahagi/nugc/

Journal ref: Astrophys.J.634:26-50,2005

The metal enrichment of elliptical galaxies in hierarchical galaxy formation models

Authors: Masahiro Nagashima, Cedric G. Lacey, Takashi Okamoto, Carlton M. Baugh, Carlos S. Frenk, Shaun Cole

Abstract: We investigate the metal enrichment of elliptical galaxies in the framework of hierarchical models of galaxy formation. The semi-analytical model we use in this Letter, which has been used to study the metal enrichment of the intracluster medium (ICM) by Nagashima et al., includes the effects of flows of gas and metals both into and out of galaxies and the processes of metal enrichment due to bo… ▽ More We investigate the metal enrichment of elliptical galaxies in the framework of hierarchical models of galaxy formation. The semi-analytical model we use in this Letter, which has been used to study the metal enrichment of the intracluster medium (ICM) by Nagashima et al., includes the effects of flows of gas and metals both into and out of galaxies and the processes of metal enrichment due to both type Ia and type II supernovae. We adopt a solar neighbourhood IMF for star formation in discs, but consider models in which starbursts have either a solar neighbourhood IMF or a top-heavy IMF. We find that the alpha-element abundance in ellipticals is consistent with observed values only if the top-heavy IMF is used. This result is consistent with our previous study on the metal enrichment of the ICM. We also discuss the abundance ratio of alpha elements to iron as a function of velocity dispersion and metallicity. We find that models with a top-heavy IMF match the alpha/Fe ratios observed in typical L_* ellipticals, but none of the models reproduce the observed increase of alpha/Fe with velocity dispersion. △ Less

Submitted 14 July, 2005; v1 submitted 27 April, 2005; originally announced April 2005.

Comments: 5 pages, 2 figures, to appear in MNRAS Letters

Gamma-Ray Background from Neutralino Annihilation in the First Cosmological Objects

Authors: Takeshi Oda, Tomonori Totani, Masahiro Nagashima

Abstract: The paradigm of the neutralino dark matter predicts that the first gravitationally bound objects are earth-mass sized microhaloes, which would emit annihilation gamma-rays. Here we show that, though the flux from individual nearest microhaloes is extremely difficult to detect, meaningful constraints on their survival probability and internal density profile can be set by requiring that the galac… ▽ More The paradigm of the neutralino dark matter predicts that the first gravitationally bound objects are earth-mass sized microhaloes, which would emit annihilation gamma-rays. Here we show that, though the flux from individual nearest microhaloes is extremely difficult to detect, meaningful constraints on their survival probability and internal density profile can be set by requiring that the galactic and extragalactic gamma-ray background flux from the microhaloes does not exceed the existing EGRET background data. Possible disruption of microhaloes by stellar encounters does not significantly reduce the background flux. If the probability for microhaloes to survive the hierarchical clustering process of dark matter is as large as indicated by a recent simulation, they could be a significant component of the observed background flux in some photon energy range, even with the standard annihilation cross section and conservative internal density profile of microhaloes. The integrated gamma-ray flux from microhaloes in the halo of the Andromeda galaxy may also be detectable by observations in the near future. △ Less

Submitted 7 September, 2005; v1 submitted 5 April, 2005; originally announced April 2005.

Comments: 5 pages including 3 figures. The version accepted for publication in ApJL. Presentation improved, main conclusions unchanged

Journal ref: Astrophys.J. 633 (2005) L65-L68

Evaporation and Condensation of HI clouds in thermally bistable interstellar media: semi-analytic description of isobaric dynamics of curved interfaces

Authors: Masahiro Nagashima, Hiroshi Koyama, Shu-ichiro Inutsuka

Abstract: We analyze the evaporation and condensation of spherical and cylindrical HI clouds of the cold neutral medium surrounded by the warm neutral medium. Because the interstellar medium including those two phases is well described as a thermally bistable fluid, it is useful to apply pattern formation theories to the dynamics of the interface between the two phases. Assuming isobaric evolution of flui… ▽ More We analyze the evaporation and condensation of spherical and cylindrical HI clouds of the cold neutral medium surrounded by the warm neutral medium. Because the interstellar medium including those two phases is well described as a thermally bistable fluid, it is useful to apply pattern formation theories to the dynamics of the interface between the two phases. Assuming isobaric evolution of fluids and a simple cubic form of the heat-loss function, we show the curvature effects of the interface. We find that approximate solutions for spherical clouds are in good agreement with numerically obtained solutions. We extend our analysis to general curved fronts taking into account the curvature effects explicitly. We find that the curvature effects always stabilise curved interfaces under assumptions such as isobaric evolution we adopt in this Letter. △ Less

Submitted 27 April, 2005; v1 submitted 7 March, 2005; originally announced March 2005.

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

Gravitational Radiation from Coalescing Supermassive Black Hole Binaries in a Hierarchical Galaxy Formation Model

Authors: Motohiro Enoki, Kaiki T. Inoue, Masahiro Nagashima, Naoshi Sugiyama

Abstract: We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts r… ▽ More We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts rates owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 μ{\rm Hz})^{-2/3}$ for $f \lesssim 1 μ{\rm Hz}$. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$. △ Less

Submitted 28 February, 2005; v1 submitted 25 February, 2005; originally announced February 2005.

Comments: 6 pages, 5 figures, To appear in Proceedings of the 22nd Texas Symposium on Relativistic Astrophysics, Stanford University, December 2004 (typos corrected)

Gamma-Rays from Large Scale Structure Formation and the Warm-Hot Intergalactic Medium: Cosmic Baryometry with Gamma-Rays

Authors: Susumu Inoue, Masahiro Nagashima

Abstract: It is shown that inverse Compton gamma-rays from electrons accelerated in large scale structure formation shocks can be crucially affected by non-gravitational effects such as radiative cooling and galaxy formation, with corresponding uncertainties by an order of magnitude in either the gamma-ray source counts or the extragalactic background contribution. However, this also implies that such gam… ▽ More It is shown that inverse Compton gamma-rays from electrons accelerated in large scale structure formation shocks can be crucially affected by non-gravitational effects such as radiative cooling and galaxy formation, with corresponding uncertainties by an order of magnitude in either the gamma-ray source counts or the extragalactic background contribution. However, this also implies that such gamma-rays may in the near future provide us with valuable information about the fraction of cosmic baryons in different forms, particularly the warm-hot intergalactic medium where the majority of the baryons in the universe are believed to reside. We address this problem in a simple way through semi-analytic modeling of structure formation shocks which self-consistently treats merger and accretion shocks. △ Less

Submitted 17 February, 2005; originally announced February 2005.

Comments: 7 pages, 2 figures; slightly modified version of article to appear in Proc. of the "International Symposium on High Energy Gamma-Ray Astronomy" Heidelberg (2004), eds. F.A. Aharonian, H.J. Voelk and D. Horns (AIP, NY)

Damped Lyman alpha Absorbing Galaxies At Low Redshifts z<1 From Hierarchical Galaxy Formation Models

Authors: Katsuya Okoshi, Masahiro Nagashima

Abstract: We investigate Damped Ly-alpha absorbing galaxies (DLA galaxies) at low redshifts z<1 in the hierarchical structure formation scenario to clarify the nature of DLA galaxies because observational data of such galaxies mainly at low redshifts are currently available. We find that our model well reproduces distributions of fundamental properties of DLA galaxies such as luminosities, column densitie… ▽ More We investigate Damped Ly-alpha absorbing galaxies (DLA galaxies) at low redshifts z<1 in the hierarchical structure formation scenario to clarify the nature of DLA galaxies because observational data of such galaxies mainly at low redshifts are currently available. We find that our model well reproduces distributions of fundamental properties of DLA galaxies such as luminosities, column densities, impact parameters obtained by optical and near-infrared imagings. Our results suggest that DLA systems primarily consist of low luminosity galaxies with small impact parameters (typical radius about 3 kpc, surface brightness from 22 to 27 mag arcsec^{-2}) similar to low surface brightness (LSB) galaxies. In addition, we investigate selection biases arising from the faintness and from the masking effect which prevents us from identifying a DLA galaxy hidden or contaminated by a point spread function of a background quasar. We find that the latter affects the distributions of DLA properties more seriously rather than the former, and that the observational data are well reproduced only when taking into account the masking effect. The missing rate of DLA galaxies by the masking effect attains 60-90 % in the sample at redshift 0<z<1 when an angular size limit is as small as 1 arcsec. Furthermore we find a tight correlation between HI mass and cross section of DLA galaxies, and also find that HI-rich galaxies with M(HI) \sim 10^{9} M_sun dominate DLA systems. These features are entirely consistent with those from the Arecibo Dual-Beam Survey which is a blind 21 cm survey. Finally we discuss star formation rates, and find that they are typically about 10^{-2} M_sun yr^{-1} as low as those in LSB galaxies. △ Less

Submitted 21 December, 2004; originally announced December 2004.

Comments: 21 pages, 13 figures, Accepted for publication in Astrophsical Journal

Journal ref: Astrophys.J. 623 (2005) 99-111

The metal enrichment of the intracluster medium in hierarchical galaxy formation models

Authors: Masahiro Nagashima, Cedric G. Lacey, Carlton M. Baugh, Carlos S. Frenk, Shaun Cole

Abstract: We investigate the metal enrichment of the intracluster medium (ICM) in the framework of hierarchical models of galaxy formation. We calculate the formation and evolution of galaxies and clusters using a semi-analytical model which includes the effects of flows of gas and metals both into and out of galaxies. For the first time in a semi-analytical model, we calculate the production of both alph… ▽ More We investigate the metal enrichment of the intracluster medium (ICM) in the framework of hierarchical models of galaxy formation. We calculate the formation and evolution of galaxies and clusters using a semi-analytical model which includes the effects of flows of gas and metals both into and out of galaxies. For the first time in a semi-analytical model, we calculate the production of both alpha and iron-peak elements based on theoretical models for the lifetimes and ejecta of type Ia and type II supernovae (SNe Ia and SNe II). It is essential to include the long lifetimes of the SNIa progenitors in order to correctly model the evolution of the iron-peak elements. We find that if all stars form with an IMF similar to that found in the solar neighbourhood, then the metallicities of O, Mg, Si and Fe in the ICM are predicted to be 2-3 times lower than observed values. In contrast, a model (also favoured on other grounds) in which stars formed in bursts triggered by galaxy mergers have a top-heavy IMF reproduces the observed ICM abundances of O, Mg, Si and Fe. The same model predicts ratios of ICM mass to total stellar luminosity in clusters which agree well with observations. According to our model, the bulk of the metals in clusters are produced by L* and brighter galaxies. [abridged] △ Less

Submitted 23 December, 2004; v1 submitted 27 August, 2004; originally announced August 2004.

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

Journal ref: Mon.Not.Roy.Astron.Soc. 358 (2005) 1247-1266

Type Ia supernovae in a hierarchical galaxy formation model: the Milky Way

Authors: Masahiro Nagashima, Takashi Okamoto

Abstract: We investigate chemical evolution in Milky Way-like galaxies based on the cold dark matter model in which cosmic structures form via hierarchical merging. We introduce chemical enrichment due to type Ia supernovae (SNe Ia) into the Mitaka semi-analytic galaxy formation model developed by Nagashima & Yoshii. For the first time we derive distributions of stellar metallicities and their ratios in M… ▽ More We investigate chemical evolution in Milky Way-like galaxies based on the cold dark matter model in which cosmic structures form via hierarchical merging. We introduce chemical enrichment due to type Ia supernovae (SNe Ia) into the Mitaka semi-analytic galaxy formation model developed by Nagashima & Yoshii. For the first time we derive distributions of stellar metallicities and their ratios in Milky Way-like galaxies treating chemical enrichment due to SNe Ia in a hierarchical galaxy formation model self-consistently. As a first attempt, we assume all SNe Ia to have the same lifetime, and assume instantaneous recycling for type II supernovae (SNe II). We find that our model reproduces well the metal abundance ratio [O/Fe] against [Fe/H] and the {iron metallicity distribution function} in the solar neighborhood. This means that the so-called G-dwarf problem is resolved by the hierarchical formation of galaxies, and a gas infall term introduced in traditional monolithic collapse models to solve this problem is well explained by the mixture of some physical processes such as hierarchical merging of dark halos, gas cooling, energy feedback and injection of gas and metals into hot gas due to supernovae. Our model predicts more oxygen-enhanced stars in bulges at [Fe/H] $\simeq 0$ than in disks. This trend seems to be supported by recent observations while they have still uncertainties. More data in number and accuracy will provide independent and important constraints on galaxy formation. (abridged) △ Less

Submitted 3 February, 2006; v1 submitted 25 April, 2004; originally announced April 2004.

Comments: 14 pages including 18 figures, to appear in ApJ

Journal ref: Astrophys.J. 643 (2006) 863-880

Hierarchical Formation of Galaxies with Dynamical Response to Supernova-Induced Gas removal

Authors: Masahiro Nagashima, Yuzuru Yoshii

Abstract: We reanalyze the formation and evolution of galaxies in the hierarchical clustering scenario. Using a semi-analytic model (SAM) of galaxy formation described in this paper, which we hereafter call the Mitaka model, we extensively investigate the observed scaling relations of galaxies among photometric, kinematic, structural and chemical characteristics. In such a scenario, spheroidal galaxies ar… ▽ More We reanalyze the formation and evolution of galaxies in the hierarchical clustering scenario. Using a semi-analytic model (SAM) of galaxy formation described in this paper, which we hereafter call the Mitaka model, we extensively investigate the observed scaling relations of galaxies among photometric, kinematic, structural and chemical characteristics. In such a scenario, spheroidal galaxies are assumed to be formed by major merger and subsequent starburst, in contrast to the traditional scenario of monolithic cloud collapse. As a new ingredient of SAMs, we introduce the effects of dynamical response to supernova-induced gas removal on size and velocity dispersion, which play an important role on dwarf galaxy formation. In previous theoretical studies of dwarf galaxies based on the monolithic cloud collapse given by Yoshii & Arimoto and Dekel & Silk, the dynamical response was treated in the extremes of a purely baryonic cloud and a baryonic cloud fully supported by surrounding dark matter. To improve this simple treatment, in our previous paper, we formulated the dynamical response in more realistic, intermediate situations between the above extremes. While the effects of dynamical response depend on the mass fraction of removed gas from a galaxy, how much amount of the gas remains just after major merger depends on the star formation history. A variety of star formation histories are generated through the Monte Carlo realization of merging histories of dark halos, and it is found that our SAM naturally makes a wide variety of dwarf galaxies and their dispersed characteristics as observed. (Abridged) △ Less

Submitted 25 April, 2004; originally announced April 2004.

Comments: 24 pages including 29 figures, using emulateapj.cls; accepted for publication in ApJ

Journal ref: Astrophys.J. 610 (2004) 23-44

Gravitational Waves from Supermassive Black Hole Coalescence in a Hierarchical Galaxy Formation Model

Authors: Motohiro Enoki, Kaiki T. Inoue, Masahiro Nagashima, Naoshi Sugiyama

Abstract: We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. When galaxies merge, the SMBHs in the host galaxies sink to the center of the new merged galaxy and form a binary system. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estim… ▽ More We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. When galaxies merge, the SMBHs in the host galaxies sink to the center of the new merged galaxy and form a binary system. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 μ{\rm Hz})^{-2/3}$ for $f \lesssim 1 μ{\rm Hz}$ just below the detection limit from measurements of the pulsar timing provided that SMBHs coalesce simultaneously when host galaxies merge. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$. △ Less

Submitted 15 July, 2004; v1 submitted 20 April, 2004; originally announced April 2004.

Comments: Accepted for publication in ApJ. 11 pages, 7 figures

Journal ref: Astrophys.J. 615 (2004) 19

The roles of ram-pressure stripping and minor mergers in evolution of galaxies

Authors: Takashi Okamoto, Masahiro Nagashima

Abstract: We investigate environmental effects on evolution of bright cluster galaxies in a Lambda-dominated cold dark matter universe using a combination of dissipationless N-body simulations and a semi-analytic galaxy formation model. We incorporate effects of ram-pressure stripping (RPS) and minor merger-induced small starburst (minor burst) into our model. By considering minor burst, observed morpholo… ▽ More We investigate environmental effects on evolution of bright cluster galaxies in a Lambda-dominated cold dark matter universe using a combination of dissipationless N-body simulations and a semi-analytic galaxy formation model. We incorporate effects of ram-pressure stripping (RPS) and minor merger-induced small starburst (minor burst) into our model. By considering minor burst, observed morphology-radius relation is successfully reproduced. When we do not consider minor burst, the RPS hardly increases the intermediate B/T population. In addition, the RPS and minor burst are not important for colours or star formation rates of galaxies in the cluster core if star formation time-scale is properly chosen, because the star formation is sufficiently suppressed by consumption of the cold gas. We also find that SF in bulge-dominated galaxies is mainly terminated by starburst induced by major mergers in all environments. △ Less

Submitted 15 April, 2004; originally announced April 2004.

Comments: To appear in the proceedings of IAU colloq. No. 195, "Outskirts of Galaxy Clusters: intense life in the suburbs", Torino, 12-16 March 2004, 5 pages, 2 figures, uses IAU macro

Mass Function of Low Mass Dark Halos

Authors: Hideki Yahagi, Masahiro Nagashima, Yuzuru Yoshii

Abstract: The mass function of dark halos in a Lambda-dominated cold dark matter (LambdaCDM) universe is investigated. 529 output files from five runs of N-body simulations are analyzed using the friends-of-friends cluster finding algorithm. All the runs use 512^3 particles in the box size of 35 h^{-1}Mpc to 140 h^{-1}Mpc. Mass of particles for 35 h^{-1} Mpc runs is 2.67 times 10^7 h^{-1} M_{solar}. Becau… ▽ More The mass function of dark halos in a Lambda-dominated cold dark matter (LambdaCDM) universe is investigated. 529 output files from five runs of N-body simulations are analyzed using the friends-of-friends cluster finding algorithm. All the runs use 512^3 particles in the box size of 35 h^{-1}Mpc to 140 h^{-1}Mpc. Mass of particles for 35 h^{-1} Mpc runs is 2.67 times 10^7 h^{-1} M_{solar}. Because of the high mass resolution of our simulations, the multiplicity function in the low-mass range, where the mass is well below the characteristic mass and $ν= δ_c / σ\lesssim 1.0$, is evaluated in the present work, and is well fitted by the functional form proposed by Sheth & Tormen (ST). However, the maximum value of the multiplicity function from our simulations at $ν\sim 1$ is smaller, and its low mass tail is shallower when compared with the ST multiplicity function. △ Less

Submitted 7 January, 2004; originally announced January 2004.

Comments: to be published in ApJ

Journal ref: Astrophys.J. 605 (2004) 709-713

Chemical Enrichment in Damped Lyman Alpha Systems From Hierarchical Galaxy Formation Models

Authors: Katsuya Okoshi, Masahiro Nagashima, Naoteru Gouda, Satoshi Yoshioka

Abstract: We investigate chemical enrichment in Damped Lyman alpha (DLA) systems in the hierarchical structure formation scenario using a semi-analytic model of galaxy formation. The model developed by Nagashima, Totani, Gouda and Yoshii takes into account various selection effects on high-redshift galaxies and can show fundamental observational properties of galaxies, such as luminosity functions and num… ▽ More We investigate chemical enrichment in Damped Lyman alpha (DLA) systems in the hierarchical structure formation scenario using a semi-analytic model of galaxy formation. The model developed by Nagashima, Totani, Gouda and Yoshii takes into account various selection effects on high-redshift galaxies and can show fundamental observational properties of galaxies, such as luminosity functions and number-magnitude/redshift relations. DLA systems offer the possibilities of measuring metal abundance more accurately than faint galaxies. For example, recent measurements of zinc abundance can provide good evidence for understanding the processes of metal pollution and star formation in DLA systems because zinc is virtually unaffected by dust depletion. Here we focus on this advantage for observation in order to explore the metallicity evolution in DLA systems at high redshifts. We can consistently show the metallicity evolution for reasonable models which also reproduce fundamental properties of local galaxy population. This result suggests that the chemical evolution of DLA systems can be consistently reconciled with the observational features of typical galaxies. We also investigate other properties of DLA systems (column density distribution and mass density of cold gas), and find that star formation in massive galaxies should be more active than that in low-mass ones. This is consistent with the results by Nagashima et al. and Cole et al. in which the star formation timescale is set by reproducing cold gas mass fraction in local spiral galaxies. Finally we discuss host galaxies associated with DLA systems. We conclude that they primarily consist of sub-L* and/or dwarf galaxies from the observations. △ Less

Submitted 25 November, 2003; originally announced November 2003.

Comments: 16 pages, 11 figures, Accepted for publication in Astrophsical Journal

Journal ref: Astrophys.J. 603 (2004) 12-22