Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Microwave-induced selective decomposition of cellulose: Computational and experimental mechanistic study

Understanding microwave-material interactions will help facilitate the utilization of microwave technology in gasification and renewable energy production. In this study, cellulose was used as a model compound to simulate the organic matter in biomass, and its catalytic decomposition under a microwave (MW) field was studied to identify structural changes from the reaction. The study was conducted using a MW source coupled to a fixed-bed gas-flow reactor, mass spectroscopic, and Fourier transform infrared spectroscopy post-reaction analysis. Zeolite 13X was chosen as a microwave absorber to study the catalytic enhancement of the decomposition of cellulose. Density functional theory (DFT) was used to gain insights into the molecular transformations occurring in the presence of a static electric field, which was used to simulate the electric field component of the microwave electromagnetic radiation. Theoretical calculations demonstrated that both the positive and negative portion of the electric field interact with the permanent dipoles of the cellulose effecting the decomposition through the glycosidic bond breaking mechanism. The theoretical result was verified using infrared spectroscopic analysis of the pure cellulose during microwave heating. The theoretical calculations suggest only the positive electric field component may be active for the cellulose decomposition in the presence of Zeolite 13X. Physically mixing Zeolite 13X with the cellulose led to a significant enhancement in the decomposition rate of the glycosidic bond. Zeolite 13X enhanced the glycosidic O-C decomposition at lower MW power (lower temperatures), whereas the O-H functional group required higher MW power (higher temperature) for its decomposition. The DFT study coupled with the reaction studies revealed that the electric field polarizability is dependent upon both the direction and the orientation of the cellulose. Gas products revealed that applying 250 W of MW power led to the production of CO, H 2 , along with some CO 2 , CH 4 , and benzene at 305 °C. Finally, reaction under 500, 750, and 1000 W of power at constant temperature (305 °C) revealed that higher power led to the complete decomposition of cellulose to mostly CO, H 2 , and CO 2 .

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

Methanol adsorption and dissociation on GaP(110) studied by ambient pressure X-ray photoelectron spectroscopy

Ambient pressure X-ray photoelectron spectroscopy (AP-XPS) was used to investigate methanol (CH 3 OH) adsorption and reaction on the GaP(110) surface. Exposure of CH 3 OH to GaP(110) at room temperature led to the formation of at least four different surface species as indicated by analysis of C 1s and O 1s XPS features. By combining AP-XPS data with density functional theory calculations, the surface species were identified as methoxy (CH 3 O*), formaldehyde (CH 2 O*), and paired methanol (p-CH 3 O*H) and methoxy (p-CH 3 O*) species, where “paired” means that they belong to a hydrogen-bonded methoxy-methanol complex. Asterisk * here indicates an adsite. The formation of CH 2 O* via the dehydrogenation of CH 3 O* was shown to be limited by the availability of vacant phosphorus (P) sites on GaP(110). With an increase in CH 3 OH pressure, the fractional coverage of CH 3 O* species reached 0.55, and the surface P sites were completely saturated with hydrogen. Under a constant CH 3 OH pressure of 0.5 Torr, the surface concentration of the paired species and of CH 2 O* remained constant until 400 K. At higher temperatures, thermally driven reactions led to a significant increase in the concentration of surface CH x * species, which suggests that C-O bond cleavage of the CH 3 O group is the dominant decomposition mechanism on GaP(110). In conclusion, based on the reactivity of GaP(110) toward CH 3 OH dehydrogenation, elevated temperatures and CH 3 OH pressures may be used to functionalize this surface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Crystal structure of halofuginone hydrobromide, C 16 H 18 BrClN 3 O 3 Br

