Search NASA⌕ Search

SEARCH · Search NASA

Results for “analysis of variance”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗
DOE OSTI2025

Boosting efficiency and reducing graph reliance: Basis adaptation integration in Bayesian multi-fidelity networks

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

42 ENGINEERING↗
DOE OSTI2023

Scrambling and quantum chaos indicators from long-time properties of operator distributions

Scrambling is a key concept in the analysis of nonequilibrium properties of quantum many-body systems. Most studies focus on its characterization via out-of-time-ordered correlation (OTOC) functions, particularly through the early-time decay of the OTOC. However, scrambling is a complex process which involves operator spreading and operator entanglement, and a full characterization requires one to access more refined information on the operator dynamics at several timescales. In this work we analyze operator scrambling by expanding the target operator in a complete basis and studying the structure of the expansion coefficients treated as a coarse-grained probability distribution in the space of operators. Here we study different features of this distribution, such as its mean, variance, and participation ratio, for the Ising model with longitudinal and transverse fields, kicked collective spin models, and random circuit models. We show that the long-time properties of the operator distribution display common features across these cases and discuss how these properties can be used as a proxy for the onset of quantum chaos. Finally, we discuss the connection with OTOCs and analyze the cost of probing the operator distribution experimentally using these correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantifying the impact of temporal analysis of products reactor initial state uncertainties on kinetic parameters

The temporal analysis of products (TAP) reactor, a transient kinetic tool, provides users with information as the catalyst state evolves. However, the state of the art for TAP uncertainty quantification only considers the experimental noise present in the outlet flow signal. Additional sources of uncertainty, including, but not limited to, initial surface coverages, the catalyst zone location, the inert void fraction, and the gas pulse intensity and pulse delay, are not included. For this reason, a framework for quantifying all uncertainty sources present in TAP experiments is presented and applied to a carbon monoxide oxidation case study. Herein, two methods for quantifying these sources of uncertainty are introduced. The first utilizes initial state sensitivities to approximate the parameter variances, as well as to provide insights into the structural certainty of the model. The second generates parameter confidence distributions through an ensemble-based sampling algorithm. This initial state covariance matrix can ultimately be merged with the experimental noise covariance matrix, providing a unified description of the parameter uncertainties for a single TAP experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2020

Data Collection and Analysis Challenges and Mitigation Strategies for Quantitative Human Factors Research Studies in Nuclear Power Plant Modernization

The United States (U.S.) Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) program Plant Modernization Pathway is conducting targeted research and development (R&D) to address aging and reliability concerns with the legacy instrumentation and control and related information systems of the U.S. LWR fleet. In this effort, the application of human factors engineering (HFE) provides an important role in ensuring new digital plant technologies enable broad innovation and business improvement with continued operational safety. Evaluation is a key activity in HFE, which often occurs iteratively through the system design lifecycle. While qualitative methods are important in collecting information of how users perform tasks through observations, quantitative methods are equally important in assessing system design based on performance. In collecting and analyzing this quantitative data, there are notable challenges in the nuclear HFE domain that may threaten the validity and reliability of the inferences made in these studies. Notable challenges include small sample size and limited resources, large error variance and small effect size, an “adding test features to losing degrees of freedom” dilemma, non-normal distribution, and heterogeneity of variance. The results in control room usability studies are often statistically non-significant, which makes it hard to interpret. This work discusses these challenges across different scientific viewpoints and provides real-world examples of these challenges in practice. Collectively, the objective of this work is to position these challenges to the larger data science community as a means of identifying future opportunities to address these issues.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Transitioning from Stage-III to Stage-IV: Cosmology from galaxy×CMB lensing and shear×CMB lensing (Final Report for Grant SC0021429 "Cosomology with Galaxy and CMB Surveys: Mitigating Systematic Effects Through Combining Datasets")

