The Study of Reactive Intermediates in Condensed Phases

Predicting Selectivity with a Bifurcating Surface: Inaccurate Model or Inaccurate Statistics of Dynamics?

Reactions on post-transition-state bifurcation (PTSB) energy surfaces are an important class of reaction in which classical rate theories, such as the transition state theory, fail to account for the selectivity. Quasiclassical trajectory molecular dynamic (QCT-MD) simulation is an important computational approach to understanding reactions mechanisms, especially for reactions that cannot be predicted from conventional rate theories. However, the applicability of direct dynamic simulations is hampered by huge computational costs for collecting a statistically meaningful set of trajectories, making it difficult to compare simulation results with theoretical or physical insights-based predictions (non-MD predictions). In this work, we examine the PTSB of Schmidt-Aubé reactions studied by Tantillo and co-workers. With machine-learning using kernel-ridge regression (KRR) to predict atomic forces, statistical reliability was enhanced by significantly increasing the number of trajectories. With KRR, the bottleneck of simulating dynamics (atomic forces in QCT-MD with density functional theory) was accelerated more than 100-fold. We found that this KRR-aided QCT-MD approach is successful in predicting branching ratios with a much larger number of trajectories. With our approach, statistical errors are greatly reduced, and hypothetical non-MD models for predicting selectivity are tested with much higher confidence. By comparison with non-MD models, dynamical properties that affect branching ratios become more clearly described.

A 21st Century View of Allowed and Forbidden Electrocyclic Reactions

In 1965, Woodward and Hoffmann proposed a theory to predict the stereochemistry of electrocyclic reactions, which, after expansion and generalization, became known as the Woodward-Hoffmann Rules. Subsequently, Longuet-Higgins and Abrahamson used correlation diagrams to propose that the stereoselectivity of electrocyclizations could be explained by the correlation of reactant and product orbitals with the same symmetry. Immediately thereafter, Hoffmann and Woodward applied correlation diagrams to explain the mechanism of cycloadditions. We describe these discoveries and their evolution. We now report an investigation of various electrocyclic reactions using DFT and CASSCF. We track the frontier molecular orbitals along the intrinsic reaction coordinate and modeled trajectories and examine the correlation between HOMO and LUMO for thermally forbidden systems. We also investigate the electrocyclizations of several highly polarized systems for which the Houk group had predicted that donor-acceptor substitution can induce zwitterionic character, thereby providing low-energy pathways for formally forbidden reactions. We conclude with perspectives on the field of pericyclic reactions, including a refinement as the meaning of Woodward and Hoffmann's "Violations. There are none!" Lastly, we comment on the burgeoning influence of computations on all fields of chemistry.

Understanding chemistry: from "heuristic (soft) explanations and reasoning by analogy" to "quantum chemistry"

"Soft theories," i.e., "heuristic models based on reasoning by analogy" largely drove chemistry understanding for 150 years or more. But soft theories have their limitations and with the expansion of chemistry in the mid-20th century, more and more inexplicable (by soft theory) experimental results were being obtained. In the past 50 years, quantum chemistry, most often in the guise of applied theoretical chemistry including computational chemistry, has provided (a) the underlying "hard evidence" for many soft theories and (b) the explanations for chemical phenomena that were unavailable by soft theories. In this publication, we define "hard theories" as "theories derived from quantum chemistry." Both soft and hard theories can be qualitative and quantitative, and the "Houk quadrant" is proposed as a helpful categorization tool. Furthermore, the language of soft theories is often used appropriately to describe quantum chemical results. A valid and useful way of doing science is the appropriate use and application of both soft and hard theories along with the best nomenclature available for successful communication of results and ideas.

Investigation of organic monoradicals reactivity using surface-enhanced Raman spectroscopy

Surface-enhanced Raman spectroscopy (SERS) and self-assembled monolayer (SAM) approaches were used to investigate the reactions of organic monoradicals with methanol. An attempt was made to generate monoradicals from thiophenols and phenylmethanethiols substituted with bromine, iodine, and nitro groups by irradiation with UV light. Monolayers of radical precursors were deposited on SERS substrates, which were then immersed in methanol and irradiated for 1 and/or 3, 6, 12 and 24 h in a UV photochemical reactor. Pre- and postreaction SERS spectra were obtained by using a confocal Raman microscope and compared with the spectra of expected products of the radical reaction with methanol. Our studies have shown that the efficiency of monoradical generation is highly dependent on the chemical structure of the precursor. In addition, it is shown that both the SERS substrate and experimental conditions used strongly influence the obtained results.

Portable Models for Entropy Effects on Kinetic Selectivity

Differences in entropies of competing transition states can direct kinetic selectivity. Understanding and modeling such entropy differences at the molecular level is complicated by the fact that entropy is statistical in nature; i.e., it depends on multiple vibrational states of transition structures, the existence of multiple dynamically accessible pathways past these transition structures, and contributions from multiple transition structures differing in conformation/configuration. The difficulties associated with modeling each of these contributors are discussed here, along with possible solutions, all with an eye toward the development of portable qualitative models of use to experimentalists aiming to design reactions that make use of entropy to control kinetic selectivity.

