Intergalactic Medium Research Articles

Intergalactic Medium Rotation Measure of Primordial Magnetic Fields

The Faraday rotation effect, quantified by the rotation measure (RM), is a powerful probe of the large-scale magnetization of the Universe—tracing magnetic fields not only on galaxy and galaxy cluster scales but also in the intergalactic medium (IGM; referred to as RMIGM). The redshift dependence of the latter has extensively been explored with observations. It has also been shown that this relation can help to distinguish between different large-scale magnetization scenarios. We study the evolution of this RMIGM for different primordial magnetogenesis scenarios to search for the imprints of primordial magnetic fields (PMFs; magnetic fields originating in the early Universe) on the redshift-dependence of RMIGM. We use cosmological magnetohydrodynamic simulations for evolving PMFs during large-scale structure formation, coupled with the light-cone analysis to produce a realistic statistical sample of mock RMIGM images. We study the predicted behavior for the cosmic evolution of RMIGM for different correlation lengths of PMFs, and provide fitting functions for their dependence on redshifts. We compare these mock RM trends with the recent analysis of the the LOw-Frequency ARray RM Grid and find that large-scale-correlated PMFs should have (comoving) strengths ≲0.75 nG, if they originated during inflation with the scale-invariant spectrum and (comoving) correlation length of ∼19 h −1 cMpc or ≲30 nG if they originated during phase-transition epochs with the comoving correlation length of ∼1 h −1 cMpc. Our findings agree with previous observations and confirm the results of semi-analytical studies, showing that upper limits on the PMF strength decrease as their coherence scales increase.

NIHAO-RiNG: A Comparison of Simulated Disk Galaxies from GASOLINE and GIZMO

We utilize the publicly available code GIZMO to re-simulate 12 galaxies selected from the Numerical Investigation of a Hundred Astronomical Object (NIHAO) simulation suite, which were run with the GASOLINE code, then compare their properties. We find that while both codes with the same initial conditions and large-scale environments can successfully produce similar galactic disks in Milky Way-mass systems, yet significant differences are still seen in many aspects, particularly the circumgalactic medium (CGM) environment they reside in. Specifically, the thermal feedback recipe used in GASOLINE results in ubiquitous, long-term, large-scale outflows, primarily driven by high-density hot interstellar medium from the galaxy center, preventing the intergalactic medium from falling efficiently. Recycled gas and inflows in the CGM appear at 104∼5 K, playing a crucial role in the formation of cold disks in the NIHAO simulations. In contrast, disk galaxies simulated by GIZMO do not exhibit prominent outflows at low redshifts, but instead display quasi-virialized hot gas halos that arise from the interaction between inflows and feedback-driven outflows. Therefore, the origins of mass and angular momentum of the cold disk in the two simulations are quite different, even though the final morphologies of the corresponding galaxies are both disky. The differences in the distribution of CGM gas are mainly due to different feedback models implemented in the two codes. Future observations of CGM could provide valuable insight into the physics governing the baryon cycle in galaxies.

Multiwavelength investigations of PKS 2300−18: S-shaped radio quasar with precessing jets and double peaked broad emission line spectrum

Abstract S-shaped radio galaxy jets are prime sources for investigating the dynamic interplay between the central active galactic nucleus, the jets, and the ambient intergalactic medium. These sources are excellent candidates for studying jet precession, as their S-shaped inversion symmetry strongly indicates underlying precession. We present a multiwavelength analysis of the giant inversion-symmetric S-shaped radio galaxy PKS 2300−18, which spans 0.76 Mpc. The host is a quasar at a redshift of 0.128, displaying disturbed optical morphology due to an ongoing merger with a companion galaxy. We conducted a broadband radio spectral study using multifrequency data ranging from 183 MHz to 6 GHz, incorporating dedicated observations with the uGMRT and JVLA alongside archival radio data. Particle injection model was fitted to the spectra of different regions of the source to perform ageing analysis, which was supplemented with a kinematic jet precession model. The ageing analysis revealed a maximum plasma age of ∼ 40 Myr, while the jet precession model indicated a precession period of ∼ 12 Myr. ROSAT data revealed an X-ray halo of Mpc size, and from Chandra the AGN X-ray spectrum was modelled using thermal and power law components. The optical spectrum displaying double peaked broad emission lines was modelled, indicating complex broad line region kinematics at the centre with the possibility of a binary SMBH. We present the results of our multiwavelength analysis of the source, spanning scales from a few light days to a few Mpc, and discuss its potential evolutionary path.

