Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calibration and fitting methods”

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 55 records · Page 3

Characterization of a hard X-ray self-seeding diamond crystal orientation

We present a method to accurately control the photon energies for hard X-ray Self-seeding schemes with a single crystal monochromator in transmissive geometry. The energy calibration is performed by measuring which pairs of the machine pitch and yaw angles for different crystallographic planes reflect the X-ray at the same wavelength. The free parameters of an analytical formula for the self-seeding energies are determined by fitting the observed intersections and the normalized derivative with respect to the pitch and yaw angles in the observed intersections. The method requires a hard X-ray spectrometer, but it does not rely on its absolute energy calibration. Instead, identifying the self-seeded energies above the SASE background or the monochromatic notches within the SASE bandwidth is sufficient for the calibration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Constitutive model development of aluminum alloy 1100 for elevated temperature forming process

Commercially pure aluminum alloy, AA1100, presents good electrical and thermal conductivity, high formability, and low cost. Those favorable characteristics have the potential to enable bipolar plates with improved economics and enhanced performance compared to current stainless steel bipolar plates for proton exchange membrane fuel cells. An accurate constitutive model is essential to develop and optimize processing parameters and effectively control the forming process. Here, the objective of this work is to develop a constitutive model of AA1100 that is able to simulate stress-strain relation, formed geometry, and predict the onset of fracture strain to avoid forming failure. Initially, a set of tensile tests at temperature between 300 and 500°C and strain rate between 0.005 and 1.0/s were conducted to examine the deformation behavior. Then, a set of damage-based unified visco-plastic constitutive equations is proposed and calibrated based on the results of stress-strain data. A genetic algorithm optimization method is applied to search for best fitting material constants in constitutive equations. The proposed model shows good predictability of both the stress-strain relation and fracture strain at low strain rate and high temperature conditions. The accuracy of proposed model is also evaluated statistically. A comparison of the proposed model with three popular models (Arrhenius-type mode, Johnson-Cook model and Zerilli-Armstrong model) was made. The proposed model shows the best experimental agreement with correlation coefficient of 0.96 in contrast to 0.25, 0.38 and 0.75 for the popular models, respectively. The proposed model can help to optimize the elevated temperature forming process and guide die design to enable optimal geometric features in the formed components.

08 HYDROGEN↗

The DEHVILS survey overview and initial data release: high-quality near-infrared Type Ia supernova light curves at low redshift

ABSTRACT While the sample of optical Type Ia supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good ‘standard candles,’ is still ≲ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from 0.002 < z < 0.09 as a part of the Dark Energy, H0, and peculiar Velocities using Infrared Light from Supernovae (DEHVILS) survey. Observations are taken using UKIRT’s WFCAM, where the median depth of the images is 20.7, 20.1, and 19.3 mag (Vega) for Y, J, and H-bands, respectively. The median number of epochs per SN Ia is 18 for all three bands (Y,J, and H) combined and 6 for each band individually. We fit 47 SN Ia LCs that pass strict quality cuts using three LC models, SALT3, snoopy, and B ayeSN and find scatter on the Hubble diagram to be comparable to or better than scatter from optical-only fits in the literature. Fitting NIR-only LCs, we obtain standard deviations ranging from 0.128 to 0.135 mag. Additionally, we present a refined calibration method for transforming 2MASS magnitudes to WFCAM magnitudes using HST CALSPEC stars that results in a 0.03 mag shift in the WFCAM Y-band magnitudes.

Astronomy & Astrophysics↗

ZTF-SEDm Type Ia supernova sample for Twins Embedding spectrophotometric standardization

