Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “databases)”

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

Estimating soil N 2 O emissions induced by organic and inorganic fertilizer inputs using a Tier-2, regression-based meta-analytic approach for U.S. agricultural lands

Consistent methods are essential for generating country and region-specific estimates of greenhouse gas (GHG) emissions used for reporting and policymaking. The estimates of direct N 2 O emissions from U.S. agricultural soils have primarily relied on the use of emission factors (EFs, Tier-1) and process-based models (Tier-3). However, Tier-1 estimates are relatively crude while Tier-3 calculations can be costly. This work addressed this gap by developing a Tier-2, regression-based approach by leveraging a meta-database containing 1883 field N 2 O observations together with environmental and management covariates from 139 studies. Our results estimated higher monthly soil N 2 O emissions (N 2 O m , kg N/ha) during the growing season (0.38) than the fallow period (0.15), highlighting the importance of considering measurement periods when utilizing meta-databases for analyzing N 2 O drivers. Significantly different N 2 O m were found for tillage practices (conventional > no-till: 0.42 > 0.27), fertilizer type (liquid > solid manure: 0.55 > 0.32), and soil texture (fine > coarse: 0.36 > 0.22). The comparisons of the influence of crop type and rotation, water management, and soil order on N 2 O emissions are complicated by regional data availability and interactions among different factors. Additionally, the finding that N 2 O emissions reported based on area (N 2 O m ), N input rate (EF), or yield can alter treatment rankings underscores the need to establish transparent criteria for rewarding or discouraging regionally-based management practices using N 2 O metrics. Finally, we show how General Linear Models (GLMs) can be used to estimate country and regional Tier-2 N 2 O m using a suite of covariates. Our GLMs identified tillage, water management, N input type and rate, soil properties, and elevation as the most influential covariates for the conterminous U.S. The limited accuracy of regional-scale GLMs, however, suggests the need to further improve the quality and availability of GHG and covariate data through concerted efforts in data collection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Critical mineral source potential from oil & gas produced waters in the United States

The volume of produced water, a by-product of oil & gas operations and other energy processes, has been growing across the United States (U.S.) along with the need to manage or recycle this wastewater. Produced water contains many naturally occurring elements of varying concentrations, including critical minerals which are essential to the clean energy transition. However, the current understanding of critical mineral concentrations in produced water and the associated volumes across the U.S. is limited. This study has assessed available databases and literature to gain insight into the presence and concentration of five high priority critical minerals, namely cobalt, lithium, magnesium, manganese, and nickel. The U.S. Geological Survey's National Produced Waters Geochemical Database was the main data source used for determining average critical mineral concentrations in produced water from the major oil and gas reservoirs in the U.S. The volumes of produced water for these major reservoirs were coupled with these concentrations to provide insights into where critical minerals are likely to have high abundance and therefore more recovery options. The analysis indicated the highest recovery potential for lithium and magnesium from produced water in the Permian basin and the Marcellus shale region. However, these assessments should be considered conservative due to the limited availability of reliable concentration data. Finally, it is expected more critical mineral recovery options could emerge with comprehensive characterization data from more recent and representative sources of produced water.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Adsorbate-induced adatom formation on Au-Cu bimetallic alloys and its possible consequences for CO 2 electroreduction

The adsorbate-induced formation of sub-nanometer clusters on transition-metal single crystals observed in previous high-pressure microscopic studies hinted at the in-situ formation of unique active sites even on large nanoparticle catalysts. We propose that the adatom formation energy can be used as an energetic descriptor for the initial step toward the adsorbate-induced metal-cluster formation process. This descriptor can be efficiently computed using density functional theory (DFT) calculations and applied for screening and identification of metal catalysts where this phenomenon may play an important role in generating active sites in-situ. As a proof of concept, here, we construct an adatom formation energy database for three Au x Cu y alloys (x:y = 3:1, 1:1, or 1:3) and eighteen adsorbates (H, C, N, O, F, S, Cl, Br, I, CH x , NH x (x = 1 – 3), CO, NO, and OH) commonly involved in catalytic reactions. The energetics of adatom formation were examined in all cases where the (111) terrace, (211) step-edge, and (874) kink were the sources of the adatom. We demonstrate that the presence of an adsorbate could alter not only the energetics for adatom formation but also the elemental nature of the preferred adatom being formed. Using our database, we identified promising systems which favor adsorbate-induced adatom formation under near-ambient conditions. Specifically, CO-induced adatom formation on all three Au-Cu alloy surfaces could occur under CO 2 electroreduction (CO 2 RR) conditions. This phenomenon offers a qualitative explanation for the experimentally observed CO 2 RR activity on Au-Cu alloy catalysts. As a result, our methodology offers an easily expandable and efficient approach for large-scale catalyst screening with regards to adatom/cluster formation under reaction conditions and provides insight into the possible nature of active sites on alloy catalysts from a novel perspective.