New Insights and Predictions into Complex Homogeneous Reactions Enabled by Computational Chemistry in Synergy with Experiments: Isotopes and Mechanisms

Homogeneous catalysis and biocatalysis have been widely applied in synthetic, medicinal, and energy chemistry as well as synthetic biology. Driven by developments of new computational chemistry methods and better computer hardware, computational chemistry has become an essentially indispensable mechanistic "instrument" to help understand structures and decipher reaction mechanisms in catalysis. In addition, synergy between computational and experimental chemistry deepens our mechanistic understanding, which further promotes the rational design of new catalysts. In this Account, we summarize new or deeper mechanistic insights (including isotope, dispersion, and dynamical effects) into several complex homogeneous reactions from our systematic computational studies along with subsequent experimental studies by different groups. Apart from uncovering new mechanisms in some reactions, a few computational predictions (such as excited-state heavy-atom tunneling, steric-controlled enantioswitching, and a new geminal addition mechanism) based on our mechanistic insights were further verified by ensuing experiments.The Zimmerman group developed a photoinduced triplet di-π-methane rearrangement to form cyclopropane derivatives. Recently, our computational study predicted the first excited-state heavy-atom (carbon) quantum tunneling in one triplet di-π-methane rearrangement, in which the reaction rates and ¹²C/¹³C kinetic isotope effects (KIEs) can be enhanced by quantum tunneling at low temperatures. This unprecedented excited-state heavy-atom tunneling in a photoinduced reaction has recently been verified by an experimental ¹²C/¹³C KIE study by the Singleton group. Such combined computational and experimental studies should open up opportunities to discover more rare excited-state heavy-atom tunneling in other photoinduced reactions. In addition, we found unexpectedly large secondary KIE values in the five-coordinate Fe(III)-catalyzed hetero-Diels-Alder pathway, even with substantial C-C bond formation, due to the non-negligible equilibrium isotope effect (EIE) derived from altered metal coordination. Therefore, these KIE values cannot reliably reflect transition-state structures for the five-coordinate metal pathway. Furthermore, our density functional theory (DFT) quasi-classical molecular dynamics (MD) simulations demonstrated that the coordination mode and/or spin state of the iron metal as well as an electric field can affect the dynamics of this reaction (e.g., the dynamically stepwise process, the entrance/exit reaction channels).Moreover, we unveiled a new reaction mechanism to account for the uncommon Ru(II)-catalyzed geminal-addition semihydrogenation and hydroboration of silyl alkynes. Our proposed key gem-Ru(II)-carbene intermediates derived from double migrations on the same alkyne carbon were verified by crossover experiments. Additionally, our DFT MD simulations suggested that the first hydrogen migration transition-state structures may directly and quickly form the key gem-Ru-carbene structures, thereby "bypassing" the second migration step. Furthermore, our extensive study revealed the origin of the enantioselectivity of the Cu(I)-catalyzed 1,3-dipolar cycloaddition of azomethine ylides with β-substituted alkenyl bicyclic heteroarenes enabled by dual coordination of both substrates. Such mechanistic insights promoted our computational predictions of the enantioselectivity reversal for the corresponding monocyclic heteroarene substrates and the regiospecific addition to the less reactive internal C═C bond of one diene substrate. These predictions were proven by our experimental collaborators. Finally, our mechanistic insights into a few other reactions are also presented. Overall, we hope that these interactive computational and experimental studies enrich our mechanistic understanding and aid in reaction development.

Heavy-Atom Kinetic Isotope Effects: Primary Interest or Zero Point?

Chemists have many options for elucidating reaction mechanisms. Global kinetic analysis and classic transition-state probes (e.g., LFERs, Eyring) inevitably form the cornerstone of any strategy, yet their application to increasingly sophisticated synthetic methodologies often leads to a wide range of indistinguishable mechanistic proposals. Computational chemistry provides powerful tools for narrowing the field in such cases, yet wholly simulated mechanisms must be interpreted with great caution. Heavy-atom kinetic isotope effects (KIEs) offer an exquisite but underutilized method for reconciling the two approaches, anchoring the theoretician in the world of calculable observables and providing the experimentalist with atomistic insights. This Perspective provides a personal outlook on this synergy. It surveys the computation of heavy-atom KIEs and their measurement by NMR spectroscopy, discusses recent case studies, highlights the intellectual reward that lies in alignment of experiment and theory, and reflects on the changes required in chemical education in the area.

Effect of solvent motions on the dynamics of the Diels-Alder reaction

How solvent motions affect the dynamics of chemical reactions in which the solute undergoes a substantial shape change is a fundamental but elusive issue. This work utilizes reactive simulation and Grote-Hynes theory to explore the effect of solvent motions on the dynamics of the Diels-Alder reaction (in the reverse direction, this reaction involves very substantial solute expansion) in aprotic solvents. The results reveal that the solvent environment is not sufficiently constraining to influence transition state passage dynamics, with the calculated transmission coefficients being close to unity. Even when solvent motions are suppressed or artificially slowed down, the solvent only affects the reaction dynamics in the transition state region to a very small extent. The only notable effect of solvent occurs far from the transition state region and corresponds to caging of the reactants within the reactant well.

