Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

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

Validation of the Cossee–Arlman mechanism for propylene oligomerization on Ni/UiO-66

Steady state rate expressions can be derived to distinguish the Cossee–Arlman and metallacycle mechanisms postulated for propylene oligomerization on nickel-based catalysts based on product selectivities, where product selectivities for the former are a function of olefin pressure because sequential coordination and insertion steps lead to independent mechanistic pathways for different hexene isomers. In contrast, the metallacycle mechanism presents pressure-independent product selectivities due to successive coordination prior to the kinetically relevant steps in each mechanism. In this work, steady state propylene oligomerization rates and selectivities were measured in the absence of an activator on nickel functionalized UiO-66 metal organic framework (MOF), Ni/UiO-66, to validate the Cossee–Arlman mechanism for light olefin oligomerization. In situ NO titrations reveal that ~5% of nickel sites were active during the reaction, and thus, not all nickel sites are relevant for catalysis. Propylene dimerization was first order in propylene pressure from 5 to 500 kPa with an apparent activation energy of ~20 kJ mol -1 from 453 to 493 K. Calculated apparent activation energies with density functional theory (DFT) calculations on cluster models of Ni/UiO-66 are in agreement with experiment to corroborate the Cossee–Arlman mechanism. Selectivities of hexene products and the ratio of hexene product selectivities on Ni/UiO-66 are in accordance with selectivity expressions derived from the Cossee–Arlman mechanism. In conclusion, analysis of product selectivities can be used more extensively to demarcate the Cossee–Arlman and metallacycle mechanisms for olefin oligomerization on metal-based catalysts.

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

Scale-invariant machine-learning model accelerates the discovery of quaternary chalcogenides with ultralow lattice thermal conductivity

We design an advanced machine-learning (ML) model based on crystal graph convolutional neural network that is insensitive to volumes (i.e., scale) of the input crystal structures to discover novel quaternary chalcogenides, AMM'Q 3 (A/M/M' = alkali, alkaline earth, post-transition metals, lanthanides, and Q = chalcogens). These compounds are shown to possess ultralow lattice thermal conductivity (κ l ), a desired requirement for thermal-barrier coatings and thermoelectrics. Upon screening the thermodynamic stability of ~1 million compounds using the ML model iteratively and performing density-functional theory (DFT) calculations for a small fraction of compounds, we discover 99 compounds that are validated to be stable in DFT. Taking several DFT-stable compounds, we calculate their κ l using Peierls–Boltzmann transport equation, which reveals ultralow κ l (<2 Wm -1 K -1 at room temperature) due to their soft elasticity and strong phonon anharmonicity. Our work demonstrates the high efficiency of scale-invariant ML model in predicting novel compounds and presents experimental-research opportunities with these new compounds.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Association Mechanisms of Unsaturated C2 Hydrocarbons with Their Cations: Acetylene and Ethylene

The ion-molecule association mechanism of acetylene and ethylene with their cations is investigated by ab initio quantum chemical methods to understand the structures, association energies, and the vibrational and electronic spectra of the products. Stable puckered cyclic isomers are found as the result of first forming less stable linear and bridge isomers. The puckered cyclic complexes are calculated to be strongly bound, by 87, 35 and 56 kcal/mol for acetylene-acetylene cation, ethylene-ethylene cation and acetylene-ethylene cation, respectively. These stable complexes may be intermediates that participate in further association reactions. There are no association barriers, and no significant inter-conversion barriers, so the initial linear and bridge encounter complexes are unlikely to be observable. However, the energy gap between the bridged and cyclic puckered isomers greatly differs from complex to complex: it is 44 kcal/mol in C4H4 +, but only 6 kcal/mol in C4H8 +. The accurate CCSD(T) calculations summarized above are also compared against less computationally expensive MP2 and density functional theory (DFT) calculations for structures, relative energies, and vibrational spectra. Calculated vibrational spectra are compared against available experiments for cyclobutadiene cation. Electronic spectra are also calculated using time-dependent DFT.

ethylene↗
DOE OSTI2022

Origin of ultralow phonon transport and strong anharmonicity in lead-free halide perovskites