Active site↗
DOE OSTI2026

Evaluation and optimization of flow boiling frictional pressure drop correlations using the data from traditional and next-generation refrigerants in a micro-fin tube

This study presents an experimental evaluation and optimization of flow boiling frictional pressure drop correlations for conventional and next-generation refrigerants in a horizontal micro-fin tube, with particular emphasis on the newly emerging refrigerant blends R-454C and R-455A, for which pressure-drop data in enhanced tubes remain limited. Experiments were conducted with R-410A, R-454C, R-455A, R-134a, R-1234yf, and R-1234ze(E) in a copper micro-fin tube with an inner diameter of 8.468 mm, over mass fluxes ranging from 100 to 300 kg/(m²·s) depending on the refrigerant, and evaporation temperatures of 7, 12, and 14 °C. Frictional pressure gradients were determined from measured total pressure drops after subtracting acceleration pressure drop, and the resulting database was used to assess four existing models: Kuo and Wang (1996), Cavallini et al. (1997), Goto et al. (2001), and Diani et al. (2014). The measured frictional pressure gradient increased with vapor quality and mass flux for all refrigerants and increased further at lower evaporation temperatures, with the overall trend strongly related to liquid viscosity. Among the four correlations, the Goto et al. (2001) model provided the best overall agreement with the measured data before optimization. To further improve prediction accuracy, the Kuo and Wang (1996) and Goto et al. (2001) models were optimized using the complete experimental database. After optimization, both models reduced the overall mean absolute deviation to below 15%, while the optimized Goto et al. (2001) model maintained the best and most consistent overall performance. The results provide new pressure-drop data for next-generation refrigerants and demonstrate that parameter optimization can significantly enhance the applicability of existing micro-fin-tube correlations.

Hu, Yifeng [ORNL] (ORCID:0000000242875185)↗
DOE OSTI2021

MCL1 inhibition enhances the efficacy of docetaxel against airway-derived squamous cell carcinoma cells

Highlights: • Authentic squamous cell carcinoma (SqCC) cell lines of the lung are sparse. • HARA lung SqCC cells depend highly on MCL1 for survival in vitro and in vivo. • MCL1-silenced HARA cells decline in number upon docetaxel (DTX) exposure. • High MCL1 mRNA expression is linked with shorter survival in lung SqCC patients. • MCL1 depletion + DTX is highly effective against Detroit 562 pharyngeal SqCC cells. MCL1 is an anti-apoptotic BCL2 family member that is often overexpressed in various malignant tumors. However, few reports have described the role of MCL1 in squamous cell carcinoma (SqCC) derived from airways including the lung. In this study, we examined whether MCL1 could be a novel druggable target for airway-derived SqCC, for which effective molecular targeted drugs are unavailable. We searched the Kaplan–Meier Plotter database and found that high MCL1 mRNA expression was significantly associated with shorter survival in patients with lower airway (lung) or upper airway (head and neck) derived SqCC. We also explored the Expression Atlas database and learned that authentic lung SqCC cell lines expressing both TP63 and KRT5 mRNA were extremely sparse among the publicly available “lung SqCC cell lines”, with an exception being HARA cells. HARA cells were highly dependent on MCL1 for survival, and MCL1-depleted cells were not able to grow, and even declined in number, upon docetaxel (DTX) exposure in vitro and in vivo. Similar in vitro experimental findings, including those in a 3D culture model, were also obtained using Detroit 562 pharyngeal SqCC cells. These findings suggested that combined treatment with MCL1 silencing plus DTX appears highly effective against airway-derived SqCC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Mast1 mediates radiation-induced gastric injury via the P38 MAPK pathway

Radiation-induced gastric injury is a serious adverse effect and reduces the efficacy of radiotherapy treatment. However, the mechanisms underlying radiation-induced stomach injury remain unclear. Here, mouse stomach and gastric epithelial cells were irradiated with different doses of X-ray radiation. The results showed that radiation induced gastric injury in vivo and in vitro. Differentially expressed functional mRNAs in irradiation-induced gastric tissues were screened from the Gene Expression Omnibus (GEO) database. We found that the expression of microtubule-associated serine/threonine kinase 1 (Mast1) was downregulated in mouse gastric tissues and gastric epithelial cells after irradiation. Furthermore, functional assays showed that knockdown of Mast1 inhibited growth and promoted apoptosis in gastric epithelial cells, while overexpression of Mast1 protected gastric epithelial cells from radiation damage. Mechanistically, Mast1 negatively regulated radiation-induced injury in gastric epithelial cells by inhibiting the activation of P38. The apoptosis caused by knockdown of Mast1 in gastric epithelial cells could be partially reversed by the P38 inhibitor SB203580. Moreover, data from several gastric cancer cell lines and online databases revealed that Mast1 was not involved in the development of gastric cancer. Collectively, our findings demonstrated that Mast1 is essential for radiation-induced gastric injury, providing a promising prognostic and therapeutic target.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Towards fast and accurate predictions of radio frequency power deposition and current profile via data-driven modelling: applications to lower hybrid current drive

