Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “point spread function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Detection of x ray sources in PROS

The problem of detecting discrete sources in x-ray images has much in common with the problem of automatic source detection at other wavelengths. In all cases, one searches for positive brightness enhancements exceeding a certain threshold, which appear consistent with what one expects for a point source, in the presence of a (possibly) spatially variable background. Multidimensional point spread functions (e.g., dependent on detector position and photon energy) are also common. At the same time, the problem in x-ray astronomy has some unique aspects. For example, for typical x-ray exposures in current or recent observatories, the number of available pixels far exceeds the number of actual x-ray events, so Poisson, rather than Gaussian statistics apply. Further, extended cosmic x-ray sources are common, and one often desires to detect point sources in the vicinity or even within bright, diffuse x-ray emission. Finally, support structures in x-ray detectors often cast sharp shadows in x-ray images making it necessary to detect sources in a region of rapidly varying exposure. We have developed a source detection package within the IRAF/PROS environment which attempts to deal with some of the problems of x-ray source detection. We have patterned our package after the successful Einstein Observatory x-ray source detection programs. However, we have attempted to improve the flexibility and accessibility of the functions and to provide a graphical front-end for the user. Our philosophy has been to use standard IRAF tasks whenever possible for image manipulation and to separate general functions from mission-specific ones. We will report on the current status of the package and discuss future developments, including simulation tasks, to allow the user to assess detection efficiency and source significance, tasks to determine source intensity, and alternative detection algorithms.

Deponte, J.↗
DOE OSTI2023

Beyond the 3rd moment: a practical study of using lensing convergence CDFs for cosmology with DES Y3

ABSTRACT Widefield surveys probe clustered scalar fields – such as galaxy counts, lensing potential, etc. – which are sensitive to different cosmological and astrophysical processes. Constraining such processes depends on the statistics that summarize the field. We explore the cumulative distribution function (CDF) as a summary of the galaxy lensing convergence field. Using a suite of N-body light-cone simulations, we show the CDFs’ constraining power is modestly better than the second and third moments, as CDFs approximately capture information from all moments. We study the practical aspects of applying CDFs to data, using the Dark Energy Survey (DES Y3) data as an example, and compute the impact of different systematics on the CDFs. The contributions from the point spread function and reduced shear approximation are $\lesssim 1~{{\ \rm per\ cent}}$ of the total signal. Source clustering effects and baryon imprints contribute 1–10 per cent. Enforcing scale cuts to limit systematics-driven biases in parameter constraints degrade these constraints a noticeable amount, and this degradation is similar for the CDFs and the moments. We detect correlations between the observed convergence field and the shape noise field at 13σ. The non-Gaussian correlations in the noise field must be modelled accurately to use the CDFs, or other statistics sensitive to all moments, as a rigorous cosmology tool.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Stellar-reddening-based Extinction Maps for Cosmological Applications

Cosmological surveys must correct their observations for the reddening of extragalactic objects by Galactic dust. Existing dust maps, however, have been found to have spatial correlations with the large-scale structure of the Universe. Errors in extinction maps can propagate systematic biases into samples of dereddened extragalactic objects and into cosmological measurements such as correlation functions between foreground lenses and background objects and the primordial non-Gaussianity parameter f NL . Emission-based maps are contaminated by the cosmic infrared background, while maps inferred from stellar reddenings suffer from imperfect removal of quasars and galaxies from stellar catalogs. Thus, stellar-reddening-based maps using catalogs without extragalactic objects offer a promising path to making dust maps with minimal correlations with large-scale structure. We present two high-latitude integrated extinction maps based on stellar reddenings, with a point-spread functions of FWHMs 6'. 1 and 15'. We employ a strict selection of catalog objects to filter out galaxies and quasars and measure the spatial correlation of our extinction maps with extragalactic structure. Our galactic extinction maps have reduced spatial correlation with large-scale structure relative to most existing stellar-reddening-based and emission-based extinction maps.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2019

