Engineering topics

Frieman, Joshua

Publications and source records attributed to Frieman, Joshua.

8 records

DOE OSTI2024

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Impact of property covariance on cluster weak lensing scaling relations

ABSTRACT We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, ΔΣ, and the ‘true’ cluster galaxy number count, Ngal, as measured within a spherical volume that is void of projection effects. By quantifying the impact of this covariance on mass calibration, this work reveals a significant source of systematic uncertainty. Using the MDPL2 simulation with galaxies traced by the SAGE semi-analytic model, we measure the intrinsic property covariance between these observables within the three-dimensional vicinity of the cluster, spanning a range of dynamical mass and redshift values relevant for optical cluster surveys. Our results reveal a negative covariance at small radial scales (R ≲ R200c) and a null covariance at large scales (R ≳ R200c) across most mass and redshift bins. We also find that this covariance results in a $2{\!-\!}3~{{\ \rm per\ cent}}$ bias in the halo mass estimates in most bins. Furthermore, by modelling Ngal and ΔΣ as multi-(log)-linear equations of secondary halo properties, we provide a quantitative explanation for the physical origin of the negative covariance at small scales. Specifically, we demonstrate that the Ngal–ΔΣ covariance can be explained by the secondary properties of haloes that probe their formation history. We attribute the difference between our results and the positive bias seen in other works with (mock)-cluster finders to projection effects. These findings highlight the importance of accounting for the covariance between observables in cluster mass estimation, which is crucial for obtaining accurate constraints on cosmological parameters.

Astronomy & Astrophysics↗

DOE OSTI2024

Project Dinos I: A joint lensing–dynamics constraint on the deviation from the power law in the mass profile of massive ellipticals

The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing–dynamics analysis of a large statistical sample with 77 galaxy–galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy–galaxy lenses. We find that the Sloan Lens ACS Survey lens galaxies with a mean redshift of 0.2 are consistent with the power-law profile within 1.1σ (2.8σ) and the Strong Lensing Legacy Survey lens galaxies with a mean redshift of 0.6 are consistent within 0.8σ (2.1σ), for a spatially constant (Osipkov–Merritt) stellar anisotropy profile. We adopted the spatially constant anisotropy profile as our baseline choice based on previous dynamical observables of local ellipticals. However, spatially resolved stellar kinematics of lens galaxies are necessary to differentiate between the two anisotropy models. Future studies will use our lens models to constrain the mass distribution individually in the dark matter and baryonic components.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Modelling galaxy cluster triaxiality in stacked cluster weak lensing analyses

Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis. We test whether redMaPPer selection biases the clusters’ shape and orientation and find that it only biases orientation, preferentially selecting clusters with their major axes oriented along the line of sight. Modelling the richness–mass relation as log-linear, we find that the log-richness amplitude ln (A) is boosted from the lowest to highest orientation bin with a significance of 14σ, while the orientation dependence of the richness-mass slope and intrinsic scatter is minimal. We also find that the weak lensing shear-profile ratios of cluster-associated dark haloes in different orientation bins resemble a ‘bottleneck’ shape that can be quantified with a Cauchy function. We test the correlation of orientation with two other leading systematics in cluster cosmology – miscentering and projection – and find a null correlation. The resulting mass bias predicted from our templates confirms the DES Y1 finding that triaxiality is a leading source of bias in cluster cosmology. However, the richness-dependence of the bias confirms that triaxiality does not fully resolve the tension at low-richness between DES Y1 cluster cosmology and other probes. Our model can be used for quantifying the impact of triaxiality bias on cosmological constraints for upcoming weak lensing surveys of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

LensingETC: A Tool to Optimize Multifilter Imaging Campaigns of Galaxy-scale Strong Lensing Systems

Imaging data is the principal observable required to use galaxy-scale strong lensing in a multitude of applications in extragalactic astrophysics and cosmology. In this paper, we develop Lensing Exposure Time Calculator (LensingETC; to optimize the efficiency of telescope-time usage when planning multifilter imaging campaigns for galaxy-scale strong lenses. This tool simulates realistic data tailored to specified instrument characteristics and then automatically models them to assess the power of the data in constraining lens model parameters. We demonstrate a use case of this tool by optimizing a two-filter observing strategy (in the IR and ultraviolet-visual (UVIS)) within the limited exposure time per system allowed by a Hubble Space Telescope (HST) Snapshot program. We find that higher resolution is more advantageous to gain constraining power on the lensing observables, when there is a trade-off between signal-to-noise ratio and resolution; for example, between the UVIS and IR filters of the HST. We also find that, whereas a point-spread function (PSF) with sub-Nyquist sampling allows the sample mean for a model parameter to be robustly recovered for both galaxy–galaxy and point-source lensing systems, a sub-Nyquist-sampled PSF introduces a larger scatter than a Nyquist-sampled one in the deviation from the ground truth for point-source lens systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Candidate Periodically Variable Quasars from the Dark Energy Survey and the Sloan Digital Sky Survey

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg 2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-year long multi-color ( griz ) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r ~23.5 mag) powered by less massive black holes (with masses ≳ 10 8.5 M ⊙ ) at high redshift for the first time. We find five candidates with significant (at >99.74% single-frequency significance in at least two bands with a global p-value of ~7 × 10 -4 –3× 10 -3 accounting for the look-elsewhere effect) periodicity with observed periods of ~3–5 years (i.e., 1–2 years in rest frame) having ~4–6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of |${\sim }0.8^{+0.5}_{-0.3}$| % or |${\sim }1.1^{+0.7}_{-0.5}$| quasar per deg 2 . Our detection rate is 4–80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. Furthermore, we discuss implications on the future direct detection of low-frequency gravitational waves. Additionally, continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

