Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “image simulation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Limited-angle effect compensation for respiratory binned cardiac SPECT

Purpose: In cardiac single photon emission computed tomography (SPECT), respiratory-binned study is used to combat the motion blur associated with respiratory motion. However, owing to the variability in respiratory patterns during data acquisition, the acquired data counts can vary significantly both among respiratory bins and among projection angles within individual bins. If not properly accounted for, such variation could lead to artifacts similar to limited-angle effect in image reconstruction. In this work, the authors aim to investigate several reconstruction strategies for compensating the limited-angle effect in respiratory binned data for the purpose of reducing the image artifacts. Methods: The authors first consider a model based correction approach, in which the variation in acquisition time is directly incorporated into the imaging model, such that the data statistics are accurately described among both the projection angles and respiratory bins. Afterward, the authors consider an approximation approach, in which the acquired data are rescaled to accommodate the variation in acquisition time among different projection angles while the imaging model is kept unchanged. In addition, the authors also consider the use of a smoothing prior in reconstruction for suppressing the artifacts associated with limited-angle effect. In our evaluation study, the authors first used Monte Carlo simulated imaging with 4D NCAT phantom wherein the ground truth is known for quantitative comparison. The authors evaluated the accuracy of the reconstructed myocardium using a number of metrics, including regional and overall accuracy of the myocardium, uniformity and spatial resolution of the left ventricle (LV) wall, and detectability of perfusion defect using a channelized Hotelling observer. As a preliminary demonstration, the authors also tested the different approaches on five sets of clinical acquisitions. Results: The quantitative evaluation results show that the three compensation methods could all, but to different extents, reduce the reconstruction artifacts over no compensation. In particular, the model based approach reduced the mean-squared-error of the reconstructed myocardium by as much as 40%. Compared to the approach of data rescaling, the model based approach further improved both the overall and regional accuracy of the myocardium; it also further improved the lesion detectability and the uniformity of the LV wall. When ML reconstruction was used, the model based approach was notably more effective for improving the LV wall; when MAP reconstruction was used, the smoothing prior could reduce the noise level and artifacts with little or no increase in bias, but at the cost of a slight resolution loss of the LV wall. The improvements in image quality by the different compensation methods were also observed in the clinical acquisitions. Conclusions: Compensating for the uneven distribution of acquisition time among both projection angles and respiratory bins can effectively reduce the limited-angle artifacts in respiratory-binned cardiac SPECT reconstruction. Direct incorporation of the time variation into the imaging model together with a smoothing prior in reconstruction can lead to the most improvement in the accuracy of the reconstructed myocardium.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MRI-Based Evaluation of the Vaginal Cuff in Brachytherapy Planning: Are We Missing the Target?

Purpose: Although recurrences and toxicity occur after vaginal cuff (VC) brachytherapy, little is known about dosimetry due to the inability to clearly visualize the VC on computed tomography (CT). T2-weighted (T2W) magnetic resonance imaging (MRI) is superior to CT in this setting, and we hypothesized that it could provide previously unascertainable dosimetric information. Methods and Materials: In a cohort of 32 patients who underwent cylinder-based brachytherapy for endometrial cancer with available MR simulation images, the VC was retrospectively contoured on T2W images, and cases were replanned to treat the upper VC to a dose of 7 Gy/fraction prescribed to 5 mm. Relevant dose-volume parameters for the VC were calculated. Results: T2W MRI identified significant underdosing not observed on CT or T1-weighted imaging. Over two-thirds (69%) of patients had at least 1 cm{sup 3} of VC that received less than 75% of the prescription dose and half (50%) of patients had a least 1 cm{sup 3} of VC that received less than 50% of the prescription dose. The mean minimum point dose to the VC was 2.4 Gy, or 34% of the intended prescription dose (range: 0.53-6.4 Gy). Conclusions: We identified previously unreported VC underdosing in over two-thirds of our patients, with most of these patients having volumes of undistended VC that received less than half of the prescription dose. The maximum dimension was along the craniocaudal axis in some patients or left-right/anterior-posterior axis in others, suggesting that suture material may be restricting access to the vaginal apex and that alternative applicators may be needed when the diameter of the apex is larger than the introitus. Additional follow-up will be needed to determine whether underdosing is associated with isolated VC failure or whether low failure rates across the cohort suggest that some patients are being exposed to excessive dose and unnecessary risk of toxicity.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2010

