Search

To date, JWST has detected the earliest known star clusters in our Universe (Adamo et al. 2024, Messa et al. 2024, Vanzella et al. 2024, Mowla et al. 2024). They appear to be relatively compact (~few pc, Adamo et al. 2024) and had only recently formed their stars. It was speculated that these clusters may be the earliest progenitors of globular clusters ever detected. Globular clusters are a relic of the initial stages of star formation in the Universe. However, because they contain little to no dark matter (e.g., Heggie & Hut 1996, Bradford et al. 2011, Conroy et al. 2011, Ibata et al. 2013), their formation mechanism poses a significant theoretical challenge. A recent suggestion pointed out that the relative velocity between the gas and the dark matter (Tseliakhovich & Hirata 2010) in the early Universe could naturally form potentially star-forming regions outside of dark matter halos. Here, for the first time, we follow the star formation process of these early Universe objects using high-resolution hydrodynamical simulations, including mechanical feedback. Our results suggest that the first dark matter-less star clusters are top-heavy, with a higher abundance of massive stars compared to today's clusters and extremely high stellar mass surface densities compared to the local Universe., Comment: 10 pages, 5 figures

The current paradigm for the co-evolution of galaxies and their supermassive black holes postulates that dust-obscured active galactic nuclei (AGNs) represent a transitional phase towards a more luminous and unobscured state. However, our understanding of dusty AGNs and their host galaxies at early cosmic times is inadequate due to observational limitations. Here, we present JWST observations of CID-931, an X-ray-detected AGN at a spectroscopic redshift of $z_{\rm spec}=4.91$. Multiband NIRCam imaging from the COSMOS-Web program reveals an unresolved red core, similar to JWST-discovered dusty AGNs. Strikingly, the red core is surrounded by at least eight massive star-forming clumps spread over $1.\!\!^{\prime\prime}6 \approx 10~{\rm kpc}$, each of which has a stellar mass of $10^9-10^{10}M_\odot$ and $\sim0.1-1~{\rm kpc}$ in radius. The whole system amounts to $10^{11}M_\odot$ in stellar mass, higher than typical star-forming galaxies at the same epoch. In this system, gas inflows and/or complex merger events may trigger clump formation and AGN activity thus leading to the rapid formation of a massive galaxy hosting a supermassive black hole. Future follow-up observations will provide new insights into the evolution of the galaxy-black hole relationship during such transitional phases in the early universe., Comment: 13 pages, 10 figures, submitted to PASJ

Population III stars are characterized by extremely low metallicities as they are thought to be formed from a pristine gas in the early Universe. Although the existence of Population III stars is widely accepted, the lack of direct observational evidence hampers the study of the nature of the putative stars. In this article, we explore the possibilities of constraining the nature of the oldest stars by using the luminosity function of their remnants -- white dwarfs. We study the formation and evolution of white dwarf populations by following star formation in a Milky Way-like galaxy using the semi-analytic model A-SLOTH. We derive the white dwarf luminosity function by applying a linear Initial-Final Mass Relation and Mestel's cooling model. The obtained luminosity function is generally in agreement with available observations and theoretical predictions -- with an exponential increase to a maximum of Mabs = 16 and a sudden drop for Mabs > 16. We explore the uncertainties of our model and compare them to the observational estimates. We adopt two different models of the initial mass function of Population III stars to show that the faint end of the luminosity function imprints the signature of Population III remnants. If the feature is detected in future observations, it would provide a clue to Population III stars and would also be an indirect evidence of low- to intermediate-mass Population III stars. We discuss the challenges and prospects for detecting the signatures.

We present an analytic description of the spherically symmetric gravitational collapse of radiatively cooling gas clouds. The approach is based on developing the "one-zone" density-temperature relationship of the gas into a full dynamical model. We convert this density-temperature relationship into a barotropic equation of state, which we use to calculate the density and velocity profiles of the gas. From these quantities we calculate the time-dependent mass accretion rate onto the center of the cloud. The approach clarifies the mechanism by which radiative cooling induces gravitational instability. In particular, we distinguish the rapid, quasi-equilibrium contraction of a cooling gas core to high central densities from the legitimate instability this contraction establishes in the envelope. We develop a refined criterion for the mass scale of this instability, based only on the chemical-thermal evolution in the core. We explicate our model in the context of a primordial mini-halo cooled by molecular hydrogen, and then provide two further examples, a delayed collapse with hydrogen deuteride cooling and the collapse of an atomic cooling halo. In all three cases, our results agree well with full hydrodynamical treatments., Comment: 11+3 pages, 14+3 figures

We report the highest-redshift detection of [O I] 63 $\mu$m from a luminous quasar, J2054-0005, at $z=6.04$ based on the Atacama Large Millimeter/sub-millimeter Array Band 9 observations. The [O I] 63 $\mu$m line luminosity is $(4.5\pm1.5) \times 10^{9}~L_{\rm \odot}$, corresponding to the [O I] 63 $\mu$m-to-far-infrared luminosity ratio of $\approx 6.7\times10^{-4}$, which is consistent with the value obtained in the local universe. Remarkably, [O I] 63 $\mu$m is as bright as [C II] 158 $\mu$m, resulting in the [O I]-to-[C II] line luminosity ratio of $1.3\pm0.5$. Based on a careful comparison of the luminosity ratios of [O I] 63 $\mu$m, [C II] 158 $\mu$m, and dust continuum emission to models of photo-dissociation regions, we find that J2054-0005 has a gas density log($n_{\rm H}$/cm$^{-3}$)$=3.7\pm0.3$ and an incident far-ultraviolet radiation field of log($G/G_{\rm 0}$)$= 3.0\pm0.1$, showing that [O I] 63 $\mu$m serves as an important coolant of the dense and warm gas in J2054-0005. A close examination of the [O I] and [C II] line profiles suggests that the [O I] line may be partially self-absorbed, however deeper observations are needed to verify this conclusion. Regardless, the gas density and incident radiation field are in a broad agreement with the values obtained in nearby star-forming galaxies and objects with [O I] 63 $\mu$m observations at $z=1-3$ with the Herschel Space Observatory. These results demonstrate the power of ALMA high-frequency observations targeting [O I] 63 $\mu$m to examine the properties of photo-dissociation regions in high-redshift galaxies., Comment: Submitted to PASJ. 5 figures and 3 tables. Comments are very welcome

