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

Photophysics of graphene quantum dot assemblies with axially coordinated cobaloxime catalysts

We report a study of chromophore-catalyst assemblies composed of light harvesting hexabenzocoronene (HBC) chromophores axially coordinated to two cobaloxime complexes. The chromophore-catalyst assemblies were prepared using bottom-up synthetic methodology and characterized using solid-state NMR, IR, and x-ray absorption spectroscopy. Detailed steady-state and time-resolved laser spectroscopy was utilized to identify the photophysical properties of the assemblies, coupled with time-dependent DFT calculations to characterize the relevant excited states. The HBC chromophores tend to assemble into aggregates that exhibit high exciton diffusion length (D = 18.5 molecule 2 /ps), indicating that over 50 chromophores can be sampled within their excited state lifetime. We find that the axial coordination of cobaloximes leads to a significant reduction in the excited state lifetime of the HBC moiety, and this finding was discussed in terms of possible electron and energy transfer pathways. By comparing the experimental quenching rate constant (1.0 × 10 9 s -1 ) with the rate constant estimates for Marcus electron transfer (5.7 × 10 8 s -1 ) and Förster/Dexter energy transfers (8.1 × 10 6 s -1 and 1.0 × 10 10 s -1 ), we conclude that both Dexter energy and Marcus electron transfer process are possible deactivation pathways in CoQD-A. No charge transfer or energy transfer intermediate was detected in transient absorption spectroscopy, indicating fast, subpicosecond return to the ground state. Finally, these results provide important insights into the factors that control the photophysical properties of photocatalytic chromophore-catalyst assemblies.

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

Self-compensation in heavily Ge doped AlGaN: A comparison to Si doping

Self-compensation in Ge- and Si-doped Al 0.3 Ga 0.7 N has been investigated in terms of the formation of III vacancy and donor-vacancy complexes. Both Ge- and Si-doped AlGaN layers showed a compensation knee behavior with impurity compensation (low doping regime), compensation plateau (medium doping regime), and self-compensation (high doping regime). A maximum free carrier concentration of 4–5 × 10 19 cm –3 was obtained by Ge doping, whereas Si doping resulted in only half of that value, ~2 × 10 19 cm –3 . A DFT calculation with the grand canonical thermodynamics model was developed to support the hypothesis that the difference in self-compensation arises from the difference in the formation energies of the V III -n•donor complexes relative to their onsite configurations. The model suggested that the V III -2•donor and V III -3•donor complexes were responsible for self-compensation for both Ge- and Si-doped AlGaN. However, a lower free carrier concentration in Si-doped samples was due to a high V III -3•Si concentration, resulting from a lower energy of formation of V III -3•Si.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Atomic scale investigation of aluminum incorporation, defects, and phase stability in β -(Al x Ga 1-x ) 2 O 3 films

The development of novel ultra-wide bandgap (UWBG) materials requires precise understanding of the atomic level structural origins that give rise to their important properties. We study the aluminum atom incorporation, defect formation, and their relationships with phase stability in β-(Al x Ga 1-x ) 2 O 3 films, a promising candidate for UWBG applications, to explain atomic scale structural characteristics and properties using a combination of quantitative scanning transmission electron microscopy (STEM) and density functional theory (DFT). Our STEM analysis indicates that ~54% of the incorporated Al substitutes on the octahedrally coordinated Ga 2 site in a series of films grown with different techniques and alloy concentrations. DFT calculations show that, while Al energetically prefers the octahedral site, surface reconstructions and kinetic limitations during the epitaxial growth are responsible for Al occupying both octahedral and tetrahedral sites in (Al x Ga 1-x ) 2 O 3 , ultimately limiting the stability of the β-phase at x < ~50%. Local heterogeneity of composition results in the formation of a planar defect, affecting the stability of the β-phase. The similarity of such inclusions to the metastable γ-phase is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structure prediction of epitaxial inorganic interfaces by lattice and surface matching with Ogre

