Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT simulations”

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

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

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

Mechanistic insights into carbon–carbon coupling on NiAu and PdAu single-atom alloys

Carbon–carbon coupling is an important step in many catalytic reactions, and performing sp 3 –sp 3 carbon–carbon coupling heterogeneously is particularly challenging. It has been reported that PdAu single-atom alloy (SAA) model catalytic surfaces are able to selectively couple methyl groups, producing ethane from methyl iodide. In this work, we extend this study to NiAu SAAs and find that Ni atoms in Au are active for C–I cleavage and selective sp 3 –sp 3 carbon–carbon coupling to produce ethane. Furthermore, we perform ab initio kinetic Monte Carlo simulations that include the effect of the iodine atom, which was previously considered a bystander species. We find that model NiAu surfaces exhibit a similar chemistry to PdAu, but the reason for the similarity is due to the role the iodine atoms play in terms of blocking the Ni atom active sites. Specifically, on NiAu SAAs, the iodine atoms outcompete the methyl groups for occupancy of the Ni sites leaving the Me groups on Au, while on PdAu SAAs, the binding strengths of methyl groups and iodine atoms at the Pd atom active site are more similar. These simulations shed light on the mechanism of this important sp 3 –sp 3 carbon–carbon coupling chemistry on SAAs. Furthermore, we discuss the effect of the iodine atoms on the reaction energetics and make an analogy between the effect of iodine as an active site blocker on this model heterogeneous catalyst and homogeneous catalysts in which ligands must detach in order for the active site to be accessed by the reactants.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Nonadiabatic couplings from a variational excited state method based on constrained DFT

Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

An OFDM System Using Polyphase Filter and DFT Architecture for Very High Data Rate Applications

This paper presents a conceptual architectural design of a four-channel Orthogonal Frequency Division Multiplexing (OFDM) system with an aggregate information throughput of 622 megabits per second (Mbps). Primary emphasis is placed on the generation and detection of the composite waveform using polyphase filter and Discrete Fourier Transform (DFT) approaches to digitally stack and bandlimit the individual carriers. The four-channel approach enables the implementation of a system that can be both power and bandwidth efficient, yet enough parallelism exists to meet higher data rate goals. It also enables a DC power efficient transmitter that is suitable for on-board satellite systems, and a moderately complex receiver that is suitable for low-cost ground terminals. The major advantage of the system as compared to a single channel system is lower complexity and DC power consumption. This is because the highest sample rate is half that of the single channel system and synchronization can occur at most, depending on the synchronization technique, a quarter of the rate of a single channel system. The major disadvantage is the increased peak-to-average power ratio over the single channel system. Simulation results in a form of bit-error-rate (BER) curves are presented in this paper.

Kifle, Muli↗
DOE OSTI2025

Recovering high-purity uranyl nitrate from simulated used nuclear fuel dissolver solutions by crystallization: rejecting technetium

The separation of U from Tc and other problematic fission product elements like Mo and Ru, along with Sr, Zr, Cs, and Nd, has been achieved via the crystallization of uranyl nitrate hexahydrate (UNH). Rejection of technetium as pertechnetate anion ( 99 TcO 4 – ) is an especially important feature of this system, as it otherwise tends to follow U (VI) in extractive separations. It also raises the salient question regarding why this oxoanion cannot replace nitrate within the crystalline lattice of UNH. Results showed high-yield (>90 %), high-purity (>99 %) recovery of U as UNH from solutions containing 99 TcO 4 – by simple reduction of temperature from 60°C to 20°C. There was no observable interaction of 99 TcO 4 – with UO 2 2+ . The addition of other cations like, Sr 2+ , Zr 4+ , Cs + , and Nd 3+ , also did not form secondary, contaminant solid phases, leaving the > 99 % of the fission product elements in the mother liquor, while the U was recovered at > 90 %. Similarly, Mo and Ru, when added to the mixture, were shown to behave as the other fission-product elements, remaining in the mother liquor during crystallization. As a result, DFT calculations showed that, despite the higher binding strength of TcO 4 – , HMoO 4 – , and BiO 3 – with the UO 2 2+ cation compared to NO 3 – , the hydrogen-bonding network of the two coordinated ions and four waters of hydration in the UNH crystal structure is the driving force for the high specificity of this separation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Introducing the embedded random phase approximation: H 2 dissociative adsorption on Cu(111) as an exemplar

