We present a new method for aperture masking interferometric (AMI) imaging at near-IR wavelengths using radio astronomical techniques. The method starts with derivation of interferometric visibilities from a Fourier transform of the interferograms. An iterative joint optimization process is then employed, using self-calibration of the interferometric element-based complex voltage gains (i.e. electric fields), and CLEAN deconvolution to obtain the source structure. We demonstrate the efficacy of the method using the NIRISS aperture masking interferometer on the James Webb Space Telescope (JWST) at 4.8~$μ$m and 3.8~$μ$m. Due to a number of effects (the large pixel size, charge migration, near-field optics), the method also requires an initial visibility-based amplitude normalization using observations of a well know point-source calibration star. We employ early science observations of the dusty binary Wolf-Rayet star WR137. Images with a dynamic range (peak/rms) of $\sim 240$ on the target, and $\sim 1000$ on the calibrator, are synthesized from a short integration. The self-calibration process determines the photon path-lengths through the optical system to each aperture using data on the target source itself, thereby representing an essentially 'real-time', precise wavefront error sensor. Four independent measures of the JWST mirror segment pistons (two wavelengths for two sources), agree to within 10~nm to 15~nm, comparable to the expected errors based on an analysis of closure phases on the calibrator star. Including a baseline-based phase correction improves the dynamic range of the final images by about 23\%.

We introduce a novel technique using closure amplitudes, inspired by radio interferometry, to determine with high angular resolution the two-dimensional profile of a light beam using an interferogram from a non-redundantly masked aperture. Previous techniques have required multiple interferograms or accurate estimates of the non-uniform illuminations across the aperture. In contrast, our method using closure amplitudes avoids the need to estimate the aperture illuminations while determining the two-dimensional beam shape from a single interferogram. The invariance of closure amplitudes to even time-varying aperture illuminations makes it suitable to longer averaging intervals, with potential to reducing data rates and computational overheads. By using data from the ALBA synchrotron light source to validate the method and its results against existing methods, this paper represents the first real-world application of closure amplitudes to directly determine the light beam's profile using optical interferometry in the high angular resolution regime.

Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems.

Relativistic jets are thought to play a crucial role in the formation and evolution of massive galaxies and supermassive black holes. Blazars, which are quasars with jets aligned along our line of sight, provide insights into the jetted population and have been observed up to redshifts of z = 6.1. Here, we report the discovery and multi-wavelength characterization of the blazar VLASS J041009.05-013919.88 at z = 7 (age of the Universe similar to 750 Myr), which is powered by a similar to 7 x 10(8) M-circle dot black hole. The presence of this high-redshift blazar implies a large population of similar but unaligned jetted sources in the early Universe. Our findings suggest two possible scenarios. In one, the jet in J0410-0139 is intrinsically low power but appears highly luminous due to relativistic beaming, suggesting that most ultraviolet-bright quasars at this redshift host jets. Alternatively, if J0410-0139 represents an intrinsically powerful radio source, there should be hundreds to thousands of radio-quiet quasars at z approximate to 7 with properties like those of J0410-0139, a prediction in tension with observed quasar densities based on their ultraviolet luminosity function. These results support the hypothesis that the rapid growth of black holes in the early Universe may be driven by jet-enhanced or obscured super-Eddington accretion, potentially playing a key role in forming massive black holes during the epoch of reionization.

Using the Karl G. Jansky Very Large Array (VLA), we have detected absorption lines due to carbon-monoxide, CO(J=0-1), and the cyano radical, CN(N=0-1), associated with radio galaxy B2 0902+34 at redshift z=3.4. The detection of millimeter-band absorption observed 1.5 Gyr after the Big Bang facilitates studying molecular clouds down to gas masses inaccessible to emission-line observations. The CO absorption in B2 0902+34 has a peak optical depth of τ ≥ 8.6 redshift as previously detected 21-cm absorption of neutral hydrogen (HI) gas. Each CO component traces an integrated H_2 column density of N(H2) ≥ 3x10^20 cm^-2. CN absorption is detected for both CO components, as well as for a blueshifted component not detected in CO, with CO/CN line ratios ranging from ≤0.4 to 2.4. We discuss the scenario that the absorption components originate from collections of small and dense molecular clouds that are embedded in a region with more diffuse gas and high turbulence, possibly within the influence of the central Active Galactic Nucleus or starburst region. The degree of reddening in B2 0902+34, with a rest-frame color B-K   4.2, is lower than the very red colors (B-K > 6) found among other known redshifted CO absorption systems at z<1. Nevertheless, when including also the many non-detections from the literature, a potential correlation between the absorption-line strength and B-K color is evident, giving weight to the argument that the red colors of CO absorbers are due to a high dust content.

