Engineering topics

Bacon, D. (ORCID:0000000225628537)

Publications and source records attributed to Bacon, D. (ORCID:0000000225628537).

17 records

DOE OSTI2026

Robust Measurement of Stellar Streams around the Milky Way: Correcting Spatially Variable Observational Selection Effects in Optical Imaging Surveys

Observations of density variations in stellar streams are a promising probe of low-mass dark matter substructure in the Milky Way. However, survey systematics such as variations in seeing and sky brightness can also induce artificial fluctuations in the observed densities of known stellar streams. These variations arise because survey conditions affect both object detection and star–galaxy misclassification rates. To mitigate these effects, we use Balrog synthetic source injections in the Dark Energy Survey (DES) Y3 data to calculate detection rate variations and classification rates as functions of survey properties. We show that these rates are nearly separable with respect to survey properties and can be estimated with sufficient statistics from the synthetic catalogs. Applying these corrections reduces the standard deviation of relative detection rates across the DES footprint by a factor of 5, and our corrections significantly change the inferred linear density of the Phoenix stream when including faint objects. Additionally, for artificial streams with DES-like survey properties we are able to recover density power spectra with reduced bias. We also find that uncorrected power-spectrum results for Legacy Survey of Space and Time (LSST)-like data can be around 5 times more biased, highlighting the need for such corrections in future ground-based surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

DELVE Milky Way Satellite Galaxy Census. I. Satellite Population and Survey Selection Function in DES, DELVE, and Pan-STARRS

The properties of Milky Way satellite galaxies have important implications for galaxy formation, reionization, and the fundamental physics of dark matter. However, the population of Milky Way satellites includes the faintest known galaxies, and current observations are incomplete. To understand the impact of observational selection effects on the known satellite population, we perform rigorous, quantitative estimates of the Milky Way satellite galaxy detection efficiency in three wide-field survey datasets: the Dark Energy Survey Year 6, the DECam Local Volume Exploration Data Release 3, and the Pan-STARRS1 Data Release 1. Together, these surveys cover ∼13,600 deg 2 to g ∼ 24.0 and ∼27,700 deg 2 to g ∼ 22.5, spanning ∼91% of the high-Galactic-latitude sky (∣b∣ ≥ 15°). We apply multiple detection algorithms over the combined footprint and recover 49 known satellites above a strict census detection threshold. To characterize the sensitivity of our census, we run our detection algorithms on a large set of simulated galaxies injected into the survey data, which allows us to develop models that predict the detectability of satellites as a function of their properties. We then fit an empirical model to our data and infer the luminosity function, radial distribution, and size–luminosity relation of Milky Way satellite galaxies. Our empirical model predicts a total of $265^{+79}_{-47}$ satellite galaxies with −20 ≤ M V ≤ 0, half-light radii of 15 ≤ r 1/2 , (pc) ≤ 3000, and galactocentric distances of 10 ≤ D GC (kpc) ≤ 300. We also identify a mild anisotropy in the angular distribution of the observed galaxies, at a significance of ∼2σ, which can be attributed to the clustering of satellites associated with the LMC.

