Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Band structure, superconductivity, and polytypism in AuSn 4

The orthorhombic compound AuSn 4 is compositionally similar to the Dirac node arc semimetal PtSn 4 . AuSn 4 is, contrary to PtSn 4 , superconducting with a critical temperature of T c = 2.35 K. Recent measurements present indications for quasi-two-dimensional superconducting behavior in AuSn 4 . Here we present measurements of the superconducting density of states and the band structure of AuSn 4 through scanning tunneling microscopy and angular resolved photoemission spectroscopy (ARPES). The superconducting gap values in different portions of the Fermi surface are spread around Δ 0 = 0.4 meV, which is close to but somewhat larger than Δ = 1.76k B T c expected from BCS theory. We observe superconducting features in the tunneling conductance at the surface up to temperatures about 20% larger than bulk T c . The band structure calculated with density functional theory follows well the results of ARPES. The crystal structure presents two possible stackings of Sn layers, giving two nearly degenerate polytypes. This makes AuSn 4 a rather unique case with a three-dimensional electronic band structure but properties ressembling those of low-dimensional layered compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Enhanced Guest@MOF Interaction via Stepwise Thermal Annealing: TCNQ@Cu 3 (BTC) 2

Confinement of guest molecules in porous materials such as metal organic frameworks (MOFs) promises to deliver emergent properties separate from those of the individual components. Understanding the confinement mechanism is therefore important for the development of new synthesis routes that adjust MOF properties for specific applications. In this work, we developed a new synthetic method to confine guest molecules into MOF pores through a stepwise thermal annealing process, wherein the confinement of 7,7,8,8-tetracyanoquinodimethane (TCNQ) guest molecules into Cu 3 (BTC) 2 (BTC = benzene-1,3,5-tricarboxylic acid) MOF host is used as an example of how novel materials can be created with new physical properties. The stepwise thermal annealing process includes 1) an activation process of pristine Cu 3 (BTC) 2 MOF to maximizes the TCNQ guest loading in the MOF host by effectively removing the residual solvents and 2) post-annealing of the TCNQ infiltrated MOF to enhances the interaction of the confined guest molecules with the MOF host. Obtained experimental results based on thermogravimetric analysis, N 2 gas adsorption, electron microscopy, X-ray diffraction and infrared absorption, combined with density functional theory calculations provide evidence that the use of a stepwise thermal annealing process yields enhancements in the guest loading, packing and interaction between the TCNQ guest and the MOF host. The new hybrid TCNQ@Cu 3 (BTC) 2 system is stable and shows no significant signs of structural degradation even after submersion in water. This is due to the presence of significantly stronger interactions of TCNQ with the frame-work metal ions compared to those of the water molecules competing for the same framework binding sites. It was also found that TCNQ@Cu 3 (BTC) 2 system maintains a significant CO 2 and CH 4 adsorption potential compared to the pristine MOF. The synthetic route developed in this work yields novel guest@MOF hybrid systems that will be useful for many MOF-based applications such as gas separations and chemical sensors performed under humid conditions.

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

Influence of surface chemistry on Li nucleation energetics on graphene-based surfaces

Lithium metal is a promising high-capacity anode material for solid-state batteries, but it typically suffers from poor cyclability. Carbon scaffold hosts have the potential to improve this performance due to their high electronic conductivity and large surface area, which facilitates lithium-ion adsorption and desorption. Scaffold surface chemistry is known to significantly influence performance outcomes, but the details of these interactions are not fully understood. Here, this study employs first-principles simulations to explore lithium transport and nucleation on graphene anodes with various surface chemistries. Using enhanced sampling techniques, ab initio molecular dynamics, and density functional theory calculations, we find that although surface chemistry has a minimal impact on lithium interfacial transport, it influences surface nucleation significantly. Both heteroatom dopants and intrinsic defects lower the nucleation barrier, creating a more favorable environment for lithium nucleation compared to pristine graphene. In addition, our results reveal a complex interplay between surface lithium concentration, lithium transport, and nucleation kinetics. These findings highlight the potential of surface modifications to precisely control nucleation processes on carbon-based anodes and provide design guidance for reducing dendrite formation and improving the cycle life of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-throughput investigation of the formation of double spinels

