Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Emphanitic anharmonicity in PbSe at high temperature and anomalous electronic properties in the Pb Q ( Q = S , Se , Te ) system

Here, the temperature dependence of the local structure of PbSe has been investigated using pair distribution function (PDF) analysis of x-ray and neutron powder diffraction data and density functional theory (DFT) calculations. Observation of non-Gaussian PDF peaks at high temperature indicates the presence of significant anharmonicity, which can be modeled as Pb off-centering along [100] directions that grows on warming similar to the behavior seen in PbTe and PbS and sometimes called emphanisis. Interestingly, the emphanitic response is smaller in PbSe than in both PbS and PbTe indicating a nonmonotonic response with chalcogen atomic number in the PbQ (Q = S, Se, Te) series. The DFT calculations indicate a correlation between band gap and the amplitude of [100] dipolar distortion, suggesting that emphanisis may be behind the anomalous composition and temperature dependencies of the band gaps in this series.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Influence of basicity on 1,3-butadiene formation from catalytic 2,3-butanediol dehydration over γ-alumina

The direct catalytic conversion of 2,3-butanediol (BDO) to 1,3-butadiene (BD) was studied over two commercial forms of alumina (denoted as F200 and SCFa) at temperatures between 240 C and 450 C. Even though these two catalysts are both high surface area forms of c-alumina, they gave remarkably different results, with SCFa giving higher BD selectivities at all experimental conditions. The difference is attributed to the higher surface area of F200, which means a greater number of acid sites that can convert BDO to methyl ethyl ketone (MEK). NH3 and CO2-TPD results supported this conclusion by showing that the two forms of alumina had different acid/base properties. Experimental results also showed that BD selectivity was improved by increasing temperature, increasing residence time and co-feeding water. The residence time study combined with density functional theory (DFT) calculations proved that 3-buten-2-ol (3B2OL) is an important intermediate in the conversion of BDO to BD. BD selectivity decreases over sodium modified alumina SCFa. It is hypothesized that on sodium-modified alumina, 3B2OL is dehydrogenated to form methyl vinyl ketone (MVK) as opposed to dehydration to BD. Basic sites catalyzed the retro-aldol condensation of MVK, which produces acetone and formaldehyde via cleavage of the C@C bond. This is in agreement with DFT calculations showing that the proposed pathway for acetone formation is more energetically favored on Na-modified c-Al2O3 (110) surface compared to the pristine (110) surface.

Zeng, Fan↗
DOE OSTI2019

Low-index surface energies, cleavage energies, and surface relaxations for crystalline NiAl from first-principles calculations

NiAl surfaces frequently serve as a platform for studying a broad range of physical and chemical phenomena including chemisorption, catalysis, oxidation, alloy growth, and surface nanostructure formation. Knowledge of precise values for low-index surface energies of NiAl, the most fundamental quantities characterizing surface thermodynamics, is often invaluable for understanding of these phenomena. In all previous analyses for NiAl(100) and NiAl(111), Ni- and Al-terminations are not distinguished, and half of the cleavage energy (or equivalently, the average of the surface energies of two differently terminated surfaces) is always identified as the “surface energy”. No values are available for individual surface energies of Ni- or Al-terminated NiAl(100) or NiAl(111) surfaces, whereas knowledge of only cleavage energy is often insufficient for analyzing surface-associated behavior. As such, in this work we perform extensive first-principles density-functional-theory (DFT) calculations for surface energies and cleavage energies of NiAl(110), NiAl(100) and NiAl(111) by considering the chemical-potential-based formulations to clarify the ambiguity in their surface energies and cleavage energies. We obtain a surface-energy phase diagram for these three low-index surfaces versus the relevant chemical potential, as well as the chemical-potential-dependent Wulff plots for NiAl crystal equilibrium shapes. We also provide the surface-relaxation information from our DFT calculations for comparison with previous experimental data.

42 ENGINEERING↗
DOE OSTI2025