Constraining the Hubble constant with scattering in host galaxies of fast radio bursts

Measuring the Hubble constant ($ H $) is one of the most important missions in astronomy. Nevertheless, recent studies exhibit differences between the employed methods. Fast radio bursts (FRBs) are coherent radio transients with large dispersion measures (DM) with a duration of milliseconds. DM$_ IGM $, the free electron column density along a line of sight in the intergalactic medium (IGM), could open a new avenue for probing $ H $. However, it has been challenging to separate DM contributions from different components (i.e., the IGM and the host galaxy plasma), and this hampers the accurate measurements of DM$_ IGM $ and hence $ H $. We adopted a method to overcome this problem by using the temporal scattering of the FRB pulses due to the propagation effect through the host galaxy plasma (scattering time). The scattering-inferred DM in a host galaxy improves the estimate of DM$_ IGM $, which in turn leads to a better constraint on $ H $. In previous studies, a certain value or distribution has conventionally been assumed of the dispersion measure in host galaxies (DM$_ h $). We compared this method with ours by generating 100 mock FRBs, and we found that our method reduces the systematic (statistical) error of H$_0$ by 9.1$<!PCT!>$ (1$<!PCT!>$) compared to the previous method. We applied our method to 30 localized FRB sources with both scattering and spectroscopic redshift measurements to constrain H$_0$. Our result is H$_0$ = 74$_ $ km s$^ $ Mpc$^ $, where the central value prefers the value obtained from local measurements over the cosmic microwave background. We also measured DM$_ h $ with a median value of $103^ $ pc cm$^ The $ DM_h $ had to be assumed in previous works to derive DM$_ IGM $. Scattering enables us to measure $ IGM $ without assuming DM$_ h $ to constrain H$_0$. The reduction in systematic error is comparable to the Hubble tension ($ Combined with the fact that more localized FRBs will become available, our result indicates that our method can be used to address the Hubble tension using future FRB samples.

The Composite Spectral Energy Distribution of Quasars Is Surprisingly Universal Since Cosmic Noon

Leveraging the photometric data of the Sloan Digital Sky Survey and the Galaxy Evolution Explorer (GALEX), we construct mean/median spectral energy distributions (SEDs) for unique bright quasars in redshift bins of 0.2 and up to z≃3, after taking the GALEX non-detection into account. Further correcting for the absorption of the intergalactic medium, these mean/median quasar SEDs constitute a surprisingly redshift-independent mean/median composite SED from the rest-frame optical down to ≃500 A˚ for quasars with bolometric luminosity brighter than 1045.5ergs−1. Moreover, the mean/median composite quasar SED is plausibly also independent of black hole mass and Eddington ratio, and suggests similar properties of dust and gas in the quasar host galaxies since cosmic noon. Both the mean and median composite SEDs are nicely consistent with previous mean composite quasar spectra at wavelengths beyond ≃1000 A˚, but at shorter wavelengths, are redder, indicating, on average, less ionizing radiation than previously expected. Through comparing the model-predicted to the observed composite quasar SEDs, we favor a simply truncated disk model, rather than a standard thin disk model, for the quasar central engine, though we request more sophisticated disk models. Future deep ultraviolet facilities, such as the China Space Station Telescope and the Ultraviolet Explorer, would prompt revolutions in many aspects, including the quasar central engine, production of the broad emission lines in quasars, and cosmic reionization.

On the Andromeda-Milky Way Future Encounter: Thrice Faster Over Time

This inquiry accurately provides an analytical solution of the orbital free-falling time for the system of galaxies Andromeda-Milky Way within a non-expanding tired-light-dominated framework. The mathematical and physical background involves a two-body problem based on the orbital dynamics in which, in the first scenario, the gravitational interaction between local galaxies dominates over the expanding space. The latter is the standard astrophysical approach undertaken by public fund-based universities all over the world. However, in the second alternative scenario, Andromeda’s blueshift has to be corrected for the apparent recession velocity provided by the photon energy loss as a result of multiple interactions between photons and crystallized electrons through the intergalactic medium provided by a Tired Light, specifically, a New Tired Light process. This leads towards a substantial temporal deviation between the classic research approach and this unconventional and independent mindset. Hence, the outcome gives out a disagreement consisting in an orbital free-falling time equal to 1.7 rather than 4.5 billion years. Accordingly, the encounter of the two main galaxies of the Local Group will occur much earlier than estimated thus far. Accordingly, this implies an urgent shift in the scientific mindset as well as a course change in the application of the boundary conditions into forthcoming computational methods.

The cooler past of the intracluster medium in TNG-Cluster

Abstract The intracluster medium (ICM) today is comprised largely of hot gas with clouds of cooler gas of unknown origin and lifespan. We analyze the evolution of cool gas (temperatures ≲ 104.5 K) in the ICM of 352 galaxy clusters from the TNG-Cluster simulations, with present-day mass ∼1014.3 − 15.4 M⊙. We follow the main progenitors of these clusters over the past ∼13 billion years (since z ≲ 7) and find that, according to TNG-Cluster, the cool ICM mass increases with redshift at fixed cluster mass, implying that this cooler past of the ICM is due to more than just halo growth. The cool cluster gas at z ≲ 0.5 is mostly located in and around satellite galaxies, while at z ≳ 2 cool gas can also accrete via filaments from the intergalactic medium. Lower-mass and higher-redshift clusters are more susceptible to cooling. The cool ICM mass correlates with the number of gaseous satellites and inversely with the central supermassive black hole (SMBH) mass. The average number of gaseous satellites decreases since z = 2, correlating with the decline in the cool ICM mass over cosmic time, suggesting a link between the two. Concurrently, kinetic SMBH feedback shifts the ICM temperature distribution, decreasing the cool ICM mass inside-out. At z ≈ 0.5, the predicted Mg ii column densities are in the ballpark of recent observations, where satellites and other halos contribute significantly to the total Mg ii column density. Suggestively, a non-negligible amount of the ICM cool gas forms stars in-situ at early times, reaching ∼102 M⊙ yr−1 and an Hα surface brightness of ∼10−17 erg s−1 cm−2 arcsec−2 at z ≈ 2, detectable with Euclid and JWST.

Fast radio bursts as standard candles for cosmology

Recently, fast radio bursts (FRBs) have become a thriving field in astronomy and cosmology. Due to their extragalactic and cosmological origin, they are useful to study the cosmic expansion and the intergalactic medium (IGM). In the literature, the dispersion measure (DM) of FRB has been considered extensively. It could be used as an indirect proxy of the luminosity distance dL of FRB. The observed DM contains the contributions from the Milky Way (MW), the MW halo, IGM, and the host galaxy. Unfortunately, IGM and the host galaxy of FRB are poorly known to date, and hence the large uncertainties of DMIGM and DMhost in DM plague the FRB cosmology. Could we avoid DM in studying cosmology? Could we instead consider the luminosity distance dL directly in the FRB cosmology? We are interested to find a way out for this problem in the present work. From the lessons of calibrating type Ia supernovae (SNIa) or long gamma-ray bursts (GRBs) as standard candles, we consider a universal subclassification scheme for FRBs, and there are some empirical relations for them. In the present work, we propose to calibrate type Ib FRBs as standard candles by using a tight empirical relation without DM. The calibrated type Ib FRBs at high redshifts can be used like SNIa to constrain the cosmological models. We also test the key factors affecting the calibration and the cosmological constraints. Notice that the results are speculative, since they are based on simulations for the future (instead of now). We stress that the point in the present work is just to propose a viable method and a workable pipeline using type Ib FRBs as standard candles for cosmology in the future.

Measurements of the thermal and ionization state of the intergalactic medium during the cosmic afternoon

Abstract We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 Lyα absorption lines fitted from 12 archival HST STIS quasar spectra. We employ the machine-learning-based inference method that uses joint Doppler parameter - column density ($b-N_{{\rm {H\,{\small I}}}{}}$) distributions obtained from Lyα forest decomposition. Our results show that the H i photoionization rates, $\Gamma _{{\rm{H\,{\small I}}}{}}$, agree with recent UV background synthesis models, with $\log (\Gamma _{{\rm{H\,{\small I}}}}/\text{s}^{-1})={-11.79}^{+0.18}_{-0.15}$, ${-11.98}^{+0.09}_{-0.09}$, and ${-12.32}^{+0.10}_{-0.12}$ at z = 1.4, 1.2, and 1 respectively. We obtain the IGM temperature at the mean density, T0, and the adiabatic index, γ, as [log (T0/K), γ] = [${4.13}^{+0.12}_{-0.10}$, ${1.34}^{+0.10}_{-0.15}$], $[{3.79}^{+0.11}_{-0.11}$, ${1.70}^{+0.09}_{-0.09}]$ and $[{4.12}^{+0.15}_{-0.25}$, ${1.34}^{+0.21}_{-0.26}]$ at z = 1.4, 1.2 and 1. Our measurements of T0 at z = 1.4 and 1.2 are consistent with the trend predicted from previous z < 3 temperature measurements and theoretical expectations, where the IGM cools down after He ii reionization in the absence of any non-standard heating. However, our T0 measurement at z = 1 shows unexpectedly high IGM temperature. Given the relatively large uncertainty in these measurements, where $\sigma _{T_0} \sim 5000$ K, mostly emanating from the limited size of our dataset, we can not conclude whether the IGM cools down as expected. Lastly, we generate mock datasets to test the constraining power of future measurement with larger datasets. The results demonstrate that, with redshift pathlength Δz ∼ 2 for each redshift bin, three times the current dataset, we can constrain the T0 of IGM within 1500K, which would be sufficient to constrain the IGM thermal history at z < 1.5 conclusively.

A SPectroscopic Survey of Biased Halos In the Reionization Era (ASPIRE): JWST Supports Earlier Reionization around [O iii] Emitters

Understanding when and how reionization happened is crucial for studying the early structure formation and the properties of the first galaxies in the Universe. At z > 5.5, the observed intergalactic medium (IGM) optical depth shows a significant scatter, indicating an inhomogeneous reionization process. However, the nature of the inhomogeneous reionization remains debated. A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) is a JWST Cycle 1 program that has spectroscopically identified >400 [O iii] emitters in 25 quasar fields at z > 6.5. Combined with deep ground-based optical spectroscopy of ASPIRE quasars, the ASPIRE program provides the current largest sample for IGM-galaxy connection studies during cosmic reionization. We present the first results of IGM effective optical depth measurements around [O iii] emitters using 14 ASPIRE quasar fields. We find the IGM transmission is tightly related to reionization era galaxies to the extent that a significant excess of Lyα transmission exists around [O iii] emitters. We measure the stacked IGM effective optical depth of IGM patches associated with [O iii] emitters and find they reach the same IGM effective optical depth at least d z ∼ 0.1 ahead of those IGM patches where no [O iii] emitters are detected, supporting earlier reionization around [O iii] emitters. Our results indicate an enhancement in IGM Lyα transmission around [O iii] emitters at scales beyond 25 h −1 cMpc, consistent with the predicted topology of reionization from fluctuating UV background models.