We examine the cosmological constraining power from two cross-correlation probes between galaxy and cosmic microwave background (CMB) surveys: the cross-correlation of lens galaxy density with CMB lensing convergence $\langle δ_g$$^κ\text{CMB}\rangle$, and source galaxy weak lensing shear with CMB lensing convergence $\langle^{γκ}$$_{\text{CMB}}\rangle$. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of $\langle δ_g$$^κ\text{CMB}\rangle + \langle^{γκ}$$_{\text{CMB}}\rangle$ changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from $\langle δ_g$$^κ\text{CMB}\rangle$ and $\langle^{γκ}$$_{\text{CMB}}\rangle$ are comparable in the Stage-III data sets; but in Stage-IV surveys, the noise in$\langle δ_g$$^κ\text{CMB}\rangle$ becomes subdominant to cosmic variance, preventing $\langle δ_g$$^κ\text{CMB}\rangle$ to further improve the constraints. This implies that to maximize the cosmological constraints from future $\langle δ_g$$^κ\text{CMB}\rangle$ + $\langle^{γκ}$$_{\text{CMB}}\rangle$ analyses, we should focus more on the requirements on$^{γκ}$$_{\text{CMB}}\rangle$ instead of$\langle δ_g$$^κ\text{CMB}\rangle$. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Flow boiling transient critical heat flux tests with stainless steel and FeCrAl: Transient correlation implementation, model calibration, and sensitivity analysis

In this study, steady-state and transient internal-flow critical heat flux (CHF) experiments were carried out under two flow conditions, at atmospheric pressure, on 316 stainless steel (316-SS) and iron-chromium-aluminum (FeCrAl) tubes. Slow power transients at a low mass flow (300 kg/m 2 - s) with very low subcooling (X e =-0.0054) generated premature CHFs, which were prevented by faster power transients and high mass flow (1000 kg/m 2 - s). The measured transient CHFs increased linearly with increasing power transients, compared with the steady-state CHF. Yet, the wall superheat at the CHFs decreased with faster power transients. Transient CHF correlations highlighting heterogeneous spontaneous nucleation were calibrated to the measured CHFs and compared with other existing correlations. Transient CHF multipliers were acquired from pool and flow boiling empirical CHF correlations that were generated. The multipliers were applied to the RELAP5-3D nuclear system code to analyze the discrepancies between the measured data and the predicted CHF and post-CHF behavior, which improved peak cladding temperature predictions by 24.6%. A variance-based global sensitivity study perturbing the experimental uncertainties and cladding material thermal properties showed the diminishing influence of flow boiling heat transfer with increasing power transients, highlighting the significance of the volumetric heat capacity for the cladding integrity during transients. Transient CHF correlations were applied to the most limiting design basis accident: a hot-zero-power reactivity-initiated accident on a generic pressurized water reactor RELAP5-3D model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Global Sensitivity Analysis of Bulk Properties of an Atomic Nucleus

We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus 16 O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose, we develop a subspace-projected coupled-cluster method using eigenvector continuation [Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)]. With this method, we compute the binding energy and charge radius of 16 O at more than 106 different values of the 16 low-energy constants in one hour on a standard laptop computer. For relatively small subspace projections, the root-mean-square error is about 1% compared to full-space coupled-cluster results. We find that 58(1)% of the variance in energy can be apportioned to a single contact term in the 3 S 1 wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segrè chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Plasmon-Induced Hot-Carrier Excited-State Dynamics in Plasmonic Semiconductor Nanocrystals

The variance of carrier relaxation pathways for WO 3–x plasmonic semiconductor nanocrystals (PSNCs) is monitored by transient absorption spectroscopy following excitation of the localized surface plasmon resonance (LSPR) versus the optical band gap (E g,opt ). Excitation of the LSPR leads to efficient hot carrier population above the Fermi level in WO 3–x via Landau damping, in analogy to noble metal LSPR relaxation mechanisms. Hot carrier depopulation occurs on the femtosecond timescale, observed as the concomitant recovery of an LSPR bleach with the appearance of discrete interband and intraband photoinduced absorption features. Further, by comparison, the direct excitation of E g,opt results in trion recombination at donor–acceptor sites within the WO 3–x NC, consistent with exciton decay dynamics observed for typical wide-band-gap semiconductor NCs. From the analysis of pump power dependency data, a hot-carrier electron–phonon coupling constant of 1.47 × 10 11 J K –1 s –1 cm –3 is extracted. The direct comparison of the decay dynamics following E g,opt versus LSPR excitation confirms that the observed plasmon in trioctylphosphine oxide passivated, spherical WO 3–x is a resonance state in which hot carriers are generated only from excitation on resonance with the LSPR frequency. This study on WO 3–x PSNCs provides a toolset that can be used to evaluate the role of hot carriers following LSPR excitation of n-type, plasmonic transparent conducting oxide NCs, where enhancement of photocatalysis, photovoltaic performance, and optical enhancement has been reported.

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

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

