Search

Seismic traveltime tomography represents a popular and useful tool for unravelling the structure of the subsurface across the scales. In this work we address the case where the forward model is represented by the eikonal equation and derive a formalism to solve the inverse problem where gradients are calculated efficiently using the discrete adjoint state method. Our approach provides gradients with respect to both velocity structure and source locations, allowing us to perform a consistent joint inversion. The forward problem is solved using a second-order fast-marching method, which provides a strategy to efficiently solve the adjoint problem. Our approach allows for arbitrary positions of both sources and receivers and for a refined grid around the source region to reduce errors in computed traveltimes. We show how gradients computed using the discrete adjoint method can be employed to perform either deterministic inversion, i.e., solving an optimization problem, or for a probabilistic (Bayesian) approach, i.e., obtaining a posterior probability density function. We show applications of our methodology on a set of synthetic examples both in 2D and 3D using the L-BFGS algorithm for the deterministic case and the Hamiltonian Monte Carlo algorithm for the probabilistic case., Comment: 28 pages, 7 figures

Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 250 RS (Schwarzschild radii) perpendicular to the jet axes. This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 300 RS and 10000 RS distances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 {\mu}as from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at about 130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 26000 Gauss, which is consistent with previous measurements. The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength., Comment: 22 pages, 10 figures, published in A&A

The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-K\"onigl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations.

We introduce FACTS Grounding, an online leaderboard and associated benchmark that evaluates language models' ability to generate text that is factually accurate with respect to given context in the user prompt. In our benchmark, each prompt includes a user request and a full document, with a maximum length of 32k tokens, requiring long-form responses. The long-form responses are required to be fully grounded in the provided context document while fulfilling the user request. Models are evaluated using automated judge models in two phases: (1) responses are disqualified if they do not fulfill the user request; (2) they are judged as accurate if the response is fully grounded in the provided document. The automated judge models were comprehensively evaluated against a held-out test-set to pick the best prompt template, and the final factuality score is an aggregate of multiple judge models to mitigate evaluation bias. The FACTS Grounding leaderboard will be actively maintained over time, and contains both public and private splits to allow for external participation while guarding the integrity of the leaderboard. It can be found at https://www.kaggle.com/facts-leaderboard.

Past millimeter-wave galaxy surveys have probed the brightest starburst galaxies only and suffered heavily from confusion. The interpretation of existing surveys has also been hindered by the lack of reliable redshift indicators for measuring distances for the entire sample. Thanks to recent advances in mm-wave detector technologies we can now overcome these limitations, and conduct the first truly volumetric surveys of star-forming galaxies at mm-wavelengths down to the L* luminosities of typical galaxies, with ~1000 redshift slices spanning most of the Cosmic star-forming volume (z ~ 1--12) with nearly uniform mass and luminosity selection. We describe an instrument concept capable of delivering such surveys with the technologies available today, which can be built and operated on a ground-based mm-wave facility in the near future. Such spectrometer cameras can resolve and redshift identify up to to 25,000 star-forming galaxies per year even when operated on a 10-m class telescope. On a larger aperture it can do the same faster or probe even deeper. We propose a loose, open-source collaboration to design, build, and operate one or several such cameras through the shared contributions of leading experts and telescopes from around the globe., Comment: 10 pages, 1 figure, SPIE Conference 13101

Using the Event Horizon Telescope (EHT), the gravitationally lensed rings around the supermassive black holes (SMBHs) in Messier 87 (M87) and Sagittarius A* (Sgr A*) have now been successfully imaged at a resolution under 10 gravitational radii (R$_{\rm g}$ $ = \rm{GM/c^2}$). To expand studies beyond M87 and Sgr A*, we have constructed the Event Horizon and Environs (ETHER) sample, a comprehensive database encompassing approximately 3.15 million SMBH mass estimates, $\sim$ 20,000 Very-Long Baseline Interferometry (VLBI) radio flux densities, and $\sim$ 36,000 hard X-ray flux densities. This database is designed to identify and optimize target selection for the EHT and its upgrades on the ground and in space. We have identified a Gold Sample (GS) of nearby low-luminosity Active Galactic Nuclei (AGNs) within it that are ideal for studying jet bases and potentially imaging black hole shadows. We observed 27 of these AGNs using the EHT from 2022 to 2024, providing an opportunity to resolve and image accretion flows and jets at resolutions of $\leq$ 100 R$_{\rm g}$. Only a few SMBHs have sufficiently high enough flux density to be imaged at scales of $\leq$ 50 R$_{\rm g}$ with the present EHT. Among these are M87, Sgr A*, NGC4594 (Sombrero/M104), NGC4261, and NGC4374 (Messier 84/M84). Of these, NGC4261, Sombrero, and M84 have been observed and/or are scheduled for deep imaging with EHT+ALMA from 2023 to 2025. Sombrero, NGC4261, M84, NGC4278, and NGC5232 are clearly detected in our EHT+ALMA observations in 2022, indicating that the 230 GHz flux density from the accretion flows is significantly high. Ongoing imaging of the ETHER GS will enable measurements of black hole mass and spin, help constrain General Relativity, and enrich our understanding of jet launching and accretion inflows across a broad multi-parameter space, including black hole mass, spin, accretion rate, and orientation., Comment: 9 pages, 6 figures, 1 table, published in Proceedings of the 16th EVN Symposium, Ed. E. Ros, P. Benke, S.A. Dzib, I. Rottmann, & J.A. Zensus, Bonn: Max-Planck-Institut f\"ur Radioastronomie, 2024, pages 75-84, https://cloud.mpifr-bonn.mpg.de/index.php/s/BkX2CC2Xjn2aKR4

