A three-dimensional model of binary chemical reactions is studied. We consider an ab initio quantum two-particle system subjected to an attractive interaction potential and to a heat bath at thermal equilibrium at absolute temperature T > 0 . Under the sole action of the attraction potential, the two particles can either be bound or unbound to each other. While at T = 0 , there is no transition between both states, such a transition is possible when T > 0 (due to the heat bath) and plays a key role as k B T approaches the magnitude of the attractive potential. We focus on a quantum regime, typical of chemical reactions, such that: (a) the thermal wavelength is shorter than the range of the attractive potential (lower limit on T) and (b) ( 3 / 2 ) k B T does not exceed the magnitude of the attractive potential (upper limit on T). In this regime, we extend several methods previously applied to analyze the time duration of DNA thermal denaturation. The two-particle system is then described by a non-equilibrium Wigner function. Under Assumptions (a) and (b), and for sufficiently long times, defined by a characteristic time scale D that is subsequently estimated, the general dissipationless non-equilibrium equation for the Wigner function is approximated by a Smoluchowski-like equation displaying dissipation and quantum effects. A comparison with the standard chemical kinetic equations is made. The time τ required for the two particles to transition from the bound state to unbound configurations is studied by means of the mean first passage time formalism. An approximate formula for τ, in terms of D and exhibiting the Arrhenius exponential factor, is obtained. Recombination processes are also briefly studied within our framework and compared with previous well-known methods.
A three-dimensional model of binary chemical reactions is studied. We consider an ab initio quantum two-particle system subjected to an attractive interaction potential and to a heat bath at thermal equilibrium at absolute temperature T > 0 . Under the sole action of the attraction potential, the two particles can either be bound or unbound to each other. While at T = 0 , there is no transition between both states, such a transition is possible when T > 0 (due to the heat bath) and plays a key role as k B T approaches the magnitude of the attractive potential. We focus on a quantum regime, typical of chemical reactions, such that: (a) the thermal wavelength is shorter than the range of the attractive potential (lower limit on T) and (b) ( 3 / 2 ) k B T does not exceed the magnitude of the attractive potential (upper limit on T). In this regime, we extend several methods previously applied to analyze the time duration of DNA thermal denaturation. The two-particle system is then described by a non-equilibrium Wigner function. Under Assumptions (a) and (b), and for sufficiently long times, defined by a characteristic time scale D that is subsequently estimated, the general dissipationless non-equilibrium equation for the Wigner function is approximated by a Smoluchowski-like equation displaying dissipation and quantum effects. A comparison with the standard chemical kinetic equations is made. The time τ required for the two particles to transition from the bound state to unbound configurations is studied by means of the mean first passage time formalism. An approximate formula for τ, in terms of D and exhibiting the Arrhenius exponential factor, is obtained. Recombination processes are also briefly studied within our framework and compared with previous well-known methods.
We treat the non-equilibrium evolution of an open one-particle statistical system, subject to a potential and to an external “heat bath” (hb) with negligible dissipation. For the classical equilibrium Boltzmann distribution, Wc,eq, a non-equilibrium three-term hierarchy for moments fulfills Hermiticity, which allows one to justify an approximate long-time thermalization. That gives partial dynamical support to Boltzmann’s Wc,eq, out of the set of classical stationary distributions, Wc;st, also investigated here, for which neither Hermiticity nor that thermalization hold, in general. For closed classical many-particle systems without hb (by using Wc,eq), the long-time approximate thermalization for three-term hierarchies is justified and yields an approximate Lyapunov function and an arrow of time. The largest part of the work treats an open quantum one-particle system through the non-equilibrium Wigner function, W. Weq for a repulsive finite square well is reported. W’s (< 0 in various cases) are assumed to be quasi-definite functionals regarding their dependences on momentum (q). That yields orthogonal polynomials, HQ,n(q), for Weq (and for stationary Wst), non-equilibrium moments, Wn, of W and hierarchies. For the first excited state of the harmonic oscillator, its stationary Wst is a quasi-definite functional, and the orthogonal polynomials and three-term hierarchy are studied. In general, the non-equilibrium quantum hierarchies (associated with Weq) for the Wn’s are not three-term ones. As an illustration, we outline a non-equilibrium four-term hierarchy and its solution in terms of generalized operator continued fractions. Such structures also allow one to formulate long-time approximations, but make it more difficult to justify thermalization. For large thermal and de Broglie wavelengths, the dominant Weq and a non-equilibrium equation for W are reported: the non-equilibrium hierarchy could plausibly be a three-term one and possibly not far from Gaussian, and thermalization could possibly be justified.
We investigate the irreversible evolution of a small system in which a chemical reaction takes place.We have two main goals: the first requires to find an equation to produce a time-irreversible behavior,the second consists in introducing a simple exactly solvable model in order to understand basic facts in chemical kinetics.Our basic tool is the transition function counting the number of paths joining two points in the reactive coordinates system.An exact quantum Smoluchowski equation is derived for the reactive system in vacuum, in the presence of a solvent in equilibrium at any time with the reactive system a new Smoluchowski equation is obtained.The transition from a quantum regime to a classical one is discussed.The case of a reactive system not in equilibrium with its neighborhood is investigated in terms of path integral and via a partial differential function.Memory effects and closure assumptions are discussed.Using a simple potential model, the chemical rate constant is exactly calculated and questions such as the meaning of the activation energy or the physical content of the so-called prefactor are investigated.
We consider non-equilibrium open statistical systems, subject to potentials and to external “heat baths” (hb) at thermal equilibrium at temperature T (either with ab initio dissipation or without it). Boltzmann’s classical equilibrium distributions generate, as Gaussian weight functions in momenta, orthogonal polynomials in momenta (the position-independent Hermite polynomialsHn’s). The moments of non-equilibrium classical distributions, implied by the Hn’s, fulfill a hierarchy: for long times, the lowest moment dominates the evolution towards thermal equilibrium, either with dissipation or without it (but under certain approximation). We revisit that hierarchy, whose solution depends on operator continued fractions. We review our generalization of that moment method to classical closed many-particle interacting systems with neither a hb nor ab initio dissipation: with initial states describing thermal equilibrium at T at large distances but non-equilibrium at finite distances, the moment method yields, approximately, irreversible thermalization of the whole system at T, for long times. Generalizations to non-equilibrium quantum interacting systems meet additional difficulties. Three of them are: (i) equilibrium distributions (represented through Wigner functions) are neither Gaussian in momenta nor known in closed form; (ii) they may depend on dissipation; and (iii) the orthogonal polynomials in momenta generated by them depend also on positions. We generalize the moment method, dealing with (i), (ii) and (iii), to some non-equilibrium one-particle quantum interacting systems. Open problems are discussed briefly.
Non-equilibrium evolution at absolute temperature T and approach to equilibrium of statistical systems in long-time (t) approximations, using both hierarchies and functional integrals, are reviewed. A classical non-relativistic particle in one spatial dimension, subject to a potential and a heat bath (hb), is described by the non-equilibrium reversible Liouville distribution (W) and equation, with a suitable initial condition. The Boltzmann equilibrium distribution Weq generates orthogonal (Hermite) polynomials Hn in momenta. Suitable moments Wn of W (using the Hn’s) yield a non-equilibrium three-term hierarchy (different from the standard Bogoliubov–Born–Green–Kirkwood–Yvon one), solved through operator continued fractions. After a long-t approximation, the Wn’s yield irreversibly approach to equilibrium. The approach is extended (without hb) to: (i) a non-equilibrium system of N classical non-relativistic particles interacting through repulsive short range potentials and (ii) a classical ϕ4 field theory (without hb). The extension to one non-relativistic quantum particle (with hb) employs the non-equilibrium Wigner function (WQ): difficulties related to non-positivity of WQ are bypassed so as to formulate approximately approach to equilibrium. A non-equilibrium quantum anharmonic oscillator is analyzed differently, through functional integral methods. The latter allows an extension to relativistic quantum ϕ4 field theory (a meson gas off-equilibrium, without hb), facing ultraviolet divergences and renormalization. Genuine simplifications of quantum ϕ4 theory at high T and large distances and long t occur; then, through a new argument for the field-theoretic case, the theory can be approximated by a classical ϕ4 one, yielding an approach to equilibrium.
Recent progress in our understanding of quantum effects on the Brownian motion in an external potential is reviewed. This problem is ubiquitous in physics and chemistry, particularly in the context of decay of metastable states, for example, the reversal of the magnetization of a single domain ferromagnetic particle, kinetics of a superconducting tunnelling junction, etc. Emphasis is laid on the establishment of master equations describing the diffusion process in phase space analogous to the classical Fokker-Planck equation. In particular, it is shown how Wigner's [E. P. Wigner, Phys. Rev., 1932, 40, 749] method of obtaining quantum corrections to the classical equilibrium Maxwell-Boltzmann distribution may be extended to the dissipative non-equilibrium dynamics governing the quantum Brownian motion in an external potential V(x), yielding a master equation for the Wigner distribution function W(x,p,t) in phase space (x,p). The explicit form of the master equation so obtained contains quantum correction terms up to o(h(4)) and in the classical limit, h --> 0, reduces to the classical Klein-Kramers equation. For a quantum oscillator, the method yields an evolution equation coinciding in all respects with that of Agarwal [G. S. Agarwal, Phys. Rev. A, 1971, 4, 739]. In the high dissipation limit, the master equation reduces to a semi-classical Smoluchowski equation describing non-inertial quantum diffusion in configuration space. The Wigner function formulation of quantum Brownian motion is further illustrated by finding quantum corrections to the Kramers escape rate, which, in appropriate limits, reduce to those yielded via quantum generalizations of reaction rate theory.
The ultimate regularity of quantum mechanics creates a tension with the assumption of classical chaos used in many of our pictures of chemical reaction dynamics. Out-of-time-order correlators (OTOCs) provide a quantum analog to the Lyapunov exponents that characterize classical chaotic motion. Maldacena, Shenker, and Stanford have suggested a fundamental quantum bound for the rate of information scrambling, which resembles a limit suggested by Herzfeld for chemical reaction rates. Here, we use OTOCs to study model reactions based on a double-well reaction coordinate coupled to anharmonic oscillators or to a continuum oscillator bath. Upon cooling, as one enters the tunneling regime where the reaction rate does not strongly depend on temperature, the quantum Lyapunov exponent can approach the scrambling bound and the effective reaction rate obtained from a population correlation function can approach the Herzfeld limit on reaction rates: Tunneling increases scrambling by expanding the state space available to the system. The coupling of a dissipative continuum bath to the reaction coordinate reduces the scrambling rate obtained from the early-time OTOC, thus making the scrambling bound harder to reach, in the same way that friction is known to lower the temperature at which thermally activated barrier crossing goes over to the low-temperature activationless tunneling regime. Thus, chemical reactions entering the tunneling regime can be information scramblers as powerful as the black holes to which the quantum Lyapunov exponent bound has usually been applied.