Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ALKANE”

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 649 records · Page 36

Synthetic and Mechanistic Implications of Chlorine Photoelimination in Nickel/Photoredox C(sp 3 )–H Cross-Coupling

In recent years, the development of light-driven reactions has contributed numerous advances in synthetic organic chemistry. A particularly active research area combines photoredox catalysis with nickel catalysis to accomplish otherwise inaccessible cross-coupling reactions. In these reactions, the photoredox catalyst absorbs light to generate an electronically excited charge-transfer state that can engage in electron or energy transfer with a substrate and the nickel catalyst. Our group questioned whether photoinduced activation of the nickel catalyst itself could also contribute new approaches to cross-coupling. Over the past 5 years, we have sought to advance this hypothesis for the development of a suite of mild and site-selective C(sp 3 )–H cross-coupling reactions with chloride-containing coupling partners via photoelimination of a Ni–Cl bond. On the basis of a report from the Nocera laboratory, we reasoned that photolysis of a Ni(III) aryl chloride species, generated by single-electron oxidation of a typical Ni(II) intermediate in cross-coupling, might allow for the catalytic generation of chlorine atoms. Combining this with the ability of Ni(II) to accept alkyl radicals, we hypothesized that photocatalytically generated chlorine atoms could mediate hydrogen atom transfer (HAT) with C(sp 3 )–H bonds to generate a substrate-derived alkyl radical that is captured by the Ni center in cross-coupling. A photoredox catalyst was envisioned to promote the necessary single-electron oxidation and reduction of the Ni catalyst to facilitate an overall redox-neutral process. Overall, this strategy would offer a visible-light-driven mechanism for chlorine radical formation enabled by the sequential capture of two photons. As an initial demonstration, we developed a Ni/photoredox-catalyzed α-oxy C(sp 3 )–H arylation of cyclic and acyclic ethers. This method was extended to a mild formylation of abundant and complex aryl chlorides through selective 2-functionalization of 1,3-dioxolane. Seeking to develop a suite of reactions that introduce carbon at all different oxidation states, we explored C(sp 3 )–H cross-coupling with trimethyl orthoformate, a common laboratory solvent. We found that trimethyl orthoformate serves as a source of methyl radical for a methylation reaction via β-scission from a tertiary radical generated upon chlorine-mediated HAT. Since chlorine radical is capable of abstracting unactivated C(sp 3 )–H bonds, our efforts have also been directed at cross-coupling with a range of feedstock chemicals, such as alkanes and toluenes, along with late-stage intermediates, using chloroformates as coupling partners. Overall, this platform enables access to valuable synthetic transformations with (hetero)aryl chlorides, which despite being the most ubiquitous and inexpensive aryl halide coupling partners, are rarely reactive in Ni/photoredox catalysis. Little is known about the photophysics and photochemistry of organometallic Ni complexes relevant to cross-coupling. We have conducted mechanistic investigations, including computational, spectroscopic, emission quenching, and stoichiometric oxidation studies, of Ni(II) aryl halide complexes common to Ni/photoredox reactions. These studies indicate that chlorine radical generation from excited Ni(III) is operative in the described C(sp 3 )–H functionalization methods. More generally, the studies illustrate that the photochemistry of cross-coupling catalysts cannot be ignored in metallaphotoredox reactions. We anticipate that further mechanistic understanding should facilitate catalyst design and lead to the development of new synthetic methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Story of 5d Metallocorroles: From Metal–Ligand Misfits to New Building Blocks for Cancer Phototherapeutics