The visibility of the Lyman-$\alpha$ (Ly$\alpha$) emission from reionization-epoch galaxies depends sensitively on the extent of the intrinsic \lya emission redwards of 1215.67~\AA. The prominent red peak resulting from resonant radiative transfer in the interstellar medium is often modelled as a single Gaussian. We use the \textsc{Azahar} simulation suite of a massive-reionization epoch galaxy to show that a significantly larger fraction of the \lya emission extends to $400$-$800$~km~s$^{-1}$, and thus significantly further to the red than predicted by a Gaussian line profile. A cycle of frequent galaxy mergers strongly modulates the \lya luminosity, the red peak velocity and its extended red wing emerging from the galaxy, which all also strongly vary with viewing angle. The \lya emission also depends sensitively on the implemented feedback, dust and star formation physics. Our simulations including cosmic rays reproduce the observed spectral properties of reionization epoch \lya emitters (LAEs) well if we assume that the \lya emission is affected by very little dust. The visibility of LAEs can be strongly underestimated if the extended red wings of the intrinsic \lya emission are not accounted for. We discuss implications for using the visibility of LAEs to constrain the evolution of the volume-averaged neutral fraction during reionization.

Recent JWST observations have allowed for the first time to obtain comprehensive measurements of the ionizing photon production efficiency $\xi_\text{ion} $ for a wide range of reionization-epoch galaxies. We explore implications for the inferred UV luminosity functions and escape fractions of ionizing sources in our suite of simulations. These are run with the GPU-based radiative transfer code ATON-HE and are calibrated to the XQR-30 Lyman-alpha forest data at $550$% at $z> 10$, disfavouring the oligarchic source model at very high redshift. The inferred effective clumping factors in our simulations are in the range of $3-6$, suggesting consistency between the observed ionizing properties of reionization-epoch galaxies and the ionizing photon budget in our simulations., Comment: 6 pages. 3 figures

The Lyman alpha (Lya) forest in the spectra of z>5 quasars provides a powerful probe of the late stages of the Epoch of Reionization (EoR). With the recent advent of exquisite datasets such as XQR-30, many models have struggled to reproduce the observed large-scale fluctuations in the Lya opacity. Here we introduce a Bayesian analysis framework that forward-models large-scale lightcones of IGM properties, and accounts for unresolved sub-structure in the Lya opacity by calibrating to higher-resolution hydrodynamic simulations. Our models directly connect physically-intuitive galaxy properties with the corresponding IGM evolution, without having to tune "effective" parameters or calibrate out the mean transmission. The forest data, in combination with UV luminosity functions and the CMB optical depth, are able to constrain global IGM properties at percent level precision in our fiducial model. Unlike many other works, we recover the forest observations without evoking a rapid drop in the ionizing emissivity from z~7 to 5.5, which we attribute to our sub-grid model for recombinations. In this fiducial model, reionization ends at $z=5.44\pm0.02$ and the EoR mid-point is at $z=7.7\pm0.1$. The ionizing escape fraction increases towards faint galaxies, showing a mild redshift evolution at fixed UV magnitude, Muv. Half of the ionizing photons are provided by galaxies fainter than Muv~-12, well below direct detection limits of optical/NIR instruments including JWST. We also show results from an alternative galaxy model that does not allow for a redshift evolution in the ionizing escape fraction. Despite being decisively disfavored by the Bayesian evidence, the posterior of this model is in qualitative agreement with that from our fiducial model. We caution however that our conclusions regarding the early stages of the EoR and which sources reionized the Universe are more model-dependent., Comment: 21 pages, 12 figures, 3 tables (submitted to PASA)

The X-ray polarization observations made possible with the Imaging X-ray Polarimetry Explorer (IXPE) offer new ways of probing high-energy emission processes in astrophysical jets from blazars. Here we report on the first X-ray polarization observation of the blazar S4 0954+65 in a high optical and X-ray state. During our multi-wavelength campaign on the source, we detected an optical flare whose peak coincided with the peak of an X-ray flare. This optical-X-ray flare most likely took place in a feature moving along the parsec-scale jet, imaged at 43 GHz by the Very Long Baseline Array. The 43 GHz polarization angle of the moving component underwent a rotation near the time of the flare. In the optical band, prior to the IXPE observation, we measured the polarization angle to be aligned with the jet axis. In contrast, during the optical flare the optical polarization angle was perpendicular to the jet axis; after the flare, it reverted to being parallel to the jet axis. Due to the smooth behavior of the optical polarization angle during the flare, we favor shocks as the main acceleration mechanism. We also infer that the ambient magnetic field lines in the jet were parallel to the jet position angle. The average degree of optical polarization during the IXPE observation was (14.3$\pm$4.1)%. Despite the flare, we only detected an upper limit of 14% (at 3$\sigma$ level) on the X-ray polarization degree; although a reasonable assumption on the X-ray polarization angle results in an upper limit of 8.8% ($3\sigma$). We model the spectral energy distribution (SED) and spectral polarization distribution (SPD) of S4 0954+65 with leptonic (synchrotron self-Compton) and hadronic (proton and pair synchrotron) models. The constraints we obtain with our combined multi-wavelength polarization observations and SED modeling tentatively disfavor hadronic models for the X-ray emission in S4 0954+65., Comment: Submitted to A&A, 16 pages, 5 figures, and 7 tables