JWST has brought us new insights into Cosmic Dawn with tentative detection of the unique signatures of metal-free Population III (Pop III) stars, such as strong HeII emission, extremely blue UV spectrum, and enhanced nitrogen abundance. Self-consistent theoretical predictions of the formation rates, sites, and masses of Pop III stars are crucial for interpreting the observations, but are challenging due to complex physical processes operating over the large range of length scales involved. One solution is to combine analytical models for the small-scale star formation process with cosmological simulations that capture the large-scale physics such as structure formation, radiation backgrounds, and baryon-dark matter streaming motion that regulate the conditions of Pop III star formation. We build an analytical model to predict the final masses of Pop III stars/clusters from the properties of star-forming clouds, based on the key results of small-scale star formation simulations and stellar evolution models. Our model for the first time considers the interplay between feedback and fragmentation and covers different modes of Pop III star formation ranging from ordinary small ($\sim 10-2000\ \rm M_\odot$) clusters in molecular-cooling clouds to massive ($\gtrsim 10^{4}\ \rm M_\odot$) clusters containing supermassive ($\sim 10^{4}-3\times 10^{5}\ \rm M_\odot$) stars under violent collapse of atomic-cooling clouds. As an example, the model is applied to the Pop III star-forming clouds in the progenitors of typical haloes hosting high-$z$ luminous quasars, which shows that formation of Pop III massive clusters is common ($\sim 20-70\%$) in such biased ($\sim4\sigma$) regions, and the resulting heavy black hole seeds from supermassive stars can account for a significant fraction of observed luminous ($\gtrsim 10^{46}\ \rm erg\ s^{-1}$) quasars at $z\sim 6$., Comment: 15+8 pages, 9 + 14 figures, accepted for publication in MNRAS

We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/NIRSpec observations have also revealed that very bright galaxy candidates at $z\sim10-13$ identified from ground-based telescope images before JWST are passive galaxies at $z\sim3-4$, emphasizing the necessity of strict screening and spectroscopy in the selection of the brightest galaxies at $z>10$. The UV luminosity functions derived from these spectroscopic results are consistent with a double power-law function, showing tensions with theoretical models at the bright end. To understand the origin of the overabundance of bright galaxies, we investigate their morphologies using JWST/NIRCam high-resolution images obtained in various surveys including PRIMER and COSMOS-Web. We find that $\sim70\%$ of the bright galaxies at $z\sim7$ exhibit clumpy morphologies with multiple sub-components, suggesting merger-induced starburst activity, which is consistent with SED fitting results showing bursty star formation histories. At $z\gtrsim10$, bright galaxies are classified into two types of galaxies; extended ones with weak high-ionization emission lines, and compact ones with strong high-ionization lines including NIV]$\lambda$1486, indicating that at least two different processes (e.g., merger-induced starburst and compact star formation/AGN) are shaping the physical properties of the brightest galaxies at $z\gtrsim10$ and are responsible for their overabundance., Comment: 32 pages, 22 figures, submitted to ApJ

We consider a variant of the stochastic gradient descent (SGD) with a random learning rate and reveal its convergence properties. SGD is a widely used stochastic optimization algorithm in machine learning, especially deep learning. Numerous studies reveal the convergence properties of SGD and its simplified variants. Among these, the analysis of convergence using a stationary distribution of updated parameters provides generalizable results. However, to obtain a stationary distribution, the update direction of the parameters must not degenerate, which limits the applicable variants of SGD. In this study, we consider a novel SGD variant, Poisson SGD, which has degenerated parameter update directions and instead utilizes a random learning rate. Consequently, we demonstrate that a distribution of a parameter updated by Poisson SGD converges to a stationary distribution under weak assumptions on a loss function. Based on this, we further show that Poisson SGD finds global minima in non-convex optimization problems and also evaluate the generalization error using this method. As a proof technique, we approximate the distribution by Poisson SGD with that of the bouncy particle sampler (BPS) and derive its stationary distribution, using the theoretical advance of the piece-wise deterministic Markov process (PDMP)., Comment: 28 pages

At high redshifts ($z\gtrsim12$), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically-induced gas objects (SIGOs) and their associated star clusters outside of dark matter halos. Here, we present a study of the population-level effects of streaming on star formation within both halos and SIGOs in a set of simulations with and without streaming. Notably, we find that streaming actually enhances star formation within individual halos of all masses at redshifts between $z=12$ and $z=20$. This is demonstrated both as an increased star formation rate per object as well as an enhancement of the Kennicutt-Schmidt relation for objects with streaming. We find that our simulations are consistent with some observations at high redshift, but on a population level, they continue to under-predict star formation relative to the majority of observations. Notably, our simulations do not include feedback, and so can be taken as an upper limit on the star formation rate, exacerbating these differences. However, simulations of overdense regions (both with and without streaming) agree with observations, suggesting a strategy for extracting information about the overdensity and streaming velocity in a given survey volume in future observations., Comment: 16 pages, 8 figures

Galaxies that cause the strong gravitational lensing of background galaxies provide us crucial information about the distribution of matter around them. Traditional modelling methods that analyse such strong lenses are both time and resource consuming, require sophisticated lensing codes and modelling expertise. To study the large lens population expected from imaging surveys such as LSST, we need fast and automated analysis methods. In this work, we build and train a simple convolutional neural network with an aim to rapidly predict model parameters of gravitational lenses. We focus on the most important lens mass model parameters, namely, the Einstein radius, the axis ratio and the position angle of the major axis of the mass distribution. The network is trained on a variety of simulated data with an increasing degree of realism and shows satisfactory performance on simulated test data. The trained network is then applied to the real sample of galaxy-scale candidate lenses from the Subaru HSC, a precursor survey to LSST. Unlike the simulated lenses, we do not have the ground truth for the real lenses. Therefore, we have compared our predictions with those from YattaLens, a lens modelling pipeline. Additionally, we also compare the parameter predictions for 10 HSC lenses that were also studied by other conventional modelling methods. These comparisons show a fair quantitative agreement on the Einstein radius, although the axis ratio and the position angle from the network as well as the individual modelling methods, seem to have systematic uncertainties beyond the quoted errors., Comment: 10 pages, 9 figures

