Showing 1–17 of 17 results for author: Lindoy, L P

Quantum subspace expansion approach for simulating dynamical response functions of Kitaev spin liquids

Authors: Chukwudubem Umeano, François Jamet, Lachlan P. Lindoy, Ivan Rungger, Oleksandr Kyriienko

Abstract: We develop a quantum simulation-based approach for studying properties of strongly correlated magnetic materials at increasing scale. We consider a paradigmatic example of a quantum spin liquid (QSL) state hosted by the honeycomb Kitaev model, and use a trainable symmetry-guided ansatz for preparing its ground state. Applying the tools of quantum subspace expansion (QSE), Hamiltonian operator appr… ▽ More We develop a quantum simulation-based approach for studying properties of strongly correlated magnetic materials at increasing scale. We consider a paradigmatic example of a quantum spin liquid (QSL) state hosted by the honeycomb Kitaev model, and use a trainable symmetry-guided ansatz for preparing its ground state. Applying the tools of quantum subspace expansion (QSE), Hamiltonian operator approximation, and overlap measurements, we simulate the QSL at zero temperature and finite magnetic field, thus moving outside of the symmetric subspace. Next, we implement a protocol for quantum subspace expansion-based measurement of spin-spin correlation functions. Finally, we perform QSE-based simulation of the dynamical structure factor obtained from Green's functions of the finite field Kitaev model. Our results show that quantum simulators offer an insight to quasiparticle properties of strongly correlated magnets and can become a valuable tool for studying material science. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 8 pages, 6 figures

Correlation Functions From Tensor Network Influence Functionals: The Case of the Spin-Boson Model

Authors: Haimi Nguyen, Nathan Ng, Lachlan P. Lindoy, Gunhee Park, Andrew J. Millis, Garnet Kin-Lic Chan, David R. Reichman

Abstract: We investigate the application of matrix product state (MPS) representations of the influence functionals (IF) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We… ▽ More We investigate the application of matrix product state (MPS) representations of the influence functionals (IF) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We examine three different IF approaches: one based on the Kadanoff-Baym contour targeting correlation functions at all times, one based on a complex contour targeting the correlation function at a single time, and a steady state formulation which avoids imaginary or complex times, while providing access to correlation functions at all times. We show that within the IF language, the steady state formulation provides a powerful approach to evaluate equilibrium correlation functions. △ Less

Submitted 22 June, 2024; originally announced June 2024.

Investigating the Collective Nature of Cavity Modified Chemical Kinetics under Vibrational Strong Coupling

Authors: Lachlan P. Lindoy, Arkajit Mandal, David R. Reichman

Abstract: In this paper we develop quantum dynamical methods capable of treating the dynamics of chemically reacting systems in an optical cavity in the vibrationally strong-coupling (VSC) limit at finite temperatures and in the presence of a dissipative solvent in both the few and many molecule limits. In the context of two simple models we demonstrate how reactivity in the {\em collective} VSC regime does… ▽ More In this paper we develop quantum dynamical methods capable of treating the dynamics of chemically reacting systems in an optical cavity in the vibrationally strong-coupling (VSC) limit at finite temperatures and in the presence of a dissipative solvent in both the few and many molecule limits. In the context of two simple models we demonstrate how reactivity in the {\em collective} VSC regime does not exhibit altered rate behavior in equilibrium, but may exhibit resonant cavity modification of reactivity when the system is explicitly out of equilibrium. Our results suggest experimental protocols that may be used to modify reactivity in the collective regime and point to features not included in the models studied which demand further scrutiny. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Unified Framework for Open Quantum Dynamics with Memory

Authors: Felix Ivander, Lachlan P. Lindoy, Joonho Lee

Abstract: Studies of the dynamics of a quantum system coupled to baths are typically performed by utilizing the Nakajima-Zwanzig memory kernel (${\mathcal{K}}$) or the influence functions ($\mathbf{I}$), especially when the dynamics exhibit memory effects (i.e., non-Markovian). Despite their significance, the formal connection between the memory kernel and the influence functions has not been explicitly mad… ▽ More Studies of the dynamics of a quantum system coupled to baths are typically performed by utilizing the Nakajima-Zwanzig memory kernel (${\mathcal{K}}$) or the influence functions ($\mathbf{I}$), especially when the dynamics exhibit memory effects (i.e., non-Markovian). Despite their significance, the formal connection between the memory kernel and the influence functions has not been explicitly made. We reveal their relation by inspecting the system propagator for a broad class of problems where an $N$-level system is linearly coupled to Gaussian baths (bosonic, fermionic, and spin.) With this connection, we also show how approximate path integral methods can be understood in terms of approximate memory kernels. For a certain class of open quantum system problems, we devised a non-perturbative, diagrammatic approach to construct ${\mathcal{K}}$ from $\mathbf{I}$ for (driven) systems interacting with Gaussian baths without the use of any projection-free dynamics inputs required by standard approaches. Lastly, we demonstrate a Hamiltonian learning procedure to extract the bath spectral density from a set of reduced system trajectories obtained experimentally or by numerically exact methods, opening new avenues in quantum sensing and engineering. The insights we provide in this work will significantly advance the understanding of non-Markovian dynamics, and they will be an important stepping stone for theoretical and experimental developments in this area. △ Less

