Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

Ring-Opening Terpolymerisation of Elemental Sulfur Waste with Propylene Oxide and Carbon Disulfide via Lithium Catalysis.

Elemental sulfur, a waste product of the oil refinement process, represents a promising raw material for the synthesis of degradable polymers. We show that simple lithium alkoxides facilitate the polymerisation of elemental sulfur S8 with industrially relevant propylene oxide (PO) and CS2 (a base chemical sourced from waste S8 itself) to give poly(monothiocarbonate-alt-Sx) in which x can be controlled by the amount of supplied sulfur. The in situ generation of thiolate intermediates obtained by a rearrangement, which follows CS2 and PO incorporation, allows to combine S8 and epoxides into one polymer sequence that would otherwise not be possible. Mechanistic investigations reveal that alkyl oligosulfide intermediates from S8 ring opening and sulfur chain length equilibration represent the better nucleophiles for inserting the next PO if compared to the trithiocarbonates obtained from the competing CS2 addition, which causes the sequence selectivity. The polymers can be crosslinked in situ with multifunctional thiols to yield reprocessable and degradable networks. Our report demonstrates how mechanistic understanding allows to combine intrinsically incompatible building blocks for sulfur waste utilisation.

Classical and quantum trial wave functions in auxiliary-field quantum Monte Carlo applied to oxygen allotropes and a CuBr2 model system.

In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representative of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and the relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which is industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow it to generate energies close to chemical accuracy compared to full configuration interaction or experimental results. In the field of material science, we study the electronic structure properties of cuprates through the quasi-1D Fermi-Hubbard model derived from CuBr2, where we find that trial wave functions with both significantly larger fidelities and lower energies over a mean-field solution do not necessarily lead to AFQMC results closer to the exact ground state energy.

Reaction of OH with Aliphatic and Aromatic Isocyanates.

Isocyanates are highly relevant industrial intermediates for polyurethane production. In this work, we used quantum chemistry and transition state theory (TST) to investigate the gas-phase reaction of isocyanates with the OH radical, which is likely one of the most significant chemical sinks for these compounds in the troposphere. para-Tolyl-isocyanate (p-tolyl-NCO) was chosen as a proxy substance for the large-volume aromatic diisocyanate species toluene diisocyanate and methylene diphenyl diisocyanate. Besides p-tolyl-NCO + OH, the model reactions CH3NCO + OH, H2C═CHNCO + OH, C6H5-NCO + OH, C6H5-CH3 + OH, and C6H6 + OH have been studied as well to analyze various substituent effects and to allow for comparison with literature. Quantum chemical computations at the CCSD(T)/cc-pV(T,Q → ∞)Z//M06-2X/def2-TZVP level were used as the basis for tunneling-corrected canonical TST calculations. For CH3NCO + OH, H abstraction by OH at the methyl group is the main reaction channel according to our calculations and predicted to be four orders of magnitude faster than OH addition at the NCO group. The calculated rate coefficient (8.8 × 10-14 cm3 molecule-1 s-1) at 298 K is in good agreement with experimental data from the literature. Likewise, for aromatic isocyanates, OH attack at the isocyanate group was found to be only a minor pathway compared to addition to the aromatic ring. In the OH + p-tolyl-NCO reaction, OH addition at the ortho-position relative to the NCO group has been identified as the main initial reaction channel (branching fraction: 53.2%), with smaller but significant branching fractions for the H abstraction at the methyl group (9.6%) as well as the other ring addition reactions (ipso: 2.3%, meta: 24.5%, para: 10.5%, all relative to the NCO group). By comparing OH addition to the aromatic ring in p-tolyl-NCO with the respective ring addition reactions of phenyl isocyanate and toluene, the site-selective reactivity trends observed for ring addition in the OH + p-tolyl-NCO could be rationalized by a dominating positive mesomeric effect of the NCO group and a positive electron-donating (inductive) effect of the CH3 group. Except for the OH ring adduct formed from OH addition in ipso-position to the NCO group, we estimate the first-generation radical intermediates in the OH + p-tolyl-NCO reaction to have sufficiently long lifetimes to react with O2 under atmospheric conditions and undergo typical oxidative reaction cascades like those known for benzene or toluene.

