Search for: All records

Cyclical MCMC is a novel MCMC framework recently proposed by Zhang et al. (2019) to address the challenge posed by high-dimensional multimodal posterior distributions like those arising in deep learning. The algorithm works by generating a nonhomogeneous Markov chain that tracks -- cyclically in time -- tempered versions of the target distribution. We show in this work that cyclical MCMC converges to the desired limit in settings where the Markov kernels used are fast mixing, and sufficiently long cycles are employed. However in the far more common settings of slow mixing kernels, the algorithm may fail to converge to the correct limit. In particular, in a simple mixture example with unequal variance where powering is known to produce slow mixing kernels, we show by simulation that cyclical MCMC fails to converge to the desired limit. Finally, we show that cyclical MCMC typically estimates well the local shape of the target distribution around each mode, even when we do not have convergence to the target. 

An important characteristic of cosmic hydrogen reionization is the growth of ionized gas bubbles surrounding early luminous objects. Ionized bubble sizes are beginning to be probed using Lyman α emission from high-redshift galaxies, and will also be probed by upcoming 21 cm maps. We present results from a study of bubble sizes using the state-of-the-art thesan radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation. We employ the mean free path method and track the evolution of the effective ionized bubble size at each point (Reff) throughout the Epoch of Reionization. We show that there is a slow growth period for regions ionized early, but a rapid ‘flash ionization’ process for regions ionized later as they immediately enter a large, pre-existing bubble. We also find that bright sources are preferentially in larger bubbles, and find consistency with recent observational constraints at z ≳ 9, but tension with idealized Lyman α damping-wing models at z ≈ 7. We find that high-overdensity regions have larger characteristic bubble sizes, but the correlation decreases as reionization progresses, likely due to runaway formation of large percolated bubbles. Finally, we compare the redshift at which a region transitions from neutral to ionized (zreion) with the time it takes to reach a given bubble size and conclude that zreion is a reasonable local probe of small-scale bubble size statistics ($R_\text{eff} \lesssim 1\, \rm {cMpc}$). However, for larger bubbles, the correspondence between zreion and size statistics weakens due to the time delay between the onset of reionization and the expansion of large bubbles, particularly at high redshifts.

 

We investigate how feedback and environment shapes the X-ray scaling relations of early-type galaxies (ETGs), especially at the low-mass end. We select central-ETGs from the TNG100 box of IllustrisTNG that have stellar masses $\log _{10}(M_{\ast }/\mathrm{M_{\odot }})\in [10.7, 11.9]$. We derive mock X-ray luminosity (LX, 500) and spectroscopic-like temperature (Tsl, 500) of hot gas within R500 of the ETG haloes using the MOCK-X pipeline. The scaling between LX, 500 and the total mass within 5 effective radii ($M_{5R_{\rm e}}$) agrees well with observed ETGs from Chandra. IllustrisTNG reproduces the observed increase in scatter of LX, 500 towards lower masses, and we find that ETGs with $\log _{10} (M_{5R_{\rm e}}/\mathrm{M_{\odot }}) \leqslant 11.5$ with above-average LX, 500 experienced systematically lower cumulative kinetic AGN feedback energy historically (vice versa for below-average ETGs). This leads to larger gas mass fractions and younger stellar populations with stronger stellar feedback heating, concertedly resulting in the above-average LX, 500. The LX, 500–Tsl, 500 relation shows a similar slope to the observed ETGs but the simulation systematically underestimates the gas temperature. Three outliers that lie far below the LX–Tsl relation all interacted with larger galaxy clusters recently and demonstrate clear features of environmental heating. We propose that the distinct location of these backsplash ETGs in the LX–Tsl plane could provide a new way of identifying backsplash galaxies in future X-ray surveys.

 

