Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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 OSTI2025

Optical Activity of Chiral Excitons

Recent activity in the area of chiroptical phenomena has been focused on the connection between structural asymmetry, electron spin configuration and light/matter interactions in chiral semiconductors. In these systems, spin‐splitting phenomena emerge due to inversion symmetry breaking and the presence of extended electronic states, yet the connection to chiroptical phenomena is lacking. Here, in this study, we develop an analytical effective mass model of chiral excitons, parameterized by density functional theory. The model accounts for parity mixing of the band edge Bloch functions resulting from polar distortions, resulting in allowed magnetic dipole transitions. Through the study of a prototypical chiral 2D hybrid perovskite semiconductor, we show that circular dichroism of the chiral exciton and its interband continuum emerges from spin‐splitting via cross‐coupling of Rashba‐like and chiral/helical spin‐texture components. To demonstrate the generality of our approach, and as a counterpoint, we apply our model to describe chiroptical properties of three‐dimensionally confined excitons in perovskite nanocrystals that occur without chiral lattice distortions.

14 SOLAR ENERGY↗
DOE OSTI2023

Development of an interatomic potential for the W–Ta system

A physics-inspired, data-driven interatomic potential framework for multi-element system is presented. This potential is based on the generalization of the embedded atom method potential and systematically incorporates two-, three- and many-body effects. Different atomic environments are described by atom-centered Gaussian basis sets which provides sufficient flexibility to capture diverse atomic environments and allows for easy optimization of the free parameters. An interatomic potential model for the tungsten–tantalum (W–Ta) system is developed using this framework by training on data from ab initio density functional theory (DFT) calculations. The potential is thoroughly tested at various compositions by comparing bulk and defect properties from experimental and DFT data. It is shown that the potential predicts the elastic constants, defect properties, such as vacancy and interstitial formation energies, core structures of dislocations, and melting points for both pure and tantalum and tungsten with an accuracy comparable to other single element machine learning based potentials. In conclusion, the potential model is used to investigate the formation energies of ordered alloys and vacancy formation energies for chemically random W–Ta alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

In Situ Phosphatizing of Triphenylphosphine Encapsulated within Metal–Organic Frameworks to Design Atomic Co 1 –P 1 N 3 Interfacial Structure for Promoting Catalytic Performance

The engineering coordination environment offers great opportunity in performance tunability of isolated metal single-atom catalysts. For the most popular metal–N x (MN x ) structure, the replacement of N atoms by some other atoms with relatively weak electronegativity has been regarded as a promising strategy for optimizing the coordination environment of an active metal center and promoting its catalytic performance, which is still a challenge. Herein, we proposed a new synthetic strategy of an in situ phosphatizing of triphenylphosphine encapsulated within metal–organic frameworks for designing atomic Co 1 –P 1 N 3 interfacial structure, where a cobalt single atom is costabilized by one P atom and three N atoms (denoted as Co-SA/P-in situ). In the acidic media, the Co-SA/P-in situ catalyst with Co 1 –P 1 N 3 interfacial structure exhibits excellent activity and durability for the hydrogen evolution reaction (HER) with a low overpotential of 98 mV at 10 mA cm –2 and a small Tafel slope of 47 mV dec –1 , which are greatly superior to those of catalyst with Co 1 –N 4 interfacial structure. We discover that the bond-length-extended high-valence Co 1 –P 1 N 3 atomic interface structure plays a crucial role in boosting the HER performance, which is supported by in situ X-ray absorption fine structure (XAFS) measurements and density functional theory (DFT) calculation. With this work, we hope to promote the development of high performance metal single-atom catalysts.

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

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗
DOE OSTI2022

Modeling Activities Related to Waste Form Degradation: Progress Report