A ``Tail" of Two Cities: Fatality-based Modeling of COVID-19 Evolution in New York City and Cook County, IL

I describe $SIR$ modeling of the COVID-19 pandemic in two U.S. urban environments, New York City (NYC) and Cook County, IL, from onset through the month of June, 2020. Since testing was not widespread early in the pandemic in the U.S., I do not use data on confirmed cases and rely solely on public fatality data to estimate model parameters. Fits to the first 20 days of data determine a degenerate combination of the basic reproduction number, $R_0$, and the mean time to removal from the infectious population, $\gamma^{-1}$, with $\gamma (R_0-1) = 0.25 (0.21)$ inverse days for NYC (Cook County). Equivalently, the initial doubling time was $t_d = 2.8 (3.4)$ days for NYC (Cook). The early fatality data suggest that both locations had infections in early February. I model the mitigation measures implemented in mid-March in both locations (distancing, quarantine, isolation, etc) via a time-dependent reproduction number $R_t$ that declines monotonically from $R_0$ to a sma ller asymptotic value, with a parameterized functional form. The timing (mid-March) and duration (several days) of the transitions in $R_t$ appear well determined by the data. However, the fatality data determine only a degenerate combination of the parameters $R_0$, the percentage reduction in social contact due to mitigation measures, $X$, and the infection fatality rate (IFR), $f$. With flat priors, based on simulations the NYC model parameters have 95.45\% credible intervals of $R_0 = 3.0-5.4$, $X= 80-99.9\%$ and $f=2-6\%$, with $5-13\%$ of the population asymptotically infected. A strong external prior indicating a lower value of $f$ or of $\gamma^{-1}$ would imply lower values of $R_0$ and $X$ and higher percentage infection of the population. For Cook County, the evolution was qualitatively different: after mitigation measures were implemented, the daily fatality counts reached a plateau for about a month before tailing off. This is consistent with an $SIR$ model that exhibits ``critical slowing-down", in which $R_t$ plateaus ! at a value just above unity. For Cook County, the 95.45\% credible intervals for the model parameters are much broader and shifted downward, $R_0=1.4-4.7$, $X=26-54\%$, and $f=0.1-0.6\%$ with $15-88\%$ of the population asymptotically infected. Despite the apparently lower efficacy of its social contact reduction measures, Cook County has had significantly fewer fatalities per population than NYC, $D_\infty/N=100$ vs. 270 per 100,000. In the model, this is attributed to the lower inferred IFR for Cook; an external prior pointing to similar values of the IFR for the two locations would instead chalk up the difference in $D/N$ to differences in the relative growth rate of the disease. I derive a model-dependent threshold, $X^{crit}_e$, for `safe' re-opening, that is, for easing of contact reduction that would not trigger a second wave; for NYC, the models predict that increasing social contact by more than 20\% from post-mitigation levels will lead to renewed spread, while for Cook County the threshold value is very uncertain, given the parameter degeneracies. The timing of 2nd-wave growth will depend on the amplitude of contact increase relative to $X^{crit}_e$ and on the asymptotic growth rate, and the impact in terms of fatalities will depend on the parameter $f$.

Frieman, Joshua↗

NASA NTRS1994

Can MACHOs probe the shape of the galaxy halo ?

Microlensing searches in our galaxy have recently discovered several candidates in the direction of the Large Magellanic Cloud (LMC). We study the prospects for such searches to yield useful information about the flattening of the Galaxy dark matter halo, using a self-consistent oblate halo model and allowing for the possibility of misalignment between the disk and halo symmetry axes. The microlensing optical depth for the LMC, tau(LMC), depends sensitively on the disk-halo tilt angle in the Milky Way, as does the ratio tau(SMC)/tau(LMC). If the tilt angle is as large as 30 deg, a much larger spread in values for tau(LMC) is consistent with rotation curve constraints than previously thought. Disk-halo tilt and halo flattening do not significantly affect the massive compact halo object (MACHO) masses inferred from event durations.

Frieman, Joshua↗