Search NASA⌕ Search

SEARCH · Search NASA

Results for “function approximation”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Fabrication of Catalytic Distillation Membranes with Atomic Layer Deposition

The integration of catalysts onto the surface of membranes enables simultaneous physical separation and catalytic transformation of constituents in a feed stream, facilitating improved contaminant removal and fouling mitigation. Distillation membranes are a particularly attractive platform for catalytic membranes because they reject nonvolatile species and exhibit exceptional resistance to oxidative and radical-driven degradation. However, imparting catalytic functionality onto hydrophobic, porous distillation membranes has proven challenging since the membranes used are chemically inert and difficult to modify. Furthermore, catalysts on the membrane surface can decrease hydrophobicity and increase the membrane’s susceptibility to pore wetting and failure. In this work, we create a catalytic distillation membrane by coating a polytetrafluoroethylene membrane surface with titanium dioxide (TiO 2 ) via plasma-assisted atomic layer deposition (ALD). By precisely tuning the ALD parameters, we demonstrate localized growth of TiO 2 near (within approximately 1 μm) the surface of polytetrafluoroethylene membranes, forming a catalytically active interface while preserving the underlying hydrophobic pore structure. Localized growth of TiO 2 is confirmed by electron microscopy and spectroscopy techniques, and membranes coated with 500 cycles of ALD show pressure tolerance up to 12.8 bar and higher than 95% salt rejection in pressure-driven distillation. Photocatalytic activity is demonstrated via the degradation of methylene blue dye under UV irradiation, where increasing TiO2 loading leads to an enhancement in dye degradation. These results establish a general strategy for integrating catalytic functionality into chemically inert, hydrophobic membranes without compromising distillation performance, providing a pathway toward multifunctional membranes that couple advanced oxidation with membrane separation for water treatment.

atomic layer deposition↗
DOE OSTI2026

Examination of Replicate Syntheses of Metal Organic Frameworks as a Window into Reproducibility in Materials Chemistry

Replicate experiments are a useful tool in understanding the repeatability of scientific measurements. In 2019, a systematic search for replicate syntheses of a collection of 130 metal–organic frameworks (MOFs) found that 89% of these materials had no reported replicate syntheses apart from the original publications identifying the material (Agrawal, M. Proc. Natl. Acad. Sci. U.S.A. 2020, 117, 877−88210.1073/pnas.1918484117). A potential weakness of that search was that only 5–11 years had elapsed since the original publication of each material. Here, this analysis is extended to all publications 11–17 years after the original publication. Although this extended time period identifies more repeat syntheses, 83% of the materials still have no reported replicate syntheses. We also consider how appropriately selected Density Functional Theory (DFT) calculations can provide corroboration for the experimentally reported crystal structures. By using data from previous high-throughput DFT studies, corroborating evidence from DFT was available for 17% of the 130 structures for which no replicate syntheses are available. In total, approximately 1/3 of the 130 MOFs have data associated with replicate synthesis experiments and/or directly corroborating DFT calculations.