It has been claimed that traditional models struggle to explain the tentative detection of the 21 cm absorption trough centered atz∼ 17 measured by the EDGES collaboration. On the other hand, it has been shown that the EDGES results are consistent with an extrapolation of a declining UV luminosity density, following a simple power law of deep Hubble Space Telescope observations of 4 <z< 9 galaxies. We here explore the conditions by which the EDGES detection is consistent with current reionization and post-reionization observations, including the neutral hydrogen fraction atz∼ 6–8, Thomson-scattering optical depth, and ionizing emissivity atz∼ 5. By coupling a physically motivated source model derived from radiative transfer hydrodynamic simulations of reionization to a Markov Chain Monte Carlo sampler, we find that it is entirely possible to reconcile existing high-redshift (cosmic dawn) and low-redshift (reionization) constraints. In particular, we find that high contributions from low-mass halos along with high photon escape fractions are required to simultaneously reproduce cosmic dawn and reionization constraints. Our analysis further confirms that low-mass galaxies produce a flatter emissivity evolution, which leads to an earlier onset of reionization with a gradual and longer duration, resulting in a higher optical depth. While the models dominated by faint galaxies successfully reproduce the measured globally averaged quantities over the first one billion years, they underestimate the late redshift-instantaneous measurements in efficiently star-forming and massive systems. We show that our (simple) physically motivated semianalytical prescription produces results that are consistent with the (sophisticated) state-of-the-artTHESANradiation-magnetohydrodynamic simulation of the reionization.

 

Modelling galaxy formation in hydrodynamic simulations has increasingly adopted various radiative transfer methods to account for photoionization feedback from young massive stars. However, the evolution of H ii regions around stars begins in dense star-forming clouds and spans large dynamical ranges in both space and time, posing severe challenges for numerical simulations in terms of both spatial and temporal resolution that depends strongly on gas density (∝n−1). In this work, we perform a series of idealized H ii region simulations using the moving-mesh radiation-hydrodynamic code arepo-rt to study the effects of numerical resolution. The simulated results match the analytical solutions and the ionization feedback converges only if the Strömgren sphere is resolved by at least 10–100 resolution elements and the size of each time integration step is smaller than 0.1 times the recombination time-scale. Insufficient spatial resolution leads to reduced ionization fraction but enhanced ionized gas mass and momentum feedback from the H ii regions, as well as degrading the multiphase interstellar medium into a diffuse, partially ionized, warm (∼8000 K) gas. On the other hand, insufficient temporal resolution strongly suppresses the effects of ionizing feedback. This is because longer time-steps are not able to resolve the rapid variation of the thermochemistry properties of the gas cells around massive stars, especially when the photon injection and thermochemistry are performed with different cadences. Finally, we provide novel numerical implementations to overcome the above issues when strict resolution requirements are not achievable in practice.

 

Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization during the first billion years of cosmic history. This is especially relevant now that a fleet of advanced observatories (e.g. James Webb Space Telescope, Nancy Grace Roman Space Telescope, SPHEREx, ELT, SKA) will soon provide an holistic picture of this defining epoch. To bridge this gap, we publicly release all simulation outputs and post-processing products generated within the thesan simulation project at www.thesan-project.com. This project focuses on the z ≥ 5.5 Universe, combining a radiation-hydrodynamics solver (arepo-rt), a well-tested galaxy formation model (IllustrisTNG) and cosmic dust physics to provide a comprehensive view of the Epoch of Reionization. The thesan suite includes 16 distinct simulations, each varying in volume, resolution, and underlying physical models. This paper outlines the unique features of these new simulations, the production and detailed format of the wide range of derived data products, and the process for data retrieval. Finally, as a case study, we compare our simulation data with a number of recent observations from the James Webb Space Telescope, affirming the accuracy and applicability of thesan. The examples also serve as prototypes for how to utilize the released data set to perform comparisons between predictions and observations.

 