The Cosmic Ultraviolet Baryon Survey (CUBS). IX. The Enriched Circumgalactic and Intergalactic Medium Around Star-forming Field Dwarf Galaxies Traced by O vi Absorption

The shallow potential wells of star-forming dwarf galaxies make their surrounding circumgalactic and intergalactic medium (CGM/IGM) sensitive laboratories for studying the inflows and outflows thought to regulate galaxy evolution. We present new absorption-line measurements in quasar sight lines, probing within projected distances of <300 kpc from 91 star-forming field dwarf galaxies with a median stellar mass of logM⋆/M⊙≈8.3 at 0.077 < z < 0.73, from the Cosmic Ultraviolet Baryon Survey (CUBS). In this redshift range, the CUBS quasar spectra cover a suite of transitions including H i, low, and intermediate metal ions (e.g., C ii, Si ii, C iii, and Si iii), and highly ionized O vi. This CUBS-Dwarfs survey enables constraints with samples nine times larger than past dwarf CGM/IGM studies with similar ionic coverage. We find that low and intermediate ionization metal absorption is rare around dwarf galaxies, consistent with previous surveys of local dwarfs. In contrast, highly ionized O vi is commonly observed in sight lines that pass within the virial radius of a dwarf, and O vi detection rates are nonnegligible at projected distances of 1−2× the virial radius. Based on these measurements, we estimate that the O vi-bearing phase of the CGM/IGM accounts for a dominant share of the metal budget of dwarf galaxies. The absorption kinematics suggest that a relatively modest fraction of the O vi-bearing gas is formally unbound. Together, these results imply that low-mass systems at z ≲ 1 effectively retain a substantial fraction of their metals within the nearby CGM and IGM.

