Search

Previously, it has been demonstrated that external electric fields may be used to exert control over chemical reactivity. In this study, the impact of a strong, nonresonant IR field (1064 nm) on the photoisomerization of cis- stilbene is investigated in cyclohexane solution. The design of a suitable reaction vessel for characterization of this effect is presented. The electric field supplied by the pulsed, near-IR radiation (ε l = 4.5 × 10 7 V/cm) enhances the cis → trans photoisomerization yield at the red edge of the absorption spectrum (wavelengths between 337 and 340 nm). Within the microliter focal volume, up to 75% of all cis- stilbene molecules undergo isomerization to trans- stilbene in the strong electric-field environment, indicating a significant increase relative to the 35% yield of trans- stilbene under field-free conditions. This result correlates with a 1-3% enhancement in the trans- stilbene concentration throughout the bulk solution. Theoretical analysis suggests that the observed change is the result of dynamic Stark shifting of the ground and first excited states, leading to a significant redshift in cis -stilbene's absorption spectrum. The predicted increase in the absorption cross section in this range of excitation wavelengths is qualitatively consistent with the experimental increase in trans -stilbene production.

Force can induce remarkable non-destructive transformations along a polymer, but we have a limited understanding of the energy transduction and product distribution in tandem mechanochemical reactions. Ladderanes consist of multiple fused cyclobutane rings and have recently been used as monomeric motifs to develop polymers that drastically change their properties in response to force. Here we show that [4]-ladderane always exhibits 'all-or-none' cascade mechanoactivations and the same stereochemical distribution of the generated dienes under various conditions and within different polymer backbones. Transition state theory fails to capture the reaction kinetics and explain the observed stereochemical distributions. Ab initio steered molecular dynamics reveals unique non-equilibrium dynamic effects: energy transduction from the first cycloreversion substantially accelerates the second cycloreversion, and bifurcation on the force-modified potential energy surface leads to the product distributions. Our findings illustrate the rich chemistry in closely coupled multi-mechanophores and an exciting potential for effective energy transduction in tandem mechanochemical reactions.

The relaxation dynamics of thymine and its derivatives thymidine and thymidine monophosphate are studied using time-resolved photoelectron spectroscopy applied to a water microjet. Two absorption bands are studied; the first is a bright ππ* state which is populated using tunable-ultraviolet light in the range 4.74-5.17 eV and probed using a 6.20 eV probe pulse. By reversing the order of these pulses, a band containing multiple ππ* states is populated by the 6.20 eV pulse and the lower energy pulse serves as the probe. The lower lying ππ* state is found to decay in ∼400 fs in both thymine and thymidine independent of pump photon energy, while thymidine monophosphate decays vary from 670 to 840 fs with some pump energy dependence. The application of a computational quantum mechanical/molecular mechanical scheme at the XMS-CASPT2//CASSCF/AMBER level of theory suggests that conformational differences existing between thymidine and thymidine monophosphate in solution account for this difference. The higher lying ππ* band is found to decay in ∼600 fs in all three cases, but it is only able to be characterized when the 5.17 eV probe pulse is used. Notably, no long-lived signal from an nπ* state can be identified in either experiment on any of the three molecules.

Ultrafast proton migration and isomerization are key processes for acetylene and its ions. However, the mechanism for ultrafast isomerization of acetylene [HCCH] 2+ to vinylidene [H 2 CC] 2+ dication remains nebulous. Theoretical studies show a large potential barrier ( > 2 eV) for isomerization on low-lying dicationic states, implying picosecond or longer isomerization timescales. However, a recent experiment at a femtosecond X-ray free-electron laser suggests sub-100 fs isomerization. Here we address this contradiction with a complete theoretical study of the dynamics of acetylene dication produced by Auger decay after X-ray photoionization of the carbon atom K shell. We find no sub-100 fs isomerization, while reproducing the salient features of the time-resolved Coulomb imaging experiment. This work resolves the seeming contradiction between experiment and theory and also calls for careful interpretation of structural information from the widely applied Coulomb momentum imaging method.The timescale of isomerization in molecules involving ultrafast migration of constituent atoms is difficult to measure. Here the authors report that sub-100 fs isomerization time on acetylene dication in lower electronic states is not possible and point to misinterpretation of recent experimental results.

We describe a flexible and broadly applicable energy refinement method, "nebterpolation," for identifying and characterizing the reaction events in a molecular dynamics (MD) simulation. The new method is applicable to ab initio simulations with hundreds of atoms containing complex and multimolecular reaction events. A key aspect of nebterpolation is smoothing of the reactive MD trajectory in internal coordinates to initiate the search for the reaction path on the potential energy surface. We apply nebterpolation to analyze the reaction events in an ab initio nanoreactor simulation that discovers new molecules and mechanisms, including a C-C coupling pathway for glycolaldehyde synthesis. We find that the new method, which incorporates information from the MD trajectory that connects reactants with products, produces a dramatically distinct set of minimum energy paths compared to existing approaches that start from information for the reaction end points alone. The energy refinement method described here represents a key component of an emerging simulation paradigm where molecular dynamics simulations are applied to discover the possible reaction mechanisms.

The development of accurate molecular mechanics force fields is a significant challenge that must be addressed for the continued success of molecular simulation. We developed the ForceBalance method to automatically derive accurate force field parameters using flexible combinations of experimental and theoretical reference data. The method is demonstrated in the parametrization of two rigid water models, yielding new parameter sets (TIP3P-FB and TIP4P-FB) that accurately describe many physical properties of water.