Solvent Effects on Ultrafast Photochemical Pathways

Photochemical reactions are increasingly being used for chemical and materials synthesis, for example, in photoredox catalysis, and generally involve photoexcitation of molecular chromophores dissolved in a liquid solvent. The choice of solvent influences the outcomes of the photochemistry because solute-solvent interactions modify the energies of and crossings between electronic states of the chromophores, and they affect the evolving structures of the photoexcited molecules. Ultrafast laser spectroscopy methods with femtosecond to picosecond time resolution can resolve the dynamics of these photoexcited molecules as they undergo structural and electronic changes, relax back to the ground state, dissipate their excess internal energy to the surrounding solvent, or undergo photochemical reactions. In this Account, we illustrate how experimental studies using ultrafast lasers can reveal the influences that different solvents or cosolutes exert on the photoinduced nonadiabatic dynamics of internal conversion and intersystem crossing in nonradiative relaxation pathways. Although the environment surrounding a solute molecule is rapidly changing, with fluctuations in the coordination to neighboring solvent molecules occurring on femtosecond or picosecond time scales, we show that it is possible to photoexcite selectively only those molecular chromophores transiently experiencing specific solute-solvent interactions such as intermolecular hydrogen bonding.The effects of different solvation environments on the photodynamics are illustrated using four selected examples of photochemical processes in which the solvent has a marked effect on the outcomes. We first consider two aromatic carbonyl compounds, benzophenone and acetophenone, which are known to undergo fast intersystem crossing to populate the first excited triplet state on time scales of a few picoseconds. We show that the nonadiabatic excited-state dynamics are modified by transient hydrogen bonding of the carbonyl group to a protic solvent or by coordination to a metal cation cosolute. We then examine how different solvents modify the competition between two alternative relaxation pathways in a photoexcited UVA-sunscreen molecule, diethylamino hydroxybenzoyl hexyl benzoate (DHHB). This relaxation back to the ground electronic state is an essential part of the effective operation of the sunscreen compound, but the dynamics are sensitive to the surrounding environment. Finally, we consider how solvents of different polarity affect the energies and lifetimes of excited states with locally excited or charge-transfer character in heterocyclic organic compounds used as excited-state electron donors for photoredox catalysis. With these and other examples, we seek to develop a molecular level understanding of how the choice of solution environment might be used to control the outcomes of photochemical reactions.

Melding of Experiment and Theory Illuminates Mechanisms of Metal-Catalyzed Rearrangements: Computational Approaches and Caveats

This review summarizes approaches and caveats in computational modeling of transition-metal-catalyzed sigmatropic rearrangements involving carbene transfer. We highlight contemporary examples of combined synthetic and theoretical investigations that showcase the synergy achievable by integrating experiment and theory.

1 Introduction

2 Mechanistic Models

3 Theoretical Approaches and Caveats

3.1 Recommended Computational Tools

3.2 Choice of Functional and Basis Set

3.3 Conformations and Ligand-Binding Modes

3.4 Solvation

4 Synergy of Experiment and Theory – Case Studies

4.1 Metal-Bound or Free Ylides?

4.2 Conformations and Ligand-Binding Modes of Paddlewheel Complexes

4.3 No Metal, Just Light

4.4 How To ‘Cope’ with Nonstatistical Dynamic Effects

5 Outlook

Molecules, the Ultimate Nanomotor: Linking Chemical Reaction Intermediates to their Molecular Diffusivity

The intellectual community focused on nanomotors has recently become interested in extending these concepts to individual molecules. Here, we study a chemical reaction according to whose mechanism some intermediate species should speed up while others slow down in predictable ways, if the nanomotor hypothesis of boosted diffusion holds. Accordingly, we scrutinize the absolute diffusion coefficient (D) during intermediate steps of the catalytic cycle for the CuAAC reaction (copper-catalyzed azide-alkyne cycloaddition click reaction), using proton pulsed field-gradient nuclear magnetic resonance to discriminate between the diffusion of various reaction intermediates. We observe time-dependent diffusion that is enhanced for some intermediate molecular species and depressed for those whose size increases owing to complex formation. These findings point to the failure of the conventional Stokes-Einstein equation to fully explain diffusivity during chemical reaction. Without attempting a firm explanation, this paper highlights aspects of the physics of chemical reactions that are imperfectly understood and presents systematic data that can be used to assess hypotheses.

On the origins of the mechanistic variants in the thermal reactions of Sx+ (x = 1-3) with benzene