Submitted 3 June, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Modelling non-Markovian noise in driven superconducting qubits

Authors: Abhishek Agarwal, Lachlan P. Lindoy, Deep Lall, Francois Jamet, Ivan Rungger

Abstract: Non-Markovian noise can be a significant source of errors in superconducting qubits. We develop gate sequences utilising mirrored pseudoidentities that allow us to characterise and model the effects of non-Markovian noise on both idle and driven qubits. We compare three approaches to modelling the observed noise: (i) a Markovian noise model, (ii) a model including interactions with a two-level sys… ▽ More Non-Markovian noise can be a significant source of errors in superconducting qubits. We develop gate sequences utilising mirrored pseudoidentities that allow us to characterise and model the effects of non-Markovian noise on both idle and driven qubits. We compare three approaches to modelling the observed noise: (i) a Markovian noise model, (ii) a model including interactions with a two-level system (TLS), (iii) a model utilising the post Markovian master equation (PMME), which we show to be equivalent to the qubit-TLS model in certain regimes. When running our noise characterisation circuits on a superconducting qubit device we find that purely Markovian noise models cannot reproduce the experimental data. Our model based on a qubit-TLS interaction, on the other hand, is able to closely capture the observed experimental behaviour for both idle and driven qubits. We investigate the stability of the noise properties of the hardware over time, and find that the parameter governing the qubit-TLS interaction strength fluctuates significantly even over short time-scales of a few minutes. Finally, we evaluate the changes in the noise parameters when increasing the qubit drive pulse amplitude. We find that although the hardware noise parameters fluctuate significantly over different days, their drive pulse induced relative variation is rather well defined within computed uncertainties: both the phase error and the qubit-TLS interaction strength change significantly with the pulse strength, with the phase error changing quadratically with the amplitude of the applied pulse. Since our noise model can closely describe the behaviour of idle and driven qubits, it is ideally suited to be used in the development of quantum error mitigation and correction methods. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Anderson impurity solver integrating tensor network methods with quantum computing

Authors: Francois Jamet, Lachlan P. Lindoy, Yannic Rath, Connor Lenihan, Abhishek Agarwal, Enrico Fontana, Fedor Simkovic IV, Baptiste Anselme Martin, Ivan Rungger

Abstract: Solving the Anderson impurity model typically involves a two-step process, where one first calculates the ground state of the Hamiltonian, and then computes its dynamical properties to obtain the Green's function. Here we propose a hybrid classical/quantum algorithm where the first step is performed using a classical computer to obtain the tensor network ground state as well as its quantum circuit… ▽ More Solving the Anderson impurity model typically involves a two-step process, where one first calculates the ground state of the Hamiltonian, and then computes its dynamical properties to obtain the Green's function. Here we propose a hybrid classical/quantum algorithm where the first step is performed using a classical computer to obtain the tensor network ground state as well as its quantum circuit representation, and the second step is executed on the quantum computer to obtain the Green's function. Our algorithm exploits the efficiency of tensor networks for preparing ground states on classical computers, and takes advantage of quantum processors for the evaluation of the time evolution, which can become intractable on classical computers. We demonstrate the algorithm using 24 qubits on a quantum computing emulator for SrVO$_3$ with a multi-orbital Anderson impurity model within the dynamical mean field theory. The tensor network based ground state quantum circuit preparation algorithm can also be performed for up to 60 qubits with our available computing resources, while the state vector emulation of the quantum algorithm for time evolution is beyond what is accessible with such resources. We show that, provided the tensor network calculation is able to accurately obtain the ground state energy, this scheme does not require a perfect reproduction of the ground state wave function on the quantum circuit to give an accurate Green's function. This hybrid approach may lead to quantum advantage in materials simulations where the ground state can be computed classically, but where the dynamical properties cannot. △ Less

