Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparing the Expense and Accuracy of Methods to Simulate Atomic Vibrations in Rubrene

Atomic vibrations can inform about materials properties from hole transport in organic semiconductors to correlated disorder in metal–organic frameworks. Currently, there are several methods for predicting these vibrations using simulations, but the accuracy–efficiency tradeoffs have not been examined in depth. In this work, rubrene is used as a model system to predict atomic vibrational properties using six different simulation methods: density functional theory, density functional tight binding, density functional tight binding with a Chebyshev polynomial-based correction, a trained machine learning model, a pretrained machine learning model called ANI-1, and a classical forcefield model. The accuracy of each method is evaluated by comparison to the experimental inelastic neutron scattering spectrum. All methods discussed here show some accuracy across a wide energy region, though the Chebyshev-corrected tight-binding method showed the optimal combination of high accuracy with low expense. We then offer broad simulation guidelines to yield efficient, accurate results for inelastic neutron scattering spectrum prediction.

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

Accurate Prediction of Adiabatic Ionization Potentials of Organic Molecules using Quantum Chemistry Assisted Machine Learning

In previous work (Dandu et al., J. Phys. Chem. A, 2022, 126, 4528–4536), we were successful in predicting accurate atomization energies of organic molecules using machine learning (ML) models, obtaining an accuracy as low as 0.1 kcal/mol compared to the G4MP2 method. In this work, we extend the use of these ML models to adiabatic ionization potentials on data sets of energies generated using quantum chemical calculations. Atomic specific corrections that were found to improve atomization energies from quantum chemical calculations have also been used in this study to improve ionization potentials. Here, the quantum chemical calculations were performed on 3405 molecules containing eight or fewer non-hydrogen atoms derived from the QM9 data set, using the B3LYP functional with the 6–31G(2df,p) basis set for optimization. Low-fidelity IPs for these structures were obtained using two density functional methods: B3LYP/6–31+G(2df,p) and ωB97XD/6–311+G(3df,2p). Highly accurate G4MP2 calculations were performed on these optimized structures to obtain high-fidelity IPs to use in ML models based on the low-fidelity IPs. Our best performing ML methods gave IPs of organic molecules within a mean absolute deviation of 0.035 eV from the G4MP2 IPs for the whole data set. This work demonstrates that ML predictions assisted by quantum chemical calculations can be used to successfully predict IPs of organic molecules for use in high throughput screening.

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

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

Machine learning to accelerate screening for Marcus reorganization energies

Understanding and predicting the charge transport properties of π-conjugated materials is an important challenge for designing new organic electronic devices, such as solar cells, plastic transistors, light-emitting devices, and chemical sensors. A key component of the hopping mechanism of charge transfer in these materials is the Marcus reorganization energy which serves as an activation barrier to hole or electron transfer. While modern density functional methods have proven to accurately predict trends in intramolecular reorganization energy, such calculations are computationally expensive. In this work, we outline active machine learning methods to predict computed intramolecular reorganization energies of a wide range of polythiophenes and their use toward screening new compounds with low internal reorganization energies. Our models have an overall root mean square error (RMSE) of ±0.113 eV, but a much smaller RMSE of only ±0.036 eV on the new screening set. Since the larger error derives from high-reorganization energy compounds, the new method is highly effective to screen for compounds with potentially efficient charge transport parameters.

14 SOLAR ENERGY↗
DOE OSTI2021

Octahedral rotations in Ruddlesden-Popper layered oxides under pressure from first principles

The combination of reduced dimensionality and tunable structural distortions in layered perovskite oxides makes these materials ideal platforms for designing novel properties and functionalities. One example is hybrid improper ferroelectricity in n = 2 Ruddlesden-Popper oxides, where the combination of a layered crystal structure and rotations of the metal-oxide octahedra break symmetry and induce a polarization. Precisely controlling the octahedral rotation distortions, for example by the application of hydrostatic pressure, provides a pathway to tune and optimize the properties of these materials. Here we combine group theoretic methods, density functional theory calculations, and Landau theory analysis to investigate how octahedral rotations respond to pressure in the hybrid improper ferroelectrics Sr 3 Zr 2 O 7 , Ca 3 Ti 2 O 7 , and Sr 3 Sn 2 O 7 . We find that factors that are known to control the pressure response of ABO 3 perovskites—the formal charge of the A- and B-site cations, tolerance factor, and B-site chemistry—also impact the pressure response of these layered perovskites. We also show that coupling between the octahedral rotation and strain order parameters plays a key role in determining the overall pressure response. Despite some similarities, we find that these layered perovskites display a distinct pressure response compared to their ABO 3 perovskite analogs. By identifying trends and underlying mechanisms that control octahedral rotations in Ruddlesden-Popper oxides under pressure, this work lays the foundation for tailoring the structure and properties of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2025