Observation of multiple flat bands and Van Hove singularities in the distorted kagome metal NdTi 3 Bi 4

Kagome materials have attracted enormous research interest recently owing to their diverse topological phases and manifestation of electronic correlation. Here, in this study, we present the electronic structure of a distorted ferromagnetic kagome metal, NdTi 3 ⁢Bi 4 , exhibiting a transition temperature of 9 K. Our investigation employs a combination of angle-resolved photoemission spectroscopy (ARPES) measurements and density functional theory (DFT) calculations. We discover the presence of two “flat” bands which are found to originate from the kagome structure formed by Ti atoms with a major contribution from Ti 𝑑 𝑥⁢𝑦 and Ti 𝑑 𝑥 2 −𝑦 2 orbitals. We also observed multiple Van Hove singularities (VHSs) in its electronic structure, with one VHS lying near the Fermi level. The ARPES data reveal the existence of Dirac cone at the $\overline{𝐾}$ point, a finding which is corroborated by our DFT calculations. These findings present a detailed electronic structure capable of hosting correlation-driven phenomenon in this novel ferromagnetic kagome metal.

Fermi surface↗
DOE OSTI2015

Al{sub 15}Ge{sub 4}Ni{sub 3}: A new intergrowth structure with Cu{sub 3}Au- and CaF{sub 2}-type building blocks

The new ternary compound Al{sub 15}Ge{sub 4}Ni{sub 3} (τ{sub 2} in the system Al–Ge–Ni) was synthesized in single crystalline form by a special annealing procedure from samples located in the three phase fields [L+Al+τ{sub 2}] and [L+Ge+τ{sub 2}]. The crystal structure of Al{sub 15}Ge{sub 4}Ni{sub 3} was determined by single-crystal X-ray diffraction. The compound crystallizes in a new structure type in space group I4-bar3m, Pearson Symbol cI88, cubic lattice parameter a=11.405(1) Å. Phase diagram investigations indicate stoichiometric composition without considerable homogeneity range; τ{sub 2} melts peritectically at T=444 °C. The crystal structure of Al{sub 15}Ge{sub 4}Ni{sub 3} shows a unique combination of simple Cu{sub 3}Au- and CaF{sub 2}-type building blocks: a three dimensional network of CaF{sub 2}-type units, formed by Ni and Al atoms, is interspaced by clusters (Al{sub 6}Ge{sub 8}) resembling unit cells of the Cu{sub 3}Au-type. Both structural motifs are connected by Al–Ge bonds. The ground state energy of the compound was obtained by DFT calculations and the densities of states were analyzed in detail. In addition, electron density maps were calculated in four different sections through the unit cell using the full potential linearized augmented plane-wave (FLAPW) method. The bonding situation in Al{sub 15}Ge{sub 4}Ni{sub 3} was discussed combining results from electronic calculations with the analysis of the coordination of atoms. - Graphical abstract: The new compound Al{sub 15}Ge{sub 4}Ni{sub 3} shows a unique combination of simple Cu{sub 3}Au- and CaF{sub 2}-type building blocks. - Highlights: • The crystal structure of Al{sub 15}Ge{sub 4}Ni{sub 3} (space group I4-bar3m) was determined. • It shows a unique combination of CaF{sub 2}- and Cu{sub 3}Au-type building blocks. • Electronic (DFT) calculations were performed to gain insight to chemical bonding.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Relativistic Effects on a Metal–Metal Bond: Osmium Corrole Dimers

