Search NASA⌕ Search

SEARCH · Search NASA

Results for “reliability estimation”

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 91 records · Page 5

Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT

Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.

58 GEOSCIENCES↗

Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes

Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.

58 GEOSCIENCES↗

Rapid Coal-Ash Characterization using Geophysical Methods & Machine Learning

Coal combustion products (CCP) are challenging to delineate in heterogeneous field settings. Conventional methods (test pits, coring, and laboratory analyses) are labor-intensive, slow, invasive, and provide sparse spatial coverage. This study evaluates whether rapid non-invasive geophysical screening methods—induced polarization (IP), magnetic susceptibility, and nuclear magnetic resonance (NMR) —combined with surface colorimetry (RGB_24), can discriminate CCP-soil mixtures and provide reliable estimates of CCP content. Laboratory measurements were collected on five CCP-soil mixtures (series) and modeled using (i) a linear baseline, (ii) a calibrated non-linear (power-mean) model, and (iii) a machine-learning (ML) Random Forest approach, with validation via leave-one-series-out and site-specific tests. Across the five series, individual signals—particularly IP and magnetic susceptibility—were strongly predictive of ash content but were consistently outperformed by combined models. The pooled calibrated non-linear and ML models captured the observed non-linearity and achieved high accuracy and precision, improving on linear fits. Colorimetry showed the weakest direct relationship with ash content for the tested samples but improved performance when included in multi-signal models. At pre-selected 3.5% decision threshold, calibrated and ML approaches yielded near-perfect classification (Matthews correlation coefficient ˜ 1), suggesting strong practical operability for field screening. Additionally, field-analog tests highlighted the role of endmembers—accuracy declined without access to end-member measurements but was largely recovered by collecting a minimal labeled pair for local recalibration. With end members, accuracy remained high. Globally trained models performed well on three operational unknowns; however, series-specific refits provided the most accurate predictions. Overall, these results highlight the potential of combining rapid geophysics and minimal local calibration for improved coal-ash delineation.

Peshtani, Klaudio↗

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗

Benchmark Data Set of Crystalline Organic Semiconductors

This work reports a Benchmark Data set of Crystalline Organic Semiconductors to test calculations of the structural and electronic properties of these materials in the solid state. The data set contains 67 crystals consisting of mostly rigid molecules with a single dominant conformer, covering the majority of known structural types. The experimental crystal structure is available for the entire data set, whereas zero-temperature unit cell volume can be reliably estimated for a subset of 28 crystals. Using this subset, we benchmark r 2 SCAN-D3 and PBE-D3 density functionals. Then, for the entire data set, we benchmark approximate density functional theory (DFT) methods, including GFN1-xTB and DFTB3(3ob-3-1), with various dispersion corrections against r 2 SCAN-D3. Our results show that r 2 SCAN-D3 geometries are accurate within a few percent, which is comparable to the statistical uncertainty of experimental data at a fixed temperature, but the unit cell volume is systematically underestimated by 2% on average. The several times faster PBE-D3 provides an unbiased estimate of the volume for all systems except for molecules with highly polar bonds, for which the volume is substantially overestimated in correlation with the underestimation of atomic charges. Considered approximate DFT methods are orders of magnitude faster and provide qualitatively correct but overcompressed crystal structures unless the dispersion corrections are fitted by unit cell volume.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduced‐Order Probabilistic Emulation of Physics‐Based Ring Current Models: Application to RAM‐SCB Particle Flux

Abstract In this work, we address the computational challenge of large‐scale physics‐based simulation models for the ring current. Reduced computational cost allows for significantly faster than real‐time forecasting, enhancing our ability to predict and respond to dynamic changes in the ring current, valuable for space weather monitoring and mitigation efforts. Additionally, it can also be used for a comprehensive investigation of the system. Thus, we aim to create an emulator for the Ring current‐Atmosphere interactions Model with Self‐Consistent magnetic field (RAM‐SCB) particle flux that not only improves efficiency but also facilitates forecasting with reliable estimates of prediction uncertainties. The probabilistic emulator is built upon the methodology developed by Licata and Mehta (2023), https://doi.org/10.1029/2022sw003345 . A novel discrete sampling is used to identify 30 simulation periods over 20 years of solar and geomagnetic activity. Focusing on a subset of particle flux, we use Principal Component Analysis for dimensionality reduction and Long Short‐Term Memory (LSTM) neural networks to perform dynamic modeling. Hyperparameter space was explored extensively resulting in about 5% median symmetric accuracy across all data sets for one‐step dynamic prediction. Using a hierarchical ensemble of LSTMs, we have developed a reduced‐order probabilistic emulator (ROPE) tailored for time‐series forecasting of particle flux in the ring current. This ROPE offers accurate predictions of omnidirectional flux at a single energy with no pitch angle information, providing robust predictions on the test set with an error score below 11% and calibration scores under 8% with bias under 2% providing a significant speed up as compared to the full RAM‐SCB run.