The COS CGM Compendium. V. The Dichotomy of O vi Associated with Low- and High-metallicity Cool Gas at z < 1

We analyze the O vi content and kinematics for 126 H i-selected absorbers at 0.14 ≲ z ≲ 0.73 for which the metallicities of their cool photoionized phase have been determined. We separate the absorbers into 100 strong Lyα forest systems (SLFSs with 15 ≲ log N(H i) < 16.2) and 26 partial Lyman Limit systems (pLLSs with 16.2 ≤ log N(H i) ≤ 17.2). The sample is drawn from the COS CGM Compendium (CCC) and has O vi coverage in signal-to-noise ratio ≥ 8 Hubble Space Telescope/COS G130M/G160M QSO spectra, yielding a 2σ completeness level of logN (O vi) ≥ 13.6. The O vi detection rates differ substantially between low-metallicity (LM; [X/H] ≤ −1.4) and high-metallicity (HM; [X/H] > −1.4) SLFSs, with 20% and 60% detection rates, respectively. The O vi detection frequency for the HM and LM pLLSs is, however, similar at ∼60%. The SLFSs and pLLSs without detected O vi are consistent with the absorbing gas being in a single phase, while those with O vi trace multiphase gas. We show that the O vi velocity widths and column densities have different distributions in LM and HM gas. We find a strong correlation between O vi column density and metallicity. The strongest ( logN (O vi) ≳ 14) and broadest O vi absorbers are nearly always associated with HM absorbers, while weaker O vi absorbers are found in both LM and HM absorbers. From comparisons with galaxy-selected and blind O vi surveys, we conclude absorbers with logN (O vi) ≳ 14 most likely arise in the CGM of star-forming galaxies. Absorbers with weak O vi likely trace the extended CGM or intergalactic medium, while those without O vi likely originate in the intergalactic medium.