We present a new version of the Ogre open source Python package with the capability to perform structure prediction of epitaxial inorganic interfaces by lattice and surface matching. In the lattice matching step, a scan over combinations of substrate and film Miller indices is performed to identify the domain-matched interfaces with the lowest mismatch. Subsequently, surface matching is conducted by Bayesian optimization to find the optimal interfacial distance and in-plane registry between the substrate and the film. For the objective function, a geometric score function is proposed based on the overlap and empty space between atomic spheres at the interface. The score function reproduces the results of density functional theory (DFT) at a fraction of the computational cost. The optimized interfaces are pre-ranked using a score function based on the similarity of the atomic environment at the interface to the bulk environment. Final ranking of the top candidate structures is performed with DFT. Ogre streamlines DFT calculations of interface energies and electronic properties by automating the construction of interface models. The application of Ogre is demonstrated for two interfaces of interest for quantum computing and spintronics, Al/InAs and Fe/InSb.

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

Intrinsic magnetism in KTaO3 heterostructures

There has been intense recent interest in the two-dimensional electron gases (2DEGs) that form at the surfaces and interfaces of KTaO3 (KTO), with the discovery of superconductivity at temperatures significantly higher than those of similar 2DEGs based on SrTiO3 (STO). Like STO heterostructures, these KTO 2DEGs are formed by depositing an overlayer on top of appropriately prepared KTO surfaces. Some of these overlayers are magnetic, and the resulting 2DEGs show signatures of this magnetism, including hysteresis in the magnetoresistance (MR). Here, we show that KTO 2DEGs fabricated by depositing AlOx on top of KTO also show hysteretic MR, indicative of long-range magnetic order, even though the samples nominally contain no intrinsic magnetic elements. The hysteresis appears in both the transverse and longitudinal resistance in magnetic fields both perpendicular to and in the plane of the 2DEG. The hysteretic MR has different characteristic fields and shapes for surfaces of different crystal orientations and vanishes above a few Kelvin. Density functional theory (DFT) calculations indicate that the magnetism likely arises from Ta4+ local moments created in the presence of oxygen vacancies.

Physics↗
DOE OSTI2024

An experimental and computational view of the photoionization of diol–water clusters

In the interstellar medium, diols and other prebiotic molecules adsorb onto icy mantles surrounding dust grains. Water in the ice may affect the reactivity and photoionization of these diols. Ethylene glycol (EG), 1,2-propylene glycol, and 1,3-propylene glycol clusters with water clusters were used as a proxy to study these interactions. The diol–water clusters were generated in a continuous supersonic molecular beam, photoionized by synchrotron-based vacuum ultraviolet light from the Advanced Light Source, and subsequently detected by reflectron time-of-flight mass spectrometry. The appearance energies for the detected clusters were determined from the mass spectra, collected at increasing photon energy. Clusters of both diol fragments and unfragmented diols with water were detected. The lowest energy geometry optimized conformers for the observed EG–water clusters and EG fragment–water clusters have been visualized using density functional theory (DFT), providing insight into hydrogen bonding networks and how these affect fragmentation and appearance energy. As the number of water molecules clustered around EG fragments (m/z 31 and 32) increased, the appearance energy for the cluster decreased, indicating a stabilization by water. This trend was supported by DFT calculations. Fragment clusters from 1,2-propylene glycol exhibited a similar trend, but with a smaller energy decrease, and no trend was observed from 1,3-propylene glycol. In conclusion, we discuss and suggest that the reactivity and photoionization of diols in the presence of water depend on the size of the diol, the location of the hydroxyl group, and the number of waters clustered around the diol.

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

Improving gas adsorption modeling for MOFs by local calibration of Hubbard U parameters

While computational screening with density functional theory (DFT) is frequently employed for the screening of metal–organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule–MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs.