Electronic and thermal properties of the phase-change memory material, Ge 2 Sb 2 Te 5 , and results from spatially resolved transport calculations

Here, we report new insights into the electronic, structural, and transport (heat and charge) properties of the phase-change memory material amorphous Ge 2 Sb 2 Te 5 . Using realistic structural models of Konstantinou et al., (2019), we analyze the topology, electronic states, and lattice dynamics with density functional methods, including hybrid-functional calculations and machine-learned interatomic potentials. The Kohn–Sham orbitals near the Fermi level display a strong electron–phonon coupling, and exhibit large energy fluctuations at room temperature. The conduction tail states exhibit larger phonon-induced fluctuations than the valence tail states. To resolve transport at the atomic scale, we employ space-projected electronic conductivity and site-projected thermal conductivity methods. Local analysis of heat transport highlights the role of filamentary networks dominated by Te, with Sb and Ge making progressively smaller contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparing the reaction profiles of single iron catalytic sites in enzymes and in reticular frameworks for methane-to-methanol oxidation

The design of synthetic inorganic catalysts mimicking the first coordination spheres of enzymatic cofactors often results in lower yields and selectivity than their biological counterparts. In this study, we exploit Kohn-Sham density functional methods to compare the reaction profiles of four single iron-based catalysts for the direct oxidation of methane to methanol: two biomimetic models based on two enzymes (cytochrome P450 and taurine dioxygenase [TauD]) and two synthetic reticular frameworks (iron-BEA zeolite and tri-iron oxo-center-based metal-organic framework). Both the biomimetic and inorganic catalysts show almost zero selectivity toward methanol for methane conversions >1% at ambient temperature. Furthermore, this study highlights that iron’s first coordination shell can influence selectivity toward methanol but to a limited extent. In the absence of methanol protection strategies, high selectivity can be reached only by mimicking the reaction microenvironment of enzymes beyond the first coordination shell of iron.

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

Hydrogen and deuterium tunneling in niobium

We use density functional methods to identify the atomic configurations of H and D atoms trapped by O impurities embedded in bulk Nb. The O atoms are located at the octahedral position in the Nb body-centered cubic (BCC) lattice, and H (D) atoms tunnel between two degenerate tetrahedral sites separated by a mirror plane. Using nudged elastic band (NEB) methods, we calculate the double-well potential for O-H and O-D impurities and the wave functions and tunnel splittings for H and D atoms. Our results agree with those obtained from analysis of heat capacity and neutron scattering measurements on Nb with low concentrations of O-H and O-D.

density functional theory↗
DOE OSTI2020

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Thermal effects on the electronic properties of sodium electride under high pressures

Sodium is a simple metal at ambient conditions, while it transits to an electride phase at pressures above ~160 GPa along the room temperature isotherm. Here we explore the thermal effects on the electronic properties of the $hP4$ phase of sodium along the $ρ$ = 5.872 g/cm 3 isochore. We quantitatively classify this phase as an insulator based on the criterion of nearsightedness of the one-particle density matrix. Ab initio calculations suggest that the band gap of this insulator decreases with increasing temperature along the isochore, primarily because of ionic distortions; culminating in an insulator-to-metal transition upon melting at $T_m \approx 2100$ $\hbox{K}$. This transition is accompanied by residual electronic localization (in d orbitals) in the form of dynamic electron bubbles and a change in hybridization form $p–d$ to $s–p$ upon melting. This transition is explored by tracking the electronic and electro-optical properties along the isochore under consideration.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Xenon Trapping in Metal-Supported Silica Nanocages

Xenon (Xe) is a valuable and scarce noble gas used in various applications, including lighting, electronics, and anesthetics, among many others. It is also a volatile byproduct of the nuclear fission of uranium. In this work, a novel material architecture consisting of silicate nanocages in contact with a metal surface and an approach for trapping single Xe atoms in these cages is presented. The trapping is done at low Xe pressures and temperatures between 400 and 600 K, and the process is monitored in situ using synchrotron-based ambient pressure X-ray photoelectron spectroscopy. Release of the Xe from the cages occurs only when heating to temperatures above 750 K. A model that explains the experimental trapping kinetics is proposed and tested using Monte Carlo methods. Density functional theory calculations show activation energies for Xe exiting the cages consistent with experiments. This work can have significant implications in various fields, including Xe production, nuclear power, nuclear waste remediation, and nonproliferation of nuclear weapons. The results are also expected to apply to argon, krypton, and radon, opening an even more comprehensive range of applications.

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

The influence of temperature on the strain-hardening behavior of Fe-22/25/28Mn-3Al-3Si TRIP/TWIP steels

