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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

Assessing the Performance of Exchange‐Correlation Density Functionals in Describing the Iron‐Catalyzed Ammonia Synthesis System

Density functional theory (DFT) has been widely employed for elucidating mechanistic aspects of heterogeneous catalysis. However, the accuracy of DFT calculations relies heavily on selecting an exchange-correlation (XC) functional that correctly describes the electronic structure of materials involved in the reactions. This study assesses the accuracy of several XC density functionals for modeling the iron-catalyzed ammonia synthesis system. In the assessment of functional accuracy, experimental references are compared to DFT-calculated values for the formation energy of gas-phase ammonia and nitrogen, bulk Fe/Fe-nitride (γ′-Fe 4 N) lattice constants and cohesive/formation energies, and nitrogen and ammonia binding energies on Fe(100), Fe(111), Fe(110), and γ′-Fe 4 N(111). It is observed that the experimental value for each of these descriptors is accurately modeled by at least one functional. RPBE alone provides reliable estimates for both the lattice constant and cohesive energy of Fe and γ′-Fe 4 N. Temperature-programmed desorption experiments led to estimates for N and NH 3 adsorption across several Fe-based facets that are best captured by RPBE. These results highlight the importance of choosing an appropriate XC functional that accurately describes experimental systems and offer insights into effectively modeling the interactions between nitrogen and ammonia on Fe-based surfaces.

adsorption↗
DOE OSTI2021

Catalytic CO Oxidation on MgAl 2 O 4 -Supported Iridium Single Atoms: Ligand Configuration and Site Geometry

Probing and understanding the local chemical environment of an active site is essential for designing high-performance single-atom catalysts (SACs). Density functional theory (DFT) calculations were performed to investigate the ligand configuration and site geometry of MgAl 2 O 4 -supported iridium single atoms (Ir-1) toward catalytic carbon monoxide (CO) oxidation. We employed MgAl 2 O 4 (111) and MgAl 2 O 4 (211) as the model substrates with adsorbed Ir single atoms of different site geometries. DFT calculations revealed that the Mg-site on MgAl 2 O 4 (111) and the step site on MgAl 2 O 4 (211) are the most stable adsorption sites for Ir single atoms. Irrespective of site choices, CO oxidation on supported Ir single atoms follows the Eley-Rideal (E-R) mechanism, in which the surface oxygen vacancies close to the Ir single atoms activate molecular O 2 with a negligible barrier and the rate-limiting step is the gas-phase CO directly attacking the O-Ir species that is modulated by a CO ligand. Furthermore, first-principles X-ray absorption near-edge spectra of reaction intermediates along with in situ/operando X-ray absorption spectroscopy (XAS) suggest that Ir single atoms adsorb primarily on the step sites of MgAl 2 O 4 . However, microkinetic modeling predicts that a higher activity can be attained on the equally stable Mg-site, maximizing the population of which in catalyst synthesis might prove fruitful in future studies. Electronic structure analysis indicates that the CO ligand increases the reactivity of adsorbed oxygen atoms bound to Ir single atoms by increasing the antibonding nature of the O-Ir bond.

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

Tuning the bandgap of Cs 2 AgBiBr 6 through dilute tin alloying

The promise of lead halide hybrid perovskites for optoelectronic applications makes finding less-toxic alternatives a priority. The double perovskite Cs 2 AgBiBr 6 ( 1 ) represents one such alternative, offering long carrier lifetimes and greater stability under ambient conditions. However, the large and indirect 1.95 eV bandgap hinders its potential as a solar absorber. Here we report that alloying crystals of 1 with up to 1 atom% Sn results in a bandgap reduction of up to ca. 0.5 eV while maintaining low toxicity. Crystals can be alloyed with up to 1 atom% Sn and the predominant substitution pathway appears to be a ~2 : 1 substitution of Sn 2+ and Sn 4+ for Ag + and Bi 3+ , respectively, with Ag + vacancies providing charge compensation. Spincoated films of 1 accommodate a higher Sn loading, up to 4 atom% Sn, where we see mostly Sn 2+ substitution for both Ag + and Bi 3+ . Density functional theory (DFT) calculations ascribe the bandgap redshift to the introduction of Sn impurity bands below the conduction band minimum of the host lattice. Using optical absorption spectroscopy, photothermal deflection spectroscopy, X-ray absorption spectroscopy, 119 Sn NMR, redox titration, single-crystal and powder X-ray diffraction, multiple elemental analysis and imaging techniques, and DFT calculations, we provide a detailed analysis of the Sn content and oxidation state, dominant substitution sites, and charge-compensating defects in Sn-alloyed Cs 2 AgBiBr 6 ( 1:Sn ) crystals and films. An understanding of heterovalent alloying in halide double perovskites opens the door to a wider breadth of potential alloying agents for manipulating their band structures in a predictable manner.

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

