Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covalency”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

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

Optimization of Random Phase Approximation Calculations for Improved Energies of Molecules, Solids, and Surfaces

We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.

Adsorption↗
DOE OSTI2020

Robust, Accurate, and Efficient: Quantum Embedding Using the Huzinaga Level-Shift Projection Operator for Complex Systems

Using wave function (WF) in density functional theory (DFT) embedding methods provides a framework for performing localized, high-accuracy WF calculations on a system, while not incurring the full computational cost of the WF calculation on the full system. To effectively partition a system into localized WF and DFT subsystems, we utilize the Huzinaga level-shift projection operator within an absolutely localized basis. In this work, we study the ability of the absolutely localized Huzinaga level-shift projection operator method to study complex WF and DFT partitions, including partitions between multiple covalent bonds, a double bond, and transition-metal–ligand bonds. We find that our methodology can accurately describe all of these complex partitions. Additionally, we study the robustness of this method with respect to the WF method, specifically where the embedded systems were described using a multiconfigurational WF method. We found that the method is systematically improvable with respect to both the number of atoms in the WF region and the size of the basis set used, with energy errors less than 1 kcal/mol. Additionally, we calculated the adsorption energy of H2 to a model of an iron metal–organic framework (Fe-MOF-74) to within 1 kcal/mol compared to CASPT2 calculations performed on the full model while incurring only a small fraction of the full computational cost. In conclusion, this work demonstrates that the absolutely localized Huzinaga level-shift projection operator method is applicable to very complex systems with difficult electronic structures.

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

Discovery of TAK-981, a First-in-Class Inhibitor of SUMO-Activating Enzyme for the Treatment of Cancer

SUMOylation is a reversible post-translational modification that regulates protein function through covalent attachment of small ubiquitin-like modifier (SUMO) proteins. The process of SUMOylating proteins involves an enzymatic cascade, the first step of which entails the activation of a SUMO protein through an ATP-dependent process catalyzed by SUMO-activating enzyme (SAE). Here, we describe the identification of TAK-981, a mechanism-based inhibitor of SAE which forms a SUMO–TAK-981 adduct as the inhibitory species within the enzyme catalytic site. Optimization of selectivity against related enzymes as well as enhancement of mean residence time of the adduct were critical to the identification of compounds with potent cellular pathway inhibition and ultimately a prolonged pharmacodynamic effect and efficacy in preclinical tumor models, culminating in the identification of the clinical molecule TAK-981.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Design, Synthesis, and Biological Activity of Novel Ornithine Decarboxylase (ODC) Inhibitors

We here describe the design, synthesis, and biological activity of novel ornithine decarboxylase (ODC) inhibitors that show significantly higher potency in vitro than α-difluoromethylornithine (DFMO), a U.S. Food and Drug Administration (FDA) approved drug. We report two X-ray structures of ODC complexed with new ODC inhibitors, computational docking, molecular dynamics, and binding free energy calculations to validate the experimental models. The X-ray structures reveal that covalent adducts with pyridoxal phosphate (PLP) are formed in the active site of the human ODC enzyme, as verified by their preparation and enzymatic testing. Finally, we verified that the cellular activity of endogenous ODC was inhibited, and polyamine levels were reduced. Given that ODC is a clinically validated target, combined with the fact that DFMO is currently the only ODC inhibitor in clinical use for several indications, the further development of more potent ODC inhibitors with superior activity and physical properties is warranted.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Balancing Charge Transfer and Frenkel Exciton Coupling Leads to Excimer Formation in Molecular Dimers: Implications for Singlet Fission

