ℏ 4 quantum corrections to semiclassical transmission probabilities.

The combination of vibrational perturbation theory with the replacement of the harmonic oscillator quantization condition along the reaction coordinate with an imaginary action to be used in the uniform semiclassical approximation for the transmission probability has been shown in recent years to be a practical method for obtaining thermal reaction rates. To date, this theory has been developed systematically only up to second order in perturbation theory. Although it gives the correct leading order term in an ℏ2 expansion, its accuracy at lower temperatures, where tunneling becomes important, is not clear. In this paper, we develop the theory to fourth order in the action. This demands developing the quantum perturbation theory up to sixth order. Remarkably, we find that the fourth order theory gives the correct ℏ4 term in the expansion of the exact thermal rate. The relative magnitude of the fourth order correction as compared to the second order term objectively indicates the accuracy of the second order theory. We also extend the previous modified second order theory to the fourth order case, creating an ℏ2 modified potential for this purpose. The resulting theory is tested on the standard examples-symmetric and asymmetric Eckart potentials and a Gaussian potential. The modified fourth order theory is remarkably accurate for the asymmetric Eckart potential.

A personal perspective of the present status and future challenges facing thermal reaction rate theory.

Reaction rate theory has been at the center of physical chemistry for well over one hundred years. The evolution of the theory is not only of historical interest. Reliable and accurate computation of reaction rates remains a challenge to this very day, especially in view of the development of quantum chemistry methods, which predict the relevant force fields. It is still not possible to compute the numerically exact rate on the fly when the system has more than at most a few dozen anharmonic degrees of freedom, so one must consider various approximate methods, not only from the practical point of view of constructing numerical algorithms but also on conceptual and formal levels. In this Perspective, I present some of the recent analytical results concerning leading order terms in an ℏ2m series expansion of the exact rate and their implications on various approximate theories. A second aspect has to do with the crossover temperature between tunneling and thermal activation. Using a uniform semiclassical transmission probability rather than the "primitive" semiclassical theory leads to the conclusion that there is no divergence problem associated with a "crossover temperature." If one defines a semiclassical crossover temperature as the point at which the tunneling energy of the instanton equals the barrier height, then it is a factor of two higher than its previous estimate based on the "primitive" semiclassical approximation. In the low temperature tunneling regime, the uniform semiclassical theory as well as the "primitive" semiclassical theory were based on the classical Euclidean action of a periodic orbit on the inverted potential. The uniform semiclassical theory wrongly predicts that the "half-point," which is the energy at which the transmission probability equals 1/2, for any barrier potential, is always the barrier energy. We describe here how augmenting the Euclidean action with constant terms of order ℏ2 can significantly improve the accuracy of the semiclassical theory and correct this deficiency. This also leads to a deep connection with and improvement of vibrational perturbation theory. The uniform semiclassical theory also enables an extension of the quantum version of Kramers' turnover theory to temperatures below the "crossover temperature." The implications of these recent advances on various approximate methods used to date are discussed at length, leading to the conclusion that reaction rate theory will continue to challenge us both on conceptual and practical levels for years to come.

ℏ2 Corrections to Semiclassical Transmission Coefficients.

The uniform semiclassical expression for the energy-dependent transmission probability through a barrier has been a staple of reaction rate theory for almost 90 years. Yet, when using the classical Euclidean action, the transmission probability is identical to 1/2 when the energy equals the barrier height since the Euclidean action vanishes at this energy. This result is generally incorrect. It also leads to an inaccurate estimate of the leading order term in an ℏ2n expansion of the thermal transmission coefficient. The central result of this paper is that adding an ℏ2 dependent correction to the uniform semiclassical expression, whether as a constant action or as a shift in the energy scale, not only corrects this inaccuracy but also leads to a theory that is more accurate than the previous one for almost any energy. Shifting the energy scale is a generalization of the vibrational perturbation theory 2 (VPT2) and is much more accurate than the "standard" VPT2 theory, especially when the potential is asymmetric. Shifting the action by a constant is a generalization of a result obtained by Yasumori and Fueki (YF) only for the Eckart barrier. The resulting modified VPT2 and YF semiclassical theories are applied to the symmetric and asymmetric Eckart barrier, a Gaussian barrier, and a tanh barrier. The one-dimensional theories are also generalized to many-dimensional systems. Their effect on the thermal instanton theory is discussed.

Multidimensional uniform semiclassical instanton thermal rate theory.

