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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Scaling of Transition State Vibrational Frequencies and Application of d -Band Theory to the Brønsted–Evans–Polanyi Relationship on Surfaces

Semiempirical energy relations provide a means of estimating thermodynamic properties. Specifically, linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) relationships correlate adsorption energies between adsorbates across surfaces and reaction energies with activation barriers, respectively. Although vibrational scaling relations (VSRs) exist between adsorbates at identical sites, scaling between vibrational frequencies of adsorbed local minima and transition states is lacking. Here, we present density functional theory calculations for AH X (A = C, N, O) diffusions on transition metal surfaces and reveal linear scaling between frequencies of local minima and the transition state between those minima. Using d-band theory and linear muffin tin orbital theory (LMTO), we derive the slopes of these transition state vibrational scaling relations (TSVSRs) and, in so doing, provide a rigorous theory extending the original BEP relations developed for solution chemistry to surface chemistry. Furthermore, with a single reference DFT calculation, we predict the slopes and quantify uncertainty in the predictions.

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

An automatically curated first-principles database of ferroelectrics

Ferroelectric materials have technological applications in information storage and electronic devices. The ferroelectric polar phase can be controlled with external fields, chemical substitution and size-effects in bulk and ultrathin film form, providing a platform for future technologies and for exploratory research. In this work, we integrate spin-polarized density functional theory (DFT) calculations, crystal structure databases, symmetry tools, workflow software, and a custom analysis toolkit to build a library of known, previously-proposed, and newly-proposed ferroelectric materials. With our automated workflow, we screen over 67,000 candidate materials from the Materials Project database to generate a dataset of 255 ferroelectric candidates, and propose 126 new ferroelectric materials. We benchmark our results against experimental data and previous first-principles results. The data provided includes atomic structures, output files, and DFT values of band gaps, energies, and the spontaneous polarization for each ferroelectric candidate. Finally, we contribute our workflow and analysis code to the open-source python packages atomate and pymatgen so others can conduct analogous symmetry driven searches for ferroelectrics and related phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Correlating binding energies of adsorbed CO and H on model surfaces with CO/H 2 selectivity from co-electrolysis of CO 2 and H 2 O over copper–palladium bimetallic catalysts

Binding energies of adsorbed CO and H are key descriptors governing the activity and selectivity of the co-electrolysis of CO 2 and H 2 O to produce syngas with desired CO/H 2 ratios. Palladium hydride (PdH), which forms in situ at negative overpotentials, has been identified as the active Pd phase for CO 2 reduction to syngas. Herein, binding energies of CO and H are determined using temperature programmed desorption (TPD) of CO and H 2 from Pd(111), PdH/Pd(111), and Cu/PdH/Pd(111) under ultra-high vacuum (UHV) conditions. TPD results reveal that desorption of H 2 from subsurface PdH occurs at 460 K, while desorption from surface PdH is more facile at 320 K. CO desorption temperatures shift 20 K lower on PdH/Pd(111) compared to on Pd(111). The presence of 0.7 ML Cu further increases the desorption temperature of H 2 by 30 K while simultaneously reducing CO desorption temperatures by 70 K. Density functional theory (DFT) calculations show that CO adsorption onto Pd sites is hindered on the 0.7 ML Cu/PdH/Pd(111) surface while the kinetic barrier for H 2 desorption is increased. The trends in the binding energies of CO and H on model surfaces are consistent with electrochemical measurements of CuPd powder catalysts in a membrane electrode assembly (MEA), where H 2 evolution is reduced while CO production is enhanced compared to unmodified Pd catalysts. Overall, the results from model surface studies (TPD and DFT) provide a prediction and explanation for the activity and CO/H 2 ratios observed in electrochemical experiments. This study also demonstrates that CuPd is a promising catalyst with reduced Pd-loading to produce CO-rich syngas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Low-Temperature Activation and Coupling of Methane on MgO Nanostructures Embedded in Cu 2 O/Cu(111)

