Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model assignment”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
DOE OSTI2023

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part II - Clustering Extremely High-Dimensional Grid-Based Data

Building an accurate surrogate model for the spatio-temporal outputs of a computer simulation is a challenging task. A simple approach to improve the accuracy of the surrogate is to cluster the outputs based on similarity and build a separate surrogate model for each cluster. This clustering is relatively straightforward when the output at each time step is of moderate size. However, when the spatial domain is represented by a large number of grid points, numbering in the millions, the clustering of the data becomes more challenging. In this report, we consider output data from simulations of a jet interacting with high explosives. These data are available on spatial domains of different sizes, at grid points that vary in their spatial coordinates, and in a format that distributes the output across multiple files at each time step of the simulation. We first describe how we bring these data into a consistent format prior to clustering. Borrowing the idea of random projections from data mining, we reduce the dimension of our data by a factor of thousand, making it possible to use the iterative k-means method for clustering. We show how we can use the randomness of both the random projections, and the choice of initial centroids in k-means clustering, to determine the number of clusters in our data set. Our approach makes clustering of extremely high dimensional data tractable, generating meaningful cluster assignments for our problem, despite the approximation introduced in the random projections.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Initial Mobility Analysis for ORNL VA-EDH Synthetic Populations

Travel burdens are a major barrier to healthcare access among US Veteran patient populations, particularly those residing in rural areas. Spatial accessibility to points of care for US Veteran populations is commonly assessed in two ways. The first approach uses open data from the US Census to represent collective travel burdens, for example the distance between population-weighted census tract centroids and VHA points of care. The second approach uses restricted-access VHA patient data to measure travel costs (e.g., distance, time) for accessing points of care with respect to geolocated patient addresses and real or approximated transportation networks. While the advantage of the open data approach lies in its reproducibility, it has notable limitations in its tendency to infer individual travel behavior from aggregate population characteristics, a problem known as ecological fallacy. Conversely, while the patient data approach is able to account for individual travel behavior, its ability to account for localized access disparities (e.g., a neighborhood with exceptionally high transportation costs) and patient demographics is limited as protecting individual patient data requires their storage in closed systems with limited capacity for adequately modeling real-world travel patterns or for supplementing patient attributes. Additionally, the patient data approach cannot account for veterans who are not enrolled in the VHA system but who may be eligible for care. These challenges limit the ability to perform “what if” analyses on the effects of place-specific interventions on veteran populations with high access barriers to healthcare. To address these challenges, we explore the application of realistic synthetic populations to examine travel burdens and spatial accessibility issues among veteran patient populations. Synthetic populations provide a virtual, individually-resolved and cross-sectional representation of the veteran patient population that enables investigation of spatial access to points of care in ways in which aggregate data and patient data do not. First, synthetic populations allow one to directly assess how individuals access points of care, from synthesized residential locations to outpatient facilities on real-world transportation networks. Modeling access to points of care at the individual scale addresses the ecological fallacy problem associated with using aggregated census data to represent veteran populations and patterns of movement. Second, synthetic populations provide a means of completely representing an area’s veteran population using only publicly available, anonymized census microdata from the American Community Survey (ACS) to ensure the privacy of real-world individuals. Generating synthetic populations from the ACS also expands descriptive characteristics beyond what patient data typically offers to include socio-demographic, economic, housing, and mobility attributes. More detailed profiles of both VHA patient populations and veterans not enrolled in the VA system will provide a comprehensive picture of groups that may benefit from interventions or outreach. As an initial exercise for using synthetic populations to measure veteran travel burdens to VA care, we apply Oak Ridge National Laboratory’s (ORNL) UrbanPop capability to generate a series of synthetic VHA patient populations for 9 Veterans Integrated Services Networks (VISN) market areas in 9 Census Divisions across the continental United States, which are listed in Table 1. We use UrbanPop to produce synthetic populations for the VISN markets selected for each US Census Division, then assign VA outpatient clinic destinations to synthetic VHA patients based on travel about each VISN market’s road network. To demonstrate using the synthetic populations to evaluate healthcare travel burdens, we compare the time-based impedance between simulated home locations and VA outpatient clinics in each VISN market. We then perform validation exercises on the synthetic populations with respect to neighborhood (block group) demographic composition as well as patient mobility, comparing aggregate origin-destination statistics for the synthetic population to outpatient visits available in restricted patient data from the VA’s Corporate Data Warehouse (CDW) database.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Numerical Optimization of Spray-Guided Spark Assistance for Cold Idle Operation in a Heavy-Duty Gasoline Compression Ignition Engine