Dark matter overdensities around black holes can be searched for by looking at the characteristic imprint they leave on the gravitational waveform of binary black hole mergers. Current theoretical predictions of the density profile of dark matter overdensities are based on highly idealised formation scenarios, in which black holes are assumed to grow adiabatically from an infinitesimal seed to their final mass, compressing dark matter cusps at the center of galactic halos into very dense `spikes'. These scenarios were suitable for dark matter indirect detection studies, since annihilating dark matter cannot reach very high densities, but they fail to capture the dark matter distribution in the innermost regions where the gravitational wave signal is produced. We present here a more realistic formation scenario where black holes form from the collapse of supermassive stars, and follow the evolution of the dark matter density as the supermassive star grows and collapses to a black hole. We show that in this case dark matter forms shallower `mounds', instead of `spikes', on scales comparable with the size of the supermassive stars originating them. We discuss the implications for the detectability of these systems., Comment: 8 pages, 4 figures

Current models of the formation of first galaxies predict low masses and faint objects at extremely high redshifts, z=9-15. However, the first observations of this epoch indicate a higher-than-expected number of bright (sometimes massive) galaxies. Numerical simulations can help to elucidate the mild evolution of the bright end of the UV luminosity function and they can provide the link between the evolution of bright galaxies and variations of the galaxy formation efficiency across different redshifts. We use the FirstLight database of 377 zoom-in cosmological simulations of a mass-complete sample of galaxies. Mock luminosities are estimated by a dust model constrained by current observations of an evolution of the beta-MUV relation at high-z. FirstLight contains a high number of bright galaxies, MUV<-20, consistent with current data at z=6-13. The evolution of the UV cosmic density is driven by the evolution of the galaxy efficiency and the relation between MUV and halo mass. The efficiency of galaxy formation increases significantly with redshift at a fixed halo mass because galactic halos at extremely high redshifts convert gas into stars at a higher rate than at lower redshifts. The high gas densities in galaxies at z>9 enable these high efficiencies. Our simulations predict higher number densities of massive galaxies, Ms=10^9 Msun, than other models with constant efficiency. Cosmological simulations of galaxy formation with self-consistent models of star formation and feedback can reproduce the different regimes of galaxy formation across cosmic history., Comment: 11 pages, 10 figures, published at A&A

We present JWST NIRCam imaging of B14-65666 ("Big Three Dragons"), a bright Lyman-break galaxy system ($M_\text{UV}=-22.5$ mag) at $z=7.15$. The high angular resolution of NIRCam reveals the complex morphology of two galaxy components: galaxy E has a compact core (E-core), surrounded by diffuse, extended, rest-frame optical emission, which is likely to be tidal tails; and galaxy W has a clumpy and elongated morphology with a blue UV slope ($\beta_\text{UV}=-2.2\pm0.1$). The flux excess, F356W$-$F444W, peaks at the E-core ($1.05^{+0.08}_{-0.09}$ mag), tracing the presence of strong [OIII] 4960,5008 \r{A} emission. ALMA archival data show that the bluer galaxy W is brighter in dust continua than the redder galaxy E, while the tails are bright in [OIII] 88 $\mathrm{\mu m}$. The UV/optical and sub-mm SED fitting confirms that B14-65666 is a major merger in a starburst phase as derived from the stellar mass ratio (3:1 to 2:1) and the star-formation rate, $\simeq1$ dex higher than the star-formation main sequence at the same redshift. The galaxy E is a dusty ($A_\text{V}=1.2\pm0.1$ mag) starburst with a possible high dust temperature ($\ge63$-$68$ K). The galaxy W would have a low dust temperature ($\le27$-$33$ K) or patchy stellar-and-dust geometry, as suggested from the infrared excess (IRX) and $\beta_\text{UV}$ diagram. The high optical-to-FIR [OIII] line ratio of the E-core shows its lower gas-phase metallicity ($\simeq0.2$ Z$_{\odot}$) than the galaxy W. These results agree with a scenario where major mergers disturb morphology and induce nuclear dusty starbursts triggered by less-enriched inflows. B14-65666 shows a picture of complex stellar buildup processes during major mergers in the epoch of reionization., Comment: 18 pages, 6 figures, 4 tables. Submitted to ApJ

Recent JWST observations with superb angular resolution have revealed the existence of clumpy galaxies at high redshift through the detection of rest-frame optical emission lines. We use the FirstLight simulation to study the properties of (sub-)galactic clumps that are bright in [OIII] 5007$\mathrm{\mathring{A}}$ line with flux greater than $\sim 10^{-18} \, {\rm erg\, s^{-1}\, cm^{-2}}$, to be detected by JWST. For 62 simulated galaxies that have stellar masses of $(0.5-6) \times 10^{10} \, M_\odot$ at $z=5$, we find clumps in 1828 snapshots in the redshift range $z = 9.5-5.5$. The clumps are identified by the surface density of star formation rate. About one-tenth of the snapshots show the existence of clumpy systems with two or more components. Most of the clumps are formed by mergers and can be characterized by their ages; central clumps dominated by stellar populations older than 50 Myr, and off-centered clumps dominated by younger stellar populations with specific star formation rates of $\sim 50 \, {\rm Gyr^{-1}}$. The latter type of young clumps is formed from gas debris in the tidal tails of major mergers with baryonic mass ratios of $1 \leq q < 4$. The merger-induced clumps are short-lived, and merge within a dynamical time of several tens million years. The number density of the clumpy systems is estimated to be $\sim 10^{-5}\, {\rm cMpc^{-3}}$, which is large enough to be detected in recent JWST surveys., Comment: Version accepted for publication in ApJ. 19 pages, 10 figures, 3 tables, references added

Fuzzy dark matter (FDM) is an intriguing candidate alternative to the standard cold dark matter (CDM). The FDM model predicts that dark halos have characteristic core structures generated by the effect of quantum pressure, which is different from the structure of CDM halos. We devise a semi-analytic model of a FDM halo density profile by assuming that the density distribution results from the redistribution of mass in a halo with the Navarro-Frenk-White profile. We calculate the mass redistribution radius by considering dynamical relaxation within the FDM halo. We adopt a concentration-halo mass relation with lower concentration compared to that in the CDM model below the half mode mass, which originates from the suppressed matter density fluctuations at small length scales. Our model reproduces the core-halo mass relation (CHMR) found in the numerical simulation of \citet{2014NatPh..10..496S} at $z<1$. We show that the CHMR is well described by a double power law, unlike previous studies that approximate it by a single power law. Our model predictions are in reasonable agreement with the results of the largest FDM simulation of \citet{2021MNRAS.506.2603M} at $z=3$. We find that the core mass for a given halo mass follows the log-normal distribution, both in our model and in the simulation results for the first time, and quantitatively compare the variance of the distribution among them. Although our model does not fully explain the scatter of the CHMR, we show the scatter of the concentration-halo mass relation sizably contributes to them., Comment: 22 pages, 17 figures, published in PRD