In this work, the influence of temperature and stacking fault energy (SFE) on the strain-hardening behavior and critical resolved shear stress for twinning was investigated for three Fe–22/25/28Mn–3Al–3Si wt.% transformation- and twinning-induced plasticity (TRIP/TWIP) steels. The SFEs were calculated by two different methods, density functional theory and statistical thermodynamic modeling. The dislocation structure, observed at low levels of plastic deformation, transitions from “planar” to “wavy” dislocation glide with an increase in temperature, Mn content, and/or SFE. The change in dislocation glide mechanisms from planar to wavy reduces the strain hardening rate, in part due to fewer planar obstacles and greater cross slip activity. Additionally, the alloys exhibit a large decrease in strength and ductility with increasing temperature from 25 to 200 °C, attributed to a substantial reduction in the thermally activated component of the flow stress, predominate suppression of TRIP and TWIP, and a significant increase in the critical resolved shear stress for mechanical twinning. Interestingly, the increase in SFE with temperature had a rather minor influence on the critical resolved shear stress for mechanical twinning, and other temperature dependent factors which likely play a more dominant role are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantum-Chemically Informed Machine Learning: Prediction of Energies of Organic Molecules with 10 to 14 Non-hydrogen Atoms

High-fidelity quantum-chemical calculations can provide accurate predictions of molecular energies, but their high computational costs limit their utility, especially for larger molecules. We have shown in previous work that machine learning models trained on high-level quantum-chemical calculations (G4MP2) for organic molecules with one to nine non-hydrogen atoms can provide accurate predictions for other molecules of comparable size at much lower costs. Here we demonstrate that such models can also be used to effectively predict energies of molecules larger than those in the training set. To implement this strategy, we first established a set of 191 molecules with 10–14 non-hydrogen atoms having reliable experimental enthalpies of formation. We then assessed the accuracy of computed G4MP2 enthalpies of formation for these 191 molecules. The error in the G4MP2 results was somewhat larger than that for smaller molecules, and the reason for this increase is discussed. Two density functional methods, B3LYP and ωB97X-D, were also used on this set of molecules, with ωB97X-D found to perform better than B3LYP at predicting energies. The G4MP2 energies for the 191 molecules were then predicted using these two functionals with two machine learning methods, the FCHL-Δ and SchNet-Δ models, with the learning done on calculated energies of the one to nine non-hydrogen atom molecules. The better-performing model, FCHL-Δ, gave atomization energies of the 191 organic molecules with 10–14 non-hydrogen atoms within 0.4 kcal/mol of their G4MP2 energies. Thus, this work demonstrates that quantum-chemically informed machine learning can be used to successfully predict the energies of large organic molecules whose size is beyond that in the training set.

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

Crystal Structure of Colloidally Prepared Metastable Ag 2 Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag 2 Se nanocrystals has traditionally been referred to as the “tetragonal” phase. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl 2 -like structure type (space group P2 1 /n) more accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag 2 Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that this anti-PbCl 2 -like Ag 2 Se polymorph is a dynamically stable, narrow-band-gap semiconductor. Here, the anti-PbCl 2 -like structure of Ag 2 Se is a low-lying metastable polymorph at 5–25 meV/atom above the ground state, depending on the exchange-correlation functional used.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

The efficacy of Lewis affinity scale metrics to represent solvent interactions with reagent salts in all-inorganic metal halide perovskite solutions

Solvents employed in the solution processing of metal halide perovskites are known to play a key role in defining the morphology and properties of the resulting thin film, and thus the performance of perovskite solar cell devices. Accurate metrics are needed that are capable of differentiating among candidates, finding solvents that adequately solubilize the various precursor species in solution and facilitate the nucleation and growth of these materials. Existing metrics such as the unsaturated Mayer bond order (UMBO) and the Gutmann donor number (DN) have been tested for lead iodide perovskite systems; but there has yet to be a comprehensive study on their transferability to lead-free perovskite solutions. Here, we use ab initio methods (density functional theory) and regression analysis tools to study the usefulness of DN and BF 3 affinity scales in this regard. We compared the relative effectiveness of these scales to describe interactions between solvents and BXn perovskite salts of lead (Pb 2+ ), tin (Sn 2+ and Sn 4+ ), germanium (Ge 2+ ), bismuth (Bi 3+ ), and antimony (Sb 3+ and Sb 5+ ). The DN proved to be a better representation than the BF 3 of such interactions, reflecting the closer similarity of these species to the “parent” SbCl 5 Lewis acid than to BF 3 . In addition, we have uncovered the usefulness of the lithium cation affinity metric (LCA) to describe the strength of interactions between solvents and A-site cations (e.g. Na + , K + , Rb + and Cs + ) in all-inorganic metal halide perovskite solutions. We find that the coordination strengths of solvents towards species in all-inorganic metal halide perovskite solutions are best described by two different metrics with distinct modes of action: DN differentiates among BX n salt complexes, and LCA among A-site cation species. This revelation can help guide the choice of solvent to optimize processing conditions. It also emphasizes the importance of selecting solvents whose DN and LCA optimize coordination to key Lewis acid species in all-inorganic perovskite solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