It is well known that the edge limit of Gaussian/Laguerre Beta ensembles is given by Airy($\beta$) point process. We prove an universality result that this also holds for a general class of additions of Gaussian and Laguerre ensembles. In order to make sense of $\beta$-addition, we introduce type A Bessel function as the characteristic function of our matrix ensemble following the line of Gorin-Marcus, Benaych Georges-Cuenca-Gorin. Then we extract its moment information through the action of Dunkl operators, a class of differential operators originated from special function theory. We do the action explicitly on the Bessel generating functions of our additions, and after the asymptotic analysis, we obtain certain limiting functional in terms of conditional Brownian bridges of the Laplace transform of Airy($\beta$), which is universal up to proper rescaling among all our additions., Comment: 69 pages, more than 10 figures;v2: fixed some minor issues

We present images of IRAS~08235--4316 with the Dark Energy Camera Plane Survey (DECaPS; spanning 0.398--1.034$\,\mu$m, at ${\sim}1''$ resolution) and the Submillimeter Array (SMA; at 1.38\,mm/217\,GHz, at ${\sim}1.9''\times1.2''$ resolution), a YSO located in the Vela constellation near to the Puppis boundary, detected in a systematic search for new large/extended emission sources. The DECaPS data show an asymmetric bi--polar morphology with a large angular extent of ${\sim}7.1''$ separated by a dark lane, characteristic of highly--inclined protoplanetary disks and less-evolved YSOs with outflows. The SMA data show an extended continuum structure along the optical dark lane with a smaller angular extent of ${\sim}4.6''$. The detected $^{12}$CO J=2--1 emission tentatively shows a velocity gradient along the position angle of the dark lane/millimeter continuum, that may trace rotating gas. Additional $^{12}$CO emission is present which could trace infalling/outflowing gas, and/or a nearby gas cloud. We estimate a distance to IRAS~08235--4316 of at least ${\sim}191\,$pc. Supported by additional SED modelling, we infer IRAS~08235--4316 to be a newly discovered class~I YSO with an outflow, host to an embedded protoplanetary disk, with a large millimeter radius of ${\sim}440$\,au and dust mass ${\gtrsim}11\,M_\oplus$., Comment: Accepted for publication in AJ 14 pages (10 page paper), 6 figures

In the past few years, the improved sensitivity and cadence of wide-field optical surveys have enabled the discovery of several afterglows without associated detected gamma-ray bursts (GRBs). We present the identification, observations, and multiwavelength modeling of a recent such afterglow (AT2023lcr), and model three literature events (AT2020blt, AT2021any, and AT2021lfa) in a consistent fashion. For each event, we consider the following possibilities as to why a GRB was not observed: 1) the jet was off-axis; 2) the jet had a low initial Lorentz factor; and 3) the afterglow was the result of an on-axis classical GRB (on-axis jet with physical parameters typical of the GRB population), but the emission was undetected by gamma-ray satellites. We estimate all physical parameters using afterglowpy and Markov Chain Monte Carlo methods from emcee. We find that AT2023lcr, AT2020blt, and AT2021any are consistent with on-axis classical GRBs, and AT2021lfa is consistent with both on-axis low Lorentz factor ($\Gamma_0 \approx 5 - 13$) and off-axis ($\theta_\text{obs}=2\theta_\text{jet}$) high Lorentz factor ($\Gamma_0 \approx 100$) jets., Comment: 40 pages, 18 figures, 20 tables

Modern deep neural networks have been observed to exhibit a simple structure in their final layer features and weights, commonly referred to as neural collapse. This phenomenon has also been noted in layers beyond the final one, an extension known as deep neural collapse. Recent findings indicate that such a structure is generally not optimal in the deep unconstrained feature model, an approximation of an expressive network. This is attributed to a low-rank bias induced by regularization, which favors solutions with lower-rank than those typically associated with deep neural collapse. In this work, we extend these observations to the cross-entropy loss and analyze how the low-rank bias influences various solutions. Additionally, we explore how this bias induces specific structures in the singular values of the weights at global optima. Furthermore, we examine the loss surface of these models and provide evidence that the frequent observation of deep neural collapse in practice, despite its suboptimality, may result from its higher degeneracy on the loss surface., Comment: 40 pages

We study a flaring activity of the HSP Mrk421 that was characterized from radio to very-high-energy (VHE; E $>0.1$TeV) gamma rays with MAGIC, Fermi-LAT, Swift, XMM-Newton and several optical and radio telescopes. These observations included, for the first time for a gamma-ray flare of a blazar, simultaneous X-ray polarization measurements with IXPE. We find substantial variability in both X-rays and VHE gamma rays throughout the campaign, with the highest VHE flux above 0.2 TeV occurring during the IXPE observing window, and exceeding twice the flux of the Crab Nebula. However, the VHE and X-ray spectra are on average softer, and the correlation between these two bands weaker that those reported in previous flares of Mrk421. IXPE reveals an X-ray polarization degree significantly higher than that at radio and optical frequencies. The X-ray polarization angle varies by $\sim$100$^\circ$ on timescales of days, and the polarization degree changes by more than a factor 4. The highest X-ray polarization degree reaches 26%, around which a X-ray counter-clockwise hysteresis loop is measured with XMM-Newton. It suggests that the X-ray emission comes from particles close to the high-energy cutoff, hence possibly probing an extreme case of the Turbulent Extreme Multi-Zone model. We model the broadband emission with a simplified stratified jet model throughout the flare. The polarization measurements imply an electron distribution in the X-ray emitting region with a very high minimum Lorentz factor, which is expected in electron-ion plasma, as well as a variation of the emitting region size up to a factor of three during the flaring activity. We find no correlation between the fluxes and the evolution of the model parameters, which indicates a stochastic nature of the underlying physical mechanism. Such behaviour would be expected in a highly turbulent electron-ion plasma crossing a shock front., Comment: Accepted in Astronomy and Astrophysics (on January 17th 2025). Corresponding authors: Axel Arbet-Engels, Lea Heckmann, David Paneque