This report represents the milestone deliverable M4SF-22SN010309092 “Modeling Activities Related to Waste Form Degradation: Progress Report” that describes the progress of R&D activities of ongoing modeling investigations specifically on nuclear waste glass degradation, Density Functional Theory (DFT) studies on clarkeite structure and stability, and electrochemical model development of spent nuclear fuel (SNF). These activities are part of the Waste Form Testing, Modeling, and Performance work package at Sandia National Laboratories (SNL). This work package is part of the “Inventory and Waste Form Characteristics and Performance” control account that includes various experimental and modeling activities on nuclear waste degradation conducted at Oak Ridge National Laboratory (ORNL), SNL, Argonne National Laboratory (ANL), and Pacific Northwest National Laboratory (PNNL).

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

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Engineering Solid Electrolyte Interphase Composition by Assessing Decomposition Pathways of Fluorinated Organic Solvents in Lithium Metal Batteries

Studies have shown fluorinated electrolyte solvents can form desirable solid electrolyte interphase (SEI) in lithium metal batteries. In this study, we develop a detailed mechanistic understanding of two high performing electrolytes, Fluoroethylene Carbonate (FEC) and Difluoroethylene Carbonate (DFEC) to demonstrate minimal structural variations can lead to different decomposition products, and thereby the nature of the SEI. Using density functional theory (DFT) calculations, we find different initial bond-breaking mechanisms between FEC and DFEC. We develop free energy diagrams for the decomposition pathways including both electrochemical and chemical steps. Using the computational Li electrode, we identify the largest limiting potential of 1.77 V for FEC decomposition, associated with the formation of lithium fluoride, lithium oxide and FEC oligomers, and 1.53 V for DFEC, which correspond to the formation of polymerized vinylene carbonate and lithium fluoride. We suggest the formation of oligomers in the case of FEC instead of long polymers may lead to better SEI compactness. We also demonstrate the SEI components of FEC and DFEC are not stable on typical cathode voltage (3.87 V). This study presents a unified electrocatalytic perspective on SEI formation and decomposition.

25 ENERGY STORAGE↗
DOE OSTI2024

Leveraging Curvature on N–Doped Carbon Materials for Hydrogen Storage

Carbon sorbent materials have shown great promise for solid-state hydrogen (H 2 ) storage. Modification of these materials with nitrogen (N) dopants has been undertaken to develop materials that can store H 2 at ambient temperatures. In this work density functional theory (DFT) calculations are used to systematically probe the influence of curvature on the stability and activity of undoped and N-doped carbon materials toward H binding. Specifically, four models of carbon materials are used: graphene, [5,5] carbon nanotube, [5,5] D 5d -C 120 , and C 60 , to extract and correlate the thermodynamic properties of active sites with varying degrees of sp 2 hybridization (curvature). From the calculations and analysis, it is found that graphitic N-doping is thermodynamically favored on more pyramidal sites with increased curvature. In contrast, it is found that the hydrogen binding energy is weakly affected by curvature and is dominated by electronic effects induced by N-doping. These findings highlight the importance of modulating the heteroatom doping configuration and the lattice topology when developing materials for H 2 storage.

08 HYDROGEN↗
DOE OSTI2023

Molecular Dynamics Simulation of Hygroscopic Aging Effects in Epoxy Polymer

The automobile industry is incorporating more lightweight content in car designs to boost fuel-economy. New structural adhesives are needed to mitigate the corrosion and thermal expansion issues associated with joining dissimilar lightweight materials, but adhesive developers lack a fundamental understanding of the chemistry that occurs in the adhesive as the joint ages. In this study, we developed structural adhesive molecular models and applied classical molecular dynamics simulations and density functional theory calculations to gain molecular insights into the influence of water molecules on the properties of epoxy-based adhesives (DGEBA + Jeffamine (JD230)). The simulations were complemented by experimental synthesis and characterization. Our work underscores the impact of water molecules on the local structure of the epoxy network as well as resulting mechanical properties. Water molecules were mainly coordinated with hydroxyls, primary amines and secondary amines, but also weakly coordinated with ether linkages, which were found most probable to be labile. Simulated stress–strain data indicates that increasing the water content deteriorates the mechanical properties. The Young’s modulus decreased by ~ 30% when the water content increased to 3 wt%. We conclude, this integration of molecular-level chemical insights with mechanical property simulations of the hydrated epoxy system and experimental validation holds the promise to advance lightweight joint technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Catalytic production of δ-valerolactone (DVL) from biobased 2-hydroxytetrahydropyran (HTHP) – Combined experimental and modeling study