Tan, Chin Yi [Univ. of Chicago, IL (United States)↗

DOE OSTI2026

Ultra-faint Milky Way Satellites Discovered in Carina, Phoenix, and Telescopium with DELVE Data Release 3

We report the discovery of three Milky Way satellite candidates: Carina IV, Phoenix III, and DELVE 7, in the third data release of the DECam Local Volume Exploration survey (DELVE). The candidate systems were identified by cross-matching results from two independent search algorithms. All three are extremely faint systems composed of old, metal-poor stellar populations (τ ≳ 10 Gyr, [Fe/H] ≲−1.4). Carina IV (M V = −2.8; r 1/2 = 40 pc) and Phoenix III (M V = −1.2; r 1/2 = 19 pc) have half-light radii that are consistent with the known population of dwarf galaxies, while DELVE 7 (M V = 1.2; r 1/2 = 2 pc) is very compact and seems more likely to be a star cluster, though its nature remains ambiguous without spectroscopic follow-up. The Gaia proper motions of stars in Carina IV ($M_{\star} = 2250^{+1180}_{-830} M_⊙$) indicate that it is unlikely to be associated with the LMC, while DECam CaHK photometry confirms that its member stars are metal poor. Phoenix III ($M_{\star} = 520^{+660}_{-290} M_⊙$) is the faintest known satellite in the extreme outer stellar halo (D GC > 100 kpc), while DELVE 7 ($M_{\star} = 60^{+120}_{-40} M_⊙$) is the faintest known satellite with D GC > 20 kpc.

Tan, Chin Yi [Univ. of Chicago, IL (United States)↗

DOE OSTI2025

Euclid I. Overview of the Euclid mission

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015–2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14 000 deg 2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

cosmology↗

DOE OSTI2024

The Dark Energy Survey Supernova Program: Cosmological Biases from Host Galaxy Mismatch of Type Ia Supernovae

Abstract Redshift measurements, primarily obtained from host galaxies, are essential for inferring cosmological parameters from type Ia supernovae (SNe Ia). Matching SNe to host galaxies using images is nontrivial, resulting in a subset of SNe with mismatched hosts and thus incorrect redshifts. We evaluate the host galaxy mismatch rate and resulting biases on cosmological parameters from simulations modeled after the Dark Energy Survey 5 Yr (DES-SN5YR) photometric sample. For both DES-SN5YR data and simulations, we employ the directional light radius method for host galaxy matching. In our SN Ia simulations, we find that 1.7% of SNe are matched to the wrong host galaxy, with redshift differences between the true and matched hosts of up to 0.6. Using our analysis pipeline, we determine the shift in the dark energy equation of state parameter (Δ w ) due to including SNe with incorrect host galaxy matches. For SN Ia–only simulations, we find Δ w = 0.0013 ± 0.0026 with constraints from the cosmic microwave background. Including core-collapse SNe and peculiar SNe Ia in the simulation, we find that Δ w ranges from 0.0009 to 0.0032, depending on the photometric classifier used. This bias is an order of magnitude smaller than the expected total uncertainty on w from the DES-SN5YR sample of ∼0.03. We conclude that the bias on w from host galaxy mismatch is much smaller than the uncertainties expected from the DES-SN5YR sample, but we encourage further studies to reduce this bias through better host-matching algorithms or selection cuts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Designing an Optimal Kilonova Search Using DECam for Gravitational-wave Events

We address the problem of optimally identifying all kilonovae detected via gravitational-wave emission in the upcoming LIGO/Virgo/KAGRA observing run, O4, which is expected to be sensitive to a factor of ~7 more binary neutron star (BNS) alerts than previously. Electromagnetic follow-up of all but the brightest of these new events will require >1 m telescopes, for which limited time is available. We present an optimized observing strategy for the DECam during O4. We base our study on simulations of gravitational-wave events expected for O4 and wide-prior kilonova simulations. We derive the detectabilities of events for realistic observing conditions. We optimize our strategy for confirming a kilonova while minimizing telescope time. For a wide range of kilonova parameters, corresponding to a fainter kilonova compared to GW170817/AT 2017gfo, we find that, with this optimal strategy, the discovery probability for electromagnetic counterparts with the DECam is ~80% at the nominal BNS gravitational-wave detection limit for O4 (190 Mpc), which corresponds to an ~30% improvement compared to the strategy adopted during the previous observing run. For more distant events (~330 Mpc), we reach an ~60% probability of detection, a factor of ~2 increase. For a brighter kilonova model dominated by the blue component that reproduces the observations of GW170817/AT 2017gfo, we find that we can reach ~90% probability of detection out to 330 Mpc, representing an increase of ~20%, while also reducing the total telescope time required to follow up events by ~20%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Synchronous Rotation in the (136199) Eris–Dysnomia System

We combine photometry of Eris from a 6 month campaign on the Palomar 60 inch telescope in 2015, a 1 month Hubble Space Telescope WFC3 campaign in 2018, and Dark Energy Survey data spanning 2013–2018 to determine a light curve of definitive period 15.771 ± 0.008 days (1σ formal uncertainties), with nearly sinusoidal shape and peak-to-peak flux variation of 3%. This is consistent at part-per-thousand precision with the P = 15.785 90 ± 0.00005 day sidereal period of Dysnomia's orbit around Eris, strengthening the recent detection of synchronous rotation of Eris by Szakáts et al. with independent data. Photometry from Gaia are consistent with the same light curve. We detect a slope of 0.05 ± 0.01 mag per degree of Eris's brightness with respect to illumination phase averaged across g, V, and r bands, intermediate between Pluto's and Charon's values. Variations of 0.3 mag are detected in Dysnomia's brightness, plausibly consistent with a double-peaked light curve at the synchronous period. The synchronous rotation of Eris is consistent with simple tidal models initiated with a giant-impact origin of the binary, but is difficult to reconcile with gravitational capture of Dysnomia by Eris. The high albedo contrast between Eris and Dysnomia remains unexplained in the giant-impact scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Timing the r -process Enrichment of the Ultra-faint Dwarf Galaxy Reticulum II

The ultra-faint dwarf galaxy Reticulum II (Ret II) exhibits a unique chemical evolution history, with ${72}_{-12}^{+10}$% of its stars strongly enhanced in r-process elements. We present deep Hubble Space Telescope photometry of Ret II and analyze its star formation history. As in other ultra-faint dwarfs, the color–magnitude diagram is best fit by a model consisting of two bursts of star formation. If we assume that the bursts were instantaneous, then the older burst occurred around the epoch of reionization, forming ~80% of the stars in the galaxy, while the remainder of the stars formed ~3 Gyr later. When the bursts are allowed to have nonzero durations, we obtain slightly better fits. The best-fitting model in this case consists of two bursts beginning before reionization, with approximately half the stars formed in a short (100 Myr) burst and the other half in a more extended period lasting 2.6 Gyr. Considering the full set of viable star formation history models, we find that 28% of the stars formed within 500 ± 200 Myr of the onset of star formation. The combination of the star formation history and the prevalence of r-process-enhanced stars demonstrates that the r-process elements in Ret II must have been synthesized early in its initial star-forming phase. We therefore constrain the delay time between the formation of the first stars in Ret II and the r-process nucleosynthesis to be less than 500 Myr. This measurement rules out an r-process source with a delay time of several Gyr or more, such as GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

DeepZipper. II. Searching for Lensed Supernovae in Dark Energy Survey Data with Deep Learning

Gravitationally lensed supernovae (LSNe) are important probes of cosmic expansion, but they remain rare and difficult to find. Current cosmic surveys likely contain 5–10 LSNe in total while next-generation experiments are expected to contain several hundred to a few thousand of these systems. We search for these systems in observed Dark Energy Survey (DES) five year SN fields—10 3 sq. deg. regions of sky imaged in the griz bands approximately every six nights over five years. To perform the search, we utilize the DeepZipper approach: a multi-branch deep learning architecture trained on image-level simulations of LSNe that simultaneously learns spatial and temporal relationships from time series of images. We find that our method obtains an LSN recall of 61.13% and a false-positive rate of 0.02% on the DES SN field data. DeepZipper selected 2245 candidates from a magnitude-limited (m i < 22.5) catalog of 3,459,186 systems. We employ human visual inspection to review systems selected by the network and find three candidate LSNe in the DES SN fields.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The Dark Energy Survey Supernova Program results: type Ia supernova brightness correlates with host galaxy dust

ABSTRACT Cosmological analyses with type Ia supernovae (SNe Ia) often assume a single empirical relation between colour and luminosity (β) and do not account for varying host-galaxy dust properties. However, from studies of dust in large samples of galaxies, it is known that dust attenuation can vary significantly. Here, we take advantage of state-of-the-art modelling of galaxy properties to characterize dust parameters (dust attenuation AV, and a parameter describing the dust law slope RV) for 1100 Dark Energy Survey (DES) SN host galaxies. Utilizing optical and infrared data of the hosts alone, we find three key aspects of host dust that impact SN cosmology: (1) there exists a large range (∼1–6) of host RV; (2) high-stellar mass hosts have RV on average ∼0.7 lower than that of low-mass hosts; (3) for a subsample of 81 spectroscopically classified SNe there is a significant (>3σ) correlation between the Hubble diagram residuals of red SNe Ia and the host RV that when corrected for reduces scatter by $\sim 13{{\ \rm per\ cent}}$ and the significance of the ‘mass step’ to ∼1σ. These represent independent confirmations of recent predictions based on dust that attempted to explain the puzzling ‘mass step’ and intrinsic scatter (σint) in SN Ia analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Measuring Cosmological Parameters with Type Ia Supernovae in redMaGiC Galaxies

Abstract Current and future cosmological analyses with Type Ia supernovae (SNe Ia) face three critical challenges: (i) measuring the redshifts from the SNe or their host galaxies; (ii) classifying the SNe without spectra; and (iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each of these challenges. In the context of the Dark Energy Survey (DES), we analyze an SN Ia sample with host galaxies in the redMaGiC galaxy catalog, a selection of luminous red galaxies. redMaGiC photo- z estimates are expected to be accurate to σ Δ z /(1+ z ) ∼ 0.02. The DES-5YR photometrically classified SN Ia sample contains approximately 1600 SNe, and 125 of these SNe are in redMaGiC galaxies. We demonstrate that redMaGiC galaxies almost exclusively host SNe Ia, reducing concerns relating to classification uncertainties. With this subsample, we find similar Hubble scatter (to within ∼0.01 mag) using photometric redshifts in place of spectroscopic redshifts. With detailed simulations, we show that the bias due to using redMaGiC photo- z s on the measurement of the dark energy equation of state w is up to Δ w ∼ 0.01–0.02. With real data, we measure a difference in w when using the redMaGiC photo- z s versus the spec- z s of Δ w = 0.005. Finally, we discuss how SNe in redMaGiC galaxies appear to comprise a more standardizable population, due to a weaker relation between color and luminosity ( β ) compared to the DES-3YR population by ∼5 σ . These results establish the feasibility of performing redMaGiC SN cosmology with photometric survey data in the absence of spectroscopic data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Dwarf AGNs from Optical Variability for the Origins of Seeds (DAVOS): insights from the dark energy survey deep fields

ABSTRACT We present a sample of 706, z < 1.5 active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg2. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 d, a 6-yr baseline, and a single-epoch imaging depth of up to g ∼ 24.5. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast }\lt 10^{9.5}\, {\rm M}_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability time-scales (∼ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole (BH) masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalogue, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower BH masses in deep fields, which may be more ‘pristine’ analogues of supermassive BH seeds.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Milky Way Satellite Census. IV. Constraints on Decaying Dark Matter from Observations of Milky Way Satellite Galaxies

Abstract We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetime τ comparable to the age of the universe to a lighter DM particle (with mass splitting ϵ ) and to a dark radiation species. These decays impart a characteristic “kick velocity,” V kick = ϵ c , on the DM daughter particles, significantly depleting the DM content of low-mass subhalos and making them more susceptible to tidal disruption. We fit the suppression of the present-day DDM subhalo mass function (SHMF) as a function of τ and V kick using a suite of high-resolution zoom-in simulations of MW-mass halos, and we validate this model on new DDM simulations of systems specifically chosen to resemble the MW. We implement our DDM SHMF predictions in a forward model that incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk using an empirical model for the galaxy–halo connection. By comparing to the observed MW satellite population, we conservatively exclude DDM models with τ < 18 Gyr (29 Gyr) for V kick = 20 kms −1 (40 kms −1 ) at 95% confidence. These constraints are among the most stringent and robust small-scale structure limits on the DM particle lifetime and strongly disfavor DDM models that have been proposed to alleviate the Hubble and S 8 tensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Optical variability of quasars with 20-yr photometric light curves

ABSTRACT We study the optical gri photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during ∼1998−2020 with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average ∼200 nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10–15 yr baselines and ≲ 100 epochs. We find that the average damping time-scale τDRW continues to rise with increased baseline, reaching a median value of ∼750 d (g band) in the rest frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained τDRW (less than 20 per cent of the baseline), we confirm a weak wavelength dependence of τDRW∝λ0.51 ± 0.20. We further quantify optical variability of these quasars over days to decades time-scales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping time-scales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on time-scales less than ∼ a month for these luminous quasars, and this second break point correlates with the longer DRW damping time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The Observed Evolution of the Stellar Mass–Halo Mass Relation for Brightest Central Galaxies

We quantify evolution in the cluster-scale stellar mass–halo mass (SMHM) relation's parameters using 2323 clusters and brightest central galaxies (BCGs) over the redshift range 0.03 ≤ z ≤ 0.60. The precision on the inferred SMHM parameters is improved by including the magnitude gap (m gap ) between the BCG and fourth-brightest cluster member (M14) as a third parameter in the SMHM relation. At fixed halo mass, accounting for m gap , through a stretch parameter, reduces the SMHM relation's intrinsic scatter. To explore this redshift range, we use clusters, BCGs, and cluster members identified using the Sloan Digital Sky Survey C4 and redMaPPer cluster catalogs and the Dark Energy Survey redMaPPer catalog. Through this joint analysis, we detect no systematic differences in BCG stellar mass, m gap , and cluster mass (inferred from richness) between the data sets. We utilize the Pareto function to quantify each parameter's evolution. We confirm prior findings of negative evolution in the SMHM relation's slope (3.5σ), and detect negative evolution in the stretch parameter (4.0σ) and positive evolution in the offset parameter (5.8σ). This observed evolution, combined with the absence of BCG growth, when stellar mass is measured within 50 kpc, suggests that this evolution results from changes in the cluster's m gap . For this to occur, late-term growth must be in the intracluster light surrounding the BCG. We also compare the observed results to IllustrisTNG 300-1 cosmological hydrodynamic simulations and find modest qualitative agreement. However, the simulations lack the evolutionary features detected in the real data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Dark energy survey year 3 results: Cosmology with peaks using an emulator approach

ABSTRACT We constrain the matter density Ωm and the amplitude of density fluctuations σ8 within the ΛCDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We take into account the uncertainty from the remaining, largely unconstrained ΛCDM parameters (Ωb, ns, and h). We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv ~\sigma _8\sqrt{\Omega _{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68 per cent confidence limits), with a precision of 1.8 per cent, an improvement of 38 per cent compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in S8 is recorded. We find a mild tension of $1.5 \, \sigma$ between our study and the results from Planck 2018, with our analysis yielding a lower S8. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment AIA and S8, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

OzDES reverberation mapping program: Lag recovery reliability for 6-yr C iv analysis

ABSTRACT We present the statistical methods that have been developed to analyse the OzDES reverberation mapping sample. To perform this statistical analysis we have created a suite of customizable simulations that mimic the characteristics of each source in the OzDES sample. These characteristics include: the variability in the photometric and spectroscopic light curves, the measurement uncertainties, and the observational cadence. By simulating the sources in the OzDES sample that contain the C iv emission line, we developed a set of criteria that rank the reliability of a recovered time-lag depending on the agreement between different recovery methods, the magnitude of the uncertainties, and the rate at which false positives were found in the simulations. These criteria were applied to simulated light curves and these results used to estimate the quality of the resulting Radius–Luminosity relation. We grade the results using three quality levels (gold, silver, and bronze). The input slope of the R–L relation was recovered within 1σ for each of the three quality samples, with the gold standard having the lowest dispersion with a recovered a R–L relation slope of 0.454 ± 0.016 with an input slope of 0.47. Future work will apply these methods to the entire OzDES sample of 771 AGN.

79 ASTRONOMY AND ASTROPHYSICS↗