X-ray polarization is a unique new probe of the particle acceleration in astrophysical jets made possible through the Imaging X-ray Polarimetry Explorer. Here we report on the first dense X-ray polarization monitoring campaign on the blazar Mrk 421. Our observations were accompanied by an even denser radio and optical polarization campaign. We find significant short-timescale variability in both X-ray polarization degree and angle, including a $\sim90^\circ$ angle rotation about the jet axis. We attribute this to random variations of the magnetic field, consistent with the presence of turbulence but also unlikely to be explained by turbulence alone. At the same time, the degree of lower-energy polarization is significantly lower and shows no more than mild variability. Our campaign provides further evidence for a scenario in which energy-stratified shock-acceleration of relativistic electrons, combined with a turbulent magnetic field, is responsible for optical to X-ray synchrotron emission in blazar jets., Comment: 23 pages, including 8 pages of appendices. 12 figures, 3 tables. Submitted to ApJ

At millimeter wavelengths, the atmospheric emission is circularly polarized owing to the Zeeman splitting of molecular oxygen by the Earth's magnetic field. We report a measurement of the signal in the 150 GHz band using 3 years of observational data with the \textsc{Polarbear} project. Non-idealities of a continuously rotating half-wave plate (HWP) partially convert circularly polarized light to linearly polarized light. While \textsc{Polarbear} detectors are sensitive to linear polarization, this effect makes them sensitive to circular polarization. Although this was not the intended use, we utilized this conversion to measure circular polarization. We reconstruct the azimuthal gradient of the circular polarization signal and measure its dependency from the scanning direction and the detector bandpass. We compare the signal with a simulation based on atmospheric emission theory, the detector bandpass, and the HWP leakage spectrum model. We find the ratio of the observed azimuthal slope to the simulated slope is $0.92 \pm 0.01\rm{(stat)} \pm 0.07\rm{(sys)}$. This ratio corresponds to a brightness temperature of $3.8\,\mathrm{m K}$ at the effective band center of $121.8\,\mathrm{GHz}$ and bandwidth of $3.5\,\mathrm{GHz}$ estimated from representative detector bandpass and the spectrum of Zeeman emission. This result validates our understanding of the instrument and reinforces the feasibility of measuring the circular polarization using the imperfection of continuously rotating HWP. Continuously rotating HWP is popular in ongoing and future cosmic microwave background experiments to modulate the polarized signal. This work shows a method for signal extraction and leakage subtraction that can help measuring circular polarization in such experiments., Comment: 23 pages, 11 figures

The first very long baseline interferometry (VLBI) detections at 870$\mu$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$\lambda$ corresponding to an angular resolution, or fringe spacing, of 19$\mu$as. The Allan deviation of the visibility phase at 870$\mu$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$\mu$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time., Comment: Corresponding author: S. Doeleman

The relative contribution of emission from stellar sources and accretion onto supermassive black holes to reionization has been brought into focus again by the apparent high abundance of faint AGN at $4\lesssim z\lesssim11$ uncovered by JWST. We investigate the contribution of these faint AGN to hydrogen and the early stages of helium reionization using the GPU-based radiative transfer code ATON-HE by post-processing a cosmological hydrodynamical simulation from the SHERWOOD-RELICS suite of simulations. We study four models of reionization: two previously studied galaxy-only late-end reionization models and two new models -- a QSO-assisted model and a QSO-only model. In the QSO-assisted model, 1\% of the haloes host AGN with a 10~Myr lifetime, and the AGN luminosities are scaled such that the AGN contribution to the hydrogen-ionizing emissivity is 20\% of that contributed by galaxies. In the QSO-only model, quasars account for all the hydrogen-ionizing emissivity, with 10\% of the haloes hosting AGN, each with a 10 Myr lifetime. All models are calibrated to the observed mean Lyman-$\alpha$ forest transmission at $5\lesssim z\lesssim6.2$. We find that the QSO-assisted model requires an emissivity factor of $1.8$ lower than the galaxy-only models towards the end of reionization and fits the observed distribution of the Lyman-$\alpha$ optical depths well. Our QSO-only model is inconsistent with the observed Lyman-$\alpha$ optical depths distribution. It also results in too high IGM temperatures at $z\lesssim 5$ due to an early onset of HeII reionization unless the escape fraction of HeII-ionizing photons is assumed to be low. Our results suggest that a modest contribution to reionization by faint AGN is in good agreement with the Lyman-$\alpha$ forest data. In contrast, a scenario dominated by faint AGN appears difficult to reconcile with these observations., Comment: 19 pages, 11 figures