Instanton-based rate theory is a powerful tool that is used to explore tunneling in many-dimensional systems. Yet, it diverges at the so-called "crossover temperature." Using the uniform semiclassical transmission probability of Kemble [Phys. Rev. 48, 549 (1935)], we showed recently that in one dimension, one might derive a uniform semiclassical instanton rate theory, which has no divergence. In this paper, we generalize this uniform theory to many-dimensional systems. The resulting theory uses the same input as in the previous instanton theory, yet does not suffer from the divergence. The application of the uniform theory to dissipative systems is considered and used to revise Wolynes' well-known analytical expression for the rate [P. G. Wolynes, Phys. Rev. Lett. 47, 968 (1981)] so that it does not diverge at the "crossover temperature."

Uniform Semiclassical Instanton Rate Theory.

The instanton expression for the thermal transmission probability through a one-dimensional barrier is derived by using the uniform semiclassical energy-dependent transmission coefficient of Kemble. The resulting theory does not diverge at the "crossover temperature" but changes smoothly. The temperature-dependent energy of the instanton is the same as the barrier height when ℏßω‡ = π and not 2π as in the "standard" instanton theory. The concept of a crossover temperature between tunneling and thermal activation, based on the divergence of the instanton rate, is obsolete. The theory is improved by assuring that at high energy when the energy-dependent transmission coefficient approaches unity the integrand decays exponentially as dictated by the Boltzmann factor and not as a Gaussian. This ensures that at sufficiently high temperatures the uniform theory reduces to the classical. Application to Eckart barriers demonstrates that the uniform theory provides a good estimate of the numerically exact result over the whole temperature range.

Recent Developments in Kramers' Theory of Reaction Rates.

In this short review, we provide an update of recent developments in Kramers' theory of reaction rates. After a brief introduction stressing the importance of this theory initially developed for chemical reactions, we briefly present the main theoretical formalism starting from the generalized Langevin equation and continue by showing the main points of the modern Pollak, Grabert and Hänggi theory. Kramers' theory is then sketched for quantum and classical surface diffusion. As an illustration the surface diffusion of Na atoms on a Cu(110) surface is discussed showing escape rates, jump distributions and diffusion coefficients as a function of reduced friction. Finally, some very recent applications of turnover theory to different fields such as nanoparticle levitation, microcavity polariton dynamics and simulation of reaction in liquids are presented. We end with several open problems and future challenges faced up by Kramers turnover theory.

Nonadiabatic Couplings Can Speed Up Quantum Tunneling Transition Path Times.

Quantum tunneling is known to play an important role in the dynamics of systems with nonadiabatic couplings. However, until recently, the time-domain properties of nonadiabatic scattering have been severely under-explored. Using numerically exact quantum methods, we study the impact that nonadiabatic couplings have on the time it takes to tunnel through a barrier. We find that the Wigner phase time is the appropriate measure to use when determining the tunneling flight time also when considering nonadiabatic systems. The central result of the present study is that in an avoided crossing system in one dimension, the nonadiabatic couplings speed up the tunneling event, relative to the adiabatic case in which all nonadiabatic coupling is ignored. This has implications for both the study of quantum tunneling times and for the field of nonadiabatic scattering and chemistry.

Microscopic origin of diffusive dynamics in the context of transition path time distributions for protein folding and unfolding.

Experimentally measured transition path time distributions are usually analyzed theoretically in terms of a diffusion equation over a free energy barrier. It is though well understood that the free energy profile separating the folded and unfolded states of a protein is characterized as a transition through many stable micro-states which exist between the folded and unfolded states. Why is it then justified to model the transition path dynamics in terms of a diffusion equation, namely the Smoluchowski equation (SE)? In principle, van Kampen has shown that a nearest neighbor Markov chain of thermal jumps between neighboring microstates will lead in a continuum limit to the SE, such that the friction coefficient is proportional to the mean residence time in each micro-state. However, the practical question of how many microstates are needed to justify modeling the transition path dynamics in terms of an SE has not been addressed. This is a central topic of this paper where we compare numerical results for transition paths based on the diffusion equation on the one hand and the nearest neighbor Markov jump model on the other. Comparison of the transition path time distributions shows that one needs at least a few dozen microstates to obtain reasonable agreement between the two approaches. Using the Markov nearest neighbor model one also obtains good agreement with the experimentally measured transition path time distributions for a DNA hairpin and PrP protein. As found previously when using the diffusion equation, the Markov chain model used here also reproduces the experimentally measured long time tail and confirms that the transition path barrier height is ∼3kBT. This study indicates that in the future, when attempting to model experimentally measured transition path time distributions, one should perhaps prefer a nearest neighbor Markov model which is well defined also for rough energy landscapes. Such studies can also shed light on the minimal number of microstates needed to unravel the experimental data.

