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

Timothy J. H. Hele

Timothy J. H. Hele contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phcond-mat.stat-mechquant-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

12 published item(s)

preprint2022arXiv

Electronic energies from coupled fermionic 'Zombie' states imaginary time evolution

Zombie States are a recently introduced formalism to describe coupled coherent Fermionic states which address the Fermionic sign problem in a computationally tractable manner. Previously it has been shown that Zombie States with fractional occupations of spin-orbitals obeyed the correct Fermionic creation and annihilation algebra and presented results for real-time evolution [Dmitrii V. Shalashilin, J. Chem. Phys. 148, 194109 (2018)]. In this work we extend and build on this formalism by developing efficient algorithms for evaluating the Hamiltonian and other operators between Zombie States and address their normalization. We also show how imaginary time propagation can be used to find the ground state of a system. We also present a biasing method, for setting up a basis set of random Zombie States, that allow much smaller basis sizes to be used while still accurately describing the electronic structure Hamiltonian and its ground state and describe a technique of wave function "cleaning" which removes the contributions of configurations with the wrong number of electrons, improving the accuracy further. We also show how low-lying excited states can be calculated efficiently using a Gram-Schmidt orthogonalization procedure.The proposed algorithm of imaginary time propagation on a biased random grids of Zombie States may present an alternative to existing Quantum Monte Carlo methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Inverse molecular design from first principles: tailoring organic chromophore spectra for optoelectronic applications

The discovery of molecules with tailored optoelectronic properties such as specific frequency and intensity of absorption or emission is a major challenge in creating next-generation organic light-emitting diodes (OLEDs) and photovoltaics. This raises the question: how can we predict a potential chemical structure from these properties? Approaches that attempt to tackle this inverse design problem include virtual screening, active machine learning and genetic algorithms. However, these approaches rely on a molecular database or many electronic structure calculations, and significant computational savings could be achieved if there was prior knowledge of (i) whether the optoelectronic properties of a parent molecule could easily be improved and (ii) what morphing operations on a parent molecule could improve these properties. In this perspective we address both of these challenges from first principles. We firstly adapt the Thomas-Reiche-Kuhn sum rule to organic chromophores and show how this indicates how easily the absorption and emission of a molecule can be improved. We then show how by combining electronic structure theory and intensity borrowing perturbation theory we can predict whether or not the proposed morphing operations will achieve the desired spectral alteration, and thereby derive widely-applicable design rules. We go on to provide proof-of-concept illustrations of this approach to optimizing the visible absorption of acenes and the emission of radical OLEDs. We believe this approach can be integrated into genetic algorithms by biasing morphing operations in favour of those which are likely to be successful, leading to faster molecular discovery and greener chemistry.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An alternative derivation of ring-polymer molecular dynamics transition-state theory

