Quasiparticle band structure and excitonic optical response in V 2 O 5 bulk and monolayer
Not Available
Garcia, Claudio (ORCID:0009000077503536)↗
SEARCH · Engineering Papers
Results for “Localized basis sets”
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.
The DECADE cosmic shear project I: A new weak lensing shape catalog of 107 million galaxies
We present the Dark Energy Camera All Data Everywhere (DECADE) weak lensing dataset: a catalog of 107 million galaxies observed by the Dark Energy Camera (DECam) in the northern Galactic cap. This catalog was assembled from public DECam data including survey and standard observing programs. These data were consistently processed with the Dark Energy Survey Data Management pipeline as part of the DECADE campaign and serve as the basis of the DECam Local Volume Exploration survey (DELVE) Early Data Release 3 (EDR3). We apply the Metacalibration measurement algorithm to generate and calibrate galaxy shapes. After cuts, the resulting cosmology-ready galaxy shape catalog covers a region of $5,\!412 \,\,{\rm deg}^2$ with an effective number density of $4.59\,\, {\rm arcmin}^{-2}$. The coadd images used to derive this data have a median limiting magnitude of $r = 23.6$, $i = 23.2$, and $z = 22.6$, estimated at ${\rm S/N} = 10$ in a 2 arcsecond aperture. We present a suite of detailed studies to characterize the catalog, measure any residual systematic biases, and verify that the catalog is suitable for cosmology analyses. In parallel, we build an image simulation pipeline to characterize the remaining multiplicative shear bias in this catalog, which we measure to be $m = (-2.454 \pm 0.124) \times10^{-2}$ for the full sample. Despite the significantly inhomogeneous nature of the data set, due to it being an amalgamation of various observing programs, we find the resulting catalog has sufficient quality to yield competitive cosmological constraints.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems
Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation
This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems
The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Greedy emulators for nuclear two-body scattering
Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.
Bayesian methods↗
Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations
Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf–O or Hf–B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet–triplet energy gap (ΔE S–T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange–correlation functionals spanning multiple rungs of “Jacob’s ladder.” When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Guide to Energy Master Planning of High-Performance Districts and Communities
This guide was developed with partners throughout the United States to demonstrate how implementing district-scale high-performance strategies can be successful and scalable approaches to achieving deep energy savings that increase affordability, improve resilience, reduce emissions, and foster economic development. This document serves as a framework for districts, campuses, and communities, illustrating an iterative process of building support for, planning, and implementing high-performance districts by engaging stakeholders, setting aggressive energy goals, completing technical and financial planning, and implementing a high-performance energy master plan. The information in this guide is based on a 3-year U.S. Department of Energy Zero Energy District Accelerator and a range of real-world examples of emerging high-performance districts. It is particularly useful for architects, planners, engineers, local government agencies, and real estate developers in the early phases of planning a district with high-performance or other deep energy goals. For the purposes of this guide, a high-performance district is a multibuilding development that achieves aggressive energy and related goals such as zero energy, carbon neutrality, sustainability, ultra-efficiency, etc. High-performance districts optimize energy efficiency to reduce energy loads and use renewable energy resources to meet the remaining loads whenever possible. The Energy Independence and Security Act of 2007 defines a high-performance building as “a building that integrates and optimizes on a life cycle basis all major high-performance attributes, including energy conservation, environment, safety, security, durability, accessibility, cost-benefit, productivity, sustainability, functionality, and operational considerations.” High-performance districts are collections of such buildings that take advantage of the synergies available when energy consumption and production are considered at a district level rather than one building at a time.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Multireference Methods for Chemistry and Materials Science: Automated Active Spaces, Efficient Dynamic Correlation, and Extended Systems
While multiconfigurational approaches have long been relegated to expert practitioners working on a case-by-case basis, recent developments have increasingly made these methods more routine and applicable to broader sets of systems. This article outlines the state-of-the-art in multiconfigurational approaches, with an emphasis on moving from delicate hand-selected pathways through configuration space toward more robust and efficient approaches to treating a host of challenging chemical systems accurately. First, we overview recent work in automated active-space selection, which has enabled increasingly large-scale applications of multireference methods to modeling vertical excitations and reactivity. Second, we highlight the increasingly efficient methods for recovering correlation energy beyond the active space, as headlined by extensions of pair-density functional theory and its role in accurate and efficient treatment of excited-state dynamics and its utilization to train machine-learned potentials. Finally, we highlight recent efforts to treat extended systems that until recently have lied beyond the traditional limits of active-space methods, giving center stage to product-form wave functions of the localized active space family of methods that allow for the computation of multiconfigurational band structures. These recent advancements point to a broader use of multireference approaches for high-impact chemical and materials science applications.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Progress in pedestal and edge physics: Chapter 3 of the special issue: on the path to tokamak burning plasma operation
This paper describes the extensive progress that has been made in the understanding of tokamak pedestal physics since the 2007 publication of ‘Progress in the ITER Physics Basis’ (Ikeda 2007 Nucl. Fusion 47 E01–S500). It serves as Chapter 3 of the 2025 Nuclear Fusion Special Issue titled ‘On the Path to Tokamak Burning Plasma Operation’ (Campbell et al 2025 Nucl. Fusion ). This review was compiled by the pedestal and edge physics (PEP) community affiliated with the International Tokamak Physics Activity organization. It attempts to collect in one place citations to the majority of published literature on the pedestal physics topics that will be most important for the operation of a future power producing burning plasma tokamak. These include citations to publications describing the physics of the pedestal plasmas in many operating tokamaks worldwide and the pedestal physics projections for several near-term future devices including ITER. Descriptions of experimental results, interpretive modeling and predictive extrapolations are integrated together and comprehensive references are provided. This review is organized around four primary technical sections, viz.: pedestal structure, edge localized mode (ELM) characteristics, ELM control and regimes without large ELMs. Key results from many of the references are described briefly and set into the tokamak burning plasma power plant context. In addition, different perspectives on pedestal physics topics that are currently under debate within the community are also described, to provide guidance on needs for future research. Finally, attempts are made to describe conclusions from all of this progress consistent with discussions by the pedestal physics community at this time. The goal of this review is to provide a useful reference document for pedestal physics researchers going forward toward operation of a burning tokamak fusion plasma.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Strategic Energy Plan: Playa de Ponce, Ponce, Puerto Rico
Un Nuevo Amanecer Inc. (UNA), a community-based NGO (501C(3)), sought technical assistance through the U.S. Department of Energy's (DOE's) Energy Technology Innovation Partnership Project (ETIPP) to develop a strategic energy plan (SEP) for the community of Playa de Ponce on the south coast of Puerto Rico. A SEP helps communities identify energy goals and priority projects, and outline strategies and actions to build a shared path toward their energy vision. The primary goal of this SEP is to collaboratively explore pathways for the community's energy future. Researchers from the National Laboratory of the Rockies (NLR) engaged with the community for a period of seven months and facilitated a series of in-person community engagement workshops along with ETIPP regional partner Puerto Rico Hispanic Federation, UNA, and community leaders. The SEP workshops were held twice, March 18-20 and June 24-26, 2025. Through facilitated discussions, the community worked together to define and articulate a clear vision for its short-, medium-, and long-term energy goals. The workshops played a key role in identifying the specific strategies and priority projects that would best suit the needs of the community. The discussions outlined the community's energy aspirations, establishing a shared vision, setting goals, and outlining the priority actions needed to bring them to fruition. During the workshops, the community shared some of the most pressing energy challenges that they face on a regular basis, such as high energy costs that disproportionately impact an aging population and unreliable electricity during extreme weather events that result in frequent power outages. These issues have also exacerbated risks to public health and safety of the community. In response, UNA requested technical assistance through ETIPP to inform the development of a SEP tailored to Playa de Ponce to support its energy vision. This document explores four community-driven energy strategies aligned with that vision. Figure 1 shows the strategies outlined in this plan, which are categorized into three focus areas as identified by the community: energy reliability; safety and security; and stakeholder participation, education, and capacity building. This document serves as a shared roadmap for developing and implementing energy solutions that meet the needs of Playa de Ponce. It is intended to serve as a public resource for the community, UNA, and local authorities - providing tools, information and a shared foundation to inform future implementation efforts. This plan is intended to complement ongoing planning frameworks at the municipal, regional, and island-wide levels. As a living document, this SEP can evolve alongside Playa de Ponce's changing energy goals, while offering a framework to help them realize their energy vision.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
General, Rigorous Approach for the Treatment of Interfragment Covalent Bonds
Here, a generalized, projection-based transformation of the method-agnostic Fock operator in various ab initio fragment-based quantum chemistry methods has been developed for the treatment of interfragment covalent bonds. This transformation freezes the relevant localized molecular orbital associated with each interfragment bond, thereby restricting the variational subspace of the fragment wave functions, in order to maintain the proper physical characteristics of the involved covalent bonds. In addition, sets of orbitals that would lead to multiple occupancy of certain orbitals are explicitly removed from the variational space. The transformation is developed for the specific case of mutually orthonormal frozen and unfrozen orbitals within each fragment. The newly developed approach is then used to study model systems with two popular ab initio fragment-based methods, and the results of these calculations are compared to those obtained by existing methodologies. Analysis is focused on both quantitative and qualitative accuracy as well as computational scalability and stability. Other methods for which the developed formalisms are appropriate are outlined, and future extensions of the methods are discussed.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint
In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.
36 MATERIALS SCIENCE↗
Analytical gradients for molecular-orbital-based machine learning
We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Low-Energy Isomers of the Magic Number H + (H 2 O) 21 Cluster
Electronic structure calculations are used to characterize low-energy isomers of H + (H 2 O) 21 . Eleven different classes of isomers, based on the (H 2 O) 20 pentagonal dodecahedron with the excess proton localized on the surface (as a hydrated hydronium ion) and the “extra” water molecule located in the interior of the cluster, are characterized. In 10 of these classes, the internal water molecule is engaged in six 5-membered rings, but in the remaining class, which is predicted to start at only 0.6 kcal/mol above the global minimum, the internal water is engaged in a 4- membered ring, an additional 6-membered ring, and four 5- membered rings. In addition, isomers with two 4-membered rings and two 6-membered rings on the cluster surface are predicted to start at only ∼1.3 kcal/mol above the lowest-energy dodecahedralbased structure.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Simulating organic aerosol in Delhi with WRF-Chem using the volatility-basis-set approach: exploring model uncertainty with a Gaussian process emulator
The nature and origin of organic aerosol in the atmosphere remain unclear. The gas–particle partitioning of semi-volatile organic compounds (SVOCs) that constitute primary organic aerosols (POAs) and the multigenerational chemical aging of SVOCs are particularly poorly understood. The volatility basis set (VBS) approach, implemented in air quality models such as WRF-Chem (Weather Research and Forecasting model with Chemistry), can be a useful tool to describe emissions of POA and its chemical evolution. However, the evaluation of model uncertainty and the optimal model parameterization may be expensive to probe using only WRF-Chem simulations. Gaussian process emulators, trained on simulations from relatively few WRF-Chem simulations, are capable of reproducing model results and estimating the sources of model uncertainty within a defined range of model parameters. In this study, a WRF-Chem VBS parameterization is proposed; we then generate a perturbed parameter ensemble of 111 model runs, perturbing 10 parameters of the WRF-Chem model relating to organic aerosol emissions and the VBS oxidation reactions. This allowed us to cover the model's uncertainty space and to compare outputs from each run to aerosol mass spectrometer observations of organic aerosol concentrations and O:C ratios measured in New Delhi, India. The simulations spanned the organic aerosol concentrations measured with the aerosol mass spectrometer (AMS). However, they also highlighted potential structural errors in the model that may be related to unsuitable diurnal cycles in the emissions and/or failure to adequately represent the dynamics of the planetary boundary layer. While the structural errors prevented us from clearly identifying an optimized VBS approach in WRF-Chem, we were able to apply the emulator in the following two periods: the full period (1–29 May) and a subperiod period of 14:00–16:00 h LT (local time) on 1–29 May. The combination of emulator analysis and model evaluation metrics allowed us to identify plausible parameter combinations for the analyzed periods. We demonstrate that the methodology presented in this study can be used to determine the model uncertainty and to identify the appropriate parameter combination for the VBS approach and hence to provide valuable information to improve our understanding of OA production.
54 ENVIRONMENTAL SCIENCES↗
Opacity Table Generation
Radiation opacity is a key property of material systems which describes the absorption and scattering of photons. Such processes are an important part of energy transport, diagnostics, and astrophysical phenomena. Recent experiments in iron opacity and continuum lowering have sparked renewed interest in and careful vetting of opacity models. Even within the local thermodynamic equilibrium approximation, opacity varies over a large parameter space: material, temperature, density, and energy. Opacity tables that span wide regimes of this parameter space are of prime importance to radiation-hydrodynamics codes. Accurate opacity calculations cannot be achieved inline during a radiation hydrodynamics simulation since they involve accounting for all possible electron transitions of all important atomic states of a thermal ensemble. Livermore’s current framework utilizes precalculated tables at many densities, temperatures, and energies. Last year, the opacity theory team completed an L2 milestone concerning the development of a new code, Opus, which served the dual purpose of creating a modern code infrastructure and version control as well as the basis for training new members of that team, a critical need as key members near or are passed retirement. For the current L2 milestone, we exercised and validated Opus on three elements: boron, carbon, and nitrogen, which have few enough electrons to allow a careful convergence study over many input parameters. Here, we report on the production of new tables, show their convergence with respect to important parameters, state some of their limitations, and document a mostly automated production process. In this section, we address the milestone completion criteria, found on the title page and the subsection headings. In section 2, we validate Opus against Tycho and compare with TabOp and Atomic. Section 3 details some of the physics implemented in Opus and relevant to this report. We then cover details of the density-temperature grid in section 4 and do a spectral comparison between different parameter settings and different codes in section 5. After summarizing these results in section 6, we give instructions how to run the opacity table scripts in appendix A and document the scripts and input/output files in appendix B.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Updates to Relevance Vector Machine: Multiclass Classification, Variable Selection, and Proof-of-Concept Application to Safeguards Fresh Fuel Verification using List-Mode Neutron Collar Data
To expand the capabilities of safeguards authorities to verify the integrity of fresh fuel assemblies, Oak Ridge National Laboratory has retrofit the existing electronics of the JCC-71 uranium neutron coincidence collar, which contains 18 3 He neutron detectors and an external 241 AmLi(α, n) neutron interrogation source arranged to surround a fresh nuclear fuel assembly. The new electronics system allows analysts to record list-mode neutron multiplicity data in addition to the singles and doubles rates that are currently measured. Based on previous proof-of-concept research, analysis of these new data will identify off-normal fuel configurations in an assembly and characterize or localize the specific partial fuel defects. The purpose of this report it to document the analysis algorithm development and then to demonstrate its capability for the safeguards verification of fresh fuel assemblies using list mode neutron collar data. To analyze the complex list-mode data collected with the upgraded uranium neutron collar, multivariate classification algorithms are being developed using a novel classification method, the relevance vector machine. This approach may be applied to multiclass problems to estimate the probability that test data belongs to one of many possible classes of data. In addition, our method identifies the most useful variables/channels for making predictions, which illuminates the basis for the model’s predictions, and this interpretability is largely unique among data analytics methods. Variable selection occurs during model training and parameter tuning and does not need any external hyperparameter tuning routines. Finally, we apply the modified relevance vector machine to a simulated dataset of list-mode neutron collar data generated with the radiation transport code MCNP. The method can correctly identify off-normal fuel configurations, categorize the data according to four fuel defect scenarios, and rank the channels in the data according to prediction utility. For nuclear safeguards applications, it is concluded that this method has the potential to increase the sensitivity and reliability to detect missing fuel rods from a standard 17 x 17 Pressurized Water Reactor (PWR) fresh fuel assembly. Within this analysis, “off-normal” (i.e., missing fuel rods) were correctly classified in 17 simulated test scenarios with one quarter (25%) of the fresh fuel rods missing using a training data set of 58 simulated measurements.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