Understanding the co-evolution of galaxies and active galactic nuclei (AGN) requires accurate modeling of dust-obscured systems. Recent surveys using the Mid Infrared Instrument (MIRI) onboard the James Webb Space Telescope (JWST) have uncovered a large population of dust obscured AGN, challenging current theoretical frameworks. We present an updated version of the Simulated Infrared Extragalactic Dusty Sky (SIDES) simulation framework. Our updates include modified star-forming and starburst galaxy spectral energy distribution (SED) templates as well as quiescent and AGN templates. We also incorporate a probabilistic assignment of the fraction of the IR emission that is due to an AGN. Our simulations successfully reproduce the observed MIRI source number counts, redshift distributions, and AGN population fractions. We find that AGN dominate at bright flux densities $(S_\nu \gtrsim 20\, \mu \rm Jy$ while main sequence galaxies dominate at the faint end. We also quantify the effects of cosmic variance, showing that surveys with areas below $25\, \rm arcmin^2$ suffer from $\sim 30 \%$ uncertainty in bright AGN counts. Finally, we provide diagnostic color-color diagrams and joint Near Infrared Camera (NIRCam) and MIRI flux distributions to aid interpretation of current and upcoming JWST surveys.

Modern computer systems store vast amounts of personal data, enabling advances in AI and ML but risking user privacy and trust. For privacy reasons, it is desired sometimes for an ML model to forget part of the data it was trained on. This paper presents a new approach to machine unlearning using forgetting neural networks (FNN). FNNs are neural networks with specific forgetting layers, that take inspiration from the processes involved when a human brain forgets. While FNNs had been proposed as a theoretical construct, they have not been previously used as a machine unlearning method. We describe four different types of forgetting layers and study their properties. In our experimental evaluation, we report our results on the MNIST handwritten digit recognition and fashion datasets. The effectiveness of the unlearned models was tested using Membership Inference Attacks (MIA). Successful experimental results demonstrate the great potential of our proposed method for dealing with the machine unlearning problem.

The circumgalactic medium (CGM) plays a vital role in galaxy evolution, however, studying the emission from CGM is challenging due to its low surface brightness and the complexities involved in interpreting resonant lines such as Ly$\alpha$. The near-infrared coverage, unprecedented sensitivity, and high spatial resolution of JWST enable us to study the optical strong lines associated with the extended Ly$\alpha$ "nebulae" at redshifts of 2--3. These lines serve as diagnostic tools to infer the physical conditions in the CGM gas reservoir of these systems. In deep medium-band images taken by the JWST, we serendipitously discovered the [O III] emission from the CGM around a massive interacting galaxy system at a redshift z~2.8, known to be embedded in a bright extended (100 kpc) Ly$\alpha$ "nebula." This is the first time that the [O III] lines have been detected from a Ly$\alpha$ "nebula." The JWST images reveal that the CGM gas actually resides in narrow (~ 2.5 kpc) filamentary structures with strong [O III] emission, tracing the same extent as the Ly$\alpha$ emission. An analysis of the [O III] suggests that the emitting CGM is fully ionized and is energetically dominated by mechanical heating. We also find that the density and pressure are higher than those commonly predicted by simulations of the CGM. We conclude that the observed CGM emission originates from the gas expelled by the episodic feedback processes, cooling down and enriching the CGM, while traveling a distance of at least 60 kpc. These observations demonstrate how intensive feedback processes shape gas distribution and properties in the CGM around massive halos. While access to such deep, high-resolution imaging opens up a new discovery space for investigating the CGM, it also challenges numerical simulations with respect to explaining and reproducing the exquisitely complex structures revealed by the observations.

Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,\nu}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results.

Galaxy-scale outflows are of critical importance for galaxy formation and evolution. Dust grains are the main sites for the formation of molecules needed for star formation but are also important for the acceleration of outflows that can remove the gas reservoir critical for stellar mass growth. Using the MIRI medium-resolution integral field spectrograph aboard the James Webb Space Telescope (JWST), we detect the 3.28 $\mu$m aromatic and the 3.4 $\mu$m aliphatic hydrocarbon dust features in absorption in a redshift 4.601 hot dust-obscured galaxy, blue-shifted by $\Delta$V=$-5250^{+276}_{-339}$ kms$^{-1}$ from the systemic redshift of the galaxy. The extremely high velocity of the dust indicates that the wind was accelerated by radiation pressure from the central quasar. These results pave a novel way for probing the physics of dusty outflows in active galaxies at early cosmic time.

