Car-Parrinello molecular dynamics (CP-MD) simulations of the benzosemiquinone radical anion in aqueous solution have been performed at ambient conditions. Analysis of the trajectory shows not only extensive hydrogen bonding to the carbonyl oxygen atoms (ca. 4-5.6 water molecules depending on distance criteria), but also relatively long-lived "T-stacked" hydrogen bonds to the semiquinone pi-system. These results are discussed in the context of recent findings on semiquinone-protein interactions in photosynthetic reaction centers, and of EPR and vibration spectroscopical data for the aqueous system. Snapshots from the CP-MD trajectory are used for the first quantum chemical analyses of dynamical effects on electronic g-tensors, using cluster models and a recently developed density functional method. In particular, the effects of intermolecular hydrogen-bond dynamics on the g-tensor components are examined, in comparison with recent EPR and ENDOR studies.

Car-Parrinello molecular dynamics simulations of benzoquinone and benzosemiquinone radical anion in both aqueous solution and the gas phase have been carried out at ambient conditions. Hydrogen bonding is considerably more extensive to the anionic than to the neutral aqueous system. In addition to the conventional hydrogen bonding to the carbonyl oxygen atoms, T-stacked hydrogen bonding to the pi-system is statistically and energetically significant for the semiquinone anion but not for the neutral quinone. EPR g-tensors have been calculated at DFT level for snapshots taken at regular intervals from the gas-phase and solution semiquinone anion trajectories. Different criteria for extraction of semiquinone/water clusters from the solution trajectory give insight into the effects of different interactions on the g-tensor, as does correlation of the g-tensor with statistically significant hydrogen-bond configurations identified along the trajectories. Comparison of gas-phase and solution results indicates opposite directions of direct electronic and indirect structural influences of hydrogen bonding on g-tensors. Short-time oscillations in g(x) along the trajectory are due mainly to C-O bond vibrations.

No abstract

Based on Car-Parrinello ab initio molecular dynamics simulations of the benzosemiquinone radical anion in both aqueous solution and the gas phase, density functional calculations provide the currently most refined EPR hyperfine coupling (HFC) tensors of semiquinone nuclei and solvent protons. For snapshots taken at regular intervals from the molecular dynamics trajectories, cluster models with different criteria for inclusion of water molecules and an additional continuum solvent model are used to analyse the HFCs. These models provide a detailed picture of the effects of dynamics and of different intermolecular interactions on the spin-density distribution and HFC tensors. Comparison with static calculations allows an assessment of the importance of dynamical effects, and of error compensation in static DFT calculations. Solvent proton HFCs depend characteristically on the position relative to the semiquinone radical anion. A point-dipolar model works well for in-plane hydrogen-bonded protons but deviates from the quantum chemical values for out-of-plane hydrogen bonding.

The calculation of the g tensor of the main + NH 3 -CH-COO À radiation-induced radical in solid-state a-glycine presents a real challenge to computational methods. Density functional calculations of this spectroscopic property struggle with its small anisotropy and the zwitterionic nature of the amino acids in the crystal of this seemingly simple system. Here, several factors influencing the calculated g tensor are examined by comparing with experimental data. The extent of the molecular environment is varied in both a cluster and a periodic approach and dynamic calculations are performed to account for temperature effects. The latter does not necessarily lead to a better agreement with experiment than a static calculation. Application of a periodic approach is straightforward, but an all-electron scheme clearly is favorable. In a cluster approach, the selected basis set and density functional are of less importance, provided a hybrid functional is used to prevent cluster boundary effects. The applied spin-orbit coupling operators and proper treatment of the gauge origin of the magnetic vector potential also seem to be less critical than in other, similar molecular systems. But a careful selection of the cluster size proves to be essential for this glycine radical system. The calculated g tensor varies significantly with increasing cluster size, yielding only a good agreement with experiment when 5-7 glycine molecules in the immediate environment of the central glycine radical are incorporated. Further expansion of the cluster size can even lead to an essentially incorrect description of the radical in the condensed phase, indicating that bigger clusters can become unbalanced.