In a previous article [J. Chem. Phys. 138, 084108 (2013)], we showed that the $t\to 0_+$ limit of ring-polymer molecular dynamics (RPMD) rate-theory is also the $t\to 0_+$ limit of a new type of quantum flux-side time-correlation function, in which the dividing surfaces are invariant to imaginary-time translation; in other words, that RPMD transition-state theory (RPMD-TST) is a $t\to 0_+$ quantum transition-state theory (QTST). Recently, Jang and Voth [J. Chem. Phys. 144, 084110 (2016)] rederived this quantum $t\to 0_+$ limit, and claimed that it gives instead the centroid-density approximation. Here we show that the $t\to 0_+$ limit derived by Jang and Voth is in fact RPMD-TST.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or classical Wigner approximation) results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e. a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads $N \to \infty$, such that the lowest normal-mode frequencies take their Matsubara values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of $\hbar^2$ at $\hbar^0$ (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting Matsubara dynamics is inherently classical (since all terms $\mathcal{O}\left(\hbar^{2}\right)$ disappear from the Matsubara Liouvillian in the limit $N \to \infty$), and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes, and gives better agreement than LSC- IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

How to obtain thermostatted ring polymer molecular dynamics from exact quantum dynamics and when to use it

We obtain thermostatted ring polymer molecular dynamics (TRPMD) from exact quantum dynamics via Matsubara dynamics, a recently-derived form of linearization which conserves the quantum Boltzmann distribution. Performing a contour integral in the complex quantum Boltzmann distribution of Matsubara dynamics, replacement of the imaginary Liouvillian which results with a Fokker-Planck term gives TRPMD. We thereby provide error terms between TRPMD and quantum dynamics and predict the systems in which they are likely to be small. Using a harmonic analysis we show that careful addition of friction causes the correct oscillation frequency of the higher ring-polymer normal modes in a harmonic well, which we illustrate with calculation of the position-squared autocorrelation function. However, no physical friction parameter will produce the correct fluctuation dynamics for a parabolic barrier. The results in this paper are consistent with previous numerical studies and advise the use of TRPMD for the computation of spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relation of centroid molecular dynamics and ring-polymer molecular dynamics to exact quantum dynamics

We recently obtained a quantum-Boltzmann-conserving classical dynamics by making a single change to the derivation of the `Classical Wigner' approximation. Here, we show that the further approximation of this `Matsubara dynamics' gives rise to two popular heuristic methods for treating quantum Boltzmann time-correlation functions: centroid molecular dynamics (CMD) and ring-polymer molecular dynamics (RPMD). We show that CMD is a mean-field approximation to Matsubara dynamics, obtained by discarding (classical) fluctuations around the centroid, and that RPMD is the result of discarding a term in the Matsubara Liouvillian which shifts the frequencies of these fluctuations. These findings are consistent with previous numerical results, and give explicit formulae for the terms that CMD and RPMD leave out.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Should Thermostatted Ring Polymer Molecular Dynamics be used to calculate thermal reaction rates?

We apply Thermostatted Ring Polymer Molecular Dynamics (TRPMD), a recently-proposed approximate quantum dynamics method, to the computation of thermal reaction rates. Its short-time Transition-State Theory (TST) limit is identical to rigorous Quantum Transition-State Theory, and we find that its long-time limit is independent of the location of the dividing surface. TRPMD rate theory is then applied to one-dimensional model systems, the atom-diatom bimolecular reactions H+H$_2$, D+MuH and F+H$_2$, and the prototypical polyatomic reaction H+CH$_4$. Above the crossover temperature, the TRPMD rate is virtually invariant to the strength of the friction applied to the internal ring-polymer normal modes, and beneath the crossover temperature the TRPMD rate generally decreases with increasing friction, in agreement with the predictions of Kramers theory. We therefore find that TRPMD is approximately equal to, or less accurate than, Ring Polymer Molecular Dynamics (RPMD) for symmetric reactions, and for certain asymmetric systems and friction parameters closer to the quantum result, providing a basis for further assessment of the accuracy of this method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum Transition-State Theory

This dissertation unifies one of the central methods of classical rate calculation, `Transition-State Theory' (TST), with quantum mechanics, thereby deriving a rigorous `Quantum Transition-State Theory' (QTST). The resulting QTST is identical to ring polymer molecular dynamics transition-state theory (RPMD-TST), which was previously considered a heuristic method, and whose results we thereby validate. The key step in deriving a QTST is alignment of the flux and side dividing surfaces in path-integral space to obtain a quantum flux-side time-correlation function with a non-zero $t\to 0_+$ limit. We then prove that this produces the exact quantum rate in the absence of recrossing by the exact quantum dynamics, fulfilling the requirements of a QTST. Furthermore, strong evidence is presented that this is the only QTST with positive-definite Boltzmann statistics and therefore the pre-eminent method for computation of thermal quantum rates in direct reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Derivation of a true (t->0+) quantum transition-state theory. II. Recovery of the exact quantum rate in the absence of recrossing

In Part I [J. Chem. Phys. 138, 084108 (2013)] we derived a quantum transition-state theory by taking the t->0+ (short-time) limit of a new form of quantum flux-side time-correlation function containing a ring-polymer dividing surface. This t->0+ limit appears to be unique in giving positive-definite Boltzmann statistics, and is identical to ring-polymer molecular dynamics (RPMD) TST. Here, we show that quantum TST (i.e. RPMD-TST) is exact if there is no recrossing (by the real-time quantum dynamics) of the ring-polymer dividing surface, nor of any surface orthogonal to it in the space describing fluctuations in the polymer-bead positions along the reaction coordinate. In practice, this means that RPMD-TST gives a good approximation to the exact quantum rate for direct reactions, provided the temperature is not too far below the cross-over to deep tunnelling. We derive these results by comparing the long-time limit of the ring-polymer flux-side time-correlation function with that of a hybrid flux-side time-correlation function (containing a ring-polymer flux operator and a Miller-Schwarz-Tromp side function), and by representing the resulting ring-polymer momentum integrals as hypercubes. Together with Part I, the results of this article validate a large number of RPMD calculations of reaction rates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Derivation of a true (t->0+) quantum transition-state theory. Part I: Uniqueness and equivalence to ring-polymer molecular dynamics transition-state-theory

Surprisingly, there exists a quantum flux-side time-correlation function which has a non-zero short-time (t->0+) limit, and thus yields a rigorous quantum generalization of classical transition-state theory (TST). In this Part I of two articles, we introduce the new time-correlation function, and derive its short-time limit. The new ingredient is a generalized Kubo transform which allows the flux and side dividing surfaces to be the same function of path-integral space. Choosing this common dividing surface to be a single point gives a short-time limit which is identical to an expression introduced on heuristic grounds by Wigner in 1932, but which does not give positive-definite quantum statistics, causing it to fail while still in the shallow-tunnelling regime. Choosing the dividing surface to be invariant to imaginary-time translation gives, uniquely, a short-time limit that gives the correct positive- definite quantum statistics at all temperatures, and which is identical to ring-polymer molecular dynamics (RPMD) TST. We find that the RPMD-TST rate is not a strict upper bound to the exact quantum rate, but a good approximation to one if real-time coherence effects are small. Part II will show that the RPMD-TST rate is equal to the exact quantum rate in the absence of recrossing.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the uniqueness of t->0+ quantum transition-state theory

It was shown recently that there exists a true quantum transition-state theory (QTST) corresponding to the t->0+ limit of a (new form of) quantum flux-side time-correlation function. Remarkably, this QTST is identical to ring-polymer molecular dynamics (RPMD) TST. Here we provide evidence which suggests very strongly that this QTST (= RPMD-TST) is unique, in the sense that the t->0+ limit of any other flux-side time-correlation function gives either non-positive-definite quantum statistics or zero. We introduce a generalized flux-side time-correlation function which includes all other (known) flux-side time-correlation functions as special limiting cases. We find that the only non-zero t->0+ limit of this function that contains positive-definite quantum statistics is RPMD-TST.

0 citations0 reviewsNo code linkedOpen access