Nucleophilic Substitution Dynamics: Comparing Wave Packet Calculations with Experiment

Vibrational mode-specificity in the dynamics of the OH- + CH3I multi-channel reaction

We report a comprehensive characterization of the vibrational mode-specific dynamics of the OH- + CH3I reaction. Quasi-classical trajectory simulations are performed at four different collision energies on our previously-developed full-dimensional high-level ab initio potential energy surface in order to examine the impact of four different normal-mode excitations in the reactants. Considering the 11 possible pathways of OH- + CH3I, pronounced mode-specificity is observed in reactivity: In general, the excitations of the OH- stretching and CH stretching exert the greatest influence on the channels. For the SN2 and proton-abstraction products, the reactant initial attack angle and the product scattering angle distributions do not show major mode-specific features, except for SN2 at higher collision energies, where forward scattering is promoted by the CI stretching and CH stretching excitations. The post-reaction energy flow is also examined for SN2 and proton abstraction, and it is unveiled that the excess vibrational excitation energies rather transfer into the product vibrational energy because the translational and rotational energy distributions of the products do not represent significant mode-specificity. Moreover, in the course of proton abstraction, the surplus vibrational energy in the OH- reactant mostly remains in the H2O product owing to the prevailing dominance of the direct stripping mechanism.

Suppressing non-radiative decay of photochromic organic molecular systems in the strong coupling regime

The lifetimes of electronic excited states have a strong influence on the efficiency of organic solar cells. However, in some molecular systems a given excited state lifetime is reduced due to the non-radiative decay through conical intersections. Several strategies may be used to suppress this decay channel. The use of the strong light-matter coupling provided in optical nano-cavities is the focus of this paper. Here, we consider the meso-tert-butyl-4,4-difluoro-4-bora-3a,4a-diaza-s-indacene molecule (meso-tert-butyl-BODIPY) as a showcase of how strong and ultrastrong coupling might help in the development of organic solar cells. The meso-tert-butyl-BODIPY is known for its low fluorescence yield caused by the non-radiative decay through a conical intersection. However, we show here that, by considering this system within a cavity, the strong coupling can lead to significant changes in the multidimensional landscape of the potential energy surfaces of meso-tert-butyl-BODIPY, suppressing almost completely the decay of the excited state wave packet back to the ground state. By means of multi configuration electronic structure calculations and nuclear wave packet dynamics, the coupling with the cavity is analyzed in-depth to provide further insight of the interaction. By fine-tuning the cavity field strength and resonance frequency, we show that one can change the nuclear dynamics in the excited state, and control the non-radiative decay. This may lead to a faster and more efficient population transfer or the suppression of it.

Vibrational mode-specific dynamics of the F- + CH3CH2Cl multi-channel reaction

We investigate the mode-specific dynamics of the ground-state, C-Cl stretching (v₁₀), CH₂ wagging (v₇), sym-CH₂ stretching (v₁), and sym-CH₃ stretching (v₃) excited F^- + CH₃CH₂Cl(v_k = 0, 1) [k = 10, 7, 1, 3] → Cl^- + CH₃CH₂F (S_N2), HF + CH₃CHCl^-, FH⋯Cl^- + C₂H₄, and Cl^- + HF + C₂H₄ (E2) reactions using a full-dimensional high-level analytical global potential energy surface and the quasi-classical trajectory method. Excitation of the C-Cl stretching, CH₂ stretching, and CH₂/CH₃ stretching modes enhances the S_N2, proton abstraction, and FH⋯Cl^- and E2 channels, respectively. Anti-E2 dominates over syn-E2 (kinetic anti-E2 preference) and the thermodynamically-favored S_N2 (wider reactive anti-E2 attack angle range). The direct (a) S_N2, (b) proton abstraction, (c) FH⋯Cl^- + C₂H₄, (d) syn-E2, and (e) anti-E2 channels proceed with (a) back-side/backward, (b) isotropic/forward, (c) side-on/forward, (d) front-side/forward, and (e) back-side/forward attack/scattering, respectively. The HF products are vibrationally cold, especially for proton abstraction, and their rotational excitation increases for proton abstraction, anti-E2, and syn-E2, in order. Product internal-energy and mode-specific vibrational distributions show that CH₃CH₂F is internally hot with significant C-F stretching and CH₂ wagging excitations, whereas C₂H₄ is colder. One-dimensional Gaussian binning technique is proved to solve the normal mode analysis failure caused by methyl internal rotation.

Q-Band relaxation in chlorophyll: new insights from multireference quantum dynamics

The ultrafast relaxation within the Q-bands of chlorophyll plays a crucial role in photosynthetic light-harvesting. We investigate this process via nuclear and electronic quantum dynamics on multireference potential energy surfaces.

Fifty years of nucleophilic substitution in the gas phase

Bimolecular nucleophilic substitution ( S N 2 ) reactions have become a model system for the investigation of structure-reactivity relationships, stereochemistry, solvent influences, and detailed atomistic dynamics. In this review, the progress during five decades of experimental and theoretical research on gas phase S N 2 reactions is discussed. Many advancements of the employed methods have led to a tremendous increase in our understanding of the properties and the dynamics of these reactions. For reactions involving six atoms a quantitative agreement of the differential reactive scattering cross sections has already been achieved, in the future it is expected that even larger polyatomic reactions systems become tractable. Furthermore, studies with higher precision, improved reactant control, and a more accurate theoretical treatment of quantum effects are envisioned.

Mode-Specific Quasiclassical Dynamics of the F- + CH3I SN2 and Proton-Transfer Reactions