The crystal structure of one form of halofuginone hydrobromide has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Halofuginone hydrobromide crystallizes in space group P2 1 (#4) with a = 8.87398(13), b = 14.25711(20), c = 15.0153(3) Å, β = 91.6867(15)°, V = 1898.87(4) Å 3 , and Z = 4. The crystal structure consists of alternating layers (parallel to the ab-plane) of planar and nonplanar portions of the cations. N–H∙∙∙Br and O–H∙∙∙Br hydrogen bonds link the protonated piperidine rings and bromide anions into a two-dimensional network parallel to the ab-plane. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal structure of oxibendazole, C 12 H 15 N 3 O 3

The crystal structure of oxibendazole has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Oxibendazole crystallizes in space group C2/c (#15) with a = 23.18673(22), b = 5.35136(5), c = 19.88932(13) Å, β = 97.0876(9)°, V = 2449.018(17) Å 3 , and Z = 8. The structure consists of hydrogen-bonded layers of planar molecules parallel to the bc-plane. Strong N–H∙∙∙N hydrogen bonds link the molecules into dimers, with a graph set R2,2(8). N–H∙∙∙O hydrogen bonds further link these dimers into layers parallel to the bc-plane. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring Ca–Ce–M–O (M = 3d Transition Metal) Oxide Perovskites for Solar Thermochemical Applications

Solar thermochemical (STC) processes hold promise as efficient ways to generate renewable fuels, fuel precursors, or chemical feedstocks using concentrated sunlight. Specifically, one actively researched approach is the two-step STC cycle, which uses a redox-active, off-stoichiometric, transition-metal oxide material to split water and/or CO 2 , generating H 2 and/or CO, respectively, or syngas (a combination of H 2 and CO). Identifying novel metal oxides that yield larger reduction extents (practically achievable off-stoichiometries) than the state-of-the-art CeO 2 is critical. Here, we explore the chemical space of Ca–Ce–M–O (M = 3d transition metal, except Cu and Zn) metal oxide perovskites, with Ca and/or Ce occupying the A site and M occupying the B site within an ABO 3 framework, as potential STC candidates. We use density functional theory (DFT)-based calculations and systematically evaluate the oxygen vacancy (VaO) formation energy (≈ enthalpy of reduction in an STC cycle), electronic properties, thermodynamic stability of CaMO 3 , CeMO 3 , and Ca 0.5 Ce 0.5 MO 3 perovskites, and the VaO formation energy within Ca 0.5 Ce 0.5 Ti 0.5 Mg 0.5 O 3 perovskite. We consider only Ca and/or Ce on the A site because of their similar size and the potential redox activity of Ce 4+ . If both Ce and M exhibit simultaneous reduction with Va O formation, the resulting perovskite could exhibit a larger entropy of reduction than a single cation reduction. The increased entropy produces increased reduction for fixed temperature, partial pressure of oxygen, and reduction enthalpy, and therefore increased STC efficiency. Importantly, we identify Ca 0.5 Ce 0.5 MnO 3 , Ca 0.5 Ce 0.5 FeO 3 , and Ca 0.5 Ce 0.5 VO 3 to be promising candidates based on their Va O formation energy and thermodynamic (meta)stability. Moreover, based on our calculated on-site magnetic moments, electron density of states, and electron density differences between pristine and defective structures, we find Ca 0.5 Ce 0.5 MnO 3 to exhibit simultaneous reduction of both Ce 4+ (A-site) and Mn 3+ (B-site), highlighting a particularly promising candidate for STC applications with a predicted higher entropy of reduction than CeO 2 . Lastly, we extract metrics that govern the trends in Va O formation energies, such as standard reduction potentials, and provide pointers for further experimental and theoretical studies, which will enable the design of improved materials for the STC cycle.

14 SOLAR ENERGY↗
DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

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

Thermochemical Data Fusion Using Graph Representation Learning

Large databases are required for “Big Data” applications in catalysis and materials science. Thermochemical databases can be created by combining data from various sources and by correcting low-fidelity datasets to higher accuracy with minimal computation. To achieve this “data fusion”, thermochemical quantities of interest, calculated at various levels of density functional theory (DFT), need to be mapped to the same, high levels of theory. In this work, a graph theoretical, statistical framework is proposed for such tasks. Subgraph frequencies are shown to provide a natural representation for learning these fusion maps. The maps are linear and are learnt with automated descriptor selection. Using a dataset of as few as ~1% from the QM9 database of 133,885 molecules, these models can predict multiple thermochemical quantities at a higher level of theory with an accuracy of 1 kcal/mol. Here, the method is explainable, generalizable, and provides a diagnostic tool for outlier identification

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

Prediction of organic homolytic bond dissociation enthalpies at near chemical accuracy with sub-second computational cost

Bond dissociation enthalpies (BDEs) of organic molecules play a fundamental role in determining chemical reactivity and selectivity. However, BDE computations at sufficiently high levels of quantum mechanical theory require substantial computing resources. In this paper, we develop a machine learning model capable of accurately predicting BDEs for organic molecules in a fraction of a second. We perform automated density functional theory (DFT) calculations at the M06-2X/def2-TZVP level of theory for 42,577 small organic molecules, resulting in 290,664 BDEs. A graph neural network trained on a subset of these results achieves a mean absolute error of 0.58 kcal mol-1 (vs DFT) for BDEs of unseen molecules. We further demonstrate the model on two applications: first, we rapidly and accurately predict major sites of hydrogen abstraction in the metabolism of drug-like molecules, and second, we determine the dominant molecular fragmentation pathways during soot formation.

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

Electronic properties of single-layer CoO 2 /Au(111)

Here we report direct measurements via angle-resolved photoemission spectroscopy (ARPES) of the electronic dispersion of single-layer (SL) CoO 2 . The Fermi contour consists of a large hole pocket centered at the $\overline{\Gamma}$ point. To interpret the ARPES results, we use density functional theory (DFT) in combination with the multi-orbital Gutzwiller Approximation (DFT+GA), basing our calculations on crystalline structure parameters derived from x-ray photoelectron diffraction and low-energy electron diffraction. Our calculations are in good agreement with the measured dispersion. We conclude that the material is a moderately correlated metal. We also discuss substrate effects, and the influence of hydroxylation on the CoO 2 SL electronic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a -1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Co(salen)-Catalyzed Oxidation of Lignin Models to Form Benzoquinones and Benzaldehydes: A Computational and Experimental Study

Lignin is a highly abundant polyphenolic polymer that imparts mechanical strength to plant biomass. Transition-metal complexes can catalyze lignin oxidation to produce value-added products, but low catalytic efficiency has hampered their use in industry. Identifying the chemical and structural factors that govern catalytic activity is a prerequisite to rational design of catalysts with improved activity. In this study, we combine computational and experimental approaches to investigate the mechanism of Co(salen)-catalyzed oxidation of the monomeric lignin models syringyl (S), vanillyl (G), and 4-hydroxybenzyl alcohol (H) to produce benzoquinone and benzaldehyde products. Experimentally, S oxidation to form dimethoxybenzoquinone proceeded efficiently with a Co(salen) catalyst coordinated by a pyridine ligand, but G and H did not undergo oxidation. Density functional theory calculations reveal that catalyst regeneration is energetically unfavorable in the presence of H, which prevents oxidation. In contrast, S readily facilitates catalyst regeneration. Formation of methoxybenzoquinone from G was achieved experimentally by adding bulky, noncoordinating bases. These findings provide a fundamental baseline for enhancing the activity of Co-Schiff base catalysts toward lignin-like molecules by adding sterically hindered nitrogenous bases or potentially by including a cocatalyst that promotes catalyst regeneration.

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

Thermally frustrated phase transition at high pressure in B 2-ordered FeV

X-ray diffraction measurements of equiatomic B2-ordered FeV were performed in a diamond-anvil cell at room temperature at several pressure points up to 80 GPa that showed the cubic phase to be stable with no indication of structural phase transitions. Density functional theory at 0 K predicts Fermi surface nesting, an electronic topological transition, and a phonon dynamical instability within the experimentally investigated pressure range. Nevertheless, the instability is absent in phonon dispersion curves extracted from ab initio molecular dynamics simulations below the critical volume at temperatures as low as 10 K, indicating that thermal atomic displacements can frustrate the phase transition by renormalizing the phonon dispersion curves. Ferrimagnetism is critical for the stability of the cubic phase at low temperature, but thermal atomic displacements are enough to support the structure at and above the Néel temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Weakly coupled alternating S = 1 2 chains in the distorted honeycomb lattice compound Na 2 Cu 2 TeO 6

Spin- 1 2 chains with alternating antiferromagnetic (AF) and ferromagnetic (FM) couplings exhibit quantum entanglement as the integer-spin Haldane chains and might be similarly utilized for quantum computations. Such alternating AF-FM chains have been proposed to be realized in the distorted honeycomb lattice compound Na 2 Cu 2 TeO 6 , but to confirm this picture a comprehensive understanding of the exchange interactions including terms outside of the idealized model is required. Here, we employ neutron scattering to study the spin dynamics in Na 2 Cu 2 TeO 6 and accurately determine the coupling strengths through the random phase approximation and density functional theory approaches. We find the AF and FM intrachain couplings are the dominant terms in the spin Hamiltonian, while the interchain couplings are AF but perturbative. This hierarchy in the coupling strengths and the alternating signs of the intrachain couplings can be understood through their different exchange paths. Our results establish Na 2 Cu 2 TeO 6 as a weakly coupled alternating AF-FM chain compound and reveal the robustness of the gapped ground state in alternating chains under weak interchain couplings.

1-dimensional spin chains↗
DOE OSTI2020

Temporally decoherent and spatially coherent vibrations in metal halide perovskites

The long carrier lifetime and defect tolerance in metal halide perovskites (MHPs) are major contributors to the superb performance of MHP optoelectronic devices. Large polarons were reported to be responsible for the long carrier lifetime. Yet microscopic mechanisms of the large polaron formation, including the so-called phonon melting, are still under debate. In this study, time-of-flight inelastic neutron scattering experiments and first-principles density-functional theory calculations were employed to investigate the lattice vibrations (or phonon dynamics) in methylammonium lead iodide, a prototypical example of MHPs. Our findings are that optical phonons lose temporal coherence gradually with increasing temperature which vanishes at the orthorhombic-to-tetragonal structural phase transition. Surprisingly, however, we found that the spatial coherence is still retained throughout the decoherence process. We argue that the temporally decoherent and spatially coherent vibrations contribute to the formation of large polarons in this metal halide perovskite.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Isotope effect and critical magnetic fields of superconducting YH 6 : A Migdal-Eliashberg theory approach

The emergence of near-ambient temperature superconductivity under pressure in metal hydride systems has motivated a desire to further understand such remarkable properties, specifically critical magnetic fields. Here, YH 6 is suggested to be a departure from conventional superconductivity, due to apparent anomalous behavior. Using density functional calculations in conjunction with Migdal-Eliashberg theory we show that in YH 6 the critical temperature and the isotope effect under pressure, as well as the high critical fields, are consistent with strong-coupling conventional superconductivity, a property anticipated to extend to other related systems. Furthermore, strong-coupling corrections occur to the expected BCS values for the isotope effect coefficient (α), Ginzburg-Landau parameter [κ 1 ⁡(T)], London penetration depth [λ L (T)], electromagnetic coherence length [ξ⁡(T)], and the energy gap (Δ 0 ).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling of aqueous species interaction energies prior to nucleation in cement-based gel systems

Arguably the most ubiquitous construction material in modern civilization, concrete is enabling the development of megacities around the globe together with increasing living standards in developing nations. However, the exact formation mechanisms of the strength-giving calcium-rich gels remain a topic of debate. Using density functional modeling, we simulate the fundamental solution-based building blocks of cement hydrates (calcium ions and silicate and aluminate monomers) and their propensity to form pair-wise complexes with bonding environments characteristic of those found in calcium-silicate-hydrate, calcium-alumino-silicate-hydrate and sodium-containing calcium-alumino-silicate-hydrate gels, as assessed from Gibbs free energies of chemical reactions. By accurately simulating the high pH pore solution chemistry in Portland cements and related systems, along with discrete solvation of the species, we hypothesize potential early age formation routes of the gels and discuss limitations and future work associated with this approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal structure of toceranib, C 22 H 25 FN 4 O 2

The crystal structure of toceranib has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Toceranib crystallizes in space group P2 1 /c (#14) with a = 10.6899(6), b = 24.5134(4), c = 7.8747(4) Å, β = 107.7737(13)°, V = 1965.04(3) Å 3 , and Z = 4. The crystal structure consists of stacks of approximately planar molecules, with N–H∙∙∙O hydrogen bonds between the layers. The commercial reagent sample was a mixture of two or more phases with toceranib being the dominant phase. The difference between the Rietveld-refined and DFT-optimized structures is larger than usual. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Designing order–disorder transformation in high-entropy ferritic steels

Order–disorder transformations hold an essential place in chemically complex high-entropy ferritic steels (HEFSs) due to their critical technological application. The chemical inhomogeneity arising from mixing of multi-principal elements of varying chemistry can drive property altering changes at the atomic scale, in particular short-range order. Using density-functional theory-based linear-response theory, we predict the effect of compositional tuning on the order–disorder transformation in ferritic steels—focusing on Cr–Ni–Al–Ti–Fe HEFSs. We show that Ti content in Cr–Ni–Al–Ti–Fe solid solutions can be tuned to modify short-range order that changes the order–disorder path from BCC-B2 (Ti atomic-fraction = 0) to BCC-B2-L2 1 (Ti atomic-fraction > 0) consistent with existing experiments. Our study suggests that tuning degree of SRO through compositional variation can be used as an effective means to optimize phase selection in technologically useful alloys.

36 MATERIALS SCIENCE↗