LiteBIRD is a planned JAXA-led cosmic microwave background (CMB) B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, δr, down to δr < 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio r uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on r after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher χ 2 sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on r increases by only about 5% for a 500 GHz configuration compared to the baseline.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Examining the changes in the spatial manifestation and the rate of arrival of large tornado outbreaks

This study presents an assessment of the spatial and temporal characteristics of large tornado outbreak (LTOs) days, in which several counties were impacted by tornadoes rated F2(EF2) or greater on the Fujita (Enhanced Fujita) scale in one day. A statistical evaluation of changes in the LTO clusters for two periods, 1950–1980 and 1989–2019, has been performed. There is a geographical shift of the nucleus (central impact location) towards the southeast United States. This spatial shift is also accompanied by reduced spatial variance, suggesting LTOs have become less dispersed (or more localized) in the recent period. The overall inter-arrival rate of LTOs, and how it changed during successive 31-year climatological blocks between 1950–2019 was investigated using an exponential probability model. The arrival rate has changed from 124 days during 1950–1980 to 164 days during 1977–2007 and remained relatively constant during later periods, indicating that LTOs are becoming less frequent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Passive heat removal in horizontally oriented micro-HTGRs

There are novel, horizontally oriented high-temperature gas-cooled micro-reactor (HTGR) designs under consideration because of their portability, simplicity, and reliable operation. These features of the micro-HTGR design can meet the fission battery attributes, such as economic, standardize, installed, unattended and reliable. Although there have been HTGR related safety studies, micro-HTGRs when horizontally oriented are expected to exhibit significantly different thermal physics. During the unavailability of helium circulation, the internal reactor core is designed to cool by block-to-block conduction and radiation, and the reactor vessel surface is cooled by the ambient air. This scenario is anticipated during the transport of the micro-HTGR in a shipping container. The conduction and radiation between the prismatic micro-HTGR blocks in the core can be influenced by variances in the thermal contacts. This work investigated the conduction within a simulated horizontal HTGR core. An experimental setup was used to validate a numerical model. The experimental setup consisted of a hexagonal assembly with scaled prismatic blocks placed within a high-temperature vacuum environment. The gaps between the blocks were well controlled and monitored. The experimental setup was designed in such a way that the temperature variation in the axial direction was minimal, such that the experiment could be observed as a 2D (r, ) heat transfer problem. The pressurized tube was equipped with Infra-Red (IR) transparent windows, and an IR camera was used to capture the spatio-temporal temperature evolution in the hexagonal block assembly during the cooldown experiments. The experimental scenario was computationally modeled with a finite element analysis (FEA) program. Once validated, the computational model was used to investigate the impact of gap conductance on overall decay heat removal. Using a conservative estimate for gap conductance value (100 W/m 2 -K) between the prismatic blocks, there is a negligible increase in temperature observed during decay heat generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Computational study of the adsorption of bimetallic clusters on alumina substrate

