The next generation of wide-field cosmic microwave background (CMB) surveys are uniquely poised to open a new window for time-domain astronomy in the millimeter band. Here we explore the discovery phase space for extragalactic transients with near-term and future CMB experiments to characterize the expected population. We use existing millimeter-band light curves of known transients (gamma-ray bursts, tidal disruption events, fast blue optical transients, neutron star mergers) and theoretical models, in conjunction with known and estimated volumetric rates. Using Monte Carlo simulations of various CMB survey designs (area, cadence, depth, duration) we estimate the detection rates and the resulting light curve characteristics. We find that existing and near-term surveys will find tens to hundreds of long-duration gamma-ray bursts (LGRBs), driven primarily by detections of the reverse shock emission, and including off-axis LGRBs. Next-generation experiments (CMB-S4, CMB-HD) will find tens of fast blue optical transients (FBOTs) in the nearby universe and will detect a few tidal disruption events. CMB-HD will additionally detect a small number of short gamma-ray bursts, where these will be discovered within the detection volume of next generation gravitational wave experiments like Cosmic Explorer.

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.

Supernovae (SNe) powered by interaction with circumstellar material provide evidence for intense stellar mass loss during the final years leading up to core collapse. We have argued that during and after core neon burning, internal gravity waves excited by core convection can tap into the core fusion power and transport a super-Eddington energy flux out to the stellar envelope, potentially unbinding up to ~ 1 solar mass of material. In this work, we explore the internal conditions of SN progenitors using the MESA 1-D stellar evolution code, in search of those most susceptible to wave-driven mass loss. We focus on simple, order of magnitude considerations applicable to a wide range of progenitors. Wave-driven mass loss during core neon and oxygen fusion happens preferentially in either lower mass (<~ 20 solar mass ZAMS) stars or massive, sub-solar metallicity stars. Roughly 20 per cent of the SN progenitors we survey can excite ~ 10^46 - 10^48 erg of energy in waves that can potentially drive mass loss within a few months to a decade of core collapse. This energy can generate a circumstellar environment with 10^-3 - 1 solar mass reaching ~ 100 AU before explosion. We predict a correlation between the energy associated with pre-SN mass ejection and the time to core collapse, with the most intense mass loss preferentially happening closer to core collapse. During silicon burning, a ~ 5 day long phase for our progenitor models, wave energy may inflate ~ 10^-3 - 1 solar mass of the stellar envelope to ~ 10 - 100s of solar radii. This suggests that some nominally compact SN progenitors (Type Ibc progenitors) will have a significantly different SN shock breakout signature than traditionally assumed. We discuss the implications of our results for the core-collapse SN mechanism, Type IIn SNe, Type IIb SNe from extended progenitors (e.g., SNe 1993j and 2011dh), and observed pre-SN outbursts.

We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to streamline the extraction, preservation, integration and cross-linking of valuable data from astrophysics literature into major astronomical databases, and consequently facilitate new modes of science discovery that will better exploit the vast quantities of panchromatic and multi-dimensional data associated with the literature. We encourage authors, journal editors, referees, and publishers to implement the best practices reviewed here, as well as related recommendations from international astronomical organizations such as the International Astronomical Union (IAU) for publication of nomenclature, data, and metadata. A convenient Checklist of Recommendations for Publishing Data in the Literature is included for authors to consult before the submission of the final version of their journal articles and associated data files. We recommend that publishers of journals in astronomy and astrophysics incorporate a link to this document in their Instructions to Authors.

Supernovae (SNe) powered by interaction with circumstellar material provide evidence for intense stellar mass loss during the final years leading up to core collapse. We have argued that during and after core neon burning, internal gravity waves excited by core convection can tap into the core fusion power and transport a super-Eddington energy flux out to the stellar envelope, potentially unbinding up to ~ 1 solar mass of material. In this work, we explore the internal conditions of SN progenitors using the MESA 1-D stellar evolution code, in search of those most susceptible to wave-driven mass loss. We focus on simple, order of magnitude considerations applicable to a wide range of progenitors. Wave-driven mass loss during core neon and oxygen fusion happens preferentially in either lower mass (<~ 20 solar mass ZAMS) stars or massive, sub-solar metallicity stars. Roughly 20 per cent of the SN progenitors we survey can excite ~ 10^46 - 10^48 erg of energy in waves that can potentially drive mass loss within a few months to a decade of core collapse. This energy can generate a circumstellar environment with 10^-3 - 1 solar mass reaching ~ 100 AU before explosion. We predict a correlation between the energy associated with pre-SN mass ejection and the time to core collapse, with the most intense mass loss preferentially happening closer to core collapse. During silicon burning, a ~ 5 day long phase for our progenitor models, wave energy may inflate ~ 10^-3 - 1 solar mass of the stellar envelope to ~ 10 - 100s of solar radii. This suggests that some nominally compact SN progenitors (Type Ibc progenitors) will have a significantly different SN shock breakout signature than traditionally assumed. We discuss the implications of our results for the core-collapse SN mechanism, Type IIn SNe, Type IIb SNe from extended progenitors (e.g., SNe 1993j and 2011dh), and observed pre-SN outbursts.

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.