Aims. This paper has two aims: the first aim is to build a large homogeneous spectrophotometric sample of Type Ia supernovae (SNe Ia) from the second data release of the Zwicky Transient Facility (ZTF DR2). We used the spectrum sample from the low-resolution ( R ∼ 100) SEDmachine (SEDm) Integral Field Spectrograph (IFS) that gathered 3069 spectra. This is one of the largest samples of such collections that can attempt to reproduce the Twins Embedding (TE) spectrophotometric standardization method. This is our second objective. The method was developed based on high-quality spectra from 200 SNe Ia of the Nearby Supernova factory (SNfactory) and led to an exceptionally low value of 0.073 mag for the intrinsic scatter. Methods. As the SEDm is not designed as a spectrophotometric instrument, we first improved the flux-calibration accuracy of the SN Ia spectrum sample using the ZTF photometric data, which were calibrated at the percent level. We corrected the spectra for second-order polynomials, fitted by comparing the synthetic photometry in the ZTF g , r , i filters with the light-curve (LC) data. We then applied the three steps of the TE parameterization to a subset of 783 ZTF SN spectra near maximum light while comparing results from SNfactory and ZTF. We finally analyzed the standardization methods based on the TE parameters. Results. The precision of the phase-correction model, which is the first step of the TE, is estimated at 0.01 mag in g band based on ZTF data. Despite the challenge posed by the spectrum-extraction pipeline associated with the SEDm (flux calibration, leftover host signal, low signal-to-noise ratio, and low resolution), we applied a first standardization in color based on the second step of the TE, called read between the lines (RBTL), to the ZTF sample. We reached a Hubble residual scatter of 0.153 mag, all in normalized median absolute deviation, which is to be compared to the ∼0.11 mag obtained with the SNfactory data. The SALT color and stretch standardization reaches a scatter of 0.164 mag for the same ZTF SN Ia sample, and its host steps are ∼0.1 mag and zero for RBTL. When considering the scatter due to the redshift error and flux calibration error, we estimated a RBTL scatter of ∼0.129 mag for this ZTF sample as an upper limit because we identified an additional contribution from a systematic error in color. We tested the standardization based on the nonlinear TE parameters, and, as expected from the low spectrum quality, it did not improve the overall dispersion. Conclusions. We release 1897 flux calibrated spectra of 1607 SNe Ia with an estimated photometric accuracy of 0.07 mag. We further demonstrate that some amount of spectrophotometric SN Ia standardization can be done with limited-quality spectra. The RBTL standardization is more efficient than that of SALT with one parameter less, and the resulting host steps are consistent with zero. This makes it less prone to astrophysical bias. For future spectroscopic surveys, targeting the extraction pipeline for a thorough flux calibration and good signal-to-noise ratio would enable us to compute the full TE standardization, which would further reduce the scatter in the distance estimate.

Ganot, C↗

Star formation rate and stellar mass calibrations based on infrared photometry and their dependence on stellar population age and extinction

The stellar mass (M $\star$ ) and the star formation rate (SFR) are among the most important features that characterize galaxies. Measuring these fundamental properties accurately is critical for understanding the present state of galaxies, their history, and future evolution. Infrared (IR) photometry is widely used to measure the M $\star$ and SFR of galaxies because the near-IR traces the continuum emission of the majority of their stellar populations (SPs), and the mid/far-IR traces the dust emission powered by star-forming activity. This work explores the dependence of the IR emission of galaxies on their extinction, and the age of their SPs. It aims to provide accurate and precise IR-photometry SFR and M $\star$ calibrations that account for SP age and extinction while providing quantification of their scatter. We used the CIGALE spectral energy distribution (SED) fitting code to create model SEDs of galaxies with a wide range of star formation histories, dust content, and interstellar medium properties. We fit the relations between M $\star$ and SFR with IR and optical photometry of the model-galaxy SEDs with the Markov chain Monte Carlo (MCMC) method. As an independent confirmation of the MCMC fitting method, we performed a machine-learning random forest (RF) analysis on the same data set. The RF model yields similar results to the MCMC fits, thus validating the latter. This work provides calibrations for the SFR using a combination of the WISE bands 1 and 3, or the JWST NIR-F200W and MIRI-F2100W. It also provides mass-to-light ratio calibrations based on the WISE band-1, the JWST NIR-F200W, and the optical u - r or g - r colors. These calibrations account for the biases attributed to the SP age, while they are given in the form of extinction-dependent and extinction-independent relations. The proposed calibrations show robust estimations while minimizing the scatter and biases throughout a wide range of SFRs and stellar masses. The SFR calibration offers better results, especially in dust-free or passive galaxies where the contributions of old SPs or biases from the lack of dust are significant. Similarly, the M $\star$ calibration yields significantly better results for dusty and high-SFR galaxies where dust emission can otherwise bias the estimations.