We performed computational investigation of the adsorption of bimetallic Pd 3 M 2 (where M changes from Ag, Au, Co, Cu, Mn, Ni, Pt, and Ru) cluster on the hydroxylated alumina surface. Previously, it was shown that small silver cluster can control the rate of discharge at the cathode of lithium-oxygen battery. The gap near the fermi energy was shown to control the oxygen reduction, an important reaction for LiO 2 formation. Controlling the gap would ultimately control the rate of LiO 2 formation. One can vary the size of the cluster to vary the gap, however, this "knob" provides limited variance in the gap. Alloying, in combination with size variation offer a much wider control of the gap, hence the LiO 2 formation. Using Density Functional Theory (DFT), we determined the most stable geometry of the bimetallic clusters Pd 3 M 2 and calculated the binding energies of these clusters on the alumina substrate which ranges from 0.2 eV to 0.25 eV depending upon the composition of the alloy-cluster, its orientation and the adsorption site. We also find that Pd atoms bind strongly with the substrate oxygen atoms with an average short bond-length of about 2.2 angstrom. We explored how the gap at the Fermi level of the system varies as a function of elemental composition and the calculated gap ranges from 0 meV to 90 meV. Charges distribution using Bader analysis was also performed to probe how charges are transferred between the cluster and the substrate. These preliminary results will open the door for more systematic studies of alloy clusters of different size and stoichiometry for Li-O 2 battery cathode design.

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

Agent-Based Model of Combined Community- and Jail-Based Take-Home Naloxone Distribution

Importance Opioid-related overdose accounts for almost 80 000 deaths annually across the US. People who use drugs leaving jails are at particularly high risk for opioid-related overdose and may benefit from take-home naloxone (THN) distribution. Objective To estimate the population impact of THN distribution at jail release to reverse opioid-related overdose among people with opioid use disorders. Design, Setting, and Participants This study developed the agent-based Justice-Community Circulation Model (JCCM) to model a synthetic population of individuals with and without a history of opioid use. Epidemiological data from 2014 to 2020 for Cook County, Illinois, were used to identify parameters pertinent to the synthetic population. Twenty-seven experimental scenarios were examined to capture diverse strategies of THN distribution and use. Sensitivity analysis was performed to identify critical mediating and moderating variables associated with population impact and a proxy metric for cost-effectiveness (ie, the direct costs of THN kits distributed per death averted). Data were analyzed between February 2022 and March 2024. Intervention Modeled interventions included 3 THN distribution channels: community facilities and practitioners; jail, at release; and social network or peers of persons released from jail. Main Outcomes and Measures The primary outcome was the percentage of opioid-related overdose deaths averted with THN in the modeled population relative to a baseline scenario with no intervention. Results Take-home naloxone distribution at jail release had the highest median (IQR) percentage of averted deaths at 11.70% (6.57%-15.75%). The probability of bystander presence at an opioid overdose showed the greatest proportional contribution (27.15%) to the variance in deaths averted in persons released from jail. The estimated costs of distributed THN kits were less than $\$$15 000 per averted death in all 27 scenarios. Conclusions and Relevance This study found that THN distribution at jail release is an economical and feasible approach to substantially reducing opioid-related overdose mortality. Training and preparation of proficient and willing bystanders are central factors in reaching the full potential of this intervention.

Tatara, Eric [Argonne National Laboratory (ANL), A↗
DOE OSTI2023

Large-scale analysis of structural brain asymmetries in schizophrenia via the ENIGMA consortium

Left–right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case–control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case–control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case–control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case–control status. Subtle case–control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Rare Events via Cross-Entropy Population Monte Carlo

Rare events are events that happen with very low frequency. Estimating rare event probabilities using Monte Carlo techniques is computationally expensive, often to the point of intractability, and special methods are required. Importance sampling (IS) is a well known technique that uses a proposal distribution in place of a target distribution to lower the variance of the estimator. Key to the success of IS methods is the choice of a proposal distribution, or the parameters governing the distribution. Adaptive importance sampling improves the parameters of a family or population of proposal distributions iteratively through trials. We present a novel cross-entropy population Monte Carlo algorithm, which adapts the parameters of proposals through the cross-entropy method. The proposed method stands apart from previous work in that we are not optimizing a mixture distribution. Instead, we leverage deterministic mixture weights and optimize the distributions individually through a reinterpretation of the typical derivation of the cross-entropy method. Demonstrations on rare event examples show that the algorithm can outperform existing resampling based population Monte Carlo methods, especially for higher-dimensional problems. Finally, we also demonstrate efficacy on a conjunction analysis problem.

97 MATHEMATICS AND COMPUTING↗