79 ASTRONOMY AND ASTROPHYSICS↗

Instantaneous tracking of earthquake growth with elastogravity signals

Rapid and reliable estimation of large earthquake magnitude (above 8) is key to mitigating the risks associated with strong shaking and tsunamis. Standard early warning systems based on seismic waves fail to rapidly estimate the size of such large earthquakes. Geodesy-based approaches provide better estimations, but are also subject to large uncertainties and latency associated with the slowness of seismic waves. Recently discovered speed-of-light prompt elastogravity signals (PEGS) have raised hopes that these limitations may be overcome, but have not been tested for operational early warning. Here we show that PEGS can be used in real time to track earthquake growth instantaneously after the event reaches a certain magnitude. We develop a deep learning model that leverages the information carried by PEGS recorded by regional broadband seismometers in Japan before the arrival of seismic waves. After training on a database of synthetic waveforms augmented with empirical noise, we show that the algorithm can instantaneously track an earthquake source time function on real data. Our model unlocks ‘true real-time’ access to the rupture evolution of large earthquakes using a portion of seismograms that is routinely treated as noise, and can be immediately transformative for tsunami early warning.

58 GEOSCIENCES↗

Moment expansion of polarized dust SED: A new path towards capturing the CMB B -modes with LiteBIRD

Accurate characterization of the polarized dust emission from our Galaxy will be decisive in the quest for the cosmic microwave background (CMB) primordial B-modes. An incomplete modeling of its potentially complex spectral properties could lead to biases in the CMB polarization analyses and to a spurious measurement of the tensor-to-scalar ratio r. It is particularly crucial for future surveys like the LiteBIRD satellite, the goal of which is to constrain the faint primordial signal leftover by inflation with an accuracy on the tensor-to-scalar ratio r of the order of 10 –3 . Variations of the dust properties along and between lines of sight lead to unavoidable distortions of the spectral energy distribution (SED) that cannot be easily anticipated by standard component-separation methods. This issue can be tackled using a moment expansion of the dust SED, an innovative parametrization method imposing minimal assumptions on the sky complexity. In the present paper, we apply this formalism to the B-mode cross-angular power spectra computed from simulated LiteBIRD polarization data at frequencies between 100 and 402 GHz that contain CMB, dust, and instrumental noise. The spatial variation of the dust spectral parameters (spectral index β and temperature T) in our simulations lead to significant biases on r (~21 σ r ) if not properly taken into account. Performing the moment expansion in β, as in previous studies, reduces the bias but does not lead to sufficiently reliable estimates of r. We introduce, for the first time, the expansion of the cross-angular power spectra SED in both β and T, showing that, at the sensitivity of LiteBIRD, the SED complexity due to temperature variations needs to be taken into account in order to prevent analysis biases on r. Thanks to this expansion, and despite the existing correlations between some of the dust moments and the CMB signal responsible for a rise in the error on r, we can measure an unbiased value of the tensor-to-scalar ratio with a dispersion as low as σ r = 8.8 × 10 –4 .

79 ASTRONOMY AND ASTROPHYSICS↗

High-redshift quasars along the main sequence

The 4D Eigenvector 1 empirical formalism (4DE1) and its main sequence (MS) for quasars has emerged as a powerful tool for organising the diversity among quasar populations, as several key observational measures and physical parameters are systematically changing along it. Trends revealed by 4DE1 are very well established to explain all the diverse characteristics seen in low-redshift quasar samples. Nevertheless, the situation is far less clear when dealing with high-luminosity and high-redshift sources. Here, we aim to evaluate the behaviour of our sample of 22 quasars at high redshift (2.2 ≤ z ≤ 3.7) and high luminosity (47.39 ≤ L bol ≤ 48.36) in the context of the 4DE1. Our approach involves studying quasar physics through a spectroscopic exploration of UV and optical emission line diagnostics. We used new observations from the ISAAC instrument at ESO-VLT and primarily from the SDSS to cover the optical and the UV rest-frames, respectively. The emission lines were characterised both via a quantitative parametrisation of the line profiles and a decomposition of the emission line profiles using multicomponent fitting routines. We provide spectrophotometric properties and line profile measurements for Hβ+[O III]λλ4959,5007, as well as Si IVλ1397+O IV]λ1402, C IVλ1549+He IIλ1640, and the 1900 Å blend (including Al IIIλ1860, Si III]λ1892, and C III]λ1909). For six out of the 22 objects, a significantly blueshifted component on the Hβ profile is present. In 14 out of 22 cases, an Hβ outflowing component associated with [O III] is detected. The majority of [O III]λλ4959,5007 emission line profiles show blueshifted velocities higher than 250 km s -1 . We find extremely broad [O III]λλ4959,5007 emission that is comparable to the width of Hβ broad profile in some highly accreting quasars. The [O III]λλ4959,5007 and C IVλ1549 blueshifts show very high amplitudes and a high degree of correlation. The line widths and shifts are correlated for both [O III]λλ4959,5007 and C IVλ1549, suggesting that emission from outflowing gas is providing a substantial broadening effect to both lines. Otherwise, the links between C IVλ1549 centroid velocity at half intensity (c(1/2)), Eddington ratio (L/L Edd ), and bolometric luminosity are found to be in agreement with previous studies of high-luminosity quasars. Our analysis suggests that the behaviour of quasars of very high luminosity all along the main sequence is strongly affected by powerful outflows involving a broad range of spatial scales. The main sequence correlations remain valid at high redshift and high luminosity even if a systematic increase in line width is observed. Scaling laws based on UV Al IIIλ1860 and Hβ emission lines are equally reliable estimators of M BH .