The random phase approximation (RPA) as a means of treating electron correlation recently has been shown to outperform standard density functional theory (DFT) approximations in a variety of cases. However, the computational cost of the RPA is substantially more than DFT, especially when aiming to study extended surfaces. Properly accounting for sufficient surface ensemble size, Brillouin zone sampling, and vacuum separation of periodic images in standard periodic-planewave-based DFT code raises the cost to achieve converged results. Here, we show that sub-system embedding schemes enable use of the RPA for modeling heterogeneous reactions at reduced computational cost. Further, we explore two different embedded RPA (emb-RPA) approaches, periodic emb-RPA and cluster emb-RPA. We use the (experimentally and theoretically) well-studied H 2 dissociative adsorption on Cu(111) as our exemplar, and first perform full periodic RPA calculations as a benchmark. The full RPA results match well the semi-empirical barrier fit to experimental observables and others derived from high-level computations, e.g., from recent embedded n-electron valence second order perturbation theory [Zhao et al., J. Chem. Theory Comput. 16(11), 7078–7088 (2020)] and quantum Monte Carlo [Doblhoff-Dier et al., J. Chem. Theory Comput. 13(7), 3208–3219 (2017)] simulations. Among the two emb-RPA approaches tested, the cluster emb-RPA accurately reproduces the energy profile (maximum error of 50 meV along the reaction pathway) while reducing the computational cost by approximately two orders of magnitude. We therefore expect that the embedded cluster approach will enable wider RPA implementation in heterogeneous catalysis.

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

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗
DOE OSTI2021

Pseudoelastic deformation in Mo-based refractory multi-principal element alloys

Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal-element alloys. We present phase and property analysis of quinary (MoW) x Zr y (TaTi) 1-x-y refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo-W-Ta-Ti-Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress-strain. Our analysis shows that only Mo-W-rich compositions of Mo-W-Ta-Ti-Zr, i.e., Mo+W, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW) 85 Zr 7.5 (TaTi) 7.5 was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning-assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. Furthermore, this alloy shows identical cyclic deformation characteristics during uniaxial < 100 > loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1,500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase Diagrams of Alloys and Their Hydrides via On-Lattice Graph Neural Networks and Limited Training Data

Efficient prediction of sampling-intensive thermodynamic properties is needed to evaluate material performance and permit high-throughput materials modeling for a diverse array of technology applications. To alleviate the prohibitive computational expense of high-throughput configurational sampling with density functional theory (DFT), surrogate modeling strategies like cluster expansion are many orders of magnitude more efficient but can be difficult to construct in systems with high compositional complexity. We therefore employ minimal-complexity graph neural network models that accurately predict and can even extrapolate to out-of-train distribution formation energies of DFT-relaxed structures from an ideal (unrelaxed) crystallographic representation. This enables the large-scale sampling necessary for various thermodynamic property predictions that may otherwise be intractable and can be achieved with small training data sets. Two exemplars, optimizing the thermodynamic stability of low-density high-entropy alloys and modulating the plateau pressure of hydrogen in metal alloys, demonstrate the power of this approach, which can be extended to a variety of materials discovery and modeling problems.

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

Structure Prediction of Ionic Epitaxial Interfaces with Ogre Demonstrated for Colloidal Heterostructures of Lead Halide Perovskites

Colloidal epitaxial heterostructures are nanoparticles composed of two different materials connected at an interface, which can exhibit properties different from those of their individual components. Combining dissimilar materials offers exciting opportunities to create a wide variety of functional heterostructures. However, assessing structural compatibility–the main prerequisite for epitaxial growth–is challenging when pairing complex materials with different lattice parameters and crystal structures. This complicates both the selection of target heterostructures for synthesis and the assignment of interface models when new heterostructures are obtained. Here, we demonstrate Ogre as a powerful tool to accelerate the design and characterization of colloidal heterostructures. To this end, we implemented developments tailored for the high-efficiency prediction of epitaxial interfaces between ionic/polar materials, which encompass most colloidal semiconductors. These include the use of pre-screening candidate models based on charge balance at the interface and the use of a classical potential for fast energy evaluations, with parameters automatically calculated based on the input bulk structures. These developments are validated for perovskite-based CsPbBr 3 /Pb 4 S 3 Br 2 heterostructures, where Ogre produces interface models in excellent agreement with density functional theory and experiments. Furthermore, we use Ogre to rationalize the templating effect of CsPbCl 3 on the growth of lead sulfochlorides, where perovskite seeds induce the formation of Pb 4 S 3 Cl 2 rather than Pb 3 S 2 Cl 2 due to better epitaxial compatibility. Finally, combining Ogre simulations with experimental data enables us to unravel the structure and composition of the hitherto unsolved CsPbBr 3 /Bi x Pb y S z interface, and to assign a structure to several other reported metal halide- and oxide-based interfaces. The Ogre package is available on GitHub or via the OgreInterface desktop application, available for Windows, Linux, and Mac.

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

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗
NASA NTRS2014