Continuing work presented in Li et al. (2021), we performed a series of tests to our high-resolution three-dimensional mass map reconstruction algorithm \splinv{}. We test the mass reconstruction accuracy against realistic mock catalogs generated using shear field produced by triaxial halos with the inner density profile of $\rho \propto r^{-1}$ and of $\rho \propto r^{-1.5}$. The galaxy shape noise is modeled based on the Year-1 Subaru Hyper Suprime-Cam (HSC) Survey. After reviewing mathematical details of our algorithm and dark matter halo models, we determine an optimal value of the coefficient of the adaptive LASSO regression penalty term for single halo reconstruction. We successfully measure halo masses for massive triaxial halos; the mass determination accuracy is 5 percent for halos with $M = 10^{14.6}~M_\odot$ at $0.0625\leq z \leq 0.2425$, and 5 percent for those with $10^{14.8}~M_\odot$ at $0.0625\leq z \leq 0.4675$, and 20 percent for $M= 10^{15.0} ~M_\odot$ and $M=10^{15.2}~M_\odot$ in the redshift range $0.0625\leq z \leq 0.4675$. The redshift estimate accuracy is consistently below $\Delta z /z \leq 0.05$ for the above halo masses in the range $0.1525\leq z \leq 0.4675$. We also demonstrate that the orientation of triaxial halos and systematic error in our halo model do not affect reconstruction result significantly. Finally, we present results from reconstruction of mass distribution using shear catalogs produced by multiple halos, to show \splinv{}'s capability using realistic shear maps from ongoing and future galaxy surveys., Comment: 20 pages, 20 figures

The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function (GSMF) down to extremely faint and low mass galaxies ($M_{UV}>\sim-15$, $10^4M_\odot<=M_*<=10^8 M_\odot)$. We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint-end of the UVLF at $z=12$. We demonstrate that JWST observations are close to this enhanced regime, and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories., Comment: 13 pages, 7 figures. Accepted ApJL

Investigating the cosmological implication of the fact that neutrino has finite mass is of importance for fundamental physics. In particular, massive neutrino affects the formation of the large-scale structure (LSS) of the universe, and, conversely, observations of the LSS can give constraints on the neutrino mass. Numerical simulations of the LSS formation including massive neutrino along with conventional cold dark matter is thus an important task. For this, calculating the neutrino distribution in the phase space by solving the Vlasov equation is a suitable approach, but it requires solving the PDE in the $(6+1)$-dimensional space and is thus computationally demanding: Configuring $n_\mathrm{gr}$ grid points in each coordinate and $n_t$ time grid points leads to $O(n_\mathrm{gr}^6)$ memory space and $O(n_tn_\mathrm{gr}^6)$ queries to the coefficients in the discretized PDE. We propose a quantum algorithm for this task. Linearizing the Vlasov equation by neglecting the relatively weak self-gravity of the neutrino, we perform the Hamiltonian simulation to produce quantum states that encode the phase space distribution of neutrino. We also propose a way to extract the power spectrum of the neutrino density perturbations as classical data from the quantum state by quantum amplitude estimation with accuracy $\epsilon$ and query complexity of order $\widetilde{O}((n_\mathrm{gr} + n_t)/\epsilon)$. Our method also reduces the space complexity to $O(\mathrm{polylog}(n_\mathrm{gr}/\epsilon))$ in terms of the qubit number, while using quantum random access memories with $O(n_\mathrm{gr}^3)$ entries. As far as we know, this is the first quantum algorithm for the LSS simulation that outputs the quantity of practical interest with guaranteed accuracy.

Recent observations by the James Webb Space Telescope (JWST) discovered unexpectedly abundant luminous galaxies at high redshift, posing possibly a severe challenge to popular galaxy formation models. We study early structure formation in a cosmological model with a blue, tilted power spectrum (BTPS) given by $P(k) \propto k^{m_{\rm s}}$ with $m_{\rm s} > 1$ at small length scales. We run a set of cosmological $N$-body simulations and derive the abundance of dark matter halos and galaxies under simplified assumptions on star formation efficiency. The enhanced small-scale power allows rapid nonlinear structure formation at $z>7$, and galaxies with stellar mass exceeding $10^{10}\,M_\odot$ can be formed by $z=9$. Because of frequent mergers, the structure of galaxies and galaxy groups appears clumpy. The BTPS model reproduces the observed stellar mass density at $z=7-9$, and thus eases the claimed tension between galaxy formation theory and recent JWST observations. The large-scale structure of the present-day Universe is largely unaffected by the modification of the small-scale power spectrum. We conduct a systematic study by varying the slope of the small-scale power spectrum to derive constraints on the BTPS model from a set of observations of high-redshift galaxies., Comment: 10 pages, 7 figures, 1 table, accepted for publication in ApJ

The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficiently despite forming outside of dark matter halos. We investigate the potential for star formation in SIGOs using cosmological hydrodynamic simulations, including the aforementioned relative motions of baryons and dark matter, molecular hydrogen cooling in primordial gas clouds, and including explicit star formation. We find that SIGOs do form stars and that the nascent star clusters formed through this process are accreted by dark matter halos on short timescales (a few hundreds of Myr). Thus, SIGOs may be found as intact substructures within these halos, analogous to many present-day GCs. From this result, we conclude that SIGOs are capable of forming star clusters with similar properties to globular clusters in the early Universe and we discuss their detectablity by upcoming JWST surveys., Comment: 12 pages, 5 figures

