SEARCH · Engineering Papers

Results for “normal measurement error”

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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Differential, integral, and momentum-transfer cross sections for elastic electron scattering by neon - 5 to 100 eV

Relative elastic-scattering differential cross sections were measured in the 5-100-eV impact energy and 10-145 deg angular ranges. Normalization of these cross sections was achieved by utilizing accurate total electron-scattering cross sections. A phase-shift analysis of the angular distributions in terms of real phase shifts has been carried out. From the differential cross sections, momentum-transfer cross sections were obtained and the values of the critical energy and angle were established (associated with the lowest value of the differential cross section) as 62.5 + or - 2.5 eV and 101.7 deg + or - 1.5 deg, respectively. The present phase shifts, the critical parameters, and differential, integral, and momentum-transfer cross sections are compared to previous experimental and theoretical results. The error associated with the present data is about 10 percent.

DOE OSTI2020

An improved method for direct incident solar radiation calculation from hourly solar insolation data in building energy simulation

Solar radiation considerably influences the energy consumption of buildings and the power production of building integrated photovoltaic (BIPV) systems. Hourly solar insolation (Wh/m 2 ), represented as the amount of solar irradiance collected on the ground during a 1-h period, is the most common solar radiation data available and widely used in weather files applied in building energy modeling programs (BEMPs). Because the solar beam and position vary over time, the use of hourly insolation data as the input might result in errors in the estimation of the direct incident solar radiation on a particular surface. In this study, methods used in BEMPs for direct incident solar radiation calculations are first analyzed, and an improved method adopting a new algorithm for estimating the solar irradiance is proposed. Furthermore, the algorithm assumes that the solar irradiance changes linearly within a 1-h period and can be estimated based on the solar irradiance at the half clock and slope. The collected direct normal solar irradiance data of 2016 from eight solar radiation stations in China were used to demonstrate the proposed method and evaluate its performance by comparing the results with those from three conventional methods used in BEMPs along with the ground truth measurements. In addition, in this study, factors affecting the accuracy of the calculation results are explored. The results of the estimated direct incident solar radiation show that the proposed method achieves the best accuracy, followed by the methods used in DOE-2, EnergyPlus, and DeST. The proposed method guarantees that the hourly direct solar insolation will remain the same and reflects the variation in the direct solar irradiance across a 1-h time frame. The proposed method can be adopted in BEMPs to improve the accuracy of the solar radiation calculation, thereby improving the accuracy of the simulated building performance and the BIPV production.

14 SOLAR ENERGY↗

DOE OSTI2025

The Second Skin: A Wearable Sensor Suite that Enables Real-Time Human Biomechanics Tracking Through Deep Learning

Objective: Real-time determination of human kinematics and kinetics could advance biomechanics research and enable valuable applications of biofeedback and generalizable exoskeleton control. Here, this work aims to investigate a taskindependent, user-independent method for obtaining precise realtime joint state estimation across lower-body joints during a wide variety of tasks. Methods: We developed a generalizable sensing approach using a suit comprised of inertial measurement units (IMUs) and pressure insoles. With the suit, we collected a dataset of 33 tasks commonly performed during construction and hazardous waste cleanup (N = 10). We then trained deep learning user-independent, task-agnostic models to estimate joint lowerbody kinematics and dynamics using only worn sensor data. We likewise computed joint kinematics and dynamics analytically from sensor data to serve as a comparison tool for model results. Results: Our models achieved overall angle estimation root-meansquared-errors (RMSE) of 6.56±.92°, 8.60±1.01°, 7.58±.89°, and 6.00±.73° compared to 13.9±.1.3°, 15.31±1.0°, 10.76±.70°, and 7.56±.48° via analytical methods at the lower back, hip, knee, and ankle, respectively. Likewise, our models achieved overall normalized moment estimation RMSEs of .207±.069 Nm/kg, .242±.044 Nm/kg, .202±.038 Nm/kg, and .193±.034 Nm/kg compared to .306±.036 Nm/kg, .407±.021 Nm/kg, 1.18 ±.022 Nm/kg, and 1.73±.071 Nm/kg via analytical methods at the lower back, hip, knee, and ankle, respectively. Conclusion: These results are comparable to other state-of-the-art wearable sensing systems, establishing deep learning as a viable sensing approach that generalizes to new users and tasks. Significance: This work shows promise for enabling accurate real-world biomechanical data collection and enhancement of biofeedback systems and wearable robot control.

