Identification of the invariant manifolds of the LiCN molecule using Lagrangian descriptors.

In this paper, we apply Lagrangian descriptors to study the invariant manifolds that emerge from the top of two barriers existing in the LiCN⇌LiNC isomerization reaction. We demonstrate that the integration times must be large enough compared with the characteristic stability exponents of the periodic orbit under study. The invariant manifolds manifest as singularities in the Lagrangian descriptors. Furthermore, we develop an equivalent potential energy surface with 2 degrees of freedom, which reproduces with a great accuracy previous results [F. Revuelta, R. M. Benito, and F. Borondo, Phys. Rev. E 99, 032221 (2019)2470-004510.1103/PhysRevE.99.032221]. This surface allows the use of an adiabatic approximation to develop a more simplified potential energy with solely 1 degree of freedom. The reduced dimensional model is still able to qualitatively describe the results observed with the original 2-degrees-of-freedom potential energy landscape. Likewise, it is also used to study in a more simple manner the influence on the Lagrangian descriptors of a bifurcation, where some of the previous invariant manifolds emerge, even before it takes place.

Bohr-Sommerfeld-like quantization in the theory of walking droplets.

Recent experiments have shown that self-propelled millimetric walking droplets bouncing on a vibrating liquid surface exhibit phenomena, such as interference or tunneling, that so far were thought to be possible only in the microscopic realm. Here we present calculations showing that the surface wave satisfies, in the long-memory limit, a Bohr-Sommerfeld quantization-like relation. This strongly suggest the possibility of a novel fundamental type of quantization in these experiments, which can simultaneously explain their emulation of the quantum behavior and, more importantly, shed light into some of the interpretational difficulties of the standard quantum theory.

Short-periodic-orbit method for excited chaotic eigenfunctions.

An alternative method for the calculation of excited chaotic eigenfunctions in arbitrary energy windows is presented. We demonstrate the feasibility of using wave functions localized on unstable periodic orbits as efficient basis sets for this task in classically chaotic systems. The number of required localized wave functions is only of the order of the ratio t_{H}/t_{E}, with t_{H} the Heisenberg time and t_{E} the Ehrenfest time. As an illustration, we present convincing results for a coupled two-dimensional quartic oscillator with chaotic dynamics.

Finite-barrier corrections for multidimensional barriers in colored noise.

The usual identification of reactive trajectories for the calculation of reaction rates requires very time-consuming simulations, particularly if the environment presents memory effects. In this paper, we develop a method that permits the identification of reactive trajectories in a system under the action of a stochastic colored driving. This method is based on the perturbative computation of the invariant structures that act as separatrices for reactivity. Furthermore, using this perturbative scheme, we have obtained a formally exact expression for the reaction rate in multidimensional systems coupled to colored noisy environments.

Unveiling the chaotic structure in phase space of molecular systems using Lagrangian descriptors.

We explore here the feasibility of using the recently introduced Lagrangian descriptors [A. M. Mancho et al., Commun. Nonlinear Sci. Numer. Simul. 18, 3530 (2013)1007-570410.1016/j.cnsns.2013.05.002] to unveil the usually rich dynamics taking place in the vibrations of molecular systems, especially if they are floppy. The principal novelty of our work is the inclusion of p norms in the definition of the descriptors in this kind of system, which greatly enhances their power to discern among the different structures existing in the phase space. As an illustration we use the LiCN molecule described by realistic potentials in two and three dimensions, which exhibits chaotic motion within a mixed phase space in the isomerization between the two wells corresponding to the linear isomer stable configurations, LiNC and LiCN. In particular, we pay special attention to the manifolds emerging from the unstable fixed point between the corresponding isomer wells, and also to the marginally stable structures around a parabolic point existing near the LiNC well.

Quantization Scheme for the Experiments with "Walking Droplets".

In this paper we explore the interest and feasibility of quantizing the macroscopic surface wave generated in the dynamics of walking droplets on a vertically vibrated liquid surface in the limit of high memory of the droplet trajectory, where an astonishing similarity with the quantum behavior has been experimentally observed.

Transition state theory for activated systems with driven anharmonic barriers.

Classical transition state theory has been extended to address chemical reactions across barriers that are driven and anharmonic. This resolves a challenge to the naive theory that necessarily leads to recrossings and approximate rates because it relies on a fixed dividing surface. We develop both perturbative and numerical methods for the computation of a time-dependent recrossing-free dividing surface for a model anharmonic system in a solvated environment that interacts strongly with an oscillatory external field. We extend our previous work, which relied either on a harmonic approximation or on periodic force driving. We demonstrate that the reaction rate, expressed as the long-time flux of reactive trajectories, can be extracted directly from the stability exponents, namely, Lyapunov exponents, of the moving dividing surface. Comparison to numerical results demonstrates the accuracy and robustness of this approach for the computation of optimal (recrossing-free) dividing surfaces and reaction rates in systems with Markovian solvation forces. The resulting reaction rates are in strong agreement with those determined from the long-time flux of reactive trajectories.

Semiclassical basis sets for the computation of molecular vibrational states.