Innovative Compact Coronagraph Approach for Balloon-Borne Investigation of Temperature and Speed of Electrons in the Corona (BITSE)

We are developing an innovative compact coronagraph for studying the physical conditions in the solar wind acceleration region. This paper presents the new development of the compact coronagraph for the investigation of temperature and speed of electrons in the solar corona. The proposed compact coronagraph is a one stage externally occulted coronagraph without internal occulter or Lyot stop mask. The key of the new idea is to set the inner field cutoff at the External Occulter (EO) much smaller than the specified inner field cutoff. A second occulter on the surface of the detector array removes the remaining diffraction. The occulter on the detector surface functions similar to an internal occulter with the Inner Field of View Cutoff (IFoVC) exactly the same as specified. For BITSE, the desired inner field cutoff is 3 R, but the cutoff at EO is only 1.5 R☉. The diffraction analysis shows that in the sensor plane, the diffraction intensity at the 3 R☉ is not sensitive to the EO cutoff, either at 1.5 R☉ or close to 3 R☉. The advantage of having a smaller EO cutoff is that the vignetting decreased for the Field of View (FoV) near 3 R☉, therefore, the signal increases. Meanwhile, the diffraction of Point Spread Function is much less in the radial direction, which not only increases the image quality around 3 R☉, but also increases the encircled energy and signal to noise ratio. In other words, the data is useful right at 3 R☉! The BITSE optical design and diffraction analysis will be presented in detail. The simulation shows the signal to noise ratio obtained from the diffraction and vignetting data enables corona temperature and speed measurement.

Gong, Qian↗
NASA NTRS2019

Innovative Compact Coronagraph Approach for Balloon-Borne Investigation of Temperature and Speed of Electrons in the Corona (BITSE)

We are developing an innovative compact coronagraph for studying the physical conditions in the solar wind acceleration region. This paper presents the new development of the compact coronagraph for the investigation of temperature and speed of electrons in the solar corona. The proposed compact coronagraph is a one stage externally occulted coronagraph without internal occulter or Lyot stop mask. The key of the new idea is to set the inner field cutoff at the External Occulter (EO) much smaller than the specified inner field cutoff. A second occulter on the surface of the detector array removes the remaining diffraction. The occulter on the detector surface functions similar to an internal occulter with the Inner Field of View Cutoff (IFoVC) exactly the same as specified. For BITSE, the desired inner field cutoff is 3 R, but the cutoff at EO is only 1.5 R☉. The diffraction analysis shows that in the sensor plane, the diffraction intensity at the 3 R☉ is not sensitive to the EO cutoff, either at 1.5 R☉ or close to 3 R☉. The advantage of having a smaller EO cutoff is that the vignetting decreased for the Field of View (FoV) near 3 R☉, therefore, the signal increases. Meanwhile, the diffraction of Point Spread Function is much less in the radial direction, which not only increases the image quality around 3 R☉, but also increases the encircled energy and signal to noise ratio. In other words, the data is useful right at 3 R☉! The BITSE optical design and diffraction analysis will be presented in detail. The simulation shows the signal to noise ratio obtained from the diffraction and vignetting data enables corona temperature and speed measurement.

Gong, Qian↗
NASA NTRS2011

Fine Guidance Sensing for Coronagraphic Observatories

Three options have been developed for Fine Guidance Sensing (FGS) for coronagraphic observatories using a Fine Guidance Camera within a coronagraphic instrument. Coronagraphic observatories require very fine precision pointing in order to image faint objects at very small distances from a target star. The Fine Guidance Camera measures the direction to the target star. The first option, referred to as Spot, was to collect all of the light reflected from a coronagraph occulter onto a focal plane, producing an Airy-type point spread function (PSF). This would allow almost all of the starlight from the central star to be used for centroiding. The second approach, referred to as Punctured Disk, collects the light that bypasses a central obscuration, producing a PSF with a punctured central disk. The final approach, referred to as Lyot, collects light after passing through the occulter at the Lyot stop. The study includes generation of representative images for each option by the science team, followed by an engineering evaluation of a centroiding or a photometric algorithm for each option. After the alignment of the coronagraph to the fine guidance system, a "nulling" point on the FGS focal point is determined by calibration. This alignment is implemented by a fine alignment mechanism that is part of the fine guidance camera selection mirror. If the star images meet the modeling assumptions, and the star "centroid" can be driven to that nulling point, the contrast for the coronagraph will be maximized.