Spinel compounds, with the general chemical formula AB 2 O 4 , are a wide class of materials, where A and B can be a variety of cations, providing this structure with a great deal of functional flexibility and giving rise to its considerable scientific interest. Recently, a spinel with the general formula ABB'O 4 has been predicted, increasing the possible usability of the spinels due to the higher cation diversity in the so-called double spinel structure. In this study, we use density functional theory calculations to predict if double spinels can be formed between experimentally synthesized normal and inverse single spinels. Our computations reveal that 49 double spinels have negative mixing enthalpies and are thus thermodynamically stable, with most of the stable compounds being formed from one of two distinct cation orderings. We show that the 17 different cations that form the different double spinels have a preferred site, tetrahedral or octahedral, except for Mn, Fe and Co which can occupy both sites interchangeably. We also study the relation between mixing enthalpies and cation-specific properties, as well as ways to classify the double spinels into distinct types and spinel groups depending on the cation ordering and cation oxidation states, respectively. Furthermore, by applying the Sure Independence Screening and Sparsifying Operator (SISSO) approach on the coordination-dependent ionic radii of the elemental constituents, we show that an interplay of local strain and electrostatic dominated terms can be used to separate the double spinels into distinct structural types depending on the cation order and their oxidation states.

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

Low-frequency Raman spectrum of 2D layered perovskites: Local atomistic motion or superlattice modes?

We report the low-frequency Raman spectrum (ω = 10 cm -1 –150 cm -1 ) of a wide variety of alkylammonium iodide based 2D lead halide perovskites (2D LHPs) as a function of A-site cation (MA = methylammonium and FA = formamidinium), octahedral layer thickness (n = 2–4), organic spacer chain length (butyl-, pentyl-, hexyl-), and sample temperature (T = 77 K–293 K). Using density functional theory calculations under the harmonic approximation for n = 2 BA:MAPbI, we assign several longitudinal/transverse optical phonon modes between 30 cm -1 and 100 cm -1 , the eigen displacements of which are analogous to that observed previously for octahedral twists/distortions in bulk MAPbI. Additionally, we propose an alternative assignment for low-frequency modes below this band (<30 cm -1 ) as zone-folded longitudinal acoustic phonons corresponding to the periodicity of the entire layered structure. Furthermore, we compare measured spectra to predictions of the Rytov elastic continuum model for zone-folded dispersion in layered structures. Our results are consistent across the various 2D LHPs studied herein, with energetic shifts of optical phonons corresponding to microscopic structural differences between materials and energetic shifts of acoustic phonons according to changes in the periodicity and elastic properties of the perovskite/organic subphases. This study highlights the importance of both the local atomic order and the superlattice structure on the vibrational properties of layered 2D materials.

2D materials↗
DOE OSTI2024

Coexistence of Interacting Charge Density Waves in a Layered Semiconductor

Coexisting orders are key features of strongly correlated materials and underlie many intriguing phenomena from unconventional superconductivity to topological orders. Here, we report the coexistence of two interacting charge-density-wave (CDW) orders in EuTe 4 , a layered crystal that has drawn considerable attention owing to its anomalous thermal hysteresis and a semiconducting CDW state despite the absence of perfect FS nesting. By accessing unoccupied conduction bands with time- and angle-resolved photoemission measurements, we find that mono- and bi-layers of Te in the unit cell host different CDWs that are associated with distinct energy gaps. The two gaps display dichotomous evolutions following photoexcitation, where the larger bilayer CDW gap exhibits less renormalization and faster recovery. Surprisingly, the CDW in the Te monolayer displays an additional momentum-dependent gap renormalization that cannot be captured by density-functional theory calculations. This phenomenon is attributed to interlayer interactions between the two CDW orders, which account for the semiconducting nature of the equilibrium state. Finally, our findings not only offer microscopic insights into the correlated ground state of EuTe 4 but also provide a general non-equilibrium approach to understand coexisting, layer-dependent orders in a complex system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Superstructures and magnetic order in heavily Cu-substituted ( Fe 1 – x Cu x ) 1 + y Te

Most iron-based superconductors exhibit stripe-type magnetism, characterized by the ordering vector Q = ($\frac{1}{2},\frac{1}{2}$). In contrast, Fe 1+y Te, the parent compound of the Fe 1+y Te 1–x Se x superconductors, exhibits double-stripe magnetic order associated with the ordering vector Q = ($\frac{1}{2},0$). Here, we use elastic neutron scattering to investigate heavily Cu-substituted (Fe 1–x Cu x ) 1+y Te compounds and reveal that (1) for x ≳ 0.4, short-range magnetic order emerges around the stripe-type vector at Q = ($\frac{1}{2}$ ± δ, $\frac{1}{2}$ ± δ, $\frac{1}{2}$) with δ ≈ 0.05; (2) the short-range magnetic order is associated with a superstructure modulation at Q = ($\frac{1}{3},\frac{1}{3},\frac{1}{2}$), with the magnetic correlation length shorter than that for the superstructure; and (3) for x ≳ 0.55, we observe an additional intergrown phase with higher Cu content, characterized by a superstructure modulation vector Q = ($\frac{1}{3},\frac{1}{3},0$) and magnetic peaks at Q = ($\frac{2}{3},\frac{1}{3},\frac{1}{2}$)/($\frac{1}{3},\frac{2}{3},\frac{1}{2}$). The positions of superstructure peaks suggest that relative to the tetragonal unit cell of Fe 1+y Te, heavy Cu substitution leads to Fe-Cu orderings that expand the unit cell by $\sqrt{2}$ × 3$\sqrt{2}$ times in the ab plane, corroborated by first-principles calculations that suggest the formation of spin chains and spin ladders. Finally, our findings show that stripe-type magnetism is common in magnetically diluted iron pnictides and chalcogenides, despite the varying associated atomic orderings

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Magnetic multilayer edges in Bernal-stacked hexagonal boron nitride