Collisional Grooming Models of the Kuiper Belt Dust Cloud

We modeled the three-dimensional structure of the Kuiper Belt (KB) dust cloud at four different dust production rates, incorporating both planet-dust interactions and grain-grain collisions using the collisional grooming algorithm. Simulated images of a model with a face-on optical depth of approximately 10 (exp -4) primarily show an azimuthally- symmetric ring at 40-47 AU in submillimeter and infrared wavelengths; this ring is associated with the cold classical KB. For models with lower optical depths (10 (exp -6) and 10 (exp-7)), synthetic infrared images show that the ring widens and a gap opens in the ring at the location of Neptune; this feature is caused by trapping of dust grains in Neptune's mean motion resonances. At low optical depths, a secondary ring also appears associated with the hole cleared in the center of the disk by Saturn. Our simulations, which incorporate 25 different grain sizes, illustrate that grain-grain collisions are important in sculpting today's KB dust, and probably other aspects of the solar system dust complex; collisions erase all signs of azimuthal asymmetry from the submillimeter image of the disk at every dust level we considered. The model images switch from being dominated by resonantly trapped small grains ("transport dominated") to being dominated by the birth ring ("collision dominated") when the optical depth reaches a critical value of r approximately v/c, where v is the local Keplerian speed.

Kuchner, Marc J.↗
DOE OSTI2024

The atomic-level structure and stability of interfaces of Pt nanoparticles in alumina: An experimental and computational evaluation

The atomic-level structure of interfaces between Pt and a transition form of Al2O3 were studied using a combination of electron microscopy and first principles calculations. A model system of Pt nanoprecipitates in Al2O3 were formed in sapphire wafers via high-energy ion implantation of Pt followed by thermal annealing at 1000 °C in air. The Pt nanoparticles took the form of tetrahedra and truncated tetrahedra primarily bound by {111}Pt facets. The high prevalence of these facets motivated the development of density functional theory (DFT) based models of (111)Pt interfaces with six different chemical terminations of ( 2 ¯ 01 ) θ-alumina. The atomic-level structure of the Pt/Al2O3 interfaces was characterized with aberration-corrected scanning transmission electron microscopy (STEM) and the experimental images were compared to STEM image simulations of the DFT models. The model interface with Pt bonded to oxygen-terminated θ-Al2O3, with the Pt located on top of the O and with an underlying layer of octahedral Al, provided the best match to the experimental images. This interfacial termination is also the most stable for the thermal annealing conditions used based on thermodynamic calculations of the interfacial energy as a function of temperature and oxygen partial pressure. This experimentally verified model provides a basis for improving models of Pt/γ-alumina interfaces.

Clauser, A L↗
DOE OSTI2015

Fine precipitation scenarios of AlZnMg(Cu) alloys revealed by advanced atomic-resolution electron microscopy study Part II: Fine precipitation scenarios in AlZnMg(Cu) alloys

Although they are among the most important precipitation-hardened materials for industry applications, the high-strength AlZnMg(Cu) alloys have thus far not yet been understood adequately about their underlying precipitation scenarios in relation with the properties. This is partly due to the fact that the structures of a number of different precipitates involved in the alloys are unknown, and partly due to the complexity that the precipitation behaviors of the alloys may be closely related to the alloy's composition. In Part I of the present study, we have determined all the unknown precipitate structures in the alloys. Here in Part II, using atomic-resolution electron microscopy in association with the first principles energy calculations, we further studied and correlated the phase/structure transformation/evolution among these hardening precipitates in relation with the alloy's composition. It is shown that there are actually two coexisting classes of hardening precipitates in these alloys: the first class includes the η′-precipitates and their early-stage Guinier–Preston (GP-η′) zones; the second class includes the precursors of the equilibrium η-phase (referred to η{sub p}, or η-precursor) and their early-stage Guinier–Preston (GP-η{sub p}) zones. The two coexisting classes of precipitates correspond to two precipitation scenarios. - Highlights: • We determine and verify all the key precipitate structures in AlMgZn(Cu) alloys. • We employ aberration-corrected scanning transmission electron microscopy (STEM). • We use aberration-corrected high-resolution TEM (HRTEM) for the investigations. • We obtain atomic-resolution images of the precipitates and model their structures. • We refine all precipitate structures with quantitative image simulation analysis. • The hardening precipitates in AlZnMg alloys shall be classified into two groups. • Two precipitation scenarios coexist in the alloys. • The precipitation behavior of such an alloy depends on the alloy's composition. • Very detailed phase/structure transformations among the precipitates are revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Dosimetric and deformation effects of image-guided interventions during stereotactic body radiation therapy of the prostate using an endorectal balloon