Photoexcitation of molecular chromophore aggregates can form excimer states that play a significant role in photophysical processes such as charge and energy transfer as well as singlet fission. An excimer state is commonly defined as a superposition of Frenkel exciton and charge transfer states. In this work, we investigate the dynamics of excimer formation and decay in π-stacked 9,10-bis(phenylethynyl)anthracene (BPEA) covalent dimers appended to a xanthene spacer, where the electronic coupling between the two BPEA molecules is adjusted by changing their longitudinal molecular slip distances. Using exciton coupling calculations, we quantify the relative contributions of Frenkel excitons and charge transfer states and find that there is an upper and lower threshold of the charge transfer contribution for efficient excimer formation to occur. In conclusion, knowing these thresholds can aid the design of molecular aggregates that optimize singlet fission.

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

Hydrolysis of Small Oxo/Hydroxo Molecules Containing High Oxidation State Actinides (Th, Pa, U, Np, Pu): A Computational Study

The energetics of hydrolysis reactions for high oxidation states of oxo/hydroxo monomeric actinide species (Th IV O 2 , Pa IV O 2 , U IV O 2 , Pa V O 2 (OH), U V O 2 (OH), U VI O 3 , Np VI O 3 , Np VII O 3 (OH), and Pu VII O 3 (OH)) were calculated at the CCSD(T) level. The first step is the formation of a Lewis acid/base adduct with H 2 O (hydration), followed by a proton transfer to form a dihydroxide molecule (hydrolysis); this process is repeated until all oxo groups are hydrolyzed. The physisorption (hydration) for each H 2 O addition was predicted to be exothermic, ca. -20 kcal/mol. The hydrolysis products are preferred energetically over the hydration products for the +IV and +V oxidation states. The compounds with An VI are a turning point in terms of favoring hydration over hydrolysis. For An VII O 3 (OH), hydration products are preferred, and only two waters can bind; the complete hydrolysis process is now endothermic, and the oxidation state for the An in An(OH) 7 is +VI with two OH groups each having one-half an electron. The natural bond order charges and the reaction energies provide insights into the nature of the hydrolysis/hydration processes. The actinide charges and bond ionicity generally decrease across the period. Here, the ionic character decreases as the oxidation state and coordination number increase so that covalency increases moving to the right in the actinide period.

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

Excited-State Bond Contraction and Charge Migration in a Platinum Dimer Complex Characterized by X-ray and Optical Transient Absorption Spectroscopy

Interactions between metal centers in dimeric transition metal complexes (TMCs) play important roles in their excited-state energetics and pathways, and, thus, affect their photophysical properties relevant to their applications, e.g., photoluminescent materials and photocatalysis. In this study, we report electronic and nuclear structural dynamics studies of two photoexcited pyrazolate-bridged [Pt(ppy)(µ-R 2 pz)] 2 -type Pt(II) dimers (ppy = 2-phenylpyridine, µ-R 2 pz = 3,5-substituted pyrazolate): [Pt(ppy)(µ-H 2 pz)] 2 (1) and [Pt(NDI-ppy)(µ-Ph 2 pz)] 2 (2, NDI = 1,4,5,8-naphthalenediimide), both of which have distinct ground state Pt-Pt distances. X-ray transient absorption spectroscopy at the Pt L III -edge revealed a new d-orbital vacancy due to the one-electron oxidation of the Pt centers in 1 and 2. However, while a transient Pt-Pt contraction was observed in 2, such an effect was completely absent in 1, demonstrating how the excited states of these complexes are determined by the overlap of the Pt (d z 2 ) orbitals, which is tuned by steric bulk of the pyrazolate R-groups in the 3- and 5-positions. In tandem with analysis of the Pt-Pt distance structural parameter, we observed photoinduced charge separation in 2 featuring a covalently linked NDI acceptor on the ppy ligand. The formation and subsequent decay of the NDI radical anion absorption signals were detected upon photoexcitation using optical transient absorption spectroscopy. The NDI radical anion decayed on the same timescale, hundreds of picoseconds, as that of the d-orbital vacancy signal of the oxidized Pt-Pt core observed in the XTA measurements. The data indicated ultrafast formation of the charge-separated state and subsequent charge recombination to the original Pt(II-II) species.

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

Multiconfiguration Pair-Density Functional Theory for Vertical Excitation Energies in Actinide Molecules

