Search

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: 68 pages, more than 10 figures

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: Submitted to Astronomy and Astrophysics. 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 observations of the \textsc{Polarbear} project. Although the detectors are sensitive to linear polarization, we can measure the circular polarization because a continuously rotating half-wave plate in the optics converts part of circular polarization into linear polarization. The atmospheric circular polarization signal appears as a modulated signal at twice the frequency of rotation of the half-wave plate. We reconstruct the azimuthal gradient of the circular polarization signal and measure the dependencies on the scanning azimuth and the detector bandpass. We compare the signal with a simulation based on atmospheric emission theory, the detector bandpass, and the half-wave plate 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)}$, which demonstrates that our measurement is consistent with theoretical prediction. This result validates our understanding of the instrument and reinforces the feasibility of measuring the circular polarization using the imperfection of the half-wave plate. Quantifying atmospheric circular polarization is the first step toward conducting a search for cosmological circular polarization at these wavelengths., Comment: 22 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: Submitted to 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: Resubmitted to ApJ after incorporating reviewer feedback

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

The vSPACE experimental proof-of-concept (PoC) on the TrueElect[Anon][Creds] protocol presents a novel approach to secure, private, and scalable elections, extending the TrueElect and ElectAnon protocols with the integration of AnonCreds SSI (Self-Sovereign Identity). Such a protocol PoC is situated within a Zero-Trust Architecture (ZTA) and leverages confidential computing, continuous authentication, multi-party computation (MPC), and well-architected framework (WAF) principles to address the challenges of cybersecurity, privacy, and trust over IP (ToIP) protection. Employing a Kubernetes confidential cluster within an Enterprise-Scale Landing Zone (ESLZ), vSPACE integrates Distributed Ledger Technology (DLT) for immutable and certifiable audit trails. The Infrastructure as Code (IaC) model ensures rapid deployment, consistent management, and adherence to security standards, making vSPACE a future-proof solution for digital voting systems.

We study the effects due to mismatches in passbands, polarization angles, and temperature and polarization calibrations in the context of the upcoming cosmic microwave background experiment Simons Observatory (SO). Using the SO multi-frequency likelihood, we estimate the bias and the degradation of constraining power in cosmological and astrophysical foreground parameters assuming different levels of knowledge of the instrumental effects. We find that incorrect but reasonable assumptions about the values of all the systematics examined here can have significant effects on cosmological analyses, hence requiring marginalization approaches at the likelihood level. When doing so, we find that the most relevant effect is due to bandpass shifts. When marginalizing over them, the posteriors of parameters describing astrophysical microwave foregrounds (such as radio point sources or dust) get degraded, while cosmological parameters constraints are not significantly affected. Marginalization over polarization angles with up to 0.25$^\circ$ uncertainty causes an irrelevant bias $\lesssim 0.05 \sigma$ in all parameters. Marginalization over calibration factors in polarization broadens the constraints on the effective number of relativistic degrees of freedom $N_\mathrm{eff}$ by a factor 1.2, interpreted here as a proxy parameter for non standard model physics targeted by high-resolution CMB measurements., Comment: 61 pages (27 of which are appendices), 26 figures, 12 tables. (Replacement to update the acknowledgments)

The Crab Nebula, also known as Tau A, is a polarized astronomical source at millimeter wavelengths. It has been used as a stable light source for polarization angle calibration in millimeter-wave astronomy. However, it is known that its intensity and polarization vary as a function of time at a variety of wavelengths. Thus, it is of interest to verify the stability of the millimeter-wave polarization. If detected, polarization variability may be used to better understand the dynamics of Tau~A, and for understanding the validity of Tau~A as a calibrator. One intriguing application of such observation is to use it for the search of axion-like particles (ALPs). Ultralight ALPs couple to photons through a Chern-Simons term, and induce a temporal oscillation in the polarization angle of linearly polarized sources. After assessing a number of systematic errors and testing for internal consistency, we evaluate the variability of the polarization angle of the Crab Nebula using 2015 and 2016 observations with the 150 GHz POLARBEAR instrument. We place a median 95% upper bound of polarization oscillation amplitude $A < 0.065^\circ$ over the oscillation frequencies from $0.75~\mathrm{year}^{-1}$ to $0.66~\mathrm{hour}^{-1}$. Assuming that no sources other than ALP are causing Tau A's polarization angle variation, that the ALP constitutes all the dark matter, and that the ALP field is a stochastic Gaussian field, this bound translates into a median 95% upper bound of ALP-photon coupling $g_{a\gamma\gamma} < 2.16\times10^{-12}\,\mathrm{GeV}^{-1}\times(m_a/10^{-21} \mathrm{eV})$ in the mass range from $9.9\times10^{-23} \mathrm{eV}$ to $7.7\times10^{-19} \mathrm{eV}$. This demonstrates that this type of analysis using bright polarized sources is as competitive as those using the polarization of cosmic microwave background in constraining ALPs., Comment: 24 pages, 19 figures, 5 tables