WISE J224607.6-052634.9 (W2246-0526) is a hot dust-obscured galaxy at $z$ = 4.601, and the most luminous obscured quasar known to date. W2246-0526 harbors a heavily obscured supermassive black hole that is most likely accreting above the Eddington limit. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) in seven bands, including band 10, of the brightest far-infrared (FIR) fine-structure emission lines of this galaxy: [OI]$_{63\mu m}$, [OIII]$_{88\mu m}$, [NII]$_{122\mu m}$, [OI]$_{145\mu m}$, [CII]$_{158\mu m}$, [NII]$_{205\mu m}$, [CI]$_{370\mu m}$, and [CI]$_{609\mu m}$. A comparison of the data to a large grid of Cloudy radiative transfer models reveals that a high hydrogen density ($n_{H}\sim3\times10^3$ cm$^{-3}$) and extinction ($A_{V}\sim300$ mag), together with extreme ionization ($log(U)=-0.5$) and a high X-ray to UV ratio ($\alpha_{ox}\geq-0.8$) are required to reproduce the observed nuclear line ratios. The values of $\alpha_{ox}$ and $U$ are among the largest found in the literature and imply the existence of an X-ray-dominated region (XDR). In fact, this component explains the a priori very surprising non-detection of the [OIII]$_{88\mu m}$ emission line, which is actually suppressed, instead of boosted, in XDR environments. Interestingly, the best-fitted model implies higher X-ray emission and lower CO content than what is detected observationally, suggesting the presence of a molecular gas component that should be further obscuring the X-ray emission over larger spatial scales than the central region that is being modeled. These results highlight the need for multiline infrared observations to characterize the multiphase gas in high redshift quasars and, in particular, W2246-0526 serves as an extreme benchmark for comparisons of interstellar medium conditions with other quasar populations at cosmic noon and beyond.

A new data format for Monte Carlo (MC) events, or any structural data, including experimental data, is discussed. The format is designed to store data in a compact binary form using variable-size integer encoding as implemented in the Google's Protocol Buffers package. This approach is implemented in the ProMC library which produces smaller file sizes for MC records compared to the existing input-output libraries used in high-energy physics (HEP). Other important features of the proposed format are a separation of abstract data layouts from concrete programming implementations, self-description and random access. Data stored in ProMC files can be written, read and manipulated in a number of programming languages, such C++, JAVA, FORTRAN and PYTHON.

We report the highest-redshift detection of [O I] 63 $\mu$m from a luminous quasar, J2054-0005, at $z=6.04$ based on the Atacama Large Millimeter/sub-millimeter Array Band 9 observations. The [O I] 63 $\mu$m line luminosity is $(4.5\pm1.5) \times 10^{9}~L_{\rm \odot}$, corresponding to the [O I] 63 $\mu$m-to-far-infrared luminosity ratio of $\approx 6.7\times10^{-4}$, which is consistent with the value obtained in the local universe. Remarkably, [O I] 63 $\mu$m is as bright as [C II] 158 $\mu$m, resulting in the [O I]-to-[C II] line luminosity ratio of $1.3\pm0.5$. Based on a careful comparison of the luminosity ratios of [O I] 63 $\mu$m, [C II] 158 $\mu$m, and dust continuum emission to models of photo-dissociation regions, we find that J2054-0005 has a gas density log($n_{\rm H}$/cm$^{-3}$)$=3.7\pm0.3$ and an incident far-ultraviolet radiation field of log($G/G_{\rm 0}$)$= 3.0\pm0.1$, showing that [O I] 63 $\mu$m serves as an important coolant of the dense and warm gas in J2054-0005. A close examination of the [O I] and [C II] line profiles suggests that the [O I] line may be partially self-absorbed, however deeper observations are needed to verify this conclusion. Regardless, the gas density and incident radiation field are in a broad agreement with the values obtained in nearby star-forming galaxies and objects with [O I] 63 $\mu$m observations at $z=1-3$ with the Herschel Space Observatory. These results demonstrate the power of ALMA high-frequency observations targeting [O I] 63 $\mu$m to examine the properties of photo-dissociation regions in high-redshift galaxies.

The paper describes a curricular framework for introducing microcontroller programming in the University Physics lab. The approach makes use of Modeling Instruction, an effective approach for teaching science at the secondary level in which student learn the standard material by developing and deploying models of the physical world. In our approach, students engage with a context-rich problem that can be solved with one or more sensors and a microcontroller. The solution path we describe then consists of developing a mathematical model for how the sensors' data can be mapped to a meaningful measurement, and further, developing an algorithmic model that will be implemented in the microcontroller. Once the system is developed and implemented, students are given an array of similar problems in which they can deploy their data collection system. Results from the implementation of this idea, in two University Physics sections, using Arduino microcontrollers, are also described.