δ-Valerolactone (DVL) is a five-carbon (C5) cyclic ester that can undergo ring-opening polymerization to yield high-performance, biocompatible polyesters. But current market prices of C5 chemicals like DVL are very high due to poor availability of C5 feedstock in petroleum. Herein, we demonstrate a novel route to DVL synthesis via dehydrogenation of biomass-derived 2-hydroxytetrahydropyran (HTHP) over Cu/SiO 2 without the use of toxic reagents. Since HTHP exists in thermal equilibrium with 3,4-dihydropyran (DHP) via dehydration, and with 2,2’-oxybis(tetrahydropyran) and 5-(tetrahydropyran-2-yloxy)pentanal via acetalization, we have also determined the thermochemistry (ΔH rxn and ΔG rxn ) of each competing reaction using density functional theory (DFT) calculations at the M06-2X/cc-pVTZ level. Finally, by developing a kinetic model of all 8 reactions involved, we have achieved 84% selectivity to DVL at 150°C in a packed bed reactor for over 72 hours of time on stream.

42 ENGINEERING↗
DOE OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2024

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗
DOE OSTI2020

Symmetry-breaking polymorphous descriptions for correlated materials without interelectronic U

Correlated materials with open-shell d- and f-ions having degenerate band edge states show a rich variety of interesting properties ranging from metal-insulator transition to unconventional superconductivity. The textbook view for the electronic structure of these materials is that mean-field approaches are inappropriate, as the interelectronic interaction U is required to open a band gap between the occupied and unoccupied degenerate states while retaining symmetry. We show that the latter scenario often defining what Mott insulators are, is in fact not needed for the 3d binary oxides MnO, FeO, CoO, and NiO. The mean-field-like band theory can indeed lift such degeneracies in the binaries when nontrivial unit cell representations (polymorphous networks) are allowed to break symmetries, in conjunction with a recently developed non-empirical exchange and correlation density-functional without an on-site interelectronic interaction U. We explain how density-functional theory (DFT) in the polymorphous representation achieves band gap opening in correlated materials through a separate mechanism to the Mott-Hubbard approach. Here, we show the method predicts magnetic moments and gaps for the four binary monoxides in both the antiferromagnetic and paramagnetic phases, offering an effective alternative to symmetry-conserving approaches for studying a range of functionalities in open d- and f-shell complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Bespoke Liquid/Liquid Interfaces (Final Technical Report)

The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Computational Study of an Iron(II) Polypyridine Electrocatalyst for CO 2 Reduction: Key Roles for Intramolecular Interactions in CO 2 Binding and Proton Transfer

A solar-driven conversion of CO 2 into fuels by artificial photosynthesis would not only mitigate the greenhouse effect but also provide an alternative to obtain fuels in a renewable fashion. To this end, the new iron polypyridine catalyst [Fe(bpy NHEt PY2Me)L 2 ] 2+ (L = H 2 O, CH 3 CN) was recently developed for the electrochemical reduction of CO 2 to CO. Here, we performed density functional theory (DFT) electronic structure calculations to shed light on a possible pathway for CO 2 reduction and the origin of the selectivity between CO 2 reduction versus the hydrogen evolution reaction. The metal center remains Lewis acidic throughout the reduction process due to ligand loss and mainly ligand-based reduction stabilized by antiferromagnetic coupling to a high-spin Fe(II) center. This results in a high barrier for hydride formation but a facile addition and activation of CO 2 via an η 2 coordination and stabilizing hydrogen bonding by the amine group. The second unoccupied equatorial coordination site opens up the possibility for an intramolecular protonation with a coordinated water ligand. This facilitates protonation because not only CO 2 but also the proton source H 2 O is activated and properly aligned for a proton transfer due to the Fe-OH 2 bond; consequently, both protonation steps are facile. The moderate ligand field allows a rapid ligand exchange for a second intramolecular protonation step and facilitates an exergonic CO release. The lower selectivity of the related [Fe(bpy OH PY2Me)L 2 ] 2+ complex can be related to its more acidic second coordination sphere, which enables an intramolecular proton transfer that is kinetically competitive with CO 2 addition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