This article describes the results of a response surface model (RSM)-based numerical optimization campaign for spray-guided spark assistance at cold operations in a heavy-duty gasoline compression ignition (GCI) engine. On the basis of an earlier work on spark-assisted GCI cold combustion, a space-filling design of experiments (DoE) method was first undertaken to investigate a multitude of hardware design variables and engine operating parameters. The main design variables included the number of injector nozzles, fuel split quantities and injection timings, and spark timing. The objective variables were engine combustion efficiency (ŋ c ), maximum pressure rise rate (MPRR), and engine-out nitrogen oxide (NO x ) emissions. A total of 150 design candidates were automatically generated using the Sobol sequence method provided by the commercial software package, CAESES. Then, closed-cycle computational fluid dynamic (CFD) spark-assisted GCI simulations under cold idling operations were performed. The outcomes from the CFD-DoE design campaign were utilized to construct high-fidelity RSMs that allowed for further design optimization of the spark plug- and fuel injector-related design variables, along with fuel injection strategy parameters. A merit function with respect to objective variables was formulated with an appropriate weight assignment on each objective variable. Finally, the best design candidate was identified from the RSM-based optimization process and further validated in the CFD analysis. The best design candidate showed the potential to significantly improve combustion efficiency (ŋ c > 90%) over the baseline at cold idle while satisfying MPRR and NOx emissions constraints (MPRR < 5 bar/CAD and NO x < 4.5 g/kWh).

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

A representation-independent electronic charge density database for crystalline materials

Abstract In addition to being the core quantity in density-functional theory, the charge density can be used in many tertiary analyses in materials sciences from bonding to assigning charge to specific atoms. The charge density is data-rich since it contains information about all the electrons in the system. With the increasing prevalence of machine-learning tools in materials sciences, a data-rich object like the charge density can be utilized in a wide range of applications. The database presented here provides a modern and user-friendly interface for a large and continuously updated collection of charge densities as part of the Materials Project. In addition to the charge density data, we provide the theory and code for changing the representation of the charge density which should enable more advanced machine-learning studies for the broader community.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spectral evidence of solar neighborhood analogs in CALIFA galaxies

We introduce a novel nonparametric method to find solar neighborhood analogs (SNAs) in extragalactic integral field spectroscopic surveys. The main ansatz is that the physical properties of the solar neighborhood (SN) should be encoded in its optical stellar spectrum. We assume that our best estimate of such a spectrum is the one extracted from the analysis performed by the Code for Stellar properties Heuristic Assignment (CoSHA) from the MaStar stellar library. It follows that finding SNAs in other galaxies consist in matching, in a χ 2 sense, the SN reference spectrum across the optical extent of the observed galaxies. We applied this procedure to a selection of CALIFA galaxies, by requiring a close to face-on projection, relative isolation, and non-active galactic nucleus. We explore how the local and global properties of the SNAs (stellar age, metallicity, dust extinction, mass-to-light ratio, stellar surface mass density, star-formation density, and galactocentric distance) and their corresponding host galaxies (morphological type, total stellar mass, star-formation rate, and effective radius) compared with those of the SN and the Milky Way (MW). We find that SNAs are located preferentially in S(B)a–S(B)c galaxies, in a ring-like structure, which radii seem to scale with the galaxy size. Despite the known sources of systematics and errors, most properties present a considerable agreement with the literature on the SN. We conclude that the solar neighborhood is relatively common in our sample of SNAs. Our results warrant a systematic exploration of correlations among the physical properties of the SNAs and their host galaxies. We reckon that our method should inform current models of the galactic habitable zone in our MW and other galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Reductive Dynamic and Static Excited State Quenching of a Homoleptic Ruthenium Complex Bearing Aldehyde Groups