Here, the efficient conversion of methane into valuable hydrocarbons such as ethane and ethylene at relatively low temperatures without deactivation issues is crucial for advancing sustainable energy solutions. Herein, AP-XPS and STM studies show that MgO nanostructures (0.2-0.5 nm wide, 0.4-0.6 Å high) embedded in a Cu 2 O/Cu(111) substrate activate methane at room temperature, mainly dissociating it into CH x (x = 2 or 3) and H adatoms, with minimal conversion to C adatoms. These MgO nanostructures in contact with Cu 2 O/Cu(111) exhibit unique reactivity, enabling C-C coupling into ethane and ethylene at 500 K, a significantly lower temperature than that required for bulk MgO catalysts (>700 K), with negligible carbon deposition and no deactivation. DFT calculations corroborate these experimental findings. The CH 4,gas → *CH 3 +*H reaction is a downhill process on MgO/Cu 2 O/Cu(111) surfaces. The activation of methane is facilitated by an electron transfer from copper to MgO and the existence of Mg and O atoms with a low coordination number in the oxide nanostructures. The formation of O-CH 3 and O-H bonds overcomes the energy necessary for the cleavage of a C-H bond in methane. DFT studies reveal that smaller Mg 2 O 2 model clusters provide stronger binding and lower activation barriers for C-H dissociation in CH 4 , while larger Mg 3 O 3 clusters promote C-C coupling due to weaker *CH 3 binding. All these results emphasize the importance of size when optimizing the catalytic performance of MgO nanostructures in the selective conversion of methane.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Final Technical Report for "Cloud-based Low-Scaling Quantum Chemistry Simulations for Materials"

The objective of phase I was to develop the basic infrastructure for performing efficient hybrid DFT calculations on solid materials with a turnaround time amenable to the use in large-scale data-driven approaches. To fulfill the goal we have developed a novel algorithm that significantly reduces the cost of exchange matrix formation. The algorithm does so by making use of three ingredients (a) interpolative decomposition of the electron integrals (b) robust pseudospectral method and (c) occ-RI approach. Our published work demonstrates that the algorithm is orders of magnitude faster than any other hybrid-DFT method and is on the order of only 3-4 times slower than pure DFT. All steps of the algorithm developed during phase I can be accelerated to the point that, for large enough systems, the diagonalization of the Fock matrix will become the most expensive step. To treat such large systems we have created an interface to the ASCR-funded PEXSI library. The PEXSI method is used to compute the density matrix directly from the Fock matrix in a manner that preserves the sparsity of the local representation.

Shiozaki, Toru↗
DOE OSTI2025

Grain boundary zirconia-modified garnet solid-state electrolyte

Here, we report a method for promoting electrochemical stability in garnet Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 solid-state electrolyte based on a composite two-phase oxide–oxide microstructure. Grain boundary precipitation of the controlled distribution of amorphous zirconium oxide microparticles is achieved through the addition of reactive tantalum carbide. During ambient-atmosphere sintering, the carbide decomposes through an in situ reaction, the ‘extra’ Ta substituting for Zr within the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 lattice. Density functional theory (DFT) calculations identify a thermodynamically favourable reaction path and show how substituting Ta 5+ at Zr 4+ sites affects the crystal structure as well as bulk ionic and electronic conductivities. Quantitative stereology highlights that zirconia also acts as a sintering aid, reducing compact porosity. Cryogenic focused-ion-beam scanning electron microscopy and fractography analysis of cycled solid-state electrolytes illustrates that near-universally observed intergranular Li-metal dendrite propagation is suppressed by the two-phase microstructure, favouring transgranular dendrites instead. Importantly, DFT demonstrates that compared with the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 surface, the zirconium oxide surface per se is less electronically conductive and does not trap excess electrons to reduce Li ions. This is a key reason for the substantial improvement in the electrochemical properties over the single-phase baseline.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

C–H Bond Activation via U(II) in the Reduction of Heteroleptic Bis(trimethylsilyl)amide U(III) Complexes