The J=1$\rightarrow$0 spectral line of carbon monoxide (CO(1-0)) is the canonical tracer of molecular gas. However, CO(2-1) is frequently used in its place, following the assumption that the higher energy line can be used to infer the CO(1-0) luminosity and molecular gas mass. The use of CO(2-1) depends on a knowledge of the ratio between CO(2-1) and CO(1-0) luminosities, r21. Here we present galaxy-integrated r21 measurements for 122 galaxies spanning stellar masses from 10$^9$ to 10$^{11.5}$ M$_\odot$ and star formation rates (SFRs) from 0.08 to 35 M$_\odot$/yr. We find strong trends between r21 and SFR, SFR surface density, star formation efficiency, and distance from the star formation main sequence (SFMS). We show that the assumption of a constant r21 can introduce biases into the molecular gas trends in galaxy population studies and demonstrate how this affects the recovery of important galaxy scaling relations, including the Kennicutt-Schmidt law and the relation between SFMS offset and star formation efficiency. We provide a prescription which accounts for variations in r21 as a function of SFR and can be used to convert between CO(2-1) and CO(1-0) when only one line is available. Our prescription matches variations in r21 for both AMISS and literature samples and can be used to derive more accurate gas masses from CO(2-1) observations., Comment: Accepted version, in press at ApJ

We study homological mirror symmetry for $(\mathbb{P}^2, \Omega)$ viewed as an object of birational geometry, with $\Omega$ the standard meromorphic volume form. First, we construct universal objects on the two sides of mirror symmetry, focusing on the exact symplectic setting: a smooth complex scheme $U_\mathrm{univ}$ and a Weinstein manifold $M_\mathrm{univ}$, both of infinite type; and we prove homological mirror symmetry for them. Second, we consider autoequivalences. We prove that automorphisms of $U_\mathrm{univ}$ are given by a natural discrete subgroup of $\operatorname{Bir} (\mathbb{P}^2, \pm \Omega)$; and that all of these automorphisms are mirror to symplectomorphisms of $M_\mathrm{univ}$. We conclude with some applications., Comment: 24 pages, 3 figures; comments welcome

The CO(1-0) line has been carefully calibrated as a tracer of molecular gas mass. However, recent studies often favor higher J transitions of the CO molecule which are brighter and accessible for redshift ranges where CO(1-0) is not. These lines are not perfect analogues for CO(1-0), owing to their more stringent excitation conditions, and must be calibrated for use as molecular gas tracers. Here we introduce the Arizona Molecular ISM Survey with the SMT (AMISS), a multi-CO line survey of z~0 galaxies conducted to calibrate the CO(2-1) and CO(3-2) lines. The final survey includes CO(2-1) spectra of 176 galaxies and CO(3-2) spectra for a subset of 45. We supplement these with archival CO(1-0) spectra from xCOLD GASS for all sources and additional CO(1-0) observations with the Kitt Peak 12m Telescope. Targets were selected to be representative of the galaxy population in the stellar mass range of $10^9$ to $10^{11.5}$ M$_\odot$. Our project emphasized careful characterization of statistical and systematic uncertainties to enable studies of trends in CO line ratios. We show that optical and CO disk sizes are on average equal, for both the CO(1-0) and CO(2-1) line. We measure the distribution of CO line luminosity ratios, finding medians (16th-84th percentile) of 0.71 (0.51-0.96) for the CO(2-1)-to-CO(1-0) ratio, 0.39 (0.24-0.53) for the CO(3-2)-to-CO(1-0) ratio, and 0.53 (0.41-0.74) for the CO(3-2)-to-CO(2-1) ratio. A companion paper presents our study of CO(2-1)'s applicability as a molecular gas mass tracer and search for trends in the CO(2-1)-to-CO(1-0) ratio. Our catalog of CO line luminosities will be publicly available with the published version of this article., Comment: Accepted version; published in ApJ

The joint moments of the derivatives of the characteristic polynomial of a random unitary matrix, and also a variant of the characteristic polynomial that is real on the unit circle, in the large matrix size limit, have been studied intensively in the past twenty five years, partly in relation to conjectural connections to the Riemann zeta-function and Hardy's function. We completely settle the most general version of the problem of convergence of these joint moments, after they are suitably rescaled, for an arbitrary number of derivatives and with arbitrary positive real exponents. Our approach relies on a hidden, higher-order exchangeable structure, that of an exchangeable array. Using these probabilistic techniques, we then give a combinatorial formula for the leading order coefficient in the asymptotics of the joint moments, when the power on the characteristic polynomial itself is a positive real number and the exponents of the derivatives are integers, in terms of a finite number of finite-dimensional integrals which are explicitly computable. Finally, we develop a method, based on a class of Hankel determinants shifted by partitions, that allows us to give an exact representation of all these joint moments, for finite matrix size, in terms of derivatives of Painlev\'e V transcendents, and then for the leading order coefficient in the large-matrix limit in terms of derivatives of solutions of the $\sigma$-Painlev\'e III' equation. Equivalently, we can represent all the joint moments of power sum linear statistics of a certain determinantal point process behind this problem in terms of derivatives of $\sigma$-Painlev\'e III' transcendents. This gives an efficient way to compute all these quantities explicitly. Our methods can be used to obtain analogous results for a number of other models sharing the same features., Comment: Typos corrected. Submitted version