Modeling actinides with electronic structure theories is challenging because these systems present a strong ligand field and metal–ligand covalency. Here, we systematically investigate the effectiveness of pair-density functional theory (PDFT) for the calculation of vertical excitation energies in An(III), [An III Cl 6 ] 3– , and [An VI O 2 ] 2+ (An = U, Np, Pu, and Am). We compare the performance of PDFT, hybrid PDFT, and multistate PDFT with traditional active-space methods followed by perturbation theory, like multistate CASPT2, and with experimental data. Overall, multistate PDFT gives quantitative agreement with multistate CASPT2 at a significantly reduced computational cost.

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

Diabatic Valence-Hole Concept

A global diabatization scheme, based on the “valence-hole” concept, has been previously applied to model webs of avoided crossings that exist in four electronic-state symmetry manifolds of C 2 ( 1 Π g , 3 Π g , 1 Σ u + , and 3 Σ u + ). Here, this model is extended to the electronically excited states of four more molecules: CN ( 2 Σ + ), N 2 ( 3 Π u ), SiC ( 3 Π), and Si 2 ( 3 Π g ). Many strangenesses in the spectroscopic observations (e.g., energy level structure, predissociation linewidths, and radiative lifetimes) for all four electronic state systems discussed here are accounted for by this unified model. The key concept of the model is valence-hole electron configurations: 3σ 2 4σ 1 1π 4 5σ 2 in CN ( 2 Σ + ), 2σ g 2 2σ u 1 1π u 4 3σ g 2 1π g 1 in N 2 ( 3 Π u ), 5σ 2 6σ 1 7σ 2 2π 3 in SiC ( 3 Π), and 4σ g 2 4σ u 1 5σ g 2 2π u 3 in Si 2 ( 3 Π g ), all of which have a triply occupied “valence-core” (i.e., 2σ g 2 2σ u 1 or the equivalent). These valence-hole configurations have a nominal bond order of three or higher and correlate with high-energy separated-atom limits with an np ← ns (n = 2, 3) promotion in one of the atomic constituents. On its way to dissociation, the strongly bound diabatic valence-hole state crosses multiple weakly bound or repulsive states, which belong to electron configurations with a completely filled valence-core. These curve crossings between diabatic potentials result in a network of many avoided crossings among multiple electronic states, analogous to the well-studied electronic structure landscape of ionic-covalent crossings in strongly ionic molecules. In conclusion, considering the unique role of valence-hole states in shaping the global electronic structure, the valence-hole concept should be added to our intuitive framework of chemical bonding.

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

Energetic and Electronic Properties of AcX and LaX (X = O and F)

The bonding and spectroscopic properties of LaX and AcX (X = O and F) diatomic molecules were studied by highlevel ab initio CCSD(T) and SO-CASPT2 electronic structure calculations. Bond dissociation energies (BDEs) were calculated at the Feller−Peterson−Dixon (FPD) level. Potential energy curves and spectroscopic constants for the lowest-lying spin−orbit Ω states were obtained at the SO-CASPT2/aQ-DK level. A dense manifold of excited states was described for the monofluorides with the ground states well separated from the excited states. The spectroscopic parameters were in good agreement with those reported experimentally for LaO and LaF. For the diatomic molecules containing actinides, no experimental data of these parameters was found, but the results were consistent with other high-level calculations. The BDEs calculated at the FPD level were 791.3 (LaO), 705.2 (AcO), 650.0 (LaF), and 678.6 (AcF) kJ/ mol. The NBO analysis showed that the monofluorides are essentially ionic, which explains why the BDE(AcF) is higher than BDE(LaF); for the monoxides, covalent contributions involving the d orbitals of the metal and the p orbitals of the oxygen are stronger for LaO than AcO, which explains the higher BDE for LaO. Finally, the bond orders are predicted to be 2 for LaF and AcF, 3 for AcO, and higher than 3 for LaO.

Basis sets↗
DOE OSTI2025