Pursuing Heteroleptic Ligand Design Principles for Photoactive Fe Complexes with Ultrafast X-ray Emission and Variable-Temperature Optical Spectroscopies

Understanding the key parameters that govern the photophysical and photochemical properties of transition metal complexes is essential for the development of efficient photosensitizers for photocatalytic applications. Achieving this objective necessitates clear and detailed investigations of their electronic excited states, for which time-resolved metal Kβ X-ray emission spectroscopy (XES) has proven highly effective. Here, we present a time-resolved Fe Kβ XES study of a heteroleptic Fe(II) polypyridyl carbene complex, [Fe(phen) 2 (C 4 H 10 N 4 )] 2+ (1; phen = 1,10-phenanthroline), utilizing both the valence-to-core and Kβ mainline spectral regions, complemented by variable-temperature transient optical absorption (VT-TA) spectroscopy. Detailed analysis of the time-resolved Kβ XES data, supported by density functional theory (DFT) calculations and an Eyring analysis of the VT-TA data, reveals parallel excited state relaxation dynamics that support an assignment of the long-lived excited state to a triplet metal-centered state. Placing these results in the context of prior studies of heteroleptic Fe(II) polypyridyl cyanide complexes motivated a series of DFT calculations to investigate the effects of ligand structural flexibility and arrangement. These calculations reinforce the experimentally derived conclusion that constraining structural flexibility with multidentate ligands significantly impacts the excited state relaxation dynamics. Furthermore, our study emphasizes that the arrangement of strong field ligands in heteroleptic complexes substantially affects the energy of Jahn–Teller active triplet metal-centered states in low-spin d 6 metal complexes. Together, these findings provide synthetic design principles for extending metal-to-ligand charge transfer excited state lifetimes of heteroleptic Fe complexes.

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

Coexistence of Dirac fermion and charge density wave in the square-net-based semimetal LaAuSb 2

Here, we report a comprehensive study of magnetotransport properties, angle-resolved photoemission spectroscopy (ARPES), and density functional theory (DFT) calculations on self-flux grown LaAuSb 2 single crystals. Resistivity and Hall measurements reveal a charge density wave (CDW) transition at 77 K. MR and de Haas-Van Alphen (dHvA) measurements indicate that the transport properties of LaAuSb 2 are dominated by Dirac fermions that arise from Sb square nets. ARPES measurements and DFT calculations reveal an electronic structure with a common feature of the square-net-based topological semimetals, which is in good agreement with the magnetotransport properties. Our results indicate the coexistence of CDW and Dirac fermion in LaAuSb 2 , both of which are linked to the bands arising from the Sb-square net, suggesting that the square net could serve as a structural motif to explore various electronic orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

In‐Plane Anisotropy in van der Waals NiTeSe Ternary Alloy

Abstract The anisotropic properties of materials profoundly influence their electronic, magnetic, optical, and mechanical behaviors and are critical for a wide range of applications. In this study, the anisotropic characteristics of Ni‐based van der Waals materials, specifically NiTe 2 and its alloy NiTeSe, utilizing a combination of comprehensive scanning tunneling microscopy (STM), angle‐resolved photoemission spectroscopy (ARPES), and density functional theory (DFT) calculations, are explored. Unlike 1T‐NiTe 2 , which exhibits trigonal in‐plane symmetry, the substitution of Te with Se in NiTe 2 (resulting in the NiTeSe alloy) induces a pronounced in‐plane anisotropy. This anisotropy is clear in the STM topographs, which reveal a distinct linear order of charge distribution. Corroborating these observations, ARPES measurements and DFT calculations reveal an anisotropic Fermi surface centered at the point, which is notably elongated along the k y direction, leading to directional variations in in‐plane carrier velocities. Consequently, the Fermi velocity is highest along the k x direction where the linear charge distribution aligns in real space and is lowest along the k y direction. These findings offer valuable insights into the tunability of anisotropic properties in ternary transition metal dichalcogenide systems, highlighting their potential applications in the development of anisotropic electronic and optoelectronic devices.

Chemistry↗
DOE OSTI2021

AlCl 3 -Dosed Si(100)-2 × 1: Adsorbates, Chlorinated Al Chains, and Incorporated Al