Understanding the optical properties of exoplanet cloud particles is a top priority. Many cloud condensates form as nonspherical particles and their optical properties can be very different from those of spheres. In this study, we focus on KCl particles, which likely form as cuboids in warm (T=500-1000K) exoplanet atmospheres. We compare the phase functions (at 532 nm wavelength) of KCl particles computed with Mie theory, the two-term Henyey-Greenstein (TTHG) approximation, laboratory data, and the discrete dipole approximation (DDA). Mie theory assumes scattering from spheres, while TTHG functions are used to approximate cloud scattering in two-stream radiative transfer models like PICASO. Laboratory measurements and DDA allow for a robust understanding of scattering from cuboid and deformed cuboid particle shapes. We input these phase functions into PICASO using cloud distributions from the cloud model Virga, to determine how different phase functions can impact the reflected-light intensities of the benchmark sub-Neptune, GJ 1214b. Simulated reflected light phase curves of GJ1214b produced using the TTHG, laboratory, and DDA phase curves differ by less than 3ppm. Our findings suggest that TTHG phase functions may be useful for approximating the scattering intensity of certain cuboid and irregular particle shapes. Future work should expand upon the wavelengths and particles considered to better determine when scattering approximations, like TTHG, may be useful in lieu of more accurate, but time-consuming laboratory measurements and/or nonspherical scattering theory.

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.

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.

Weak lensing surveys, along with most other late-Universe probes, have consistently measured a lower amplitude of the matter fluctuation spectrum, denoted by the parameter $S_8$, compared to predictions from early-Universe measurements in cosmic microwave background data. Improper modelling of nonlinear scales may partially explain these discrepancies in lensing surveys. This study investigates whether the conventional approach to addressing small-scale biases remains optimal for Stage-IV lensing surveys. We demonstrate that conventional weak lensing estimators are affected by scale leakage from theoretical biases at nonlinear scales, which influence all observed scales. Using the BNT transform, we propose an $\ell$-cut methodology that effectively controls this leakage. The Bernardeau-Nishimichi-Taruya (BNT) transform reorganises weak lensing data in $\ell$ space, aligning it with $k$ space, thereby reducing the mixing of nonlinear scales and providing a more accurate interpretation of the data. We evaluate the BNT approach by comparing HMcode, Halofit, Baryon Correction Model and AxionHMcode mass power spectrum models using Euclid-like survey configurations. Additionally, we introduce a new estimator to quantify scale leakage in both the BNT and noBNT approaches. Our findings show that BNT outperforms traditional methods, preserving cosmological constraints while significantly mitigating theoretical biases.

A primary goal of many undergraduate- and graduate-level courses in the physical sciences is to prepare students to engage in scientific research, or to prepare students for careers that leverage skillsets similar to those used by research scientists. Even for students who may not intend to pursue a career with these characteristics, exposure to the context of applications in modern research can be a valuable tool for teaching and learning. However, a persistent barrier to student participation in research is familiarity with the technical language, format, and context that academic researchers use to communicate research methods and findings with each other: the literature of the field. Astrobites, an online web resource authored by graduate students, has published brief and accessible summaries of more than 1,300 articles from the astrophysical literature since its founding in 2010. This article presents three methods for introducing students at all levels within the formal higher education setting to approaches and results from modern research. For each method, we provide a sample lesson plan that integrates content and principles from Astrobites, including step-by-step instructions for instructors, suggestions for adapting the lesson to different class levels across the undergraduate and graduate spectrum, sample student handouts, and a grading rubric.

On 2019 April 25.346 and 26.640 UT the LIGO and Virgo gravitational wave (GW) observatories announced the detection of the first candidate events in Observing Run 3 that contain at least one neutron star. S190425z is a likely binary neutron star (BNS) merger at $d_L = 156 \pm 41$ Mpc, while S190426c is possibly the first NS-BH merger ever detected, at $d_L = 377 \pm 100$ Mpc, although with marginal statistical significance. Here we report our optical follow-up observations for both events using the MMT 6.5-m telescope, as well as our spectroscopic follow-up of candidate counterparts (which turned out to be unrelated) with the 4.1-m SOAR telescope. We compare to publicly reported searches, explore the overall areal coverage and depth, and evaluate those in relation to the optical/NIR kilonova emission from the BNS merger GW170817, to theoretical kilonova models, and to short GRB afterglows. We find that for a GW170817-like kilonova, the partial volume covered spans up to about 40% for S190425z and 60% for S190426c. For an on-axis jet typical of short GRBs, the search effective volume is larger, but such a configuration is expected in at most a few percent of mergers. We further find that wide-field $\gamma$-ray and X-ray limits rule out luminous on-axis SGRBs, for a large fraction of the localization regions, although these searches are not sufficiently deep in the context of the $\gamma$-ray emission from GW170817 or off-axis SGRB afterglows. The results indicate that some optical follow-up searches are sufficiently deep for counterpart identification to about 300 Mpc, but that localizations better than 1000 deg$^2$ are likely essential.

Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to <10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model of Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data; both are feasible with a modern server. The algorithm is broadly applicable to other interferometers and is valid for wide-field and non-coplanar arrays.

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as "redundant-baseline calibration" resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions--both in data and in simulations--and present strategies for mitigating that structure.