The protoclusters in the epoch of reionization, traced by galaxies overdensity regions, are ideal laboratories for studying the process of stellar assembly and cosmic reionization. We present the spectroscopic confirmation of the core of the most distant protocluster at $z = 7.88$, A2744-z7p9OD, with the James Webb Space Telescope NIRSpec integral field unit spectroscopy. The core region includes as many as 4 galaxies detected in [OIII] 4960 \AA\ and 5008 \AA\ in a small area of $\sim 3\arcsec \times 3\arcsec$, corresponding to $\sim$ 11 kpc $\times$ 11 kpc, after the lensing magnification correction. Three member galaxies are also tentatively detected in dust continuum in Atacama Large Millimeter/submillimeter Array Band 6, which is consistent with their red ultraviolet continuum slopes, $\beta \sim -1.3$. The member galaxies have stellar masses in the range of log($M_{*}/M_{\rm \odot}$) $\sim 7.6-9.2$ and star formation rates of $\sim 3-50$ $M_{\rm \odot}$ yr$^{-1}$, showing a diversity in their properties. FirstLight cosmological simulations reproduce the physical properties of the member galaxies including the stellar mass, [OIII] luminosity, and dust-to-stellar mass ratio, and predict that the member galaxies are on the verge of merging in a few to several tens Myr to become a large galaxy with $M_{\rm *}\sim 6\times10^{9} M_{\rm \odot}$. The presence of a multiple merger and evolved galaxies in the core region of A2744-z7p9OD indicates that environmental effects are already at work 650 Myr after the Big Bang., Comment: 12 pages, 4 figures, 1 table, published in The Astrophysical Journal Letters

The collisionless Boltzmann equation (CBE) is a fundamental equation that governs the dynamics of a broad range of astrophysical systems from space plasma to star clusters and galaxies. It is computationally expensive to integrate the CBE directly in a multi-dimensional phase space, and thus the applications to realistic astrophysical problems have been limited so far. Recently, Todorova & Steijl (2020) proposed an efficient quantum algorithm to solve the CBE with significantly reduced computational complexity. We extend the algorithm to perform quantum simulations of self-gravitating systems, incorporating the method to calculate gravity with the major Fourier modes of the density distribution extracted from the solution-encoding quantum state. Our method improves the dependency of time and space complexities on Nv , the number of grid points in each velocity coordinate, compared to the classical simulation methods. We then conduct some numerical demonstrations of our method. We first run a 1+1 dimensional test calculation of free streaming motion on 64*64 grids using 13 simulated qubits and validate our method. We then perform simulations of Jeans collapse, and compare the result with analytic and linear theory calculations. It will thus allow us to perform large-scale CBE simulations on future quantum computers., Comment: 12 pages, 10 figures

An array of large observational programs using ground-based and space-borne telescopes is planned in the next decade. The forthcoming wide-field sky surveys are expected to deliver a sheer volume of data exceeding an exabyte. Processing the large amount of multiplex astronomical data is technically challenging, and fully automated technologies based on machine learning and artificial intelligence are urgently needed. Maximizing scientific returns from the big data requires community-wide efforts. We summarize recent progress in machine learning applications in observational cosmology. We also address crucial issues in high-performance computing that are needed for the data processing and statistical analysis., Comment: 55 pages, 8 figures, accepted for publication in Reports on Progress in Physics

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

We present integral field, far-infrared (FIR) spectroscopy of Mrk 54, a local Lyman Continuum Emitter (LCE), obtained with FIFI-LS on the Stratospheric Observatory for Infrared Astronomy. This is only the second time, after Haro 11, that [C II] 158 $\mu$m and [O III] 88 $\mu$m spectroscopy of the known LCEs have been obtained. We find that Mrk 54 has a strong [C II] emission that accounts for $\sim1$% of the total FIR luminosity, whereas it has only moderate [O III] emission, resulting in the low [O III]/[C II] luminosity ratio of $0.22\pm0.06$. In order to investigate whether [O III]/[C II] is a useful tracer of $f_{\rm esc}$ (LyC escape fraction), we examine the correlations of [O III]/[C II] and (i) the optical line ratio of $\rm O_{32} \equiv$ [O III] 5007 \AA/[O II] 3727 \AA, (ii) specific star formation rate, (iii) [O III] 88 $\mu$m/[O I] 63 $\mu$m ratio, (iv) gas phase metallicity, and (v) dust temperature based on a combined sample of Mrk 54 and the literature data from the Herschel Dwarf Galaxy Survey and the LITTLE THINGS Survey. We find that galaxies with high [O III]/[C II] luminosity ratios could be the result of high ionization (traced by $\rm O_{32}$), bursty star formation, high ionized-to-neutral gas volume filling factors (traced by [O III] 88 $\mu$m/[O I] 63 $\mu$m), and low gas-phase metallicities, which is in agreement with theoretical predictions. We present an empirical relation between the [O III]/[C II] ratio and $f_{\rm esc}$ based on the combination of the [O III]/[C II] and $\rm O_{32}$ correlation, and the known relation between $\rm O_{32}$ and $f_{\rm esc}$. The relation implies that high-redshift galaxies with high [O III]/[C II] ratios revealed by ALMA may have $f_{\rm esc}\gtrsim0.1$, significantly contributing to the cosmic reionization., Comment: 14 pages, 5 figures, Accepted for publication in ApJ

Tensor decomposition is now being used for data analysis, information compression, and knowledge recovery. However, the mathematical property of tensor decomposition is not yet fully clarified because it is one of singular learning machines. In this paper, we give the upper bound of its real log canonical threshold (RLCT) of the tensor decomposition by using an algebraic geometrical method and derive its Bayesian generalization error theoretically. We also give considerations about its mathematical property through numerical experiments.

Recent observations by James Webb Space Telescope discovered a number of high-redshift galaxies with strong emission lines from doubly ionized oxygen. Combined with ALMA observations of far-infrared lines, multi-line diagnostics can be applied to the high-redshift galaxies in order to probe the physical conditions of the inter-stellar medium. We study the formation and evolution of galaxies using the FirstLight simulation suite, which provides outputs of 62 high-resolution, zoom-in galaxy simulations. We devise a physical model of HII regions and calculate spatially resolved [OIII] line emission. We show that massive galaxies with stellar masses of $M_* > 10^9 M_\odot$ chemically evolve rapidly to $z=9$. Young stellar populations in the star-forming galaxies boost the [OIII] line emission, rendering the ratio of line luminosity to star formation rate larger than that for low-redshift galaxies, which is consistent with recent observations. Measuring the flux ratios of rest-frame optical and far-infrared lines allows us to estimate the physical conditions such as density and metallicity of the star-forming gas in high-redshift [OIII] emitters., Comment: 10 pages, 5 figures, 1 table, ApJ in press