Porphyrin chemistry is Shakespearean: over a century of study has not withered the field’s apparently infinite variety. Heme proteins continually astonish us with novel molecular mechanisms, while new porphyrin analogues bowl us over with unprecedented optical, electronic, and metal-binding properties. Within the latter domain, corroles occupy a special place, exhibiting a unique and rich coordination chemistry. The 5d metallocorroles are arguably the icing on that cake. New Zealand chemist Penny Brothers has used the word “misfit” to describe the interactions of boron, a small atom with a predilection for tetrahedral coordination, and porphyrins, classic square-planar ligands. Steve Jobs lionized misfits as those who see things differently and push humanity forward. Both perspectives have inspired us. The 5d metallocorroles are misfits in that they encapsulate a large 5d transition metal ion within the tight cavity of a contracted porphyrin ligand. Given the steric mismatch inherent in their structures, the syntheses of some 5d metallocorroles are understandably capricious, proceeding under highly specific conditions and affording poor yields. Three broad approaches may be distinguished. (a) In the metal–alkyl approach, a free-base corrole is exposed to an alkyllithium and the resulting lithio-corrole is treated with an early transition metal chloride; a variant of the method eschews alkyllithium and deploys a transition metal–alkyl instead, resulting in elimination of the alkyl group as an alkane and insertion of the metal into the corrole. This approach is useful for inserting transition metals from groups 4, 5, and, to some extent, 6, as well as lanthanides and actinides. (b) In our laboratory, we have often deployed a low-valent organometallic approach for the middle transition elements (groups 6, 7, 8, and 9). The reagents are low-valent metal–carbonyl or −olefin complexes, which lose one or more carbon ligands at high temperature, affording coordinatively unsaturated, sticky metal fragments that are trapped by the corrole nitrogens. (c) Finally, a metal acetate approach provides the method of choice for gold and platinum insertion (groups 10 and 11). This Account provides a first-hand perspective of the three approaches, focusing on the last two, which were largely developed in our laboratory. In general, the products were characterized with X-ray crystallography, electrochemistry, and a variety of spectroscopic methods. The physicochemical data, supplemented by relativistic DFT calculations, have provided fascinating insights into periodic trends and relativistic effects. An unexpected feature of many 5d metallocorroles, given their misfit character, is their remarkable stability under thermal, chemical, and photochemical stimulation. Many of them also exhibit long triplet lifetimes on the order of 100 μs and effectively sensitize singlet oxygen formation. Many exhibit phosphorescence in the near-infrared under ambient conditions. Furthermore, water-soluble ReO and Au corroles exhibit impressive photocytotoxicity against multiple cancer cell lines, promising potential applications as cancer phototherapeutics. We thus envision a bright future for the compounds as rugged building blocks for new generations of therapeutic and diagnostic (theranostic) agents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bimetallic-Derived Catalytic Structures for CO 2 -Assisted Ethane Activation

In recent years, the simultaneous upgrading of CO 2 and ethane has emerged as a promising approach for generating valuable gaseous (CO, H 2 , and ethylene) and liquid (aromatics and C 3 oxygenates) chemicals from the greenhouse gas CO 2 and large-reserved shale gas. The key challenges for controlling product selectivity lie in the selective C-H and C-C bond cleavage of ethane with the assistance of CO 2 . Bimetallic-derived catalysts likely undergo alloying or oxygen-induced segregation under reaction conditions, thus providing diverse types of interfacial sites, e.g., metal/support (M/M'O x ) interface and metal oxide/metal (M'O x /M) inverse interface, that are beneficial for selective CO 2 -assisted ethane upgrading. The alloying extent can be initially predicted by cohesive energy and atomic radius (or Wigner-Seitz radius), while the preference for segregation to form the on-top suboxide can be approximated using the work function, electronegativity, and binding strength of adsorbed oxygen. Furthermore, bimetallic-derived catalysts are typically supported on high surface area oxides. Modifying the reducibility and acidity/basicity of the oxide supports and introducing surface defects facilitate CO 2 activation and oxygen supplies for ethane activation. Using in situ synchrotron characterization and density functional theory (DFT) calculations, we found that the electronic properties of oxygen species influence the selective cleavage of C-H/C-C bonds in ethane, with electron-deficient oxygen over the metal (or alloy) surface promoting nonselective bond scission to produce syngas and electron-enriched oxygen over the metal oxide/metal interface enhancing selective C-H scission to yield ethylene. Here, we further demonstrate that the preferred structures of the catalyst surfaces, either alloy surfaces or metal oxide/metal inverse interfaces, can be controlled through the appropriate choice of metal combinations and their atomic ratios. Through a comprehensive comparison of experimental results and DFT calculations, the selectivity of C-C/C-H bond scission is correlated with the thermodynamically favorable bimetallic-derived structures (i.e., alloy surfaces or metal oxide/metal inverse interfaces) under reaction conditions over a wide range of bimetallic catalysts. These findings not only offer structural and mechanistic insights into bimetallic-derived catalysts but also provide design principles for selective catalysts for CO 2 -assisted activation of ethane and other light alkanes. This Account concludes by discussing challenges and opportunities in designing advanced bimetallic-derived catalysts, incorporating new reaction chemistries for other products, employing precise synthesis strategies for well-defined structures with optimized site densities, and leveraging time/spatial/energy-resolved in situ spectroscopy/scattering/microscopy techniques for comprehensive structural analysis. The research methodologies established here are helpful for the investigation of dynamic alloy and interfacial structures and should inspire more efforts toward the simultaneous upgrading of CO 2 and shale gas.

