BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jeremy O. Richardson

Jeremy O. Richardson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phquant-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2024arXiv

Understanding the Cavity Born-Oppenheimer Approximation

Experiments have demonstrated that vibrational strong coupling between molecular vibrations and light modes can significantly change molecular properties, such as ground-state reactivity. Theoretical studies towards the origin of this exciting observation can roughly be divided in two categories, with studies based on Hamiltonians that simply couple a molecule to a cavity mode via its ground-state dipole moment on the one hand, and on the other hand ab initio calculations that self-consistently include the effect of the cavity mode on the electronic ground state within the cavity Born-Oppenheimer (CBO) approximation; these approaches are not equivalent. The CBO approach is more rigorous, but unfortunately it requires the rewriting of electronic-structure code, and gives little physical insight. In this work, we exploit the relation between the two approaches and demonstrate on a real molecule (hydrogen fluoride) that for realistic coupling strengths, we can recover CBO energies and spectra to high accuracy using only out-of-cavity quantities from standard electronic-structure calculations. In doing so, we discover what the physical effects underlying the CBO results are. Our methodology can aid in incorporating more, possibly important features in models, play a pivotal role in demystifying CBO results and provide a practical and efficient alternative to full CBO calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonadiabatic instanton rate theory beyond the golden-rule limit

Fermi's golden rule describes the leading-order behaviour of the reaction rate as a function of the diabatic coupling. Its asymptotic $(\hbar \rightarrow 0)$ limit is the semiclassical golden-rule instanton rate theory, which rigorously approximates nuclear quantum effects, lends itself to efficient numerical computation and gives physical insight into reaction mechanisms. However the golden rule by itself becomes insufficient as the strength of the diabatic coupling increases, so higher-order terms must be additionally considered. In this work we give a first-principles derivation of the next-order term beyond the golden rule, represented as a sum of three components. Two of them lead to new instanton pathways that extend the golden-rule case and, among other factors, account for the effects of recrossing on the full rate. The remaining component derives from the equilibrium partition function and accounts for changes in potential energy around the reactant and product wells due to diabatic coupling. The new semiclassical theory demands little computational effort beyond a golden-rule instanton calculation. It makes it possible to rigorously assess the accuracy of the golden-rule approximation and sets the stage for future work on general semiclassical nonadiabatic rate theories.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum tunnelling driven H$_2$ formation on graphene

It is commonly believed that it is unfavourable for adsorbed H atoms on carbonaceous surfaces to form H$_2$ without the help of incident H atoms. Using ring-polymer instanton theory to describe multidimensional tunnelling effects, combined with ab initio electronic structure calculations, we find that these quantum-mechanical simulations reveal a qualitatively different picture. Recombination of adsorbed H atoms, which was believed to be irrelevant at low temperature due to high barriers, is enabled by deep tunnelling, with reaction rates enhanced by tens of orders of magnitude. Furthermore, we identify a new path for H recombination that proceeds via multidimensional tunnelling, but would have been predicted to be unfeasible by a simple one-dimensional description of the reaction. The results suggest that hydrogen molecule formation at low temperatures are rather fast processes that should not be ignored in experimental settings and natural environments with graphene, graphite and other planar carbon segments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Transfer learning for affordable and high quality tunneling splittings from instanton calculations

The combination of transfer learning (TL) a low level potential energy surface (PES) to a higher level of electronic structure theory together with ring-polymer instanton (RPI) theory is explored and applied to malonaldehyde. The RPI approach provides a semiclassical approximation of the tunneling splitting and depends sensitively on the accuracy of the PES. With second order Møller-Plesset perturbation theory (MP2) as the low-level (LL) model and energies and forces from coupled cluster singles, doubles and perturbative triples (CCSD(T)) as the high-level (HL) model, it is demonstrated that CCSD(T) information from only 25 to 50 judiciously selected structures along and around the instanton path suffice to reach HL-accuracy for the tunneling splitting. In addition, the global quality of the HL-PES is demonstrated through a mean average error of 0.3 kcal/mol for energies up to 40 kcal/mol above the minimum energy structure (a factor of 2 higher than the energies employed during TL) and $< 2 $ cm$^{-1}$ for harmonic frequencies compared with computationally challenging normal mode calculations at the CCSD(T) level.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Instanton theory of tunnelling in molecules with asymmetric isotopic substitutions

We consider quantum tunnelling in asymmetric double-well systems for which the local minima in the two wells have the same energy, but the frequencies differ slightly. We derive a generalization of instanton theory for these asymmetric systems, leading to a semiclassical expression for the tunnelling matrix element and hence the energy level splitting. We benchmark the method using a set of one- and two-dimensional models, for which the results compare favourably with numerically exact quantum calculations. Using the ring-polymer instanton approach, we apply the method to compute the level splittings in various isotopic-substituted versions of malonaldehyde in full dimensionality and analyse the relative contributions from the zero-point energy difference and tunnelling effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revisiting nuclear tunnelling in the aqueous ferrous-ferric electron transfer