79 ASTRONOMY AND ASTROPHYSICS↗

Increased accuracy and signal-to-noise ratio through recent improvements in infra-red video bolometer fabrication and calibration

The infra-red video bolometer (IRVB) is a diagnostic equipped with an infra-red camera that measures the total radiated power in thousands of lines of sight within a large field of view. Recently validated in MAST-U [Fderici et al., Rev. Sci. Instrum. 94, 033502 (2023)], it offers a high spatial resolution map of the radiated power in the divertor region, where large gradients are expected. The IRVB’s sensing element comprises a thin layer of high Z absorbing material, typically platinum, usually coated with carbon to reduce reflections [Peterson et al., Rev. Sci. Instrum. 79, 10E301 (2008)].Here, the possibility of using a relatively inert material such as titanium, is explored that can be produced in layers up to 1 μm compared to 2.5 μm for Pt and then coat it with Pt of the desired thickness (0.3 μm per side here) and carbon. This leads to a higher temperature signal (about 3 times) and better spatial resolution (about 4 times), resulting in higher accuracy in the measured power [Peterson et al., Rev. Sci. Instrum. 79, 10E301 (2008)]. This assembly is also expected to improve foil uniformity, as the Pt layer is obtained via deposition rather than mechanical processes [Mukai et al., Rev. Sci. Instrum. 87, 2014 (2016)].Given its multi-material composition, measuring the thermal properties of the foil assembly is vital. Various methods using a calibrated laser as a heat source have been developed, analyzing the temperature profile shape [Sano et al., Plasma and Fusion Res. 7, 2405039 (2012)] and [Mukai et al., Rev. Sci. Instrum. 89, 10E114 (2018)] or fitting the calculated laser power for different intensities and frequencies [Fderici et al., Rev. Sci. Instrum. 94, 033502 (2023)]. Here, a simpler approach is presented, which relies on analyzing the separate components of the foil heat equation for a single laser exposure in a given area. This can then be iterated over the entire foil to capture local deviations.

Magnetic Fusion Energy↗

Harmonizing solar induced fluorescence across spatial scales, instruments, and extraction methods using proximal and airborne remote sensing: A multi-scale study in a soybean field

Solar-induced chlorophyll fluorescence (SIF) has been widely used to track vegetation photosynthesis at different scales ranging from in-situ measurements to satellite products. Airborne platforms sample SIF data at a spatial scale intermediate between in-situ and satellite, matching that of ground measurement (e.g. flux tower footprints and other field sampling), enabling us to explore causes of SIF variation and validate satellite-based SIF products. However, harmonizing SIF across sensors and platforms (correcting for systematic errors to yield a consistent, comparable SIF product) is challenging because SIF can be retrieved in different absorption windows, with different instruments and methods complicating the comparison between different observational levels (i.e., ground, airborne, satellites) and between sites equipped with different instruments with varying optical properties (spectral resolution and sampling intervals, spatial resolution). Additionally, the spatial and temporal variability of atmospheric properties can influence the retrieval of the weak SIF signal. Because of these complications, direct comparisons of airborne and ground SIF across scales are rarely attempted. Here, in this study, we combined airborne SIF data with simultaneous ‘ground truth’ data collected by stationary and mobile platforms in a soybean field in Nebraska, USA. In this effort, we tested several SIF extraction methods, including Fraunhofer Line Discrimination (FLD), improved Fraunhofer Line Discrimination (iFLD), Spectral Fitting Method (SFM), SpecFit, and a Singular Vector Decomposition (SVD) method. The SpecFit method was sensitive to the 715–740 nm water bands and removing the water bands in the fitting process yielded better agreement between the airborne and ground SIF spectra. Accurate estimation of the ground level downwelling irradiance obtained by ground measurements over a calibration target improved agreement between airborne and ground SIF retrievals at the O 2 A band, and allowed us to derive a SIF dataset with improved agreement across platforms and sampling scales. This experimental approach provided a method for generating comparable SIF signals across instruments, methods and platforms, which is critical to understanding the SIF-GPP relationship at different scales and to cross-validate the diversity of platforms used for satellite products calibration and validation.

