Showing papers by "David E. Manolopoulos published in 2018"

Spin-selective electron transfer reactions of radical pairs: Beyond the Haberkorn master equation

Abstract: 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.

Analytic continuation of Wolynes theory into the Marcus inverted regime.

Abstract: The Wolynes theory of electronically nonadiabatic reaction rates [P. G. Wolynes, J. Chem. Phys. 87, 6559 (1987)] is based on a saddle point approximation to the time integral of a reactive flux autocorrelation function in the nonadiabatic (golden rule) limit. The dominant saddle point is on the imaginary time axis at tsp=iλspℏ, and provided λsp lies in the range −β/2≤λsp≤β/2, it is straightforward to evaluate the rate constant using information obtained from an imaginary time path integral calculation. However, if λsp lies outside this range, as it does in the Marcus inverted regime, the path integral diverges. This has led to claims in the literature that Wolynes theory cannot describe the correct behaviour in the inverted regime. Here we show how the imaginary time correlation function obtained from a path integral calculation can be analytically continued to λsp<−β/2, and the continuation used to evaluate the rate in the inverted regime. Comparison with exact golden rule results for a spin-boson model ...

On the low magnetic field effect in radical pair reactions

Abstract: Radical pair recombination reactions are known to be sensitive to the application of both low and high magnetic fields. The application of a weak magnetic field reduces the singlet yield of a singlet-born radical pair, whereas the application of a strong magnetic field increases the singlet yield. The high field effect arises from energy conservation: when the magnetic field is stronger than the sum of the hyperfine fields in the two radicals, S → T± transitions become energetically forbidden, thereby reducing the number of pathways for singlet to triplet interconversion. The low field effect arises from symmetry breaking: the application of a weak magnetic field lifts degeneracies among the zero field eigenstates and increases the number of pathways for singlet to triplet interconversion. However, the details of this effect are more subtle and have not previously been properly explained. Here we present a complete analysis of the low field effect in a radical pair containing a single proton and in a radical pair in which one of the radicals contains a large number of hyperfine-coupled nuclear spins. We find that the new transitions that occur when the field is switched on are between S and T0 in both cases, and not between S and T± as has previously been claimed. We then illustrate this result by using it in conjunction with semiclassical spin dynamics simulations to account for the observation of a biphasic-triphasic-biphasic transition with increasing magnetic field strength in the magnetic field effect on the time-dependent survival probability of a photoexcited carotenoid-porphyrin-fullerene radical pair.

Spin-selective electron transfer reactions of radical pairs: beyond the Haberkorn master equation.

Abstract: 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 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 (HEOM) 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.

Simple and Accurate Method for Central Spin Problems.

Abstract: We describe a simple quantum mechanical method that can be used to obtain accurate numerical results over long timescales for the spin correlation tensor of an electron spin that is hyperfine coupled to a large number of nuclear spins. This method does not suffer from the statistical errors that accompany a Monte Carlo sampling of the exact eigenstates of the central spin Hamiltonian obtained from the algebraic Bethe ansatz, or from the growth of the truncation error with time in the time-dependent density matrix renormalization group (TDMRG) approach. As a result, it can be applied to larger central spin problems than the algebraic Bethe ansatz, and for longer times than the TDMRG algorithm. It is therefore an ideal method to use to solve central spin problems, and we expect that it will also prove useful for a variety of related problems that arise in a number of different research fields.

On the low magnetic field effect in radical pair reactions

Abstract: Radical pair recombination reactions are known to be sensitive to the application of both low and high magnetic fields. The application of a weak magnetic field reduces the singlet yield of a singlet-born radical pair, whereas the application of a strong magnetic field increases the singlet yield. The high field effect arises from energy conservation: when the magnetic field is stronger than the sum of the hyperfine fields in the two radicals, ${\rm S}\to {\rm T}_{\pm}$ transitions become energetically forbidden, thereby reducing the number of pathways for singlet to triplet interconversion. The low field effect arises from symmetry breaking: the application of a weak magnetic field lifts degeneracies among the zero field eigenstates and increases the number of pathways for singlet to triplet interconversion. However, the details of this effect are more subtle, and have not previously been properly explained. Here we present a complete analysis of the low field effect in a radical pair containing a single proton, and in a radical pair in which one of the radicals contains a large number of hyperfine-coupled nuclear spins. We find that the new transitions that occur when the field is switched on are between ${\rm S}$ and ${\rm T}_0$ in both cases, and not between ${\rm S}$ and ${\rm T}_{\pm}$ as has previously been claimed. We then illustrate this result by using it in conjunction with semiclassical spin dynamics simulations to account for the observation of a biphasic--triphasic--biphasic transition with increasing magnetic field strength in the magnetic field effect on the time-dependent survival probability of a photoexcited carotenoid-porphyrin-fullerene radical pair.