Search NASA⌕ Search

SEARCH · Search NASA

Results for “Classical density functional theory”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Electronic Transport Through Carbon Nanotubes: Effects of Structural Deformation and the Tube Chirality

Atomistic simulations using a combination of classical force field and Density-Functional-Theory (DFT) show that carbon atoms remain essentially sp2 coordinated in either bent tubes or tubes pushed by an atomically sharp AFM tip. Subsequent Green's-function-based transport calculations reveal that for armchair tubes there is no significant drop in conductance, while for zigzag tubes the conductance can drop by several orders of magnitude in AFM-pushed tubes. The effect can be attributed to simple stretching of the tube under tip deformation, which opens up an energy gap at the Fermi surface.

Maiti, Amitesh↗
NASA NTRS2002

Effects of Structural Deformation and Tube Chirality on Electronic Conductance of Carbon Nanotubes

A combination of large scale classical force-field (UFF), density functional theory (DFT), and tight-binding Green's function transport calculations is used to study the electronic properties of carbon nanotubes under the twist, bending, and atomic force microscope (AFM)-tip deformation. We found that in agreement with experiment a significant change in electronic conductance can be induced by AFM-tip deformation of metallic zigzag tubes and by twist deformation of armchair tubes. The effect is explained in terms of bandstructure change under deformation.

Svizhenko, Alexei↗
DOE OSTI2023

Quantum simulation of exact electron dynamics can be more efficient than classical mean-field methods

Quantum algorithms for simulating electronic ground states are slower than popular classical mean-field algorithms such as Hartree–Fock and density functional theory but offer higher accuracy. Accordingly, quantum computers have been predominantly regarded as competitors to only the most accurate and costly classical methods for treating electron correlation. However, here we tighten bounds showing that certain first-quantized quantum algorithms enable exact time evolution of electronic systems with exponentially less space and polynomially fewer operations in basis set size than conventional real-time time-dependent Hartree–Fock and density functional theory. Although the need to sample observables in the quantum algorithm reduces the speedup, we show that one can estimate all elements of the k-particle reduced density matrix with a number of samples scaling only polylogarithmically in basis set size. We also introduce a more efficient quantum algorithm for first-quantized mean-field state preparation that is likely cheaper than the cost of time evolution. We conclude that quantum speedup is most pronounced for finite-temperature simulations and suggest several practically important electron dynamics problems with potential quantum advantage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2022

First-principles ionized-impurity scattering and charge transport in doped materials

Scattering of carriers with ionized impurities governs charge transport in doped semiconductors. However, electron interactions with ionized impurities cannot be fully described with quantitative first-principles calculations, so their understanding relies primarily on simplified models. Here we show an ab initio approach to compute the interactions between electrons and ionized impurities or other charged defects. It includes the short- and long-range electron-defect (e-d) interactions on equal footing and allows for efficient interpolation of the e-d matrix elements. Here we combine the e-d and electron-phonon interactions in the Boltzmann transport equation to compute the carrier mobilities in doped silicon over a wide range of temperature and doping concentrations, seamlessly spanning the defect- and phonon-limited transport regimes. The individual contributions of the defect- and phonon-scattering mechanisms to the carrier relaxation times and mean-free paths are analyzed. Our method provides a powerful tool to study electronic interactions in doped materials. It broadens the scope of first-principles transport calculations, enabling studies of a wide range of doped semiconductors and oxides with application to electronics, energy, and quantum technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Modeling of Battery Materials

This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.

25 ENERGY STORAGE↗
DOE OSTI2024

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

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

Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study

Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.

Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗
DOE OSTI2022

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl 2 from T m + 85 K to T m + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl 2 . We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. Furthermore, this result may contribute to new understandings of the structural relaxation of other network liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanism of Solid Electrolyte Interphase Formation on Lithium Anodes Studied by Molecular Dynamics Simulations

We investigated the initial formation of the SEI on lithium-metal anodes in presence of several electrolyte compositions using ReaxFF, in absence of a bias potential. We observed the formation of well-defined regions with different density and oxidation states of Li. Dissolution of Li metal atoms results in an uneven distribution of Li in the liquid phase. The electrolyte decomposition mechanism was studied for various electrolyte mixtures. Overall, the reaction mechanisms observed with ReaxFF are in good agreement with previous reports using DFT. However, we found some discrepancies in specific predictions by the ReaxFF that should be addressed in future work.

25 ENERGY STORAGE↗
DOE OSTI2022

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

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

Engineering local strain for single-atom nuclear acoustic resonance in silicon

Mechanical strain plays a key role in the physics and operation of nanoscale semiconductor systems, including quantum dots and single-dopant devices. In this report we describe the design of a nanoelectronic device, where a single nuclear spin is coherently controlled via nuclear acoustic resonance (NAR) through the local application of dynamical strain. The strain drives spin transitions by modulating the nuclear quadrupole interaction. We adopt an AlN piezoelectric actuator compatible with standard silicon metal–oxide–semiconductor processing and optimize the device layout to maximize the NAR drive. We predict NAR Rabi frequencies of order 200 Hz for a single 123Sb nucleus in a wide region of the device. Spin transitions driven directly by electric fields are suppressed in the center of the device, allowing the observation of pure NAR. Using electric field gradient-elastic tensors calculated by the density-functional theory, we extend our predictions to other high-spin group-V donors in silicon and to the isoelectronic 73Ge atom.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2024

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

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

Structure prediction of porous organic crystals

In this work, we explore the possibility of applying automated crystal structure prediction to reproduce the experimentally identified metastable porous polymorphs. Using our recently developed High-Throughput Organic Crystal Structure Prediction ( HTOCSP ) framework, we conducted a systematic study on five representative organic crystalline systems including hydrogen-bonded frameworks (HOFs), featured by the presence of significant porosity, in conjunction with different choices of energy models from classical, machine learning force fields, tight binding to density functional theory. Our results suggest that the current structure generation framework, with careful selection of symmetry conditions, is likely to generate rather complex and abundant metastable crystal candidates for porous crystals. In conjunction with the recent advance in universal machine learning force fields, it becomes possible to identify experimental structures as the energetically favorable candidates from a simple energy versus density analysis, thus paving the way for computational design of complex porous materials with the target systems prior to the experimental synthesis and characterization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Liquid to crystal Si growth simulation using machine learning force field

Machine learning force field (ML-FF) has emerged as a potential promising approach to simulate various material phenomena for large systems with ab initio accuracy. However, most ML-FFs have been used to study the phenomena relatively close to the equilibrium ground states. In this work, we have studied a far from equilibrium system of liquid to crystal Si growth using ML-FF. Here, we found that our ML-FF based on ab initio decomposed atomic energy can reproduce all the aspects of ab initio simulated growth, from local energy fluctuations to transition temperatures, to diffusion constant, and growth rates. We have also compared the growth simulation with the Stillinger-Weber classical force field and found significant differences. A procedure is also provided to correct a systematic fitting bias in the ML-FF training process, which exists in all training models, otherwise critical results like transition temperature will be wrong.

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

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