Rates of Molecular Vibrational Energy Transfer in Organic Solutions

Synthesis of Bis(amino acids) Containing the Styryl-cyclobutane Core by Photosensitized [2+2]-Cross-cycloaddition of Allylidene-5(4H)-oxazolones

The irradiation of 2-aryl-4-(E-3'-aryl-allylidene)-5(4H)-oxazolones 1 with blue light (456 nm) in the presence of [Ru(bpy)₃](BF₄)₂ (bpy = 2,2'-bipyridine, 5% mol) gives the unstable cyclobutane-bis(oxazolones) 2 by [2+2]-photocycloaddition of two oxazolones 1. Each oxazolone contributes to the formation of 2 with a different C=C bond, one of them reacting through the exocyclic C=C bond, while the other does so through the styryl group. Treatment of unstable cyclobutanes 2 with NaOMe/MeOH produces the oxazolone ring opening reaction, affording stable styryl-cyclobutane bis(amino acids) 3. The reaction starts with formation of the T₁ excited state of the photosensitizer ³[Ru*(bpy)₃]²⁺, which reacts with S₀ of oxazolones 1 through energy transfer to give the oxazolone T₁ state ³(oxa*)-1, which is the reactive species and was characterized by transient absorption spectroscopy. Measurement of the half-life of ³(oxa*)-1 for 1a, 1b and 1d shows large values for 1a and 1b (10-12 μs), while that of 1d is shorter (726 ns). Density functional theory (DFT) modeling displays strong structural differences in the T₁ states of the three oxazolones. Moreover, study of the spin density of T₁ state ³(oxa*)-1 provides clues to understanding the different reactivity of 4-allylidene-oxazolones described here with respect to the previously reported 4-arylidene-oxazolones.

Influence of a Methyl Group on the Unidirectional Flow of Vibrational Energy in an Adenine-Thymine Base Pair

The role of a methyl group in intramolecular vibrational energy redistribution (IVR) of the hydrogen-bonded adenine-thymine base pair has been studied using classical dynamics procedures. Energy transferred to the doorway bond thymine-NH from the vibrationally excited H₂O(v) efficiently redistributes among various bonds of the base pair through vibration-to-vibration coupling, depositing a large fraction of the available energy in the terminal bond adenine-NH. On the other hand, the extent of energy flow in the reverse direction from the excited adenine-NH to thymine-NH is insignificant, indicating IVR in adenine-thymine resulting from the intermolecular interaction with a vibrationally excited H₂O molecule, is direction-specific. The unidirectional flow is due to the coupling of stretch-torsion vibrations of a methyl group with conjugated bonds on the thymine ring, when the methyl rotor is present and is adjacent to the vibrationally excited thymine-NH. The insignificance of energy flow from the terminal-to-terminal bond in the reverse direction is attributed to the absence of a methyl group on the adenine moiety, even though the molecule has many CC and CN bonds coupled to their neighbors.

Intramolecular vibrational energy redistribution in nucleobases: Excitation of NH stretching vibrations in adenine–uracil + H2O

Redistribution of vibrational energy in the adenine–uracil base pair is studied when the base pair undergoes an intermolecular interaction with an overtone-bending vibration excited H2O(2[Formula: see text]bend) molecule. Energy transfer is calculated using the structural information obtained from density functional theory in the solution of the equations of motion. Intermolecular vibrational energy transfer (VET) from H2O(2[Formula: see text]bend) to the uracil–NH stretching mode is efficient and rapidly followed by intramolecular vibrational energy redistribution (IVR) resulting from coupling between vibrational modes. An important pathway is IVR carrying energy to the NH-stretching mode of the adenine moiety in a subpicosecond scale, the energy build-up being sigmoidal, when H2O interacts with the uracil–NH bond. The majority of intermolecular hydrogen bonds between the base pair and H2O are weakened but unbroken during the ultrafast energy redistribution period. Lifetimes of intermolecular HB are on the order of 0.5 ps. The efficiency of IVR in the base pair is due to near-resonance between coupled CC and CN vibrations. The resonance also exists between the frequencies of H2O bend and NH stretch, thus facilitating VET. When H2O interacts with the NH bond at the adenine end of the base pair, energy flow in the reverse direction to the uracil–NH stretch is negligible, the unidirectionality discussed in terms of the effects of uracil CH stretches. The energy distributed in the CH bonds is found to be significant. The IVR process is found to be nearly temperature independent between 200 and 400 K.

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.