The S-π interaction between sulfur atom(s) and aromatic ring prevails in chemical and biochemical processes. The thermal gas-phase reactions of the S_n⁺ (n = 1-3) ions with benzene have been explored by using Quadrupole-Ion Trap (Q-IT) mass spectrometry complemented by quantum chemical calculations. Charge transfer was found to be the only reaction channel for S₂⁺/C₆H₆, while both charge transfer and bond activation are available for the S⁺/C₆H₆ and S₃⁺/C₆H₆ couples. Upon interrogating the associated electronic origins, multiple factors were found to matter for these processes. In contrast to the σ-type two-center three-electron (2c-3e) S-π hemibond as reported previously, unusual S-π hemibonds were addressed for the S_n⁺/C₆H₆ couples, i.e. the 2c-3e π(S061Eπ) and the three-center three-electron (3c-3e) σ(S₂061Eπ) hemibonds. Such S-π interaction was found to be responsible for the charge transfer processes in S⁺/C₆H₆ and S₂⁺/C₆H₆, but uninvolved in any transformation for S₃⁺/C₆H₆.

Machine learning classification of disrotatory IRC and conrotatory non-IRC trajectory motion for cyclopropyl radical ring opening

Quasiclassical trajectory analysis is now a standard tool to analyze non-minimum energy pathway motion of organic reactions. However, due to the large amount of information associated with trajectories, quantitative analysis of the dynamic origin of reaction selectivity is complex. For the electrocyclic ring opening of cyclopropyl radical, more than 4000 trajectories were run showing that allyl radicals are formed through a mixture of disrotatory intrinsic reaction coordinate (IRC) motion as well as conrotatory non-IRC motion. Geometric, vibrational mode, and atomic velocity transition-state features from these trajectories were used for supervised machine learning analysis with classification algorithms. Accuracy >80% with a random forest model enabled quantitative and qualitative assessment of transition-state trajectory features controlling disrotatory IRC versus conrotatory non-IRC motion. This analysis revealed that there are two key vibrational modes where their directional combination provides prediction of IRC versus non-IRC motion.

Direct observation of coherent femtosecond solvent reorganization coupled to intramolecular electron transfer

It is well known that the solvent plays a critical role in ultrafast electron-transfer reactions. However, solvent reorganization occurs on multiple length scales, and selectively measuring short-range solute-solvent interactions at the atomic level with femtosecond time resolution remains a challenge. Here we report femtosecond X-ray scattering and emission measurements following photoinduced charge-transfer excitation in a mixed-valence bimetallic (FeⁱⁱRuⁱⁱⁱ) complex in water, and their interpretation using non-equilibrium molecular dynamics simulations. Combined experimental and computational analysis reveals that the charge-transfer excited state has a lifetime of 62 fs and that coherent translational motions of the first solvation shell are coupled to the back electron transfer. Our molecular dynamics simulations identify that the observed coherent translational motions arise from hydrogen bonding changes between the solute and nearby water molecules upon photoexcitation, and have an amplitude of tenths of ångströms, 120-200 cm^-1 frequency and ~100 fs relaxation time. This study provides an atomistic view of coherent solvent reorganization mediating ultrafast intramolecular electron transfer.

Metal-Micelle Cooperativity: Phosphine Ligand-Free Ultrasmall Palladium(II) Nanoparticles for Oxidative Mizoroki-Heck-type Couplings in Water at Room Temperature

The amphiphile PS-750-M generates stable, phosphine ligand-free, and catalytically active ultrasmall Pd(II) nanoparticles (NPs) from Pd(OAc)₂, preventing their precipitation, polymerization, and oxidation state changes. PS-750-M directly interacts with Pd(II) NP surfaces, as confirmed by high-resolution mass spectrometry and IR spectroscopy, resulting in their high stability. The Pd cations in NPs are most likely held together by hydroxides and acetate ions. The NPs were characterized by HRTEM, revealing their morphology and particle size distribution, and by HRMS and IR, providing evidence for NP-amphiphile interaction. The NP catalytic activity was examined in the context of oxidative Mizoroki-Heck-type couplings in water at room temperature. Hot filtration, hot extraction, and three-phase tests indicate heterogeneous catalysis occurring at the micellar interface rather than homogeneous catalysis occurring in the solution. NMR studies indicate that the catalytic activity stems from metal cation-π interactions of the styrene along with transmetalation by the arylboronic acid, followed by insertion and β-H elimination to furnish the coupled product along with the reoxidation of Pd by benzoquinone to complete the catalytic cycle. This method is very mild and sustainable, both in terms of NP synthesis and subsequent catalysis, and shows broad substrate scope while circumventing the need for organic solvents for this important class of couplings.

Vibrational Energy Relaxation of Deuterium Fluoride in d-Dichloromethane: Insights from Different Potentials

