NSF’s NOIRLab

Abstract We have combined Wide-field Infrared Survey Explorer (WISE) and NEOWISE science exposures spanning from 2010 January through 2024 July to create the deepest ever full-sky maps at 3.4 μ m (W1) and 4.6 μ m (W2). These maps, which we publicly release, are based on ∼11.5 yr of WISE and NEOWISE observations in each of these two bands. As our new all-sky set of full-depth unWISE coadds uniformly incorporates observations from the entirety of the WISE and NEOWISE missions, it represents a valuable legacy dataset for Galactic and extragalactic astronomy in the infrared.

Abstract The detection of circularly polarized, low-frequency radio emission offers the tantalizing possibility of observing interactions between stars and their possible substellar companions, as well as direct emission from exoplanets. Additional follow-up of systems with radio emission is key to understanding the true origin of the emission, since multiple astrophysical mechanisms can plausibly lead to such signals. While 19 M dwarfs were detected by the LOw-Frequency ARray in circular polarization as part of the V-LoTSS survey, HD 220242 is the only F star to have a circularly polarized, low-frequency radio detection in the same survey. We conducted radial velocity (RV) follow-up with the Habitable-zone Planet Finder and combined these observations with additional archival RVs and Hipparcos-Gaia proper motion accelerations to determine that HD 220242 has a stellar companion with P = 16.79 ± 0.04 yr and a mass of 0.619 ± 0.014 M ⊙ . We use spectral energy distribution fitting and lack of any UV excess to rule out a coevolved white dwarf companion, and confirm that the companion is an M dwarf star. Given that F stars lack the coronal properties to produce such coherent emission, and the companion mass and lack of UV excess are consistent with an M dwarf, the radio emission is most plausibly associated with the companion.

The dusty and molecular torus is an elusive structure surrounding supermassive black holes, yet its importance is unequivocal for understanding feedback and accretion mechanisms. The torus and accretion disk feed the inspiraling gas onto the active nucleus, launching outflows that fundamentally connect the active nucleus's activity to the host galaxy. In this work, we utilize the aperture-masking interferometric mode onboard the JWST to achieve a resolution of 0.08" at 4.3 μm and bring out the fainter features in the central 10 pc of the Circinus galaxy. We show that most of the dust mass is located along the equatorial axis in the form of a 5×3 pc disk feeding the active nucleus. Only <1% of the dust emission arises from an arc structure composed of hot dust entrained in a molecular and ionized outflow, while the extended emission is associated with dust heated by the active galaxy at large scales.

Abstract We report observations of main-belt asteroid 2025 VZ 8 using multiple telescopes from UT 2025 November 24 to UT 2025 December 23 that confirm the presence of visible comet-like activity. These observations were motivated by online reports that an abrupt brightening event may have led to the object’s discovery on UT 2025 November 9. In all observations reported here, the object displays a tapered tail approximately 3″ long aligned with the antisolar direction. The object has an asteroidal Tisserand parameter value with respect to Jupiter ( T J = 3.371), and is thus considered an active asteroid. We measure an average apparent r ′ -band magnitude of m r ∼ 21.6 mag over the reported observing period. Due to the sudden appearance of activity far from perihelion, we suggest that the activity is more likely due to an impact or rotational destabilization, rather than sublimation.

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑mν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 wa CDM relaxes the constraints past 59 meV again, but only w 0 wa CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑mν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Abstract Among the most puzzling early discoveries of JWST are “little red dots” (LRDs), compact red sources that host broad Balmer emission lines, and in many cases exhibit a “V-shaped” change in slope in the rest-optical. The physical properties of LRDs currently have order-of-magnitude uncertainties, because models to explain the continuum of these sources differ immensely. Here, we leverage the complete selection of red sources in the RUBIES program, supplemented with public PRISM spectra, to study the origin of this V shape. By fitting a broken power law with a flexible inflection point, we find that a large fraction of red H α emitters at 2 &lt; z &lt; 6 exhibit a strong change in slope, and that all strong inflections appear associated with the Balmer limit (0.3645 μ m). Using a simple model of a reddened active galactic nucleus (AGN) with an unobscured scattered-light component, we demonstrate that the observed V shape in LRDs is unlikely to occur at any specific wavelength if the entire continuum is dominated by light from a power-law AGN continuum. In contrast, models with an intrinsic feature at the Balmer limit, such as those that are dominated by an evolved stellar population, can produce the observed spectral shapes, provided that a reddened component picks up sufficiently redward of the break. While no model can comfortably explain the full LRD spectral energy distribution, the common inflection location suggests that a single component consistently dominates the rest-frame UV optical in LRDs, and that this component is associated with T ∼ 10 4 K hydrogen.