A new homoleptic Ru polypyridyl complex bearing two aldehyde groups on each bipyridine ligand, [Ru(dab) 3 ](PF 6 ) 2 , where dab is 4,4′-dicarbaldehyde-2,2′-bipyridine, was synthesized, characterized, and utilized for iodide photo-oxidation studies. In acetonitrile (CH 3 CN) solution, the complex displayed an intense metal-to-ligand charge transfer (MLCT) absorbance maximum at 475 nm (ε = 22,000 M –1 cm –1 ) and an infrared (IR) band at 1712 cm –1 assigned to the pendent aldehyde groups. Visible light excitation in air-saturated solution resulted in room temperature photoluminescence (PL) with a maximum at 675 nm, a quantum yield, ϕ PL = 0.048, and an excited state lifetime, τ ο = 440 ns, from which radiative and nonradiative relaxation rate constants were extracted, k r = 9.1 × 10 4 s –1 and knr = 1.8 × 10 6 s –1 . Pulsed visible light excitation yielded transient UV–vis and IR absorption spectra consistent with an MLCT excited state; relaxation occurred with the maintenance of two isosbestic points in the visible region, and a lifetime that agreed with that measured by time-resolved PL. Cyclic voltammetry studies in a CH 3 CN solution with 0.1 M TBAPF 6 electrolyte revealed a quasi-reversible oxidation, E°(Ru III/II ) = +1.25 V vs. Fc +/0 , and three sequential one-electron reductions at −1.10, −1.25, and −1.54 V vs. Fc +/0 . Here, an excited state reduction potential of E°(Ru *2+/+ ) = +0.89 V vs. Fc +/0 was estimated with the Rehm–Weller expression. Titration of tetrabutylammonium iodide, TBAI, into a CD 3 CN solution of [Ru(dab) 3 ](PF 6 ) 2 resulted in significant shifts in the aldehyde H atom and 3,3′-biypridyl resonances that were analyzed with a 1:1 equilibrium model, from which K eq = 460 M –1 was extracted, increasing to 5800 M –1 when the solvent was changed to acetone-d 6 . Iodide titrations resulted in a significant quenching of the [Ru(dab) 3 ] *2+ lifetime and quantum yield in both CH 3 CN and acetone solvents. In CH 3 CN, the quenching was mainly dynamic and well described by the Stern–Volmer model, from which a quenching rate constant, k q , of 4.5 × 10 10 M –1 s –1 and an equilibrium constant, K eq , of 8.3 × 10 3 M –1 were obtained. In acetone, the static quenching pathway by iodide was greatly enhanced, with a K eq of 1.2 × 10 4 M –1 and a higher k q of 9.2 × 10 10 M –1 s –1 .

Aldehydes↗
DOE OSTI2025

Wavelength dependent photochemistry of an iron dinitrogen hydride complex via multiple spectroscopies – competing ejection of axial ligands

Nitrogenase (N 2 ase) is a critical enzyme which catalyzes the reaction of N 2 → NH 3 in nature. Studies on the spectroscopy and photochemistry of trans-[Fe II (DMeOPrPE) 2 (N 2 )H][BPh 4 ] (1) and its isotopologues (2–6) provide a possible first step to evaluate the geometries and properties of the real N 2 ase–N 2 structure(s). In this article, we have used FT-IR, FT-Raman, synchrotron-based nuclear resonant vibrational spectroscopy (NRVS) and DFT calculations to examine and assign the normal modes of these complexes. In addition, we have monitored their wavelength dependent photochemistry using mid-IR, near-IR, NRVS, and Mössbauer spectroscopies. Two distinct photolysis pathways are observed with mid-IR at (nominal) 4 K – (1) the cleavage of Fe–N 2 bond in UV or visible light photolyses, which presents a unipolar disappearance of the N 2 peak at 2094 cm −1 and is recombinable; (2) the ejection of trans hydrogen atom with UV irradiation, which has a pair of bipolar peaks with the disappearance of N 2 at 2094 cm −1 and the appearance of a new species at 2056 cm −1 and is non-recombinable. The latter peak is well aligned with the N 2 peak in an Fe I reference complex (7). The combination of mid IR monitored photolysis/recombination and NRVS monitored photolysis form the central evidence for the conclusions in this article. In particular, the Fe II –N 2 and Fe I –H·dissociations are in competition with each other in UV or UV-inclusive photolyses of this dinitrogen hydride complex. In addition, near-IR and Mössbauer also provide consistent evidence about Fe I . This wavelength dependent photochemical work is the first one on a reaction active N 2 ase–N 2 model complex and it also demonstrates the competition ejection between two axial ligands (H· and N 2 ). It offers valuable information for future studies on real N 2 ase–N 2 and its photolysis products.