Automated input structure generation for single-ended reaction path optimizations.

The exploration of a reaction network requires highly automated workflows to avoid error-prone and time-consuming manual steps. In this respect, a major bottleneck is the search for transition-state (TS) structures, which frequently fails and, therefore, makes (manual) revision necessary. In this work, we present a technique for obtaining suitable input structures for automated TS searches based on single-ended reaction path optimization algorithms, which makes subsequent TS searches via this method significantly more robust. First, possible input structures are generated based on the spatial alignment of the reactants. The appropriate orientation of reacting groups is achieved via stepwise rotations along selected torsional degrees of freedom. Second, a ranking of the obtained structures is performed according to selected geometric criteria. The main goals are to properly align the reactive atoms, to avoid hindrance within the reaction channel and to resolve steric clashes between the reactants. The developed procedure has been carefully tested on a variety of examples and provides suitable input structures for TS searches within seconds. The method is in daily use in an industrial setting.

Small-Molecule Investigation of Diels-Alder Complexes for Thermoreversible Crosslinking in Polymeric Applications.

Combinations of dienes and dienophiles were examined in order to elicit possible combinations for thermoreversible crosslinking units. Comparison of experimental results and quantum calculations indicated that reaction kinetics and activation energy were much better prediction factors than change in enthalpy for the prediction of successful cycloaddition. Further testing on diene-dienophile pairs that underwent successful cycloaddition determined the feasibility of thermoreversibility/retro-reaction of each of the Diels-Alder compounds. Heating and testing of the compounds in the presence of a trapping agent allowed for experimental determination of reverse kinetics and activation energy for the retro-reaction. The experimental values were in good agreement with quantum calculations. The combination of chemical calculations with experimental results provided a strong insight into the structure-property relationships and how quantum calculations can be used to examine the feasibility of the thermoreversibility of new Diels-Alder complexes in potential polymer systems or to fine-tune thermoreversible Diels-Alder systems already in use.

Role of Acetate Anions in the Catalytic Formation of Isocyanurates from Aromatic Isocyanates.

The formation of isocyanurates via cyclotrimerization of aromatic isocyanates is widely used to enhance the physical properties of a variety of polyurethanes. The most commonly used catalysts in industries are carboxylates for which the exact catalytically active species have remained controversial. We investigated how acetate and other carboxylates react with aromatic isocyanates in a stepwise manner and identified that the carboxylates are only precatalysts in the reaction. The reaction of carboxylates with an excess of aromatic isocyanates leads to irreversible formation of corresponding deprotonated amide species that are strongly nucleophilic and basic. As a result, they are active catalysts during the nucleophilic anionic trimerization, but can also deprotonate urethane and urea species present, which in turn catalyze the isocyanurate formation. The current study also shows how quantum chemical calculations can be used to direct spectroscopic identification of reactive intermediates formed during the active catalytic cycle with predictive accuracy.

Towards a converged strategy for including microsolvation in reaction mechanism calculations.

A major part of chemical conversions is carried out in the fluid phase, where an accurate modeling of the involved reactions requires to also take into account solvation effects. Implicit solvation models often cover these effects with sufficient accuracy but can fail drastically when specific solvent-solute interactions are important. In those cases, microsolvation, i.e., the explicit inclusion of one or more solvent molecules, is a commonly used strategy. Nevertheless, microsolvation also introduces new challenges-a consistent workflow as well as strategies how to systematically improve prediction performance are not evident. For the COSMO and COSMO-RS solvation models, this work proposes a simple protocol to decide if microsolvation is needed and how the corresponding molecular model has to look like. To demonstrate the improved accuracy of the approach, specific application examples are presented and discussed, i.e., the computation of aqueous pKa values and a mechanistic study of the methanol mediated Morita-Baylis-Hillman reaction.

