Search NASA⌕ Search

SEARCH · Search NASA

Results for “functions and functionals”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Function-based Taxonomy Implementation

Function-based taxonomies are relational mapping tools used to clearly illustrate how an organization delivers Critical Functions by using various people, processes, technologies, information, and infrastructure (PPTII). Critical Functions are the actions or activities that make up the organization’s primary purpose. Enabling Functions are the combination of PPTII used by the organization to deliver their Critical Functions. A function-based taxonomy provides a framework to categorize information and related artifacts that document an organization’s unique implementation of PPTII for Critical Function delivery. The ideal state of a function-based taxonomy is to organize the existing knowledge the organization already has throughout their numerous systems, policies, people, procedures, and configurations. A well-organized taxonomy helps organizations to effectively leverage their knowledge by clearly identifying dependencies and connections. This document is the result of combined engineering and analytical experience and describes a repeatable method for producing function-based taxonomies.

42 ENGINEERING↗
DOE OSTI2020

Safety Functions and Features, Events and Processes for the E-Area Performance Assessment

The DOE Technical Standard, “Disposal Authorization Statement and Tank Closure Documentation,” (DOE 2017) recommends the use of safety functions and features, events and processes (FEPs) to support development of conceptual models and identification of scenarios to be considered in a performance assessment (PA). The FEP process provides a means to describe how a PA considers and addresses the factors that could influence the performance of key barriers. Understanding the roles of barriers in terms of limiting migration helps to focus on how changes in the system could lead to a situation where those roles cannot be fulfilled and there is the potential for compromised performance. The FEPs screening and review process was used to identify FEPs that are relevant for the EArea Low-Level Waste Facility (LLWF) and specifically those FEPs that could have a detrimental impact on the effectiveness of a given safety function. For this PA, a default list of FEPs developed at the International Atomic Energy Agency (IAEA 2004) and an approach implemented for PAs at the Hanford and Idaho sites (Mehta et al. 2016, DOE-ID 2019) are used to identify processes and events that could influence the effectiveness of a given safety function for the E-Area LLWF (e.g., subsidence can impact the safety function of the cover system and lead to increased infiltration). The Hanford and Idaho PAs represent two of the most recent applications of this approach. The PA evaluates the potential impacts of changes in performance of different features of the system and demonstrates that the safety functions represent multiple and redundant barriers. Barrier analyses, assuming a safety function is not present, also test the robustness of the system in the event of the loss of one or more safety functions. Such evaluations also support a qualitative illustration of the concept of defense in depth. The safety concept for closure of the E-Area LLWF (generically referred to as “E-Area”) encompasses a variety of different features (i.e., administrative controls, natural site features, and engineered barriers) that reduce the potential impacts on human health and the environment from the residual waste that will remain after closure. These features can be represented as a collection of safety functions acting independently and as a system to provide for overall safety. In some applications, there have been attempts to assign numerical expectations to specific safety functions, but that is not the intent in this case. The concept of safety functions is used more qualitatively in two ways for this PA: 1. To illustrate the robustness of the E-Area design, operational practices and closure approach by documenting features that are and are not credited in different modeling cases. 2. To identify the roles of the different features and potential processes and events that could compromise the performance of safety features and need to be considered when developing the modeling approach. This report addresses both safety functions and FEPs for the E-Area PA.

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

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
DOE OSTI2025

Fully thermal meta-GGA exchange correlation free-energy density functional

The application of density functional theory to materials in the warm dense matter regime has motivated the development of exchange-correlation functionals which incorporate proper, explicit temperature dependence. Previous work has yielded fully-thermal exchange-correlation free energy functionals at the local density approximation (LDA) and generalized gradient approximation (GGA) levels of refinement. Recently an additive thermal correction scheme was utilized to construct a meta-GGA exchange-correlation (XC) functional in which thermal effects are treated at the GGA level. Here, the f TSCAN free-energy XC functional presented here includes thermal effects through the meta-GGA level in the context of the SCAN (strongly constrained and appropriately normed) ground-state functional. The f TSCAN functional provides generality while achieving similar performance to a thermal GGA functional at high temperatures, e.g. pressures within 1% of path integral Monte Carlo simulations of warm dense hydrogen, and a significant improvement over ground-state functionals. At low temperatures, f TSCAN demonstrates improvements in accuracy relative to lower-level and deorbitalized functionals, indicating that calculations using f TSCAN may be expected to perform well across experimentally relevant densities and pressures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional

The absence of a reliable formulation of kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions which have been proposed in the literature and the appropriate coefficients are calculated via a linear regression framework. To model the nonlocal contributions, we explore two new nonlocal functionals - a functional that captures fluctuations in electronic density and a functional that incorporates gradient information. Since, the analytical functional forms of the kernels present in these nonlocal terms are not known from theory, we propose a basis function expansion to model these seemingly difficult nonlocal quantities. This allows us to easily reconstruct kernels for any system using only a few structures. The proposed method is able to learn kinetic energy densities and total kinetic energies of molecular and periodic systems, such as H 2 , LiH, LiF and a one-dimensional chain of 8 hydrogens using data from Kohn-Sham density functional theory calculations for only a few structures. For the ease of reproduction, codes used to generate the models are provided in the supporting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

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

funcX: Federated Function as a Service for Science

Here, funcX is a distributed function as a service (FaaS) platform that enables flexible, scalable, and high performance remote function execution. Unlike centralized FaaS systems, funcX decouples the cloud-hosted management functionality from the edge-hosted execution functionality. funcX's endpoint software can be deployed, by users or administrators, on arbitrary laptops, clouds, clusters, and supercomputers, in effect turning them into function serving systems. funcX's cloud-hosted service provides a single location for registering, sharing, and managing both functions and endpoints. It allows for transparent, secure, and reliable function execution across the federated ecosystem of endpoints-enabling users to route functions to endpoints based on specific needs. funcX uses containers (e.g., Docker, Singularity, and Shifter) to provide common execution environments across endpoints. funcX implements various container management strategies to execute functions with high performance and efficiency on diverse funcX endpoints. funcX also integrates with an in-memory data store and Globus for managing data that may span endpoints. We motivate the need for funcX, present our prototype design and implementation, and demonstrate, via experiments on two supercomputers, that funcX can scale to more than 130000 concurrent workers. We show that funcX's container warming-aware routing algorithm can reduce the completion time for 3,000 functions by up to 61% compared to a randomized algorithm and the in-memory data store can speed up data transfers by up to 3x compared to a shared file system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Predation probabilities and functional responses: How piscivorous waterbirds respond to pulses in fish abundance

Abstract How predators respond to changes in prey abundance (i.e., functional responses) is foundational to consumer–resource interactions, predator–prey dynamics, and the stability of predator–prey systems. Predation by piscivorous waterbirds on out‐migrating juvenile steelhead trout ( Oncorhynchus mykiss ) is considered a factor affecting the recovery of multiple Endangered Species Act‐listed steelhead populations in the Columbia River basin. Waterbird functional responses, however, may vary by predator species and location, with important implications to predator management strategies. We used a 13‐year dataset on waterbird abundance across seven breeding colonies (three Caspian tern [ Hydroprogne caspia ], two double‐crested cormorant [ Nannopterum auritum ], and two California and ring‐billed gull [ Larus californicus and Larus delawarensis ] colonies) and steelhead tag‐recovery data (>645,000 tagged and >32,000 recovered steelhead) to quantify weekly predation probabilities and functional responses across waterbird species, colonies, and years. Weekly predation probabilities were highly variable, ranging from 0.01 to 0.30 at tern colonies, 0.01 to 0.20 at cormorant colonies, and 0.03 to 0.13 at gull colonies. Per capita predation probabilities were an order of magnitude higher at inland tern and cormorant colonies relative to estuary colonies of the same species. Terns displayed Type II functional responses across colonies and years, where predation probabilities peaked at low steelhead abundances and declined as steelhead abundance increased (i.e., predator swamping). Cormorants nesting at the large estuary colony (several thousand birds) displayed a Type III functional response, but cormorants nesting at the smaller inland colony (several hundred birds) displayed a Type II response. Consumption probabilities of steelhead by gulls remained consistent across a large range of steelhead availability, suggesting a Type I or a Type III functional response, but a lack of colony abundance data prevented quantifying functional responses. The level of tern predation combined with Type II functional responses indicate possible population‐level impacts that could destabilize small or declining prey populations. Conversely, the apparent Type III functional responses of gulls and estuary nesting cormorants are indicative of prey switching behaviors targeted at periods of high steelhead abundance. Our results illustrate the complexity of predator–prey interactions and the importance of quantifying predator‐ and location‐specific functional responses when predicting the efficacy of management strategies to enhance prey populations.