Sholl, David S. [Oak Ridge National Laboratory (OR↗
DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗
DOE OSTI2025

Search for the jet-induced diffusion wake in the quark-gluon plasma via measurements of jet-track correlations in photon-jet events in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector

This paper presents a measurement of jet-track correlations in photon-jet events, using 1.72 nb −1 of Pb+Pb data at $\sqrt{s_{NN}}$ = 5.02TeV recorded with the ATLAS detector at the LHC. Events with energetic photon-jet pairs are selected, where the photon and jet are approximately back-to-back in azimuth. The angular correlation between jets and charged-particle tracks with transverse momentum (𝑝 T ) in the range 0.5–2.0 GeV in the hemisphere opposite to the jet, |Δ⁢𝜙(jet,track)|>𝜋/2, is measured as a function of their relative pseudorapidity difference, |Δ⁢𝜂(jet,track)|. In central Pb+Pb collisions, these correlations are predicted to be sensitive to the diffusion wake in the quark-gluon plasma resulting from the lost energy of high-𝑝T partons traversing the plasma, with a characteristic modification as a function of |Δ⁢𝜂(jet,track)|. The correlations are examined with different selections on the jet-to-photon 𝑝 T ratio to select events with different degrees of energy loss. No diffusion wake signal is observed within the current sensitivity and upper limits at 95% confidence level on the diffusion wake amplitude are reported.

heavy flavor physics↗
DOE OSTI2025

Conditional Pseudo-Reversible Normalizing Flow for Surrogate Modeling in Quantifying Uncertainty Propagation

We introduce a conditional pseudo-reversible normalizing flow (PR-NF) that directly learns conditional probability distributions from noisy physical models to efficiently quantify both forward and inverse uncertainty propagation. Traditional surrogate modeling approaches approximate only the deterministic component of physical models, requiring separate noise characterization and computationally expensive sampling methods for inverse problems. Here, in this work, we develop the conditional PR-NF model to directly learn and efficiently generate samples from the conditional probability density functions (PDFs). The training process utilizes dataset consisting of input-output pairs without requiring prior knowledge about the noise and the function. Once trained, our model efficiently generates samples from conditional PDFs for any input within the training domain. Moreover, the pseudo-reversibility feature allows for the use of fully connected neural network architectures, which simplifies the implementation and enables theoretical analysis. We provide a rigorous convergence analysis of the conditional PR-NF model, showing its ability to converge to the target conditional PDF using the Kullback−Leibler divergence. To demonstrate the effectiveness of our method, we apply it to several benchmark tests and a real-world geologic carbon storage problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗
DOE OSTI2025

Universality in diffusion-controlled nucleation and growth

Nucleation and growth are studied in a system that undergoes diffusion-controlled condensation under gradual changes in parameters, such as cooling. Here, it is demonstrated that when the Gibbs–Thompson effect becomes negligible, the system falls into a universal regime. i.e., the final droplet size distribution remains invariant under certain rescaling of system parameters. An approximate yet very accurate analytic form is obtained for the final droplet size distribution in this regime.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Plutonium Oxidation State Distribution in the Presence of WIPP-Relevant Organics and Iron Corrosion Products

The oxidation state and solubility of plutonium (Pu) in high ionic strength synthetic WIPP (Waste Isolation Pilot Plant) brines as a function of pC H+ in the presence and absence of WIPP-relevant organic ligands (EDTA [Ethylenediaminetetraacetic acid], oxalate, citrate, acetate) and iron corrosion products (magnetite and metallic iron) at 𝑇 = 23 ± 2 ∘C was thoroughly studied by long-term batch solubility experiments (between approximately 800-1,100 days) from an undersaturation approach. The oxidation state of Pu in the WIPP environment has been a topic of interest since the initial Compliance Certification Application (CCA). This study aims to investigate the solubility of Pu under the expected WIPP conditions and to determine the oxidation state of the solid phase that will control the solubility. One of the most important results of this study is that the Pu oxidation state was analyzed both from the surface area of the corrosion products and in the precipitated solid. The analysis of the Pu oxidation state on the surface of the iron mineral from the ongoing undersaturated experiments is more relevant to the performance assessment of nuclear waste disposal than short term batch (plutonium-iron phase) experiments. The X-ray Absorption Near-Edge Spectroscopy (XANES) analysis showed that Pu oxidation state is different on the metal surface and in the precipitated solid. Pu(III) is the dominant oxidation state in the metallic iron (Fe 0 ) system in the presence and absence of organics. Pu(IV) is the dominant oxidation state in the magnetite system in the presence and absence of organics. Also, organics stabilize Pu(IV) in the magnetite system. Analysis of Pu in the precipitated solid by Extended X-ray Absorption Fine Structure (EXAFS) analysis showed that Pu formed an inner-sphere complex with iron with minor amounts of PuO 2 present. X-ray diffraction (XRD) results indicate that metallic iron and magnetite did not oxidize in three years in the alkaline and high ionic strength system. Under these conditions (8 < pC H+ < 10 at T = (22 ± 2) °C under nitrogen atmosphere, the solubility of Pu changes by up to three orders of magnitude (10 -5 and 10 -8 M). The spread in solubility is highest at pC H+ = 9. Pu(III) and Pu(IV) showed different solubility behavior in the synthetic WIPP brine. A summary of the data collected in this report will be submitted to Sandia National Laboratories as part of a parameter update report which will outline the changes to the OXSTAT parameter for the 2026 Compliance Recertification Application (CRA-2026). The experiments performed were done according to the U.S. Department of Energy (DOE) approved Test Plan entitled “Effects of Radiolysis, Organic Complexation, and Redox Conditions on the Speciation and Oxidation State Distribution of Pu(III/IV)” (LCO-ACP-25). All data reported were obtained under the Los Alamos National Laboratory-Carlsbad Office (LANL-CO) Quality Assurance Program, which is compliant with the DOE Carlsbad Field Office, Quality Assurance Program Document (CBFO/QAPD).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Quasielastic $\overrightarrow{^{3}He}$ $(\overrightarrow{e},e′)$ asymmetry in the threshold region with a longitudinally polarized target

We report on the measurements of the double-spin asymmetry from electron- 3 He scattering in the threshold region of two- and three-body breakup of 3 He for 4-momentum transfer values of 0.1 and 0.2 (GeV/c) 2 . The measurement was done at Jefferson Lab using a polarized 3 He target with its spin aligned parallel to the beam direction. This asymmetry is sensitive to a combination of spin-dependent response functions not previously measured and may highlight differences in two body forces, meson-exchange currents and final-state interactions. The results of this measurement serve as a test of our understanding of few-body systems. When compared with calculations from plane wave impulse approximation and Faddeev calculations, we found that the latter, which use modern nuclear potentials and prescriptions for meson-exchange currents, demonstrate an overall good agreement with data.

Few-body nuclei↗
DOE OSTI2026

Design Trade-Offs in Composite Fuel Cell Membranes: Effects of Reinforcement and Chemical Additives

Perfluorosulfonic acid (PFSA) membranes are critical components in proton exchange membrane fuel cells, where performance depends on balancing ionic conductivity, mechanical durability, and chemical stability. This study characterizes a composite membrane (NC700) featuring PFSA-impregnated expanded polytetrafluoroethylene (ePTFE) reinforcement and cerium-based radical scavengers, benchmarked against unreinforced NR211. Complementary techniques, including electron microscopy, X-ray scattering, infrared spectroscopy, thermogravimetric analysis, and dynamic mechanical analysis, identify the structural and compositional strategies employed in NC700. Water sorption isotherms reveal lower water uptake for NC700 across all conditions, attributed to reinforcement and cerium incorporation. Reinforcement reduces in-plane swelling from 11% to 2.1% at 90% RH, confirming strong swelling anisotropy, while maintaining mechanical properties at elevated temperatures. While the ionic conductivity of NC700 is approximately 10% lower than that of NR211, the reduced thickness yields a 40% decrease in calculated area-specific resistance, suggesting the composite architecture can favorably shift the conductivity-stability trade-off. The composite structure also reduces gas permeability, indicating potential for improved separator function alongside favorable transport properties. Systematic deconvolution of reinforcement and additive contributions shows that conductivity losses from cerium incorporation are largely offset by gains from the lower equivalent-weight polymer, providing quantitative relationships that may guide composite membrane design for fuel cells and other electrochemical applications.

25 ENERGY STORAGE↗
DOE OSTI2024

Computational Exploration of High Entropy Alloys as Promising Materials for Future Beam Windows

With the ever-increasing demand for high beam power, the currently used beam-intercepting devices (BIDs) such as targets, and beam windows may not be able to handle the high power required for future accelerator complexes or the lifetime may be reduced drastically. As beam power increases, the damage incurred by BIDs, including thermal shock, fatigue, and irradiation damage, also rises. Therefore, it is imperative to design materials that can withstand high beam power for longer lifetimes. High entropy alloys (HEAs) have emerged as potential alternative materials for designing next-generation BIDs. In this study, we primarily focus on materials for developing beam windows for next-generation accelerator complexes. We propose an integrated approach that combines various computational techniques to study and design new materials. Specifically, we use CALPHAD, density functional theory (DFT), and molecular dynamics (MD) to comprehensively investigate the defect properties of suitable HEAs, offering potential alternatives for future beam windows. We begin by scanning the extensive phase space provided by Cr-Mn-V-Ti-Al-Co HEAs, selecting 8 compositions after evaluating approximately 120,000 unique compositions using CALPHAD. We, then employ DFT-informed machine learning techniques to develop force-field parameters. Finally, MD simulations using these developed force-field parameters will be used to study the effects of radiation damage on the defect and mechanical properties of the selected alloys. This research explains the use of the CALPHAD approach and shows how critical modeling (DFT and MD) is in developing novel material such as HEAs. It also highlights the promising role of machine learning in this field. The results from this study will greatly improve the novel materials development to be used in next-generation accelerator components, leading to higher beam power and longer operational times of BIDs.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

The AMACStar ASIC for the HL-LHC ATLAS ITk Strip detector: design, verification, testing, and quality assurance

For the high-luminosity upgrade to the LHC (HL-LHC), the ATLAS detector at CERN requires an all-new inner detector, the Inner Tracker (ITk). The ITk Strip subdetector is made up of silicon modules, which include three types of radiation-hard ASICs. One of these is the Autonomous Monitor and Control (AMAC). The AMAC is manufactured by Global Foundries using 130 nm CMOS8RF DM technology and is approximately 3 by 5 mm in size. This ASIC autonomously monitors the temperatures, voltages, and currents in the module components while controlling critical values in order to prevent these quantities from reaching dangerous levels. The final design, AMACStar, was verified, tested, and ensured to perform all necessary functions. A quality control procedure using an in-house probe station set-up was developed in order to ensure that every individual chip on each wafer of AMACStars required by the ITk Strip project met performance requirements. The average per wafer yield of usable AMACStars is 92.33%, exceeding the design-specific 90% yield estimated for project costing. This estimate was based on actual yields for similar designs in this process.

Analogue electronic circuits↗
DOE OSTI2024

Impact of penetrating collisions of plasma ions on spectral line shapes

Spectral-line-broadening models have been moving towards including full Coulomb interactions between the atom and plasma, replacing the commonly used dipole approximation. The effects of the full Coulomb interaction have been thoroughly explored for plasma electrons, resulting in redshifts of spectral lines in high-energy-density plasmas. We explore the impact of a full Coulomb treatment on ion broadening. Penetrating collisions due to ions do not significantly impact the linewidths. The most significant aspect is the appearance of quasimolecular resonances in the far line wings, such as those previously observed in white-dwarf spectra. We identify several problems with the existing models. The most direct implementation causes the quasimolecular features to appear at the wrong photon energies. It is possible to include the detailed molecular structure, but this currently requires the Anderson-Talman approximation, which ignores the N-body properties of the plasma. We find that N-body effects can substantially broaden these quasimolecular features and even shift them, depending on the plasma conditions. We conclude that penetration of ions into the spatial extent of the radiator wave function does not strongly affect the usual diagnostics and may have a moderate to weak effect on opacity. Furthermore, the Rosseland mean opacity is weighted towards low-opacity regions of a spectrum, where these quasimolecular features are found.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nonequilibrium steady-state thermoelectrics of Kondo-correlated quantum dots

The transport across a Kondo-correlated quantum dot coupled to two leads with independent temperatures and chemical potentials is studied using a controlled nonperturbative, and in this sense numerically exact, treatment based on a hybrid numerical renormalization group combined with time-dependent density matrix renormalization group (NRG-tDMRG). In the Kondo regime, for sufficiently large fixed voltage bias V ≳ T K , with T K the Kondo temperature, we find a peak in the conductance vs the temperature gradient Δ⁢T = T R - T L across left and right lead. Focusing then on zero voltage bias but finite ΔT far beyond linear response, we reveal the dependence of the characteristic zero-bias conductance on the individual lead temperatures. Here, we find that the finite-Δ⁢T data behaves quantitatively similar to linear response with an effective equilibrium temperature derived from the different lead temperatures. The regime of sign changes in the Seebeck coefficient, signaling the presence of Kondo correlations, and its dependence on the individual lead temperatures provide a complete picture of the Kondo regime in the presence of finite-temperature gradients. The results from the zero-bias conductance and Seebeck coefficient studies unveil an approximate “Kondo circle” in the T L /T R plane as the regime within which the Kondo correlations dominate. We also study the heat current and the corresponding heat conductance vs finite Δ⁢T. We provide a polynomial fit for our numerical results for the thermocurrent as a function of the individual lead temperatures, which may be used to fit experimental data in the Kondo regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Anomalous elastic softening in ferroelectric hafnia under pressure

his study employs first-principles density-functional theory (DFT) calculations to explore the elastic and mechanical properties of ferroelectric hafnia (HfO 2 ) in its polar orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase under varying hydrostatic pressure conditions up to 30 GPa. Utilizing a plane-wave basis set and Perdew-Burke-Ernzerhof generalized-gradient approximation for solids in our DFT calculations, we investigate both pure and yttrium-substituted HfO 2 . Our findings reveal an anomalous reduction in the 𝐶 33 component of the elastic tensor with increasing pressure, which becomes significant above 15 GPa and signals a potential pressure-driven structural phase transition at higher pressure. The analysis of atomic displacements under pressure sheds light on the unusual mechanical behavior and phase stability of this material. Additionally, we observe a transition from an indirect band gap to a direct band gap with increasing pressure, which could have significant implications for optical applications. Here, the effects of yttrium substitution on the mechanical and electronic properties are further examined, revealing that yttrium substitution softens the elastic response of this material and reduces the electronic band gap. These results enhance our understanding of elastic and mechanical responses of ferroelectric hafnia and its potential for applications in microelectronics, piezoelectric devices, and nonvolatile ferroelectric random-access memories. Further experimental validation is recommended to confirm our predictions and explore the practical implications of the observed phase transitions and electronic behavior of the ferroelectric hafnia under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗
DOE OSTI2025

Precision Reconstruction of Rational Conformal Field Theory from Exact Fixed-Point Tensor Network

The novel concept of entanglement renormalization and its corresponding tensor network renormalization technique have been highly successful in developing a controlled real-space renormalization group (RG) scheme. Numerically approximate fixed-point (FP) tensors are widely used to extract the conformal data of the underlying conformal field theory (CFT) describing critical phenomena. In this paper, we present an explicit analytical construction of the FP tensor for 2D rational CFT. We define it as a correlation function between the “boundary-changing operators” (BCO) on triangles. Our construction fully captures all the real-space RG conditions. We also provide concrete examples, such as Ising, Yang-Lee, and tricritical Ising models, to compute the scaling dimensions explicitly based on the corresponding FP tensor. The BCO descendants turn out to be an optimal basis such that truncation in bond dimensions naturally produces comparable accuracies with the leading existing FP algorithms. Interestingly, our construction of FP tensors is closely related to a strange correlator, where the holographic picture naturally emerges. Our results also open a new door toward understanding CFT in higher dimensions. Published by the American Physical Society 2025

Cheng, Gong (ORCID:0009000891587404)↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