Modeling Materials: Design for Planetary Entry, Electric Aircraft, and Beyond

NASA missions push the limits of what is possible. The development of high-performance materials must keep pace with the agency's demanding, cutting-edge applications. Researchers at NASA's Ames Research Center are performing multiscale computational modeling to accelerate development times and further the design of next-generation aerospace materials. Multiscale modeling combines several computationally intensive techniques ranging from the atomic level to the macroscale, passing output from one level as input to the next level. These methods are applicable to a wide variety of materials systems. For example: (a) Ultra-high-temperature ceramics for hypersonic aircraft-we utilized the full range of multiscale modeling to characterize thermal protection materials for faster, safer air- and spacecraft, (b) Planetary entry heat shields for space vehicles-we computed thermal and mechanical properties of ablative composites by combining several methods, from atomistic simulations to macroscale computations, (c) Advanced batteries for electric aircraft-we performed large-scale molecular dynamics simulations of advanced electrolytes for ultra-high-energy capacity batteries to enable long-distance electric aircraft service; and (d) Shape-memory alloys for high-efficiency aircraft-we used high-fidelity electronic structure calculations to determine phase diagrams in shape-memory transformations. Advances in high-performance computing have been critical to the development of multiscale materials modeling. We used nearly one million processor hours on NASA's Pleiades supercomputer to characterize electrolytes with a fidelity that would be otherwise impossible. For this and other projects, Pleiades enables us to push the physics and accuracy of our calculations to new levels.

Supercomputing↗
DOE OSTI2026

Oxygen Vacancy Evolution at Li x V 2 O 5 /LiPON Solid State Electrochemical Interfaces Using Depth Resolved Cathodoluminescence Spectroscopy

The formation of oxygen vacancies at buried LiPON/ Li x V 2 O 5 interfaces has been observed on a near-nanometer scale and nondestructively using depth-resolved cathodoluminescence spectroscopy (DRCLS) and interfacial markers. Before electrochemical cycling, as-deposited LiPON/Li x V 2 O 5 exhibits a 1.6 eV defect optical emission, which density functional theory calculations identify as originating from oxygen vacancies. This defect appears first within a few nanometers of the buried LiPON/Li x V 2 O 5 interface without cycling, indicating that spontaneous O diffusion from the Li x V 2 O 5 lattice into LiPON may have caused these interface-localized oxygen vacancy defects. DRCLS measured the intensity and spatial distribution of this oxygen vacancy signal as a function of electrochemical cycling in a LiPON/Li x V 2 O 5 half-cell, showing oxygen vacancy signal increasing and moving deeper into the electrode with increased cycle number. Significant electrochemical irreversibility was also observed, with poor Coulombic efficiency and a 15% drop in capacity over 50 cycles. Theoretical simulations predict that the presence of oxygen vacancies increases the energy barrier for lithium diffusion significantly, indicating that this aggregation of oxygen vacancies could be another battery degradation mechanism accompanying lithiation induced phase changes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

New Dimensions in the Theory of Excited States and X-ray Spectra (Final Report)

This Final Technical Report briefly summarizes the achievements during the lifetime of our DOE BES grant DE-FG02-97ER45623. The long-term goal of this project has been the development of quantitative theories of the interaction between radiation and matter, with a focus on x-ray spectroscopies. X-ray spectra have long been among the most important probes of atomic-scale structure and properties of matter, ranging from atoms and molecular systems to condensed matter and exotic states. These spectroscopies are widely used in investigations at the major DOE synchrotron x-ray facilities and related centers world-wide. In addition to fundamental theory, a major goal of our project has been the development of computational software that implements the theory for calculations of x-ray spectra of various materials throughout the periodic table. Due to the complex nature of x-ray spectra, quantitative theory is essential for its interpretation. The theory is challenging since it involves excited state electronic structure and many-body correlation effects that go beyond independent particle approximations like DFT or Hartree-Fock. Moreover, the experimental investigations typically involve a broad range of energy, time, and temperature scales, from the UV-Vis to hard x-ray energies of order 10 4 eV, and temperatures T from ambient up to the warm-dense-matter regime where the Fermi energy kBTF is of order a few eV, i.e., temperatures of order 105 Kelvin. This broad range of experimental conditions has fostered many novel theoretical approaches and computational techniques, many of which we have developed systematically over the duration of the grant. In contrast to the traditional wave-function approach of quantum theory and electronic structure methods, our theoretical approach is based on modern Green's function techniques. This approach is better suited for aperiodic structures, excited states, and broad spectral ranges, since it avoids the computational bottlenecks of sum-over-states approaches, as in the Fermi golden rule. This theoretical framework has been incorporated into efficient, user-friendly x-ray spectroscopy software which is now used routinely worldwide to simulate and analyze spectra. These theoretical tools provide an essential complement to synchrotron and next-generation light sources, which are used to investigate complex materials with ever increasing precision. Moreover, the synergism between theory, computation and experiment contributed by our research enhances scientific understanding and creates opportunities for innovations in materials and energy science and in many fields. As documented in this Report, this research grant has been remarkably successful in achieving these goals. In particular, this grant has supported the development of the x-ray spectroscopy software suite known as FEFF (named for an effective scattering amplitude f eff in the theory). The FEFF codes have become one of the premier tools for quantitative simulations of x-ray spectra as documented by many thousands of citations in the Web of Science and Google-Scholar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Local cation ordering in compositionally complex Ruddlesden–Popper n = 1 oxides