Vibrationally excited deuterium fluoride (DF) formed by fluorine atom reaction with a solvent was found (Science, 2015, 347, 530) to relax rapidly (less than 10 ps) in acetonitrile-d₃ (CD₃CN) and dichloromethane-d₂ (CD₂Cl₂). However, insights into how CD₂Cl₂ facilitates this energy relaxation have so far been lacking, given the weak interaction between DF and a single CD₂Cl₂. In this work, we report the results of reactive simulations with a two-state reactive empirical valence bond (EVB) potential to study the energy deposited into nascent DF after transition-state passage and of nonequilibrium molecular dynamics simulations using multiple different potential energy functions to model the relaxation dynamics. For these second simulations, we used the standard Merck molecular force field (MMFF) potential, an MMFF-based covalent-ionic empirical valence bond (EVB) potential (EVB^CI), a newly developed potential [referred to as MMFF(rDF)] which extends upon the MMFF potential by making the DF/CD₂Cl₂ interaction depend on the value of the D-F bond stretching coordinate and by taking the anisotropic charge distribution of the solvent molecules into account, the polarizable atomic multipole optimized energetics for biomolecular applications (AMOEBA) potential, and the quantum mechanics/molecular mechanics (QM/MM) potential. The relaxation is revealed to be highly sensitive to the potential used. Neither standard MMFF nor EVB^CI reproduces the experimentally observed rapid relaxation dynamics, and they also fail to provide a good description of the interaction potential between DF and CD₂Cl₂ as calculated using CCSD(T)-F12. This is attributed to the use of a point-charge model for the solute and to failing to model the anisotropic electrostatic properties of CD₂Cl₂. The MMFF(rDF), AMOEBA, and QM/MM potentials all reproduce the CCSD(T)-F12 two-body DF---CD₂Cl₂ interaction potential rather well but only with the QM/MM approach is fast vibrational relaxation obtained (lifetimes of ∼288, ∼186, and ∼8 ps, respectively), which we attribute to differences in the solute-solvent local structure. With QM/MM, a unique "many-body" interaction pattern in which DF is in close contact with two solvent Cl atoms and more than three solvent D atoms is found, but this structure is not seen with other potentials. The QM/MM dynamics also display enhanced solute-solvent interactions with vibrationally excited DF that induce a DF band redshift and hence a resonant overlap with solvent C-D modes, which facilitate the intermolecular energy transfer. Our work also suggests that potentials used to model energy relaxation need to capture the fine structure of solute-solvent interactions and not just the two-body part.

Synthesis and characterization of supported stabilized palladium nanoparticles for selective hydrogenation in water at low temperature

Zirconia supported vacant phosphotungstate stabilized Pd nanoparticles (Pd-PW11/ZrO2) were prepared using a simple impregnation and post reduction method, characterized and their efficiency for selective C[double bond, length as m-dash]C hydrogenation of unsaturated compounds explored. The establishment of a hydrogenation strategy at low temperature using water as solvent under mild conditions makes the present system environmentally benign and green. The catalyst shows outstanding activity (96% conversion) with just a small amount of Pd(0) (0.0034 mol%) with high substrate/catalyst ratio (29 177/1), TON (28 010) and TOF (14 005 h-1) for cyclohexene (as a model substrate) hydrogenation. The catalyst was recovered by simple centrifugation and reused for up to five catalytic cycles without alteration in its activity. The present catalyst was found to be viable towards different substrates with excellent activity and TON (18 000 to 28 800). A study on the effect of addenda atom shows that the efficiency of the catalyst can be enhanced greatly by increasing the number of counter protons. This challenging strategy would greatly benefit sustainable development in chemistry as it diminishes the use of organic solvents and offers economic and environmental benefits as water is cheap and non-toxic.

Real-Time Reaction Monitoring with In Operando Flow NMR and FTIR Spectroscopy: Reaction Mechanism of Benzoxazole Synthesis

In operando observation of reaction intermediates is crucial for unraveling reaction mechanisms. To address the sensitivity limitations of commercial ReactIR, a flow cell was integrated with a Fourier transform infrared (FTIR) spectrometer yielding a "flow FTIR" device coupled with an NMR spectrometer for the elucidation of reaction mechanisms. The former device detects the low-intensity IR peaks of reaction intermediates by adjusting the path length of the FTIR sample cell, whereas the flow NMR allows the quantitative analysis of reaction species, thus offsetting the limitations of IR spectroscopy resulting from different absorption coefficients of the normal modes. Using the flow NMR and FTIR device, the controversial mechanism of benzoxazole synthesis was conclusively determined by spectroscopic evaluation of the reaction intermediates. This system enabled the accurate acquisition of previously elusive kinetic data, such as the reaction time and rate-determining step. The implementation of reaction flow cells into NMR and FTIR systems could be widely applied to study various reaction mechanisms, including dangerous and harsh reactions, thus avoiding contact with potentially harmful reaction intermediates.

Ab Initio Molecular Dynamics Simulations of the SN1/SN2 Mechanistic Continuum in Glycosylation Reactions

We report a computational approach to evaluate the reaction mechanisms of glycosylation using ab initio molecular dynamics (AIMD) simulations in explicit solvent. The reaction pathways are simulated via free energy calculations based on metadynamics and trajectory simulations using Born-Oppenheimer molecular dynamics. We applied this approach to investigate the mechanisms of the glycosylation of glucosyl α-trichloroacetimidate with three acceptors (EtOH, i-PrOH, and t-BuOH) in three solvents (ACN, DCM, and MTBE). The reactants and the solvents are treated explicitly using density functional theory. We show that the profile of the free energy surface, the synchronicity of the transition state structure, and the time gap between leaving group dissociation and nucleophile association can be used as three complementary indicators to describe the glycosylation mechanism within the SN1/SN2 continuum for a given reaction. This approach provides a reliable means to rationalize and predict reaction mechanisms and to estimate lifetimes of oxocarbenium intermediates and their dependence on the glycosyl donor, acceptor, and solvent environment.