Osaka Feedback Model. III. Cosmological Simulation CROCODILE

We introduce our new cosmological simulation data set CROCODILE, executed using the GADGET4-Osaka smoothed particle hydrodynamics code. This simulation incorporates an updated supernova (SN) feedback model of Y. Oku et al. and an active galactic nuclei (AGN) feedback model. A key innovation in our SN feedback model is the integration of a metallicity- and redshift-dependent, top-heavy initial mass function. Our SN model introduces a new consideration that results in an order of magnitude difference in the energy injection rate per unit stellar mass formed at high redshift. The CROCODILE data set is comprehensive, encompassing a variety of runs with diverse feedback parameters. This allows for an in-depth exploration of the relative impacts of different feedback processes in galactic evolution. Our initial comparisons with observational data, spanning the galaxy stellar mass function, the star formation main sequence, and the mass–metallicity relation, show promising agreement, especially for the Fiducial run. These results establish a solid foundation for our future work. We find that SN feedback is a key driver in the chemical enrichment of the intergalactic medium (IGM). Additionally, the AGN feedback creates metal-rich, bipolar outflows that extend and enrich the circumgalactic medium and IGM over a few Mpc scales.

QSO MUSEUM

Extended Lyα emission is routinely found around single quasars across cosmic time. However, few studies have investigated how such emission changes in fields with physically associated quasar pairs, which should reside in dense environments and are predicted to be linked through intergalactic filaments. We present VLT/MUSE snapshot observations (45 minutes/source) to unveil extended Lyα emission on scales of the circumgalactic medium (CGM) around the largest sample of physically associated quasar pairs to date, encompassing eight pairs (14 observed quasars) at z ∼ 3 with an i-band magnitude between 18 and 22.75, corresponding to absolute magnitudes Mi(z = 2) between −29.6 and −24.9. The pairs are either at close (∼50–100 kpc, five pairs) or wide (∼450–500 kpc, three pairs) angular separation and have velocity differences of Δv ≤ 2000 km s−1. We detected extended emission around 12 of the 14 targeted quasars and investigated the luminosity, size, kinematics, and morphology of these Lyα nebulae. On average, they span about 90 kpc and are 2.8 × 1043 erg s−1 bright. Irrespective of the quasars’ projected distance, the nebulae often (∼45%) extend toward the other quasar in the pair, leading to asymmetric emission whose flux-weighted centroid is at an offset position from any quasar location. We show that large nebulae are preferentially aligned with the large-scale structure, as traced by the direction between the two quasars, and conclude that the cool gas (104 K) in the CGM traces well the direction of cosmic web filaments. Additionally, the radial profile of the Lyα surface brightness around quasar pairs can be described by a power law with a shallower slope (∼−1.6) with respect to single quasars (∼−2), indicative of increased CGM densities out to large radii and/or an enhanced contribution from the intergalactic medium (IGM) due to the dense environments expected around quasar pairs. The sample presented in this study contains excellent targets for ultra-deep observations to directly study filamentary IGM structures in emission. This work demonstrates that a large snapshot survey of quasar pairs will pave the way to direct statistical study of the IGM.