Brugarolas, Paul↗
NASA NTRS2023

JWST NIRCam Defocused Imaging Photometric Stability Performance and How It Can Sense Mirror Tilts

We use JWST NIRCam short-wavelength photometry to capture a transit lightcurve of the exoplanet HAT-P-14 b to assess performance as part of instrument commissioning. The short-wavelength precision is 152 ppm per 27 s integration as measured over the full time series compared to a theoretical limit of 107 ppm, after corrections to spatially correlated 1/f noise. Persistence effects from charge trapping are well fit by an exponential function with short characteristic timescales, settling on the order of 5–15 minutes. The short-wavelength defocused photometry is also uniquely well suited to measure the real-time wave-front error of JWST. Analysis of the images and reconstructed wave-front maps indicates that two different hexagonal primary mirror segments exhibited “tilt events,” where they changed orientation rapidly in less than ∼1.4 s. In some cases, the magnitude and timing of the flux jumps caused by tilt events can be accurately predicted with a telescope model. These tilt events can be sensed by simultaneous longer-wavelength NIRCam grism spectral images alone in the form of changes to the point-spread function, diagnosed from the full width at half maximum. They can also be sensed with the fine guidance sensor instrument from difference images. Tilt events possibly from sudden releases of stress in the backplane structure behind the mirrors were expected during the commissioning period because they were found in ground-based testing. Tilt events have shown signs of decreasing in frequency but have not disappeared completely. The detectors exhibit some minor (less than 1%) deviations from linear behavior in the first few groups of each integration, potentially impacting absolute fluxes and transit depths on bright targets, where only a handful of groups are possible. Overall, the noise is within 50% of the theoretical photon noise and read noise. This bodes well for high-precision measurements of transiting exoplanets and other time variable targets.

Everett Schlawin↗
NASA NTRS1995

Wavelet transform analysis of the small-scale X-ray structure of the cluster Abell 1367

We have developed a new technique based on a wavelet transform analysis to quantify the small-scale (less than a few arcminutes) X-ray structure of clusters of galaxies. We apply this technique to the ROSAT position sensitive proportional counter (PSPC) and Einstein high-resolution imager (HRI) images of the central region of the cluster Abell 1367 to detect sources embedded within the diffuse intracluster medium. In addition to detecting sources and determining their fluxes and positions, we show that the wavelet analysis allows a characterization of the sources extents. In particular, the wavelet scale at which a given source achieves a maximum signal-to-noise ratio in the wavelet images provides an estimate of the angular extent of the source. To account for the widely varying point response of the ROSAT PSPC as a function of off-axis angle requires a quantitative measurement of the source size and a comparison to a calibration derived from the analysis of a Deep Survey image. Therefore, we assume that each source could be described as an isotropic two-dimensional Gaussian and used the wavelet amplitudes, at different scales, to determine the equivalent Gaussian Full Width Half-Maximum (FWHM) (and its uncertainty) appropriate for each source. In our analysis of the ROSAT PSPC image, we detect 31 X-ray sources above the diffuse cluster emission (within a radius of 24 min), 16 of which are apparently associated with cluster galaxies and two with serendipitous, background quasars. We find that the angular extents of 11 sources exceed the nominal width of the PSPC point-spread function. Four of these extended sources were previously detected by Bechtold et al. (1983) as 1 sec scale features using the Einstein HRI. The same wavelet analysis technique was applied to the Einstein HRI image. We detect 28 sources in the HRI image, of which nine are extended. Eight of the extended sources correspond to sources previously detected by Bechtold et al. Overall, using both the PSPC and the HRI observations, we detect 16 extended features, of which nine have galaxies coincided with the X-ray-measured positions (within the positional error circles). These extended sources have luminosities lying in the range (3 - 30) x 10(exp 40) ergs/s and gas masses of approximately (1 - 30) x 10(exp 9) solar mass, if the X-rays are of thermal origin. We confirm the presence of extended features in A1367 first reported by Bechtold et al. (1983). The nature of these systems remains uncertain. The luminosities are large if the emission is attributed to single galaxies, and several of the extended features have no associated galaxy counterparts. The extended features may be associated with galaxy groups, as suggested by Canizares, Fabbiano, & Trinchieri (1987), although the number required is large.