Let $K$ be a local field of characteristic 0 with residue characteristic $p$. Let $G$ be an extraspecial $p$-group and let $L/K$ be a totally ramified $G$-extension. In this paper we find sufficient conditions for $L/K$ to admit a Galois scaffold. This leads to sufficient conditions for the ring of integers $\mathfrak{O}_L$ to be free of rank 1 over its associated order $\mathfrak{A}_{L/K}$, and to stricter conditions which imply that $\mathfrak{A}_{L/K}$ is a Hopf order in the group ring $K[G]$.

We present multiwavelength polarization measurements of the luminous blazar Mrk~501 over a 14-month period. The 2--8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100-ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical-infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at $>3\sigma$ above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions., Comment: Accepted for publication in ApJ

Galaxy clusters are the most massive, gravitationally-bound structures in the Universe, emerging through hierarchical structure formation of large-scale dark matter and baryon overdensities. Early galaxy ``proto-clusters'' are believed to be important physical drivers of the overall cosmic star-formation rate density and serve as ``hotspots'' for the reionization of the intergalactic medium. Our understanding of the formation of these structures at the earliest cosmic epochs is, however, limited to sparse observations of their galaxy members, or based on phenomenological models and cosmological simulations. Here we report the detection of a massive neutral, atomic hydrogen (HI) gas reservoir permeating a galaxy proto-cluster at redshift $z=5.4$, observed one billion years after the Big Bang. The presence of this cold gas is revealed by strong damped Lyman-$\alpha$ absorption features observed in several background galaxy spectra taken with JWST/NIRSpec in close on-sky projection. While overall the sightlines probe a large range in HI column densities, $N_{\rm HI} = 10^{21.7}-10^{23.5}$ cm$^{-2}$, they are similar across nearby sightlines, demonstrating that they probe the same dense, neutral gas. This observation of a massive, large-scale overdensity of cold neutral gas challenges current large-scale cosmological simulations and has strong implications for the reionization topology of the Universe., Comment: Submitted

Recent evidence supporting reionisation ending at redshift z<6 includes the rapid redshift evolution of the mean free path for Lyman-limit photons through the intergalactic medium at 510 discovered with JWST., Comment: 17 pages, 14 figures, accepted to MNRAS

We present multi-frequency (5-345 GHz) and multi-resolution radio observations of 1ES 1927+654, widely considered one of the most unusual and extreme changing-look active galactic nuclei (CL-AGN). The source was first designated a CL-AGN after an optical outburst in late 2017 and has since displayed considerable changes in X-ray emission, including the destruction and rebuilding of the X-ray corona in 2019-2020. Radio observations prior to 2023 show a faint and compact radio source typical of radio-quiet AGN. Starting in February 2023, 1ES 1927+654 began exhibiting a radio flare with a steep exponential rise, reaching a peak 60 times previous flux levels, and has maintained this higher level of radio emission for over a year to date. The 5-23 GHz spectrum is broadly similar to gigahertz-peaked radio sources, which are understood to be young radio jets less than ~1000 years old. Recent high-resolution VLBA observations at 23.5 GHz now show resolved extensions on either side of the core, with a separation of ~0.15 pc, consistent with a new and mildly relativistic bipolar outflow. A steady increase in the soft X-ray band (0.3-2 keV) concurrent with the radio may be consistent with jet-driven shocked gas, though further observations are needed to test alternate scenarios. This source joins a growing number of CL-AGN and tidal disruption events which show late-time radio activity, years after the initial outburst., Comment: Submitted to ApJ Letters 24 June 2024; Accepted 14 October 2024

The Submillimeter Array (SMA) is an array of 8 antennas operating at millimeter and submillimeter wavelengths on Maunakea, Hawaii, operated by the Smithsonian Astrophysical Observatory and Academia Sinica Institute of Astronomy and Astrophysics, Taiwan. Over the past several years, we have been preparing a major upgrade to the SMA that will replace the aging original receiver cryostats and receiver cartridges with all new cryostats and new 230 and 345 GHz receiver designs. This wideband upgrade (wSMA) will also include significantly increased instantaneous bandwidth, improved sensitivity, and greater capabilities for dual frequency observations. In this paper, we will describe the wSMA receiver upgrade and status, as well as the future upgrades that will be enabled by the deployment of the wSMA receivers., Comment: To be published in the proceedings of SPIE Astronomical Telescopes + Instrumentation 2024, paper number 13096-122

In the 1980s, Viennot developed a combinatorial approach to studying mixed moments of orthogonal polynomials using Motzkin paths. Recently, an alternative combinatorial model for these mixed moments based on lecture hall paths was introduced in arXiv:2311.12761. For sequences of orthogonal polynomials, we establish here a bijection between the Motzin paths and the lecture hall paths via a novel symmetric lecture hall graph. We use this connection to calculate the moments of the position operator in various separable quantum systems, such as the quantum harmonic oscillator and the hydrogen atom, showing that they may be expressed as generating functions of Motzkin paths and symmetric lecture hall paths, thereby providing a quantum interpretation for these paths. Our approach can be extended to other quantum systems where the wavefunctions are expressed in terms of orthogonal polynomials., Comment: 13 pages, 9 figures