A series of metal–metal bonded osmium corrole dimers, {Os[T p XPC]} 2 , were synthesized in reasonably good yields (35–46%) via the interaction of the corresponding free-base meso -tris( p -X-phenyl)corroles (H 3 [T p XPC], X = CF 3 , H, CH 3 , and OCH 3 ), Os 3 (CO) 12 , and potassium carbonate in 1,2,4-trichlorobenzene under an inert atmosphere at 180 °C over several hours. The complexes are only the second class of Os corroles reported to date (the first being Os VI N corroles) and also the second class of metal–metal bonded metallocorrole dimers (the other being Ru corrole dimers). Comparison of the X-ray structures, redox potentials, and optical spectra of analogous Ru and Os corrole dimers, along with scalar-relativistic DFT calculations, has provided an experimentally calibrated account of relativistic effects in these complexes. Three of the Os corrole dimers (X = CF 3 , H, and OCH 3 ) were analyzed with single-crystal X-ray diffraction analysis, revealing inversion-related corrole rings with eclipsed Os–N bonds and Os–Os distances of ~2.24 Å that are ~0.06 Å longer than the Ru–Ru distances in the analogous Ru corrole dimers. Interestingly, a comparison of scalar-relativistic and nonrelativistic DFT calculations indicates that this difference in metal–metal bond distance does not, in fact, reflect a differential relativistic effect. Additionally, for a given corrole ligand, the Ru and Os corrole dimers exhibit nearly identical oxidation potentials but dramatically different reduction potentials, with the Os values ~0.5 V lower relative to Ru, suggesting that whereas oxidation occurs in a ligand-centered manner, reduction is substantially metal-centered, which indeed was confirmed by scalar-relativistic calculations. The calculations further indicate that approximately a third of the ~0.5 V difference in reduction potentials can be ascribed to relativity. The somewhat muted value of this relativistic effect appears to be related to the finding that reduction of an Os corrole dimer is not exclusively metal-based but that a significant amount of spin density is delocalized over to the corrole ligand; in contrast, reduction of an Ru corrole dimer occurs exclusively on the Ru–Ru linkage. For isoelectronic complexes, the Ru and Os corrole dimers exhibit substantially different UV–vis spectra. Furthermore, a key difference is a strong near-UV feature of the Os series, which in energy terms is blue-shifted by ~0.55 V relative to the analogous feature of the Ru series. TDDFT calculations suggest that this difference may be related to higher-energy Os(5d)-based LUMOs in the Os case relative to analogous MOs for Ru.

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

A phase width for CaGaSn. Crystal structure of mixed intermetallic Ca{sub 4}Ga{sub 4+x}Sn{sub 4−x} and SmGa{sub x}Sn{sub 3−x}, stability, geometry and electronic structure

X-ray single-crystal structure has been established for new compositions in intermetallic systems of tin and gallium. Crystals were successfully obtained in alloys prepared from elements. The structure of SmGaSn{sub 2} (cubic Pm3̄m, a=4.5778(8) Å, Z=1, R1=0.012) is described with atomic disorder at all Sn/Ga positions and the structure of Ca{sub 4}Ga{sub 4.9}Sn{sub 3.1} (hexagonal, P6{sub 3}/mmc, a=4.2233(9), c=17.601(7) Å, Z=1, R1=0.062) raises an interesting question about existence of a composition domain for CaGaSn. Finally, Ca{sub 4}Ga{sub 4.9}Sn{sub 3.1} should be considered as a particular composition of Ca{sub 4}Ga{sub 4+x}Sn{sub 4−x}, a compound assumed to exist in the range x ~ 0−1. Partial atomic ordering characterizes the Sn/Ga puckered layers of hexagons whose geometries are analyzed and discussed comparatively with analogous arrangements in AlB{sub 2} related hexagonal compounds. The study is supported by rigid band model and DFT calculations performed for different experimental and hypothetic arrangements. - Graphical abstract: A phase width for Ca{sub 4}Ga{sub 4+x}Sn{sub 4−x} belonging to the hexagonal YPtAs structure-type. - Highlights: • Single crystals of mixed tin gallium ternary intermetallics were obtained. • Partial ordering at metal sites and phase width are evidenced for Ca{sub 4}Ga{sub 4+x}Sn{sub 4−x}. • Layer deviation to flatness is studied comparatively with related structures. • Geometry and stability analyses based on DFT calculations are provided.

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

Fathoming the anisotropic magnetoelasticity and magnetocaloric effect in GdNi