We describe a procedure in which elementary education students fabricate 3-cell, $\approx5 ~\mathrm{Watt}$, $1.6~\mathrm{V}$ solar panels. The panels are manufactured at a cost of about \6$ each. The paper also describes a series of characterization activities that we have had our students complete with the panels once fabrication is complete. Specifically, our students have determined: the dependence of panel power on load resistance, the variation of open-circuit panel voltage between different students' work, and energy capture efficiency. The activity also allows opportunities for discussion of manufacturing and process improvement which can provide unusual connections to engineering and technology learning outcomes.

Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration.

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 &lt;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.

We report integral field spectroscopy observations with the Near-Infrared Spectrograph on board JWST targeting the 60 kpc environment surrounding the most luminous obscured quasar known at $z=4.6$. We detect ionized gas filaments on 40 kpc scales connecting a network of merging galaxies likely to form a cluster. We find regions of low ionization consistent with large-scale shock excitation surrounding the central dust-obscured quasar, out to distances nearly eight times the effective stellar radius of the quasar host galaxy. In the nuclear region, we find an ionized outflow driven by the quasar with velocities reaching 13,000 km s$^{-1}$, one of the fastest discovered to date with an outflow rate of 2000 M$_\odot$ yr$^{-1}$ and a kinetic luminosity of 6$\times10^{46}$ erg s$^{-1}$ resulting in coupling efficiency between the bolometric luminosity of the quasar and the outflow of 5%. The kinetic luminosity of the outflow is sufficient to power the turbulent motion of the gas on galactic and circumgalactic scales and is likely the primary driver of the radiative shocks on interstellar medium and circumgalactic medium scales. This provides compelling evidence supporting long-standing theoretical predictions that powerful quasar outflows are a main driver in regulating the heating and accretion rate of gas onto massive central cluster galaxies.

Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required.

We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $\Delta^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $\Delta^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars.

Measurements of the one-point probability distribution function and higher-order moments (variance, skewness, and kurtosis) of the high-redshift 21 cm fluctuations are among the most direct statistical probes of the non-Gaussian nature of structure formation and evolution during reionization. However, contamination from astrophysical foregrounds and instrument systematics pose significant challenges in measuring these statistics in real observations. In this work, we use forward modelling to investigate the feasibility of measuring 21 cm one-point statistics through a foreground avoidance strategy. Leveraging the characteristic wedge-shape of the foregrounds in k-space, we apply a wedge-cut filter that removes the foreground contaminated modes from a mock data set based on the Hydrogen Epoch of Reionization Array (HERA) instrument, and measure the one-point statistics from the image-space representation of the remaining non-contaminated modes. We experiment with varying degrees of wedge-cutting over different frequency bandwidths and find that the centre of the band is the least susceptible to bias from wedge-cutting. Based on this finding, we introduce a rolling filter method that allows reconstruction of an optimal wedge-cut 21~cm intensity map over the full bandwidth using outputs from wedge-cutting over multiple sub-bands. We perform Monte Carlo simulations to show that HERA should be able to measure the rise in skewness and kurtosis near the end of reionization with the rolling wedge-cut method if foreground leakage from the Fourier transform window function can be controlled.

Accurately determining gas-phase metal-abundances within galaxies is critical as metals strongly affect the physics of the interstellar medium (ISM). To date, the vast majority of widely-used gas-phase abundance-indicators rely on emission from bright optical-lines, whose emissivities are highly sensitive to the electron temperature. Alternatively, direct-abundance methods exist that measure the temperature of the emitting gas directly, though these methods usually require challenging observations of highly-excited auroral lines. Low-lying far-infrared (FIR) fine-structure lines are largely insensitive to electron temperature and thus provide an attractive alternative to optically-derived abundances. Here, we introduce the far-infrared abundances (FIRA) project, which employs these FIR transitions, together with both radio free-free emission and hydrogen recombination-lines, to derive direct, absolute gas-phase oxygen-abundances. Our first target is M101, a nearby spiral-galaxy with a relatively steep abundance gradient. Our results are consistent with the O$^{++}$ electron-temperatures and absolute oxygen-abundances derived using optical direct-abundance methods by the CHemical Abundance Of Spirals (CHAOS) program, with a small difference ($\sim$ 1.5$\sigma$) in the radial abundance-gradients derived by the FIR/free-free-normalized vs. CHAOS/direct-abundance techniques. This initial result demonstrates the validity of the FIRA methodology $-$ with the promise of determining absolute metal-abundances within dusty star-forming galaxies, both locally and at high redshift.