A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because of a physical offset between the centers of mass of DM and baryonic overdensities. We present an updated numerical investigation of DM GHOSts and Supersonically Induced Gas Objects (SIGOs), including the effects of molecular cooling, in high resolution hydrodynamic simulations using the AREPO code. Supplemented by an analytical understanding of their ellipsoidal gravitational potentials, we study the population-level properties of these objects, characterizing their morphology, spin, radial mass, and velocity distributions in comparison to classical structures in non-streaming regions. The stream velocity causes deviations from sphericity in both the gas and DM components and lends greater rotational support to the gas. Low mass ($<\sim 10^{5.5}$ M$_\odot$) objects in regions of streaming demonstrate core-like rotation and mass profiles. Anomalies in the rotation and morphology of DM GHOSts could represent an early Universe analogue to observed ultra-faint dwarf galaxies with variations in DM content and unusual rotation curves., Comment: 27 pages, 20 figures. Accepted ApJ

Supersonically Induced Gas Objects (SIGOs) are a class of early Universe objects that have gained attention as a potential formation route for globular clusters. SIGOs have only recently begun to be studied in the context of molecular hydrogen cooling, which is key to characterizing their structure and evolution. Studying the population-level properties of SIGOs with molecular cooling is important for understanding their potential for collapse and star formation, and central for addressing whether SIGOs can survive to the present epoch. Here, we investigate the evolution of SIGOs before they form stars, using a combination of numerical and analytical analysis. For example, we study various timescales important to the evolution of SIGOs at a population level in the presence of molecular cooling. Revising the previous formulation for the critical density of collapse for SIGOs allows us to show that their prolateness tends to act as an inhibiting factor to collapse. We find that simulated SIGOs are limited by artificial two-body relaxation effects that tend to disperse them, an effect of their limited resolution. We expect that SIGOs in nature will be longer-lived compared to our simulations. Further, the fall-back timescale on which SIGOs fall into nearby dark matter halos, potentially producing a globular-cluster-like system, is frequently longer than their cooling timescale and the collapse timescale on which they shrink through gravity. Therefore, some SIGOs have time to cool and collapse outside of halos despite initially failing to exceed the critical density, even without considering metal line cooling. From this analysis we conclude that SIGOs should form stars outside of halos in non-negligible stream velocity patches in the Universe., Comment: 18 pages, 11 figures

In this series of papers, we present a simulation-based model for the non-linear clustering of galaxies based on separate modelling of clustering in real space and velocity statistics. In the first paper, we present an emulator for the real-space correlation function of galaxies, whereas the emulator of the real-to-redshift space mapping based on velocity statistics is presented in the second paper. Here, we show that a neural network emulator for real-space galaxy clustering trained on data extracted from the Dark Quest suite of N-body simulations achieves sub-per cent accuracies on scales $1 < r < 30 $ $h^{-1} \,\mathrm{Mpc}$, and better than $3\%$ on scales $r < 1$ $h^{-1}\mathrm{Mpc}$ in predicting the clustering of dark-matter haloes with number density $10^{-3.5}$ $(h^{-1}\mathrm{Mpc})^{-3}$, close to that of SDSS LOWZ-like galaxies. The halo emulator can be combined with a galaxy-halo connection model to predict the galaxy correlation function through the halo model. We demonstrate that we accurately recover the cosmological and galaxy-halo connection parameters when galaxy clustering depends only on the mass of the galaxies' host halos. Furthermore, the constraining power in $\sigma_8$ increases by about a factor of $2$ when including scales smaller than $5$ $h^{-1} \,\mathrm{Mpc}$. However, when mass is not the only property responsible for galaxy clustering, as observed in hydrodynamical or semi-analytic models of galaxy formation, our emulator gives biased constraints on $\sigma_8$. This bias disappears when small scales ($r < 10$ $h^{-1}\mathrm{Mpc}$) are excluded from the analysis. This shows that a vanilla halo model could introduce biases into the analysis of future datasets.

We present a deep neural network Real/Bogus classifier that improves classification performance in the Tomo-e Gozen transient survey by handling label errors in the training data. In the wide-field, high-frequency transient survey with Tomo-e Gozen, the performance of conventional convolutional neural network classifier is not sufficient as about $10^6$ bogus detections appear every night. In need of a better classifier, we have developed a new two-stage training method. In this training method, label errors in the training data are first detected by normal supervised learning classification, and then they are unlabeled and used for training of semi-supervised learning. For actual observed data, the classifier with this method achieves an area under the curve (AUC) of 0.9998 and a false positive rate (FPR) of 0.0002 at true positive rate (TPR) of 0.9. This training method saves relabeling effort by humans and works better on training data with a high fraction of label errors. By implementing the developed classifier in the Tomo-e Gozen pipeline, the number of transient candidates was reduced to $\sim$40 objects per night, which is $\sim$1/130 of the previous version, while maintaining the recovery rate of real transients. This enables more efficient selection of targets for follow-up observations., Comment: 14 pages, 17 figures, 2 tables. Published in PASJ. The source code is available at https://github.com/ichiro-takahashi/tomoe-realbogus

The semi-analytical model A-SLOTH (Ancient Stars and Local Observables by Tracing Halos) is the first public code that connects the formation of the first stars and galaxies to observables. After several successful projects with this model, we publish the source code and describe the public version in this paper. The model is based on dark matter merger trees that can either be generated based on Extended Press-Schechter theory or that can be imported from dark matter simulations. On top of these merger trees, A-SLOTH applies analytical recipes for baryonic physics to model the formation of both metal-free and metal-poor stars and the transition between them with unprecedented precision and fidelity. A-SLOTH samples individual stars and includes radiative, chemical, and mechanical feedback. It is calibrated based on six observables, such as the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely-metal poor stars, and the cosmic star formation rate density at high redshift. A-SLOTH has versatile applications with moderate computational requirements. It can be used to constrain the properties of the first stars and high-z galaxies based on local observables, predicts properties of the oldest and most metal-poor stars in the Milky Way, can serve as a subgrid model for larger cosmological simulations, and predicts next-generation observables of the early Universe, such as supernova rates or gravitational wave events., Comment: Published in ApJ, 35 pages, video summary can be watched at https://www.youtube.com/watch?v=M_JRT2AsOqg and git repository is available at https://gitlab.com/thartwig/asloth

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

