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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Pb:Sn Ratio-Driven Polytypism and Band Modulation in Photoresponsive Hexagonal Perovskitoids

Three-dimensional (3D) perovskites of the formula AMX 3 are known for their excellent optoelectronic properties, but their design is limited by the narrow range of A-site cations that can template the 3D corner-sharing structure, and many of the viable options have already been explored. These materials also face structural instability under environmental conditions. In contrast, 3D hexagonal perovskitoids, with the same chemical formula, may offer enhanced stability and richer structural diversity through a range of corner- and face-sharing octahedral configuration options, providing greater opportunities for structural design; however, challenges such as synthesis complexity and wide band gaps have hindered optoelectronic performance achieved to date. Herein, we synthesized a structural homologous series of mixed-metal hexagonal perovskitoids with the formula APb 1-x Sn x I 3 (x = 0, 0.25, 0.50, 0.75, 1; A = ethylammonium, guanidinium) and identified three polytypes (9R, 12R, 6H) using single-crystal X-ray diffraction (SCXRD). These structures exhibited increasing corner-sharing connectivity with higher Sn content, revealing a previously unobserved relationship between metal composition and structural evolution in perovskitoid materials. The incorporation of Sn reduced the band gap (tunable from 2.51 eV to 1.87 eV) and drove structural transformations, a trend seen also in density functional theory (DFT) calculations, which suggest a thermodynamic preference for Pb at face-sharing sites and Sn at corner-sharing sites. DFT band structure calculations and optical spectroscopy also reveal an anomalous behavior in these materials, which we term polytypic band modulation. This phenomenon combines conventional band bowing with the structural transformations that occur as Sn content increases, as observed in the band gap and in photoluminescence spectra. Photodiodes fabricated from thin films of these materials exhibited stable and pronounced photoresponses across various light intensities over time. Furthermore, the combination of templating 3D perovskitoids with multiple cations and alloying Pb and Sn, suggests a vastly underexplored phase space that offers new parameters to tune perovskitoids.

Cations↗
DOE OSTI2023

Dominant Role of Hole Transport Pathway in Achieving Record High Photoconductivity in Two‐Dimensional Metal–Organic Frameworks

Abstract Metal–organic frameworks (MOFs) with mobile charges have attracted significant attention due to their potential applications in photoelectric devices, chemical resistance sensors, and catalysis. However, fundamental understanding of the charge transport pathway within the framework and the key properties that determine the performance of conductive MOFs in photoelectric devices remain underexplored. Herein, we report the mechanisms of photoinduced charge transport and electron dynamics in the conductive 2D M−HHTP (M=Cu, Zn or Cu/Zn mixed; HHTP=2,3,6,7,10,11‐hexahydroxytriphenylene) MOFs and their correlation with photoconductivity using the combination of time‐resolved terahertz spectroscopy, optical transient absorption spectroscopy, X‐ray transient absorption spectroscopy, and density functional theory (DFT) calculations. We identify the through‐space hole transport mechanism through the interlayer sheet π–π interaction, where photoinduced hole state resides in HHTP ligand and electronic state is localized at the metal center. Moreover, the photoconductivity of the Cu−HHTP MOF is found to be 65.5 S m −1 , which represents the record high photoconductivity for porous MOF materials based on catecholate ligands.

14 SOLAR ENERGY↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2020

Ternary PtIrNi Catalysts for Efficient Electrochemical Ammonia Oxidation