We present [NII] 205 $\mu$m fine structure line observations of three submillimeter galaxies (SMGs) and three quasar host galaxies at 4$\lesssim$z$\lesssim$6 using the Institut de radioastronomie millim\'etrique (IRAM) interferometer. The [NII] emission is detected in three sources, and we report detections of the underlying dust continuum emission in all sources. The observed [NII]-to-infrared luminosity ratio spans at least 0.5 dex for our sources. Comparing our estimates with sources detected in the [NII] 205 $\mu$m at similar redshifts shows that the overall [NII]-to-IR luminosity ratio spans over a dex in magnitude from L$_{[NII]}$/L$_{IR}$ ~ 10$^{-4}$ - 10$^{-5}$ and follows the trend of the so-called [NII] fine structure line deficit observed in (ultra)-luminous infrared galaxies in the local Universe. The [CII]-to-[NII] luminosity ratio is >10 for most of our sources, indicating that the bulk of the [CII] 158 $\mu$m line emission (f([CII]$^{PDR}$)>75%) arises from the neutral medium. From our analysis, we do not find significant differences in the [NII] 205 $\mu$m emission and the respective ratios between SMGs and QSOs, suggesting a negligible contribution to the boosting of [NII] 205 $\mu$m emission due to the active galactic nucleus (AGN) photoionization. Future investigations involving other fine structure lines and optical diagnostics will provide further insight into a suite of ionized medium properties and reveal the diversity between AGN and non-AGN environments.

We present spatially resolved dust-continuum ALMA observations from rest-frame $\sim$60 to $\sim$600 $\mu$m (bands 3-10) of the hyperluminous hot dust-obscured galaxy (hot DOG) WISE J224607.6-052634.9 (W2246-0526), at redshift $z=4.6$. W2246-0526 is interacting with at least three companion galaxies, forming a system connected by tidal streams. We model the multiwavelength ALMA observations of the dust continuum using a modified blackbody, from which we derive the dust properties (mass, emissivity index, area of the emitting region, and temperature) in the hot DOG and resolved structures across a region of nearly $\sim$50 kpc. The peak temperature at the location of the hot DOG, $\sim$110 K, is likely the consequence of heating by the central quasar. The dust temperature drops to $\sim$40 K at a radius of $\sim$8 kpc, suggesting that heating by the quasar beyond that distance is nondominant. The dust in the connecting streams between the host and companion galaxies is at temperatures between 30-40 K, typical of starburst galaxies, suggesting it is most likely heated by recent, in-situ star formation. This is the first time dust properties are spatially resolved over several tens of kpc in a galaxy system beyond Cosmic Noon --this is more than six times the scales previously probed in galaxies at those redshifts.

We present spatially resolved Alma Band-9 observations of the [CII] 158 $\mu$m fine structure line from an optically selected quasar, SDSS J100038.01+020822.4 (J1000), at z=1.8275. By utilizing [OI] 63 $\mu$m line observations from Herschel/PACS and constructing a detailed dust SED using Herschel and Spitzer archival imaging data, we show that the [CII] line emission is well explained by a photodissociation region (PDR) model, in which the emission arises from the surfaces of molecular clouds exposed to far-UV radiation fields $\sim 5\cdot10^3$ times the local interstellar radiation field (G$_0$). We find a factor of 30 variation in spatially resolved [CII]/Far-IR continuum across the source which is explained by the reduced fraction of cooling via [CII] line emission at such high far-UV field strengths. By matching derived PDR parameters to the observed far-IR line and continuum intensities we derive cloud size-scales and find that typical cloud radii in J1000 are $\sim$3.5 pc perhaps indicating an ISM that is highly fractured due to intense star formation activity. We model the galaxy dynamically and find that the [CII] emission is contained within a compact, dynamically cold disk with v/$\sigma$=6.2, consistent with cosmological simulations. We also report the discovery of a companion galaxy to j1000 confirmed by the detection of [CII] and use recently obtained JWST/NirCAM imaging of the system to argue for J1000 being an interacting system. With total stellar mass $\sim 1.5 \times 10^{10}$ M$_\odot$ and main-component dynamical mass $\gtrsim 10^{11}$ M$_\odot$, the J1000 system is a progenitor to the most massive galaxies seen in the local Universe.