Grebeney, S. A.↗
DOE OSTI2022

Joint Survey Processing. I. Compact Oddballs in the COSMOS Field—Low-luminosity Quasars at z > 6?

The faint-end slope of the quasar luminosity function at z ~ 6 and its implication on the role of quasars in reionizing the intergalactic medium at early times has been an outstanding problem for some time. The identification of faint high-redshift quasars with luminosities of <10 44.5 erg s –1 is challenging. They are rare (few per square degree), and the separation of these unresolved quasars from late-type stars and compact star-forming galaxies is difficult from ground-based observations alone. In addition, source confusion becomes significant at >25 mag, with ~30% of sources having their flux contaminated by foreground objects when the seeing resolution is ~0''.7. We mitigate these issues by performing a pixel-level joint processing of ground and space-based data from Subaru/Hyper-SuprimeCam (HSC) and Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). We create a deconfused catalog over the 1.64 deg 2 of the COSMOS field, after accounting for spatial varying point-spread functions and astrometric differences between the two data sets. We identify twelve low-luminosity (M UV ~ –21 mag) z > 6 quasar candidates through (i) their red color measured between ACS/F814W and HSC/i band and (ii) their compactness in the space-based data. Nondetections of our candidates in Hubble DASH data argues against contamination from late-type stars. Our constraints on the faint end of the quasar luminosity function at z ~ 6.4 suggest a negligibly small contribution to reionization compared to the star-forming galaxy population. The confirmation of our candidates and the evolution of number density with redshift could provide better insights into how supermassive galaxies grew in the first billion years of cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Modeling and Correcting the Time-Dependent ACS PSF

The ability to accurately measure the shapes of faint objects in images taken with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) depends upon detailed knowledge of the Point Spread Function (PSF). We show that thermal fluctuations cause the PSF of the ACS Wide Field Camera (WFC) to vary over time. We describe a modified version of the TinyTim PSF modeling software to create artificial grids of stars across the ACS field of view at a range of telescope focus values. These models closely resemble the stars in real ACS images. Using 10 bright stars in a real image, we have been able to measure HST s apparent focus at the time of the exposure. TinyTim can then be used to model the PSF at any position on the ACS field of view. This obviates the need for images of dense stellar fields at different focus values, or interpolation between the few observed stars. We show that residual differences between our TinyTim models and real data are likely due to the effects of Charge Transfer Efficiency (CTE) degradation. Furthermore, we discuss stochastic noise that is added to the shape of point sources when distortion is removed, and we present MultiDrizzle parameters that are optimal for weak lensing science. Specifically, we find that reducing the MultiDrizzle output pixel scale and choosing a Gaussian kernel significantly stabilizes the resulting PSF after image combination, while still eliminating cosmic rays/bad pixels, and correcting the large geometric distortion in the ACS. We discuss future plans, which include more detailed study of the effects of CTE degradation on object shapes and releasing our TinyTim models to the astronomical community.

Rhodes, Jason↗
NASA NTRS2026