The breaking of surface gravity waves is a key process contributing to air-sea fluxes and turbulent ocean mixing. The highly nonlinear nature of wave breaking, combined with the challenges of observing this process in a laboratory or field setting, leaves our understanding of the energetic processes underpinning wave breaking incomplete. Progress towards refining this understanding was made in a recent study (D. G. Boettger et. al., An energetic signature for breaking inception in surface gravity waves, Journal of Fluid Mechanics 959, A33 (2023)), which identified an energetic signature in the wave kinetic energy evolution that preceded breaking onset and correlated with the strength of the breaking event. In this study, we examine the influence of wind forcing on this energetic signature. We develop a numerical wave tank that simulates wind flowing over mechanically generated waves and construct an ensemble of cases with varying wave steepness and wind forcing speed. The wind is shown to modulate the wave geometry and elevate kinetic energy at crest tip by up to 35 %. Despite these influences, the energetic inception signature was found to robustly indicate breaking inception in all cases examined, with a threshold value in the kinetic energy growth rate at this instant separating breaking and non-breaking waves. Under wind forcing, the timing of the energetic inception point occurred slightly earlier that unforced breaking waves, giving advance warning of breaking 0.3 wave periods prior to breaking onset., Comment: Accepted for publication

The Simons Observatory (SO) experiment is a cosmic microwave background (CMB) experiment located in the Atacama Desert, Chile. The SO' s small aperture telescopes (SATs) consist of three telescopes designed for precise CMB polarimetry at large angular scales. Each SAT uses a cryogenic rotating half-wave plate (HWP) as a polarization modulator to mitigate atmospheric 1/f noise and other systematics. To realize efficient polarization modulation over the observation bands, we fabricated an achromatic HWP (AHWP) consisting of three sapphire plates with anti-reflection coatings. The AHWP is designed to have broadband modulation efficiency and transmittance. This paper reports on the design and the preliminary characterization of the AHWPs for SATs.

We present resolved images of IRAS 23077+6707 ("Dracula's Chivito") in 1.3 mm/225 GHz thermal dust and CO gas emission with the Submillimeter Array (SMA) and optical (0.5-$0.8\,\mu\mathrm{m}$) scattered light with the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). The Pan-STARRS data show a bipolar distribution of optically scattering dust that is characteristic for disks observed at high inclinations. Its scattered light emission spans ${\sim}14''$, with two highly asymmetric filaments extending along the upper bounds of each nebula by ${\sim}9''$. The SMA data measure 1.3 mm continuum dust as well as $^{12}$CO, $^{13}$CO and C$^{18}$O $J$=2$-$1 line emission over $12''$-$14''$ extents, with the gas presenting the typical morphology of a disk in Keplerian rotation, in both position-velocity space and in each CO line spectrum. IRAS 23077+6707 has no reported distance estimate, but if it is located in the Cepheus star-forming region (180-800 pc), it would have a radius spanning thousands of au. Taken together, we infer IRAS 23077+6707 to be a giant and gas-rich edge-on protoplanetary disk, which to our knowledge is the largest in extent so far discovered., Comment: published in The Astrophysical Journal Letters 967:L2 (9pp), 2024 May 20

We present evidence of variable 1.3 millimeter emission from the 1-3 Myr, SpT G2-G5 class III YSO, HD~283572. HD~283572 was observed on 8 dates with the Submillimeter Array between 2021 December and 2023 May, a total on-source time of 10.2 hours, probing a range of timescales down to 5.2 seconds. Averaging all data obtained on 2022 Jan 17 shows a 4.4 mJy ($8.8\sigma$) point source detection with a negative spectral index ($\alpha{=}{-2.7}{\pm}1.2$), with peak emission rising to 13.8 mJy in one 3 minute span, and 25 mJy in one 29.7 second integration ($L_\nu=4.7\times10^{17}$ erg s$^{-1}$ Hz$^{-1}$). Combining our data for the other 7 dates shows no detection, with an rms noise of 0.24 mJy beam$^{-1}$. The stochastic millimeter enhancements on time frames of seconds--minutes--hours with negative spectral indices are most plausibly explained by synchrotron or gyro-synchrotron radiation from stellar activity. HD 283572's 1.3 mm light-curve has similarities with variable binaries, suggesting HD 283572's activity may have been triggered by interactions with an as-yet undetected companion. We additionally identify variability of HD 283572 at 10 cm, from VLASS data. This study highlights the challenges of interpreting faint mm emission from evolved YSOs that may host tenuous disks, and suggests that a more detailed temporal analysis of spatially unresolved data is generally warranted. The variability of class III stars may open up new ground for understanding the physics of flares in the context of terrestrial planet formation., Comment: Accepted for publication in ApJL. 12 pages, inc. 5 figures, 2 tables

3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $\nu_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored., Comment: 15 pages, 6 figures, published in A&A

The spectra of high-redshift ($z\gtrsim 6$) quasars contain valuable information on the progression of the Epoch of Reionization (EoR). At redshifts $z<6$, the observed Lyman-series forest shows that the intergalactic medium (IGM) is nearly ionized, while at $z>7$ the observed quasar damping wings indicate high neutral gas fractions. However, there remains a gap in neutral gas fraction constraints at $6\lesssim z \lesssim 7$ where the Lyman series forest becomes saturated but damping wings have yet to fully emerge. In this work, we use a sample of 18 quasar spectra at redshifts $6.0