03 NATURAL GAS↗

Using a Chromatographic Pseudophase Model To Elucidate the Mechanism of Olefin Separation by Silver(I) Ions in Ionic Liquids

Silver(I) ions undergo selective olefin complexation and have been utilized in various olefin/paraffin separation techniques such as argentation chromatography and facilitated transport membranes. Ionic liquids (ILs) are solvents known for their low vapor pressure, high thermal stability, low melting points, and ability to promote a favorable solvation environment for silver(I) ion–olefin interactions. To develop highly selective separation systems, a fundamental understanding of analyte partitioning to the stationary phase and the thermodynamic driving forces behind solvation is required. In this study, a chromatographic model treating silver(I) ions as a pseudophase is constructed and employed for the first time to investigate the olefin separation mechanism in silver(I) salt/IL mixtures. Stationary phases containing varying amounts of noncoordinated silver(I) salt ([Ag + ][NTf 2 – ]) dissolved in the 1-decyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 10 MIM + ][NTf 2 – ]) IL are utilized to determine the partition coefficients of various analytes including alkanes, alkenes, alkynes, aromatics, aldehyde, esters, and ketones. As ligand coordination to silver(I) ions is known to lower its olefin complexation capability, this study also examines two different types of coordinated silver(I) ion pseudophases, namely, monocoordinated silver(I) salt ([Ag + (1-decyl-2-methylimidazole, DMIM)][NTf 2 – ]) and dicoordinated silver(I) salt ([Ag + (1-methylimidazole, MIM)(DMIM)][NTf 2 – ]). The extent of olefin partitioning to the coordinated silver(I) ion pseudophases over the carrier gas and IL decreased by up to two orders of magnitude. Values for enthalpy, entropy, and free energy of solvation were determined for the three silver(I) ion-containing systems. Olefin retention was observed to be enthalpically dominated, while ligand coordination to the silver(I) ion pseudophase resulted in variations for both enthalpic and entropic contributions to the free energy of solvation. Furthermore, the developed model can be used to study chemical changes that occur in silver(I) ions over time as well as identify optimal silver(I) salt/IL mixtures that yield high olefin selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Confocal Raman Microscopy Investigation of Self-Assembly of Hybrid Phospholipid Bilayers within Individual Porous Silica Chromatographic Particles

Hybrid-supported phospholipid bilayers are a model structure utilized for measurement of molecular interactions that typically occur at cell membranes. These membrane models are prepared by adsorption of a lipid monolayer onto a stable n -alkyl chain layer that is covalently-bound to a support surface. Hybrid bilayers have been adapted to chromatographic retention measurements of lipophilicity through the assembly of a phospholipid monolayer onto n -alkane-modified silica surfaces in reversed-phase chromatographic particles. Recent Raman microscopy studies of these particles have shown that the acyl chains of the phospholipid interact with the C 18 -alkyl chains immobilized on the silica surface, where both lipid and C 18 alkyl chains become ordered due to chain interdigitation. Confocal Raman microscopy has also been used to investigate the association of small molecules with hybrid-lipid bilayers in C 18 chromatographic silica particles; the partitioning of model solutes compares favorably to that in lipid vesicle membranes with similar changes in acyl-chain structure (disordering) with solute partitioning. The present study seeks information about how these membrane-mimetic bilayers assemble onto the C 18 –derivatized silica surfaces of reversed-phase chromatographic silica particles. Confocal Raman microscopy is capable of interrogating the time-dependent internal composition and structure within individual silica particles. The Raman scattering data can be resolved into component Raman spectra and corresponding composition vectors that describe the time-dependent changes in intensity of the component spectra. This analysis provides insight into how the structures of both the lipid and C 18 alkyl chains of hybrid lipid bilayers evolve during deposition and organization on the internal surfaces of reversed-phase chromatographic silica particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catalytic Applications of Vanadium: A Mechanistic Perspective