Abstract We present an efficient estimator for higher-order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter halos than individual galaxies. When comparing the observed clustering statistics with the mock catalogs generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point auto-correlation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher-order statistic on galaxy-halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy halos in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy-dark matter connection.

ABSTRACT In preparation for the first cosmological measurements from the full shape of the Lyman-$\alpha$ (Ly $\alpha$) forest from Dark Energy Spectroscopic Instrument (DESI), we must carefully model all relevant systematics that might bias our analysis. It has been shown that random quasar redshift errors produce a smoothing effect on the mean quasar continuum in the Ly $\alpha$ forest region. This, in turn, gives rise to spurious features in the Ly $\alpha$ autocorrelation and its cross-correlation with quasars. Using synthetic data sets based on the DESI survey, we confirm that the impact on baryon acoustic oscillation measurements is small, but that a bias is introduced to parameters which depend on the full shape of our correlations. We combine a model of this contamination in the cross-correlation with a new model we introduce here for the autocorrelation. These are parametrized by three parameters, which, when included in a joint fit to both correlation functions, successfully eliminate any impact of redshift errors on our full-shape constraints. We also present a strategy for removing this contamination from real data, by removing $\sim$0.3 per cent of correlating pairs.

Abstract Much of the carbonaceous dust observed in the early Universe may originate from colliding wind binaries (CWBs) hosting hot, luminous Wolf–Rayet (W-R) stars. Downstream of the shock between the stellar winds there exists a suitable environment for dust grain formation, and the orbital motions of the stars wrap this dust into richly structured spiral geometries. The Apep system is the most extreme W-R CWB in our Milky Way: two W-R stars produce a complex spiral dust nebula, whose slow expansion has been linked to a gamma-ray burst progenitor. It has been unclear whether the O-type supergiant 0 . ″ 7 distant from the W-R + W-R binary is physically associated with the system, and whether it affects the dusty nebula. Multiepoch Very Large Telescope/VISIR and JWST/MIRI observations show that this northern companion star routinely carves a cavity in the dust nebula—the first time such an effect has been observed in a CWB—which unambiguously associates the O star as a bound component to the Apep system. These observations are used together with a new geometric model to infer the cavity geometry and the orbit of the W-R + W-R binary, yielding the first strong constraints on wind and orbital parameters. We confirm an orbital period of over 190 yr for the inner binary—nearly an order of magnitude longer than the next longest-period dust-producing W-R CWB. This, together with the confirmed classification as a hierarchical triple, cements Apep as a singular astrophysical laboratory for studying colliding winds and the terminal life stages of the most massive star systems.

Abstract The Milky Way Survey of the Dark Energy Spectroscopic Instrument (DESI) has so far observed three classical dwarf spheroidal galaxies (dSphs): Draco, Sextans, and Ursa Minor. Based on the observed line-of-sight velocities and metallicities of their member stars, we apply the axisymmetric Jeans Anisotropic Multi-Gaussian Expansion modeling (JAM) approach to recover their inner dark matter distributions. In particular, both the traditional single-population Jeans model and the multiple population chemodynamical model are adopted. With the chemodynamical model, we divide member stars of each dSph into metal-rich and metal-poor populations. The metal-rich populations are more centrally concentrated and dynamically colder, featuring lower velocity dispersion profiles than the metal-poor populations. We find a diversity of the inner density slopes γ of dark matter halos, with the best constraints by the single-population or chemodynamical models consistent with each other. The inner density slopes are <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mn>0.7</mml:mn> <mml:msubsup> <mml:mrow> <mml:mn>1</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.35</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.34</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> , <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mn>0.2</mml:mn> <mml:msubsup> <mml:mrow> <mml:mn>6</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.12</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.22</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> , and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mn>0.3</mml:mn> <mml:msubsup> <mml:mrow> <mml:mn>3</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.16</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.20</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> for Draco, Sextans, and Ursa Minor, respectively. We also present the measured astrophysical J and D factors of the three dSphs. Our results indicate that the study of the dark matter content of dSphs through stellar kinematics is still subject to uncertainties behind both the methodology and the observed data, through comparisons with previous measurements and datasets.