Single-layer h-BN is known to have edges with unique magnetism; however, in the commonly fabricated multilayer AA'-h-BN, edge relaxations occur that create interlayer bonds and eliminate the unpaired electrons at the edge. Recently, a robust method of growing the unconventional Bernal-stacked h-BN (AB-h-BN) has been reported. Here, we use theoretical approaches to investigate the nitrogen-terminated zigzag edges in AB-h-BN that can be formed in a controlled fashion using a high-energy electron beam. We find that these "open"edges remain intact in bilayer and multilayer AB-h-BN, enabling researchers potentially to investigate these edge states experimentally. We also investigate the thermodynamics of the spin configurations at the edge by constructing a lattice model that is based on parameters extracted from a set of first-principles calculations. We find that the edge spins in neighboring layers interact very weakly, resulting in a sequence of independent spin chains in multilayer samples. By solving this model using Monte Carlo simulations, we can determine nm-scale correlation lengths at liquid-N 2 temperatures and lower. At low temperatures, these edges may be utilized in magnetoresistance and spintronics applications.

1-dimensional spin chains↗
DOE OSTI2023

High-frequency phonons drive large phonon-drag thermopower in semiconductors at high carrier density

It has been well established that (i) the thermopower of semiconductors can be enhanced through a phe- nomenon known as the drag effect, and (ii) the drag enhancement involves only low-frequency acoustic phonons and benefits from low electron densities and low temperatures. Using first-principles calculations we show that large drag enhancements to the thermopower are possible at high carrier density even at room temperature and arise from high-frequency acoustic phonons. A fascinating example is cubic boron arsenide (BAs) for which the calculated room temperature drag enhancement of the thermopower exceeds an order of magnitude at a high hole density of 10 21 cm –3 . This remarkable behavior stems from the simultaneously weak phonon-phonon and phonon-hole scattering of the high-frequency phonons in BAs that become drag active at high carrier densities through electron-phonon interactions. Furthermore, this work advances our understanding of coupled electron-phonon nanoscale transport and introduces an unexpected paradigm for achieving large thermopowers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deep Understanding of Strong Metal Interface Confinement: A Journey of Pd/FeO x Catalysts

Tuning the atomic interface configuration of noble metals (NMs) and transition-metal oxides is an effective straightforward yet challenging strategy to modulate the activity and stability of heterogeneous catalysts. Herein, Pd supported on mesoporous Fe 2 O 3 with a high specific surface area was rationally designed and chosen to construct the Pd/iron oxide interface. As a versatile model, the physicochemical environments of Pd nanoparticles (NPs) could be precisely controlled by taming the reduction temperature. The experimental and density functional theory calculation results unveiled that the catalyst in the support–metal interface confinement (SMIC) state showed significantly enhanced catalytic activity and sintering resistance for CO oxidation. The constructed Fe sites at the interfaces between FeO x overlayers and Pd NPs not only provided additional coordinative unsaturated ferrous sites for the adsorption and activation of O 2 , thereby facilitating the activation efficiency of O 2 , but also impressively changed the reaction pathway of CO oxidation. As a result, the catalyst followed the Pd/Fe dual-site mechanism instead of the classical Mars–van Krevelen mechanism. For the catalyst in the strong metal–support interaction (SMSI) state, its catalytic activity was seriously suppressed because of the excessive encapsulation of the active Pd sites by FeO x overlayers. Thus, the present study therefore provides detailed insights into the SMIC and SMSI in ferric oxide-supported Pd catalysts, which could guide the preparation of highly efficient supported catalysts for practical applications.

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

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗
DOE OSTI2024

Accelerating the Discovery of New, Single Phase High Entropy Ceramics via Active Learning

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. Here we then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiple Reaction Pathways for the Oxygen Evolution Reaction May Contribute to IrO 2 (110)’s High Activity