Due to high energy density, carbon-free feature, and easiness to be liquefied, ammonia (NH 3 ) has proved to be an effective alternative to hydrogen in low temperature fuel cells via its direct the ammonia oxidation reaction (AOR) for electricity generation. However, the kinetically sluggish AOR has prohibitively hindered the attractive direct ammonia fuel cell (DAFC) applications. Here we report an efficient AOR catalyst, in which ternary PtIrNi alloy nanoparticles well dispersed on a binary composite support consisting of porous silicon dioxide (SiO 2 ) and carboxyl-functionalized carbon nanotube (PtIrNi/SiO 2 -CNT-COOH) through a sonochemical-assisted synthesis strategy. The PtIrNi alloy nanoparticles, with the aid of abundant OHad provided by porous SiO 2 , and the improved electrical conductivity by CNTs, exhibit remarkable catalytic activity for the AOR in alkaline media. It is evidenced by a lower onset potential (~0.40 V vs. RHE) at room temperature, than that of commercial PtIr/C (ca. 0.43 V vs. RHE). Increasing NH 3 concentrations and operation temperatures significantly improve AOR performance of this catalyst. Specifically, AOR activity of the optimal PtIrNi nanoparticle catalyst can be significantly enhanced by elevating the temperature to 80ºC, with a much lower onset potential (~0.32 V vs. RHE), indicating that DAFC can be operated at higher temperature for increased performance. Constant-potential density functional theory (DFT) calculations showed that the Pt-Ir ensembles on {100}-terminated surfaces serve as the active site. Importantly, the introduction of Ni raises the center energy of the density of states projected onto the group d-orbitals of surface sites and thus lowers the theoretical onset potential for *NH 2 dehydrogenation to *NH when compared to Pt and Pt 3 Ir alloy.

25 ENERGY STORAGE↗
DOE OSTI2023

Au/Pt Bimetallic Nanowires with Stepped Pt Sites for Enhanced C–C Cleavage in C2+ Alcohol Electro-oxidation Reactions

Efficient C–C bond cleavage and oxidation of alcohols to CO 2 is the key to developing highly efficient alcohol fuel cells for renewable energy applications. In this work, we report the synthesis of core/shell Au/Pt nanowires (NWs) with stepped Pt clusters deposited along the ultrathin (2.3 nm) stepped Au NWs as an active catalyst to effectively oxidize alcohols to CO 2 . The catalytic oxidation reaction is dependent on the Au/Pt ratios, and the Au 1.0 /Pt 0.2 NWs have the largest percentage (~75%) of stepped Au/Pt sites and show the highest activity for ethanol electro-oxidation, reaching an unprecedented 196.9 A/mg Pt (32.5 A/mg Pt+Au ). This NW catalyst is also active in catalyzing the oxidation of other primary alcohols, such as methanol, n-propanol, and ethylene glycol. In situ X-ray absorption spectroscopy and infrared spectroscopy are used to characterize the catalyst structure and to identify key reaction intermediates, providing concrete evidence that the synergy between the low-coordinated Pt sites and the stepped Au NWs is essential to catalyze the alcohol oxidation reaction, which is further supported by DFT calculations that the C–C bond cleavage is indeed enhanced on the undercoordinated Pt–Au surface. Here, our study provides important evidence that a core/shell structure with stepped core/shell sites is essential to enhance electrochemical oxidation of alcohols and will also be central to understanding electro-oxidation reactions and to the future development of highly efficient direct alcohol fuel cells for renewable energy applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Reversible dual anionic-redox chemistry in NaCrSSe with fast charging capability

Utilizing the anionic redox reaction opens new approaches for the development of new cathode materials with extra capacities. Although, it suffers from several obstacles such as voltage hysteresis and sluggish kinetics. In this paper, a new layered chalcogenide-based on dual anionic-redox reaction is reported. The newly designed layered NaCrSSe exhibits the capacity of almost all Na + intercalation/deintercalation (137 mAh/g -1 at 50 mA/g), and a unique charge/discharge feature with a small polarization of 0.15 V and high energy efficiency of ~92% in initial cycles. Furthermore, a superior high-rate charge capacity of 115.5mAh/g -1 (83.7% retention) was achieved at 27.8 C (4000 mA/g -1 ). Systematic characterization studies on structure evolution and DFT calculation show the charge compensation of S and Se anions during cycling. Overall, these results will enrich the anion redox chemistry and provide valuable information for developing new anion-redox based cathode materials with high capacity and fast kinetics.

25 ENERGY STORAGE↗
DOE OSTI2023