Purpose: During stereotactic body radiation therapy (SBRT) for the treatment of prostate cancer, an inflatable endorectal balloon (ERB) may be used to reduce motion of the target and reduce the dose to the posterior rectal wall. This work assessed the dosimetric impact of manual interventions on ERB position in patients receiving prostate SBRT and investigated the impact of ERB interventions on prostate shape. Methods: The data of seven consecutive patients receiving SBRT for the treatment of clinical stage T1cN0M0 prostate cancer enrolled in a multi-institutional, IRB-approved trial were analyzed. The SBRT dose was 50 Gy in five fractions to a planning target volume (PTV) that included the prostate (implanted with three fiducial markers) with a 3-5 mm margin. All plans were based on simulation images that included an ERB inflated with 60 cm{sup 3} of air. Daily kilovoltage cone-beam computed tomography (CBCT) imaging was performed to localize the PTV, and an automated fusion with the planning images yielded displacements required for PTV relocalization. When the ERB volume and/or position were judged to yield inaccurate repositioning, manual adjustment (ERB reinflation and/or repositioning) was performed. Based on all 59 CBCT image sets acquired, a deformable registration algorithm was used to determine the dose received by, displacement of, and deformation of the prostate, bladder (BLA), and anterior rectal wall (ARW). This dose tracking methodology was applied to images taken before and after manual adjustment of the ERB (intervention), and the delivered dose was compared to that which would have been delivered in the absence of intervention. Results: Interventions occurred in 24 out of 35 (69%) of the treated fractions. The direct effect of these interventions was an increase in the prostate radiation dose that included 95% of the PTV (D95) from 9.6 {+-} 1.0 to 10.0 {+-} 0.2 Gy (p = 0.06) and an increase in prostate coverage from 94.0% {+-} 8.5% to 97.8% {+-} 1.9% (p = 0.03). Additionally, ERB interventions reduced prostate deformation in the anterior-posterior (AP) direction, reduced errors in the sagittal rotation of the prostate, and increased the similarity in shape of the prostate to the radiotherapy plan (increased Dice coefficient from 0.76 {+-} 0.06 to 0.80 {+-} 0.04, p = 0.01). Postintervention decreases in prostate volume receiving less than the prescribed dose and decreases in the voxel-wise displacement of the prostate, bladder, and anterior rectal wall were observed, which resulted in improved dose-volume histogram (DVH) characteristics. Conclusions: Image-guided interventions in ERB volume and/or position during prostate SBRT were necessary to ensure the delivery of the dose distribution as planned. ERB interventions resulted in reductions in prostate deformations that would have prevented accurate localization of patient anatomy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

archi : pipeline for light curve extraction of CHEOPS background stars

ABSTRACT High precision time-series photometry from space is being used for a number of scientific cases. In this context, the recently launched CHaracterizing ExOPlanet Satellite (CHEOPS) (ESA) mission promises to bring 20 ppm precision over an exposure time of 6 h, when targeting nearby bright stars, having in mind the detailed characterization of exoplanetary systems through transit measurements. However, the official CHEOPS (ESA) mission pipeline only provides photometry for the main target (the central star in the field). In order to explore the potential of CHEOPS photometry for all stars in the field, in this paper, we present archi, an additional open-source pipeline module1 to analyse the background stars present in the image. As archi uses the official data reduction pipeline data as input, it is not meant to be used as an independent tool to process raw CHEOPS data but, instead, to be used as an add-on to the official pipeline. We test archi using CHEOPS simulated images, and show that photometry of background stars in CHEOPS images is only slightly degraded (by a factor of 2–3) with respect to the main target. This opens a potential for the use of CHEOPS to produce photometric time-series of several close-by targets at once, as well as to use different stars in the image to calibrate systematic errors. We also show one clear scientific application where the study of the companion light curve can be important for the understanding of the contamination on the main target.