Three machine learning techniques (multilayer perceptron, random forest and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ( $n_{e0}$ , $T_{e0}$ , $I_p$ , $B_t$ , $R_0$ , $n_{\|}$ , $Z_{{\rm eff}}$ , $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$ ) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data augmentation for disruption prediction via robust surrogate models

The goal of this work is to generate large statistically representative data sets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student t process regression. We apply Student t process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via colouring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics and classic machine learning clustering algorithms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Discovery of Energy Storage Molecular Materials Using Quantum Chemistry-Guided Multiobjective Bayesian Optimization

Redox flow batteries (RFBs) are a promising technology for stationary energy storage applications due to their flexible design, scalability, and low cost. In RFBs, energy is carried in flowable redoxactive materials (redoxmers) which are stored externally and pumped to the cell during operation. Further improvements in the energy density of RFBs necessitates redoxmer designs with wider redox potential windows and higher solubility. Additionally, designing redoxmers with a fluorescence-enabled self-reporting functionality allows monitoring of the state of health of RFBs. To accelerate the discovery of redoxmers with desired properties, state-of-the- art machine learning ( ML) methods, such as multiobjective Bayesian optimization (MBO), are useful. Here, we first employed density functional theory calculations to generate a database of reduction potentials, solvation free energies, and absorption wavelengths for 1400 redoxmer molecules based on a 2,1,3-benzothiadiazole (BzNSN) core structure. From the computed properties, we identified 22 Pareto-optimal molecules that represent best trade-off among all of the desired properties. We further utilized these data to develop and benchmark an MBO approach to identify candidates quickly and efficiently with multiple targeted properties. With MBO, optimal candidates from the 1400-molecule data set can be identified at least 15 times more efficiently compared to the brute force or random selection approach. Importantly, we utilized this approach for discovering promising redoxmers from an unseen database of 1 million BzNSN-based molecules, where we discovered 16 new Pareto-optimal molecules with significant improvements in properties over the initial 1400 molecules. Furthermore, we anticipate that this active learning technique is general and can be utilized for the discovery of any class of functional materials that satisfies multiple desired property criteria.

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

Bubble Point Measurements of Three Binary Mixtures of Refrigerants: R-32/1234yf, R-32/1234ze(E), and R-1132a/1234yf

In this article, bubble point pressures of three binary mixtures at two compositions each have been measured utilizing a static method. The mixtures studied were difluoromethane + 2,3,3,3-tetrafluoroprop-1-ene (R-32/1234yf), difluoromethane + (E)-1,3,3,3-tetrafluoroprop-1-ene (R-32/1234ze(E)), and 1,1-difluoroethene + 2,3,3,3-tetrafluoroprop-1-ene (R-1132a/1234yf). For all of the mixtures, the minimum temperature measured was 270 K. The maximum temperatures measured for the R-32/1234ze(E), R-32/1234yf, and R-1132a/1234yf mixtures were 360, 345, and 325 K, respectively. A total of 100 bubble point pressures are reported and compared to available literature data relative to the REFPROP database. Lemmon, E. W. [ NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, 10.0; National Institute of Standards and Technology, Standard Reference Data Program: Gaithersburg, 2018]. Additionally, the bubble point pressures obtained in this study are modeled using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state. The PC-SAFT model parameters, as well as deviations of this work and literature data from the model, are reported.

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

Reaction Mechanism Generator v3.0: Advances in Automatic Mechanism Generation

In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.

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

Acceleration of Graph Neural Network-Based Prediction Models in Chemistry via Co-Design Optimization on Intelligence Processing Units

Atomic structure prediction and associated property calculations are the bedrock of chemical physics. Since high-fidelity ab initio modeling techniques for computing the structure and properties can be prohibitively expensive, this motivates the development of machine-learning (ML) models that make these predictions more efficiently. Training graph neural networks over large atomistic databases introduces unique computational challenges such as the need to process millions of small graphs with variable size and support communication patterns that are distinct from learning over large graphs such as social networks. We demonstrate a novel hardware-software co-design approach to scale up the training of atomistic graph neural networks (GNN) for structure and property prediction. First, to eliminate redundant computation and memory associated with alternative padding techniques and to improve throughput via minimizing communication, we formulate the effective coalescing of the batches of variable-size atomistic graphs as the bin packing problem and introduce a hardware-agnostic algorithm to pack these batches. In addition, we propose hardware-specific optimizations including a planner and vectorization for the gather-scatter operations targeted for Graphcore’s Intelligence Processing Unit (IPU), as well as model-specific optimizations such as merged communication collectives and optimized softplus. Putting these all together, we demonstrate the effectiveness of the proposed co-design approach by providing an implementation of a well-established atomistic GNN on the Graphcore IPUs. We evaluate the training performance on multiple atomistic graph databases with varying degrees of graph counts, sizes and sparsity. Here, we demonstrate that such a co-design approach can reduce the training time of atomistic GNNs and can improve the performance by up to 1.5× compared to the baseline implementation of the model on the IPUs. Additionally, we compare our IPU implementation with a Nvidia GPU-based implementation and show that our atomistic GNN implementation on the IPUs can run 1.8× faster on average compared to the execution time on the GPUs.

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

Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal–Organic Frameworks

Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO 2 . Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal–Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. In conclusion, a strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.

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

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2024

Robust Tensor Hypercontraction of the Particle–Particle Ladder Term in Equation-of-Motion Coupled Cluster Theory

One method of representing a high-rank tensor as a (hyper-)product of lower-rank tensors is the tensor hypercontraction (THC) method of Hohenstein et al. This strategy has been found to be useful for reducing the polynomial scaling of coupled-cluster methods by representation of a four-dimensional tensor of electron-repulsion integrals in terms of five two-dimensional matrices. Pierce et al. have already shown that the application of a robust form of THC to the particle–particle ladder (PPL) term reduces the cost of this term in couple-cluster singles and doubles (CCSD) from O(N 6 ) to O(N 5 ) with negligible errors in energy with respect to the density-fitted variant. In this work, we have implemented the least-squares variant of THC (LS-THC) which does not require a nonlinear tensor factorization, including the robust form (R-LS-THC), for the calculation of the excitation and electron attachment energies using equation-of-motion coupled cluster methods EOMEE-CCSD and EOMEA-CCSD, respectively. We have benchmarked the effect of the R-LS-THC-PPL approximation on excitation energies using the comprehensive QUEST database and the accuracy of electron attachment energies using the NAB22 database. Here, we find that errors on the order of 1 meV are achievable with a reduction in total calculation time of approximately 5x.

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

Confirmation of PNNL Quantitative Infrared Cross-Sections for Isobutane

The Pacific Northwest National Laboratory (PNNL) gas-phase database is a compilation of quantitative experimental (5, 25, 50 °C) infrared spectra of ca. 500 molecules, designed for in situ, standoff or remote sensing of gases and vapors at or near atmospheric pressure. The data are characterized by calibration on both the wavenumber and intensity axes. Recent papers have called into question the PNNL intensity values for isobutane, (2-methylpropane, HC(CH 3 ) 3 ) suggesting discrepancies of 30 to 40%. In this study we re-measure and re-examine the intensity values of isobutane using both similar and alternate methods to those used to generate the original PNNL database spectra. Here, indirect confirmation from literature data of homologous molecules and direct confirmation from new results confirm that for many band integrals across the isobutane spectrum, the original PNNL data are indeed accurate to within the reported 3% experimental uncertainty.

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

R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

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

Modeling Nanoconfinement Effects Using Active Learning

Predicting the spatial configuration of gas in nanopores of is relevant in applications such as fluid flow forecasting and hydrocarbon reserves estimation. For example, shale reservoirs have suffered from computationally intractable multiscale problems, since fluid properties such as viscosity, density, and adsorption must be calculated by using expensive molecular dynamics (MD) simulations within each nanopore, whereas flow through these connected nanopores must be simulated at the micrometer scale. We utilize machine learning techniques to quickly and accurately model nanoscale confinement effects as an important step toward bridging the nano and micro scales. Our workflow is based on building and training physics-based deep-neural-networks models by learning from a database of MD calculations. The model accounts for the adsorption phenomenon by predicting the statistical distribution of gas inside nanopores. Because large databases of MD calculations are expensive to create, we investigate active learning (AL) as a data set construction strategy. In this workflow, new data are selected based on the model uncertainty via the query-by-committee approach. We show that our workflow obtains accurate models that generalize to real scanning electron microscopy geometries with 1/10th of the number of MD calculations required vs random data set generation. Our method enables the possibility of modeling nanoconfinement effects at the mesoscale, where complex connected sets of nanopores affect flow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