BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David E. Manolopoulos

David E. Manolopoulos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phcond-mat.mtrl-sciphysics.comp-phcond-mat.stat-mechcond-mat.softquant-ph

Catalog footprint

What is connected

16works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A general non-adiabatic quantum instanton approximation

We present a general quantum instanton approach to calculating reaction rates for systems with two electronic states and arbitrary values of the electronic coupling. This new approach, which we call the non-adiabatic quantum instanton (NAQI) approximation, reduces to Wolynes theory in the golden rule limit and to a recently proposed projected quantum instanton (PQI) method in the adiabatic limit. As in both of these earlier theories, the NAQI approach is based on making a saddle point approximation to the time integral of a reactive flux autocorrelation function, although with a generalised definition of the projection operator onto the product states. We illustrate the accuracy of the approach by comparison with exact rates for one dimensional scattering problems and discuss its applicability to more complex reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum mechanical spin dynamics of a molecular magnetoreceptor

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An efficient quantum mechanical method for radical pair recombination reactions

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(Z^2\log 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(MZ\log Z)$, where $M$ is the number of Monte Carlo samples needed for convergence. Example calculations on a strongly-coupled radical pair with $Z>10^6$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Floquet theory of radical pairs in radiofrequency magnetic fields

We present a new method for calculating the product yield of a radical pair recombination reaction in the presence of a weak time-dependent magnetic field. This method successfully circumvents the computational difficulties presented by a direct solution of the Liouville-von Neumann equation for a long-lived radical pair containing many hyperfine-coupled nuclear spins. Using a modified formulation of Floquet theory, treating the time-dependent magnetic field as a perturbation, and exploiting the slow radical pair recombination, we show that one can obtain a good approximation to the product yield by considering only nearly-degenerate sub-spaces of the Floquet space. Within a significant parameter range, the resulting method is found to give product yields in good agreement with exact quantum mechanical results for a variety of simple model radical pairs. Moreover it is considerably more efficient than the exact calculation, and it can be applied to radical pairs containing significantly more nuclear spins. This promises to open the door to realistic theoretical investigations of the effect of radiofrequency electromagnetic radiation on the photochemically induced radical pair recombination reactions in the avian retina which are believed to be responsible for the magnetic compass sense of migratory birds.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nuclear quantum effects in water exchange around lithium and fluoride ions

We employ classical and ring polymer molecular dynamics simulations to study the effect of nuclear quantum fluctuations on the structure and the water exchange dynamics of aqueous solutions of lithium and fluoride ions. While we obtain reasonably good agreement with experimental data for solutions of lithium by augmenting the Coulombic interactions between the ion and the water molecules with a standard Lennard-Jones ion-oxygen potential, the same is not true for solutions of fluoride, for which we find that a potential with a softer repulsive wall gives much better agreement. A small degree of destabilization of the first hydration shell is found in quantum simulations of both ions when compared with classical simulations, with the shell becoming less sharply defined and the mean residence time of the water molecules in the shell decreasing. In line with these modest differences, we find that the mechanisms of the exchange processes are unaffected by quantization, so a classical description of these reactions gives qualitatively correct and quantitatively reasonable results. We also find that the quantum effects in solutions of lithium are larger than in solutions of fluoride. This is partly due to the stronger interaction of lithium with water molecules, partly due to the lighter mass of lithium, and partly due to competing quantum effects in the hydration of fluoride, which are absent in the hydration of lithium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An improved semiclassical theory of radical pair recombination reactions

We present a practical semiclassical method for computing the electron spin dynamics of a radical in which the electron spin is hyperfine coupled to a large number of nuclear spins. This can be used to calculate the singlet and triplet survival probabilities and quantum yields of radical recombination reactions in the presence of magnetic fields. Our method differs from the early semiclassical theory of Schulten and Wolynes [J. Chem. Phys. 68, 3292 (1978)] in allowing each individual nuclear spin to precess around the electron spin, rather than assuming that the hyperfine coupling-weighted sum of nuclear spin vectors is fixed in space. The downside of removing this assumption is that one can no longer obtain a simple closed-form expression for the electron spin correlation tensor: our method requires a numerical calculation. However, the computational effort increases only linearly with the number of nuclear spins, rather than exponentially as in an exact quantum mechanical calculation. The method is therefore applicable to arbitrarily large radicals. Moreover, it approaches quantitative agreement with quantum mechanics as the number of nuclear spins increases and the environment of the electron spin becomes more complex, owing to the rapid quantum decoherence in complex systems. Unlike the Schulten-Wolynes theory, the present semiclassical theory predicts the correct long-time behaviour of the electron spin correlation tensor, and it therefore correctly captures the low magnetic field effect in the singlet yield of a radical recombination reaction with a slow recombination rate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Asymmetric recombination and electron spin relaxation in the semiclassical theory of radical pair reactions

We describe how the semiclassical theory of radical pair recombination reactions recently introduced by two of us [D. E. Manolopoulos and P. J. Hore, J. Chem. Phys. 139, 124106 (2013)] can be generalised to allow for different singlet and triplet recombination rates. This is a non-trivial generalisation because when the recombination rates are different the recombination process is dynamically coupled to the coherent electron spin dynamics of the radical pair. Furthermore, because the recombination operator is a two-electron operator, it is no longer sufficient simply to consider the two electrons as classical vectors: one has to consider the complete set of 16 two-electron spin operators as independent classical variables. The resulting semiclassical theory is first validated by comparison with exact quantum mechanical results for a model radical pair containing 12 nuclear spins. It is then used to shed light on the spin dynamics of a carotenoid-porphyrin-fullerene (CPF) triad containing considerably more nuclear spins which has recently been used to establish a 'proof of principle' for the operation of a chemical compass [K. Maeda et al., Nature 453, 387 (2008)]. We find in particular that the intriguing biphasic behaviour that has been observed in the effect of an Earth-strength magnetic field on the time-dependent survival probability of the photo-excited C+PF- radical pair arises from a delicate balance between its asymmetric recombination and the relaxation of the electron spin in the carotenoid radical.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei

Light nuclei at room temperature and below exhibit a kinetic energy which significantly deviates from the predictions of classical statistical mechanics. This quantum kinetic energy is responsible for a wide variety of isotope effects of interest in fields ranging from chemistry to climatology. It also furnishes the second moment of the nuclear momentum distribution, which contains subtle information about the chemical environment and has recently become accessible to deep inelastic neutron scattering experiments. Here we show how, by combining imaginary time path integral dynamics with a carefully designed generalized Langevin equation, it is possible to dramatically reduce the expense of computing the quantum kinetic energy. We also introduce a transient anisotropic Gaussian approximation to the nuclear momentum distribution which can be calculated with negligible additional effort. As an example, we evaluate the structural properties, the quantum kinetic energy, and the nuclear momentum distribution for a first-principles simulation of liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How to remove the spurious resonances from ring polymer molecular dynamics

Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the desirable features of RPMD: thermostatted RPMD (TRPMD) is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

i-PI: A Python interface for ab initio path integral molecular dynamics simulations

Recent developments in path integral methodology have significantly reduced the computational expense of including quantum mechanical effects in the nuclear motion in ab initio molecular dynamics simulations. However, the implementation of these developments requires a considerable programming effort, which has hindered their adoption. Here we describe i-PI, an interface written in Python that has been designed to minimise the effort required to bring state-of-the-art path integral techniques to an electronic structure program. While it is best suited to first principles calculations and path integral molecular dynamics, i-PI can also be used to perform classical molecular dynamics simulations, and can just as easily be interfaced with an empirical forcefield code. To give just one example of the many potential applications of the interface, we use it in conjunction with the CP2K electronic structure package to showcase the importance of nuclear quantum effects in high pressure water.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct Measurement of Competing Quantum Effects on the Kinetic Energy of Heavy Water upon Melting

Even at room temperature, quantum mechanics plays a major role in determining the quantitative behaviour of light nuclei, changing significantly the values of physical properties such as the heat capacity. However, other observables appear to be only weakly affected by nuclear quantum effects (NQEs): for instance, the melting temperatures of light and heavy water differ by less than 4 K. Recent theoretical work has attributed this to a competition between intra and inter molecular NQEs, which can be separated by computing the anisotropy of the quantum kinetic energy tensor. The principal values of this tensor change in opposite directions when ice melts, leading to a very small net quantum mechanical effect on the melting point. This paper presents the first direct experimental observation of this phenomenon, achieved by measuring the deuterium momentum distributions n(p) in heavy water and ice using Deep Inelastic Neutron Scattering (DINS), and resolving their anisotropy. Results from the experiments, supplemented by a theoretical analysis, show that the anisotropy of the quantum kinetic energy tensor can also be captured for heavier atoms such as oxygen.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accelerating the convergence of path integral dynamics with a generalized Langevin equation

The quantum nature of nuclei plays an important role in the accurate modelling of light atoms such as hydrogen, but it is often neglected in simulations due to the high computational overhead involved. It has recently been shown that zero-point energy effects can be included comparatively cheaply in simulations of harmonic and quasi-harmonic systems by augmenting classical molecular dynamics with a generalized Langevin equation (GLE). Here we describe how a similar approach can be used to accelerate the convergence of path integral (PI) molecular dynamics to the exact quantum mechanical result in more strongly anharmonic systems exhibiting both zero point energy and tunnelling effects. The resulting PI-GLE method is illustrated with applications to a double-well tunnelling problem and to liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The inefficiency of re-weighted sampling and the curse of system size in high order path integration

Computing averages over a target probability density by statistical re-weighting of a set of samples with a different distribution is a strategy which is commonly adopted in fields as diverse as atomistic simulation and finance. Here we present a very general analysis of the accuracy and efficiency of this approach, highlighting some of its weaknesses. We then give an example of how our results can be used, specifically to assess the feasibility of high-order path integral methods. We demonstrate that the most promising of these techniques -- which is based on re-weighted sampling -- is bound to fail as the size of the system is increased, because of the exponential growth of the statistical uncertainty in the re-weighted average.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Competing quantum effects in the dynamics of a flexible water model

Numerous studies have identified large quantum mechanical effects in the dynamics of liquid water. In this paper, we suggest that these effects may have been overestimated due to the use of rigid water models and flexible models in which the intramolecular interactions were described using simple harmonic functions. To demonstrate this, we introduce a new simple point charge model for liquid water, q-TIP4P/F, in which the O--H stretches are described by Morse-type functions. We have parameterized this model to give the correct liquid structure, diffusion coefficient, and infra-red absorption frequencies in quantum (path integral-based) simulations. By comparing classical and quantum simulations of the liquid, we find that quantum mechanical fluctuations increase the rates of translational diffusion and orientational relaxation in our model by a factor of around 1.15. This effect is much smaller than that observed in all previous simulations of simple empirical water models, which have found a quantum effect of at least 1.4 regardless of the quantum simulation method or the water model employed. The small quantum effect in our model is a result of two competing phenomena. Intermolecular zero point energy and tunneling effects destabilize the hydrogen bonding network, leading to a less viscous liquid with a larger diffusion coefficient. However this is offset by intramolecular zero point motion, which changes the average water monomer geometry resulting in a larger dipole moment, stronger intermolecular interactions, and slower diffusion. We end by suggesting, on the basis of simulations of other potential energy models, that the small quantum effect we find in the diffusion coefficient is associated with the ability of our model to produce a single broad O-H stretching band in the infra-red absorption spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient stochastic thermostatting of path integral molecular dynamics

The path integral molecular dynamics (PIMD) method provides a convenient way to compute the quantum mechanical structural and thermodynamic properties of condensed phase systems at the expense of introducing an additional set of high-frequency normal modes on top of the physical vibrations of the system. Efficiently sampling such a wide range of frequencies provides a considerable thermostatting challenge. Here we introduce a simple stochastic path integral Langevin equation (PILE) thermostat which exploits an analytic knowledge of the free path integral normal mode frequencies. We also apply a recently-developed colored-noise thermostat based on a generalized Langevin equation (GLE), which automatically achieves a similar, frequency-optimized sampling. The sampling efficiencies of these thermostats are compared with that of the more conventional Nosé-Hoover chain (NHC) thermostat for a number of physically relevant properties of the liquid water and hydrogen-in-palladium systems. In nearly every case, the new PILE thermostat is found to perform just as well as the NHC thermostat while allowing for a computationally more efficient implementation. The GLE thermostat also proves to be very robust delivering a near-optimum sampling efficiency in all of the cases considered. We suspect that these simple stochastic thermostats will therefore find useful application in many future PIMD simulations.

0 citations0 reviewsNo code linkedOpen access