Accessing the Activation Mechanisms of Ethylene Photo-Polymerization under Pressure by Transient Infrared Absorption Spectroscopy

The ambient temperature photoinduced polymerization of compressed (P < 1 GPa) fluid ethylene was characterized by transient infrared absorption spectroscopy with a resolution of few nanoseconds, 3 orders of magnitude higher than previously reported. The reaction has been studied under both one- and two-photon excitation evidencing in the latter case its occurrence only in the presence of different transition metal oxides. Their photocatalytic activity is ascribed to the stabilization of the excited biradicals through electron density exchange between the d orbitals of the metal and the π antibonding orbitals of ethylene which lengthens the lifetime of the biradicals. In both one- and two-photon activation cases the polymerization is characterized by an initial step distinguished by a molecularity of 0.15 ± 0.02 identified as the activation step of the reaction lasting, in the one-photon excitation case, a few hundreds of nanoseconds. Using pulsed excitation the reaction evolves toward a free radical polymerization only under one-photon excitation whereas the critical concentration of radicals required to propagate the reaction is never achieved in the two-photon excitation case. Comparison with continuous wave excitation unambiguously identifies in the average power released to the sample the key factor to drive quantitatively and qualitatively the polymerization.

Role of Polar Protic Solvents in the Dissociation and Reactivity of Photogenerated Radical Ion Pairs

The UV photolysis of bimolecular charge transfer complexes is employed to yield reactive radical ions in their solvent-equilibrated electronic ground state. In polar protic media, noncovalent complexes of 1,2,4,5-tetracyanobenzene and toluene undergo efficient, ultrafast dissociation to ion pairs and equilibrate with their solvent environment before the resulting radical ions engage in electron transfer and proton abstraction on subnanosecond time scales. Solvent molecules play a critical role in these reactive pathways and in the dissociation and relaxation processes that precede them. We report a clear separation of time scales for these relaxation and reactive processes, which implies that solvent-solute interactions can be used as a tool for tuning the reaction pathways of equilibrated radical ions in solution.

Unusual KIE and dynamics effects in the Fe-catalyzed hetero-Diels-Alder reaction of unactivated aldehydes and dienes

Hetero-Diels-Alder (HDA) reaction is an important synthetic method for many natural products. An iron(III) catalyst was developed to catalyze the challenging HDA reaction of unactivated aldehydes and dienes with high selectivity. Here we report extensive density-functional theory (DFT) calculations and molecular dynamics simulations that show effects of iron (including its coordinate mode and/or spin state) on the dynamics of this reaction: considerably enhancing dynamically stepwise process, broadening entrance channel and narrowing exit channel from concerted asynchronous transition states. Also, our combined computational and experimental secondary KIE studies reveal unexpectedly large KIE values for the five-coordinate pathway even with considerable C-C bond forming, due to equilibrium isotope effect from the change in the metal coordination. Moreover, steric and electronic effects are computationally shown to dictate the C=O chemoselectivity for an α,β-unsaturated aldehyde, which is verified experimentally. Our mechanistic study may help design homogeneous, heterogeneous and biological catalysts for this challenging reaction.

Stable chaos and delayed onset of statisticality in unimolecular dissociation reactions

Statistical models provide a powerful and useful class of approximations for calculating reaction rates by bypassing the need for detailed, and often difficult, dynamical considerations. Such approaches invariably invoke specific assumptions about the extent of intramolecular vibrational energy flow in the system. However, the nature of the transition to the statistical regime as a function of the molecular parameters is far from being completely understood. Here, we use tools from nonlinear dynamics to study the transition to statisticality in a model unimolecular reaction by explicitly visualizing the high dimensional classical phase space. We identify generic features in the phase space involving the intersection of two or more independent anharmonic resonances and show that the presence of correlated, but chaotic, intramolecular dynamics near such junctions leads to nonstatisticality. Interestingly, akin to the stability of asteroids in the Solar System, molecules can stay protected from dissociation at the junctions for several picoseconds due to the phenomenon of stable chaos.

Intramolecular vibrational energy redistribution and the quantum ergodicity transition: a phase space perspective

The onset of facile intramolecular vibrational energy flow can be related to features in the connected network of anharmonic resonances in the classical phase space.

Elucidating cation effects in homogeneously catalyzed formic acid dehydrogenation