Dong, Weibing [Liaoning Normal University (China)]↗
DOE OSTI2024

Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2 ) n , n = 1,2

Pyrimidine has two in-plane CH(δ+)/N̈(δ–)/CH(δ+) binding sites that are complementary to the (δ–/2δ+/δ–) quadrupole moment of CO 2 . For this study, we recorded broadband microwave spectra over the 7.5–17.5 GHz range for pyrimidine-(CO 2 ) n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14 N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO 2 complex is planar, with CO 2 occupying one of the quadrupolar binding sites, forming a structure in which the CO 2 is stabilized in the plane by interactions with the C–H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO 2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191–195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C 2v symmetry in which the second CO 2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO 2 is –13.7 kJ mol –1 , including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO 2 ) 2 is almost exactly double that size, indicating little interaction between the two CO 2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

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

A decision support tool for e-waste recycling operations using the hen-and-chicks bio-inspired optimization metaheuristic

E-waste from end-of-life electrical and electronic devices is one of the fastest growing waste streams from households and businesses. E-waste recycling yields environmental sustainability and economic benefits. Due to continuous changes in e-waste types and compositions, recycling businesses face challenges to optimize their operational configuration to achieve better economic and environmental performance. To help e-waste recyclers mitigate this problem, we have developed a modular decision support tool called the Comprehensive Manufacturing Assessment Tool (CMAT) that can simulate both e-waste recycling operations and economics. This tool can give valuable insights regarding the profitability of the entire operation and different e-waste types. In addition, a new bio-inspired metaheuristic optimization algorithm, hen-and-chicks optimization (HACO), was developed to assign manpower to different workstations to maximize operational efficiencies. According to the results of our case study, laptops, desktops, and computer peripherals are the three electronic waste products that produce the most profit. Our examination of the sensitivity of material prices shows that the price of steel has the most significant influence on total profit, because it is the most widely used material in the majority of electronic devices. We have released the decision support tool as open-source software under a general public license. It could be customized for other recycling industries beyond e-waste to achieve business sustainability by making their operations more efficient.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

High-resolution infrared spectroscopy of supersonically cooled singlet carbenes: Bromomethylene (HCBr) in the CH stretch region

First high-resolution spectra of cold (~35 K) singlet bromomethylene HCBr in the CH stretching (v 1 ) region from 2770 to 2850 cm-1 are reported using near quantum shot-noise limited laser absorption methods in a slit jet supersonic discharge expansion source. Three rovibrational bands are identified at high S/N (20:1 - 40:1) and rotationally assigned to i) the CH stretch fundamental (v 1 ) band $\tilde{X}$(1,0,0) ← $\tilde{X}$ (0,0,0) and ii) vibrational hot bands ($\tilde{X}$(1,1,0) ← $\tilde{X}$(0,1,0) and $\tilde{X}$(1,0,1) ← $\tilde{X}$(0,0,1)) arising from vibrationally excited HCBr populated in the discharge with single quanta in either the H–C–Br bend (v 2 ) or C–Br stretch (v 3 ) modes. Precision rotational constants are reported for a total of six states, with an experimentally determined CH stretch vibrational frequency (2799.38 cm -1 ) in good agreement with previous low-resolution fluorescence studies. Detailed analysis of the fundamental v 1 band highlights the presence of perturbations in the $\tilde{X}$(1,0,0) level, which we tentatively attribute to arise from the nearby triplet state $\tilde{a}$(0,0,1) through spin-orbit interaction or the multiple quanta $\tilde{X}$(0,2,1) singlet state via c-type Coriolis coupling. Reduced-Doppler resolution (60 MHz) in the slit-jet IR spectrometer permits clear observation of nuclear spin hyperfine structure, with experimental line shapes well reproduced by nuclear quadrupole/spin-rotation coupling constants from microwave studies. Lastly, the a-type to b-type transition intensity ratio for the fundamental CH stretch band is notably larger than predicted by a bond-dipole model, which from high level ab initio quantum calculations (CCSD(T)/PVQZ) can be attributed to vibrationally induced “charge-sloshing” of electron density along the polar C–Br bond.

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