The Simons Observatory (SO) is a next-generation ground-based telescope located in the Atacama Desert in Chile, designed to map the cosmic microwave background (CMB) with unprecedented precision. The observatory consists of three small aperture telescopes (SATs) and one large aperture telescope (LAT), each optimized for distinct but complementary scientific goals. To achieve these goals, optimized scan strategies have been defined for both the SATs and LAT. This paper describes a software system deployed in SO that effectively translates high-level scan strategies into realistic observing scripts executable by the telescope, taking into account realistic observational constraints. The data volume of SO also necessitates a scalable software infrastructure to support its daily data processing needs. This paper also outlines an automated workflow system for managing data packaging and daily data reduction at the site., Comment: 16 pages, 7 figures, 6 tables, to be presented at SPIE Astronomical Telescopes + Instrumentation 2024

We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars., Comment: 17 pages, 10 figures, 4 tables, Accepted for publication in A&A

Recent quasar absorption line observations suggest that reionization may end as late as $z \approx 5.3$. As a means to search for large neutral hydrogen islands at $z<6$, we revisit long dark gaps in the Ly$\beta$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$\alpha$ forest corresponding to both edges of these Ly$\beta$ dark gaps and identify a damping wing-like extended absorption profile. The average redshift of the stacked forest is $z=5.8$. By comparing these observations with reionization simulations, we infer that such a damping wing-like feature can be naturally explained if these gaps are at least partially created by neutral islands. Conversely, simulated dark gaps lacking neutral hydrogen struggle to replicate the observed damping wing features. Furthermore, this damping wing-like profile implies that the volume-averaged neutral hydrogen fraction must be $\langle x_{\rm HI} \rangle \geq 6.1 \pm 3.9\%$ at $z = 5.8$. Our results offer robust evidence that reionization extends below $z=6$., Comment: 8 pages, 5 figures, 1 table; accepted for publication in MNRAS Letters

Context. The blazar AO 0235+164, located at redshift $z=0.94$, has displayed interesting and repeating flaring activity in the past, the latest episodes occurring in 2008 and 2015. In 2020, the source brightened again, starting a new flaring episode that peaked in 2021. Aims. We study the origin and properties of the 2021 flare in relation to previous studies and the historical behavior of the source, in particular to the 2008 and 2015 flaring episodes. Methods. We analyze the multi-wavelength photo-polarimetric evolution of the source. From Very Long Baseline Array images, we derive the kinematic parameters of new components associated with the 2021 flare. We use this information to constrain a model for the spectral energy distribution of the emission during the flaring period. We propose an analytical geometric model to test whether the observed wobbling of the jet is consistent with precession. Results. We report the appearance of two new components that are ejected in a different direction than previously, confirming the wobbling of the jet. We find that the direction of ejection is consistent with that of a precessing jet.The derived period independently agrees with the values commonly found in the literature. Modeling of the spectral energy distribution further confirm that the differences between flares can be attributed to geometrical effects., Comment: Accepted 15 May 2024

The Simons Observatory (SO) is a cosmic microwave background (CMB) survey experiment that includes small-aperture telescopes (SATs) observing from an altitude of 5,200 m in the Atacama Desert in Chile. The SO SATs will cover six spectral bands between 27 and 280 GHz to search for primordial B-modes to a sensitivity of $\sigma(r)=0.002$, with quantified systematic errors well below this value. Each SAT is a self-contained cryogenic telescope with a 35$^\circ$ field of view, 42 cm diameter optical aperture, 40 K half-wave plate, 1 K refractive optics, and $<0.1$ K focal plane that holds $>12,000$ TES detectors. We describe the nominal design of the SATs and present details about the integration and testing for one operating at 93 and 145 GHz.

The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling., Comment: 46 pages, 23 figures, accepted by Astronomy & Astrophysics on August. 29, 2024

The attenuation of Ly$\alpha$ photons by neutral hydrogen in the intergalactic medium (IGM) at $z\gtrsim5$ continues to be a powerful probe for studying the epoch of reionisation. Given a framework to estimate the intrinsic (true) Ly$\alpha$ emission of high-$z$ sources, one can infer the ionisation state of the IGM during reionisation. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high-SNR QSO spectra between $5.8\lesssim\,z\lesssim\,6.6$ obtained with VLT/X-Shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Ly$\alpha$ emission profiles (obtained via covariance matrix reconstruction of the Ly$\alpha$ and N V emission lines from unattenuated high-ionisation emission line profiles; C IV, Si IV + O IV] and C III]) and (ii) distribution of ionised regions within the IGM using synthetic damping wing profiles drawn from a $1.6^3$ Gpc$^3$ reionisation simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionisation. Our median and 68th percentile constraints on the IGM neutral fraction are: $0.20\substack{+0.14\\-0.12}$ and $0.29\substack{+0.14\\-0.13}$ at $z = 6.15$~and 6.35. Further, we also report 68th percentile upper-limits of $\bar{x}_{\mathrm{H\,{\scriptscriptstyle I}}} < 0.21$, 0.20, 0.21 and 0.18 at $z = 5.8, 5.95, 6.05$~and 6.55. These results imply reionisation is still ongoing at $5.8\lesssim\,z\lesssim\,6.55$, consistent with previous results from XQR-30 (dark fraction and Ly$\alpha$ forest) along with other observational probes considered in the literature., Comment: 21 pages, 9 figures and 2 tables. Accepted for publication in MNRAS