In this work, we use density functional theory based molecular dynamics with an explicit description of methanol solvent to study the effect of cations on formic acid dehydrogenation catalyzed by a ruthenium PNP pincer complex (RuPNP). Formic acid dehydrogenation is a two step process that involves the reorientation of the formate moiety bound via its oxygen to a H bound intermediate, followed by the hydride transfer step to form CO2 and the hydrogenated catalyst. We find the reorientation step to proceed with a low barrier in methanol solvent and in the presence of a Li+ cation, while the hydride transfer is significantly hindered by the presence of cations (Li+ and K+). The cation seems to strongly stabilize the negatively charged formate moiety, hindering complete hydride transfer and resulting in a high barrier for this step. This study is a first step towards addressing the exact role of cations in formic acid dehydrogenation reactions.

Ultrafast trans → cis Photoisomerization Dynamics of Alkyl-Substituted Stilbenes in a Supramolecular Capsule

Ultrafast spectroscopy reveals the effects of confinement on the excited-state photoisomerization dynamics for a series of alkyl-substituted trans-stilbenes encapsulated in the hydrophobic cavity of an aqueous supramolecular organic host-guest complex. Compared with the solvated compounds, encapsulated trans-stilbenes have broader excited-state absorption spectra, excited-state lifetimes that are 3-4 times longer, and photoisomerization quantum yields that are 1.7-6.5 times lower in the restricted environment. The organic capsule disrupts the equilibrium structure and restricts torsional rotation around the central C═C double bond in the excited state, which is an important motion for the relaxation of trans-stilbene from S₁ to S₀. The location and identity of alkyl substituents play a significant role in determining the excited-state dynamics and photoisomerization quantum yields by tuning the relative crowding inside the capsule. The results are discussed in terms of distortions of the ground- and excited-state potential energy surfaces, including the topology of the S₁-S₀ conical intersection.

The full dynamics of energy relaxation in large organic molecules: from photo-excitation to solvent heating

In some molecular systems, such as nucleobases, polyenes or sunscreens, substantial amounts of photo-excitation energy are dissipated on a sub-picosecond time scale. Where does this energy go or among which degrees of freedom it is being distributed at such early times?

In some molecular systems, such as nucleobases, polyenes or the active ingredients of sunscreens, substantial amounts of photo-excitation energy are dissipated on a sub-picosecond time scale, raising questions such as: where does this energy go or among which degrees of freedom it is being distributed at such early times? Here we use transient absorption spectroscopy to track excitation energy dispersing from the optically accessible vibronic subsystem into the remaining vibrational subsystem of the solute and solvent. Monitoring the flow of energy during vibrational redistribution enables quantification of local molecular heating. Subsequent heat dissipation away from the solute molecule is characterized by classical thermodynamics and molecular dynamics simulations. Hence, we present a holistic approach that tracks the internal temperature and vibronic distribution from the act of photo-excitation to the restoration of the global equilibrium. Within this framework internal vibrational redistribution and vibrational cooling are emergent phenomena. We demonstrate the validity of the framework by examining a highly controversial example, carotenoids. We show that correctly accounting for the local temperature unambiguously explains their energetically and temporally congested spectral dynamics without the ad hoc postulation of additional ‘dark’ states. An immediate further application of this approach would be to monitor the excitation and thermal dynamics of pigment–protein systems.

Unified Mechanism of Oxygen Atom Transfer and Hydrogen Atom Transfer Reactions with a Triflic Acid-Bound Nonheme Manganese(IV)-Oxo Complex via Outer-Sphere Electron Transfer

Outer-sphere electron transfer from styrene, thioanisole, and toluene derivatives to a triflic acid (HOTf)-bound nonheme Mn(IV)-oxo complex, [(N4Py)Mn^IV(O)]²⁺-(HOTf)₂ (N4Py = N, N-bis(2-pyridylmethyl)- N-bis(2-pyridyl)methylamine), has been shown to be the rate-determining step of different types of redox reactions such as epoxidation, sulfoxidation, and hydroxylation of styrene, thioanisole, and toluene derivatives, respectively, by [(N4Py)Mn^IV(O)]²⁺-(HOTf)₂. The rate constants of HOTf-promoted epoxidation of all styrene derivatives with [(N4Py)Mn^IV(O)]²⁺ and electron transfer from electron donors to [(N4Py)Mn^V(O)]²⁺ exhibit a remarkably unified correlation with the driving force of outer-sphere electron transfer in light of the Marcus theory of electron transfer. The same electron-transfer driving force dependence is observed in the oxygen atom transfer from [(N4Py)Mn^IV(O)]²⁺-(HOTf)₂ to thioanisole derivatives as well as in the hydrogen atom transfer from toluene derivatives to [(N4Py)Mn^IV(O)]²⁺-(HOTf)₂. Thus, mechanisms of oxygen atom transfer (epoxidation and sulfoxidation) reactions of styrene and thioanisole derivatives and hydrogen atom transfer (hydroxylation) reactions of toluene derivatives by [(N4Py)Mn^IV(O)]²⁺-(HOTf)₂ have been unified for the first time as the same reaction pathway via outer-sphere electron transfer, followed by the fast bond-forming step, which exhibits the singly unified electron-transfer driving force dependence of the rate constants as outer-sphere electron-transfer reactions. In the case of the epoxidation of cis-stilbene by [(N4Py)Mn^IV(O)]²⁺-(HOTf)₂, the isomerization of cis-stilbene radical cation to trans-stilbene radical cation occurs after outer-sphere electron transfer from cis-stilbene to [(N4Py)Mn^IV(O)]²⁺-(HOTf)₂ to yield trans-stilbene oxide selectively, which is also taken as evidence for the occurrence of electron transfer in the acid-catalyzed epoxidation.