Ultrafast Correlation Energy Estimator

A virtually no-cost method is proposed that can compute the correlation energies of general, covalently bonded, organic, and inorganic molecules (including conjugated π-electron systems) with a well-defined dominant Lewis structure at the accuracy of 99.5% of the near-exact values determined by the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] in the complete-basis-set (CBS) limit. This Correlation Energy Per Bond (CEPB) method assigns a partial correlation energy to each bond type (characterized by the identities of the two atoms forming the bond and its integer bond order) and to a lone pair, regardless of the bond length, bond angle, sp-hybridization, π-electron conjugation, ionicity, noncovalent interactions, etc. At its current stage, the method is mainly suitable for near-equilibrium geometries. The correlation energies per bond are determined by a fit to the CCSD(T)/CBS benchmarks. It can neither improve the equilibrium structures nor discern conformers or positional isomers, yet its accuracy for reaction energies rivals that of the second-order Møller–Plesset perturbation theory, which is far more expensive. Its promising performance underscores the possibility that surprisingly compact, chemically intuitive molecular fragments exist into which correlation energies can be partitioned, leading to various ultrafast correlation-energy estimators tailored to different purposes.

Binding energy↗
DOE OSTI2022

Influence of the Nature of Aliphatic Hydrophobic Physical Crosslinks on Water Crystallization in Copolymer Hydrogels

The local environment within a hydrogel influences the properties of water, including the propensity for ice crystallization. Water-swollen amphiphilic copolymers produce tunable nanoscale environments, which are defined by hydrophobic associations, for the water molecules. In this work, the antifreeze properties for equilibrium-swollen amphiphilic copolymers with a common hydrophilic component, hydroxyethyl acrylate (HEA), but associated through crystalline (octadecyl acrylate, ODA) or rubbery (ethylhexyl acrylate, EHA) hydrophobic segments, are examined. Differences in the efficacy of the associations can be clearly enunciated from compositional solubility limits for the copolymers in water (<2.6 mol % ODA vs ≤14 mol % EHA), and these differences can be attributed to the strength of the association. The equilibrium-swollen HEA-ODA copolymers are viscoelastic solids, while the swollen HEA-EHA copolymers are viscoelastic liquids. Cooling these swollen copolymers to nearly 200 K induces some crystallization of the water, where the fraction of water frozen depends on the details of the nanostructure. Decreasing the mean free path of water by increasing the ODA composition from 10 to 25 mol % leads to fractionally more unfrozen water (66–87%). The swollen HEA-EHA copolymers only marginally inhibit ice (<13%) except with 45 mol % EHA, where nearly 60% of the water remains amorphous on cooling to 200 K. In general, the addition of the EHA leads to less effective ice inhibition than analogous covalently crosslinked HEA hydrogels (19.9 ± 1.8%). These results illustrate that fluidity of confining surfaces can provide pathways for critical nuclei to form and crystal growth to proceed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerated Sequence Design of Star Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning

Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatibilization, chemical transformations, and separations. s-BCPs have chain architectures where three or more linear diblock copolymer arms comprised of two chemically distinct linear polymers, e.g., solvophobic and solvophilic chains, are covalently joined at one point. The chemical composition of each of the subunit polymer chains comprising the arms, their molecular weights, and the number of arms can be varied to tailor the surface and interfacial activity of these architecturally unique molecules. Further, this makes identification of the optimal s-BCP design nontrivial as the total number of plausible s-BCP architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with a reinforcement learning-based Monte Carlo tree search (MCTS) to identify s-BCP designs that minimize the interfacial tension between polar and nonpolar solvents. We first validate the MCTS approach for the design of small- and medium-sized s-BCPs and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified by using the configurations obtained from MD simulations. Chemical insights into the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provides important groundwork for future experimental investigation of s-BCPs for various applications.

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

RbB 3 Si 3 : An Alkali Metal Borosilicide that is Metastable and Superconducting at 1 atm