Chemistry↗
DOE OSTI2024

Improving the precision of forces in real-space pseudopotential density functional theory

The high-order finite difference real-space pseudopotential density functional theory (DFT) approach is a valuable method for large-scale, massively parallel DFT calculations. A significant challenge in the approach is the oscillating “egg-box” error introduced by aliasing associated with a coarse grid spacing. To address this issue while minimizing computational cost, we developed a finite difference interpolation (FDI) scheme [Roller et al., J. Chem. Theory Comput. 19, 3889 (2023)] as a means of exploiting the high resolution of the pseudopotential to reduce egg-box effects systematically. Here, we show an implementation of this method in the PARSEC code and examine the practical utility of the combination of FDI with additional methods for improving force precision and/or reducing its computational cost, including orbital-based forces, compensating charges (namely, adding and subtracting a judiciously chosen charge density such that the total density is unaltered), and a modified spatial domain in which the real-space grid is defined. Using selected small molecules, as well as metallic Li, as test cases, we show that a combination of all four aspects leads to a significant reduction in computational cost while retaining a high level of precision that supports accurate structures and vibrational spectra, as well as stable and accurate molecular dynamics runs.

Chemistry↗
DOE OSTI2024

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗
DOE OSTI2025

Machine learning interatomic potential for predicting the thermal properties of uranium nitride

We present a combined computational and experimental investigation of the thermal properties of uranium nitride (UN), focusing on the development of a machine learning interatomic potential (MLIP) using the moment tensor potential framework. The MLIP was trained on density functional theory (DFT) data and validated against various quantities including energies, forces, elastic constants, phonon dispersion, and defect formation energies, achieving excellent agreement with DFT calculations, prior experimental results, and our thermal conductivity measurement. The potential was then employed in molecular dynamics simulations to predict key thermal properties such as melting point, thermal expansion, specific heat, and lattice thermal conductivity. To further assess model accuracy, we fabricated a UN sample and performed new thermal conductivity measurements representative of single-crystal properties, which showed strong agreement with the MLIP predictions. This work confirms the reliability and predictive capability of the developed potential for determining the thermal properties of UN.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Shear strength of a refractory high entropy alloy MoNbTaVW under high pressure

Radial X-ray diffraction (R-XRD) was performed in situ using a Panoramic Diamond Anvil Cell on the refractory high entropy alloy MoNbTaVW. The lattice parameters were determined through a Le Bail fit using a body centered cubic (BCC) lattice of symmetry $Im$$\overline{3}$$m$ (international space group number 229). Upper and lower bounds to the shear strength were determined up to 80 GPa nonhydrostatically using copper as a pressure standard. The equation of state was derived at the magic angle ψ = 54.7° and yielded a bulk modulus of K 0 = 220.8 ± 1.65 GPa. The experimental lattice parameters and bulk modulus match closely with corresponding density functional theory (DFT) calculations. The BCC phase remains stable up to the highest pressure of 80 GPa studied and is shown to be elastically anisotropic. The shear strength was found to saturate around 70 GPa with a value of τ = 1.75 GPa, and the shear moduli are presented in different limits of iso-strain and iso-stress.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Superconductivity enhancement in phase-engineered molybdenum carbide/disulfide vertical heterostructures

Stacking layers of atomically thin transition-metal carbides and two-dimensional (2D) semiconducting transition-metal dichalcogenides, could lead to nontrivial superconductivity and other unprecedented phenomena yet to be studied. In this work, superconducting α-phase thin molybdenum carbide flakes were first synthesized, and a subsequent sulfurization treatment induced the formation of vertical heterolayer systems consisting of different phases of molybdenum carbide—ranging from α to γ' and γ phases—in conjunction with molybdenum sulfide layers. These transition-metal carbide/disulfide heterostructures exhibited critical superconducting temperatures as high as 6 K, higher than that of the starting single-phased α-Mo 2 C (4 K). We analyzed possible interface configurations to explain the observed moiré patterns resulting from the vertical heterostacks. Our density-functional theory (DFT) calculations indicate that epitaxial strain and moiré patterns lead to a higher interfacial density of states, which favors superconductivity. Such engineered heterostructures might allow the coupling of superconductivity to the topologically nontrivial surface states featured by transition-metal carbide phases composing these heterostructures potentially leading to unconventional superconductivity. Moreover, we envisage that our approach could also be generalized to other metal carbide and nitride systems that could exhibit high-temperature superconductivity.