This study explores the metal enrichment signatures attributed to the first generation of stars (PopIII) in the Universe, focusing on the E-XQR-30 sample. We aim to identify traces of Pop III metal enrichment by analyzing neutral gas in the interstellar medium of primordial galaxies and their satellite clumps, detected in absorption. To chase the chemical signature of PopIII stars, we studied metal absorption systems in the E-XQR-30 sample, selected through the detection of the OI absorption line at 1302A. The OI line is a reliable tracer of HI and allowed us to overcome the challenges posed by the Lyman-$\alpha$ forest's increasing saturation at redshifts above $\sim5$ to identify Damped Lyman-$\alpha$ systems (DLA). We detected and analyzed 29 OI systems at $z\geq5.4$, differentiating between proximate DLAs (PDLA) and intervening DLAs. Voigt function fits were applied to obtain ionic column densities, and relative chemical abundances were determined for 28 systems. These were then compared with the predictions of theoretical models. Our findings expand the study of OI systems at $z\geq5.4$ fourfold. No systematic differences were observed in the average chemical abundances between PDLAs and intervening DLAs. The chemical abundances in our sample align with literature systems at $z>4.5$, suggesting a similar enrichment pattern for this class of absorption systems. A comparison between these DLA-analogues at $4.5

We present a new suite of late-end reionization simulations performed with ATON-HE, a revised version of the GPU-based radiative transfer code ATON that includes helium. The simulations are able to reproduce the Ly$\alpha$ flux distribution of the E-XQR-30 sample of QSO absorption spectra at $5 \lesssim z \lesssim 6.2$, and show that a large variety of reionization models are consistent with these data. We explore a range of variations in source models and in the early-stage evolution of reionization. Our fiducial reionization history has a midpoint of reionization at $z = 6.5$, but we also explore an `Early' reionization history with a midpoint at $z = 7.5$ and an `Extremely Early' reionization history with a midpoint at $z = 9.5$. Haloes massive enough to host observed Ly$\alpha$ emitters are highly biased. The fraction of such haloes embedded in ionized bubbles that are large enough to allow high Ly$\alpha$ transmission becomes close to unity much before the volume filling factor of ionized regions. For our fiducial reionization history this happens at $z = 8$, probably too late to be consistent with the detection by JWST of abundant Ly$\alpha$ emission out to $z = 11$. A reionization history in our `Early' model or perhaps even our `Extremely Early' model may be required, suggesting a Thomson scattering optical depth in tension with that reported by Planck, but consistent with recent suggestions of a significantly higher value., Comment: 25 pages, 15 figures

Modern deep neural networks have achieved high performance across various tasks. Recently, researchers have noted occurrences of low-dimensional structure in the weights, Hessian's, gradients, and feature vectors of these networks, spanning different datasets and architectures when trained to convergence. In this analysis, we theoretically demonstrate these observations arising, and show how they can be unified within a generalized unconstrained feature model that can be considered analytically. Specifically, we consider a previously described structure called Neural Collapse, and its multi-layer counterpart, Deep Neural Collapse, which emerges when the network approaches global optima. This phenomenon explains the other observed low-dimensional behaviours on a layer-wise level, such as the bulk and outlier structure seen in Hessian spectra, and the alignment of gradient descent with the outlier eigenspace of the Hessian. Empirical results in both the deep linear unconstrained feature model and its non-linear equivalent support these predicted observations., Comment: 35 pages, 14 figures

One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$\alpha$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$\alpha$ absorption (DLA) wings from high column densities of neutral atomic hydrogen (HI), signifying major gas accretion events in the formation of these galaxies. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at $z=5.5-13.4$. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly$\alpha$ damping parameter $D_{\rm Ly\alpha}$ to measure and quantify the Ly$\alpha$ emission strength, HI fraction in the IGM, or local HI column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories. The fraction of strong galaxy DLAs are found to be in the range $65-95\%$ at $z>5.5$. The fraction of strong Ly$\alpha$ emitters (LAEs) is found to increase with decreasing redshift, in qualitative agreement with previous observational results, and are predominantly associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine HI gas accretion. [abridged], Comment: 20 pages, 18 figures, 3 tables. Submitted to A&A. Comments welcome! All data and catalogs are available through the DAWN JWST Archive (DJA): https://dawn-cph.github.io/dja/ and https://github.com/keheintz/jwst-primal

Intervening metal absorbers in quasar spectra at $z > 6$ can be used as probes to study the chemical enrichment of the Universe during the Epoch of Reionization (EoR). This work presents the comoving line densities ($dn/dX$) of low ionisation absorbers, namely, Mg II (2796\r{A}), C II (1334\r{A}) and O I (1302\r{A}) across $2 0.03$\r{A}. For the first time, we present the comoving line densities of 131 weak ($W < 0.3$\r{A}) intervening Mg II absorbers at $1.9 < z < 6.4$ which exhibit constant evolution with redshift similar to medium ($0.3 < W < 1.0$\r{A}) absorbers. However, the cosmic mass density of Mg II - dominated by strong Mg II systems - traces the evolution of global star formation history from redshift 1.9 to 5.5. E-XQR-30 also increases the absorption path length by a factor of 50% for C II and O I whose line densities show a rising trend towards $z > 5$, in agreement with previous works. In the context of a decline in metal enrichment of the Universe at $z > 5$, the overall evolution in the incidence rates of absorption systems can be explained by a weak - possibly soft fluctuating - UV background. Our results, thereby, provide evidence for a late reionization continuing to occur in metal-enriched and therefore, biased regions in the Universe., Comment: 20 pages, 21 figures