The stability, electronic structure and potential superconductivity in AB 3 Si 3 (A = Na, K, Rb, and Cs) compounds that assume the clathrate-based sodalite structure whose frameworks consists of covalent B-Si bonds is investigated via first-principles calculations. This structure type has recently been predicted in a number of high temperature superconducting hydrides, but these are only stable under megabar pressures. Herein, we predict a novel superconducting phase, RbB 3 Si 3 , that could be synthesized under pressures that are a factor of ten smaller, ~10 GPa, and quenched to atmospheric conditions. Electron-phonon coupling calculations predict that RbB 3 Si 3 possesses a superconducting critical temperature, T c , of 14 K at 1 atm. Here, the dynamic stability of RbB 3 Si 3 and CsB 3 Si 3 at ambient pressure can be explained by considering the chemical pressure exerted on the B-Si framework that is caused by the size effect of the alkali metal atom.

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

Interplay between Electrostatic Properties of Molecular Adducts and Their Positions at Carbon Nanotubes

Formation of the sp 3 -defect due to covalent functionalization of a single-walled carbon nanotube (SWCNT) results in a new red-shifted emissive excitons. Here, using density functional theory, we study the impact of the intrinsic molecular dipole moments and the local charges induced by various molecular adducts on the energy and localization properties of the optically active defect-based exciton in a (10,5) SWCNT. The interplay of two effects plays a role in the localization of the exciton and, therefore, its red shift, but at different scales: The source of the leading order contribution is the defect conformation, resulting in the red shift of the defect-associated exciton with respect to the E 11 band of the pristine SWCNT on the order of ~100 meV, while the individual dipoles and polarization properties of molecular adducts lead to significantly smaller red shifts on the order of ~10 meV. While the species-dependent trends in defect-induced charges do not directly correlate to the exciton red shift, the charge at the sp 3 -defect exhibits a distinct behavior between ortho- and para-defect configurations and is relevant to the chemical reactivity of the defect position depending on the adduct type. Overall, our computational results may be helpful for diverse synthetic strategies to fine-tune emission of SWCNTs toward desired applications.

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

Suppressing Auger Recombination in Multiply Excited Colloidal Silicon Nanocrystals with Ligand-Induced Hole Traps

Nonradiative Auger recombination in multiply excited nanocrystals is a dominant efficiency loss pathway for nanocrystal-containing optoelectronic devices that rely on high-rate emission and absorption operating conditions. Overcoming Auger recombination in these quantum-confined systems is therefore a longstanding challenge to the synthetic nanocrystal as well as device manufacturing communities. Several successful strategies have been realized to reduce Auger recombination, but they rely on complex and time-consuming nanocrystal core/shell synthesis. Alternatively, controlling Auger rates by varying the nanocrystal–ligand-binding chemistry is a promising route to obtain functional tunability, which reduces the barrier to large-scale manufacturing. The covalent surface chemistry and the extremely long-lived photoexcited lifetimes of silicon nanocrystals (Si NCs) make them a unique system among colloidal semiconductor NCs to study the intersection between surface chemistry and photoexcited carrier dynamics. Here, we show that changing the functional group that binds a saturated dodecyl ligand to the surface of nonthermal plasma-synthesized Si NCs from alkyl to thiolate slows Auger recombination rates within multiply excited Si NCs. This reduction in Auger rate persists across Si NC sizes ranging from 3.5 to 8 nm in diameter, but the expected linear dependence of Auger rates on the NC volume is retained for both alkyl and alkylthiolate surface terminations. To understand the origin behind this elongation, we carry out steady-state and time-resolved photoluminescence measurements as well as time-resolved terahertz spectroscopy measurements. These measurements reveal that thiolate groups introduce mid-gap surface states, which, we argue, reduces the photoexcited electron–hole overlap and elongates Auger recombination times. These results highlight how a typically detrimental chemical species—mid-band gap NC surface states—can be beneficial under high-rate absorption/emission conditions.

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

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