The chemistry of vanadium has seen remarkable activity in the past SO years. In the present review, reactions catalyzed by homogeneous and supported vanadium complexes from 2008 to 2018 are summarized and discussed. Here particular attention is given to mechanistic and kinetics studies of vanadium-catalyzed reactions including oxidations of alkanes, alkenes, arenes, alcohols, aldehydes, ketones, and sulfur species, as well as oxidative C-C and C-O bond cleavage, carbon carbon bond formation, deoxydehydration, haloperoxidase, cyanation, hydrogenation, dehydrogenation, ring-opening metathesis polymerization, and oxo/imido heterometathesis. Additionally, insights into heterogeneous vanadium catalysis are provided when parallels can be drawn from the homogeneous literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Selected High-Performance Fuel Blends on Three-Way Catalyst Light Off under Synthetic Spark-Ignition Engine-Exhaust Conditions

The U.S. Department of Energy funded Co-Optimization of Fuels and Engines initiative aims to simultaneously develop advanced engines along with high-performance fuels to reduce petroleum consumption. The engine exhaust of spark-ignited light-duty vehicles contains pollutants such as nitrogen oxides, carbon monoxide (CO), and non-methane organic gases. When operated above their “light-off” temperature, three-way catalysts (TWCs) efficiently control the emissions of these pollutants from the vehicle exhaust. However, below the catalyst light-off temperature, during cold start, the TWCs are not effective. Thus, the stringent environmental regulations necessitate cold-start compliance of advanced engines operating on novel fuels for commercialization. Exhaust composition strongly impacts the effectiveness of TWCs. Hence, ensuring that the high-performance fuels under consideration do not have detrimental effects on current emissions control technology is necessary. To mitigate cold-start emissions, a low light-off temperature of the fuel on the TWC is desirable. As real-world fuels are multicomponent blends, we conducted investigations into the light-off behavior of representative fuel mixtures on a three-mode redox-aged commercial TWC under synthetic engine-exhaust conditions. The high-performance fuels in this study included 10–30% volumetric levels of ethanol, isobutanol, diisobutylene, and an aromatic mixture. Each of these high-performance fuel components was mixed into a gasoline surrogate blendstock for oxygenate blending (BOB). Our results showed that aromatics and alkenes in the surrogate BOB inhibit low-temperature reactivity of alkanes, alcohols, and CO on the TWC and dominate the blend light-off behavior. All the high-performance fuel blends had a very similar light-off behavior to the surrogate gasoline BOB, indicating that blending up to 30% (by vol.) of high-performance blendstocks in a gasoline base fuel can potentially reduce greenhouse gas emissions through improved engine efficiency and petroleum displacement without jeopardizing the ability to meet emissions regulations. While some high-performance blendstocks demonstrated lower light-off temperatures than a surrogate gasoline blend, taking advantage of the higher catalytic reactivity of these blendstocks to reduce cold-start emissions would require reducing the aromatic content in petroleum-based market fuels.

33 ADVANCED PROPULSION SYSTEMS↗

Bio-Oil Impact on Water Diffusion and Durability of Bitumen: Influence of Aging and Salinity