The aqueous ferrous-ferric system provides a classic example of an electron-transfer process in solution. There has been a long standing argument spanning more than three decades around the importance of nuclear tunnelling in this system, with estimates based on Wolynes theory suggesting a quantum correction factor of 65, while estimates based on a related spin-boson model suggest a smaller factor of 7-36. Recently, we have shown that Wolynes theory can break down for systems with multiple transition states leading to an overestimation of the rate, and we suggest that a liquid system such as the one investigated here may be particularly prone to this. We re-investigate this old yet interesting system with the first application of the recently developed golden-rule quantum transition-state theory (GR-QTST). We find that GR-QTST can be applied to this complex system without apparent difficulties and that it gives a prediction for the quantum rate 6 times smaller than that from Wolynes theory. The fact that these theories give different results suggests that although it is well known that the system can be treated using linear response and therefore resembles a spin-boson model in the classical limit, this approximation is questionable in the quantum case. It also intriguingly suggests the possibility that the previous predictions were overestimating the rate due to a break down of Wolynes theory.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Generalized spin mapping for quantum-classical dynamics

We recently derived a spin-mapping approach for treating the nonadiabatic dynamics of a two-level system in a classical environment [J. Chem. Phys. 151, 044119 (2019)] based on the well-known quantum equivalence between a two-level system and a spin-1/2 particle. In the present paper, we generalize this method to describe the dynamics of $N$-level systems. This is done via a mapping to a classical phase space that preserves the $SU(N)$-symmetry of the original quantum problem. The theory reproduces the standard Meyer--Miller--Stock--Thoss Hamiltonian without invoking an extended phase space, and we thus avoid leakage from the physical subspace. In contrast with the standard derivation of this Hamiltonian, the generalized spin mapping leads to an $N$-dependent value of the zero-point energy parameter that is uniquely determined by the Casimir invariant of the $N$-level system. Based on this mapping, we derive a simple way to approximate correlation functions in complex nonadiabatic molecular systems via classical trajectories, and present benchmark calculations on the seven-state Fenna--Matthews--Olson complex. The results are significantly more accurate than conventional Ehrenfest dynamics, at a comparable computational cost, and can compete in accuracy with other state-of-the-art mapping approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Instanton formulation of Fermi&#39;s golden rule in the Marcus inverted regime

Fermi&#39;s golden rule defines the transition rate between weakly coupled states and can thus be used to describe a multitude of molecular processes including electron-transfer reactions and light-matter interaction. However, it can only be calculated if the wave functions of all internal states are known, which is typically not the case in molecular systems. Marcus theory provides a closed-form expression for the rate constant, which is a classical limit of the golden rule, and indicates the existence of a normal regime and an inverted regime. Semiclassical instanton theory presents a more accurate approximation to the golden-rule rate including nuclear quantum effects such as tunnelling, which has so far been applicable to complex anharmonic systems in the normal regime only. In this paper we extend the instanton method to the inverted regime and study the properties of the periodic orbit, which describes the tunnelling mechanism via two imaginary-time trajectories, one of which now travels in negative imaginary time. It is known that tunnelling is particularly prevalent in the inverted regime, even at room temperature, and thus this method is expected to be useful in studying a wide range of molecular transitions occurring in this regime.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nonadiabatic quantum transition-state theory in the golden-rule limit. II. Overcoming the pitfalls of the saddle-point and semiclassical approximations

We describe a path-integral molecular dynamics implementation of our recently developed golden-rule quantum transition-state theory (GR-QTST). The method is applied to compute the reaction rate in various models of electron transfer and benchmarked against exact results. We demonstrate that for systems exhibiting two or more transition states, rates computed using Wolynes theory [P. G. Wolynes, J.\ Chem.\ Phys.\ 87, 6559 (1987)] can be overestimated by orders of magnitude, whereas the GR-QTST predictions are numerically accurate. This is the case both at low temperature, where nuclear tunneling makes a considerable contribution, and also in the classical limit, where only GR-QTST rigorously tends to the correct result. Analysis shows that the saddle-point approximation employed by Wolynes theory is not valid in this case, which results in predictions of unphysical reaction pathways, whilst the energy constraint employed by GR-QTST resolves this problem. The GR-QTST method is also seen to give accurate results for a strongly anharmonic system by sampling configurations around the instanton pathway without making the semiclassical approximation. These promising results indicate that the GR-QTST method could be an efficient and accurate approach for simulating electron-transfer reactions in complex molecular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Spin-mapping approach for nonadiabatic molecular dynamics

We propose a trajectory-based method for simulating nonadiabatic dynamics in molecular systems with two coupled electronic states. Employing a quantum-mechanically exact mapping of the two-level problem to a spin-1/2 coherent state, we construct a classical phase space of a spin vector constrained to a spherical surface with a radius consistent with the quantum magnitude of the spin. In contrast with the singly-excited harmonic oscillator basis used in Meyer-Miller-Stock-Thoss (MMST) mapping, the theory requires no additional projection operators onto the space of physical states. When treated under a quasiclassical approximation, we show that the resulting dynamics is equivalent to that generated by the MMST Hamiltonian. What differs is the value of the zero-point energy parameter as well as the initial distribution and the measurement operators. For various spin-boson models the results of our method are seen to be a significant improvement compared to both standard Ehrenfest dynamics and linearized semiclassical MMST mapping, without adding any computational complexity.

0 citations0 reviewsNo code linkedOpen access