The Effect of Time Resolution on Apparent Transition Path Times Observed in Single-Molecule Studies of Biomolecules.

Single-molecule experiments have now achieved a time resolution allowing observation of transition paths, the brief trajectory segments where the molecule undergoing an unfolding or folding transition enters the energetically or entropically unfavorable barrier region from the folded/unfolded side and exits to the unfolded/folded side, thereby completing the transition. This resolution, however, is yet insufficient to identify the precise entrance/exit events that mark the beginning and the end of a transition path: the nature of the diffusive dynamics is such that a molecular trajectory will recross the boundary between the barrier region and the folded/unfolded state, multiple times, at a time scale much shorter than that of the typical experimental resolution. Here we use theory and Brownian dynamics simulations to show that, as a result of such recrossings, the apparent transition path times are generally longer than the true ones. We quantify this effect using a simple model where the observed dynamics is a moving average of the true dynamics and discuss experimental implications of our results.

ℏ2 expansion of the transmission probability through a barrier.

Ninety years ago, Wigner derived the leading order expansion term in ℏ2 for the tunneling rate through a symmetric barrier. His derivation included two contributions: one came from the parabolic barrier, but a second term involved the fourth-order derivative of the potential at the barrier top. He left us with a challenge, which is answered in this paper, to derive the same but for an asymmetric barrier. A crucial element of the derivation is obtaining the ℏ2 expansion term for the projection operator, which appears in the flux-side expression for the rate. It is also reassuring that an analytical calculation of semiclassical transition state theory (TST) reproduces the anharmonic corrections to the leading order of ℏ2. The efficacy of the resulting expression is demonstrated for an Eckart barrier, leading to the conclusion that especially when considering heavy atom tunneling, one should use the expansion derived in this paper, rather than the parabolic barrier approximation. The rate expression derived here reveals how the classical TST limit is approached as a function of ℏ and, thus, provides critical insights to understand the validity of popular approximate theories, such as the classical Wigner, centroid molecular dynamics, and ring polymer molecular dynamics methods.

Transition Path Flight Times and Nonadiabatic Electronic Transitions.

Transition path flight times are studied for scattering on two electronic surfaces with a single crossing. These flight times reveal nontrivial quantum effects such as resonance lifetimes and nonclassical passage times and reveal that nonadiabatic effects often increase flight times. The flight times are computed using numerically exact time propagation and compared with results obtained from the Fewest Switches Surface Hopping (FSSH) method. Comparison of the two methods shows that the FSSH method is reliable for transition path times only when the scattering is classically allowed on the relevant adiabatic surfaces. However, where quantum effects such as tunneling and resonances dominate, the FSSH method is not adequate to accurately predict the correct times and transition probabilities. These results highlight limitations in methods which do not account for quantum interference effects, and suggest that measuring flight times is important for obtaining insights from the time-domain into quantum effects in nonadiabatic scattering.

Coherent state representation of thermal correlation functions with applications to rate theory.

A coherent state phase space representation of operators, based on the Husimi distribution, is used to derive an exact expression for the symmetrized version of thermal correlation functions. In addition to the time and temperature independent phase space representation of the two operators whose correlation function is of interest, the integrand includes a non-negative distribution function where only one imaginary time and one real time propagation are needed to compute it. The methodology is exemplified for the flux side correlation function used in rate theory. The coherent state representation necessitates the use of a smeared Gaussian flux operator whose coherent state phase space representation is identical to the classical flux expression. The resulting coherent state expression for the flux side correlation function has a number of advantages as compared to previous formulations. Since only one time propagation is needed, it is much easier to converge it with a semiclassical initial value representation. There is no need for forward-backward approximations, and in principle, the computation may be implemented on the fly. It also provides a route for analytic semiclassical approximations for the thermal rate, as exemplified by a computation of the transmission factor through symmetric and asymmetric Eckart barriers using a thawed Gaussian approximation for both imaginary and real time propagations. As a by-product, this example shows that one may obtain "good" tunneling rates using only above barrier classical trajectories even in the deep tunneling regime.

Perturbation theory of scattering for grazing-incidence fast-atom diffraction.