79 ASTRONOMY AND ASTROPHYSICS↗

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of tungsten filament in gas discharge in H - ion source

Sources of negative ions such as H - are essential elements of proton accelerators and tokamaks. They have limited lifetime. The replacement of an ion source is a costly process causing delays in the operation of the entire machine. The hot tungsten filament is the key element that limits the lifetime of the H - ion source at the Los Alamos Neutron Science Center (LANSCE) facility. An accurate model, which describes the filament physics and provides a reliable estimate on its lifetime, may significantly improve the operation of the facility. In this report we describe such a model, including a comprehensive list of relevant physical processes from first principles. The model can be used to describe different regimes of operation using diagnostics data as input parameters. It has been benchmarked against the data collected during the production cycle at LANSCE and shows good agreement with experimental data.

47 OTHER INSTRUMENTATION↗

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elimination of remnant phases in low-temperature growth of wurtzite ScAlN by molecular-beam epitaxy

We report growth of wurtzite Sc x Al 1-x N (x < 0.23) by plasma-assisted molecular-beam epitaxy on c-plane GaN at high temperatures significantly alters the extracted lattice constants of the material due to defects likely associated with remnant phases. In contrast, ScAlN grown below a composition-dependent threshold temperature exhibits uniform alloy distribution, reduced defect density, and atomic-step surface morphology. The c-plane lattice constant of this low-temperature ScAlN varies with composition as expected from previous theoretical calculations and can be used to reliably estimate alloy composition. Moreover, lattice-matched Sc 0.18 Al0 .82 N/GaN multi-quantum wells grown under these conditions display strong and narrow near-infrared intersubband absorption lines that confirm advantageous optical and electronic properties.

36 MATERIALS SCIENCE↗

Cross-correlation image analysis for real-time single particle tracking

Accurately measuring the translations of objects between images is essential in many fields, including biology, medicine, chemistry, and physics. One important application is tracking one or more particles by measuring their apparent displacements in a series of images. Popular methods, such as the center of mass, often require idealized scenarios to reach the shot noise limit of particle tracking and, therefore, are not generally applicable to multiple image types. More general methods, such as maximum likelihood estimation, reliably approach the shot noise limit, but are too computationally intense for use in real-time applications. These limitations are significant, as real-time, shot-noise-limited particle tracking is of paramount importance for feedback control systems. To fill this gap, we introduce a new cross-correlation-based algorithm that approaches shot-noise-limited displacement detection and a graphics processing unit-based implementation for real-time image analysis of a single particle.

Instruments & Instrumentation↗

Water dissociation at the water–rutile TiO 2 (110) interface from ab initio-based deep neural network simulations

The interaction of water with TiO 2 surfaces is of crucial importance in various scientific fields and applications, from photocatalysis for hydrogen production and the photooxidation of organic pollutants to self-cleaning surfaces and bio-medical devices. In particular, the equilibrium fraction of water dissociation at the TiO 2 –water interface has a critical role in the surface chemistry of TiO 2 , but is difficult to determine both experimentally and computationally. Among TiO 2 surfaces, rutile TiO 2 (110) is of special interest as the most abundant surface of TiO 2 ’s stable rutile phase. While surface-science studies have provided detailed information on the interaction of rutile TiO 2 (110) with gas-phase water, much less is known about the TiO 2 (110)–water interface, which is more relevant to many applications. In this work, we characterize the structure of the aqueous TiO 2 (110) interface using nanosecond timescale molecular dynamics simulations with ab initio-based deep neural network potentials that accurately describe water/TiO 2 (110) interactions over a wide range of water coverages. Simulations on TiO 2 (110) slab models of increasing thickness provide insight into the dynamic equilibrium between molecular and dissociated adsorbed water at the interface and allow us to obtain a reliable estimate of the equilibrium fraction of water dissociation. We find a dissociation fraction of 22 ± 6% with an associated average hydroxyl lifetime of 7.6 ± 1.8 ns. These quantities are both much larger than corresponding estimates for the aqueous anatase TiO 2 (101) interface, consistent with the higher water photooxidation activity that is observed for rutile relative to anatase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetic structure, excitations and short-range order in honeycomb Na 2 Ni 2 TeO 6