Hostetter, Nathan J.↗
DOE OSTI2021

Covalent Functionalization of Nickel Phosphide Nanocrystals with Aryl-Diazonium Salts

Covalent functionalization of Ni 2 P nanocrystals was demonstrated using aryl-diazonium salts. Spontaneous adsorption of aryl functional groups was observed, with surface coverages ranging from 20 to 96% depending on the native reactivity of the salt as determined by the aryl substitution pattern. Increased coverage was possible for low reactivity species using a sacrificial reductant. Functionalization was confirmed using thermogravimetric analysis, Fourier transform infrared spectroscopy, and X-ray photoelectron spectroscopy. The structure and energetics of this nanocrystal electrocatalyst system, as a function of ligand coverage, were explored with density functional theory calculations. The Hammett parameter of the surface functional group was found to linearly correlate with the change in Ni and P core-electron binding energies and the nanocrystal's experimentally and computationally determined work function. The electrocatalytic activity and stability of the functionalized nanocrystals for hydrogen evolution were also improved when compared to the unfunctionalized material, but a simple trend based on electrostatics was not evident. Density functional theory was used to understand this discrepancy, revealing that H adsorption energies on the covalently functionalized Ni 2 P also do not follow the electrostatic trend and are predictive descriptors of the experimental results.

14 SOLAR ENERGY↗
DOE OSTI2025

A compendium of human gene functions derived from evolutionary modelling

A comprehensive, computable representation of the functional repertoire of all macromolecules encoded within the human genome is a foundational resource for biology and biomedical research. The Gene Ontology Consortium has been working towards this goal by generating a structured body of information about gene functions, which now includes experimental findings reported in more than 175,000 publications for human genes and genes in experimentally tractable model organisms 1,2 . Here, we describe the results of a large, international effort to integrate all of these findings to create a representation of human gene functions that is as complete and accurate as possible. Specifically, we apply an expert-curated, explicit evolutionary modelling approach to all human protein-coding genes. This approach integrates available experimental information across families of related genes into models that reconstruct the gain and loss of functional characteristics over evolutionary time. The models and the resulting set of 68,667 integrated gene functions cover approximately 82% of human protein-coding genes. The functional repertoire reveals a marked preponderance of molecular regulatory functions, and the models provide insights into the evolutionary origins of human gene functions. We show that our set of descriptions of functions can improve the widely used genomic technique of Gene Ontology enrichment analysis. The experimental evidence for each functional characteristic is recorded, thereby enabling the scientific community to help review and improve the resource, which we have made publicly available.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Linear-depth quantum circuits for loading Fourier approximations of arbitrary functions

Abstract The ability to efficiently load functions on quantum computers with high fidelity is essential for many quantum algorithms, including those for solving partial differential equations and Monte Carlo estimation. In this work, we introduce the Fourier series loader (FSL) method for preparing quantum states that exactly encode multi-dimensional Fourier series using linear-depth quantum circuits. Specifically, the FSL method prepares a (Dn)-qubit state encoding the 2 Dn -point uniform discretization of aD-dimensional function specified by aD-dimensional Fourier series. A free parameter,m, which must be less thann, determines the number of Fourier coefficients, 2 D ( m + 1 ) , used to represent the function. The FSL method uses a quantum circuit of depth at most 2 ( n − 2 ) + ⌈ log 2 ( n − m ) ⌉ + 2 D ( m + 1 ) + 2 − 2 D ( m + 1 ) , which is linear in the number of Fourier coefficients, and linear in the number of qubits (Dn) despite the fact that the loaded function’s discretization is over exponentially many (2 Dn ) points. The FSL circuit consists of at most D n + 2 D ( m + 1 ) + 1 − 1 single-qubit and D n ( n + 1 ) / 2 + 2 D ( m + 1 ) + 1 − 3 D ( m + 1 ) − 2 two-qubit gates; we present a classical compilation algorithm with runtime O ( 2 3 D ( m + 1 ) ) to determine the FSL circuit for a given Fourier series. The FSL method allows for the highly accurate loading of complex-valued functions that are well-approximated by a Fourier series with finitely many terms. We report results from noiseless quantum circuit simulations, illustrating the capability of the FSL method to load various continuous 1D functions, and a discontinuous 1D function, on 20 qubits with infidelities of less than 10 −6 and 10 −3 , respectively. We also demonstrate the practicality of the FSL method for near-term quantum computers by presenting experiments performed on the Quantinuum H1-1 and H1-2 trapped-ion quantum computers: we loaded a complex-valued function on 3 qubits with a fidelity of over 95 % , as well as various 1D real-valued functions on up to 6 qubits with classical fidelities ≈99%, and a 2D function on 10 qubits with a classical fidelity ≈94%.