Mechanisms and Beyond: Elucidation of Fluxional Dynamics by Exchange NMR Spectroscopy

Detailed mechanistic information is crucial to our understanding of reaction pathways and selectivity. Dynamic exchange NMR techniques, in particular 2D exchange spectroscopy (EXSY) and its modifications, provide indispensable intricate information on the mechanisms of organic and inorganic reactions and other phenomena, for example, the dynamics of interfacial processes. In this Review, key results from exchange NMR studies of small molecules over the last few decades are systemised and discussed. After a brief introduction to the theory, the key types of dynamic processes are identified and fundamental examples given of intra- and intermolecular reactions, which, in turn, could involve, or not, bond-making and bond-breaking events. Following that logic, internal molecular rotation, intramolecular stereomutation and molecular recognition will first be considered because they do not typically involve bond breaking. Then, rearrangements, substitution-type reactions, cyclisations, additions and other processes affecting chemical bonds will be discussed. Finally, interfacial molecular dynamics and unexpected combinations of different types of fluxional processes will also be highlighted. How exchange NMR spectroscopy helps to identify conformational changes, coordination and molecular recognition processes as well as quantify reaction energy barriers and extract detailed mechanistic information by using reaction rate theory in conjunction with computational techniques will be shown.

Assessment of reactivities with explicit and implicit solvent models: QM/MM and gas-phase evaluation of three different Ag-catalysed furan ring formation routes

A representative reaction illustrates cases where strong solvent–solute interactions can be sufficiently well captured by continuum solvation model rendering QM/MM unnecessary.

EVB and polarizable MM study of energy relaxation in fluorine–acetonitrile reactions

Many-body effects can impact on rates of energy transfer from a ‘hot’ DF solute to acetonitrile solvent.

Post-transition state bifurcations induce dynamical detours in Pummerer-like reactions

A post-transition state bifurcation (PTSB) involved in a Pummerer-type rearrangement is characterized using density functional theory (DFT) calculations on potential energy stationary points and direct dynamics simulations. A sensitivity of the ratio of products produced via this PTSB to solvent dielectric constant is revealed and implications of such a dependence for selectivity control of organic reactions are discussed.

Relevance of the Resonance Junctions on the Arnold Web to Dynamical Tunneling and Eigenstate Delocalization

We study the competition and correspondence between the classical and quantum routes to intramolecular vibrational energy redistribution (IVR) in a three degrees of freedom model effective Hamiltonian. Specifically, we focus on the classical and the quantum dynamics near the resonance junctions on the Arnold web that are formed by an intersection of independent resonances. The regime of interest models the IVR dynamics from highly excited initial states near dissociation thresholds of molecular systems wherein both classical and purely quantum, involving dynamical tunneling, routes to IVR coexist. In the vicinity of a resonance junction, classical chaos is inevitably present, and hence one expects the quantum IVR pathways to have a strong classical component as well. We show that with increasing resonant coupling strengths the classical component of IVR leads to a transition from coherent dynamical tunneling to incoherent dynamical tunneling. Furthermore, we establish that the quantum IVR dynamics can be predicted based on the structures on the classical Arnold web. In addition, we investigate the nature of the highly excited eigenstates to identify the quantum signatures of the multiplicity-2 junctions. For the parameter regimes studies herein, by projecting the eigenstates onto the Arnold web, we find that eigenstates in the vicinity of the junctions are primarily delocalized due to dynamical tunneling.

Uncertainty quantification for quantum chemical models of complex reaction networks

For the quantitative understanding of complex chemical reaction mechanisms, it is, in general, necessary to accurately determine the corresponding free energy surface and to solve the resulting continuous-time reaction rate equations for a continuous state space. For a general (complex) reaction network, it is computationally hard to fulfill these two requirements. However, it is possible to approximately address these challenges in a physically consistent way. On the one hand, it may be sufficient to consider approximate free energies if a reliable uncertainty measure can be provided. On the other hand, a highly resolved time evolution may not be necessary to still determine quantitative fluxes in a reaction network if one is interested in specific time scales. In this paper, we present discrete-time kinetic simulations in discrete state space taking free energy uncertainties into account. The method builds upon thermo-chemical data obtained from electronic structure calculations in a condensed-phase model. Our kinetic approach supports the analysis of general reaction networks spanning multiple time scales, which is here demonstrated for the example of the formose reaction. An important application of our approach is the detection of regions in a reaction network which require further investigation, given the uncertainties introduced by both approximate electronic structure methods and kinetic models. Such cases can then be studied in greater detail with more sophisticated first-principles calculations and kinetic simulations.