Comparative Performance of Three Eye-Tracking Devices in Detection of Mild Traumatic Brain Injury in Acute Versus Chronic Subject Populations

ABSTRACT Introduction Presently, traumatic brain injury (TBI) triage in field settings relies on symptom-based screening tools such as the updated Military Acute Concussion Evaluation. Objective eye-tracking may provide an alternative means of neurotrauma screening due to sensitivity to neurotrauma brain-health changes. Previously, the US Army Medical Research and Development Command Non-Invasive NeuroAssessment Devices (NINAD) Integrated Product Team identified 3 commercially available eye-tracking devices (SyncThink EYE-SYNC, Oculogica EyeBOX, NeuroKinetics IPAS) as meeting criteria toward being operationally effective in the detection of TBI in service members. We compared these devices to assess their relative performance in the classification of mild traumatic brain injury (mTBI) subjects versus normal healthy controls. Materials and Methods Participants 18 to 45 years of age were assigned to Acute mTBI, Chronic mTBI, or Control group per study criteria. Each completed a TBI assessment protocol with all 3 devices counterbalanced across participants. Acute mTBI participants were tested within 72 hours following injury whereas time since last injury for the Chronic mTBI group ranged from months to years. Discriminant analysis was undertaken to determine device classification performance in separating TBI subjects from controls. Area Under the Curves (AUCs) were calculated and used to compare the accuracy of device performance. Device-related factors including data quality, the need to repeat tests, and technical issues experienced were aggregated for reporting. Results A total of 63 participants were recruited as Acute mTBI subjects, 34 as Chronic mTBI subjects, and 119 participants without history of TBI as controls. To maximize outcomes, poorer quality data were excluded from analysis using specific criteria where possible. Final analysis utilized 49 (43 male/6 female, mean [x̅] age = 24.3 years, SD [s] = 5.1) Acute mTBI subjects, and 34 (33 male/1 female, x̅ age = 38.8 years, s = 3.9) Chronic mTBI subjects were age- and gender-matched as closely as possible with Control subjects. AUCs obtained with 80% of total dataset ranged from 0.690 to 0.950 for the Acute Group and from 0.753 to 0.811 for the Chronic mTBI group. Validation with the remaining 20% of dataset produced AUCs ranging from 0.600 to 0.750 for Acute mTBI group and 0.490 to 0.571 for the Chronic mTBI group. Conclusions Potential eye-tracking detection of mTBI, per training model outcomes, ranged from acceptable to excellent for the Acute mTBI group; however, it was less consistent for the Chronic mTBI group. The self-imposed targeted performance (AUC of 0.850) appears achievable, but further device improvements and research are necessary. Discriminant analysis models differed for the Acute versus Chronic mTBI groups, suggesting performance differences in eye-tracking. Although eye-tracking demonstrated sensitivity in the Chronic group, a more rigorous and/or longitudinal study design is required to evaluate this observation. mTBI injuries were not controlled for this study, potentially reducing eye-tracking assessment sensitivity. Overall, these findings indicate that while eye-tracking remains a viable means of mTBI screening, device-specific variability in data quality, length of testing, and ease of use must be addressed to achieve NINAD objectives and DoD implementation.

General & Internal Medicine↗
DOE OSTI2020

The European wood pellets for heating market - Price developments, trade and market efficiency