Physics↗
DOE OSTI2023

An expectation–maximization framework for comprehensive prediction of isoform-specific functions

Advances in RNA sequencing technologies have achieved an unprecedented accuracy in the quantification of mRNA isoforms, but our knowledge of isoform-specific functions has lagged behind. There is a need to understand the functional consequences of differential splicing, which could be supported by the generation of accurate and comprehensive isoform-specific gene ontology annotations. We present isoform interpretation, a method that uses expectation–maximization to infer isoform-specific functions based on the relationship between sequence and functional isoform similarity. We predicted isoform-specific functional annotations for 85 617 isoforms of 17 900 protein-coding human genes spanning a range of 17 430 distinct gene ontology terms. Comparison with a gold-standard corpus of manually annotated human isoform functions showed that isoform interpretation significantly outperforms state-of-the-art competing methods. We provide experimental evidence that functionally related isoforms predicted by isoform interpretation show a higher degree of domain sharing and expression correlation than functionally related genes. We also show that isoform sequence similarity correlates better with inferred isoform function than with gene-level function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Functional Redundancy in Soil Microbial Community Based on Metagenomics Across the Globe

Understanding the contribution of soil microbial communities to ecosystem processes is critical for predicting terrestrial ecosystem feedbacks under changing climate. Our current understanding lacks a consistent strategy to formulate the linkage between microbial systems and ecosystem processes due to the presumption of functional redundancy in soil microbes. Here we present a global soil microbial metagenomic analysis to generalize patterns of microbial taxonomic compositions and functional potentials across climate and geochemical gradient. Our analyses show that soil microbial taxonomic composition varies widely in response to climate and soil physicochemical gradients, while microbial functional attributes based on metagenomic gene abundances are redundant. Among 17 climate zones, microbial taxonomic compositions were more distinct than functional potentials, as climate and edaphic properties showed more significant influence on microbial taxonomic compositions than on functional potentials. Microbial taxonomies formed a larger and more complex co-occurrence network with more module structures than functional potentials. Functional network was strongly inter-connected among different categories, whereas taxonomic network was more positively interactive in the same taxonomic groups. This study provides strong evidence to support the hypothesis of functional redundancy in soil microbes, as microbial taxonomic compositions vary to a larger extent than functional potentials based on metagenomic gene abundances in terrestrial ecosystems across the globe.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Self‐Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces

The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as “Jacob's ladder.” This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBEα framework with mathematical equation being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.

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

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Differential structure and functional gene response to geochemistry associated with the suspended and attached shallow aquifer microbiomes from the Illinois Basin, IL

Despite the clear ecological significance of the microbiomes inhabiting groundwater and connected ecosystems, our current understanding of their habitats, functionality, and the ecological processes controlling their assembly have been limited. In this study, an efficient pipeline combining geochemistry, high-throughput Fluidigm TM functional gene amplification and sequencing was developed to analyze the suspended and attached microbial communities inhabiting five groundwater monitoring wells in the Illinois Basin, USA. The dominant taxa in the suspended and the attached microbial communities exhibited significantly different spatial and temporal changes in both alpha- and beta-diversity. Further analyses of representative functional genes affiliated with N 2 fixation (nifH), methane oxidation (pmoA), and sulfate reduction (dsrB, and aprA), suggested functional redundancy within the shallow aquifer microbiomes. While more diversified functional gene taxa were observed for the suspended microbial communities than the attached ones except for pmoA, different levels of changes over time and space were observed between these functional genes. Notably, deterministic and stochastic ecological processes shaped the assembly of microbial communities and functional gene reservoirs differently. While homogenous selection was the prevailing process controlling assembly of microbial communities, the neutral processes (e.g., dispersal limitation, drift and others) were more important for the functional genes. The results suggest complex and changing shallow aquifer microbiomes, whose functionality and assembly vary even between the spatially proximate habitats and fractions. As a result, this research underscored the importance to include all the interface components for a more holistic understanding of the biogeochemical processes in aquifer ecosystems, which is also instructive for practical applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