Source author record

Rigoberto Hernandez

Rigoberto Hernandez appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

11 published item(s)

preprint2022arXiv

Mean first-passage times for solvated LiCN isomerization at intermediate to high temperatures

The behavior of a particle in a solvent has been framed using stochastic dynamics since the early theory of Kramers. A particle in a chemical reaction reacts slower in a diluted solvent because of the lack of energy transfer via collisions. The flux-over-population reaction rate constant rises with increasing density before falling again for very dense solvents. This Kramers turnover is observed in this paper at intermediate and high temperatures in the backward reaction of the LiNC $\rightleftharpoons$ LiCN isomerization via Langevin dynamics and mean first-passage times (MFPTs). It is in good agreement with the Pollak-Grabert-Hänggi (PGH) reaction rates at lower temperatures. Furthermore, we find a square root behavior of the reaction rate at high temperatures and have made direct comparisons of the methods in the intermediate- and high- temperature regimes; all suggesting increased ranges in accuracy of both the PGH and MFPT approaches.

preprint2022arXiv

Transition state theory characterizes thin film macrospin dynamics driven by an oscillatory magnetic field: Inertial effects

Understanding the magnetization switching process in ferromagnetic thin films is essential for many technological applications. We investigate the effects of periodic driving via magnetic fields on a macrospin system under explicit consideration of inertial dynamics. This is usually achieved by extending the Landau-Lifshitz-Gilbert equation with a term including the second time derivative of the magnetization. The dynamics of the magnetization switching can then be characterized by its switching rate. We apply methods from transition state theory for driven systems to resolve the rate of magnetization switching in this general case. In doing so, we find that magnetization exhibits resonance-like behavior under certain driving conditions, and it can be affected strongly by the system's relaxation rate.

preprint2021arXiv

Dynamics and decay rates of a time-dependent two-saddle system

The framework of transition state theory (TST) provides a powerful way for analyzing the dynamics of physical and chemical reactions. While TST has already been successfully used to obtain reaction rates for systems with a single time-dependent saddle point, multiple driven saddles have proven challenging because of their fractal-like phase space structure. This paper presents the construction of an approximately recrossing-free dividing surface based on the normally hyperbolic invariant manifold in a time-dependent two-saddle model system. Based on this, multiple methods for obtaining instantaneous (time-resolved) decay rates of the underlying activated complex are presented and their results discussed.

preprint2021arXiv

Identifying Reaction Pathways in Phase Space via Asymptotic Trajectories

In this paper, we revisit the concepts of the reactivity map and the reactivity bands as an alternative to the use of perturbation theory for the determination of the phase space geometry of chemical reactions. We introduce a reformulated metric, called the asymptotic trajectory indicator, and an efficient algorithm to obtain reactivity boundaries. We demonstrate that this method has sufficient accuracy to reproduce phase space structures such as turnstiles for a 1D model of the isomerization of ketene in an external field. The asymptotic trajectory indicator can be applied to higher dimensional systems coupled to Langevin baths as we demonstrate for a 3D model of the isomerization of ketene.

preprint2021arXiv

Thermal decay rates of an activated complex in a driven model chemical reaction

Recent work has shown that in a non-thermal, multidimensional system, the trajectories in the activated complex possess different instantaneous and time-averaged reactant decay rates. Under dissipative dynamics, it is known that these trajectories, which are bound on the normally hyperbolic invariant manifold (NHIM), converge to a single trajectory over time. By subjecting these dissipative systems to thermal noise, we find fluctuations in the saddle-bound trajectories and their instantaneous decay rates. Averaging over these instantaneous rates results in the decay rate of the activated complex in a thermal system. We find, that the temperature dependence of the activated complex decay in a thermal system can be linked to the distribution of the phase space resolved decay rates on the NHIM in the non-dissipative case. By adjusting the external driving of the reaction, we show that it is possible to influence how the decay rate of the activated complex changes with rising temperature.

preprint2020arXiv

Neural network approach for the dynamics on the normally hyperbolic invariant manifold of periodically driven systems

Chemical reactions in multidimensional systems are often described by a rank-1 saddle, whose stable and unstable manifolds intersect in the normally hyperbolic invariant manifold (NHIM). Trajectories started on the NHIM in principle never leave this manifold when propagated forward or backward in time. However, the numerical investigation of the dynamics on the NHIM is difficult because of the instability of the motion. We apply a neural network to describe time-dependent NHIMs and use this network to stabilize the motion on the NHIM for a periodically driven model system with two degrees of freedom. The method allows us to analyze the dynamics on the NHIM via Poincaré surfaces of section (PSOS) and to determine the transition state (TS) trajectory as a periodic orbit with the same periodicity as the driving saddle, viz. a fixed point of the PSOS surrounded by near-integrable tori. Based on Transition State Theory and a Floquet analysis of a periodic TS trajectory we compute the rate constant of the reaction with significantly reduced numerical effort compared to the propagation of a large trajectory ensemble.