In-Flight Calibration of ESA Hera’s Hyperscout-H Imager

ESA’s Hera space mission is on its way to the mission target, the binary asteroid (65803) Didymos. HyperScout-H, one of the instruments onboard Hera, is a hyperspectral imager operating in the visible and near-infrared regions between 0.65 and 0.95 μm. HyperScout-H will enable a detailed assessment of the composition of both objects, Didymos and its satellite Dimorphos, the characterization of space weathering effects, and the possible presence of exogenous material on their surfaces. To monitor instrument functionality, calibration exposures are acquired regularly. This article describes the in-flight calibrations carried out for HyperScout-H during the commissioning and cruise phases. Bias and dark exposures, as well as stellar field observations, were acquired several times after launch. We update the calibration data and monitor instrument performance in the space environment. In addition, images of Earth and Moon were acquired from distances of 1.5 × 10 6 to 2.0 × 10 6 km, and Mars and its satellite were imaged during the flyby. In five images, the surface of Mars fills the entire field of view, enabling cross-validation of HyperScout-H results with those reported by other Mars missions. We characterize the detector under in-flight operational conditions. The calibration data indicate that the bias pattern is stable, the dark current remains negligible for short exposures, and the detector response is highly linear. We quantify the field-of-view alignment and geometric distortion, and evaluate the point spread function based on the stellar field observations. Stellar observations and Mars swing-by data provide updated radiometric calibration constants, suggesting that in-flight conditions have slightly modified the detector’s spectral response. In-flight calibrations are essential to ensure data quality and reliability.

Image Processing↗
DOE OSTI2025

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

A plasma-based non-intrusive point source for acoustic beamforming applications

A laser-generated plasma acoustic point source is used to directly measure the point spread function (PSF) of a microphone phased array. In beamforming analysis of microphone phased array data, the true acoustic field is convolved with the array's PSF. By directly measuring the PSF, corrections to the array analysis can be computed and applied. The acoustic source is measured in an open-jet aeroacoustic facility to evaluate the effects of sampling rate, microphone installation, source shift, reflections, shear layer refraction and model presence. Results show that measurements exhibit behavior consistent with theory with regard to source shift and shear layer refraction. Application of a measured PSF in beamforming analysis shows the process provides an effective in-situ method for array calibration both with and without flow and allows for corrections to incorporate reflections and scattering. The technique improves the agreement of beamforming results with the true spectrum of a known source, especially in the presence of reflections.

phased array↗
DOE OSTI2024

Mitigation of the Brighter-fatter Effect in the LSST Camera

Abstract Thick, fully depleted charge-coupled devices are known to exhibit nonlinear behavior at high signal levels due to the dynamic behavior of charges collecting in the potential wells of pixels, called the brighter-fatter effect (BFE). The effect results in distorted images of bright calibration stars, creating a flux-dependent point-spread function that if left unmitigated, could make up a large fraction of the error budget in Stage IV weak-lensing (WL) surveys such as the Legacy Survey of Space and Time (LSST). In this paper, we analyze image measurements of flat fields and artificial stars taken at different illumination levels with the LSST Camera (LSSTCam) at SLAC National Accelerator Laboratory in order to quantify this effect in the LSSTCam before and after a previously introduced correction technique. We observe that the BFE evolves anisotropically as a function of flux due to higher-order BFEs, which violates the fundamental assumption of this correction method. We then introduce a new method based on a physically motivated model to account for these higher-order terms in the correction, and then we test the modified correction on both data sets. We find that the new method corrects the effect in flat fields better than it corrects the effect in artificial stars, which we suggest is the result of sub-pixel physics not included in this correction model. We use these results to define a new metric for the full-well capacity of our sensors and advise image processing strategies to further limit the impact of the effect on LSST WL science pathways.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Integrated Strategies for Overcoming Resolution Limits in Electron Beam Lithography of Chemically Amplified Resists

