Recently, there has been much interest in the chirality-induced spin selectivity effect, whereby electron spin polarization, which is dependent on molecular chirality, is produced in electrode− molecule electron transfer processes. Naturally, one might consider if a similar effect can be observed in simple molecular charge transfer reactions, for example, in light-induced electron transfer from an electron donor to an electron acceptor. In this work, I explore the effect of electron transfer on spins in chiral single radicals and chiral radical pairs using Nakajima−Zwanzig theory. In these cases, chirality, in conjuction with spin−orbit coupling, does not lead to spin polarization, but instead, the electron transfer generates quantum coherence between spins states. In principle, this chirality-induced spin coherence could manifest in a range of experiments, and in particular, I demonstrate that the out of phase electron spin echo envelope modulation pulse electron paramagnetic resonance experiment would be able to detect this effect in oriented radical pairs.

Radical pair recombination reactions are normally described using a quantum mechanical master equation for the electronic and nuclear spin density operator. The electron spin state selective (singlet and triplet) recombination processes are described with a Haberkorn reaction term in this master equation. Here we consider a general spin state selective electron transfer reaction of a radical pair and use Nakajima-Zwanzig theory to derive the master equation for the spin density operator, thereby elucidating the relationship between non-adiabatic reaction rate theory and the Haberkorn reaction term. A second order perturbation theory treatment of the diabatic coupling naturally results in the Haberkorn master equation with an additional reactive scalar electron spin coupling term. This term has been neglected in previous spin chemistry calculations, but we show that it will often be quite significant. We also show that beyond the second order in perturbation theory, i.e., beyond the Fermi golden rule limit, an additional reactive singlet-triplet dephasing term appears in the master equation. A closed form expression for the reactive scalar electron spin coupling in terms of the Marcus theory parameters that determine the singlet and triplet recombination rates is presented. By performing simulations of radical pair reactions with the exact hierarchical equations of motion method, we demonstrate that our master equations provide a very accurate description of radical pairs undergoing spin-selective non-adiabatic electron transfer reactions. The existence of a reactive electron spin coupling may well have implications for biologically relevant radical pair reactions such as those which have been suggested to play a role in avian magnetoreception.

Semiclassical methods cannot accurately simulate magnetic field effects relevant to avian magnetoreception, which may therefore deserve the label “quantum biology”.

Here we propose a mechanism by which spin-polarization can be generated dynamically in chiral molecular systems undergoing photoinduced electron transfer. The proposed mechanism explains how spin-polarization emerges in systems where charge transport is dominated by incoherent hopping, mediated by spin–orbit and electronic exchange couplings through an intermediate charge transfer state. We derive a simple expression for the spin-polarization that predicts a nonmonotonic temperature dependence, consistent with recent experiments, and a maximum spin-polarization that is independent of the magnitude of the spin–orbit coupling. We validate this theory using approximate quantum master equations and the numerically exact hierarchical equations of motion. The proposed mechanism of chirality induced spin selectivity should apply to many chiral systems, and the ideas presented here have implications for the study of spin transport at temperatures relevant to biology and provide simple principles for the molecular control of spins in fluctuating environments.

The standard quantum mechanical expressions for the singlet and triplet survival probabilities and product yields of a radical pair recombination reaction involve a trace over the states in a combined electronic and nuclear spin Hilbert space. If this trace is evaluated deterministically, by performing a separate time-dependent wavepacket calculation for each initial state in the Hilbert space, the computational effort scales as O(Zlog⁡Z), where Z is the total number of nuclear spin states. Here we show that the trace can also be evaluated stochastically, by exploiting the properties of spin coherent states. This results in a computational effort of O(MZlog⁡Z), where M is the number of Monte Carlo samples needed for convergence. Example calculations on a strongly coupled radical pair with Z>10 show that the singlet yield can be converged to graphical accuracy using just M=200 samples, resulting in a speed up by a factor of >5000 over a standard deterministic calculation. We expect that this factor will greatly facilitate future quantum mechanical simulations of a wide variety of radical pairs of interest in chemistry and biology.