We infer the expected detection number of pair instability supernovae (PISNe) during the operation of the Euclid space telescope based on binary population models. Our models reproduce the global maximum of the rate at the primary BH mass of $\sim 9-10$ $M_\odot$, and the overall gradient of the primary BH mass distribution in the binary BH merger rate consistent with recent observations. We consider different PISN conditions depending on the $^{12}$C$(\alpha, \gamma)^{16}$O reaction rate. The fiducial and $3\sigma$ models adopt the standard and $3\sigma$-smaller reaction rate, respectively. Our fiducial model predicts that Euclid detects several hydrogen-poor PISNe. For the $3\sigma$ model, detection of $\sim 1$ hydrogen-poor PISN by Euclid is expected if the stellar mass distribution extends to $M_{\max} = 600 M_\odot$, but the expected number becomes significantly smaller if $M_{\max} = 300 M_\odot$. We may be able to distinguish the fiducial and $3\sigma$ models by the observed PISN rate. This will help us to constrain the origin of binary BHs and the reaction rate, although there remains degeneracy between $M_{\max}$ and the reaction rate. PISN ejecta mass estimates from light curves and spectra obtained by follow-up observations would be important to disentangle the degeneracy., Comment: 7 pages, 4 figures, 1 table, accepted for MNRAS letters

We perform radiation hydrodynamics simulations to study the structure and evolution of a photoevaporating protoplanetary disk. Ultraviolet and X-ray radiation from the host star heats the disk surface, where H$_2$ pumping also operates efficiently. We run a set of simulations with varying the amount of dust grains, or the dust-to-gas mass ratio, which determines the relative importance between photoelectric heating and H$_2$ pumping. We show that H$_2$ pumping and X-ray heating contribute stronger to the mass-loss of the disk if the dust-to-gas mass ratio is $\mathcal{D}\leq10^{-3}$.The disk mass loss rate decreases with a lower dust amount, but remains around $10^{-10-11} M_{\odot} {\rm yr}^{-1}$. In such dust-deficient disks, H$_2$ pumping enhances photoevaporation from the inner disk region and shapes the disk mass-loss profile. We thus argue that the late-stage disk evolution is affected by the ultra-violet H$_2$ pumping effect. The mass-loss rates derived from our simulations can be used in the study of long-term disk evolution., Comment: 12 pages, 3 figures, 1 table

We study the formation and gravitational collapse of supersonically induced gas objects (SIGOs) in the early universe. We run cosmological hydrodynamics simulations of SIGOs, including relative streaming motions between baryons and dark matter. Our simulations also follow nonequilibrium chemistry and molecular hydrogen cooling in primordial gas clouds. A number of SIGOs are formed in the run with fast-streaming motions of 2 times the rms of the cosmological velocity fluctuations. We identify a particular gas cloud that condensates by H$_2$ cooling without being hosted by a dark matter halo. The SIGO remains outside the virial radius of its closest halo, and it becomes Jeans unstable when the central gas-particle density reaches $\sim 100~{\rm cm}^{-3}$ with a temperature of $\sim$ 200 K. The corresponding Jeans mass is $\sim 10^5 M_{\odot}$, and thus the formation of primordial stars or a star cluster is expected in the SIGO., Comment: 8 pages, 4 figures, Accepted for publication in ApJL

We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($\Delta\!\Sigma$) and projected galaxy clustering ($w_{\rm p}$) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range $0.15

Atmospheric escape of close-in exoplanets can be driven by high energy radiation from the host star. The planetary outflows interacting with the stellar wind may generate observable transit signals that depend on the strength of the stellar wind. We perform detailed radiation-hydrodynamics simulations of the atmospheric escape of hot Jupiters with including the wind from the host star in a self-consistent, dynamically coupled manner. We show that the planetary outflow is shaped by the balance between its thermal pressure and the ram pressure of the stellar wind. We use the simulation outputs to calculate the Lyman-$\alpha$ and H$\alpha$ transit signatures. Strong winds can confine the outflow and decrease the Lyman-$\alpha$ transit depth. Contrastingly, the wind effect on H$\alpha$ is weak because of the small contribution from the uppermost atmosphere of the planet. Observing both of the lines is important to understand the effect of the UV radiation and wind from the host. The atmospheric mass-loss rate is approximately independent of the strength of the wind. We also discuss the effect of the coronal mass ejections on the signatures. We argue that around M dwarfs the effect can be significant in every transit., Comment: 5 pages, 4 figures, accepted for publication in MNRAS

We report a Vlasov simulation of cosmic relic neutrinos combined with N-body simulation of cold dark matter in the context of large-scale structure formation in the Universe performed on Fugaku supercomputer. Gravitational dynamics of the neutrinos is followed, for the first time, by directly integrating the Vlasov equation in a six-dimensional phase space. Our largest simulation combines the Vlasov simulation on 400 trillion grids and 330 billion-body calculations in a self-consistent manner, and reproduces accurately the nonlinear dynamics of neutrinos in the Universe. The novel high-order Vlasov solver is optimized by combining an array of state-of-the-art numerical schemes and fully utilizing the SIMD instructions on the A64FX processors. Time-To-Solution of our simulation is an order of magnitude shorter than the largest N-body simulations. The performance scales excellently with up to 147,456 nodes (7 million CPU cores) on Fugaku; the weak and strong scaling efficiencies are 82% - 96% and 82% - 93%, respectively., Comment: 11 pages, 8 figures. 2021 ACM Gordon-Bell Prize Finalist. A simulation movie is available at https://youtu.be/p2fHt8f7chQ

Line intensity mapping is emerging as a novel method that can measure the collective intensity fluctuations of atomic/molecular line emission from distant galaxies. Several observational programs with various wavelengths are ongoing and planned, but there remains a critical problem of line confusion; emission lines originating from galaxies at different redshifts are confused at the same observed wavelength. We devise a generative adversarial network that extracts designated emission line signals from noisy three-dimensional data. Our novel network architecture allows two input data, in which the same underlying large-scale structure is traced by two emission lines of H$\rm \alpha$ and [OIII], so that the network learns the relative contributions at each wavelength and is trained to decompose the respective signals. After being trained with a large number of realistic mock catalogs, the network is able to reconstruct the three-dimensional distribution of emission-line galaxies at $z = 1.3-2.4$. Bright galaxies are identified with a precision of 84%, and the cross-correlation coefficients between the true and reconstructed intensity maps are as high as 0.8. Our deep-learning method can be readily applied to data from planned space-borne and ground-based experiments., Comment: 7 pages, 5 figures, ApJL in press