Impact of dimensionality and confinement on reaction dynamics and thermodynamics within 1D and 2D nanostructures

Confinement has been shown to contribute to the dynamics of small molecules within nanoscale hydrophobic or hydrophilic cavities. Enclosure within a confined space can also influence energy transfer pathways, such as the enhancement of fluorescence over thermal relaxation. In this paper, the effect of confinement on the thermodynamic properties and reaction kinetics of small hydrophobic molecules confined in a soft polymeric template is detailed. A quasi-elastic neutron scattering experiment identified a substantial decrease in translational diffusion of pyrrole after solubilization within a hydrophobic cavity. This decrease in mobility is due to pyrrole's closer packing and increased density under confinement vs the bulk liquid. The decreased mobility and increased density explain the spontaneous polymerization reaction of pyrrole observed within the cavity. The precise characterization of the polymerization kinetics under confinement found that the reaction is independent of pyrrole concentration, consistent with the close packing density. Kinetic data also show that confinement dimensionality finds a thermodynamic expression in the transition state entropy. The dynamics and kinetics experiments reported here offer rare empirical insight into the important influence that cavity geometry places on the reactions they host.

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.

Understanding and Predicting Post H-Atom Abstraction Selectivity through Reactive Mode Composition Factor Analysis

The selective functionalization of C-H bonds is one of the Grails of synthetic chemistry. In this work, we demonstrate that the selectivity toward fast hydroxylation or radical diffusion (known as the OH-rebound and dissociation mechanisms) following H-atom abstraction (HAA) from a substrate C-H bond by high-valent iron-oxo oxidants is already encoded in the HAA step when the post-HAA barriers are much lower than the preceding one. By applying the reactive mode composition factor (RMCF) analysis, which quantifies the kinetic energy distribution (KED) at the reactive mode (RM) of transition states, we show that reactions following the OH-rebound coordinate concentrate the RM kinetic energy on the motion of the reacting oxygen atom and the nascent substrate radical, whereas reactions following the dissociation channel localize most of their kinetic energy in H-atom motion. These motion signatures serve to predict the post-HAA selectivity, and since KED is affected by the free energy of reaction and asynchronicity (factor η) of HAA, we show that bimolecular HAA reactions in solution that are electron transfer-driven and highly exergonic have the lowest fraction of KED on the transferred H-atom and the highest chance to follow rebound hydroxylation. Finally, the RMCF analysis predicts that the H/D primary kinetic isotope effect can serve as a probe for these mechanisms, as confirmed in virtually all reported examples in the literature.

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.

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.

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.

The effect of coordination of alkanes, Xe and CO2 (η1-OCO) on changes in spin state and reactivity in organometallic chemistry: a combined experimental and theoretical study of the photochemistry of CpMn(CO)3

A combined experimental and theoretical study is presented of several ligand addition reactions of the triplet fragment 3CpMn(CO)2 formed upon photolysis of CpMn(CO)3. Experimental data are provided for reactions in n-heptane and perfluoromethylcyclohexane (PFMCH), as well as in PFMCH doped with C2H6, Xe and CO2. In PFMCH we find that the conversion of 3CpMn(CO)2 to 1CpMn(CO)2(PFMCH) is much slower (τ = 18 (±3) ns) than the corresponding reactions in conventional alkanes (τ = 111 (±10) ps). We measure the effect of the coordination ability by doping PFMCH with alkane, Xe and CO2; these doped ligands form the corresponding singlet adducts with significantly variable formation rates. The reactivity as measured by the addition timescale follows the order 1CpMn(CO)2(C5H10) (τ = 270 (±10) ps) > 1CpMn(CO)2Xe (τ = 3.9 (±0.4) ns) ∼ 1CpMn(CO)2(CO2) (τ = 4.7 (±0.5) ns) > 1CpMn(CO)2(C7F14) (τ = 18 (±3) ns). Electronic structure theory calculations of the singlet and triplet potential energy surfaces and of their intersections, together with non-adiabatic statistical rate theory, reproduce the observed rates semi-quantitatively. It is shown that triplet adducts of the ligand and 3CpMn(CO)2 play a role in the kinetics, and account for the variable timescales observed experimentally.