New Platinum Complexes from Salen- and Hydroxy-Substituted Salpn-Naphthalene Ligands with CO 2 Reduction Activity

The electrocatalytic reduction of carbon dioxide (CO 2 ) into added-value products is a promising alternative to completing the cycle of atmospheric CO 2 . We report two new platinum complexes—a salen-like naphthalene (PtL1) and a hydroxy-substituted salpn naphthalene (PtL2)—that are capable of activating CO 2 to produce carbon monoxide (CO). The predominant keto tautomer of the non-innocent ligands was determined using DFT calculations and UV-Vis spectroscopy. The PtL2 complex has a CO Faradaic efficiency >40% in the presence of water as a sacrificial proton source at -2.5 V vs. Fc/Fc + . The addition of the hydroxy group in combination with water as a proton source decreased the reduction potential and increased the CO formation tenfold when compared to PtL1.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Excited-State Charge Transfer in Push–Pull Platinum(II) π-Extended Porphyrins Fused with Pentacenequinone

Platinum(II) π-extended porphyrins fused with pentacenequinone and dihydropentacene have been successfully synthesized. These porphyrins were investigated using various techniques including absorption, steady-state, and time-resolved phosphorescence spectroscopy and differential pulse voltammetry. UV–vis absorption spectra of pentacenequinone-fused porphyrins (SW-Pt1 and SW-Pt2) showed unusually broad and nontypical absorption patterns. Phosphorescence spectra of SW-Pt1, SW-Pt2, and SW-Pt3 displayed similar emissions in the 704–706 nm region indicating electronic transitions of similar origin; however, the triplet lifetimes were found to be quenched in the case of both SW-Pt1 and SW-Pt2, suggesting the occurrence of excited-state events. Facile reductions were obtained for both the pentacene-quinone-fused monomer (SW-Pt2) and dimer (SW-Pt1) and were identified to be located at the pentacenequinone components. Here, the observed orbital segregations for SW-Pt2 and SW-Pt1 from DFT calculations supported the possibility of charge transfer in these push–pull systems. Interestingly, the established energy level diagram revealed that the charge transfer from the triplet excited Pt porphyrin is thermodynamically an uphill process. Systematic studies involving both femtosecond and nanosecond transient absorption techniques revealed that the singlet excited Pt porphyrins undergo an intermediate charge transfer state prior to populating the triplet state, providing a plausible explanation for phosphorescence quenching. The lifetime of the intermediate charge transfer states was found to be 25.9 and 5.68 ps, respectively, for SW-Pt1 and SW-Pt2.

Charge transfer↗
DOE OSTI2023

New Zintl Phase Yb 10 MgSb 9 with High Thermoelectric Performance

Abstract Yb 10 MgSb 9 is a new Zintl compound (with a composition closer to Yb 10.5 MgSb 9 ) and a promising thermoelectric material first reported in this work. Undoped Yb 10 MgSb 9 has an ultralow thermal conductivity due to crystallographic complexity and exhibits a relatively high peak p‐type Seebeck coefficient and high electrical resistivity. This is consistent with Zintl counting and density functional theory (DFT) calculations that the composition Yb 10.5 MgSb 9 should be a semiconductor. Na is found experimentally to be an effective p‐type dopant potentially due to the replacement of Na + for Yb 2+ , allowing for a significant decrease in electrical resistivity. With doping, a dramatic improvement of electrical conductivity is observed and the glass‐like thermal conductivity remains low, allowing for a significant enhancement of the thermoelectric figure of merit, zT . Doping increases the zT from 0.23 in undoped Yb 10 MgSb 9 to 1.06 in 7 at% Na‐doped Yb 10 MgSb 9 at 873K. This high thermoelectric performance found through Na‐doping places this material amongst the leading p‐type Zintl thermoelectrics, making it a promising candidate for future studies and high‐temperature thermoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A low-cost sulfate-based all iron redox flow battery