The member stars in globular cluster M15 show a substantial spread in the abundances of r-process elements. We argue that a rare and prolific r-process event enriched the natal cloud of M15 in an inhomogeneous manner. To critically examine the possibility, we perform cosmological galaxy formation simulations and study the physical conditions for the inhomogeneous enrichment. We explore a large parameter space of the merger event time and the site. Our simulations reproduce the large r-process abundance spread if a neutron-star merger occurs at \sim 100 pc away from the formation site of the cluster and in a limited time range of a few tens million years before the formation. Interestingly, a bimodal feature is found in the Eu abundance distribution in some cases, similarly to that inferred from recent observations. M15 member stars do not show clear correlation between the abundances of Eu and light elements such as Na that is expected in models with two stellar populations. We thus argue that a majority of the stars in M15 are formed in a single burst. The ratio of heavy to light r-process element abundance [Eu/Y] \sim 1.0 is consistent with that of the so-called r-II stars, suggesting that a lanthanide-rich r-process event dominantly enriched M15., Comment: 9 pages, 5 figures, accepted for publication in ApJ Letters

Studies of low redshift galaxy clusters suggest the intra-cluster medium (ICM) has experienced non-gravitational heating during the formation phase of the clusters. Using simple phenomenological heating prescriptions, we simulate the effect of this preheating of the nascent ICM in galaxy proto-clusters and examine its effect on Lyman-$\alpha$ (Ly$\alpha$) forest tomographic maps. We analyse a series of cosmological zoom-in simulations of proto-clusters within the framework of the Ly$\alpha$ transmission-dark matter (DM) density distribution. We find that the more energy is injected into the proto-ICM at $z$ = 3, the more the distribution at high DM density tilts towards higher Ly$\alpha$ transmission. This effect has been confirmed in both low-resolution simulations adopting a preheating scheme based on entropy floors, as well as in higher-resolution simulations with another scheme based on energy floors. The evolution of the slope of this distribution is shown to vary with redshift. The methodology developed here can be applied to current and upcoming Ly$\alpha$ forest tomographic survey data to help constrain feedback models in galaxy proto-clusters., Comment: Accepted for publication in ApJ, 24 pages, 16 figures, 1 table

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|\delta m|<9\times 10^{-3}$. The bias residual $\delta m$ shows no dependence on redshift in the range $01$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements., Comment: 39 pages, 31 figures, 5 tables; PASJ (in press)

We study the fragmentation of collisional ring galaxies (CRGs) using a linear perturbation analysis that computes the physical conditions of gravitational instability, as determined by the balance of self-gravity of the ring against pressure and Coriolis forces. We adopt our formalism to simulations of CRGs and show that the analysis can accurately characterise the stability and onset of fragmentation, although the linear theory appears to under-predict the number of fragments of an unstable CRG by a factor of 2. In addition, since the orthodox `density-wave' model is inapplicable to such self-gravitating rings, we devise a simple approach that describes the rings propagating as material waves. We find that the toy model can predict whether the simulated CRGs fragment or not using information from their pre-collision states. We also apply our instability analysis to a CRG discovered at a high redshift, $z=2.19$. We find that a quite high velocity dispersion is required for the stability of the ring, and therefore the CRG should be unstable to ring fragmentation. CRGs are rarely observed at high redshifts, and this may be because CRGs are usually too faint. Since the fragmentation can induce active star formation and make the ring bright enough to observe, the instability could explain this rarity. An unstable CRG fragments into massive clumps retaining the initial disc rotation, and thus it would evolve into a clumpy galaxy with a low surface density in an inter-clump region., Comment: 7 pages, 5 figures (plus Appendix). Accepted for publication in MNRAS

Ultra-faint dwarf galaxies (UFDs) are promising observable proxies to building blocks of galaxies formed in the early Universe. We study the formation and evolution of UFDs using cosmological hydrodynamic simulations. In particular, we show that a major merger of two building block galaxies with 3,900 Msun and 7,500 Msun at the cosmic age of 510 Myr results in a system with an extended stellar distribution consistent with the de Vaucouleurs profile. The simulated galaxy has an average stellar metallicity of [Fe/H]=-2.7 and features a metallicity gradient. These results closely resemble the properties of a recently discovered UFD, Tucana II, which is extremely metal-poor and has a spatially extended stellar halo with the more distant stars being more metal-poor. Our simulation suggests that the extended stellar halo of Tucana II may have been formed through a past major merger. Future observational searches for spatially extended structures around other UFDs, combined with further theoretical studies, will provide tangible measures of the evolutionary history of the ancient, surviving satellite galaxies., Comment: 9 pages, 4 figures, ApJ Letters in press

We propose to use relative strengths of far-infrared fine structure lines from galaxies to characterise early phases of the inside-out quenching by massive black holes (BHs). The BH feedback is thought to quench star formation by evacuating the ambient gas. In order to quantify the feedback effect on the gas density in the galactic centres, we utilise the outputs of IllustrisTNG and Illustris simulations, which implement different BH feedback models. We devise a physical model of H$_{\rm ~II}$ regions and compute the intensities of [O$_{\rm ~III}$] $52$ and $88~{\rm \mu m}$ lines. The line intensity ratio is sensitive to the local electron density, and thus can be used to measure the strength and physical extent of the BH quenching. If the BH feedback abruptly operates and expel the gas when it grows to a certain mass, as modelled in IllustrisTNG, the low-density gas yields relatively weak [O$_{\rm ~III}$] $52$ line with respect to $88~{\rm \mu m}$. In contrast, if the feedback strength and hence the local gas density are not strongly correlated with the BH mass, as in Illustris, the line ratio is not expected to vary significantly among galaxies with different evolutionary stages. We find these features are reproduced in the simulations. We also show that the line ratios are not sensitive to the aperture size for measurement, and thus observations do not need to resolve the galactic centres. We argue that the integrated line ratios can be used to capture the onset of the inside-out quenching by BHs., Comment: 13 pages, 15 figures (plus Appendix), submitted to MNRAS