Intermetallic GdNi adopts a CrB type of crystal structure (space group Cmcm), and it orders ferromagnetically via a second-order phase transition at 70 K, exhibiting unusually strong spontaneous striction along the three independent crystallographic axes in the ferromagnetically ordered state. Here we introduce a microscopic model to describe anisotropic changes of lattice parameters and elastic contribution to magnetocaloric effect of GdNi. In the model, results of density functional theory (DFT) calculations are used as inputs into a Hamiltonian that includes elastic energy of an anisotropic crystal lattice, exchange interactions, and Zeeman effect. The magnetic and elastic Hamiltonians are coupled through an anisotropic Bean-Rodbell model of magnetoelastic interactions. This coupling gives rise to anisotropic changes in the lattice parameters observed experimentally, and the model reveals good to reasonable agreements between the current theoretical results and earlier experimental data, thus validating the model within the limits of assumptions made. We also show that DFT calculations with 4f electrons of Gd treated as core electrons lead to a more adequate estimate of elastic constants of GdNi in comparison with the LDA + U method where 4f electrons are treated as valence electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Visualizing Redox Dynamics of a Single Ag/AgCl Heterogeneous Nanocatalyst at Atomic Resolution

Operando characterization of gas solid reactions at the atomic scale is of great importance for determining the mechanism of catalysis. This is especially true in the study of heterostructures because of structural correlation between the different parts. However, such experiments are challenging and have rarely been accomplished. In this work, atomic scale redox dynamics of Ag/AgCl heterostructures have been studied using in situ environmental transmission electron microscopy (ETEM) in combination with density function theory (DFT) calculations. The reduction of Ag/AgCl to Ag is likely a result of the formation of Cl vacancies while Ag+ ions accept electrons. The oxidation process of Ag/AgCl has been observed: rather than direct replacement of Cl by O, the Ag/AgCl nanocatalyst was first reduced to Ag, and then Ag was oxidized to different phases of silver oxide under different O-2 partial pressures. Ag2O formed at low O-2 partial pressure, whereas AgO formed at atmospheric pressure. By combining in situ ETEM observation and DFT calculations, this structural evolution is characterized in a distinct nanoscale environment.

Ag/AgCl heterogeneous nanocatalyst↗
DOE OSTI2020

Elucidation of Active Sites in Aldol Condensation of Acetone over Single-Facet Dominant Anatase TiO 2 (101) and (001) Catalysts

Aldol condensations of carbonyl compounds for C–C bond formation are a very important class of reactions in organic synthesis and upgrading of biomass-derived feedstocks. However, the atomic level understanding of reaction mechanisms and structure-activity correlation on widely used transition metal oxide catalysts are limited due to the high degree of struc-tural heterogeneity of catalysts such as commercial TiO 2 powders. Here, we provide a deep understanding of the reaction mechanisms, kinetics, and structure–function relationships for vapor phase acetone aldol condensation through the con-trolled synthesis of two catalysts with high surface areas and clean, dominant facets, coupled with detailed characterization and kinetic studies that are further assisted by density functional theory (DFT) calculations. Temperature-dependent diffuse reflectance infrared Fourier transform spectroscopy showed the existence of abundant acetone bonded to surface hydroxyl groups (acetone-O s H) and acetone bonded to Lewis acid sites (acetone-Ti 5c ) on surface of both {101} and {001} facet dom-inant TiO 2 . Intermolecular C–C coupling of enolate intermediate from acetone-Ti 5c and a vicinal acetone-O s H is a kinetically relevant step, which is consistent with kinetic and isotopic studies and DFT calculations. The {001} facet showed a lower apparent activation energy (or higher reactivity) than the {101} facet. This is likely caused by the weaker Lewis acid and Brønsted base strengths of the {001} facet which favors the reprotonation–desorption of the coupled intermediate, making the C–C coupling step more exothermic on the {001} facet and resulting in an earlier transition state with a lower activation barrier. It is also possible that the {001} facet has a smoother surface configuration and less steric hindrance during intermo-lecular C–C bond formation than the {101} facet.

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