All-inorganic lead-free halide double perovskites offer a promising avenue toward non-toxic, stable optoelectronic materials, properties that are missing in their prominent lead-containing counterparts. Their large thermopowers and high carrier mobilities also make them promising for thermoelectric applications. Here, in this work, we present a first-principles study of the lattice vibrations and thermal transport behaviors of Cs 2 SnI 6 and γ-CsSnI 3 , two prototypical compounds in this materials class. We show that conventional static zero temperature density functional theory (DFT) calculations severely underestimate the lattice thermal conductivities (κ l ) of these compounds, indicating the importance of dynamical effects. By calculating anharmonic renormalized phonon dispersions, we show that some optic phonons significantly harden with increasing temperature (T), which reduces the scattering of heat carrying phonons and enhances calculated κ l values when compared with standard zero temperature DFT. Furthermore, we demonstrate that coherence contributions to κ l , arising from wave like phonon tunneling, are important in both compounds. Overall, calculated κ l with temperature-dependent interatomic force constants, built from particle and coherence contributions, are in good agreement with available measured data, for both magnitude and temperature dependence. Large anharmonicity combined with low phonon group velocities yield ultralow values, with room temperature values of 0.26 W/m-K and 0.72 W/m-K predicted for Cs 2 SnI 6 and γ-CsSnI 3 , respectively. We further show that the lattice dynamics of these compounds are highly anharmonic, largely mediated by rotation of the SnI6 octahedra and localized modes originating from Cs rattling motion. These thermal characteristics combined with their previously computed excellent electronic properties make these perovskites promising candidates for optoelectronic and room temperature thermoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine-Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

To cope with legacy greenhouse gas emissions and to achieve net-zero emissions by 2050, the U.S. Department of Energy (DOE) is funding efforts to develop direct air capture (DAC), a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been well studied as DAC sorbent materials due to their tunable structural and compositional properties. Thermodynamic simulations using force fields are often used to provide predictions of a material’s performance in many separations. However, these force fields often make assumptions about bonds and the physics of the adsorption process. A new class of force fields called machine-learned force fields (MLFFs) use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). In this work, models were developed to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using MLFFs. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗
DOE OSTI2023

Interpretable, Extensible Linear Regression Model for Electron Density Prediction

This code allows us to generate a model to predict ground state electron densities directly from atomic structure thereby bypassing the need to perform direct DFT calculations. It is based on many body correlation descriptors defined in https://arxiv.org/abs/2004.14442 (IM release numbers: LLNL-JRNL- 808791 and LLNL-JRNL-840698). The code reads in electron densities, in CHGCAR format used by the DFT code VASP, calculates many-body correlation descriptors based on atoms surrounding a grid point and then uses linear regression to obtain the fitting parameters.

Kumar, Shashikant↗
DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2021

Elimination of the linearization error in APW/LAPW basis set: Dirac-Kohn-Sham equations