Silva, André M.↗
DOE OSTI2015

MO-F-CAMPUS-I-03: GPU Accelerated Monte Carlo Technique for Fast Concurrent Image and Dose Simulation

Purpose: To develop an accurate and fast Monte Carlo (MC) method of simulating CT that is capable of correlating dose with image quality using voxelized phantoms. Methods: A realistic voxelized phantom based on patient CT data, XCAT, was used with a GPU accelerated MC code for helical MDCT. Simulations were done with both uniform density organs and with textured organs. The organ doses were validated using previous experimentally validated simulations of the same phantom under the same conditions. Images acquired by tracking photons through the phantom with MC require lengthy computation times due to the large number of photon histories necessary for accurate representation of noise. A substantial speed up of the process was attained by using a low number of photon histories with kernel denoising of the projections from the scattered photons. These FBP reconstructed images were validated against those that were acquired in simulations using many photon histories by ensuring a minimal normalized root mean square error. Results: Organ doses simulated in the XCAT phantom are within 10% of the reference values. Corresponding images attained using projection kernel smoothing were attained with 3 orders of magnitude less computation time compared to a reference simulation using many photon histories. Conclusion: Combining GPU acceleration with kernel denoising of scattered photon projections in MC simulations allows organ dose and corresponding image quality to be attained with reasonable accuracy and substantially reduced computation time than is possible with standard simulation approaches.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Towards emulating cosmic shear data: revisiting the calibration of the shear measurements for the Kilo-Degree Survey

Exploiting the full statistical power of future cosmic shear surveys will necessitate improvements to the accuracy with which the gravitational lensing signal is measured. We present a framework for calibrating shear with image simulations that demonstrates the importance of including realistic correlations between galaxy morphology, size, and more importantly, photometric redshifts. This realism is essential to ensure that selection and shape measurement biases can be calibrated accurately for a tomographic cosmic shear analysis. We emulate Kilo-Degree Survey (KiDS) observations of the COSMOS field using morphological information from Hubble Space Telescope imaging, faithfully reproducing the measured galaxy properties from KiDS observations of the same field. We calibrate our shear measurements from lens fit, and find through a range of sensitivity tests that lens fit is robust and unbiased within the allowed two per cent tolerance of our study. Our results show that the calibration has to be performed by selecting the tomographic samples in the simulations, consistent with the actual cosmic shear analysis, because the joint distributions of galaxy properties are found to vary with redshift. Ignoring this redshift variation could result in misestimating the shear bias by an amount that exceeds the allowed tolerance. To improve the calibration for future cosmic shear analyses, it will also be essential to correctly account for the measurement of photometric redshifts, which requires simulating multi-band observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fast Neutron Scintillator Screens for Neutron Imaging Using a Layered Polymer-Phosphor Architecture

Fast neutrons enable a nondestructive examination of dense, large, and highly attenuating samples due to their lower interaction probability compared to thermal neutrons. However, this also creates a challenge in fast neutron imaging, as the thicker sensors necessary to detect fast neutrons degrade an image’s spatial resolution due to scattering within the sensor and the indeterminate depth of interaction in the sensor. This work explores the advantages of a fast neutron imaging screen with a layered polymer-phosphor screen approach as opposed to a mixed polymer-phosphor screen typically used in fast neutron imaging. Proton recoil is the primary conversion mechanism for fast neutron imaging. Simulations showed that the recoil proton range of typical fast neutrons is approximately 200 µm, however, tests at Idaho National Laboratory revealed that the light output of these screens increased at much greater polymer thicknesses. The NECTAR fast neutron beamline at FRM II was used to test the imaging performance of layered fast neutron imaging screens. Distinguishing between the fast-neutron and γ-ray signals is a major challenge in fast neutron imaging because all fast neutron sources also produce γ-rays. A relative comparison between a control plate and the fast neutron screen was made to distinguish between a γ-ray and fast neutron signals. MCNP modeling quantified the γ-ray and fast neutron contributions to the images measured at NECTAR, which were approximately a 75% γ-ray image.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