Bio-oils derived from inexpensive biomass offer environmentally friendly options for modifying bitumen for improved durability or rejuvenating aged material, although their impact on bitumen moisture resistance can be mixed. A better understanding of water diffusion into bitumen and its interactions with bio-oil compounds would aid the development of effective bio-oils. In this study, the effect of water exposure on bio-oil-modified bitumen blends was examined by de-wetting contact angle measurements and differential FTIR spectroscopy. The bio-oils improved the anti-stripping behavior of the bitumen but also increased water absorption which could weaken cohesive strength. Short-term thermal aging increased water diffusion for all bitumen blends including the control, probably due to the presence of more oxidized compounds. Negative peaks of alkanes and polar groups in the differential FTIR data suggest the co-diffusion of surfactant molecules towards the bitumen-water interface. Basic pH increased de-wetting of some bitumen blends from silica, possibly by attacking the silica surface itself. The presence of salt in solution altered the bitumen surface composition through the formation of salt complexes with bitumen or bio-oil compounds. In particular, the formation of calcium-carboxylate complexes appeared to greatly improve anti-stripping effects of bio-oils. Water diffusion into most of the bitumen blends was insensitive to pH or salt concentration except for a few notable outliers. Identifying the compounds responsible for increasing or decreasing water diffusion in these outlier cases would be valuable developing future bio-oil formations that avoid or promote those compounds.

Hung, Albert M.↗

Hierarchical Galloaluminosilicate MFI Catalysts for Ethane Nonoxidative Dehydroaromatization

This study presented an approach to overcome diffusion limitation and imbalance between Brønsted and Lewis acidities that affect conversion and aromatic yield in alkane nonoxidative dehydroaromatization (DHA) conversion. Three different Ga-promoted ZSM-5 catalysts were synthesized for the ethane DHA reaction. The hierarchical galloaluminosilicate Meso-GaZSM-5 was prepared by the solid-state crystallization technique wherein Ga was introduced through the in situ synthesis process. The microporous galloaluminosilicate Micro-GaZSM-5 was synthesized by the hydrothermal technique wherein Ga was introduced in situ during zeolite formation. The hierarchically structured Ga-Meso-ZSM-5 was synthesized by the solid-state crystallization method, but the Ga was introduced by the incipient wetness method. Additional Pt was incorporated by the incipient wetness technique for all three catalysts to enhance the DHA activity. The solid-state crystallization technique created the hierarchical structure without using a template. With the in situ technique, Ga species replaced the Al framework, reducing the strong Brønsted acidity. Meanwhile, they created the extraframework GaO + , introducing the strong Lewis acid sites into catalysts. Due to better synergistic interaction between the oxidation state of metals, acidic properties, and catalyst porosity, Pt-Meso-GaZSM-5 showed higher C 2 H 6 conversion and aromatic selectivity. The improvement on diffusivity not only reduces the effect of channel blockage but also improves the utilization of active sites. As a result, this work sheds light on the direct, nonoxidative conversion of shale gas without going through the syngas route.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Natural Attenuation in Streambed Sediment Receiving Chlorinated Solvents from Underlying Fracture Networks

Contaminant discharge from fractured bedrock formations remains a remediation challenge. Here, we applied an integrated approach to assess the natural attenuation potential of sediment that forms the transition zone between upwelling groundwater from a chlorinated solvent-contaminated fractured bedrock aquifer and the receiving surface water. In situ measurements demonstrated that reductive dechlorination in the sediment attenuated chlorinated compounds before reaching the water column. Microcosms established with creek sediment or in situ incubated Bio-Sep beads degraded C 1 -C 3 chlorinated solvents to less-chlorinated or innocuous products. Quantitative PCR and 16S rRNA gene amplicon sequencing revealed the abundance and spatial distribution of known dechlorinator biomarker genes within the creek sediment and demonstrated that multiple dechlorinator populations degrading chlorinatedC 1 -C 3 alkanes and alkenes co-inhabit the sediment. Phylogenetic classification of bacterial and archaeal sequences indicated a relatively uniform distribution over spatial (300 m horizontally) scale, but Dehalococcoides and Dehalobacter were more abundant in deeper sediment, where 5.7 ± 0.4 × 10 5 and 5.4 ± 0.9 × 10 6 16S rRNA gene copies per g of sediment, respectively, were measured. The microbiological and hydrogeological characterization demonstrated that microbial processes at the fractured bedrock-sediment interface were crucial for preventing contaminants reaching the water column, emphasizing the relevance of this critical zone environment for contaminant attenuation.

58 GEOSCIENCES↗

