Search

In their recent study, Labb\'e et al. used multi-band infrared images captured by the James Webb Space Telescope (JWST) to discover a population of red massive galaxies that formed approximately 600 million years after the Big Bang. The authors reported an extraordinarily large density of these galaxies, with stellar masses exceeding $10^{10}$ solar masses, which, if confirmed, challenges the standard cosmological model as suggested by recent studies. However, this conclusion is disputed. We contend that during the early epochs of the universe the stellar mass-to-light ratio could not have reached the values reported by Labb\'e et al. A model of galaxy formation based on standard cosmology provides support for this hypothesis, predicting the formation of massive galaxies with higher ultraviolet (UV) luminosity, which produce several hundred solar masses of stars per year and containing significant dust. These forecasts are consistent with the abundance of JWST/HST galaxies selected photometrically in the rest-frame UV wavelengths and with the properties of the recent spectroscopically-confirmed JWST/HST galaxies formed during that era. Discrepancies with Labb\'e et al. may arise from overestimation of the stellar masses, systematic uncertainties, absence of JWST/MIRI data, heavy dust extinction affecting UV luminosities, or misidentification of faint red AGN galaxies at closer redshifts. The current JWST/HST results, combined with a realistic galaxy formation model, provide strong confirmation of the standard cosmology., Comment: Submitted to Nature, matters arising

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$, that has sufficient mass resolution to resolve dwarf-size systems, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$\nu^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$\nu^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, in particular for small survey areas (0.01 and 0.1 deg$^2$). 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 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). In particular, we predict that the Euclid deep survey will find 120--240 (16--80) AGN -- depending on the model -- with rest-frame UV absolute magnitude brighter than $-20$ ($-20.5$) at $z=6.3$ ($z=7$) in the Euclid H-band deep survey (abridged)., Comment: 17 pages, 17 figures, submitted to MNRAS, corrected the following link, All Uchuu-$\nu^2$GC data used for this paper is publicly available in our Skies & Universes website http://skiesanduniverses.org/Simulations/Uchuu/

The cluster mass-richness relation (MRR) is an observationally efficient and potentially powerful cosmological tool for constraining the mean matter density of the universe and the amplitude of fluctuations using the cluster abundance technique. We derive the MRR relation using GalWCat19, a publicly available galaxy cluster catalog we created from the Sloan Digital Sky Survey-DR13 spectroscopic dataset. The MRR shows a tail at the low-richness end. Using the Illustris-TNG and mini-Uchuu cosmological numerical simulations, we demonstrate that this tail is caused by systematical uncertainties. We show that, by means of a judicious cut, identified by the use of the Hinge function, it is possible to determine a richness threshold above which the MRR is linear i.e., where cluster mass scales with richness as logM_200 = alpha + beta logN_200. We derive the MRR and show it is consistent with both sets of simulations with a slope of beta ~ 1. We use our MRR to estimate cluster masses from the GalWCat19 catalog which we then use to set constraints on omega_m and sigma_8. Utilizing the all-member MRR, we obtain constraints of omega_m = 0.31 (+0.04-0.03) and sigma_8 = 0.82 (+0.05-0.04), and utilizing the red-member MRR, we obtain omega_m = 0.31 (+0.04-0.03) and sigma_8 = 0.81 (+0.05-0.04). Our constraints on omega_m and sigma_8 are consistent and very competitive with the Planck 2018 results.

Elucidating dark matter density profiles in the Galactic dwarf satellites is essential to understanding not only the quintessence of dark matter, but also the evolution of the satellites themselves. In this work, we present the current constraints on dark matter densities in the Galactic ultra-faint dwarf (UFD) and diffuse galaxies. Applying our constructed non-spherical mass models to the currently available kinematic data of the 25 UFDs and 2 diffuse satellites, we find that whereas most of the galaxies have huge uncertainties on the inferred dark matter density profiles, Eridanus~II, Segue~I, and Willman~1 favor cuspy central profiles even considering effects of a prior bias. We compare our results with the simulated subhalos on the plane between the dark matter density at 150~pc and the pericenter distance. We find that the most observed satellites and the simulated subhalos are similarly distributed on this plane, except for Antlia~2, Crater~2, and Tucana~3, which are less than one tenth of the density. Despite considerable tidal effects, the subhalos detected by commonly-used subhalo finders have difficulty in explaining such a huge deviation. We also estimate the dynamical mass-to-light ratios of the satellites and confirm the ratio is linked to stellar mass and metallicity. Tucana~3 deviates largely from these relations, while it follows the mass-metallicity relation. This indicates that Tucana~3 has a cored dark matter halo, despite a significant uncertainty in its ratios., 19 pages, 13 figures, accepted for publication in ApJ

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., 13 pages, 11 figures, accepted for publication in the ApJ