We present millimeter observations of the host galaxy of the most distant blazar known, VLASS J041009.05-013919.88 (hereafter J0410-0139) at z = 7, using Atacama Large Millimeter/submillimeter Array (ALMA) and NOrthern Extended Millimeter Array (NOEMA) observations. The ALMA data reveal a (2.02 +/- 0.36) x 10(42) erg s(-1) [C II] 158 mu m emission line at z = 6.9964 with a [C II]-inferred star formation rate (SFR) of 58 +/- 9 M-circle dot yr(-1). We estimate a dynamical mass of M-dyn,[C II(]) = (4.6 +/- 2.0) x 10(9)M(circle dot), implying a black hole mass to host a dynamical mass ratio of 0.15(-0.05)(+0.08). The 238 GHz continuum (rest-frame IR) decreased by similar to 33% from the NOEMA to the ALMA observations taken similar to 10 months apart. The Very Large Array 3-10 GHz radio flux densities showed a similar to 37% decrease in a similar time frame, suggesting a causal connection. At face value, J0410-0139 would have the lowest [CII]-to-IR luminosity ratio of a z > 5.7 quasar reported to date (similar to 10(-4)). However, if only <20% of the measured IR luminosity was due to thermal emission from dust, the [C II]-to-IR luminosity ratio would be typical of (U)LIRGs, and the SFRs derived from [C II] and IR luminosities would be consistent. These results provide further evidence that synchrotron emission significantly contributes to the observed rest-frame IR emission of J0410-0139, similar to what has been reported in some radio-loud active galactic nuclei at z < 1.

We present millimeter observations of the host galaxy of the most distant blazar known, VLASSJ041009.05-013919.88 (hereafter J0410-0139) at z=7, using ALMA and NOEMA observations. The ALMA data reveal a 2e42 erg/s [CII] 158um emission line at z=6.9964 with a [CII]-inferred star-formation rate of 58 Msun/yr. We estimate a dynamical mass of 4.6e9 Msun, implying a black hole mass to host a dynamical mass ratio of 0.15. The 238 GHz continuum (rest-frame IR) decreased by 33 apart. The VLA 3-10 GHz radio flux densities showed a 37 similar time frame, suggesting a causal connection. At face value, J0410-0139 would have the lowest [CII]-to-IR luminosity ratio of a z>5.7 quasar reported to date ( 1e-4). However, if only <20 to thermal emission from dust, the [CII]-to-IR luminosity ratio would be typical of (U)LIRGS, and the star formation rates derived from [CII] and IR luminosities would be consistent. These results provide further evidence that synchrotron emission significantly contributes to the observed rest-frame IR emission of J0410-0139, similar to what has been reported in some radio-loud AGN at z<1.

We present JWST/NIRSpec integral field spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at z = 6.2342 caught as its host galaxy interacts with companion galaxies. We detect the spatially extended emission of several emission lines (H alpha, H beta, [O III], [N II], [S II], and He II), which we used to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detected continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with photoionization conditions typical of an Active Galactic Nucleus, the western companion shows minimal dust extinction, low metallicity (Z similar to 0.4 Z(circle dot)), and star formation driven photoionization. The eastern companion shows higher extinction and metallicity (Z similar to 0.8 Z(circle dot)) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive (M-dyn > 10(11) M-circle dot), it is still rapidly building up by accreting two relatively massive (M-star similar to 10(10) M-circle dot) companion sources. This dataset showcases the power of JWST in exposing the buildup of massive galaxies in the first gigayear of the Universe.

Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method allows us to form statistical uncertainty estimates for both our models and the recovered visibilities, which is an important ingredient in establishing robust upper limits on the Epoch of Reionisation (EoR) power spectrum. In cases where the noise is large compared to the EoR signal, this approach can constrain the systematics well enough to mitigate them down to the noise level for both systematics studied. Where the noise is smaller than the EoR, our modelling can mitigate the majority of the reflections with there being only a minor level of residual systematics, while cross-coupling sees essentially complete mitigation. Our approach performs similarly to existing filtering/fitting techniques used in the HERA pipeline, but with the added benefit of rigorously propagating uncertainties. In all cases it does not significantly attenuate the underlying signal.