Redox flow batteries (RFBs) are promising choices for stationary electric energy storage. Nevertheless, commercialization is impeded by high-cost electrolyte and membrane materials. Here, we report a low-cost all-iron RFB that features inexpensive FeSO 4 electrolytes, microporous membrane along with a glass fiber separator. The addition of 0.1 м 1-ethyl-3-methylimidazolium chloride (EMIC) overcomes the low solubility of FeSO 4 in water which is raised to 2.2 м. DFT calculations demonstrate that EMI + can strengthen the interaction between sulfate anions and water molecules. This electrolyte composition also allows both anode and cathode reactions to operate without actively maintaining a pH gradient between them, thus eliminating the need for expensive ion exchange membranes. The all-iron RFB demonstrates stable operation at a current density of 20 mA cm –2 for more than 800 cycles via a simple, periodic regeneration process. Furthermore, the system cost of FeSO 4 /EMIC RFBs is projected to be $ 50 per kWh due to its low-cost active materials and the inexpensive microporous membrane. Finally, this low-cost, high-concentration all-iron RFB is a promising stationary energy–storage system for storing renewable energy.

25 ENERGY STORAGE↗
DOE OSTI2021

Crystal growth, density functional theory, and scintillation properties of Tl 3 LnCl 6 :Ce 3+ and TlLn 2 Cl 7 :Ce 3+ (Ln = Y, Gd)

In this paper we report on the crystal growth, density functional theory (DFT) calculations and scintillation properties of Tl 3 LnCl 6 :Ce and TlLn 2 Cl 7 :Ce(Ln = Y, Gd). Crystals were grown by the Vertical Bridgman technique up to 16 mm in diameter and 25 mm long. Crystals of Tl 3 YCl 6 :Ce and Tl 3 GdCl 6 :Ce belong to the family of the Nesohalides which have the monoclinic crystal structure. Crystals of TlY 2 Cl 7 :Ce and TlGd 2 Cl 7 :Ce belong to the family of the Phyllohalides which have either the monoclinic or the orthorhombic crystal structure. We report the light yields of these Nesohalides and Phyllohalides are typically on the order of 35,000–40,000 ph/MeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Nuclear–electronic orbital approach to quantization of protons in periodic electronic structure calculations

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

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

S K-edge XAS of Cu Ⅱ , Cu Ⅰ , and Zn Ⅱ oxidized Dithiolene complexes: Covalent contributions to structure and the Jahn-Teller effect

We report reduced dithiolene ligands are bound to high valent Mo centers in the active site of the oxotransferase family of enzymes. Related model complexes have been studied with great insight by Prof. Holm and his colleagues. This study focuses on the other limit of dithiolene chemistry: an investigation of the 2-electron oxidized dithiolene bound to low-valent late transition metal (TM) ions (Zn II , Cu I , and Cu II ). The bonding descriptions of the oxidized dithiolene [N,N-dimethyl piperazine 2,3-dithione (Me 2 Dt 0 )] complexes are probed using S K-edge X-ray absorption spectroscopy (XAS) and the results are correlated to density functional theory (DFT) calculations. These experimentally supported calculations are then extended to explain the different geometric structures of the three complexes. The Zn II (Me 2 Dt 0 ) 2 complex has only ligand-ligand repulsion so it is stabilized at the D 2d symmetry limit. The Cu I (Me 2 Dt 0 ) 2 complex has additional weak backbonding thus distorts somewhat from D 2d toward D 2h symmetry. The Cu II (Me 2 Dt 0 ) 2 complex has a strong σ donor bond that leads to both a large Jahn-Teller stabilization to D 2h and an additional covalent contribution to the geometry. The combined strong stabilization results in the square planar, D 2h structure. This study quantifies the competition between the ligand-ligand repulsion and the change in electronic structures in determining the final geometric structures of the oxidized dithiolene complexes, and provides quantitative insights into the Jahn-Teller stabilization energy and its origin.

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

Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)

Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.

36 MATERIALS SCIENCE↗
DOE OSTI2023

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

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

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

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2022

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2024

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

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

Findley, John↗