CO 2 -assisted ethane aromatization over zinc and phosphorous modified ZSM-5 catalysts

The tandem reactions of CO 2 -assisted oxidative dehydrogenation and aromatization (CO 2 -ODA) of ethane were investigated over ZSM-5 modified with zinc (Zn) and phosphorus (P) by experimental studies and DFT calculations. The nature of active Zn species and their interaction with P were elucidated for understanding Zn/P-ZSM-5 catalyzed CO 2 -ODA reaction of ethane. Compared to the direct dehydrogenation and aromatization (DDA) reaction, the Zn/P-ZSM-5 catalyst possessed higher ethane conversion, enhanced aromatics yield, and better stability in the CO 2 -ODA reaction. Here, the presence of CO 2 reduced the carbon deposition, stabilized the Zn species, and thereby facilitated the generation of more liquid aromatics in CO 2 -assisted ethane aromatization. In situ characterization and DFT calculations provided insight into the mechanisms of CO 2 -ODA over the Zn/P-ZSM-5 catalyst.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spin-informed universal graph neural networks for simulating magnetic ordering

The screening and discovery of magnetic materials are hindered by the computational cost of first-principles density-functional theory (DFT) calculations required to find the ground state magnetic ordering. Although universal machine-learning interatomic potentials (uMLIPs), also known as atomistic foundation models, offer high-fidelity models of many atomistic systems with significant speedup, they currently lack the inputs required for predicting magnetic ordering. In this work, we present a data-efficient, spin-informed graph neural network framework that incorporates spin degrees of freedom as inputs and preserves physical symmetries, extending the functionality of uMLIPs to simulate magnetic orderings. This framework speeds up DFT calculations through better initial guesses for magnetic moments, determines the ground-state ordering of bulk materials and even generalizes to magnetic ordering in surfaces. Furthermore, we implement a closed-loop anomaly detection approach that effectively addresses the classic "chicken-and-egg" problem of creating a high-quality dataset while developing a uMLIP, unearthing anomalies in large benchmark datasets and boosting model accuracy.

Xu, Wenbin↗
DOE OSTI2018

Understanding the Intrinsic Surface Reactivity of Single-Layer and Multilayer PdO(101) on Pd(100)

In this work, we investigated the intrinsic reactivity of CO on single and multilayer PdO(101) grown on Pd(100) using temperature programmed reaction spectroscopy (TPRS) and reflection absorption infrared spectroscopy (RAIRS) experiments as well as density functional theory (DFT) calculations. We find that CO binds more strongly on multilayer than single-layer PdO(101) (~119 vs. 43 kJ/mol), and that CO oxidizes negligibly on single-layer PdO(101) whereas nearly 90% of a saturated layer of CO oxidizes on multilayer PdO(101) during TPRS experiments. RAIRS further shows that CO molecules adsorb on both bridge and atop-Pd cus sites (coordinatively-unsaturated Pd sites) of single-layer PdO(101)/Pd(100), while CO binds exclusively on atop-Pd cus sites of multilayer PdO(101). The DFT calculations reproduce the much stronger binding of CO on multilayer PdO(101) as well as the observed binding site preferences, and reveal that the stronger binding is entirely responsible for the higher CO oxidation activity of multilayer PdO(101)/Pd(100). We show that the O-atom below the Pd cus site, present only on multi-layer PdO(101), modifies the electronic states of the Pd cus atom in a way that enhances the CO-Pd cus bonding. Lastly, we show that a precursor-mediated kinetic model, with energetics determined from the present study, predicts that the intrinsic CO oxidation rates achieved on both single and multilayer PdO(101)/Pd(100) can be expected to exceed the gaseous CO diffusion rate to the surface during steady-state CO oxidation at elevated pressures, even though the intrinsic reaction rates are 4-5 orders of magnitude lower on single layer than on multilayer PdO(101)/Pd(100).

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