The twisted jets and magnetic fields of the extended radio galaxy 4C 70.19

Context. The appearance of the jets and lobes of some radio galaxies makes it difficult to assign them to a known class of objects. This is often due to the activity of the central engine and/or interactions with the environment, as well as projection effects. Aims. We analyse the radio data for an apparently asymmetric radio source 4C 70.19, which is associated with the giant elliptical galaxy NGC 6048. The source shows distorted radio jets and lobes, one of which bends by 180°. The aim of our study is to explain the nature of the observed distortions. Methods. We used LOFAR, Effelsberg, and VLA radio data in a wide range of frequencies. At high frequencies, we also used radio polarimetry to study the properties of the magnetic fields. Additionally, we made use of optical, infra-red, and X-ray data. Results. Polarisation data suggest shearing of the magnetic fields at points where the jets bend. The low-frequency LOFAR map at 145 MHz, as well as the sensitive single-dish Effelsberg map at 8.35 GHz, reveal previously undetected diffuse emission around the source. The rotation measure (RM) derived from the polarimetric data allowed us to estimate the density of the medium surrounding the source, which agrees with typical densities of the intergalactic medium or the outer parts of insterstellar halos. Conclusions. We propose that the southern jet is bent in the same manner as the northern one, but that it is inclined to the sky plane. Both these bends are likely caused by the orbital motion within the galaxy group, as well as interactions with the intergalactic medium. Our analyses suggest that, despite its complex morphology, 4C 70.19 seems to be intrinsically symmetric with a physical extent of up to 600 kpc, and that the diffuse emission detected in our high-sensitivity maps is related to radio plumes that are expanding behind the source.

Dark photon constraints from CMB temperature anisotropies

The resonant conversion, within the inter-galactic medium, of regular photons into dark photons amplifies the anisotropy observed in the CMB, thereby imposing stringent constraints on the existence of light dark photons. In this study, we investigate the impact of light dark photons, with masses in the range 3 × 10-15 eV < mA ' < 3 × 10-12 eV on the power spectrum of temperature anisotropies within the cosmic microwave background (CMB) radiation utilizing the state-of-the-art large-volume FLAMINGO cosmological simulations. Our results show that using full Planck data, one can expect the existing constraints on the dark photon mixing parameter in this mass range to improve by an order of magnitude.