heterostructure↗
DOE OSTI2020

Enhanced and stabilized hydrogen production from methanol by ultrasmall Ni nanoclusters immobilized on defect-rich h-BN nanosheets

Employing liquid organic hydrogen carriers (LOHCs) to transport hydrogen to where it can be utilized relies on methods of efficient chemical dehydrogenation to access this fuel. Therefore, developing effective strategies to optimize the catalytic performance of cheap transition metal-based catalysts in terms of activity and stability for dehydrogenation of LOHCs is a critical challenge. Here, we report the design and synthesis of ultrasmall nickel nanoclusters (∼1.5 nm) deposited on defect-rich boron nitride (BN) nanosheet (Ni/BN) catalysts with higher methanol dehydrogenation activity and selectivity, and greater stability than that of some other transition-metal based catalysts. The interface of the two-dimensional (2D) BN with the metal nanoparticles plays a strong role both in guiding the nucleation and growth of the catalytically active ultrasmall Ni nanoclusters, and further in stabilizing these nanoscale Ni catalysts against poisoning by interactions with the BN substrate. We provide detailed spectroscopy characterizations and density functional theory (DFT) calculations to reveal the origin of the high productivity, high selectivity, and high durability exhibited with the Ni/BN nanocatalyst and elucidate its correlation with nanocluster size and support–nanocluster interactions. This study provides insight into the role that the support material can have both regarding the size control of nanoclusters through immobilization during the nanocluster formation and also during the active catalytic process; this twofold set of insights is significant in advancing the understanding the bottom-up design of high-performance, durable catalytic systems for various catalysis needs.

08 HYDROGEN↗
DOE OSTI2023

A first principles analysis of potential-dependent structural evolution of active sites in Fe-N-C catalysts

Fe-N-C (iron–nitrogen–carbon) electrocatalysts have emerged as potential alternatives to precious metal-based materials for the oxygen reduction reaction (ORR). However, the structure of these materials under electrochemical conditions is not well understood, and their poor stability in acidic environments poses a formidable challenge for successful adoption in commercial fuel cells. To provide molecular-level insights into these complex phenomena, we combine periodic density functional theory (DFT) calculations, exhaustive treatment of coadsorption effects for ORR reaction intermediates, including O and OH, and comprehensive analysis of solvation stabilization effects to construct voltage-dependent ab initio thermodynamic phase diagrams that describe the in situ structure of the active sites. These structures are further linked to activity and stability descriptors that can be compared with experimental parameters such as the half-wave potential for ORR and the onset potential for carbon corrosion and CO 2 evolution. The results indicate that pyridinic Fe sites at zigzag carbon edges, as well as other edge sites, exhibit high activity for ORR compared to sites in the bulk. However, edges neighboring the active sites are prone to instability via overoxidation and consequent site loss. The results suggest that it could be beneficial to synthesize Fe-N-C catalysts with small sizes and large perimeter edge lengths to enhance ORR activity, while voltage fluctuations should be limited during fuel cell operation to prevent carbon corrosion of overoxidized edges.

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

Secondary structure determines electron transport in peptides