Submitted 21 October, 2024; v1 submitted 13 April, 2023; originally announced April 2023.

Quantum Dynamics of Vibrational Polariton Chemistry

Authors: Lachlan P. Lindoy, Arkajit Mandal, David R. Reichman

Abstract: We employ an exact quantum mechanical simulation technique to investigate a model of cavity-modified chemical reactions in the condensed phase. The model contains the coupling of the reaction coordinate to a generic solvent, cavity coupling to either the reaction coordinate or a non-reactive mode, and the coupling of the cavity to lossy modes. Thus, many of the most important features needed for r… ▽ More We employ an exact quantum mechanical simulation technique to investigate a model of cavity-modified chemical reactions in the condensed phase. The model contains the coupling of the reaction coordinate to a generic solvent, cavity coupling to either the reaction coordinate or a non-reactive mode, and the coupling of the cavity to lossy modes. Thus, many of the most important features needed for realistic modeling of the cavity modification of chemical reactions are included. We find that when a molecule is coupled to an optical cavity it is essential to treat the problem quantum mechanically in order to obtain a quantitative account of alterations to reactivity. We find sizable and sharp changes in the rate constant that are associated with quantum mechanical state splittings and resonances. The features that emerge from our simulations are closer to those observed in experiments than are previous calculations, even for realistically small values of coupling and cavity loss. This work highlights the importance of a fully quantum treatment of vibrational polariton chemistry. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Two-Dimensional Spectroscopy of Two-Dimensional Materials

Authors: Lachlan P Lindoy, Yao-Wen Chang, David R Reichman

Abstract: In this work we provide an exact and efficient numerical approach to simulate multi-time correlation functions in the Mahan-Nozières-De Dominicis model, which crudely mimics the spectral properties of doped two-dimensional semiconductors such as monolayer transition metal dichalcogenides. We apply this approach to study the coherent two-dimensional electronic spectra of the model. We show that sev… ▽ More In this work we provide an exact and efficient numerical approach to simulate multi-time correlation functions in the Mahan-Nozières-De Dominicis model, which crudely mimics the spectral properties of doped two-dimensional semiconductors such as monolayer transition metal dichalcogenides. We apply this approach to study the coherent two-dimensional electronic spectra of the model. We show that several experimentally observed phenomena, such as peak asymmetry and coherent oscillations in the waiting-time dependence of the trion-exciton cross peaks of the two-dimensional rephasing spectrum, emerge naturally in our approach. Additional features are also present which find no correspondence with experimentally expected behavior. We trace these features to the infinite hole mass property of the model. We use this understanding to construct an efficient approach which filters out configurations associated with the lack of exciton recoil, enabling the connection to previous work and providing a route to the construction of realistic two-dimensional spectra over a broad doping range in two-dimensional semiconductors. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 8 Pages, 6 Figures

Resonant Cavity Modification of Ground State Chemical Kinetics

Authors: Lachlan P. Lindoy, Arkajit Mandal, David R. Reichman

Abstract: Recent experiments have suggested that ground state chemical kinetics can be suppressed or enhanced by coupling the vibrational degrees of freedom of a molecular system with a radiation mode inside an optical cavity. Experiments show that the chemical rate is strongly modified when the photon frequency is close to characteristic vibrational frequencies. The origin of this remarkable effect remains… ▽ More Recent experiments have suggested that ground state chemical kinetics can be suppressed or enhanced by coupling the vibrational degrees of freedom of a molecular system with a radiation mode inside an optical cavity. Experiments show that the chemical rate is strongly modified when the photon frequency is close to characteristic vibrational frequencies. The origin of this remarkable effect remains unknown. In this work, we develop an analytical rate theory for cavity-modified ground state chemical kinetics based on the Pollak-Grabert-Hänggi rate theory. Unlike previous work, our theory covers the complete range of solvent friction values, from the energy-diffusion limited to the spatial-diffusion limited regimes. We show that the chemical reaction rate can either be enhanced or suppressed depending on the bath friction; when bath friction is weak chemical kinetics is enhanced as opposed to the case of strong bath friction, where chemical kinetics is suppressed. Further, we show that the photon frequency at which maximum modification of chemical rate is achieved is close to the reactant well, and hence resonant rate modification occurs. In the strong friction limit the {\it resonant} photon frequency is instead close to the barrier frequency, as obtained using the Grote-Hynes rate theory. Finally, we observe that the rate changes (as a function of photon frequency) are much sharper and more sizable in the weak friction limit than in the strong friction limit, and become increasingly sharp with decreasing well frequency. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Time Evolution of ML-MCTDH Wavefunctions II: Application of the Projector Splitting Integrator