Mode-specific quasiclassical trajectory computations are performed for the F^- + CH₃I( v _k = 0, 1) S_N2 and proton-transfer reactions at nine different collision energies in the range of 1.0-35.3 kcal/mol using a full-dimensional high-level ab initio analytical potential energy surface with ground-state and excited CI stretching ( v₃), CH₃ rocking ( v₆), CH₃ umbrella ( v₂), CH₃ deformation ( v₅), CH symmetric stretching ( v₁), and CH asymmetric stretching ( v₄) initial vibrational modes. Millions of trajectories provide statistically definitive mode-specific cross sections, opacity functions, scattering angle distributions, and product internal energy distributions. The excitation functions reveal slight vibrational S_N2 inversion inhibition/enhancement at low/high collision energies ( E_coll), whereas large decaying-with- E_coll vibrational enhancement effects for the S_N2 retention (double inversion) and proton-transfer channels. The most efficient vibrational enhancement is found by exciting the CI stretching (high E_coll) for S_N2 inversion and the CH stretching modes (low E_coll) for double inversion and proton transfer. Mode-specific effects do not show up in the scattering angle distributions and do blue-shift the hot/cold S_N2/proton-transfer product internal energies.

Stretching vibration is a spectator in nucleophilic substitution

How chemical reactions are influenced by reactant vibrational excitation is a long-standing question at the core of chemical reaction dynamics. In reactions of polyatomic molecules, where the Polanyi rules are not directly applicable, certain vibrational modes can act as spectators. In nucleophilic substitution reactions, CH stretching vibrations have been considered to be such spectators. While this picture has been challenged by some theoretical studies, experimental insight has been lacking. We show that the nucleophilic substitution reaction of F^- with CH₃I is minimally influenced by an excitation of the symmetric CH stretching vibration. This contrasts with the strong vibrational enhancement of the proton transfer reaction measured in parallel. The spectator behavior of the stretching mode is supported by both quasi-classical trajectory simulations and the Sudden Vector Projection model.

Conservation of direct dynamics in sterically hindered SN2/E2 reactions

Nucleophilic substitution (S_N2) and base-induced elimination (E2), two indispensable reactions in organic synthesis, are commonly assumed to proceed under stereospecific conditions. Understanding the way in which the reactants pre-orient in these reactions, that is its stereodynamics, is essential in order to achieve a detailed atomistic picture and control over such processes. Using crossed beam velocity map imaging, we study the effect of steric hindrance in reactions of Cl^- and CN^- with increasingly methylated alkyl iodides by monitoring the product ion energy and scattering angle. For both attacking anions the rebound mechanism, indicative of a direct S_N2 pathway, is found to contribute to the reaction at high relative collision energies despite being increasingly hindered. An additional forward scattering mechanism, ascribed to a direct E2 reaction, also contributes at these energies. Inspection of the product energy distributions confirms the direct and fast character of both mechanisms as opposed to an indirect reaction mechanism which leads to statistical energy redistribution in the reaction complex. This work demonstrates that nonstatistical dynamics and energetics govern S_N2 and E2 pathways even in sterically hindered exchange reaction systems.

Two New Methods To Generate Internal Coordinates for Molecular Wave Packet Dynamics in Reduced Dimensions

The curse of dimensionality still remains as the central challenge of molecular quantum dynamical calculations. Either compromises on the accuracy of the potential landscape have to be made or methods must be used that reduce the dimensionality of the configuration space of molecular systems to a low dimensional one. For dynamic approaches such as grid-based wave packet dynamics that are confined to a small number of degrees of freedom this dimensionality reduction can become a major part of the overall problem. A common strategy to reduce the configuration space is by selection of a set of internal coordinates using chemical intuition. We devised two methods that increase the degree of automation of the dimensionality reduction as well as replace chemical intuition by more quantifiable criteria. Both methods reduce the dimensionality linearly and use the intrinsic reaction coordinate as guidance. The first one solely relies on the intrinsic reaction coordinate (IRC), whereas the second one uses semiclassical trajectories to identify the important degrees of freedom.

Mode-Specific SN2 Reaction Dynamics

Despite its importance in chemistry, the microscopic dynamics of bimolecular nucleophilic substitution (SN2) reactions is still not completely elucidated. In this publication, the dynamics of a prototypical SN2 reaction (F(-) + CH3Cl → CH3F + Cl(-)) is investigated using a high-dimensional quantum mechanical model on an accurate potential energy surface (PES) and further analyzed by quasi-classical trajectories on the same PES. While the indirect mechanism dominates at low collision energies, the direct mechanism makes a significant contribution. The reactivity is found to depend on the specific reactant vibrational mode excitation. The mode specificity, which is more prevalent in the direct reaction, is rationalized by a transition-state-based model.

Quantum Dynamics of a Photochemical Bond Cleavage Influenced by the Solvent Environment: A Dynamic Continuum Approach

In every day chemistry, solvents are used to influence the outcome of chemical synthesis. Electrostatic effects stabilize polar configurations during the reaction and in addition dynamic solvent effects can emerge. How the dynamic effects intervene on the ultrafast time scale is in the focus of this theoretical study. We selected the photoinduced bond cleavage of Ph2CH-PPh3(+) for which the electrostatic interactions are negligible. Elaborate ultrafast pump-probe studies already exist and serve as a reference. We compared quantum dynamical simulations with and without environment and noticed the necessity to model the influence of the solvent cage on the reactive motions of the solute. The frictional force induced by the dynamic viscosity of the solvent is implemented in the quantum mechanical formalism with a newly developed approach called the dynamic continuum ansatz. Only when the environment is included are the experimentally observed products reproduced on the subpicosecond time scale.