Predicting Aerosol Reactivity Across Scales: from the Laboratory to the Atmosphere

To fully utilize the results of laboratory-based studies of the chemistry of model atmospheric aerosol reactions, it is important to understand how to relate them to the conditions found in nature. In this study, we have taken a validated reaction-diffusion mechanism for oxidation of C 30 H 62 aerosol by OH under flow tube conditions and examined its predictions for another experimental regime (continuous flow stirred tank reactor) and for the atmosphere, spanning alkane aerosol viscosities from liquid to semisolid. The results show that under OH-concentration-limited and aerosol-mixing-limited conditions, it should be possible to select laboratory experimental conditions where many aspects of the particle phase and volatile product chemistry under atmospheric conditions can be revealed. If the OH collision and organic diffusion rates are comparable, however, reactivity is highly sensitive to the details of both OH concentration and internal mixing. Furthermore the characteristics of the transition between limiting conditions provide key insights into which parts of the reaction mechanism dominate in the various kinetic regimes.

54 ENVIRONMENTAL SCIENCES↗

Analysis of Tar and Oil Derived from Pyrolysis and Copyrolysis of Waste Plastics and Biomass

Pyrolysis has been proposed as a potential technology for managing the growing volume of plastic waste generated worldwide. Co-pyrolysis of plastic waste with biomass is a promising technology for generating fuel and chemical products. However, this process generates tar as a waste product. The chemical properties of this tar have yet to be thoroughly analyzed. Further, this study presents the results of gas chromatography–mass spectrometry (GC–MS), Fourier-transform infrared spectroscopy (FTIR), and thermogravimetric analysis (TGA) of oil and tar obtained from the pyrolysis of pure plastics including high-density polyethylene (HDPE), low-density polyethylene (LDPE), polyethylene (PE), polystyrene (PS), and plastic-biomass mixtures. GC–MS analysis revealed the presence of C 7 –C 37 carbon-containing hydrocarbons, which include alkanes and alkenes as the dominant products. FTIR data revealed the presence of various functional groups, including alcohols, aldehydes, ketones, and carboxylic acids, indicating the complexity of the pyrolysis and copyrolysis oil obtained from waste plastics and biomass. TGA data show that tar from all four plastics has a higher decomposition rate, suggesting the presence of heavier hydrocarbons compared with their corresponding oils. This research will be of interest to researchers looking to advance the study of plastic and biomass waste management.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Trickier than It Looks: Isomerization between Five- and Six-Coordinated Zinc in Heterometallic Li 2 Zn 2 Molecule

This report describes the synthesis and characterization of two hetero bi metallic Li–Zn coordination isomers [Li 2 Zn 2 (tbaoac) 6 ] (tbaoac = tert -butyl acetoacetato) that have been isolated separately by the same stoichiometric reaction run in different organic solvents. The 6-coordinated zinc isomer ( 6-Zn ) was synthesized in acetone with high yield, while the 5-coordinated one ( 5-Zn ) was readily obtained from ethanol. The 5-Zn isomer has a low solubility in organic solvents such as alkanes and haloalkanes, while its 6-Zn counterpart exhibits a good solubility in almost all common solvents. Two isomeric molecules feature similar centrosymmetric tetranuclear cyclic assemblies, which are different in their arrangement of tbaoac ligands. While all ligands act as μ 2 -type in the structure of 5-Zn , the two tbaoac groups chelating Li appear as μ 3 -type in 6-Zn , thus providing an additional coordination for Zn ions. However, the real structural transformation between these isomers was shown to be more complex than simply making or breaking a couple of Zn–O bonds. X-ray single-crystal structure analysis, powder X-ray diffraction, multinuclear NMR, DART mass spectrometry, ICP-OES analysis, and TGA have been employed for the characterization of the isomers. The combination of powder X-ray diffraction and 1 H NMR investigation revealed that 6-Zn isomer can be quantitatively transformed to 5-Zn in ethanol, while the reverse conversion instantly takes place in acetone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Simulation of the Catalytic Regeneration of n BuLi through a Hydrometalation Route