Authors: Lachlan P. Lindoy, Benedikt Kloss, David R. Reichman

Abstract: The multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) approach suffers from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of a regularization parameter. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MC… ▽ More The multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) approach suffers from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of a regularization parameter. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to $10^6$ bath modes, and find that for these problems the PSI requires roughly 3-4 orders of magnitude fewer Hamiltonian evaluations and 2-3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH, and 2-3/1-2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to $\sim 1.3\times 10^9$ parameters to obtain accurate dynamics. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: Main paper: 17 pages, 7 figures. Supplemental information (for paper I and paper II of the series): 22 pages, 2 figures

Time Evolution of ML-MCTDH Wavefunctions I: Gauge Conditions, Basis Functions, and Singularities

Authors: Lachlan P. Lindoy, Benedikt Kloss, David R. Reichman

Abstract: We derive a family of equations of motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than that used for standard ML-MCTDH. These alternative conditions re… ▽ More We derive a family of equations of motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than that used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOMs approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 10 + $ε$ pages, 2 figures

Spin relaxation in radical pairs from the stochastic Schrödinger equation

Authors: Thomas P. Fay, Lachlan P. Lindoy, David E. Manolopoulos

Abstract: We show that the stochastic Schrödinger equation (SSE) provides an ideal way to simulate the quantum mechanical spin dynamics of radical pairs. Electron spin relaxation effects arising from fluctuations in the spin Hamiltonian are straightforward to include in this approach, and their treatment can be combined with a highly efficient stochastic evaluation of the trace over nuclear spin states that… ▽ More We show that the stochastic Schrödinger equation (SSE) provides an ideal way to simulate the quantum mechanical spin dynamics of radical pairs. Electron spin relaxation effects arising from fluctuations in the spin Hamiltonian are straightforward to include in this approach, and their treatment can be combined with a highly efficient stochastic evaluation of the trace over nuclear spin states that is required to compute experimental observables. These features are illustrated in example applications to a flavin-tryptophan radical pair of interest in avian magnetoreception, and to a problem involving spin-selective radical pair recombination along a molecular wire. In the first of these examples, the SSE is shown to be both more efficient and more widely applicable than a recent stochastic implementation of the Lindblad equation, which only provides a valid treatment of relaxation in the extreme-narrowing limit. In the second, the exact SSE results are used to assess the accuracy of a recently-proposed combination of Nakajima-Zwanzig theory for the spin relaxation and Schulten-Wolynes theory for the spin dynamics, which is applicable to radical pairs with many more nuclear spins. An appendix analyses the efficiency of trace sampling in some detail, highlighting the particular advantages of sampling with SU(N) coherent states. △ Less

Submitted 26 February, 2021; originally announced February 2021.

Comments: 14 pages, 4 figures

Journal ref: J. Chem. Phys. 154, 084121 (2021)

Quantum mechanical spin dynamics of a molecular magnetoreceptor

Authors: Lachlan P. Lindoy, Thomas P. Fay, David E. Manolopoulos

Abstract: Radical pair recombination reactions are known to be sensitive to extremely weak magnetic fields, and can therefore be said to function as molecular magnetoreceptors. The classic example is a carotenoid-porphyrin-fullerene (C+PF-) radical pair that has been shown to provide a "proof-of-principle" for the operation of a chemical compass [K. Maeda et al., Nature 453, 387 (2008)]. Previous simulation… ▽ More Radical pair recombination reactions are known to be sensitive to extremely weak magnetic fields, and can therefore be said to function as molecular magnetoreceptors. The classic example is a carotenoid-porphyrin-fullerene (C+PF-) radical pair that has been shown to provide a "proof-of-principle" for the operation of a chemical compass [K. Maeda et al., Nature 453, 387 (2008)]. Previous simulations of this radical pair have employed semiclassical approximations, which are routinely applicable to its 47 coupled electronic and nuclear spins. However, calculating the exact quantum mechanical spin dynamics presents a significant challenge, and has not been possible before now. Here we use a recently developed method to perform numerically converged simulations of the C+PF- quantum mechanical spin dynamics, including all coupled spins. Comparison of these quantum mechanical simulations with various semiclassical approximations reveals that, while it is not perfect, the best semiclassical approximation does capture essentially all of the relevant physics in this problem. △ Less