54 ENVIRONMENTAL SCIENCES↗

A proposed high-intensity radiometer calibration method using concentrated solar radiation

Accurate calibration of irradiance measurement devices, or radiometers, is essential for ensuring the reliability of measurements in high heat applications such as concentrating solar power (CSP), aerospace, defense, and pulsed power systems. Despite the critical need, existing calibration standards and service providers are limited to irradiance levels below 100 kW/m 2 and specific radiation sources, which is insufficient for many applications. For instance, CSP technologies, particularly those under the Department of Energy’s Solar Energy Technologies Office (SETO) Gen 3 program, require accurate measurements of broadband irradiance at levels exceeding 2000 kW/m 2 . In even more extreme scenarios, such as re-entry vehicles, heat levels can surpass 10000 kW/m 2 . Current ISO standards, specifically ISO 14934–2 and ISO 14934–3, are constrained to lower irradiance levels and dependent on black body heat sources, limiting their applicability for high-intensity broadband irradiance measurements, particularly in concentrated solar applications. Here, to address this shortfall, the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories (SNL) proposes a calibration method and facility capable of characterizing radiometers up to 2750 kW/m 2 using concentrated solar irradiance. Calibrating with concentrated sunlight is important for solar applications as it aligns the calibration process with the solar spectrum. This alignment is crucial for minimizing systematic errors and avoiding the need for additional corrections that may arise when radiometers designed for solar applications are calibrated using black-body or electrical sources. This paper presents the present day NSTTF characterization facility and procedure, detailing the proposed calibration method and uncertainty quantification. The presented method builds upon 1980′s NSTTF methodology and involves both theoretical and empirical methods to establish a robust relationship between gauge voltage output and irradiance intensity, quantifying both measurement and fitting errors. By addressing the limitations of existing standards and extending the characterization range, this work provides an advancement in the field of high-intensity irradiance measurement and instrumentation characterization.

Gardon gauge↗

SIF and Vegetation Indices in the US Midwestern Agroecosystems, 2016-2021, Carbon Monitoring System (CMS)

This dataset contains half-hourly ground solar-induced chlorophyll fluorescence (SIF) and vegetation indices including NDVI, EVI, Red edge chlorophyll index, green chlorophyll index, and photochemical reflectance index at seven crop sites in Nebraska and Illinois for the period 2016-2021. Four sites were located at Eddy Covariance (EC) tower sites (sites US-Ne2, US-Ne3, US-UiB, and US-UiC), and three sites were located on private farms (sites Reifsteck, Rund, and Reinhart). The sites were either miscanthus, corn-soybean rotation or corn-corn-soybean rotation. The spectral data for SIF retrieval and hyperspectral reflectance for vegetation index calculation were collected by the FluoSpec2 system, installed near planting, and uninstalled after harvest to collect whole growing-season data. Raw nadir SIF at 760 nm from different algorithms (sFLD, 3FLD, iFLD, SFM) are included. SFM_nonlinear and SFM_linear represent the Spectral fitting method (SFM) with the assumption that fluorescence and reflectance change with wavelength non-linearly and linearly, respectively. Additional data include two SIF correction factors including calibration coefficient adjustment factor (f_cal_corr_QEPRO) and upscaling nadir SIF to eddy covariance footprint factor (ratio_EC footprint, SIF pixel), and measured FPAR from quantum sensors and Rededge NDVI calculated FPAR. The data are provided in comma-separated values (CSV) format.