SIMULATION OF ASTRONOMICAL IMAGES FROM OPTICAL SURVEY TELESCOPES USING A COMPREHENSIVE PHOTON MONTE CARLO APPROACH

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons s{sup −1}, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to (1) construct kilometer scale phase screens necessary for wide-field telescopes, (2) reproduce atmospheric point-spread function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, (3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and (4) recover system effective area predicted from analytic photometry integrals. This new code, the Photon Simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Structure of a zinc oxide ultra-thin film on Rh(100)

The structural parameters of ultra-thin zinc oxide films on Rh(100) are investigated using low-energy electron diffraction intensity (LEED I–V) curves, scanning tunneling microscopy (STM), and first-principles density functional theory (DFT) calculations. From the analysis of LEED I–V curves and DFT calculations, two optimized models A and B are determined. Their structures are basically similar to the planer h-BN ZnO(0001) structure, although some oxygen atoms protrude from the surface, associated with an in-plane shift of Zn atoms. From a comparison of experimental STM images and simulated STM images, majority and minority structures observed in the STM images represent the two optimized models A and B, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Gamma Imaging and Mapping Advancements

The SRNL developed GrayQb{sup TM} gamma mapping system has been successfully deployed in field situations. Currently, to get the image from the GrayQb{sup TM} the internal Phosphor Storage Plate (PSP) used to capture radiation images must be recovered and scanned by a ScanX Scanner. Upgrades are being implemented which will enable the GrayQB{sup TM} to read the PSP using a camera and red LEDs then transmit results to a remote location. This will greatly reduce personnel interaction, increasing safety as well as eliminating the need for the ScanX. Upgrades have been based off of the MiniMax, a Los Alamos National Laboratory developed PSP reading device constructed with a Nikon DSLR and a red LED flash to obtain the images from the PSP. The electronics in the upgraded GrayQb{sup TM} need to be both functional and environmentally rugged. A PC/104 was chosen for its wider operational temperature range and its physical durability. The Sony α7 ii and the Mightex are the candidate cameras for the upgrades based on test results. PSPs were exposed to gamma radiation sources at SRS's Health Physics Instrumentation Calibration Laboratory (HPICL). PSPs were tested with both Am-241 and Cs-137. Linear intensity increased with dose verifying UV image simulations. Data collected enabled the PSP intensity to be referenced to a total dose. Further work includes design and implementation of a 3D printed housing for the new camera system. The housing will include the camera, light flash, PSP holders, tungsten dual pinhole, PC/104, dosimeter, and a battery or tethering option. Once the housing is finalized further testing at the HPICL facility is planned for early August. During this testing, the final housing design will capture both Am-241 and Cs-137 images.

3D PRINTING↗
DOE OSTI2015

Aberration‐Corrected Transmission Electron Microscopy and In Situ XAFS Structural Characterization of Pt/γ‐Al 2 O 3 Nanoparticles

Abstract Aberration‐corrected (AC) STEM, AC TEM and in situ X‐ray absorption fine structure spectroscopy (XAFS) were used to characterize the Pt clusters present on a 0.35 wt % Pt on γ‐alumina support after reduction in hydrogen at 700 °C. STEM high‐angle annular dark field imaging shows that cluster formation takes place at temperatures up to approximately 350 °C, and this is followed by gradual growth in cluster size for heat treatments in hydrogen up to 700 °C. The STEM data show that after 700 °C reduction the Pt clusters are present in a narrow size distribution centered at 0.88 nm, and using a method involving a redistribution of the Pt atoms using a high electron dosage in the STEM, it is shown that the clusters are present in two‐dimensional morphology. This conclusion is verified using intensity line scans. The in situ extended X‐ray absorption fine structure data are in good agreement with these observations. High‐resolution AC–TEM, which uses a broad coherent electron beam, and can thus offer advantages relative to STEM for structure determination of fine clusters, supported by image simulations of through‐focus series, were used to analyze the structures of Pt particles. The structures determined by using AC–TEM are consistent with STEM and EXAFS data in having a flat two‐dimensional morphology. Comparison of AC–STEM and AC TEM data for the same 700 °C reduced sample suggests that parallel‐beam TEM mode of imaging may be advantageous because of the less pronounced beam‐induced structural rearrangements that occur when imaging with a fine STEM probe.

Sinkler, Wharton↗
DOE OSTI2016

Sci-Fri PM: Radiation Therapy, Planning, Imaging, and Special Techniques - 07: Transitioning from extended-distance total body irradiation to optimized VMAT total marrow irradiation

Purpose: TMI targets only the bone marrow, with the intent of sparing normal tissues. The NSCC has recently implemented a TMI protocol which includes VMAT fields to treat the bone marrow from head to mid-thigh and extended SSD POP fields to treat the lower legs. This work describes the commissioning and initial clinical results of the first reported VMAT TMI treatments in Canada. Methods: Detailed CT simulation, imaging, planning and treatment procedures were developed by a multi-disciplinary team. Patients have 1 cm of bolus over the lower legs and 0.5 cm of bolus around the lower arms. The PTV includes all bone, except mandible, facial bones and hands, with the objective of V(12 Gy) > 90%. Detailed analysis of the influence of field overlap was performed to determine optimal field placement and image-guidance tolerances. Results: PTV coverage was achieved for all cases as V(12 Gy) ranged from 90.4–96.3%. The minimum dose to the PTV, D(99%), ranged from 91.4–97.87% and V(90%Rx=10.8 Gy) ranged from 99.1–100.0%. The lungs, liver and heart had an average D{sub mean} of (7.8±0.3)Gy/(65±2)%, (7.6±0.7)Gy/(63±5)%, and (6.8±0.4)Gy/(56±4)% respectively. Conclusions: Commissioning required input and collaboration from all team members. Transitioning from TBI to TMI requires additional time for contouring, treatment planning, QA, and treatment. Patient benefit can however be seen in the quality of OAR sparing.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Approximate Bayesian computation for forward modeling in cosmology

Bayesian inference is often used in cosmology and astrophysics to derive constraints on model parameters from observations. This approach relies on the ability to compute the likelihood of the data given a choice of model parameters. In many practical situations, the likelihood function may however be unavailable or intractable due to non-gaussian errors, non-linear measurements processes, or complex data formats such as catalogs and maps. In these cases, the simulation of mock data sets can often be made through forward modeling. We discuss how Approximate Bayesian Computation (ABC) can be used in these cases to derive an approximation to the posterior constraints using simulated data sets. This technique relies on the sampling of the parameter set, a distance metric to quantify the difference between the observation and the simulations and summary statistics to compress the information in the data. We first review the principles of ABC and discuss its implementation using a Population Monte-Carlo (PMC) algorithm and the Mahalanobis distance metric. We test the performance of the implementation using a Gaussian toy model. We then apply the ABC technique to the practical case of the calibration of image simulations for wide field cosmological surveys. We find that the ABC analysis is able to provide reliable parameter constraints for this problem and is therefore a promising technique for other applications in cosmology and astrophysics. Our implementation of the ABC PMC method is made available via a public code release.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Roman View of Strong Gravitational Lenses

Galaxy–galaxy strong gravitational lenses can constrain dark matter models and the Lambda cold dark matter cosmological paradigm at subgalactic scales. Currently, there is a dearth of images of these rare systems with high signal-to-noise ratio (SNR) and angular resolution. The Nancy Grace Roman Space Telescope (hereafter Roman), scheduled for launch in late 2026, will play a transformative role in strong-lensing science with its planned wide-field surveys. With its remarkable 0.281 square degree field of view and diffraction-limited angular resolution of ~0$^{''}_.$1, Roman is uniquely suited to characterizing dark matter substructure from a robust population of strong lenses. We present a yield simulation of detectable strong lenses in Roman’s planned High Latitude Wide Area Survey (HLWAS). We simulate a population of galaxy–galaxy strong lenses across cosmic time with cold dark matter subhalo populations, select those detectable in the HLWAS, and generate simulated images accounting for realistic Wide Field Instrument detector effects. For a fiducial case of single 146 s exposures, we predict around 160,000 detectable strong lenses in the HLWAS, of which about 500 will have sufficient SNR to be amenable to detailed substructure characterization. We investigate the effect of variation of the point-spread function across Roman’s field of view on detecting individual subhalos and the suppression of the subhalo mass function at low masses. Our simulation products are available to support strong-lens science with Roman, such as training neural networks and validating dark matter substructure analysis pipelines.

79 ASTRONOMY AND ASTROPHYSICS↗