Efficient regeneration of organolithium compounds is a challenging aspect in the process of novel organometathetical catalytic cycles. A new method for Mg production from seawater has been proposed that capitalizes on the rich chemistry of Grignard reagents. Based on a 3-step catalytic cycle with Cp2MClL catalyst (M=Ti, Zr and L select organic ligands), it requires the regeneration of nBuLi. The potential of the proposed approach is evaluated with density functional theory based molecular simulations. The results reveal that the high affinity of Li towards Cl and N results in the formation of alkanes, and the strong coupling of catalyst and BuLi leads to catalyst deactivation. To improve its catalytic performance, we proposed the use of a diamine co-catalyst and a modified catalyst with a ligand that does not contain N, which would help release BuLi from the vicinity of the catalyst. Ab inito molecular dynamics simulations at 298 K in explicit solvent (THF) were used to estimate the Gibb’s free energetics and equilibrium constants obtained from the vibrational density of states. The results show a marked improvement in the free-energetics with lower barriers towards the completion of the catalytic cycle and suppression of de-activation channels. M.-S. L. and B. P. M. were supported by the US Department of Energy (DOE), Advanced Research Projects Agency for Energy performed at the Pacific Northwest National Laboratory (PNNL). V.-A. G. and R. R. were supported by DOE, Office of Science, Office of Basic Energy Sciences-CGBS Catalysis Program. Computational resources were provided by PNNL Research Computing Cluster and the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231.

Lee, Mal-Soon (ORCID:000000016851177X)↗

Open-Shell Tensor Hypercontraction

The extension of least-squares tensor hypercontracted second- and third-order Møller–Plesset perturbation theory (LS-THC-MP2 and LS-THC-MP3) to open-shell systems is an important development due to the scaling reduction afforded by THC and the ubiquity of molecular ions, radicals, and other open-shell reactive species. The complexity of wavefunction-based quantum-chemical methods such as Møller–Plesset and coupled cluster theory is reflected in the steep scaling of the computational costs with the molecular size. The least-squares tensor hypercontraction (LS-THC) method is an efficient, single-step factorization for the two-electron integral tensor but can also be used to factorize the double excitation amplitudes, leading to significant scaling reduction. Herein we extend this promising method to open-shell variants of LS-THC-MP2 and -MP3 by using diagrammatic techniques and explicit spin summation. The accuracy of the resulting methods for open-shell species is benchmarked on standard test systems such as regular alkanes as well as realistic systems involving bond breaking, radical stabilization, and other effects. We find that open-shell LS-THC-MPn methods exhibit errors highly comparable to those produced by closed-shell LS-THC-MPn and are highly insensitive to particular chemical interactions, geometries, or even moderate spin contamination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

The first isomerization reaction of an alkylperoxy (RO 2 ) radical holds significant importance in low-temperature oxidation, as it governs the branching ratios of the hydroperoxyalkyl (QOOH) radicals, which influence the competition between the chain-propagation and chain-branching reactions. In this study, we systematically calculated high-pressure rate rules for the RO 2 isomerization reaction of monoethers, exploring 5-, 6-, 7-, and 8-membered ring transition states. Primary, secondary, and tertiary carbon sites, where both the abstracting peroxy group and the abstracted hydrogen are located, were considered, with particular emphasis on distinguishing between secondary carbons adjacent (alpha) and nonadjacent to the ether functional group. Using the G4//B3LYP/6-311++G(2df,2pd) level of theory and the transition state theory, we estimated the rate constants and the Arrhenius coefficient for over 120 possible isomerization reactions. We examined the effect of ring size and ring atoms, revealing that 6- and 7-membered ring isomerizations were generally the fastest. The impact of the ether functional group on transition states was investigated by comparing reactions with identical ring size, peroxy, and radical positions, but with the ether functional group positioned either outside (i.e., out) or inside (i.e., in) the transition state ring, leading to differences in the rate constants. When comparing to analogous alkane rate constants, differences of up to an order of magnitude were observed, underscoring the need for caution when assigning rate rules by analogy. We applied our rate constants in the di-iso-butyl ether kinetic model and evaluated their influence on low-temperature chemistry finding that they altered the branching ratios by up to a factor of 9, highlighting the significance of site-specific rate constants for more accurate low-temperature modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