Pentamethylcyclopentadienyl-substituted hypersilylsilylene: reversible and irreversible activation of C[double bond, length as m-dash]C double bonds and dihydrogen.

Recent studies of low-valent main group species underscore their resemblance to transition metal complexes with regards to the ability to activate small molecules. Here, we report synthesis and full characterisation of the persistent (hypersilyl)(pentamethylcyclopentadienyl)silylene Cp*[(Me3Si)3Si]Si: as well as its unique reactivity. Silylene Cp*[(Me3Si)3Si]Si: activates dihydrogen to give the corresponding dihydrosilane Cp*[(Me3Si)3Si]SiH2 at particularly mild conditions as well as ethylene to afford the three-membered cyclic silirane c-Cp*[(Me3Si)3Si]Si(H2CCH2). The addition of N-heterocyclic carbene NHC (NHC = 1,3,4,5-tetramethyl-imidazol-2-ylidene) to dihydrosilane Cp*[(Me3Si)3Si]SiH2 induces the reductive elimination of Cp*H, which according to DFT calculations is thermodynamically preferred over H2 elimination. With NHC, Cp*[(Me3Si)3Si]Si: forms a typical donor-acceptor complex with concomitant change in hapticity of the Cp* ligand from η2 to η1 (σ). In contrast, the reaction with the N-heterocyclic silylene c-[(CH[double bond, length as m-dash]CH(tBuN)2]Si: leads to an unusual "masked" disilene with the former Cp* ligand bridging the two silicon centres. The heterodimer is stable in the solid state, but dissociates reversibly to the constituting silylene fragments in solution.

Mechanistic Investigation into the Acetate-Initiated Catalytic Trimerization of Aliphatic Isocyanates: A Bicyclic Ride.

The acetate-initiated aliphatic isocyanate trimerization to isocyanurate was investigated by state-of-the-art analytical and computational methods. Although the common cyclotrimerization mechanism assumes the consecutive addition of three equivalents of isocyanate to acetate prior to product formation, we found that the underlying mechanism is more complex. In this work, we demonstrate that the product, in fact, is formed via the connection of two unexpected catalytic cycles, with acetate being only the precatalyst. The initial discovery of a precatalyst activation by quantum chemical computations and the resulting first catalysis cycle were corroborated by mass spectrometric and NMR experiments, thereby additionally revealing a catalyst migration to the second catalytic cycle. These results were further confirmed by computations, completing the full mechanistic understanding of this catalytic system. Identification of a side product with undesired properties for final coating applications allows for process optimization in the chemical industry.

Accuracy and Resource Estimations for Quantum Chemistry on a Near-Term Quantum Computer.

One of the most important application areas of molecular quantum chemistry is the study and prediction of chemical reactivity. Large-scale, fully error-tolerant quantum computers could provide exact or near-exact solutions to the underlying electronic structure problem with exponentially less effort than a classical computer thus enabling highly accurate predictions for comparably large molecular systems. In the nearer future, however, only "noisy" devices with a limited number of qubits that are subject to decoherence will be available. For such near-term quantum computers the hybrid quantum-classical variational quantum eigensolver algorithm in combination with the unitary coupled-cluster ansatz (UCCSD-VQE) [ Peruzzo et al. Nat. Commun. 2014 , 5 , 4213 and McClean et al. New J. Phys. 2016 , 18 , 023023 ] has become an intensively discussed approach that could provide accurate results before the dawn of error-tolerant quantum computing. In this work we present an implementation of UCCSD-VQE that allows for the first time to treat both open- and closed-shell molecules. We study the accuracy of the obtained energies for nine small molecular systems as well as for four exemplary chemical reactions by comparing to well-established electronic structure methods like (nonunitary) coupled-cluster and density functional theory. Finally, we roughly estimate the required quantum hardware resources to obtain "useful" results for practical purposes.

Design Principles for Rational Polyurethane Catalyst Development.

Tertiary amine catalysts are essential components in manufacturing polyurethane materials. The low-emission requirements for indoor applications are typically achieved by employing tertiary amines with catalytically active N, N-dimethyl groups as the base catalyst and a longer alkyl substituent with a reactive end, that is, alcohol or amine, to incorporate it in the polyurethane matrix. N, N-dimethyl groups are, however, oxidized when exposed to air and lead to undesired formaldehyde emissions. Here, we employ modern quantum chemical methods to understand design principles how the structure of tertiary amine catalysts having N, N-dimethyl groups can be modified to avoid this source of formaldehyde formation but still preserve their catalytic activity. We found the pyrrolidine derivative of commonly used N, N-dimethylated catalysts to be the most promising candidate and developed design principles to rationalize why longer alkyl chains or larger ring sizes inhibit the catalytic activity. The computationally predicted catalyst performances were confirmed experimentally in model polyurethane systems for selected amine catalysts, and emission measurements showed that the formaldehyde emission was completely suppressed when pyrrolidine derivative was used as a catalyst. Our results further illustrate how condensed phase reactions can be predicted using quantum chemical methods and that to account for steric hindrance near the reaction center, it was also necessary to include conformational energy contributions in the calculated activation free energies.

Exploring Solvation Effects in Ligand-Exchange Reactions via Static and Dynamic Methods.

We investigate ligand-exchange reactions of a biomimetic Co(II)-based heterocubane complex in aqueous solution by means of various approaches for consideration of solvent effects. Static calculations based on geometry optimizations carried out in vacuum, with solvent continuum models, or with several explicit solvent molecules have been carried out as well as density functional theory (DFT)-based molecular dynamics simulations. In addition, reaction pathways and barriers have been elucidated via nudged elastic band calculations and metadynamics. The results show that static approaches with approximate consideration of the solvent environment lead to reaction energies, which may change drastically depending on the method employed. A more sophisticated approach is DFT-molecular dynamics at ambient conditions with full solvation, i.e. enough solvent molecules to retain bulk water properties far from the solute, which, however, comes with a much higher computational cost. The investigated example of the exchange of an acetate ligand by a hydroxide demonstrates that entropic contributions can be vital and consideration of electronic energies alone may be a rather rough approximation.

Influence of Antimony-Halogen Additives on Flame Propagation.

A kinetic model for flame inhibition by antimony-halogen compounds in hydrocarbon flames is developed. Thermodynamic data for the relevant species are assembled from the literature, and calculations are performed for a large set of additional species of Sb-Br-C-H-O system. The main Sb- and Br-containing species in the combustion products and reaction zone are determined using flame equilibrium calculations with a set of possible Sb-Br-C-H-O species, and these are used to develop the species and reactions in a detailed kinetic model for antimony flame inhibition. The complete thermodynamic data set and kinetic mechanism are presented. Laminar burning velocity simulations are used to validate the mechanism against available data in the literature, as well as to explore the relative performance of the antimony-halogen compounds. Further analysis of the premixed flame simulations has unraveled the catalytic radical recombination cycle of antimony. It includes (primarily) the species Sb, SbO, SbO2, and HOSbO, and the reactions: Sb+O+M=SbO+M; Sb+O2+M=SbO2+M; SbO+H=Sb+OH; SbO+O=Sb+O2; SbO+OH+M=HOSbO+M; SbO2+H2O=HOSbO+OH; HOSbO+H=SbO+H2O; SbO+O+M=SbO2+M. The inhibition cycles of antimony are shown to be more effective than those of bromine, and intermediate between the highly effective agents CF3Br and trimethylphosphate. Preliminary examination of a Sb/Br gas-phase system did not show synergism in the gas-phase catalytic cycles (i.e., they acted essentially independently).

Haloperoxidase Mimicry by CeO2-x Nanorods Combats Biofouling.

CeO2-x nanorods are functional mimics of natural haloperoxidases. They catalyze the oxidative bromination of phenol red to bromophenol blue and of natural signaling molecules involved in bacterial quorum sensing. Laboratory and field tests with paint formulations containing 2 wt% of CeO2-x nanorods show a reduction in biofouling comparable to Cu2 O, the most typical biocidal pigment.

CMIRS Solvation Model for Methanol: Parametrization, Testing, and Comparison with SMD, SM8, and COSMO-RS.

The new continuum solvation model, composite method for implicit representation of solvent (CMIRS), proposed by Pomogaeva and Chipman and implemented in GAMESS was parametrized for methanol solvent, with the aim of using it for ionic reactions in solution. The model was tested for predicting single-ion solvation free energy, pKa of acids and protonated bases, and the activation free-energy barriers of SN2 and SNAr reactions in methanol. A comparison was performed with other continuum models, such as SMD, SM8, and COSMO-RS. For a prediction of pKa and free-energy barriers, the order of performance was CMIRS > COSMO-RS > SMD > SM8. In particular, the CMIRS model is much superior to the other continuum models for predicting pKa of acids (without empirical corrections) and is able to evenly describe hard ions like methoxide and charge-dispersed ions like 2,4,6-trinitrophenol. On the basis of our results, we suggest that the field-extremum contribution, present in CMIRS, should be included in continuum solvation models, which can result in substantial improvement in the modeling of ionic reactions in solution.

Multiscale Materials Modeling in an Industrial Environment.

In this review, we sketch the materials modeling process in industry. We show that predictive and fast modeling is a prerequisite for successful participation in research and development processes in the chemical industry. Stable and highly automated workflows suitable for handling complex systems are a must. In particular, we review approaches to build and parameterize soft matter systems. By satisfying these prerequisites, efficiency for the development of new materials can be significantly improved, as exemplified here for formulation polymer development. This is in fact in line with recent Materials Genome Initiative efforts sponsored by the US government. Valuable contributions to product development are possible today by combining existing modeling techniques in an intelligent fashion, provided modeling and experiment work hand in hand.

Mechanistic Investigations of the Stereoselective Rare Earth Metal-Mediated Ring-Opening Polymerization of ß-Butyrolactone.

Poly(3-hydroxybutyrate) (PHB) is produced by numerous bacteria as carbon and energy reserve storage material. Whereas nature only produces PHB in its strictly isotactic (R) form, homogeneous catalysis, when starting from racemic (rac) ß-butyrolactone (BL) as monomer, can in fact produce a wide variety of tacticities. The variation of the metal center and the surrounding ligand structure enable activity as well as tacticity tuning. However, no homogeneous catalyst exists to date that is easy to modify, highly active, and able to produce PHB with high isotacticities from rac-ß-BL. Therefore, in this work, the reaction kinetics of various 2-methoxyethylamino-bis(phenolate) lanthanide (Ln=Sm, Tb, Y, Lu) catalysts are examined in detail. The order in monomer and catalyst are determined to elucidate the reaction mechanism and the results are correlated with DFT calculations of the catalytic cycle. Furthermore, the enthalpies and entropies of the rate-determining steps are determined through temperature-dependent in situ IR measurements. Experimental and computational results converge in one specific mechanism for the ring-opening polymerization of BL and even allow us to rationalize the preference for syndiotactic PHB.

Mechanistic Aspects of a Highly Active Dinuclear Zinc Catalyst for the Co-polymerization of Epoxides and CO2.

The dinuclear zinc complex reported by us is to date the most active zinc catalyst for the co-polymerization of cyclohexene oxide (CHO) and carbon dioxide. However, co-polymerization experiments with propylene oxide (PO) and CO2 revealed surprisingly low conversions. Within this work, we focused on clarification of this behavior through experimental results and quantum chemical studies. The combination of both results indicated the formation of an energetically highly stable intermediate in the presence of propylene oxide and carbon dioxide. A similar species in the case of cyclohexene oxide/CO2 co-polymerization was not stable enough to deactivate the catalyst due to steric repulsion.