The Ruddlesden–Popper (RP) layered perovskite structure is of great interest due to its inherent tunability, and the emergence and growth of the compositionally complex oxide (CCO) concept endows the RP family with further possibilities. Here, a comprehensive assessment of thermodynamic stabilization, local order/disorder, and lattice distortion was performed in the first two reported examples of lanthanum-deficient La n+1 B n O 3n+1 (n = 1, B = Mg, Co, Ni, Cu, Zn) obtained via various processing conditions. Chemical short-range order (CSRO) at the B-site and the controllable excess interstitial oxygen (δ) in RP-CCOs are uncovered by neutron pair distribution function analysis. Reverse Monte Carlo analysis of the data, Metropolis Monte Carlo simulations, and extended x-ray absorption fine structure analysis implies a modest degree of magnetic element segregation on the local scale. Further, ab initio molecular dynamics simulations results obtained from special quasirandom structure disagree with experimentally observed CSRO but confirm Jahn–Teller distortion of CuO 6 octahedra. These findings highlight potential opportunities to control local order/disorder and excess interstitial oxygen in layered RP-CCOs and demonstrate a high degree of freedom for tailoring application-specific properties. They also suggest a need for expansion of theoretical and data modeling approaches in order to meet the innate challenges of CCO and related high-entropy phases.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Computational Study of Field Initiated Surface Reactions for Synthesis of Diamond and Silicon

This project involves using quantum chemistry to simulate surface chemical reactions in the presence of an electric field for nanofabrication of diamond and silicon. A field delivered by a scanning tunneling microscope (STM) to a nanometer scale region of a surface affects chemical reaction potential energy surfaces (PES) to direct atomic scale surface modification to fabricate sub-nanometer structures. Our original hypothesis is that the applied voltage polarizes the charge distribution of the valence electrons and that these distorted molecular orbitals can be manipulated with the STM so as to change the relative stabilities of the electronic configurations over the reaction coordinates and thus the topology of the PES and reaction kinetics. Our objective is to investigate the effect of applied bias on surface reactions and the extent to which STM delivered fields can be used to direct surface chemical reactions on an atomic scale on diamond and silicon. To analyze the fundamentals of field induced chemistry and to investigate the application of this technique for the fabrication of nanostructures, we have employed methods capable of accurately describing molecular electronic structure. The methods we employ are density functional theory (DFT) quantum chemical (QC) methods. To determine the effect of applied bias on surface reactions we have calculated the QC PESs in various applied external fields for various reaction steps for depositing or etching diamond and silicon. We have chosen reactions which are thought to play a role in etching and the chemical vapor deposition growth of Si and diamond. The PESs of the elementary reaction steps involved are then calculated under the applied fields, which we vary in magnitude and configuration. We pay special attention to the change in the reaction barriers, and transition state locations, and search for low energy reaction channels which were inaccessible without the applied bias.

Musgrave, Charles Bruce↗
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 OSTI2021

Building Toward the Future in Chemical and Materials Simulation with Accessible and Intelligently Designed Web Applications

Over the last few decades, significant progress has been made in the development and use of electronic structure and other molecular simulation methods. As these methods become more mature and are able to simulate larger and more complex chemical simulations, the need for improvement in scientific visualization, molecular builders, simplified input to simulation methods, and the development of new approaches and languages to describe simulations, along with workflows to carry them out, becomes more apparent. In this chapter, we describe our recent efforts in developing a prototype open-source computational tool called Arrows that combines NWChem, SQL and NoSQL databases, email, web APIs, and web applications in a way that make molecular and materials modeling accessible to all scientists and engineers. At the same time, because of its simplified input, it provides a framework for expert users to carry out large numbers of calculations and run complex workflows.

2↗