Density functional theory calculations in conjunction with statistical mechanical arguments are performed on the rutile IrO 2 (110) facet in order to characterize multiple reaction pathways on the surface at the highest active limit (the stoichiometric surface with all metal sites available) and at the lowest active limit (the oxygen-terminated surface). Alternative pathways to the oxygen evolution reaction (OER) are found, with multiple pathways determined at each step of the four proton-coupled electron transfer reaction. Of particular interest is the detailed characterization of a co-adsorption pathway utilizing neighboring, adsorbed O, OH species in order to evolve oxygen; activation energies of this pathway are <0.5 eV and therefore easily surmountable at the high operating potentials of OER. We also determined that surface Ir atoms can potentially participate in deprotonating an OOH* intermediate; the activation energy to this is 0.67 eV on the oxygen-terminated surface. These theoretical findings explain in part the high activity present in iridium oxide catalysts and also provide insight into the mechanistic pathways available on metal oxide catalysts, which may require the concerted interaction of nearest neighbor co-adsorbates to produce chemicals of interest.

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

Phase transitions and potential ferroelectricity in noncentrosymmetric KNaNbOF 5

Increasing the number of known ferroelectrics requires expanding the chemical design space and mechanisms producing ferroelectricity. To that end, we examine the displacive, noncentrosymmetric-to-centrosymmetric phase transition in the oxyfluoride KNaNbOF 5 using ab initio calculations and Landau theory. Here we predict an intermediate Pnma phase occurs in the transition from the known low- and high-temperature phases, Pna2 1 and Cmcm, respectively, which we subsequently confirm using high-resolution and in situ 19 F magic-angle spinning nuclear magnetic resonance. Using the intermediate structure, we then calculate the monodomain polarization switching barrier in KNaNbOF 5 to be ~ 93 meV per formula unit, comparable to existing ferroelectrics. The reversal of the polarization is due to uncompensated antipolar displacements of oxygen and fluorine and does not require oxygen and fluorine site exchange as in other nonswitchable heteroanionic materials, which makes KNaNbOF 5 a viable oxyfluoride for a switchable electric polarization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanism of the Accelerated Water Formation Reaction under Interfacial Confinement

Subnanoscale spaces at the interface between weakly coupled thin films and their metal supports offer exciting opportunities for studying chemical reactions under confinement. Here, we investigated the kinetics of water formation (from hydrogen and chemisorbed oxygen) in the confined space at the interface between bilayer (BL) silica and a Ru(0001) support, compared to the reaction on the bare Ru(0001) surface. Ambient pressure X-ray photoelectron spectroscopy (AP-XPS) experiments were carried out at different temperatures at elevated pressures of H 2 to follow the reaction kinetics. The apparent activation energy at the BL-silica/Ru(0001) interface was found to be 0.38 eV lower than that on bare Ru(0001), consistent with a recent report by Prieto et al. (Angew. Chem., Int. Ed.2018, 57(28), 8749–8753) carried out at lower H 2 pressures using low-energy electron microscopy. Density functional theory calculations revealed that the rate-limiting step in the direct hydrogenation pathway on the Ru(0001) surface is the first hydrogen addition step (*H + *O ↔ *OH). The confinement at the BL-silica/Ru(0001) interface only marginally affects the energy barrier of the first hydrogen addition. Instead, it activates an alternative disproportionation reaction pathway (*H 2 O + *O ↔ 2*OH). On the bare Ru(0001) surface, the disproportionation pathway can only occur at cryogenic temperatures or under high water vapor pressures. However, the presence of the BL-silica increases the desorption barrier for water molecules at the interface. The increased residence time allows trapped water molecules to react with chemisorbed oxygen to produce two *OH per H 2 O with an activation energy 0.25 eV lower than that of the first hydrogen addition step. Finally, this work reveals the origin of the observed accelerated water formation reaction at the BL-silica/Ru(0001) interface in the low-temperature regime (T < 350 K) and points to a route to engineer chemical reaction pathways by leveraging subnanoscale confined spaces at metal–oxide interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamic structural and microstructural responses of a metal–organic framework type material to carbon dioxide under dual gas flow and supercritical conditions

The structural and microstructural responses of a model metal–organic framework material, Ni(3-methyl-4,4'-bipyridine)[Ni(CN) 4 ] (Ni-BpyMe or PICNIC-21), to CO 2 adsorption and desorption are reported for in situ small-angle X-ray scattering and X-ray diffraction measurements under different gas pressure conditions for two technologically important cases. These conditions are single or dual gas flow (CO 2 with N 2 , CH 4 or H 2 at sub-critical CO 2 partial pressures and ambient temperatures) and supercritical CO 2 (with static pressures and temperatures adjusted to explore the gas, liquid and supercritical fluid regimes on the CO 2 phase diagram). The experimental results are compared with density functional theory calculations that seek to predict where CO 2 and other gas molecules are accommodated within the sorbent structure as a function of gas pressure conditions, and hence the degree of swelling and contraction in the associated structure spacings and void spaces. Furthermore, these predictions illustrate the insights that can be gained concerning how such sorbents can be designed or modified to optimize the desired gas sorption properties relevant to enhanced gas recovery or to addressing carbon dioxide reduction through carbon mitigation, or even direct air capture of CO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