Spectral reconstruction for radiation hydrodynamic simulations of galaxy evolution

Radiation from stars and active galactic nuclei (AGN) plays an important role in galaxy formation and evolution, and profoundly transforms the intergalactic, circumgalactic, and interstellar medium (IGM, CGM, and ISM). On-the-fly radiative transfer (RT) has started being incorporated in cosmological simulations, but the complex evolving radiation spectra are often crudely approximated with a small number of broad bands with piece-wise constant intensity and a fixed photo-ionisation cross-section. Such a treatment is unable to capture the changes to the spectrum as light is absorbed while it propagates through a medium with non-zero opacity. This can lead to large errors in photo-ionisation and heating rates. In this work we present a novel approach of discretising the radiation field at discrete photon energies, at the edges of the typically used photo-ionising bands, in order to capture the power-law slope of the radiation field. In combination with power-law approximations for the photo-ionisation cross-sections, this model allows us to self-consistently combine radiation from sources with different spectra and accurately follow the ionisation states of primordial and metal species through time. The method is implemented in GASOLINE2 in connection with TREVR2, a fast reverse ray tracing algorithm with 𝒪(Nactive log2 N) scaling. We compare our new piece-wise power-law reconstruction to the piece-wise constant method in calculating the primordial chemistry photo-ionisation and heating rates under an evolving UV background (UVB) and stellar spectrum, and find that our method reduces errors significantly, by up to two orders of magnitude in the case of HeII ionisation. We apply our new spectral reconstruction method in RT post-processing of a cosmological zoom-in simulation, MUGS2 g1536, including radiation from stars and a live UVB, and find a significant increase in total neutral hydrogen (HI) mass in the ISM and the CGM due to shielding of the UVB and a low escape fraction of the stellar radiation. This demonstrates the importance of RT and an accurate spectral approximation in simulating the CGM-galaxy ecosystem.

The MUSE eXtremely Deep Field: Detections of circumgalactic Si II* emission at z ≳ 2

Context. The circumgalactic medium (CGM) serves as a baryon reservoir that connects galaxies to the intergalactic medium and fuels star formation. The spatial distribution of the metal-enriched cool CGM has not yet been directly revealed at cosmic noon (z ≃ 2–4), as bright emission lines at these redshifts are not covered by optical integral field units. Aims. To remedy this situation, we performed the first-ever detections and exploration of extended Si II* emission in the low-ionization state (LIS), referred to as Si II* halos, at redshifts ranging from z = 2 to 4 as a way to trace the metal-enriched cool CGM. Methods. We used a sample of 39 galaxies with systemic redshifts of z = 2.1–3.9 measured with the [C III] doublet in the MUSE Hubble Ultra Deep Field catalog, whose integration times span from ≃30 to 140 hours. We searched for extended Si II* λ1265, 1309, 1533 emission (fluorescent lines) around individual galaxies. We also stacked a subsample of 14 UV-bright galaxies. Results. We report five individual detections of Si II* λ1533 halos. We also confirm the presence of Si II* λ1533 halos in stacks for the subsample containing UV-bright sources. The other lines do not show secure detections of extended emission in individual or in stacking analyses. These detections may imply that the presence of metal-enriched CGM is a common characteristic for UV-bright galaxies. To investigate whether the origin of Si II* is continuum pumping, as suggested in previous studies, we checked the consistency of the equivalent width (EW) of Si II* emission and the EW of Si II absorption for the individual halo object with the most reliable detection. We confirm the equivalence, suggesting that photon conservation works for this object and points toward continuum pumping as the source of Si II*. We also investigated Si II* lines in a RAMSES-RT zoom-in simulation including continuum pumping, and find the ubiquitous presence of extended halos.

JWST/NIRSpec Observations of Lyman α Emission in Star-forming Galaxies at 6.5 ≲ z ≲ 13