Competitive international markets imply adjustments towards competitive spatial equilibrium in which excess from one market is transferred to another and prices are equilibrated except for remaining differences that can be assigned to transfer costs. The European market for wood pellets used in small-scale heating systems has been expanding significantly over the past decade. Small scale pellet heating is arguably a mature technology, but whether the market is mature is another question. In this paper we analyse recent data on trade flows and price developments between Italy, Austria, Germany and France to understand the developments of wood pellet market efficiency and to draw conclusions about market function. The objective of this study is to establish a framework to test the European residential wood pellet market for competitive spatial equilibrium using modern trade theory. We find mainly inefficiently integrated markets with remaining positive marginal profits and detectable arbitrageurs’ activity. Based on a thorough discussion of these findings and the underlying data we outline possible methodology advancements and list policy recommendations to secure access and affordability of this renewable heating commodity in the long run.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

HarDWR - Harmonized Water Rights Records

A dataset within the Harmonized Database of Western U.S. Water Rights (HarDWR). For a detailed description of the database, please see the meta-record v2.0. Changelog v2.0 - Recalculated based on data sourced from WestDAAT - Changed using a Site ID column to identify unique records to using aa combination of Site ID and Allocation ID - Removed the Water Management Area (WMA) column from the harmonized records. The replacement is a separate file which stores the relationship between allocations and WMAs. This allows for allocations to contribute to water right amounts to multiple WMAs during the subsequent cumulative process. - Added a column describing a water rights legal status - Added "Unspecified" was a water source category - Added an acre-foot (AF) column - Added a column for the classification of the right's owner v1.02 - Added a .RData file to the dataset as a convenience for anyone exploring our code. This is an internal file, and the one referenced in analysis scripts as the data objects are already in R data objects. v1.01 - Updated the names of each file with an ID number less than 3 digits to include leading 0s v1.0 - Initial public release Description Here we present an updated database of Western U.S. water right records. This database provides consistent unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of seven broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of inter-sectoral water allocation changes though water markets (Grogan et al., *in review*). Specifically, the data were formatted for use as input to a process-based hydrologic model, Water Balance Model (WBM), with a water rights module (Grogan et al., *in review*). While this specific study motivated the development of the database presented here, water management in the U.S. West is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al., 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. West. We produced the water rights database presented here in four main steps: (1) data collection, (2) data quality control, (3) data harmonization, and (4) generation of cumulative water rights curves. Each of steps (1)-(3) had to be completed in order to produce (4), the final product that was used in the modeling exercise in Grogan et al. (*in review*). All data in each step is associated with a spatial unit called a Water Management Area (WMA), which is the unit of water right administration utilized by the state in which the right came from. Steps (2) and (3) required use to make assumptions and interpretation, and to remove records from the raw data collection. We describe each of these assumptions and interpretations below so that other researchers can choose to implement alternative assumptions an interpretation as fits their research aims. Motivation for Changing Data Sources The most significant change has been a switch from collecting the raw water rights directly from each state to using the water rights records presented in WestDAAT, a product of the Water Data Exchange (WaDE) Program under the Western States Water Council (WSWC). One of the main reasons for this is that each state of interest is a member of the WSWC, meaning that WaDE is partially funded by these states, as well as many universities. As WestDAAT is also a database with consistent categorization, it has allowed us to spend less time on data collection and quality control and more time on answering research questions. This has included records from water right sources we had previously not known about when creating v1.0 of this database. The only major downside to utilizing the WestDAAT records as our raw data is that further updates are tied to when WestDAAT is updated, as some states update their public water right records daily. However, as our focus is on cumulative water amounts at the regional scale, it is unlikely most records updates would have a significant effect on our results. The structure of WestDAAT led to several important changes to how HarWR is formatted. The most significant change is that WaDE has calculated a field known as `SiteUUID`, which is a unique identifier for the Point of Diversion (POD), or where the water is drawn from. This separate from `AllocationNativeID`, which is the identifier for the allocation of water, or the amount of water associated with the water right. It should be noted that it is possible for a single site to have multiple allocations associated with it and for an allocation to be able to be extracted from multiple sites. The site-allocation structure has allowed us to adapt a more consistent, and hopefully more realistic, approach in organizing the water right records than we had with HarDWR v1.0. This was incredibly helpful as the raw data from many states had multiple water uses within a single field within a single row of their raw data, and it was not always clear if the first water use was the most important, or simply first alphabetically. WestDAAT has already addressed this data quality issue. Furthermore, with v1.0, when there were multiple records with the same water right ID, we selected the largest volume or flow amount and disregarded the rest. As WestDAAT was already a common structure for disparate data formats, we were better able to identify sites with multiple allocations and, perhaps more importantly, allocations with multiple sites. This is particularly helpful when an allocation has sites which cross WMA boundaries, instead of just assigning the full water amount to a single WMA we are now able to divide the amount of water between the number of relevant WMAs. As it is now possible to identify allocations with water used in multiple WMAs, it is no longer practical to store this information within a single column. Instead the stAllocationToWMATab.csv file was created, which is an allocation by WMA matrix containing the percent Place of Use area overlap with each WMA. We then use this percentage to divide the allocation's flow amount between the given WMAs during the cumulation process to hopefully provide more realistic totals of water use in each area. However, not every state provides areas of water use, so like HarDWR v1.0, a hierarchical decision tree was used to assign each allocation to a WMA. First, if a WMA could be identified based on the allocation ID, then that WMA was used; typically, when available, this applied to the entire state and no further steps were needed. Second was the spatial analysis of Place of Use to WMAs. Third was a spatial analysis of the POD locations to WMAs, with the assumption that allocation's POD is within the WMA it should belong to; if an allocation still had multiple WMAs based on its POD locations, then the allocation's flow amount would be divided equally between all WMAs. The fourth, and final, process was to include water allocations which spatially fell outside of the state WMA boundaries. This could be due to several reasons, such as coordinate errors / imprecision in the POD location, imprecision in the WMA boundaries, or rights attached with features, such as a reservoir, which crosses state boundaries. To include these records, we decided for any POD which was within one kilometer of the state's edge would be assigned to the nearest WMA. Other Changes WestDAAT has Allowed In addition to a more nuanced and consistent method of assigning water right's data to WMAs, there are other benefits gained from using the WestDAAT dataset. Among those is a consistent categorization of a water right's legal status. In HarDWR v1.0, legal status was effectively ignored, which led to many valid concerns about the quality of the database related to the amounts of water the rights allowed to be claimed. The main issue was that rights with legal status' such as "application withdrawn", "non-active", or "cancelled" were included within HarDWR v1.0. These, and other water rights status' which were deemed to not be in use have been removed from this version of the database. Another major change has been the addition of the "unspecified water source category. This is water that can come from either surface water or groundwater, or the source of which is unknown. The addition of this source category brings the total number of categories to three. Due to reviewer feedback, we decided to add the acre-foot (AF) column so that the data may be more applicable to a wider audience. We added the ownerClassification column so that the data may be more applicable to a wider audience. File Descriptions The dataset is a series of various files organized by state sub-directories. In addition, each file begins with the state's name, in case the file is separate from its sub-directory for some reason. After the state name is the text which describes the contents of the file. Here is each file described in detail. Note that st is a placeholder for the state's name. stFullRecords_HarmonizedRights.csv: A file of the complete water records for each state. The column headers for each of this type of file are: state - The name of the state to which the allocations belong to. FIPS - The two digit numeric state ID code. siteID - The site location ID for POD locations. A site may have multiple allocations, which are the actual amount of water which can be drawn. In a simplified hypothetical, a farm stead may have an allocation for "irrigation" and an allocation for "domestic" water use, but the water is drawn from the same pumping equipment. It should be noted that many of the site ID appear to have been added by WaDE, and therefore may not be recognized by a given state's water rights database. allocationID - The allocation ID for the water right. For most states this is the water right ID, and what is recommended to use should a right be looked up on a given state's water rights database. The water amounts associated with these IDs tend to be finer scaled than those associated with siteID. It should be noted that some allocations may be extracted from multiple sites, particularly for larger Places of Use. ownerClassification - A classification of the types of owners for water rights. The most common is `Private` which incorporates a wide range of entities. Several classifications would be grouped into a government category, most of which are for the U.S. Federal Government. These allocations could be listed as "Federal", "United States of America", or as the names of any number of federal agencies. The last major grouping of entities is for "Native American"s. priorityDate - The date we use as the water right priority date for our modeling analysis. This is the legal priority date when it is available. However, for some rights, specifically from California and New Mexico, we used a pseudo priority date (e.g. well completion date or start of well drilling date) when a legal priority date was not available. The most questionable dates come from New Mexico, where the only date associated with certain water right records was the date the allocation was recorded in the database. As the allocation record creation tended to be within a few months of the filing of the application of the water right, from manually double checking the water rights, and our analysis focuses on aggregating water rights on the timescale of years, we determined it was acceptable to use such dates to include as many records as possible. primaryBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories WestDAAT. This column is the original WaDE category for the primary water use at the PoD site. allocationBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories for WestDAAT. This column is the original WaDE category