Surface chemistry of aromatic reactants on Pt- and Mo-modified Pt catalysts

Supported catalysts containing an oxophilic metal such as Mo and a noble metal such as Pt have shown promising activity and selectivity for deoxygenation of biomass-derived compounds. Here, we report that PtMo catalysts also promote hydrogenolysis of the model compound benzyl alcohol, while decarbonylation is most prevalent over unmodified Pt. A combination of single crystal surface science studies, density functional theory (DFT) calculations, and vapor phase upgrading experiments using supported catalysts was carried out to better understand the mechanism by which Mo promotes deoxygenation. Molybdenum was deposited in submonolayer quantities on a Pt(111) surface and reduced at high temperature. Temperature-programmed desorption (TPD) experiments using benzyl alcohol as a reactant showed greatly enhanced yields of the deoxygenation product toluene at moderate Mo coverages. To understand how the interaction of the aromatic group with the surface influenced this reactivity, we investigated the adsorption of toluene as a probe molecule. We found that the addition of Mo to Pt(111) resulted in a significant decrease in toluene decomposition. DFT calculations indicated that this decrease was consistent with decreased aromatic adsorption strengths that accompany incorporation of Mo into the Pt subsurface. The weaker aromatic-surface interaction on Pt/Mo surfaces led to a tilted adsorption geometry for benzyl alcohol, which presumably promotes hydrogenolysis to produce toluene instead of decarbonylation to produce benzene and CO. Alumina-supported Pt and PtMo catalysts were also tested for benzyl alcohol deoxygenation. PtMo catalysts had a higher rate of toluene production and lower rates of benzene and benzaldehyde production. Additionally, when benzaldehyde was used as the reactant to measure decarbonylation activity the mass-normalized rate of benzene production was 2.5 times higher on Pt than PtMo. Altogether, the results of TPD, DFT, and supported catalyst experiments suggest that subsurface Mo sites weaken the binding of aromatic rings on PtMo surfaces; the weakened aromatic-surface interaction is correlated with an improvement in selectivity to C-O bond scission.

09 BIOMASS FUELS↗
DOE OSTI2016

Highly Preorganized Ligand 1,10-Phenanthroline-2,9-dicarboxylic Acid for the Selective Recovery of Uranium from Seawater in the Presence of Competing Vanadium Species

Studies of the complexation of new promising ligands with uranyl (UO 2 2+ ) and other seawater cations can aid the development of more efficient, selective, and robust sorbents for the recovery of uranium from seawater. Here, we propose that the ligand design principles based on structural preorganization can be successfully applied to obtain a dramatic enhancement in UO 2 2+ ion binding affinity and selectivity. This concept is exemplified through the investigation of the com-plexes of UO 2 2+ , VO 2+ , and VO 2+ with the highly preorganized ligand PDA (1,10-phenanthroline-2,9-dicarboxylic acid) using a combination of fluores-cence and absorbance techniques, along with den-sity functional theory (DFT) calculations. Moreover, the measured stability constant value, log K1, of 16.5 for the UO 2 2+ /PDA complex is very high compared to uranyl complexes with other dicarboxylic ligands. Moreover, PDA exhibits strong selectivity for uranyl over vanadium ions, since the determined sta-bility constant values of the PDA complexes of the vanadium ions are quite low (V(IV) log K1 = 7.4, V(V) = 7.3). Finally, the structures of the corresponding UO 2 2+ , VO 2+ , and VO 2+ complexes with PDA were identified by systematic DFT calculations, and helped to interpret the stronger binding affinity for uranium over the vanadium ions. Due to its high chemical stability, selectivity, and structural preor-ganization for UO 2 2+ complexation, PDA is a very promising candidate that can be potentially used in the development of novel adsorbent materials for the selective extraction of uranium from sea-water.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Chemical Pressure-Driven Enhancement of the Hydrogen Evolving Activity of Ni 2 P from Nonmetal Surface Doping Interpreted via Machine Learning

