Competitive reaction pathways in vibrationally induced photodissociation of H2SO4

Vibrationally induced photodissociation of sulfuric acid into H2O + SO3 is investigated based on reactive molecular dynamics (MD) simulations. Multisurface adiabatic reactive MD simulations allow us to follow both, H-transfer and water elimination after excitation of the ν9 OH-stretching mode. Analysis of several thousand trajectories finds that the H2O and SO3 fragments have distinct final state distributions with respect to translational, rotational, and vibrational degrees of freedom. Rotational distributions peak at quantum numbers j ≤ 5 for water and j ≈ 60 for SO3. The final state distributions should be useful in identifying products in forthcoming experiments. Based on the MD trajectories, a kinetic scheme has been developed which is able to explain most of the trajectory data and suggests that IVR is very rapid. Typical lifetimes of the excited complex range from several 10 picoseconds to hundreds of nanoseconds, depending on the excitation level. Including temperature and pressure profiles characteristic for the stratosphere in the kinetic model shows that excitations higher than ν9 = 4 can significantly contribute to the photolysis rate. This extends and specifies earlier work in that multi-level modeling is required to understand the significance of vibrationally induced decomposition pathways of sulfuric acid in the middle atmosphere.

CO-dynamics in the active site of cytochrome c oxidase

The transfer of CO from heme a3 to the Cu(B) site in Cytochrome c oxidase (CcO) after photolysis is studied using molecular dynamics simulations using an explicitly reactive, parametrized potential energy surface based on density functional theory calculations. After photodissociation from the heme-Fe, the CO ligand rebinds to the Cu(B) site on the sub-picosecond time scale. Depending on the simulation protocol the characteristic time ranges from 260 fs to 380 fs which compares with an estimated 450 fs from experiment based on the analysis of the spectral changes as a function of time delay after the photodissociating pulse. Following photoexcitation ≈90% of the ligands are found to rebind to either the Cu(B) (major component, 85%) or the heme-Fe (minor component, 2%) whereas about 10% remain in an unbound state. The infrared spectra of unbound CO in the active site is broad and featureless and no appreciable shift relative to gas-phase CO is found, which is in contrast to the situation in myoglobin. These observations explain why experimentally, unbound CO in the binuclear site of CcO has not been found as yet.

Vibrationally induced dissociation of sulfuric acid (H2SO4)

One of the important reactive steps in Earth's atmosphere is the decomposition of H(2)SO(4) to H(2)O and SO(3). However, because the UV spectrum of H(2)SO(4) was not found up to 140 nm, alternative mechanisms, including vibrationally induced dissociation, were proposed. Using adiabatic reactive molecular dynamics (ARMD) simulations with validated force fields for the product and educt channels, it is shown through explicit atomistic simulation that by exciting the ν(9) (OH-stretching-) mode, photodissociation can occur on the picosecond time scale. With the potential energy surfaces used in the present work, ν(9) = 4 is sufficient for this process. From a statistically significant number of trajectories (several thousands), vibrationally induced dissociation times are found to follow Gamma-distributions with most likely reaction times between 40 and 200 ps by depositing energies ranging from 40 to 60 kcal/mol, corresponding to 4 and 6 vibrational quanta in the OH stretching vibration. Because ARMD simulations allow multiple and long-time simulations, both nonstatistical, impulsive H-transfer and statistical, IVR-regimes of the decomposition reaction can be discussed in detail at an atomistic level.