The adsorption of AlCl 3 on Si(100) and the effect of annealing the AlCl 3 -dosed substrate were studied to reveal key surface processes for the development of atomic-precision, acceptor-doping techniques. This investigation was performed via scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. At room temperature, AlCl 3 readily adsorbed to the Si substrate dimers and dissociated to form a variety of species. Annealing the AlCl 3 -dosed substrate at temperatures below 450 °C produced unique chlorinated aluminum chains (CACs) elongated along the Si(100) dimer row direction. An atomic model for the chains is proposed with supporting DFT calculations. Al was incorporated into the Si substrate upon annealing at 450 °C and above, and Cl desorption was observed for temperatures beyond 450 °C. Al-incorporated samples were encapsulated in Si and characterized by secondary ion mass spectrometry (SIMS) depth profiling to quantify the Al atom concentration, which was found to be in excess of 10 20 cm –3 across a ~2.7 nm-thick δ-doped region. Here, the Al concentration achieved here and the processing parameters utilized promote AlCl 3 as a viable gaseous precursor for novel acceptor-doped Si materials and devices for quantum computing.

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

Phase Stability in U-6Nb Alloy Doped with Ti from the First Principles Theory

First-principles calculations within the density-functional-theory (DFT) approach are conducted in order to explore and explain the effect of small amounts of titanium on phase stability in the U-6Nb alloy. During rapid quenching from high to room temperature, metastable phases α' (orthorhombic), α" (monoclinic), and γ 0 (tetragonal) can form, depending on Nb concentration. Important mechanical properties depend on the crystal structure and, therefore, an understanding of the effect of impurities on phase stability is essential. Insights on this issue are obtained from quantum-mechanical DFT calculations. The DFT framework does not rely on any material-specific assumptions and is therefore ideal for an unbiased investigation of the U-Nb system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗
DOE OSTI2025

Stacking polymorphism of PtSe 2 : its implication to layer-dependent metal-insulator transitions

Using diffusion Monte Carlo (DMC) and density functional theory (DFT) calculations, we examine the structural stability and interlayer binding properties of PtSe 2 . Our DMC study reveals that AA and AB-r bilayer stacking modes are nearly degenerate, highlighting the significant role of interlayer hybridization in offsetting the energy cost due to larger interlayer separations in the AB-r mode. Additionally, our DMC-benchmarked DFT calculations with the r 2 SCAN+rVV10 functional uncover pronounced stacking polymorphism in few-layer PtSe 2 , driven by degenerate AA and AB-r interfaces, which leads to substantial band gap variations across different stacking configurations. This polymorphism, along with selenium vacancies, influences a layer-dependent metal-insulator transition observed in few-layer PtSe 2 . Our findings emphasize the importance of both van der Waals interactions and interlayer hybridization in determining the phase stability and electronic properties of PtSe 2 , advancing our understanding of its fundamental properties and refining theoretical models for practical applications in nanoelectronic devices.

Electronic materials↗
DOE OSTI2021

Intrinsically Honeycomb-Patterned Hydrogenated Graphene

Since the advent of graphene ushered the era of 2D materials, many forms of hydrogenated graphene have been reported, exhibiting diverse properties ranging from a tunable bandgap to ferromagnetic ordering. Patterned hydrogenated graphene with micron-scale patterns has been fabricated by lithographic means. We report successful millimeter-scale synthesis of an intrinsically honeycomb-patterned form of hydrogenated graphene on Ru(0001) by epitaxial growth followed by hydrogenation is reported. Combining scanning tunneling microscopy observations with density-functional-theory (DFT) calculations, it is revealed that an atomic-hydrogen layer intercalates between graphene and Ru(0001). The result is a hydrogen honeycomb structure that serves as a template for the final hydrogenation, which converts the graphene into graphane only over the template, yielding honeycomb-patterned hydrogenated graphene (HPHG). In effect, HPHG is a form of patterned graphane. DFT calculations find that the unhydrogenated graphene regions embedded in the patterned graphane exhibit spin-polarized edge states. This type of growth mechanism provides a new pathway for the fabrication of intrinsically patterned graphene-based materials.

2D material↗
DOE OSTI2024

Mixed Valence {Ni 2+ Ni 1+ } Clusters as Models of Acetyl Coenzyme A Synthase Intermediates

Acetyl coenzyme A synthase (ACS) catalyzes the formation and deconstruction of the key biological metabolite, acetyl coenzyme A (acetyl-CoA). The active site of ACS features a {NiNi} cluster bridged to a [Fe4S4] n+ cubane known as the A-cluster. The mechanism by which the A-cluster functions is debated, with few model complexes able to replicate the oxidation states, coordination features, or reactivity proposed in the catalytic cycle. In this work, we isolate the first bimetallic models of two hypothesized intermediates on the paramagnetic pathway of the ACS function. The heteroligated {Ni 2+ Ni 1+ } cluster, [K(12-crown-4) 2 ][1], effectively replicates the coordination number and oxidation state of the proposed “A red ” state of the A-cluster. Addition of carbon monoxide to [1] - allows for isolation of a dinuclear {Ni 2+ Ni 1+ (CO)} complex, [K(12-crown-2) n ][2] (n = 1–2), which bears similarity to the “A NiFeC ” enzyme intermediate. Structural and electronic properties of each cluster are elucidated by X-ray diffraction, nuclear magnetic resonance, cyclic voltammetry, and UV/vis and electron paramagnetic resonance spectroscopies, which are supplemented by density functional theory (DFT) calculations. Calculations indicate that the pseudo-T-shaped geometry of the three-coordinate nickel in [1] – is more stable than the Y-conformation by 22 kcal mol –1 , and that binding of CO to Ni 1+ is barrierless and exergonic by 6 kcal mol –1 . UV/vis absorption spectroscopy on [2] - in conjunction with time-dependent DFT calculations indicates that the square-planar nickel site is involved in electron transfer to the CO π*-orbital. Further, we demonstrate that [2] - promotes thioester synthesis in a reaction analogous to the production of acetyl coenzyme A by ACS.

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