Recent grazing-incidence, fast atom diffraction (GIFAD) experiments have highlighted the well known observation that the distance between classical rainbow angles depends on the incident energy. The GIFAD experiments imply an incident vertical scattering angle, facilitating an analytic analysis using classical perturbation theory, which leads to the conclusion that the so called "dynamic corrugation" amplitude, as defined by Bocan et al., Phys. Rev. Lett., 2020 125, 096101 is, within first-order perturbation theory, proportional to the tangent of the rainbow angle. Therefore it provides no further information about the interaction than is gleaned from the rainbow angle and its energy dependence. Perhaps more importantly, the resulting analytic theory reveals how the energy dependence of rainbow angles may be inverted into information on the force field governing the interaction of the incident projectile with the surface.

Lower Bounds for Nonrelativistic Atomic Energies.

A recently developed lower bound theory for Coulombic problems (E. Pollak, R. Martinazzo, J. Chem. Theory Comput. 2021, 17, 1535) is further developed and applied to the highly accurate calculation of the ground-state energy of two- (He, Li+, and H-) and three- (Li) electron atoms. The method has been implemented with explicitly correlated many-particle basis sets of Gaussian type, on the basis of the highly accurate (Ritz) upper bounds they can provide with relatively small numbers of functions. The use of explicitly correlated Gaussians is developed further for computing the variances, and the necessary modifications are here discussed. The computed lower bounds are of submilli-Hartree (parts per million relative) precision and for Li represent the best lower bounds ever obtained. Although not yet as accurate as the corresponding (Ritz) upper bounds, the computed bounds are orders of magnitude tighter than those obtained with other lower bound methods, thereby demonstrating that the proposed method is viable for lower bound calculations in quantum chemistry applications. Among several aspects, the optimization of the wave function is shown to play a key role for both the optimal solution of the lower bound problem and the internal check of the theory.

Comparison of an improved self-consistent lower bound theory with Lehmann's method for low-lying eigenvalues.

Ritz eigenvalues only provide upper bounds for the energy levels, while obtaining lower bounds requires at least the calculation of the variances associated with these eigenvalues. The well-known Weinstein and Temple lower bounds based on the eigenvalues and variances converge very slowly and their quality is considerably worse than that of the Ritz upper bounds. Lehmann presented a method that in principle optimizes Temple's lower bounds with significantly improved results. We have recently formulated a Self-Consistent Lower Bound Theory (SCLBT), which improves upon Temple's results. In this paper, we further improve the SCLBT and compare its quality with Lehmann's theory. The Lánczos algorithm for constructing the Hamiltonian matrix simplifies Lehmann's theory and is essential for the SCLBT method. Using two lattice Hamiltonians, we compared the improved SCLBT (iSCLBT) with its previous implementation as well as with Lehmann's lower bound theory. The novel iSCLBT exhibits a significant improvement over the previous version. Both Lehmann's theory and the SCLBT variants provide significantly better lower bounds than those obtained from Weinstein's and Temple's methods. Compared to each other, the Lehmann and iSCLBT theories exhibit similar performance in terms of the quality and convergence of the lower bounds. By increasing the number of states included in the calculations, the lower bounds are tighter and their quality becomes comparable with that of the Ritz upper bounds. Both methods are suitable for providing lower bounds for low-lying excited states as well. Compared to Lehmann's theory, one of the advantages of the iSCLBT method is that it does not necessarily require the Weinstein lower bound for its initial input, but Ritz eigenvalue estimates can also be used. Especially owing to this property the iSCLBT method sometimes exhibits improved convergence compared to that of Lehmann's lower bounds.

The Influence of the Symmetry of Identical Particles on Flight Times.

In this work, our purpose is to show how the symmetry of identical particles can influence the time evolution of free particles in the nonrelativistic and relativistic domains as well as in the scattering by a potential δ-barrier. For this goal, we consider a system of either two distinguishable or indistinguishable (bosons and fermions) particles. Two sets of initial conditions have been studied: different initial locations with the same momenta, and the same locations with different momenta. The flight time distribution of particles arriving at a 'screen' is calculated in each case from the density and flux. Fermions display broader distributions as compared with either distinguishable particles or bosons, leading to earlier and later arrivals for all the cases analyzed here. The symmetry of the wave function seems to speed up or slow down the propagation of particles. Due to the cross terms, certain initial conditions lead to bimodality in the fermionic case. Within the nonrelativistic domain, and when the short-time survival probability is analyzed, if the cross term becomes important, one finds that the decay of the overlap of fermions is faster than for distinguishable particles which in turn is faster than for bosons. These results are of interest in the short time limit since they imply that the well-known quantum Zeno effect would be stronger for bosons than for fermions. Fermions also arrive earlier and later than bosons when they are scattered by a δ-barrier. Although the particle symmetry does affect the mean tunneling flight time, in the limit of narrow in momentum initial Gaussian wave functions, the mean times are not affected by symmetry but tend to the phase time for distinguishable particles.