preprint2020arXiv

Phase-space resolved rates in driven multidimensional chemical reactions

Chemical reactions in multidimensional driven systems are typically described by a time-dependent rank-1 saddle associated with one reaction and several orthogonal coordinates (including the solvent bath). To investigate reactions in such systems, we develop a fast and robust method ---viz., local manifold analysis (LMA)--- for computing the instantaneous decay rate of reactants. Specifically, it computes the instantaneous decay rates along saddle-bound trajectories near the activated complex by exploiting local properties of the stable and unstable manifold associated with the normally hyperbolic invariant manifold (NHIM). The LMA method offers substantial reduction of numerical effort and increased reliability in comparison to direct ensemble integration. It provides an instantaneous flux that can be assigned to every point on the NHIM and which is associated with a trajectory ---regardless of whether it is periodic, quasi-periodic, or chaotic--- that is bound on the NHIM. The time average of these fluxes in the driven system corresponds to the average rate through a given local section containing the corresponding point on the NHIM. We find good agreement between the results of the LMA and direct ensemble integration obtained using numerically constructed, recrossing-free dividing surfaces.

preprint2016arXiv

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$\rightleftharpoons$LiCN isomerization.

preprint2015arXiv

Chemical reactions induced by oscillating external fields in weak thermal environments

Chemical reaction rates must increasingly be determined in systems that evolve under the control of external stimuli. In these systems, when a reactant population is induced to cross an energy barrier through forcing from a temporally varying external field, the transition state that the reaction must pass through during the transformation from reactant to product is no longer a fixed geometric structure, but is instead time-dependent. For a periodically forced model reaction, we develop a recrossing-free dividing surface that is attached to a transition state trajectory [T. Bartsch, R. Hernandez, and T. Uzer, Phys. Rev. Lett. 95, 058301 (2005)]. We have previously shown that for single-mode sinusoidal driving, the stability of the time-varying transition state directly determines the reaction rate [G. T. Craven, T. Bartsch, and R. Hernandez, J. Chem. Phys. 141, 041106 (2014)]. Here, we extend our previous work to the case of multi-mode driving waveforms. Excellent agreement is observed between the rates predicted by stability analysis and rates obtained through numerical calculation of the reactive flux. We also show that the optimal dividing surface and the resulting reaction rate for a reactive system driven by weak thermal noise can be approximated well using the transition state geometry of the underlying deterministic system. This agreement persists as long as the thermal driving strength is less than the order of that of the periodic driving. The power of this result is its simplicity. The surprising accuracy of the time-dependent noise-free geometry for obtaining transition state theory rates in chemical reactions driven by periodic fields reveals the dynamics without requiring the cost of brute-force calculations.

preprint2015arXiv

Transition state trajectory stability determines barrier crossing rates in chemical reactions induced by time-dependent oscillating fields

When a chemical reaction is driven by an external field, the transition state that the system must pass through as it changes from reactant to product -for example, an energy barrier- becomes time-dependent. We show that for periodic forcing the rate of barrier crossing can be determined through stability analysis of the non-autonomous transition state. Specifically, strong agreement is observed between the difference in the Floquet exponents describing stability of the transition state trajectory, which defines a recrossing-free dividing surface [G. T. Craven, T. Bartsch, and R. Hernandez, Phys. Rev. E 89, 040801(R) (2014)], and the rates calculated by simulation of ensembles of trajectories. This result opens the possibility to extract rates directly from the intrinsic stability of the transition state, even when it is time-dependent, without requiring a numerically-expensive simulation of the long-time dynamics of a large ensemble of trajectories.

preprint2014arXiv

Persistence of transition state structure in chemical reactions driven by fields oscillating in time

Chemical reactions subjected to time-varying external forces cannot generally be described through a fixed bottleneck near the transition state barrier or dividing surface. A naive dividing surface attached to the instantaneous, but moving, barrier top also fails to be recrossing-free. We construct a moving dividing surface in phase space over a transition state trajectory. This surface is recrossing-free for both Hamiltonian and dissipative dynamics. This is confirmed even for strongly anharmonic barriers using simulation. The power of transition state theory is thereby applicable to chemical reactions and other activated processes even when the bottlenecks are time-dependent and move across space.