The survey of Lyman $��$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined $2^3$ models depending on the presence or absence of dispersion of halo Ly$��$ emissivity, dispersion of the halo Ly$��$ optical depth, $��_��$, and halo mass dependence of $��_��$. The unique free parameter in our model, a pivot value of $��_��$, is calibrated so as to reproduce the $z=5.7$ Ly$��$ luminosity function (LF). We compare our model predictions with Ly$��$ LFs at $z=6.6$ and $7.3$, LAE angular auto-correlation functions (ACFs) at $z=5.7$ and $6.6$, and LAE fractions in Lyman break galaxies at $5, PASJ in press, 30 pages, 20 figures

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., 31 pages and 29 figures. Accepted for publication in MNRAS

Cosmic reionization is known to be a major phase transition of the gas in the Universe. Since astronomical objects formed in the early Universe, such as the first stars, galaxies and black holes, are expected to have caused cosmic reionization, the formation history and properties of such objects are closely related to the reionization process. In spite of the importance of exploring reionization, our understandings regarding reionization is not sufficient yet. Square Kilometre Array (SKA) is a next-generation large telescope that will be operated in the next decade. Although several programs of next-generation telescopes are currently scheduled, the SKA will be the unique telescope with a potential to directly observe neutral hydrogen up to z~30, and provide us with valuable information on the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The early science with the SKA will start in a few years; it is thus the time for us to elaborate a strategy for CD/EoR Science with the SKA. The purpose of this document is to introduce Japanese scientific interests in the SKA project and to report results of our investigation., 14 pages, 12 figures, Japan SKA consortium (SKA-JP) Cosmic Dawn/Epoch of Reionization Science book

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

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., 18 pages, 16 figures, accepted, simulation data can be downloaded on http://hpc.imit.chiba-u.jp/~nngc/

As an entry for the 2012 Gordon-Bell performance prize, we report performance results of astrophysical N-body simulations of one trillion particles performed on the full system of K computer. This is the first gravitational trillion-body simulation in the world. We describe the scientific motivation, the numerical algorithm, the parallelization strategy, and the performance analysis. Unlike many previous Gordon-Bell prize winners that used the tree algorithm for astrophysical N-body simulations, we used the hybrid TreePM method, for similar level of accuracy in which the short-range force is calculated by the tree algorithm, and the long-range force is solved by the particle-mesh algorithm. We developed a highly-tuned gravity kernel for short-range forces, and a novel communication algorithm for long-range forces. The average performance on 24576 and 82944 nodes of K computer are 1.53 and 4.45 Pflops, which correspond to 49% and 42% of the peak speed., Comment: 10 pages, 6 figures, Proceedings of Supercomputing 2012 (http://sc12.supercomputing.org/), Gordon Bell Prize Winner. Additional information is http://www.ccs.tsukuba.ac.jp/CCS/eng/gbp2012

Understanding the universe is hampered by the elusiveness of its most common constituent, cold dark matter. Almost impossible to observe, dark matter can be studied effectively by means of simulation and there is probably no other research field where simulation has led to so much progress in the last decade. Cosmological N-body simulations are an essential tool for evolving density perturbations in the nonlinear regime. Simulating the formation of large-scale structures in the universe, however, is still a challenge due to the enormous dynamic range in spatial and temporal coordinates, and due to the enormous computer resources required. The dynamic range is generally dealt with by the hybridization of numerical techniques. We deal with the computational requirements by connecting two supercomputers via an optical network and make them operate as a single machine. This is challenging, if only for the fact that the supercomputers of our choice are separated by half the planet, as one is located in Amsterdam and the other is in Tokyo. The co-scheduling of the two computers and the 'gridification' of the code enables us to achieve a 90% efficiency for this distributed intercontinental supercomputer., Accepted for publication in IEEE Computer