Submitted 23 April, 2020; originally announced April 2020.

Comments: 8 pages, 4 figures

Journal ref: J. Chem. Phys. 152, 164107 (2020)

Electron spin relaxation in radical pairs: Beyond the Redfield approximation

Authors: Thomas P. Fay, Lachlan P. Lindoy, David E. Manolopoulos

Abstract: Relaxation processes can have a large effect on the spin selective electron transfer reactions of radical pairs. These processes are often treated using phenomenological relaxation superoperators or with some model for the microscopic relaxation mechanism treated within Bloch-Redfield-Wangsness theory. Here, we demonstrate that an alternative perturbative relaxation theory, based on the Nakajima-Z… ▽ More Relaxation processes can have a large effect on the spin selective electron transfer reactions of radical pairs. These processes are often treated using phenomenological relaxation superoperators or with some model for the microscopic relaxation mechanism treated within Bloch-Redfield-Wangsness theory. Here, we demonstrate that an alternative perturbative relaxation theory, based on the Nakajima-Zwanzig equation, has certain advantages over Redfield theory. In particular, the Nakajima-Zwanzig equation does not suffer from the severe positivity problem of Redfield theory in the static disorder limit. Combining the Nakajima-Zwanzig approach consistently with the Schulten-Wolynes semiclassical method, we obtain an efficient method for modeling the spin dynamics of radical pairs containing many hyperfine-coupled nuclear spins. This is then used to investigate the spin-dependent electron transfer reactions and intersystem crossing of dimethyljulolidine-naphthalenediimide (DMJ-NDI) radical ion pairs. By comparing our simulations with experimental data, we find evidence for a field-independent contribution to the triplet quantum yields of these reactions which cannot be explained by electron spin relaxation alone. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Journal ref: J. Chem. Phys. 151, 154117 (2019)

On the calculation of quantum mechanical electron transfer rates

Authors: Joseph E. Lawrence, Theo Fletcher, Lachlan P. Lindoy, David E. Manolopoulos

Abstract: We present a simple interpolation formula for the rate of an electron transfer reaction as a function of the electronic coupling strength. The formula only requires the calculation of Fermi Golden Rule and Born-Oppenheimer rates and so can be combined with any methods that are able to calculate these rates. We first demonstrate the accuracy of the formula by applying it to a one dimensional scatte… ▽ More We present a simple interpolation formula for the rate of an electron transfer reaction as a function of the electronic coupling strength. The formula only requires the calculation of Fermi Golden Rule and Born-Oppenheimer rates and so can be combined with any methods that are able to calculate these rates. We first demonstrate the accuracy of the formula by applying it to a one dimensional scattering problem for which the exact quantum mechanical, Fermi Golden Rule, and Born-Oppenheimer rates are readily calculated. We then describe how the formula can be combined with the Wolynes theory approximation to the Golden Rule rate, and the ring polymer molecular dynamics (RPMD) approximation to the Born-Oppenheimer rate, and used to capture the effects of nuclear tunnelling, zero point energy, and solvent friction on condensed phase electron transfer reactions. Comparison with exact hierarchical equations of motion (HEOM) results for a demanding set of spin-boson models shows that the interpolation formula has an error comparable to that of RPMD rate theory in the adiabatic limit, and that of Wolynes theory in non-adiabatic limit, and is therefore as accurate as any method could possibly be that attempts to generalise these methods to arbitrary electronic coupling strengths. △ Less

Submitted 21 September, 2019; originally announced September 2019.

Comments: 17 pages, 5 figures

Journal ref: J. Chem. Phys. 151, 114119 (2019)

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

Authors: Thomas P. Fay, Lachlan P. Lindoy, David E. Manolopoulos

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… ▽ More 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. △ Less

Submitted 9 August, 2018; originally announced August 2018.

Simple and accurate method for central spin problems

Authors: Lachlan P. Lindoy, David E. Manolopoulos

Abstract: We describe a simple quantum mechanical method that can be used to obtain accurate numerical results over long time scales 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 fro… ▽ More We describe a simple quantum mechanical method that can be used to obtain accurate numerical results over long time scales 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 (t-DMRG) approach. As a result, it can be applied to larger central spin problems than the algebraic Bethe ansatz, and for longer times than the t-DMRG 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. △ Less

Submitted 10 August, 2018; v1 submitted 5 June, 2018; originally announced June 2018.

Comments: Article: 5 pages, 5 figures Supplementary material: 15 pages, 6 figures

Journal ref: Phys. Rev. Lett. 120, 220604 (2018)