In this paper, we extend a method recently reported [F. Revuelta et al., Phys. Rev. E 87, 042921 (2013)] for the calculation of the eigenstates of classically highly chaotic systems to cases of mixed dynamics, i.e., those presenting regular and irregular motions at the same energy. The efficiency of the method, which is based on the use of a semiclassical basis set of localized wave functions, is demonstrated by applying it to the determination of the vibrational states of a realistic molecular system, namely, the LiCN molecule.

Transition state theory for solvated reactions beyond recrossing-free dividing surfaces.

The accuracy of rate constants calculated using transition state theory depends crucially on the correct identification of a recrossing-free dividing surface. We show here that it is possible to define such optimal dividing surface in systems with non-Markovian friction. However, a more direct approach to rate calculation is based on invariant manifolds and avoids the use of a dividing surface altogether, Using that method we obtain an explicit expression for the rate of crossing an anharmonic potential barrier. The excellent performance of our method is illustrated with an application to a realistic model for LiNC⇌LiCN isomerization.

Transition state geometry of driven chemical reactions on time-dependent double-well potentials.

Reaction rates across time-dependent barriers are difficult to define and difficult to obtain using standard transition state theory approaches because of the complexity of the geometry of the dividing surface separating reactants and products. Using perturbation theory (PT) or Lagrangian descriptors (LDs), we can obtain the transition state trajectory and the associated recrossing-free dividing surface. With the latter, we are able to determine the exact reactant population decay and the corresponding rates to benchmark the PT and LD approaches. Specifically, accurate rates are obtained from a local description regarding only direct barrier crossings and to those obtained from a stability analysis of the transition state trajectory. We find that these benchmarks agree with the PT and LD approaches for obtaining recrossing-free dividing surfaces. This result holds not only for the local dynamics in the vicinity of the barrier top, but also for the global dynamics of particles that are quenched at the reactant or product wells after their sojourn over the barrier region. The double-well structure of the potential allows for long-time dynamics related to collisions with the outside walls that lead to long-time returns in the low-friction regime. This additional global dynamics introduces slow-decay pathways that do not result from the local transition across the recrossing-free dividing surface associated with the transition state trajectory, but can be addressed if that structure is augmented by the population transfer of the long-time returns.

Scar Functions, Barriers for Chemical Reactivity, and Vibrational Basis Sets.

The performance of a recently proposed method to efficiently calculate scar functions is analyzed in problems of chemical interest. An application to the computation of wave functions associated with barriers relevant for the LiNC ⇄ LiCN isomerization reaction is presented as an illustration. These scar functions also constitute excellent elements for basis sets suitable for quantum calculation of vibrational energy levels. To illustrate their efficiency, a calculation of the LiNC/LiCN eigenfunctions is also presented.

Using basis sets of scar functions.

We present a method to efficiently compute the eigenfunctions of classically chaotic systems. The key point is the definition of a modified Gram-Schmidt procedure which selects the most suitable elements from a basis set of scar functions localized along the shortest periodic orbits of the system. In this way, one benefits from the semiclassical dynamical properties of such functions. The performance of the method is assessed by presenting an application to a quartic two-dimensional oscillator whose classical dynamics are highly chaotic. We have been able to compute the eigenfunctions of the system using a small basis set. An estimate of the basis size is obtained from the mean participation ratio. A thorough analysis of the results using different indicators, such as eigenstate reconstruction in the local representation, scar intensities, participation ratios, and error bounds, is also presented.

Reaction rate calculation with time-dependent invariant manifolds.

The identification of trajectories that contribute to the reaction rate is the crucial dynamical ingredient in any classical chemical reactivity calculation. This problem often requires a full scale numerical simulation of the dynamics, in particular if the reactive system is exposed to the influence of a heat bath. As an efficient alternative, we propose here to compute invariant surfaces in the phase space of the reactive system that separate reactive from nonreactive trajectories. The location of these invariant manifolds depends both on time and on the realization of the driving force exerted by the bath. These manifolds allow the identification of reactive trajectories simply from their initial conditions, without the need of any further simulation. In this paper, we show how these invariant manifolds can be calculated, and used in a formally exact reaction rate calculation based on perturbation theory for any multidimensional potential coupled to a noisy environment.

Computationally efficient method to construct scar functions.

The performance of a simple method [E. L. Sibert III, E. Vergini, R. M. Benito, and F. Borondo, New J. Phys. 10, 053016 (2008)] to efficiently compute scar functions along unstable periodic orbits with complicated trajectories in configuration space is discussed, using a classically chaotic two-dimensional quartic oscillator as an illustration.

Communication: transition state theory for dissipative systems without a dividing surface.

Transition state theory is a central cornerstone in reaction dynamics. Its key step is the identification of a dividing surface that is crossed only once by all reactive trajectories. This assumption is often badly violated, especially when the reactive system is coupled to an environment. The calculations made in this way then overestimate the reaction rate and the results depend critically on the choice of the dividing surface. In this Communication, we study the phase space of a stochastically driven system close to an energetic barrier in order to identify the geometric structure unambiguously determining the reactive trajectories, which is then incorporated in a simple rate formula for reactions in condensed phase that is both independent of the dividing surface and exact.