The activity of Ni 2 P catalysts for the hydrogen evolution reaction (HER) is currently limited by strong H adsorption at the Ni 3 -hollow site. We investigate the effect of surface nonmetal doping on the HER activity of the Ni 3 P 2 termination of Ni 2 P(0001), which is stable at modest electrochemical conditions. Using density functional theory (DFT) calculations, we find that both 2p nonmetals and heavier chalcogens provide nearly thermoneutral H adsorption at moderate surface doping concentrations. We also find, however, that only chalcogen substitution for surface P is exergonic. For intermediate surface concentrations of S, the free energy of H adsorption at the Ni 3 -hollow site is -0.11 eV, which is significantly more thermoneutral than the undoped surface (-0.45 eV). We use the regularized random forest machine learning algorithm to discover the relative importance of structure and charge descriptors, extracted from the DFT calculations, in determining the HER activity of Ni 2 P(0001) under different doping concentrations. Here, we discover that the Ni–Ni bond length is the most important descriptor of HER activity, which suggests that the nonmetal dopants induce a chemical pressure-like effect on the Ni 3 -hollow site, changing its reactivity through compression and expansion.

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

Molecular Level Understanding of Polyethylene Terephthalate (PET) Depolymerization in Base/Alcohol Hybrid Systems

Polyethylene terephthalate (PET) depolymerization in base/alcohol hybrid systems represents a promising low-energy approach for chemically recycling PET waste into valuable monomers. This study investigates the mechanistic pathways of PET depolymerization in NaOH/alcohol solutions, emphasizing the competing roles of hydroxide and alkoxide species. Utilizing a combination of experimental techniques, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations, we explore how factors such as base concentration, alcohol chain length, and pKa values of alcohols influence PET depolymerization efficiency and pathways. Our findings indicate that alkoxide ions (RO⁻) exhibit notably higher reactivity than hydroxide ions (HO⁻), favoring an alcoholysis pathway in the base/alcohol hybrid system. Experimental results across a series of C1 to C5 alcohols show that longer-chain alcohols, particularly 1-butanol, achieve higher PET conversion, although this does not align solely with simple nucleophilicity trends of alkoxides. While DFT calculations reveal comparable activation energies for various alkoxides in PET depolymerization, MD simulations underscore the significant role of alcohol chain length, with longer-chain alcohols forming more stable or frequent interactions with PET. Additionally, the alkoxide concentration, influenced by the alcohol’s pK a , directly impacts PET conversion. These suggest that PET depolymerization is governed by a balance between alkoxide concentration and alkoxide-PET interactions, rather than activation energies or nucleophilicity alone. From a practical perspective, incorporating long-chain alcohols as cosolvents may enhance process efficiency but increases raw material costs by approximately 30%. However, long-chain alcohols present a safer and more sustainable alternative to hazardous cosolvents such as dichloromethane. This work offers a molecular-level understanding of PET depolymerization in base/alcohol systems and provides insights into optimizing these systems for more efficient and sustainable PET recycling processes.

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

Hydrogen desorption from the surface and subsurface of cobalt

The influence of coverage on the diffusion of hydrogen into the subsurface of cobalt was studied using density functional theory (DFT) and temperature programmed desorption (TPD). DFT calculations show that as the hydrogen coverage on Co(0001) increases, the barrier for hydrogen diffusion into the bulk decreases by 20%. Additionally, subsurface hydrogen on a hydrogen covered surface was found to be more stable when compared to a clean cobalt surface. To test these theoretical findings experimentally, excited hydrogen was used in an ultra-high vacuum environment to access higher hydrogen coverages. Here, our TPD studies showed that at high hydrogen coverages, a sharp low temperature feature appeared, indicating the stabilization of subsurface hydrogen. Further DFT calculations indicate that this sharp low temperature feature results from associative hydrogen desorption from a hydrogen saturated surface with a population of subsurface hydrogen. Microkinetic modelling was used to model the TPD spectra for hydrogen desporption from cobalt with and without subsurface hydrogen, showing reasonable agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