We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm \mu m$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$\beta$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,H\beta}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$\alpha$-based black hole mass from the same spectrum ($M_{\rm BH, H\alpha}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $\lambda 2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $\lambda_{\rm Edd,H\beta}\sim 0.67$ ($\lambda_{\rm Edd,H\alpha}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$\beta$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $\lambda 5007$ emission with a velocity offset $\Delta v_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $\Delta v_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$\beta$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail.

Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different time-scales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time. Applications of FR filtering include separating effects that are correlated with the rotating sky versus those relative to the ground, down-weighting emission in the primary beam sidelobes, and suppressing noise. FR filtering causes the noise contributions to the visibility data to become correlated in time, however, making interpretation of subsequent averaging and error estimation steps more subtle. In this paper, we describe fringe-rate filters that are implemented using discrete prolate spheroidal sequences, and designed for two different purposes-beam sidelobe/horizon suppression (the 'mainlobe' filter), and ground-locked systematics removal (the 'notch' filter). We apply these to simulated data, and study how their properties affect visibilities and power spectra generated from the simulations. Included is an introduction to fringe-rate filtering and a demonstration of fringe-rate filters applied to simple situations to aid understanding.

Context. We present a detailed analysis of the thermal, diffuse emission of the proto-intracluster medium (proto-ICM) detected in the halo of the Spiderweb Galaxy at z = 2.16, within a radius of ∼150 kpc. Aims. Our main goal is to derive the thermodynamic profiles of the proto-ICM, establish the potential presence of a cool core and constrain the classical mass deposition rate (MDR) that may feed the nuclear and the star formation (SF) activity, and estimate the available energy budget of the ongoing feedback process. Methods. We combined deep X-ray data from Chandra and millimeter observations of the Sunyaev–Zeldovich (SZ) effect obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Results. Thanks to independent measurements of the pressure profile from the ALMA SZ observation and the electron density profile from the available X-ray data, we derived, for the first time, the temperature profile in the ICM of a z > 2 protocluster. It reveals the presence of a strong cool core (comparable to local ones) that may host a significant mass deposition flow, consistent with the measured local SF values. We also find mild evidence of an asymmetry in the X-ray surface brightness distribution, which may be tentatively associated with a cavity carved into the proto-ICM by the radio jets. In this case, the estimated average feedback power would be in excess of ∼1043 erg s−1. Alternatively, the asymmetry may be due to the young dynamical status of the halo. Conclusions. The cooling time of baryons in the core of the Spiderweb protocluster is estimated to be ∼0.1 Gyr, implying that the baryon cycle in the first stages of protocluster formation is characterized by a high-duty cycle and a very active environment. In the case of the Spiderweb protocluster, we are witnessing the presence of a strongly peaked core that is possibly hosting a cooling flow with a MDR up to 250–1000 M⊙ yr−1, responsible for feeding both the central supermassive black hole (SMBH) and the high star formation rate (SFR) observed in the Spiderweb Galaxy. This phase is expected to be rapidly followed by active galactic nucleus (AGN) feedback events, whose onset may have already left an imprint in the radio and X-ray appearance of the Spiderweb protocluster, eventually driving the ICM into a self-regulated, long-term evolution in less than one Gyr.

The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based image power spectrum and its window functions computed from the DOM images. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array Phase I configuration, to study the image power spectrum properties. The window functions show <10(-11) of the integrated power leaks from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR.

This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified.

Double-aperture young interferometry is widely used in accelerators to provide a one-dimensional beam measurement. We improve this technique by combining and further developing techniques of nonredundant, two-dimensional, aperture masking, and self-calibration from astronomy. Using visible synchrotron radiation, tests at the ALBA synchrotron show that this method provides an accurate twodimensional beam transverse characterization, even from a single 1 ms interferogram. The nonredundancy of the aperture mask in the technique enables it to be resistant to spatial phase fluctuations that might be introduced by vibration of optical components or in the laboratory atmosphere.

We develop a process of image-plane self-calibration for interferometric imaging data. The process is based on shape-orientation-size (SOS) conservation for the principal triangle in an image generated from the three fringes made from a triad of receiving elements, in situations where interferometric phase errors can be factorized into element-based terms. The basis of the SOS conservation principle is that, for a three-element array, the only possible image corruption due to an element-based phase screen is a tilt of the aperture plane, leading to a shift in the image plane. Thus, an image made from any three-element interferometer represents a true image of the source brightness, modulo an unknown translation. Image-plane self-calibration entails deriving the unknown translations for each triad image via cross-correlation of the observed triad image with a model image of the source brightness. After correcting for these independent shifts, and summing the aligned triad images, a good image of the source brightness is generated from the full array, recovering source structure at diffraction-limited resolution. The process is iterative, using improved source models based on previous iterations. We demonstrate the technique in a high signal-to-noise context, and include a configuration based on radio astronomical facilities, and simple models of double sources. We show that the process converges for the simple models considered, although convergence is slower than for aperture-plane self-calibration for large-N arrays. As currently implemented, the process is most relevant for arrays with a small number of elements. More generally, the technique provides geometric insight into closure phase and the self-calibration process. The technique is generalizable to non-astronomical interferometric imaging applications across the electromagnetic spectrum.

We explore non-redundant aperture masking to derive the size and shape of the ALBA synchrotron light source at optical wavelengths using synchrotron radiation interferometry. We show that non-redundant masks are required due to phase fluctuations arising within the experimental set-up. We also show, using closure phase, that the phase fluctuations are factorizable into element-based errors. We employ multiple masks, including 2, 3, 5, and 6 hole configurations. We develop a process for self-calibration of the element-based amplitudes (square root of flux through the aperture), which corrects for non-uniform illumination over the mask, in order to derive visibility coherences and phases, from which the source size and shape can be derived. We explore the optimal procedures to obtain the most reliable results with the 5-hole mask, based on the temporal scatter in measured coherences and closure phases. We find that the closure phases are very stable, and close to zero (within 2^o). Through uv-modeling, we consider the noise properties of the experiment and conclude that our visibility measurements per frame are likely accurate to an rms scatter of ∼ 1%.

Using the Karl G. Jansky Very Large Array, we have detected absorption lines due to carbon monoxide, CO(J = 0 -> 1), and the cyano radical, CN(N = 0 -> 1), associated with radio galaxy B2 0902+34 at redshift z = 3.4. The detection of millimeter-band absorption observed 1.5 Gyr after the Big Bang facilitates studying molecular clouds down to gas masses inaccessible to emission-line observations. The CO absorption in B2 0902+34 has a peak optical depth of tau >= 8.6% and consists of two components, one of which has the same redshift as previously detected 21 cm absorption of neutral hydrogen (H I) gas. Each CO component traces an integrated H-2 column density of N-H2 greater than or similar to 3 x 10(20) cm(-2). CN absorption is detected for both CO components, as well as for a blueshifted component not detected in CO, with CO/CN line ratios ranging from less than or similar to 0.4 to 2.4. We discuss the scenario that the absorption components originate from collections of small and dense molecular clouds that are embedded in a region with more diffuse gas and high turbulence, possibly within the influence of the central active galactic nucleus or a starburst region. The degree of reddening in B2 0902+34, with rest-frame color B - K similar to 4.2, is lower than the very red colors (B - K > 6) found among other known redshifted CO absorption systems at z < 1. Nevertheless, when including the many nondetections from the literature, a potential correlation between the absorption-line strength and B - K color is evident, giving weight to the argument that the red colors of CO absorbers are due to a high dust content.

Emittance measurements using synchrotron radiation are usually performed using x-rays to avoid diffraction limits. Interferometric techniques using visible light are also used to measure either the horizontal or the vertical beam projection. Several measurements rotating the interferometry axis are needed to obtain a full beam reconstruction. In this report we present a new interferometric multi-aperture masking technique and data analysis, inspired by astronomical methods, that are able to provide a full 2-D transverse beam reconstruction in a single acquisition. Results of beam characterization obtained at ALBA synchrotron light source will also been shown.

We present results of dust continuum and [CII]$\,158\,{\rm \mu m}$ emission line observations of a remarkably UV-luminous ($M_{\rm UV}=-21.6$) galaxy at $z=10.603$: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $\lambda_{\rm rest}=160\,{\rm \mu m}$ continuum $1\,\sigma$ depth of $13.0\,\rm{\mu Jy/beam}$ and a [CII] emission line $1\,\sigma$ sensitivity of $31\,\rm{mJy/beam\,km/s}$ using $50\,\rm{km/s}$ binning with a $\sim 2\,{\rm arcsec}$ synthesized beam. Neither dust continuum nor [CII]$\,158\,{\rm \mu m}$ line emission are detected at the expected frequency of $\nu_{\rm [CII]} = 163.791\,\rm{GHz}$ and the sky location of GN-z11. The upper limits show that GN-z11 is neither luminous in $L_{\rm IR}$ nor $L_{\rm [CII]}$, with a dust mass $3\,\sigma$ limit of ${\rm log}(M_{\rm dust}/{\rm M_{\odot}}) < 6.5-6.9$ and with a [CII] based molecular gas mass $3\,\sigma$ limit of ${\rm log}(M_{\rm mol,[CII]}/{\rm M_{\odot}}) < 9.3$. Together with radiative transfer calculations, we also investigated the possible cause of the dust poor nature of the GN-z11 showed by the blue color in the UV continuum of GN-z11 ($\beta_{\rm UV}=-2.4$), and found that $\gtrsim3\times$ deeper observations are crucial to study dust production at very high-redshift. Nevertheless, our observations show the crucial role of deep mm/submm observations of very high redshift galaxies to constrain multiple phases in the interstellar medium.