The 21 cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation (EoR). This has led to the construction of low-frequency radio interferometric arrays, such as the Hydrogen Epoch of Reionization Array (HERA), aimed at systematically mapping this emission for the first time. Precision calibration, however, is a requirement in 21 cm radio observations. Due to the spatial compactness of HERA, the array is prone to the effects of mutual coupling, which inevitably lead to non-smooth calibration errors that contaminate the data. When unsmooth gains are used in calibration, intrinsically spectrally-smooth foreground emission begins to contaminate the data in a way that can prohibit a clean detection of the cosmological EoR signal. In this paper, we show that the effects of mutual coupling on calibration quality can be reduced by applying custom time-domain filters to the data prior to calibration. We find that more robust calibration solutions are derived when filtering in this way, which reduces the observed foreground power leakage. Specifically, we find a reduction of foreground power leakage by 2 orders of magnitude at k=0.5.

SNAP is a candidate for the Joint Dark Energy Mission (JDEM) that seeks to place constraints on the dark energy using two distinct methods. The first, Type Ia SN, is discussed in a companion white paper. The second method is weak gravitational lensing, which relies on the coherent distortions in the shapes of background galaxies by foreground mass structures. The excellent spatial resolution and photometric accuracy afforded by a 2-meter space-based observatory are crucial for achieving the high surface density of resolved galaxies, the tight control of systematic errors in the telescope's Point Spread Function (PSF), and the exquisite redshift accuracy and depth required by this project. These are achieved by the elimination of atmospheric distortion and much of the thermal and gravity loads on the telescope. The SN and WL methods for probing dark energy are highly complementary and the error contours from the two methods are largely orthogonal.

The SuperNova/Acceleration Probe (SNAP) will measure precisely the cosmological expansion history over both the acceleration and deceleration epochs and thereby constrain the nature of the dark energy that dominates our universe today. The SNAP focal plane contains equal areas of optical CCDs and NIR sensors and an integral field spectrograph. Having over 150 million pixels and a field-of-view of 0.34 square degrees, the SNAP NIR system will be the largest yet constructed. With sensitivity in the range 0.9-1.7 microns, it will detect Type Ia supernovae between z = 1 and 1.7 and will provide follow-up precision photometry for all supernovae. HgCdTe technology, with a cut-off tuned to 1.7 microns, will permit passive cooling at 140 K while maintaining noise below zodiacal levels. By dithering to remove the effects of intrapixel variations and by careful attention to other instrumental effects, we expect to control relative photometric accuracy below a few hundredths of a magnitude. Because SNAP continuously revisits the same fields we will be able to achieve outstanding statistical precision on the photometry of reference stars in these fields, allowing precise monitoring of our detectors. The capabilities of the NIR system for broadening the science reach of SNAP are discussed.

Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

We have used data from ADS, AAS, and astro-ph, to study the publishing, preprint posting, and citation patterns for papers published in the ApJ in 1999 and 2002. This allowed us to track statistical trends in author demographics, preprint posting habits, and citation rates for ApJ papers as a whole and across various subgroups and types of ApJ papers. The most interesting results are the frequencies of use of the astro-ph server across various subdisciplines of astronomy, and the impact that such posting has on the citation history of the subsequent ApJ papers. By 2002 72% of ApJ papers were posted as astro-ph preprints, but this fraction varies from 22-95% among the subfields studied. A majority of these preprints (61%) were posted after the papers were accepted at ApJ, and 88% were posted or updated after acceptance. On average, ApJ papers posted on astro-ph are cited more than twice as often as those that are not posted on astro-ph. This difference can account for a number of other, secondary citation trends, including some of the differences in citation rates between journals and different subdisciplines. Preprints clearly have supplanted the journals as the primary means for initially becoming aware of papers, at least for a large fraction of the ApJ author community. Publication in a widely-recognized peer-reviewed journal remains as the primary determinant of the impact of a paper, however. For example, conference proceedings papers posted on astro-ph are also cited twice as frequently as those that are not posted, but overall such papers are still cited 20 times less often than the average ApJ paper. These results provide insights into how astronomical research is currently disseminated by authors and ingested by readers.