Na 2 Ni 2 TeO 6 has a layered hexagonal structure with a honeycomb lattice constituted by Ni 2+ and a chiral charge distribution of Na + that resides between the Ni layers. In the present work, the antiferromagnetic transition temperature of Na 2 Ni 2 TeO 6 is confirmed at TN˜ 27 K, and further, it is found to be robust up to 8 T magnetic field and 1.2 GPa external pressure; and, without any frequency-dependence. Slight deviations from nominal Na-content (up to 5%) does not seem to influence the magnetic transition temperature, TN. Isothermal magnetization curves remain almost linear up to 13 T. Our analysis of neutron diffraction data shows that the magnetic structure of Na 2 Ni 2 TeO 6 is faithfully described by a model consisting of two phases described by the commensurate wave vectors vec kc, (0.5 0 0) and (0.5 0 0.5), with an additional short-range order component incorporated in to the latter phase. Consequently, a zig-zag long-range ordered magnetic phase of Ni 2+ results in the compound, mixed with a short-range ordered phase, which is supported by our specific heat data. Theoretical computations based on density functional theory (DFT) predict predominantly in-plane magnetic exchange interactions that conform to a J1-J2-J3 model with a strong J3 term. The computationally predicted parameters lead to a reliable estimate for TN and the experimentally observed zig-zag magnetic structure. A spin wave excitation in Na 2 Ni 2 TeO 6 at E ˜ 5 meV at T = 5 K is mapped out through inelastic neutron scattering experiments, which is reproduced by linear spin wave theory calculations using the J values from our computations. Our specific heat data and inelastic neutron scattering data strongly indicate the presence of short-range spin correlations, at T > TN, stemming from incipient antiferromagnetic clusters.

36 MATERIALS SCIENCE↗

Femtosecond pulse parameter estimation from photoelectron momenta using machine learning

Abstract Deep learning models have provided huge interpretation power for image-like data. Specifically, convolutional neural networks (CNNs) have demonstrated incredible acuity for tasks such as feature extraction or parameter estimation. Here we test CNNs on strong-field ionization photoelectron spectra, training on theoretical data sets to ‘invert’ experimental data. Pulse characterization is used as a ‘testing ground’, specifically we retrieve the laser intensity, where ‘traditional’ measurements typically lead to 20% uncertainty. We report on crucial data augmentation techniques required to successfully train on theoretical data and return consistent results from experiments, including accounting for detector saturation. The same procedure can be repeated to apply CNNs in a range of scenarios for strong-field ionization. Using a predictive uncertainty estimation, reliable laser intensity uncertainties of a few percent can be extracted, which are consistently lower than those given by traditional techniques. Using interpretability methods can reveal parts of the distribution that are most sensitive to laser intensity, which can be directly associated with holographic interferences. The CNNs employed provide an accurate and convenient ways to extract parameters, and represent a novel interpretational tool for strong-field ionization spectra.

Szołdra, Tomasz (ORCID:0000000228970506)↗

Deep inference of simulated strong lenses in ground-based surveys

The large number of strong lenses discoverable in future astronomical surveys will likely enhance the value of strong gravitational lensing as a cosmic probe of dark energy and dark matter. However, leveraging the increased statistical power of such large samples will require further development of automated lens modeling techniques. We show that deep learning and simulation-based inference (SBI) methods produce informative and reliable estimates of parameter posteriors for strong lensing systems in ground-based surveys. We present the examination and comparison of two approaches to lens parameter estimation for strong galaxy-galaxy lenses — Neural Posterior Estimation (NPE) and Bayesian Neural Networks (BNNs). We perform inference on 1-, 5-, and 12-parameter lens models for ground-based imaging data that mimics the Dark Energy Survey (DES). We find that NPE outperforms BNNs, producing posterior distributions that are more accurate, precise, and well-calibrated for most parameters. For the 12-parameter NPE model, the calibration is consistently within <10% of optimal calibration for all parameters, while the BNN is rarely within 20% of optimal calibration for any of the parameters. Similarly, residuals for most of the parameters are smaller (by up to an order of magnitude) with the NPE model than the BNN model. This work takes important steps in the systematic comparison of methods for different levels of model complexity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