Here, a detailed account of the implementation of equations of the relativistic density functional theory (RDFT) using basis sets of APW/LAPW type with flexible extensions provided by local orbitals is given. Earlier discoveries of the importance of the high derivative local orbital (HDLO) extension of the APW/LAPW basis set for enhancing the accuracy of DFT calculations are confirmed using a fully relativistic approach and α – U as an example. High energy local orbitals (HELO's), however indispensable for GW calculations, are considerably less efficient in enhancing the accuracy of DFT applications. It is shown that a simplified approach to the relativistic effects, namely considering them only inside the muffin-tin (MT) spheres, produces basically identical results (as compared to fully relativistic approach) for the electronic free energy of the five materials considered in this work. By comparing the effect of the simplified approach on the electronic free energy with its effect on the electronic kinetic energy we conclude that the insensitivity of the free energy to the way we describe the relativistic effects in the interstitial region is related to the variational property of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Surfactant inhibition mechanisms of carbonate mineral dissolution in shale

Surfactants are common additives to hydraulic fracturing and enhanced oil recovery (EOR) fluids, and are under consideration for amendment to supercritical carbon dioxide for geological carbon sequestration (GCS). The effect of a common anionic surfactant, internal olefin sulfonate (IOS), on mineral dissolution from shale into brine was evaluated. When added to brine at concentrations exceeding the critical micelle concentration (94 mg/L), IOS inhibited carbonate mineral dissolution in an Eagle Ford shale, as well as dissolution of optical quality calcite (the dominant carbonate in the shale). Laser profilometry images provide spatial resolution across > 3 orders of magnitude, and indicate that IOS addition to brine both enhances the formation of new etch pits in calcite, and impedes their further growth. Time-of-flight secondary ion mass spectrometry surface profiles show for the first time that IOS preferentially adsorbs at calcite pit edges versus flat calcite surfaces (i.e., terraces). Surface pressure calculations, sulfur K-edge near edge X-ray absorption fine structure (NEXAFS) spectroscopy results, and density functional theory (DFT) calculations support this observation; the DFT results indicate that the sulfonate head group of the IOS molecule binds strongly to the calcite step site as compared to the terrace site. The S K-edge NEXAFS results indicate that IOS adsorbed more to etched calcite surfaces compared to smooth calcite surfaces. Overall, the results indicate that weak adsorption on flat calcite surfaces (i.e., terraces) disrupts water structure and enhances mass transfer of dissolution, while strong adsorption on calcite pit edges displaces adsorbed water and inhibits further etch pit growth. This work provides the first direct evidence of preferential adsorption of IOS to etched calcite surfaces and links it to macroscopic dissolution kinetics. Finally, this work has implications for surfactant-containing fluids used in hydraulic fracturing, EOR and potentially GCS for subsurface injection into carbonate rich reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unexpected features in the optical vibrational spectra of δ-UO 3

Uranium trioxide displays a complex chemical phase space, with at least six structurally distinct polymorphs accessible via different synthetic routes. Remarkably, despite its technological importance, full structural and electronic characterization of these polymorphs remains an open area of study. δ-UO 3 in particular has attracted significant theoretical attention due to its high point group and space group symmetries, having U (VI) in octahedral coordination with polyhedra interconnected through corner-sharing to build a 3-D cubic lattice with space group symmetry Pm-3m and Z = 1. Critical experimental information, such as its optical vibrational spectra, are not known. Here, we study the Raman and infrared (IR) spectra of δ-UO 3 together with the support of density functional theory (DFT) calculations for spectral interpretation. A symmetry analysis of the DFT-predicted phonon eigenmodes indicates that δ-UO 3 should have two IR active modes and no Raman active modes. Experimental results, however, indicate significant Raman scattering from δ-UO 3 . We therefore propose four potential explanations for this apparent contradiction: a possible tetragonal distortion to the cubic cell, the existence of a surface impurity layer, vacancy scattering, and structural activation of Raman signal. We use powder X-ray diffraction and confocal Raman spectroscopy with depth profiling to investigate these possibilities and suggest future experiments to explore this phenomenon in more detail. Understanding the lattice dynamics of δ-UO 3 is important for identification of technogenic U phases via Raman and infrared spectroscopy and our results indicate that the simple understanding of δ-UO 3 as a high-symmetry cubic structure should be reconsidered.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Single crystal structure, vibrational spectroscopy, gas sorption and antimicrobial properties of a new inorganic acidic diphosphates material (NH 4 ) 2 Mg(H 2 P 2 O 7 ) 2 • 2 H 2 O

We report on the successful synthesis of diammonium magnesium dihydrogendiphosphate (V) dihydrate compound (NH 4 ) 2 Mg(H 2 P 2 O 7 ) 2 •2H 2 O using a wet chemical route. Single crystal X-ray diffraction analysis and micro Raman spectroscopy are employed to characterize the compound. We demonstrate, using a multidisciplinary approach, that this compound is ideal for carbon dioxide (CO 2 ) capture in addition to other anthropogenic gasses. We show here -from both an experimental as well as from a density functional theory (DFT) calculations routes- the potential for adopting this compound into domestic air-conditioning units (ACUs). From these experiments, the resistance to bacterial growth is also investigated, which is critical for the adoption of this compound in ACUs. Our compound exhibits a higher methane (CH 4 ) sorptivity as compared to CO 2 at 25 °C and 45 °C under pressures up to 50 bars. Furthermore, DFT electronic structure calculations are used to compute the main structural and electronic properties of the compound, taking into consideration the characteristics of the identified pores as a function of the progressive CO 2 vs. CH 4 loadings. Finally, the antibacterial assay reveals a strong antibacterial activity against the tested Gram-positive and Gram-negative bacteria, with a large zone of inhibition against the tested E. Coli, S. Aureus and K. Pneumonia

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

Rhenium-Selenido Corroles: Reflections on 5d Metalloporphyrins and Metallocorroles as Triplet Emitters and Photosensitizers

Following the successful synthesis of rhenium-oxido, rhenium-imido, and rhenium-sulfido corroles, a series of para-X-substituted rhenium-selenido triarylcorroles (X = OCH 3 , CH 3 , H, F, CF 3 ) have been prepared in 70–76% yields. A one-pot, two-step procedure was used, analogous to that used for ReS corroles, involving rhenium insertion via high-temperature reaction of Re 2 (CO) 10 and free-base corroles in refluxing decalin, followed by exposure to PCl 3 (which is thought to deoxygenate ReO corroles formed under the reaction) and elemental selenium. An analogous procedure, however, failed to yield rhenium-tellurido corroles, presumably reflecting, according to DFT calculations, the limited stability of these species. Unlike ReO corroles, ReS and ReSe corroles were found not to exhibit phosphorescence in the NIR-I regime (600–1000 nm); nor did they sensitize singlet oxygen formation. In the hope of obtaining a broader perspective of this negative result, we also examined ReN porphyrins and found them not to phosphoresce or to sensitize singlet oxygen formation. These results were explained by DFT and TDDFT calculations in terms of low-energy triplet states, which are not energetic enough to excite molecular oxygen to its lowest singlet state. Whether some of the new Re corroles exhibit phosphorescence in the NIR-II regime remains an interesting question for future investigation.

Aromatic compounds↗
DOE OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simple self-consistent method for excited states in density functional theory to characterize defect-derived behavior in wide-band-gap-based microelectronic materials

This final report summarizes the results of the Laboratory Direct Research and Development (LDRD) Project Number 229740. Wide band gap semiconductors such as gallium nitride (GaN) have features highly desirable for multiple mission electronic applications. Realization of their potential requires atomic-scale understanding of electronic behavior. The principal experimental tools for electronically probing defects in GaN are chemically undifferentiating and lack a practical theoretical counterpart needed to identify and characterize specific defects. This project investigated whether a simple idea for modeling defect excited states and their associated photoluminescence (PL) energies is viable, as a path to accelerate the understanding of defect behavior and gain valuable insights into engineering new electronic materials and devices. The research implemented a non-self-consistent total-energy evaluation of a Koopmans-type estimation of an excited electronic state energy in density functional theory (DFT) calculations, and proceeded to design, implement, and assess a self-consistent method for computing excited states based upon an OCcupation-Constrained-DFT (occ-DFT). The occ-DFT was verified in test calculations of defect excited states and validated against well-characterized PL data for 3d transition metal defects in GaN. The method proved stable and robust in computing excited states and gave accurate predictions compared to experimental PL data. The combined ground state/excited-state capability proved capable of chemically differentiating defect species in GaN. In application to 3d dopants in GaN, we reinterpreted extensive experimental literature, proposed new defects as prospective candidates for use in quantum information applications, and outlined design strategies to create and exploit these potentially useful functional defects in GaN.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Quantum Monte Carlo Benchmarking of Molecular Adsorption on Graphene-Supported Single Pt Atom

The precise understanding of adsorption energetics and molecular geometry at catalytic sites is fundamental for advancing catalysis, particularly under the constraints of resource efficiency and environmental sustainability. Here, this study benchmarks the performance of density functional theory (DFT) calculations against diffusion Monte Carlo (DMC) calculations for adsorption properties of small gas molecules relevant to CO oxidation—namely O 2 , CO, CO 2 , and atomic oxygen—on a single Pt atom supported by pristine graphene. Our findings reveal that DMC calculations provide a significantly different landscape of adsorption energetics compared to DFT results. Notably, DFT predicts different lowest-energy configurations and spin states, particularly for O 2 , which suggests potential discrepancies in predicting the catalytic behavior. Furthermore, this study identifies the critical issue of CO poisoning, highlighted by the large disparity between the DMC adsorption energies of O 2 (−1.23(2) eV) and CO (−3.37(1) eV), which can inhibit the catalytic process. These results emphasize the necessity for more sophisticated computational approaches in catalysis research, aiming to refine the prediction accuracy of reaction mechanisms and to enhance the design of more effective catalysts.

Ahn, Jeonghwan [University of Illinois at Urbana-C↗
DOE OSTI2022

Bandgap and electronic structure of CaSiN 2 : Experiment and theory

Abstract We interrogate the electronic properties of the nitridosilicate CaSiN 2 as determined using soft X‐ray absorption and emission spectroscopy, which probe directly the partial electronic density of states. We compare our measurements to full potential, all electron density functional theory (DFT) calculations. Good agreement between experiment and theory is obtained, and the electronic bandgap of CaSiN 2 is determined to be 3.95 ± 0.3 eV, in contrast with our predicted value of 4.78 eV that was calculated using the modified Becke–Johnson exchange‐correlation functional with the DFT framework.

Materials Science↗
DOE OSTI2023

Potential Lifshitz transition at optimal substitution in nematic pnictide Ba 1 x Sr x Ni 2 As 2

BaNi 2 As 2 is a structural analog of the pnictide superconductor BaFe 2 As 2 , which, like the iron-based superconductors, hosts a variety of ordered phases including charge density waves (CDWs), electronic nematicity, and superconductivity. Upon isovalent Sr substitution on the Ba site, the charge and nematic orders are suppressed, followed by a sixfold enhancement of the superconducting transition temperature (T c ). To understand the mechanisms responsible for enhancement of T c , we present high-resolution angle-resolved photoemission spectroscopy (ARPES) measurements of the Ba 1-x Sr x Ni 2 As 2 series, which agree well with our density functional theory (DFT) calculations throughout the substitution range. Analysis of our ARPES-validated DFT results indicates a Lifshitz transition and reasonably nested electron and hole Fermi pockets near optimal substitution where T c is maximum. These nested pockets host Ni d xz /d yz orbital compositions, which we associate with the enhancement of nematic fluctuations, revealing unexpected connections to the iron-pnictide superconductors. This gives credence to a scenario in which nematic fluctuations drive an enhanced T c .

36 MATERIALS SCIENCE↗