Ameriflux↗

Bump Morphology of the CMAGIC Diagram

Abstract We apply the color–magnitude intercept calibration method (CMAGIC) to the Nearby Supernova Factory SNe Ia spectrophotometric data set. The currently existing CMAGIC parameters are the slope and intercept of a straight line fit to the linear region in the color–magnitude diagram, which occurs over a span of approximately 30 days after maximum brightness. We define a new parameter, ω XY , the size of the “bump” feature near maximum brightness for arbitrary filters X and Y . We find a significant correlation between the slope of the linear region, β XY , in the CMAGIC diagram and ω XY . These results may be used to our advantage, as they are less affected by extinction than parameters defined as a function of time. Additionally, ω XY is computed independently of templates. We find that current empirical templates are successful at reproducing the features described in this work, particularly SALT3, which correctly exhibits the negative correlation between slope and “bump” size seen in our data. In 1D simulations, we show that the correlation between the size of the “bump” feature and β XY can be understood as a result of chemical mixing due to large-scale Rayleigh–Taylor instabilities.

79 ASTRONOMY AND ASTROPHYSICS↗

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗

Sources of error and methods to improve accuracy in interface state density analysis using quasi-static capacitance–voltage measurements in wide bandgap semiconductors

Characterizing interface trap states in commercial wide bandgap devices using frequency-based measurements requires unconventionally high probing frequencies to account for both fast and slow traps associated with wide bandgap materials. The C - ψs technique has been suggested as a viable quasi-static method for determining the interface trap state densities in wide bandgap systems, but the results are shown to be susceptible to errors in the analysis procedure. Here this work explores the primary sources of errors present in the C - ψs technique using an analytical model that describes the apparent response for wide bandgap MOS capacitor devices. Measurement noise is shown to greatly impact the linear fitting routine of the 1/C$^{*2}_s$ vs Ψ s plot to calibrate the additive constant in the surface potential/gate voltage relationship, and an inexact knowledge of the oxide capacitance is also shown to impede interface trap state analysis near the band edge. In addition, a slight nonlinearity that is typically present throughout the 1/C$^{*2}_s$ vs Ψ s plot hinders the accurate estimation of interface trap densities, which is demonstrated for a fabricated n-SiC MOS capacitor device. Methods are suggested to improve quasi-static analysis, including a novel method to determine an approximate integration constant without relying on a linear fitting routine.

42 ENGINEERING↗

Risk Ratio and Risk Difference Estimation in Case-cohort Studies

Background: In case-cohort studies with binary outcomes, ordinary logistic regression analyses have been widely used because of their computational simplicity. However, the resultant odds ratio estimates cannot be interpreted as relative risk measures unless the event rate is low. The risk ratio and risk difference are more favorable outcome measures that are directly interpreted as effect measures without the rare disease assumption. Methods: We provide pseudo-Poisson and pseudo-normal linear regression methods for estimating risk ratios and risk differences in analyses of case-cohort studies. These multivariate regression models are fitted by weighting the inverses of sampling probabilities. Also, the precisions of the risk ratio and risk difference estimators can be improved using auxiliary variable information, specifically by adapting the calibrated or estimated weights, which are readily measured on all samples from the whole cohort. Finally, we provide computational code in R (R Foundation for Statistical Computing, Vienna, Austria) that can easily perform these methods. Results: Through numerical analyses of artificially simulated data and the National Wilms Tumor Study data, accurate risk ratio and risk difference estimates were obtained using the pseudo-Poisson and pseudo-normal linear regression methods. Also, using the auxiliary variable information from the whole cohort, precisions of these estimators were markedly improved. Conclusion: The ordinary logistic regression analyses may provide uninterpretable effect measure estimates, and the risk ratio and risk difference estimation methods are effective alternative approaches for case-cohort studies. These methods are especially recommended under situations in which the event rate is not low.