We present an analysis of JWST Lyα spectroscopy of z ≳ 6.5 galaxies, using observations in the public archive covering galaxies in four independent fields: Great Observatories Origins Deep Survey (GOODS)-N, GOODS-S, A2744, and the Extended Groth Strip (EGS). We measure the Lyα emission line properties for a sample of 210 z ≃ 6.5–13 galaxies, with redshifts confirmed independently of Lyα in all cases. We present three new detections of Lyα emission in JWST spectra, including a large equivalent width (EW; =143 Å) Lyα emitter (LAE) with strong C iv emission (EW = 21 Å) at z = 7.1 in GOODS-N. We measure the redshift-dependent Lyα EW distribution across our sample. We find that strong Lyα emission (EW > 25 Å) becomes increasingly rare at earlier epochs, suggesting that the transmission of Lyα photons decreases by 4× between z ≃ 5 and z ≃ 9. We describe potential implications for the intergalactic medium neutral fraction. There is significant field-to-field variance in the LAE fraction. In contrast to the three other fields, the EGS shows no evidence for reduced transmission of Lyα photons at z ≃ 7–8, suggesting a significantly ionized sight line may be present in the field. We use available NIRCam grism observations from the First Reionization Epoch Spectroscopically Complete Observations survey to characterize overdensities on large scales around known LAEs in the GOODS fields. The strongest overdensities appear linked with extremely strong Lyα detections (EW > 50 Å) in most cases. Future Lyα spectroscopy with JWST has the potential to constrain the size of ionized regions around early galaxy overdensities, providing a new probe of the reionization process.

Evidence of Extreme Ionization Conditions and Low Metallicity in GHZ2/GLASS-Z12 from a Combined Analysis of NIRSpec and MIRI Observations

GHZ2/GLASS-z12, one of the most distant galaxies found in JWST observations, has been recently observed with both the NIRSpec and MIRI spectrographs, establishing a spectroscopic redshift z spec = 12.34 and making it the first system at z > 10 with complete spectroscopic coverage from rest-frame UV to optical wavelengths. This galaxy is identified as a strong C iv λ1549 emitter (EW = 46 Å) with many other detected emission lines, such as N iv] λ1488, He ii λ1640, O iii] λ λ1661,1666, N iii] λ1750, C iii] λ λ1907,1909, [O ii] λ λ3726,3729, [Ne iii] λ3869, [O iii] λ λ4959,5007, and Hα, including a remarkable detection of the O iii Bowen fluorescence line at rest frame λ = 3133 Å. We analyze in this paper the joint NIRSpec + MIRI spectral data set. Combining six optical strong-line diagnostics (namely R2, R3, R23, O32, Ne3O2, and Ne3O2Hd), we find extreme-ionization conditions, with log10 ([O III] λ λ4959,5007/[O II] λ λ3726,3729) = 1.39 ± 0.19 and log10 ([Ne III] λ3869/[O II] λ λ3726,3729) = 0.37 ± 0.18 in stark excess compared to typical values in the interstellar medium (ISM) at lower redshifts. These line properties are compatible either with an active galactic nucleus (AGN) or with a compact, very dense star-forming environment (ΣSFR ≃ 102–103 M ⊙ yr−1 kpc−2 and ΣM* ≃ 104–105 M ⊙ pc−2), with a high ionization parameter (log10(U) =−1.75 ± 0.16), a high ionizing photon production efficiency log(ξion)=25.7−0.1+0.3 , and a low gas-phase metallicity (also confirmed by the direct, T e method) ranging between 4% and 11% Z ⊙, indicating a rapid chemical enrichment of the ISM in the past few megayears. These properties also suggest that a substantial amount of ionizing photons (∼10%) are leaking outside of GHZ2 and starting to reionize the surrounding intergalactic medium, possibly due to strong radiation-driven winds. The general lessons learned from GHZ2 are the following: (i) the UV-to-optical combined nebular indicators are broadly in agreement with UV-only or optical-only indicators; (ii) UV+optical diagnostics fail to discriminate between an AGN and star formation in a low-metallicity, high-density, and extreme-ionization environment; and (iii) comparing the nebular line ratios with local analogs may be approaching its limits at z ≳ 10, as this approach is potentially challenged by the unique conditions of star formation experienced by galaxies at these extreme redshifts.