Reduction of (C 5 Me 5 ) 2 U III (NR 2 ) and (C 5 Me 5 )U III (NR 2 ) 2 (R = SiMe 3 ) with potassium graphite in the presence of 2.2.2-cryptand (crypt) generates dark solutions that have UV-visible spectra consistent with time-dependent density functional theory (TDDFT) calculations on the U(II) products, [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - . However, the solutions quickly change color and form the U(III) C–H bond activation products [K(crypt)][(C 5 Me 5 ) 2 U III (CH 2 SiMe 2 NSiMe 3 –κC,κN)], 1, and [K(crypt)][(C 5 Me 5 )U III (NR 2 )(CH 2 SiMe 2 NSiMe 3 –κC,κN)], 2, that were identified by X-ray crystallography. DFT calculations on the putative [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - complexes revealed 5f 3 6d 1 ground state electron configurations as previously found in isolable [(C 5 H 4 SiMe 3 ) 3 U II ] - , which indicated that these low symmetry heteroleptic complexes are reasonable precursors for new U(II) complexes.

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

Magnetization reversal driven by electron localization-delocalization crossover in the inverse spinel Co 2 VO 4

Neutron diffraction, magnetization, and muon spin relaxation measurements, supplemented by density functional theory (DFT) calculations are employed to unravel temperature-driven magnetization reversal in inverse spinel Co 2 VO 4 . All measurements show a second-order magnetic phase transition at T C = 168 K to a collinear ferrimagnetic phase. Neutron diffraction measurements reveal two antiparallel ferromagnetic (FM) sublattices, belonging to magnetic ions on two distinct crystal lattice sites, where the relative balance between the two sublattices determine the net FM moment in the unit cell. As the evolution of the ordered moment with temperature differs between the two sublattices, the net magnetic moment reaches a maximum at T NC = 138 K and reverses its sign at T MR = 65 K. The DFT results suggest that the underlying microscopic mechanism for the reversal is a delocalization of the unfilled 3d -shell electrons on one sublattice just below T C , followed by a gradual localization as the temperature is lowered. Furthermore, this delocalized-localized crossover is supported by muon spectroscopy results, as strong T 1 relaxation observed below T C indicates fluctuating internal fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Solution-Based, Anion-Doping of Li 4 Ti 5 O 12 Nanoflowers for Lithium-Ion Battery Applications

Solution–based, anionic doping represents a convenient strategy with which to improve upon the conductivity of candidate anode materials such as Li 4 Ti 5 O 12 (LTO). As such, novel synthetic hydrothermally–inspired protocols have primarily been devised herein, aimed at the large–scale production of unique halogen–doped, micron–scale, three–dimensional, hierarchical LTO flower–like motifs. Although fluorine (F) doping has been explored, the use of chlorine (Cl) dopants is the primary focus here. Several experimental variables, such as dopant amount, lithium hydroxide concentration, and titanium butoxide purity, were probed and perfected. Furthermore, the Cl doping process did not damage the intrinsic LTO morphology. The analysis, based on interpreting a compilation of SEM, XRD, XPS, and TEM–EDS results, was used to determine an optimized dopant concentration of Cl. Electrochemical tests demonstrated an increased capacity via cycling of 12 % for a Cl–doped sample as compared with pristine LTO. Moreover, the Cl–doped LTO sample described in this study exhibited the highest discharge capacity yet reported at an observed rate of 2C for this material at 143mAh g –1 . Overall, these data suggest that the Cl dopant likely enhances not only the ion transport capabilities, but also the overall electrical conductivity of our as–prepared structures. Furthermore, to help explain these favorable findings, theoretical DFT calculations were used to postulate that the electronic conductivity and Li diffusion were likely improved by the presence of increased Ti 3+ ion concentration coupled with widening of the Li migration channel.

25 ENERGY STORAGE↗
DOE OSTI2021

Non-catalytic gas phase NO oxidation in the presence of decane

In this study, experiments without a catalyst revealed that NO was oxidized in a diesel exhaust gas phase mixture due to the presence of n -C 10 H 22 (decane). This reaction was observed to occur following the low temperature oxidation catalyst test protocol (LTC-D) defined by U.S. DRIVE. The purpose of LTC-D conditions is to simulate an aftertreatment diesel combustion gas mixture in order to test candidate catalyst materials. 100% NO conversion was observed, without a catalyst, after beginning to react at 330 °C accompanied by consumption of decane. After experiments that isolated hydrocarbons, ethylene was also observed to facilitate NO oxidation to a lesser degree (>470 °C). Density functional theory (DFT) calculations were conducted to investigate thermodynamically-possible initiating elementary steps during n -C 10 H 22 consumption and NO oxidation. Two feasible intermediate radicals to oxidize NO to NO 2 are ·C 10 H 21 O 2 and ·HO

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

Enhancing perovskite electrocatalysis through synergistic functionalization of B-site cation for efficient water splitting

The family of perovskite oxides is a promising class of catalysts for diverse energy conversion processes including water splitting. In this work, a facile two-step manipulation (in-situ exsolution and post-sulfurization) strategy was proposed and applied to LaCo 0.2 Fe 0.8 O 3 (LCF) perovskite parent, through which, the electronic state, spatial immersion and intrinsic activity of B-site cobalt (Co) were stepwise tuned at nanoscale proximity accordingly (i.e., lattice Co ions segregated Co 0 → embedded CoS 2 ). Impressively, the as-prepared catalyst (S-LCF) obtains an emergent oxygen deficient microstructure seamlessly pinned with uniformly distributed CoS 2 nanoparticles (NPs), which demonstrates enhanced performance toward both oxygen evolution reaction (OER) and hydrogen evolution reaction (HER), and shows good stability in overall water splitting. The density functional theory (DFT) calculations illustrate the optimized metal-oxygen covalency and hydrogen adsorption Gibbs free energy (Δ G H *) on S-LCF, which further buttresses the prominence of our B-site cation engineering tactics.

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

Synthesis and characterization of Pt(Cu 0.67 Sn 0.33 )

Pt(Cu 0.67 Sn 0.33 ) has recently been found in a natural sample. In order to be able to characterize this new ternary compound, we synthesized it from the elements. Samples were characterized by X-ray powder diffraction, differential scanning calorimetry, thermal relaxation calorimetry, and scanning electron microscopy studies. Density functional theory-based model calculations complemented the experimental studies. Pt(Cu 0.67 Sn 0.33 )was already formed at a relatively low temperature of 773 K. Rietveld refinement of Pt(Cu 0.67 Sn 0.33 ) has been carried out in CuAu-type or L1 0 -type structure, space group P4/mmm, with Pt on 0,0,0 and disordered Cu and Sn on 1/2, 1/2, 1/2 and Z = 1. The lattice parameters are a = 2.823(1) Å, c = 3.64(1) Å, and V = 29.00(4) Å which are in good agreement with values obtained earlier on the natural sample and with the results of DFT calculations. The vibrational entropy for Pt(Cu 0.67 Sn 0.33 ) is $S_{298.15}^{vib}$ = 79.9(7) J mol -1 K -1 . The pressure dependence up to 36(2) GPa of the unit-cell volume and the lattice parameters and unit-cell volume have been obtained by synchrotron based powder diffraction using a diamond anvil cell. A fit of a 3rd-order Birch–Murnaghan equation of state to the Pt(Cu 0.67 Sn 0.33 )) (p,V)-data results in a bulk modulus of B 0 = 215(27) GPa and B' = 5(2).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solar enhanced oxygen evolution reaction with transition metal telluride

The photo-enhanced electrocatalytic method of oxygen evolution reaction (OER) shows promise for enhancing the effectiveness of clear energy generation through water splitting by using renewable and sustainable source of energy. However, despite benefits of photoelectrocatalytic (PEC) water splitting, its uses are constrained by its low efficiency as a result of charge carrier recombination, a large overpotential, and sluggish reaction kinetics. Here, we illustrate that Nickel telluride (NiTe) synthesized by hydrothermal methods can function as an extremely effective photo-coupled electrochemical oxygen evolution reaction (POER) catalyst. In this study, NiTe was synthesized by hydrothermal method at 145°C within just an hour of reaction time. In dark conditions, the NiTe deposited on carbon cloth substrate shows a small oxygen evolution reaction overpotential (261 mV) at a current density of 10 mA cm –2 , a reduced Tafel slope (65.4 mV dec -1 ), and negligible activity decay after 12 h of chronoamperometry. By virtue of its enhanced photo response, excellent light harvesting ability, and increased interfacial kinetics of charge separation, the NiTe electrode under simulated solar illumination displays exceptional photoelectrochemical performance exhibiting overpotential of 165 mV at current density of 10 mA cm -2 , which is about 96 mV less than on dark conditions. In addition, Density Functional Theory investigations have been carried out on the NiTe surface, the results of which demonstrated a greater adsorption energy for intermediate -OH on the catalyst site. Since the -OH adsorption on the catalyst site correlates to catalyst activation, it indicates the facile electrocatalytic activity of NiTe owing to favorable catalyst activation. DFT calculations also revealed the facile charge density redistribution following intermediate -OH adsorption on the NiTe surface. This work demonstrates that arrays of NiTe elongated nanostructure are a promising option for both electrochemical and photoelectrocatalytic water oxidation and offers broad suggestions for developing effective PEC devices.

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

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

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

Multiple redox mechanisms for water-gas shift reaction on Fe 3 O 4 (1 1 1) surface: A density functional theory and mean-field microkinetic modeling study

In this work, we investigated redox mechanisms for the high temperature (HT) water–gas shift reaction (WGSR) on the Fe 3 O 4 (1 1 1) surface using density functional theory (DFT) calculations and mean-field microkinetic modeling (MKM). The redox pathways branched into three Langmuir-Hinshelwood processes (LH1, LH2 and LH3) and one Mars-van-Krevelen (MVK) process (in the presence of oxygen vacancy) from multidentate binding of CO and CO 2 over four available reactive sites (Fe oct2 , Fe tet1 , Fe bridge, and O1). We found that the LH1 and LH2 processes have CO adsorption at a single iron site (Fe oct2 or Fe tet1 ), while the LH3 and MVK processes have stronger chemisorption of CO or CO 2 by both the Fe oct2 site and the O 1 site. From the mechanistic study of these reaction path, we recognized that availability of O 1 sites was key to proceed to either CO oxidation by the single site (LH1, LH2) or by the dual site (LH3, MVK). We observed that the reaction energetics were significantly different in the CO oxidation steps, where the single or dual site results in markedly different apparent activation energy (73–281 kJ/mol) and reaction rates (10 -3 –10 -10 mol∙m -2 ∙s -1 ) among the four reaction mechanisms. The utilization of mean field MKM with DFT reaction energetics helps explain the experimental debates for the catalytic reaction details as well as providing a possible direction to engineer the catalyst with higher activity.

42 ENGINEERING↗
DOE OSTI2026

Quantifying the Influences of Epoxide Binding in Epoxide/CO 2 Ring Opening Copolymerization Catalysis

Understanding and predicting the effect of epoxide structure on the rate of polymerization in epoxide/CO 2 ring opening copolymerization catalysis is a long-standing challenge. Here, a known highly active Co(III)K(I) catalyst is used to investigate the influences of six different epoxides' binding strengths on their rates of copolymerization. Since calculations and experiments indicate that studying the catalytically relevant Co(III)−epoxide adduct directly is experimentally challenging, epoxide−catalyst binding interactions are quantified using a Co(II)K(I) complex to model the key catalytic intermediate. Epoxide−catalyst coordination is investigated using UV−vis spectroscopy titrations which provide fast and effective determination of association or binding constants. The epoxide−catalyst equilibrium constants show a clear exponential correlation with copolymerization rates and a new catalyst performance linear free energy relationship is revealed. Epoxides exhibiting stronger catalyst binding constants show higher copolymerization rates. The structure−activity correlation is consistent with the polymerization kinetics, mechanism and DFT calculations. Both the methods to investigate epoxide−catalyst coordination and the linear free energy relationship are shown to apply to the series of six epoxides and a second Co(III)K(I) catalyst. These structure− performance relationships are likely applicable to other transition metal catalysts and should expedite future epoxide and catalyst selection to make useful poly(carbonate) materials.

catalysts↗
DOE OSTI2025

Atomate2: modular workflows for materials science

High-throughput density functional theory (DFT) calculations have become a vital element of computational materials science, enabling materials screening, property database generation, and training of “universal” machine learning models. While several software frameworks have emerged to support these computational efforts, new developments such as machine learned force fields have increased demands for more flexible and programmable workflow solutions. This manuscript introduces atomate2, a comprehensive evolution of our original atomate framework, designed to address existing limitations in computational materials research infrastructure. Key features include the support for multiple electronic structure packages and interoperability between them, along with generalizable workflows that can be written in an abstract form irrespective of the DFT package or machine learning force field used within them. Our hope is that atomate2's improved usability and extensibility can reduce technical barriers for high-throughput research workflows and facilitate the rapid adoption of emerging methods in computational material science.

97 MATHEMATICS AND COMPUTING↗