What can we learn from transition path time distributions for protein folding and unfolding?

Recent advances in experimental measurements of transition path time distributions have raised intriguing theoretical questions. The present interpretation of the experimental data indicates a small value of the fitted transition path barrier height as compared to the barrier height of the unfolded to folded transition. Secondly, as shown in this paper, it is essential to analyse the experimental data using absorbing boundary conditions at the end points used to determine the transition paths. Such an analysis reveals long time tails that have thus far eluded quantitative theoretical interpretation. Is this due to uncertainty in the experimental data or does it call for a rethinking of the theoretical interpretation? A detailed study of the transition path time distribution using a diffusive model leads to the following conclusions. a. The present experimental data is not accurate enough to discern between functional forms of bell shaped free energy barriers. b. Long time tails indicate the possible existence of a "trap" in the transition path region. c. The "trap" may be considered as a well in the free energy surface. d. The long time tail is quite sensitive to the form of the trap so that future measurements of the long time tail as a function of the location of the end points of the transition path may make it possible to not only determine the well depth but also to distinguish between different functional forms for the free energy surface. e. Introduction of a well along the transition path leads to good fits with the experimental data provided that the transition path barrier height is ∼3kBT, substantially higher than the estimates of ∼1kBT based on bell shaped functions. The results presented here negate the need of introducing multi-dimensional effects, free energy barrier asymmetry, sub-diffusive memory kernels or systematic ruggedness to explain the experimentally measured data.

Lower Bounds for Coulombic Systems.

As of the writing of this paper, lower bounds are not a staple of quantum chemistry computations and for good reason. All previous attempts at applying lower bound theory to Coulombic systems led to lower bounds whose quality was inferior to the Ritz upper bounds so that their added value was minimal. Even our recent improvements upon Temple's lower bound theory were limited to Lanczos basis sets and these are not available to atoms and molecules due to the Coulomb singularity. In the present paper, we overcome these problems by deriving a rather simple eigenvalue equation whose roots, under appropriate conditions, give lower bounds which are competitive with the Ritz upper bounds. The input for the theory is the Ritz eigenvalues and their variances; there is no need to compute the full matrix of the squared Hamiltonian. Along the way, we present a Cauchy-Schwartz inequality which underlies many aspects of lower bound theory. We also show that within the matrix Hamiltonian theory used here, the methods of Lehmann and our recent self-consistent lower bound theory (J. Chem. Phys. 2020, 115, 244110) are identical. Examples include implementation to the hydrogen and helium atoms.

Self-consistent theory of lower bounds for eigenvalues.

A rigorous practically applicable theory is presented for obtaining lower bounds to eigenvalues of Hermitian operators, whether the ground state or excited states. Algorithms are presented for computing "residual energies" whose magnitude is essential for the computation of the eigenvalues. Their practical application is possible due to the usage of the Lanczos method for creating a tridiagonal representation of the operator under study. The theory is self-consistent, in the sense that a lower bound for one state may be used to improve the lower bounds for others, and this is then used self-consistently until convergence. The theory is exemplified for a toy model of a quartic oscillator, where with only five states the relative error in the lower bound for the ground state is reduced to 6 · 10-6, which is the same as the relative error of the least upper bound obtained with the same basis functions. The lower bound method presented in this paper suggests that lower bounds may become a staple of eigenvalue computations.

Lower bounds to eigenvalues of the Schrödinger equation by solution of a 90-y challenge.

The Ritz upper bound to eigenvalues of Hermitian operators is essential for many applications in science. It is a staple of quantum chemistry and physics computations. The lower bound devised by Temple in 1928 [G. Temple, Proc. R. Soc. A Math. Phys. Eng. Sci. 119, 276-293 (1928)] is not, since it converges too slowly. The need for a good lower-bound theorem and algorithm cannot be overstated, since an upper bound alone is not sufficient for determining differences between eigenvalues such as tunneling splittings and spectral features. In this paper, after 90 y, we derive a generalization and improvement of Temple's lower bound. Numerical examples based on implementation of the Lanczos tridiagonalization are provided for nontrivial lattice model Hamiltonians, exemplifying convergence over a range of 13 orders of magnitude. This lower bound is typically at least one order of magnitude better than Temple's result. Its rate of convergence is comparable to that of the Ritz upper bound. It is not limited to ground states. These results complement Ritz's upper bound and may turn the computation of lower bounds into a staple of eigenvalue and spectral problems in physics and chemistry.