Proteins play a key role in biological electron transport, but the structure–function relationships governing the electronic properties of peptides are not fully understood. Despite recent progress, understanding the link between peptide conformational flexibility, hierarchical structures, and electron transport pathways has been challenging. Here, we use single-molecule experiments, molecular dynamics (MD) simulations, nonequilibrium Green’s function-density functional theory (NEGF-DFT), and unsupervised machine learning to understand the role of secondary structure on electron transport in peptides. Our results reveal a two-state molecular conductance behavior for peptides across several different amino acid sequences. MD simulations and Gaussian mixture modeling are used to show that this two-state molecular conductance behavior arises due to the conformational flexibility of peptide backbones, with a high-conductance state arising due to a more defined secondary structure (beta turn or 3 10 helices) and a low-conductance state occurring for extended peptide structures. These results highlight the importance of helical conformations on electron transport in peptides. Conformer selection for the peptide structures is rationalized using principal component analysis of intramolecular hydrogen bonding distances along peptide backbones. Molecular conformations from MD simulations are used to model charge transport in NEGF-DFT calculations, and the results are in reasonable qualitative agreement with experiments. Projected density of states calculations and molecular orbital visualizations are further used to understand the role of amino acid side chains on transport. Overall, our results show that secondary structure plays a key role in electron transport in peptides, which provides broad avenues for understanding the electronic properties of proteins.

Science & Technology - Other Topics↗
DOE OSTI2024

Modeling prebiotic chemistries with quantum accuracy at classical costs

Molecular Dynamics (MD) simulations using classical force-fields are commonly employed in numerous scientific investigations. However, many natural processes involve bond breaking and quantum forces. This complexity is compounded by the presence of multiple competing length and timescales. For example, accurately modeling the thermodynamics and dynamics of a chemical reaction requires accounting for the concerted movements of numerous solvent molecules and ions with their own fast or slow timescales. While widely used static Density Functional Theory (DFT) calculations at 0 temperature can be beneficial for such investigations, they do not account for dynamics, and lack precision in describing the molecular environments. They particularly fail at correct, rigorous treatments of finite-temperature fluctuations, and thus generalization to experimentally relevant conditions. In PNAS Benayad et al develop a scalable, generalizable approach for designing Neural Network Potentials (NNPs) that can handle chemical reactivity in solvated systems with quantum accuracy at classical costs. Specifically, they study phosphoester bond formation and rupture, which is fundamentally relevant to the Phosphorus-Oxygen bond formation central to life, and especially for the RNA world hypothesis. The framework developed here has the potential to generalize to different chemical reactions of energy and biological relevance.

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

Electronic structure and effective mass of pristine and Cl-doped CsPbBr 3

Organic–inorganic lead halide perovskites (LHPs) have attracted great interest owing to their outstanding optoelectronic properties. Typically, the underlying electronic structure would determinate the physical properties of materials. But as for now, limited studies have been done to reveal the underlying electronic structure of this material system, comparing to the huge amount of investigations on the material synthesis. The effective mass of the valance band is one of the most important physical parameters which plays a dominant role in charge transport and photovoltaic phenomena. In pristine CsPbBr 3 , the Fröhlich polarons associated with the Pb–Br stretching modes are proposed to be responsible for the effective mass renormalization. In this regard, it would be very interesting to explore the electronic structure in doped LHPs. Here, we report high-resolution angle-resolved photoemission spectroscopy (ARPES) studies on both pristine and Cl-doped CsPbBr 3 . The experimental band dispersions are extracted from ARPES spectra along both $\overline{\Gamma}$–$\overline{M}$–$\overline{\Gamma}$ and $\overline{X}$–$\overline{M}$–$\overline{X}$ high symmetry directions. DFT calculations are performed and directly compared with the ARPES data. Our results have revealed the band structure of Cl-doped CsPbBr 3 for the first time, which have also unveiled the effective mass renormalization in the Cl-doped CsPbBr 3 compound. Doping dependent measurements indicate that the chlorine doping could moderately tune the renormalization strength. These results will help understand the physical properties of LHPs as a function of doping.

36 MATERIALS SCIENCE↗