Abstract The darkling beetle tribe Adesmiini (Tenebrionidae: Pimeliinae) is a prominent part of African and western Palearctic desert faunas, with most species being day-active fast-running detritivores. Taxonomic diversity within the tribe is highest in the southern Afrotropical realm (where all genera are present); only 1 genus, the species-rich Adesmia, occurs north of the Sahara. Despite notable species, such as the fog-basking beetle Onymacris unguicularis (a focal taxon in desert ecological research), Adesmiini has undergone few modern taxonomic or phylogenetic studies. Hence, generic concepts and pronounced diurnal activity, rare in the primarily nocturnal family Tenebrionidae, remain poorly explored. To investigate evolutionary relationships and diurnal origins within the tribe, we generated a genomic dataset of 529 protein-coding genes across 43 species spanning 10 of 11 Adesmiini genera. Our resulting phylogeny for the tribe rejects the monophyly of 5 currently recognized Adesmiini genera (i.e., Adesmia, Metriopus, Onymacris, Physadesmia, and Stenocara). Ancestral state reconstruction of diurnal activity using eye shape as a proxy supports the hypothesis that Adesmiini were primitively diurnal, followed by at least 4 shifts to nocturnal or crepuscular activity. 

Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (thesan-hr), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matter (FDM), and interacting dark matter (IDM) with strong dark acoustic oscillations (sDAO). In altDM scenarios, both the halo mass functions and the ultraviolet luminosity functions at z ≳ 6 are suppressed at the low-mass/faint end, leading to delayed global star formation and reionization histories. However, strong non-linear effects enable altDM models to ‘catch up’ with cold dark matter (CDM) in terms of star formation and reionization. The specific star formation rates are enhanced in halos below the half-power mass in altDM models. This enhancement coincides with increased gas abundance, reduced gas depletion times, more compact galaxy sizes, and steeper metallicity gradients at the outskirts of the galaxies. These changes in galaxy properties can help disentangle altDM signatures from a range of astrophysical uncertainties. Meanwhile, it is the first time that altDM models have been studied in RHD simulations of galaxy formation. We uncover significant systematic uncertainties in reionization assumptions on the faint-end luminosity function. This underscores the necessity of accurately modeling the small-scale morphology of reionization in making predictions for the low-mass galaxy population. Upcoming James Webb Space Telescope imaging surveys of deep lensed fields hold potential for uncovering the faint low-mass galaxy population, which could provide constraints on altDM models.

 

The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of mb ≈ 104 M⊙), to additionally enable the study of ‘mini-haloes’ with virial temperatures Tvir < 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M1500 =   − 10) by z = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.

 

Abstract A fundamental requirement for reionizing the Universe is that a sufficient fraction of the ionizing photons emitted by galaxies successfully escapes into the intergalactic medium. However, due to the scarcity of high-redshift observational data, the sources driving reionization remain uncertain. In this work, we calculate the ionizing escape fractions (fesc) of reionization-era galaxies from the state-of-the-art thesan simulations, which combine an accurate radiation-hydrodynamic solver (arepo-rt) with the well-tested IllustrisTNG galaxy formation model to self-consistently simulate both small-scale galaxy physics and large-scale reionization throughout a large patch of the universe ($L_\text{box} = 95.5\, \text{cMpc}$). This allows the formation of numerous massive haloes ($M_\text{halo} \gtrsim 10^{10}\, {\text{M}_{\odot }}$), which are often statistically underrepresented in previous studies but are believed to be important to achieving rapid reionization. We find that low-mass galaxies ($M_\text{stars} \lesssim 10^7\, {\text{M}_{\odot }}$) are the main drivers of reionization above z ≳ 7, while high-mass galaxies ($M_\text{stars} \gtrsim 10^8\, {\text{M}_{\odot }}$) dominate the escaped ionizing photon budget at lower redshifts. We find a strong dependence of fesc on the effective star formation rate (SFR) surface density defined as the SFR per gas mass per escape area, i.e. $\bar{\Sigma }_\text{SFR} = \text{SFR}/M_\text{gas}/R_{200}^2$. The variation in halo escape fractions decreases for higher mass haloes, which can be understood from the more settled galactic structure, SFR stability, and fraction of sightlines within each halo significantly contributing to the escaped flux. Dust is capable of reducing the escape fractions of massive galaxies, but the impact on the global fesc depends on the dust model. Finally, active galactic nuclei are unimportant for reionization in thesan and their escape fractions are lower than stellar ones due to being located near the centres of galaxy gravitational potential wells.