Economics↗
DOE OSTI2025

Ordinary and exotic mesons in the extended Linear Sigma Model

The extended Linear Sigma Model (eLSM) is a hadronic model based on the global symmetries of QCD and the corresponding explicit, anomalous, and spontaneous breaking patterns. In its basic three-flavor form, its mesonic part contains the dilaton/glueball as well as the nonets of pseudoscalar, scalar, vector, and axial–vector mesons, thus chiral symmetry is linearly realized. In the chiral limit and neglecting the chiral anomaly, only one term – within the dilaton potential – breaks dilatation invariance, and all terms are chirally symmetric. Spontaneous symmetry breaking is implemented by a generalization of the Mexican-hat potential, with explicit symmetry breaking responsible for its tilting. The overall mesonic phenomenology up to ~2 GeV is in agreement with the PDG compilation of masses and partial and total decay widths. The eLSM was enlarged in a straightforward way to include other conventional quark–antiquark nonets (pseudovector and orbitally excited vector mesons, tensor and axial-tensor mesons, radially excited (pseudo)scalar mesons, etc.), as well as two nonets of hybrid mesons, the lightest one with exotic quantum numbers J $\mathscr{PC}$ = 1 -+ not allowed for $\bar{q}q$ objects, such as the resonance π 1 (1600) and the recently discovered η 1 (1855) . In doing so, different types of chiral multiplets are introduced: heterochiral and homochiral multiplets, which differ in the way they transform under chiral transformations. Moreover, besides the scalar glueball that is present from the beginning as dilaton, other glueballs, the tensor, the pseudoscalar and the vector glueballs were coupled to the eLSM: the scalar resonance ƒ 0 (1710) turns out to be mostly gluonic, the tensor glueball couples strongly to vector mesons, and the pseudoscalar glueball couples sizably to $ππη'$ and can be assigned to X (2370) or X (2600). In all cases above, masses and decays can be analyzed allowing for a better understanding of both conventional and non-conventional mesons: whenever data are available, a comparison is performed and, when this is not the case, predictions of decay widths and decay ratios are outlined. The eLSM contains chiral partners on an equal footing and is therefore well suited for studies of chiral symmetry restoration at nonzero temperature and densities: this is done by coupling it to the Polyakov loop. Finally, the QCD phase diagram and the location of the critical endpoint were investigated within this framework.

Chiral symmetry↗
DOE OSTI2021

Level structure of the T z = –1 nucleus 34 Ar and its relevance for nucleosynthesis in ONe novae

The 24 Mg+ 12 C fusion reaction was used to perform a detailed γ-ray spectroscopy study of the astrophysically important nucleus 34 Ar. In particular, an experimental setup, coupling the advanced γ-ray tracking array GRETINA with the well-established Argonne fragment mass analyzer (FMA), was employed to obtain excitation energies and spin-parity assignments for excited states in 34 Ar, both above and below the proton separation energy. For the first time, an angular distribution analysis of in-beam γ rays from fusion-evaporation reactions, using a tracking array, has been performed and Coulomb energy differences of analog states in the T = 1, A = 34 mirror system, explored from 0 to 6 MeV. Lastly, we present a comprehensive discussion of the astrophysical 33 Cl(p,γ) stellar reaction rate, together with implications for the identification of nova presolar grains from sulfur isotopic abundances.

20 ≤ A ≤ 38↗
DOE OSTI2022

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