Electron beam lithography (EBL) of chemically amplified resists (CARs) faces fundamental challenges, including stochastic electron scattering and acid diffusion, that limit resolution and reproducibility. Using SU-8 as a model CAR, this study systematically investigated complementary strategies to address these challenges, combining multipass exposure, proximity effect correction (PEC) with midrange correction factors, base quencher incorporation, and post-exposure bake (PEB) suppression. Monte Carlo simulations and calibrated PEC modeling revealed that extending the point spread function to include a midrange scattering component significantly improved critical dimension (CD) control across varying pattern densities, correcting deviations that conventional two-term PEC failed to capture. Multipass exposure, particularly 4-pass writing with a 25% offset, redistributed the dose to average stochastic beam and scattering fluctuations, reducing line-width roughness by more than 50% and yielding more uniform nanoscale features. Photoacid confinement was investigated by adding urea as a base quencher, which successfully reduced acid diffusion but introduced substantial sensitivity penalties without improving ultimate resolution or Z-factor performance, underscoring the trade-offs of chemical versus physical confinement. Suppressing PEB most directly minimized acid diffusion, resulting in improved Z-factors and reproducible 30 nm half-pitch dense line/space patterns. Overall, these results demonstrated that PEC with midrange correction, multipass strategies, quencher additives, and PEB-free processing addresses different aspects of the EBL process window and that their integration provides a comprehensive framework for managing stochastic scattering, diffusion, and chemical amplification effects. This framework advances dense nanoscale patterning in CARs and establishes guiding principles for optimizing resist design and process strategies in high-resolution EBL and potentially other advanced lithographies, such as extreme ultraviolet (EUV) lithography.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Strong Upper-Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

We use earthquake-based adjoint tomography to invert for three-dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land-based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L-BFGS optimization algorithm, improving data fit and introducing P- and S-wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. Here, the inversion reveals increased along-strike heterogeneity on the margin. In Cook Strait we observe a low-velocity zone interpreted as deep sedimentary basins. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone.

58 GEOSCIENCES↗
DOE OSTI2022

TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellar kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within 1.2σ, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a ~4σ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within ~0.6σ<, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

TDCOSMO. X. Automated modeling of nine strongly lensed quasars and comparison between lens-modeling software

When strong gravitational lenses are to be used as an astrophysical or cosmological probe, models of their mass distributions are often needed. We present a new, time-efficient automation code for the uniform modeling of strongly lensed quasars with GLEE, a lens-modeling software for multiband data. By using the observed positions of the lensed quasars and the spatially extended surface brightness distribution of the host galaxy of the lensed quasar, we obtain a model of the mass distribution of the lens galaxy. We applied this uniform modeling pipeline to a sample of nine strongly lensed quasars for which images were obtained with the Wide Field Camera 3 of the Hubble Space Telescope. The models show well-reconstructed light components and a good alignment between mass and light centroids in most cases. We find that the automated modeling code significantly reduces the input time during the modeling process for the user. The time for preparing the required input files is reduced by a factor of 3 from ~3 h to about one hour. The active input time during the modeling process for the user is reduced by a factor of 10 from ~ 10 h to about one hour per lens system. This automated uniform modeling pipeline can efficiently produce uniform models of extensive lens-system samples that can be used for further cosmological analysis. A blind test that compared our results with those of an independent automated modeling pipeline based on the modeling software Lenstronomy revealed important lessons. Quantities such as Einstein radius, astrometry, mass flattening, and position angle are generally robustly determined. Other quantities, such as the radial slope of the mass density profile and predicted time delays, depend crucially on the quality of the data and on the accuracy with which the point spread function is reconstructed. Better data and/or a more detailed analysis are necessary to elevate our automated models to cosmography grade. Nevertheless, our pipeline enables the quick selection of lenses for follow-up and further modeling, which significantly speeds up the construction of cosmography-grade models. This important step forward will help us to take advantage of the increase in the number of lenses that is expected in the coming decade, which is an increase of several orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