Hydrogen trapping energetics at BCC iron-helium interfaces

Density functional theory (DFT) calculations have been performed to assess the trapping and segregation strength of hydrogen (H) to helium (He) interfaces in BCC (Fe) as a function of He density and surface orientation. The He density ranges from 0 to 2 He/V which corresponds to equilibrium bubbles that are expected to form in the structural components in a fusion reactor environment. The DFT calculations consist of a slab of BCC metal and an initial lattice of FCC He. Here, the report on the binding energy of H to these interfaces as well as the migration barriers into and away from the surface, which together provide information on the de-trapping energy for H at He bubbles. It was found that the binding energy of H to He-Fe interfaces decreases with increasing He density. Migration barriers are also modified due to the He density, which are sensitive to the surface orientation as well as the diffusion pathway. These results provide valuable first-principles energetics that are necessary for higher-order, mesoscale models, and provide insight into the extent to which He bubbles may trap tritium in fusion structural materials.

42 ENGINEERING↗
DOE OSTI2021

Low barrier height in a ZnO nanorods/NbSe 2 heterostructure prepared by van der Waals epitaxy

Two-dimensional (2D) materials as contacts for semiconductor devices have attracted much attention due to minimizing Fermi level pinning. Schottky–Mott physics has been widely employed to design 2D material-based electrodes and to elucidate their contact behavior. In this study, we revealed that charge transfer across a 2D/semiconductor heterointerface and materials characteristics besides work function should be accounted for in fabrication of electrodes based on 2D materials. Our density functional theory (DFT) calculations predicted that charge transfer between ZnO and NbSe 2 lowers the barrier height at the heterojunction and that conductive surface states of ZnO provide an additional conduction channel in the ZnO/NbSe 2 heterostructures. Crystalline ZnO/NbSe 2 heterostructures were prepared by the hydrothermal method. Electrical characterizations of the ZnO/NbSe 2 heterostructures showed Ohmic-like behavior as predicted by the DFT calculations, opposed to the prediction based on the Schottky–Mott model.

2D materials↗
DOE OSTI2021

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

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

Effects of Sn addition on precipitation of Al–Zn–Mg alloy at early stage of natural aging

The effect of Sn addition on the natural aging behavior of Al–Zn–Mg alloy was investigated using atom probe tomography (APT) experiments and density functional theory (DFT) calculations. The results demonstrate that Sn retards the natural aging behavior by facilitating the preferential formation of Mg–Sn clusters, which suppress the formation of Guinier-Preston zones (I) during the earliest aging stage. Quantitative solute analysis conducted by APT has suggested increasing amount of Mg–Sn clusters in Sn-contained alloy due to their high binding energy, supported by DFT calculations on the Sn.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Di- and Tetrameric Molybdenum Sulfide Clusters Activate and Stabilize Dihydrogen as Hydrides

NaY zeolite-encapsulated dimeric (Mo 2 S 4 ) and tetrameric (Mo 4 S 4 ) molybdenum sulfide clusters stabilize hydrogen as hydride binding to Mo atoms. Density functional theory (DFT) calculations and adsorption measurements suggest that stabilization of hydrogen as sulfhydryl (SH) groups, as typical for layered MoS 2 , is thermodynamically disfavored. Competitive adsorption of H 2 and ethene on Mo was probed by quantifying adsorbed CO on partly hydrogen and/or ethene covered samples with IR spectroscopy. During hydrogenation, experiment and theory suggest that Mo is covered predominately with ethene and sparsely with hydride. DFT calculations further predict that under reaction conditions, each Mo x S y cluster can activate only one H 2 , suggesting that the entire cluster (irrespective of its nuclearity) acts as one active site for hydrogenation. The nearly identical turnover frequencies (24.7 ± 3.3 mol ethane ·h -1 ·mol cluster -1 ), apparent activation energies (31-32 kJ·mol -1 ), and reaction orders (~0.5 in ethene and ~1.0 in H 2 ) show that the active sites in both clusters are catalytically indistinguishable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