60 APPLIED LIFE SCIENCES↗

Strong Lensing Cosmology with Population-level Calibrated Neural Ratio Estimation

Strong gravitational lensing contains key information about cosmic acceleration. Modern and next-generation galaxy imaging surveys are expected to provide high-quality data on $\mathcal{O}(10^5)$ galaxy-galaxy lensing systems. The plethora and complexity of the data are likely to present computational challenges for parameter inference methods for fitting high-dimensional likelihoods, which are often analytically intractable. Neural Ratio Estimation (NRE) efficiently computes individual likelihood ratios that can be combined into population-level posteriors. We use simulations to study the capacity of NRE to jointly predict the dark energy equation-of-state parameter $w$ and the total matter density $Ω_{m}$ from lensing images and companion spectroscopic information. We also introduce a post hoc posterior coverage calibration procedure that mitigates the model overconfidence that is typically found in neural density estimation applications. Our experiments show that the errors on both parameters decrease with increasing inference population sizes. In particular, for 100 lenses in a standard $Λ$CDM Universe, our calibrated NRE model achieves median fractional uncertainty of $22.8\%$ in $w$ and $2.9\%$ in $Ω_{m}$. This proof of concept demonstrates a potentially scalable approach for efficient cosmological parameter inference with large populations of galaxy-scale lenses observed in future surveys.

Jarugula, Sreevani [Fermilab] (ORCID:0000000253867↗

XCal: model-based approach to X-ray CT spectral calibration

Transmission X-ray computed tomography (CT) is widely used to quantitatively reconstruct 3D objects composed of multiple materials. However, accurate CT reconstruction requires the system to be calibrated to account for the effective X-ray spectrum. Unfortunately, measurement of the effective spectrum is ill-posed, and existing calibration methods require that the system be recalibrated when the system parameters are changed. In this paper, we propose XCal, a multi-energy model-based spectral calibration approach for X-ray CT. The XCal approach models the effective spectrum using a separable physics-based model of the CT system. The model parameters are then estimated by fitting calibration data with known objects at multiple energies. An important advantage of XCal is that it allows the user to change scanner settings, such as the source voltage or X-ray filters, without the need for recalibration. Evaluations on simulated and measured datasets demonstrate that XCal significantly improves the accuracy of the estimated spectrum as compared to existing calibration methods.

Li, Wenrui [Purdue Univ., West Lafayette, IN (Unit↗

Calibration of the light-flavour jet mistagging efficiency of the b -tagging algorithms with Z+jets events using 139 $\textrm{fb}^{-1}$ of ATLAS proton–proton collision data at $\sqrt{s} = 13$ TeV

The identification of b -jets, referred to as b -tagging, is an important part of many physics analyses in the ATLAS experiment at the Large Hadron Collider and an accurate calibration of its performance is essential for high-quality physics results. This publication describes the calibration of the light-flavour jet mistagging efficiency in a data sample of proton–proton collision events at $\sqrt{s} = 13$ TeV corresponding to an integrated luminosity of 139 fb -1 . The calibration is performed in a sample of Z bosons produced in association with jets. Due to the low mistagging efficiency for light-flavour jets, a method which uses modified versions of the b-tagging algorithms referred to as flip taggers is used in this work. A fit to the jet-flavour-sensitive secondary-vertex mass is performed to extract a scale factor from data, to correct the light-flavour jet mistagging efficiency in Monte Carlo simulations, while simultaneously correcting the b -jet efficiency. With this procedure, uncertainties coming from the modeling of jets from heavy-flavour hadrons are considerably lower than in previous calibrations of the mistagging scale factors, where they were dominant. The scale factors obtained in this calibration are consistent with unity within uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