Casey, Ryan T. F. [Georgia Institute of Technology↗

NASA NTRS2010

Solar Array Disturbances to Spacecraft Pointing During the Lunar Reconnaissance Orbiter (LRO) Mission

The Lunar Reconnaissance Orbiter (LRO), the first spacecraft to support NASA s return to the Moon, launched on June 18, 2009 from the Cape Canaveral Air Force Station aboard an Atlas V launch vehicle. It was initially inserted into a direct trans-lunar trajectory to the Moon. After a five day transit to the Moon, LRO was inserted into the Lunar orbit and successfully lowered to a low altitude elliptical polar orbit for spacecraft commissioning. Successful commissioning was completed in October 2009 when LRO was placed in its near circular mission orbit with an approximate altitude of 50km. LRO will spend at least one year orbiting the Moon, collecting lunar environment science and mapping data, utilizing a suite of seven instruments to enable future human exploration. The objective is to provide key science data necessary to facilitate human return to the Moon as well as identification of opportunities for future science missions. LRO's instrument suite will provide the high resolution imaging data with sub-meter accuracy, highly accurate lunar cartographic maps, mineralogy mapping, amongst other science data of interest. LRO employs a 3-axis stabilized attitude control system (ACS) whose primary control mode, the "Observing Mode", provides Lunar nadir, off-nadir, and inertial fine pointing for the science data collection and instrument calibration. This controller combines the capability of fine pointing with on-demand large angle full-sky attitude reorientation. It provides simplicity of spacecraft operation as well as additional flexibility for science data collection. A conventional suite of ACS components is employed in the Observing Mode to meet the pointing and control objectives. Actuation is provided by a set of four reaction wheels developed in-house at NASA Goddard Space Flight Center (GSFC). Attitude feedback is provided by a six state Kalman filter which utilizes two SELEX Galileo Star Trackers for attitude updates, and a single Honeywell Miniature Inertial Measurement Unit (MIMU) to provide body rates for attitude propagation. Rate is computed by differentiating accumulated angle provided by the MIMU. The Observing Mode controller is required to maintain fine pointing while a large fully-articulated solar array (SA) maintains its panel normal to the solar incidence. This paper describes the disturbances to the attitude control resulting from the SA articulation. Observing Mode performance in the presence of this disturbance was assessed while the spacecraft was in an initial elliptical low altitude orbit during the commissioning phase, which started about two weeks after launch and lasted for 90 days. LRO demonstrated excellent pointing performance during Observing Mode nadir and inertial attitude target operations during this phase. Transient LRO attitude errors observed during commissioning resulted primarily from three sources, Diviner instrument calibrations, RW zero crossings, and SA articulation. Even during times of considerable disturbance from SA articulation, the attitude errors were maintained below the statistical attitude error requirement level of 15 arc-sec (3 sigma).

Calhoun, Philip↗

DOE OSTI2024

Sequential Kalman tuning of the t -preconditioned Crank-Nicolson algorithm: efficient, adaptive and gradient-free inference for Bayesian inverse problems

Ensemble Kalman Inversion (EKI) has been proposed as an efficient method for the approximate solution of Bayesian inverse problems with expensive forward models. However, when applied to the Bayesian inverse problem EKI is only exact in the regime of Gaussian target measures and linear forward models. Here, in this work we propose embedding EKI and Flow Annealed Kalman Inversion, its normalizing flow (NF) preconditioned variant, within a Bayesian annealing scheme as part of an adaptive implementation of the t-preconditioned Crank-Nicolson (tpCN) sampler. The tpCN sampler differs from standard pCN in that its proposal is reversible with respect to the multivariate t-distribution. The more flexible tail behaviour allows for better adaptation to sampling from non-Gaussian targets. Within our Sequential Kalman Tuning (SKT) adaptation scheme, EKI is used to initialize and precondition the tpCN sampler for each annealed target. The subsequent tpCN iterations ensure particles are correctly distributed according to each annealed target, avoiding the accumulation of errors that would otherwise impact EKI. We demonstrate the performance of SKT for tpCN on three challenging numerical benchmarks, showing significant improvements in the rate of convergence compared to adaptation within standard SMC with importance weighted resampling at each temperature level, and compared to similar adaptive implementations of standard pCN. The SKT scheme applied to tpCN offers an efficient, practical solution for solving the Bayesian inverse problem when gradients of the forward model are not available. Code implementing the SKT schemes for tpCN is available at https://github.com/RichardGrumitt/KalmanMC.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1960

A Homing Missile Control System to Reduce the Effects of Radome Diffraction

The problem of radome diffraction in radar-controlled homing missiles at high speeds and high altitudes is considered from the point of view of developing a control system configuration which will alleviate the deleterious effects of the diffraction. It is shown that radome diffraction is in essence a kinematic feedback of body angular velocities which causes the radar to sense large apparent line-of-sight angular velocities. The normal control system cannot distinguish between the erroneous and actual line-of-sight rates, and entirely wrong maneuvers are produced which result in large miss distances. The problem is resolved by adding to the control system a special-purpose computer which utilizes measured body angular velocity to extract from the radar output true line-of-sight information for use in steering the missile. The computer operates on the principle of sampling and storing the radar output at instants when the body angular velocity is low and using this stored information for maneuvering commands. In addition, when the angular velocity is not low the computer determines a radome diffraction compensation which is subtracted from the radar output to reduce the error in the sampled information. Analog simulation results for the proposed control system operating in a coplanar (vertical plane) attack indicate a potential decrease in miss distance to an order of magnitude below that for a conventional system. Effects of glint noise, random target maneuvers, initial heading errors, and missile maneuverability are considered in the investigation.

Smith, Gerald L.↗

NASA NTRS2011

Multi-Segment Radius Measurement Using an Absolute Distance Meter Through a Null Assembly

This system was one of the test methods considered for measuring the radius of curvature of one or more of the 18 segmented mirrors that form the 6.5 m diameter primary mirror (PM) of the James Webb Space Telescope (JWST). The assembled telescope will be tested at cryogenic temperatures in a 17-m diameter by 27-m high vacuum chamber at the Johnson Space Center. This system uses a Leica Absolute Distance Meter (ADM), at a wavelength of 780 nm, combined with beam-steering and beam-shaping optics to make a differential distance measurement between a ring mirror on the reflective null assembly and individual PM segments. The ADM is located inside the same Pressure-Tight Enclosure (PTE) that houses the test interferometer. The PTE maintains the ADM and interferometer at ambient temperature and pressure so that they are not directly exposed to the telescope s harsh cryogenic and vacuum environment. This system takes advantage of the existing achromatic objective and reflective null assembly used by the test interferometer to direct four ADM beamlets to four PM segments through an optical path that is coincident with the interferometer beam. A mask, positioned on a linear slide, contains an array of 1.25 mm diameter circular subapertures that map to each of the 18 PM segments as well as six positions around the ring mirror. A down-collimated 4 mm ADM beam simultaneously covers 4 adjacent PM segment beamlets and one ring mirror beamlet. The radius, or spacing, of all 18 segments can be measured with the addition of two orthogonally-oriented scanning pentaprisms used to steer the ADM beam to any one of six different sub-aperture configurations at the plane of the ring mirror. The interferometer beam, at a wavelength of 687 nm, and the ADM beamlets, at a wavelength of 780 nm, pass through the objective and null so that the rays are normally incident on the parabolic PM surface. After reflecting off the PM, both the ADM and interferometer beams return to their respective instruments on nearly the same path. A fifth beamlet, acting as a differential reference, reflects off a ring mirror attached to the objective and null and returns to the ADM. The spacings between the ring mirror, objective, and null are known through manufacturing tolerances as well as through an in situ null wavefront alignment of the interferometer test beam with a reflective hologram located near the caustic of the null. Since total path length between the ring mirror and PM segments is highly deterministic, any ADM-measured departures from the predicted path length can be attributed to either spacing error or radius error in the PM. It is estimated that the path length measurement between the ring mirror and a PM segment is accurate to better than 100 m. The unique features of this invention include the differential distance measuring capability and its integration into an existing cryogenic and vacuum compatible interferometric optical test.

Merle, Cormic↗

NASA NTRS2022

Thermal Hydraulic Disaggregation of SMAP Soil Moisture Over the Continental United States

Thermal Hydraulic disaggregation of Soil Moisture (THySM) algorithm was implemented to downscale NASA’s Soil Moisture Active Passive (SMAP) Enhanced soil moisture (SM) product to 1 km over the continental United States (CONUS). This algorithm was developed by combining thermal inertia theory with a soil hydraulic-based approach that considers fine-scale SM spatial distribution driven by both heat fluxes and hydraulic conductivity in soils. Relative soil wetness values were estimated using land surface temperature and normalized difference vegetation index for the thermal inertia model and using soil properties for the hydraulic model. The relative soil wetness values at 1 km from both models were then combined by using weighting functions whereby the spatial distribution of SM was governed more by thermal fluxes during times of strong heat transport and infiltration during moisture abundant soil conditions. THySM values were evaluated using in situ SM measurements from SMAP Core Validation Sites (CVS), the USDA Soil Climate Analysis Network, and the NOAA Climate Reference Network over CONUS. THySM shows higher accuracy than the SMAP / Sentinel-1 (SPL2SMAP_S) 1 km SM product when compared to in situ measurements. The accuracy of THySM is 0.048 m3/m3 based on unbiased root mean square error (ubRMSE), outperforming SPL2SMAP_S by 0.01-0.02 m3/m3. The ubRMSE of THySM 1km SM over the SMAP grassland/rangeland-dominated CVS sites is better than 0.04 m3/m3, which meets the SMAP mission SM accuracy requirement applied at 9 and 36 km.

SMAP↗

NASA NTRS1986

The abundance of interstellar sulphur and zinc in high density sight-lines

On the basis of early absorption line studies of individual lines of sight with the Copernicus satellite, chlorine, sulphur and zinc were classed together as elements which showed little or no depletion, relative to hydrogen, in the interstellar medium. The abundances of other less volatile elements, such as Fe and Mg were found to vary widely from one sight-line to another with gas-phase abundances in some cases being orders of magnitude below their solar counterparts. Detailed studies are reported of the depletion/density behavior of two other volatile elements which were previously considered to be virtually undepleted, S and Zn, using equivalent width data from both Copernicus and IUE observations. The results provide further evidence that the established dependence of depletion on n bar (H) extends to volatile elements and show that their use as tracers of metallicity, or for estimating hydrogen column densities, may lead to large errors in sight-lines through dense regions. It now appears that such elements may take part in the surface chemistry of grains and be important constituents of grain mantle material, although they probably do not contribute significantly to the bulk mass of grains. Due to the very similar atomic masses and ionization potentials of sulphur and phosphorous, the thermal velocity distributions of the singly ionized species of these elements in interstellar clouds should be very similar. However, a comparison of Doppler widths (b-values) derived for SIT and PIT in the same sight-lines from the Bohlin et al Copernicus equivalent width measurements has revealed an unexpected systematic discrepancy of a factor of approx. 1.7. This Discrepancy indicates that the normally adopted oscillators strengths of the PII lambda lambda 1153 and 1302 A lines may require revision.

Harris, A. W.↗

NASA NTRS2001

A Diagnostic Analysis of the Kennedy Space Center LDAR Network: Data Characteristics - 1

An analytic framework is developed in which to analyze climatological VHF (66 MHz) radiation measurements taken by the Kennedy Space Center Lightning Detection and Ranging (LDAR) network. A 19 month noise-filtered sample of LDAR observations is examined using this framework. It is found that the climatological impulsive VHF source density as observed by LDAR falls off approximately 10 dB every 71 km of ground range away from the network centroid (a 31 km e-folding scale). The underlying vertical distribution of impulsive VHF sources is approximately normally distributed with a mean altitude of 9 km and a standard deviation of 2.7 km; this implies that the loss of below-horizon sources has a negligible effect on column-integrated source densities within a 200 km ground range. At medium to far ranges, location errors are primarily radial and have a slightly asymmetric distribution whose standard deviation increases as r(exp 2). Error moments estimated from observed lightning are significantly higher than those from aircraft-based signal generator or analytic estimates. LDAR bulk flash detection efficiency is predicted to be above 90% to 94-113 km range from the network centroid and to fall below 10% at ranges greater than 200-240 km.

Boccippio, D. J.↗

NASA NTRS2011

Cryogenic Scan Mechanism for Fourier Transform Spectrometer

A compact and lightweight mechanism has been developed to accurately move a Fourier transform spectrometer (FTS) scan mirror (a cube corner) in a near-linear fashion with near constant speed at cryogenic temperatures. This innovation includes a slide mechanism to restrict motion to one dimension, an actuator to drive the motion, and a linear velocity transducer (LVT) to measure the speed. The cube corner mirror is double-passed in one arm of the FTS; double-passing is required to compensate for optical beam shear resulting from tilting of the moving cube corner. The slide, actuator, and LVT are off-the-shelf components that are capable of cryogenic vacuum operation. The actuator drives the slide for the required travel of 2.5 cm. The LVT measures translation speed. A proportional feedback loop compares the LVT voltage with the set voltage (speed) to derive an error signal to drive the actuator and achieve near constant speed. When the end of the scan is reached, a personal computer reverses the set voltage. The actuator and LVT have no moving parts in contact, and have magnetic properties consistent with cryogenic operation. The unlubricated slide restricts motion to linear travel, using crossed roller bearings consistent with 100-million- stroke operation. The mechanism tilts several arc seconds during transport of the FTS mirror, which would compromise optical fringe efficiency when using a flat mirror. Consequently, a cube corner mirror is used, which converts a tilt into a shear. The sheared beam strikes (at normal incidence) a flat mirror at the end of the FTS arm with the moving mechanism, thereby returning upon itself and compensating for the shear

Brasunas, John C.↗

NASA NTRS2019

Progress in Atmospheric Carbon Monitoring Using NASA's GEOS Model and Data from the OCO and GOSAT Missions

NASA's Global Modeling and Assimilation Office (GMAO) produces a variety of carbon products based the synthesis of satellite remote sensing data and outputs of the Goddard Earth Observing System (GEOS). This includes bottom-up surface fluxes due to fossil fuel emissions, biomass burning, terrestrial biospheric exchange, and ocean exchangeconstrained by measurements of nighttime lights, fire radiative power, normalized difference vegetation index, and ocean color. These fluxes are the basis of top-down estimates of carbon concentrations and fluxes. In particular, the GMAO system processes retrievals of column carbon dioxide (XCO2) from GOSAT and OCO-2 to produce a high-resolution, long-term global analysis of CO2 in three dimensions every 6 hours. Here, we discuss the potential applications of such products for satellite intercomparison and evaluation against independent, non-coincident data. We also highlight the ability to provide monthly global atmospheric growth rates inferred from the assimilated CO2 concentration product. Finally, we discuss the challenges facing such products including bias correction and the estimation and analysis of model transport errors.

Weir, B.↗

DOE OSTI2024

Superconducting microwave cavities and qubits for quantum information systems

Superconducting microwave cavities featuring ultrahigh Q-factors, which measure the efficiency of energy storage in relation to energy loss in a system, are revolutionizing quantum computing by providing long coherence times exceeding 1 ms, crucial for the development of scalable multi-qubit quantum systems with low error rates. Here, in this work, we provide an in-depth analysis of recent advances in ultrahigh Q-factor cavities, integration of Josephson junction-based qubits, and bosonic-encoded qubits in 3D cavities. We examine the sources of quantum state dephasing caused by damping and noise mechanisms in cavities and qubits, highlighting the critical challenges that need to be addressed to achieve even higher coherence times. We critically survey the latest progress made in implementing single 3D qubits using superconducting materials, normal metals, and multi-qubit and multi-state quantum systems. Our work sheds light on the promising future of this research area, including novel materials for cavities and qubits, modes with nontrivial topological properties, error correction techniques for bosonic qubits, and new light–matter interaction effects.

3D printing↗

NASA NTRS1994

Large-scale structure effects on the gravitational lens image positions and time delay

We compute the fluctuations in gravitational lens image positions and time delay caused by large-scale structure correlations. We show that these fluctuations can be expressed as a simple integral over the density power spectrum. Using the Cosmic Background Explorer (COBE) normalization we find that positions of objects at cosmological distances are expected to deviate from their true positions by few arcminutes. These deflections are not directly observable. The positions of the images relative to one another fluctuate by a few percent of the relative separation, implying that one does not expect multiple images to be produced by large-scale structure. Nevertheless, the fluctuations are larger than the observational errors on the positions and affect reconstructions of the lens potential. The time delay fluctuations have a geometrical and a gravitational contribution. Both are much larger than the expected time delay from the primary lens, but partially cancel each other. We find that large-scale structure weakly affects the time delay and time delay measurements can be used as a probe of the distance scale in the universe.

Seljak, Uros↗

DOE OSTI2023

SKS Polarization Anomalies Due to the Coriolis Force

The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

58 GEOSCIENCES↗

NASA NTRS2003

Stochastic Models of Tropical Rain-Rate Fields

Because of the extreme variability of rain rate in space and time and the difficulties with remote sensing methods of measuring rain rates, accurate determination of rainfall over large areas and time periods has long been a problem for hydrologists, meteorogists, and climatologists. A number of statistical models of rain have been developed in order to investigate the impact of rain variability on satellite remote sensing methods, validation of satellite rain products, and generation of rain maps with accompanying error estimates. These models may be useful in examining 'sub-grid scale' issues in representing precipitation in numerical mdoels. A stochastic model will first be described which can generate time-dependent high-resolution spatial rain fields with space and time correlations similar to those seen in rain data, as well as representing the presence of areas with zero rain rate and log-normally distributed rain rates where there is rain. A simpler model derived from this, formulated in the spectral domain, seems to imply fractal-like rain statistics at small scales when fit to rain data.

Bell, Thomas L.↗

NASA NTRS2004

GPS FOM Chimney Analysis using Generalized Extreme Value Distribution

Many a time an objective of a statistical analysis is to estimate a limit value like 3-sigma 95% confidence upper limit from a data sample. The generalized Extreme Value Distribution method can be profitably employed in many situations for such an estimate. . .. It is well known that according to the Central Limit theorem the mean value of a large data set is normally distributed irrespective of the distribution of the data from which the mean value is derived. In a somewhat similar fashion it is observed that many times the extreme value of a data set has a distribution that can be formulated with a Generalized Distribution. In space shuttle entry with 3-string GPS navigation the Figure Of Merit (FOM) value gives a measure of GPS navigated state accuracy. A GPS navigated state with FOM of 6 or higher is deemed unacceptable and is said to form a FOM 6 or higher chimney. A FOM chimney is a period of time during which the FOM value stays higher than 5. A longer period of FOM of value 6 or higher causes navigated state to accumulate more error for a lack of state update. For an acceptable landing it is imperative that the state error remains low and hence at low altitude during entry GPS data of FOM greater than 5 must not last more than 138 seconds. I To test the GPS performAnce many entry test cases were simulated at the Avionics Development Laboratory. Only high value FoM chimneys are consequential. The extreme value statistical technique is applied to analyze high value FOM chimneys. The Maximum likelihood method is used to determine parameters that characterize the GEV distribution, and then the limit value statistics are estimated.

Ott, Rick↗

DOE OSTI2026

Probabilistic Inference of Low-Surface-Brightness Galaxy Morphological Parameters Using Simulation-Based Inference

Low-surface-brightness galaxies (LSBGs) are diffuse, often dark-matter-dominated systems whose faintness makes their structural parameters difficult to measure reliably in wide-field imaging surveys. Robust parameter inference, including uncertainty quantification, is important for population studies and for comparisons with models of galaxy formation, as future surveys are expected to produce increasingly large samples of diffuse galaxies. In practice, LSBG profile modeling is sensitive to sky- background errors, masking choices, contaminating background sources, and the computational cost of obtaining posterior-level uncertainties for large samples. Motivated by these questions, we develop a simulation-based inference (SBI) framework for estimating posterior distributions of LSBG morphological parameters from simulated galaxy images. Using PyImfit, we generate DES-like single-Sersic profile LSBG images with known position angle, ellipticity, Sersic index, effective surface brightness, and effective radius. We then train a normalizing-flow-based neural posterior estimator using the sbi package to infer these parameters from the simulated images. For isolated simulated galaxies, the SBI posterior recovers the true input parameters, produces posterior predictive residuals consistent with the assumed noise model, and shows good empirical calibration in a DES-motivated test regime. We also compare SBI with PyImfit-based MCMC inference and find broadly comparable posterior constraints, while SBI enables substantially faster posterior sampling after training. Finally, we test robustness to compact background contaminants. A model trained only on isolated galaxies produces undercovered posteriors on contaminated images, whereas training on simulations with variable contaminant positions and fluxes improves calibration across contaminated test sets. These results demonstrate the promise of SBI for scalable, uncertainty-aware LSBG morphology inference, while emphasizing that posterior reliability strongly depends on whether training simulations include relevant observational complications.

Batbayar, Bilguun [U. Chicago (main)]↗