Spin Dynamics of Radical Pairs Using the Stochastic Schrödinger Equation in MolSpin.

The chemical reactivity of radical pairs is strongly influenced by the interactions of electronic and nuclear spins. A detailed understanding of these effects requires a quantum description of the spin dynamics that considers spin-dependent reaction rates, interactions with external magnetic fields, spin-spin interactions, and the loss of spin coherence caused by coupling to a fluctuating environment. Modeling real chemical and biochemical systems, which frequently involve radicals with multinuclear spin systems, poses a severe computational challenge. Here, we implement a method based on the stochastic Schrödinger equation in the software package MolSpin. Large electron-nuclear spin systems can be simulated efficiently, with asymmetric spin-selective recombination reactions, anisotropic hyperfine interactions, and nonzero exchange and dipolar couplings. Spin-relaxation can be modeled using the stochastic time-dependence of spin interactions determined by molecular dynamics and quantum chemical calculations or by allowing rate coefficients to be explicitly time-dependent. The flexibility afforded by this approach opens new avenues for exploring the effects of complex molecular motions on the spin dynamics of chemical transformations.

Isotope Substitution Effects on the Magnetic Compass Properties of Cryptochrome-Based Radical Pairs: A Computational Study.

The biophysical mechanism of the magnetic compass sense of migratory songbirds is thought to rely on the photochemical reactions of flavin-containing radical pairs in cryptochrome proteins located in the birds' eyes. A consequence of this hypothesis is that the effect of the Earth's magnetic field on the quantum yields of reaction products should be sensitive to isotopic substitutions that modify the hyperfine interactions in the radicals. In this report, we use spin dynamics simulations to explore the effects of 1H → 2H, 12C → 13C, and 14N → 15N isotopic substitutions on the functioning of cryptochrome 4a as a magnetic direction sensor. Two main conclusions emerge. (1) Uniform deuteration of the flavin chromophore appears to be the best way to boost the anisotropy of the magnetic field effect and to change its symmetry. (2) 13C substitution of three of the 12 flavin carbons, in particular C4, C4a, and C8α, seems to be the best recipe for attenuating the anisotropy. These predictions should give insight into the factors that control the magnetic sensitivity once spectroscopic techniques are available for measuring magnetic field effects on oriented protein samples.