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.

Correspondence between classical and quantum resonances

Bifurcations take place in molecular Hamiltonian nonlinear systems as the excitation energy increases, leading to the appearance of different classical resonances. In this paper, we study the quantum manifestations of these classical resonances in the isomerizing system CN-Li⇆Li-CN. By using a correlation diagram of eigenenergies versus Planck constant, we show the existence of different series of avoided crossings, leading to the corresponding series of quantum resonances, which represent the quantum manifestations of the classical resonances. Moreover, the extrapolation of these series to ℏ=0 unveils the correspondence between the bifurcation energy of classical resonances and the energy of the series of quantum resonances in the semiclassical limit ℏ→0. Additionally, in order to obtain analytical expressions for our results, a semiclassical theory is developed.

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.

Using correlation diagrams to study the vibrational spectrum of highly nonlinear floppy molecules: The K-CN case

The correlation diagrams of vibrational energy levels considering the Planck constant as a variable parameter have proven as a very useful tool to study vibrational molecular states, and more specifically in relation to the quantum manifestations of chaos in such dynamical systems. In this paper, we consider the highly nonlinear K-CN molecule, showing how the regular classical structures, i.e., Kolmogorov-Arnold-Moser tori, existing in the mixed classical phase space appear in the quantum levels correlation diagram as emerging diabatic states, something that remains hidden when only the actual value of the Planck constant is considered. Additionally, a quantum transition from order to chaos is unveiled with the aid of these correlation diagrams, where it appears as a frontier of scarred functions.

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.

Shannon entropy at avoided crossings in the quantum transition from order to chaos

Shannon entropy is studied for the series of avoided crossings that characterize the transition from order to chaos in quantum mechanics. In order to be able to study jointly this entropy for discrete and continuous probability, calculations have been performed on a quantized map, the kicked Harper map, resulting in a different behavior, as order-chaos transition takes place, for the discrete (position representation) and continuous (coherent state representation) cases. This different behavior is analyzed in terms of the distribution of zeros of the Husimi function.

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.

Role of short periodic orbits in quantum maps with continuous openings

We apply a recently developed semiclassical theory of short periodic orbits to the continuously open quantum tribaker map. In this paradigmatic system the trajectories are partially bounced back according to continuous reflectivity functions. This is relevant in many situations that include optical microresonators and more complicated boundary conditions. In a perturbative regime, the shortest periodic orbits belonging to the classical repeller of the open map-a cantor set given by a region of exactly zero reflectivity-prove to be extremely robust in supporting a set of long-lived resonances of the continuously open quantum maps. Moreover, for steplike functions a significant reduction in the number needed is obtained, similarly to the completely open situation. This happens despite a strong change in the spectral properties when compared to the discontinuous reflectivity case. In order to give a more realistic interpretation of these results we compare with a Fresnel-type reflectivity function.

Water phase transitions from the perspective of hydrogen-bond network analysis

Analysis of the water phase transitions from the perspective of hydrogen bond networks.

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.

Frequency analysis of the laser driven nonlinear dynamics of HCN

We study the vibrational dynamics of a model for the HCN molecule in the presence of a monochromatic laser field. The variation of the structural behavior of the system as a function of the laser frequency is analyzed in detail using the smaller alignment index, frequency maps, and diffusion coefficients. It is observed that the ergodicity of the system depends on the frequency of the excitation field, especially in its transitions from and into chaos. This provides a roadmap for the possibility of bond excitation and dissociation in this molecule.

Theory of short periodic orbits for partially open quantum maps

We extend the semiclassical theory of short periodic orbits [M. Novaes et al., Phys. Rev. E 80, 035202(R) (2009)PLEEE81539-375510.1103/PhysRevE.80.035202] to partially open quantum maps, which correspond to classical maps where the trajectories are partially bounced back due to a finite reflectivity R. These maps are representative of a class that has many experimental applications. The open scar functions are conveniently redefined, providing a suitable tool for the investigation of this kind of system. Our theory is applied to the paradigmatic partially open tribaker map. We find that the set of periodic orbits that belongs to the classical repeller of the open map (R=0) is able to support the set of long-lived resonances of the partially open quantum map in a perturbative regime. By including the most relevant trajectories outside of this set, the validity of the approximation is extended to a broad range of R values. Finally, we identify the details of the transition from qualitatively open to qualitatively closed behavior, providing an explanation in terms of short periodic orbits.

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.

Agricultural activity shapes the communication and migration patterns in Senegal

The communication and migration patterns of a country are shaped by its socioeconomic processes. The economy of Senegal is predominantly rural, as agriculture employs over 70% of the labor force. In this paper, we use mobile phone records to explore the impact of agricultural activity on the communication and mobility patterns of the inhabitants of Senegal. We find two peaks of phone calls activity emerging during the growing season. Moreover, during the harvest period, we detect an increase in the migration flows throughout the country. However, religious holidays also shape the mobility patterns of the Senegalese people. Hence, in the light of our results, agricultural activity and religious holidays are the primary drivers of mobility inside the country.

Solvated molecular dynamics of LiCN isomerization: All-atom argon solvent versus a generalized Langevin bath

The reaction rate rises and falls with increasing density or friction when a molecule is activated by collisions with the solvent particles. This so-called Kramers turnover has recently been observed in the isomerization reaction of LiCN in an argon bath. In this paper, we demonstrate by direct comparison with those results that a reduced-dimensional (generalized) Langevin description gives rise to similar reaction dynamics as the corresponding (computationally expensive) full molecular dynamics calculations. We show that the density distributions within the Langevin description are in direct agreement with the full molecular dynamics results and that the turnover in the reaction rates is reproduced qualitatively and quantitatively at different temperatures.

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

The minimum contour in the forward Lagrangian descriptor overlaps the invariant manifold (in green) dividing reactant and product regions.

Using the small alignment index chaos indicator to characterize the vibrational dynamics of a molecular system: LiNC-LiCN

A study of the dynamical characteristics of the phase space corresponding to the vibrations of the LiNC-LiCN molecule using an analysis based on the small alignment index (SALI) is presented. SALI is a good indicator of chaos that can easily determine whether a given trajectory is regular or chaotic regardless of the dimensionality of the system, and can also provide a wealth of dynamical information when conveniently implemented. In two-dimensional (2D) systems SALI maps are computed as 2D phase space representations, where the SALI asymptotic values are represented in color scale. We show here how these maps provide full information on the dynamical phase space structure of the LiNC-LiCN system, even quantifying numerically the volume of the different zones of chaos and regularity as a function of the molecule excitation energy.

Measuring political polarization: Twitter shows the two sides of Venezuela

We say that a population is perfectly polarized when divided in two groups of the same size and opposite opinions. In this paper, we propose a methodology to study and measure the emergence of polarization from social interactions. We begin by proposing a model to estimate opinions in which a minority of influential individuals propagate their opinions through a social network. The result of the model is an opinion probability density function. Next, we propose an index to quantify the extent to which the resulting distribution is polarized. Finally, we apply the proposed methodology to a Twitter conversation about the late Venezuelan president, Hugo Chávez, finding a good agreement between our results and offline data. Hence, we show that our methodology can detect different degrees of polarization, depending on the structure of the network.

Geometrical analysis of the LiCN vibrational dynamics: a stability geometrical indicator

The vibrational dynamics of the LiNC/LiCN molecular system is examined making use of the Riemannian geometry. Stability and chaoticity are analyzed, in this context, by means of the Jacobi-Levi-Civita equations, derived from the Jacobi metric, and its solutions. A dynamical indicator, called stability geometrical indicator, is introduced in order to ascertain the dynamical characteristics of stability and chaos in the molecule under study.

Effect of the local morphology in the field emission properties of conducting polymer surfaces

In this work, we present systematic theoretical evidence of a relationship between the point local roughness exponent (PLRE) (which quantifies the heterogeneity of an irregular surface) and the cold field emission properties (indicated by the local current density and the macroscopic current density) of real polyaniline (PANI) surfaces, considered nowadays as very good candidates in the design of field emission devices. The latter are obtained from atomic force microscopy data. The electric field and potential are calculated in a region bounded by the rough PANI surface and a distant plane, both boundaries held at distinct potential values. We numerically solve Laplace's equation subject to appropriate Dirichlet's condition. Our results show that local roughness reveals the presence of specific sharp emitting spots with a smooth geometry, which are the main ones responsible (but not the only) for the emission efficiency of such surfaces for larger deposition times. Moreover, we have found, with a proper choice of a scale interval encompassing the experimentally measurable average grain length, a highly structured dependence of local current density on PLRE, considering different ticks of PANI surfaces.

Onset of quantum chaos in molecular systems and the zeros of the Husimi function

The distribution of zeros of the Husimi function has proven a useful tool in the characterization of regular and chaotic quantum states of dynamical systems. In a previous paper [Phys. Rev. E 82, 026201 (2010)], we showed how the quantum transition from order to chaos in three different molecular systems (LiCN, HCN, and HO(2)) can be understood by means of the correlation diagram of eigenenergies versus the Planck's constant. An order-chaos frontier of scars (eigenstates localized over unstable and stable complementary, in the sense of the Poincaré-Birkhoff theorem, periodic orbits) was observed. In this paper, we show how the distribution of zeros of the Husimi function can be related to the onset of chaos in these molecular systems, showing a very interesting feature at the frontier of scars; namely, some zeros of the Husimi function localize over the stable and unstable fixed points corresponding to the two complementary periodic orbits, this representing the quantum equivalent to the Poincaré-Birkhoff theorem.

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.

Classical transients and the support of open quantum maps

The basic ingredients in a semiclassical theory are the classical invariant objects serving as a support for quantization. Recent studies, mainly obtained on quantum maps, have led to the commonly accepted belief that the classical repeller-the set of nonescaping orbits in the future and past evolution-is the object that suitably plays this role in open scattering systems. In this paper we present numerical evidence warning that this may not always be the case. For this purpose we study recently introduced families of tribaker maps [L. Ermann, G. G. Carlo, J. M. Pedrosa, and M. Saraceno, Phys. Rev. E 85, 066204 (2012)], which share the same asymptotic properties but differ in their short-time behavior. We have found that although the eigenvalue distribution of the evolution operator of these maps follows the fractal Weyl law prediction, the theory of short periodic orbits for open maps fails to describe the resonance eigenfunctions of some of them. This is a strong indication that new elements must be included in the semiclassical description of open quantum systems. We provide an interpretation of the results in order to have hints about them.

Characterizing and modeling an electoral campaign in the context of Twitter: 2011 Spanish Presidential election as a case study

Transmitting messages in the most efficient way as possible has always been one of politicians' main concerns during electoral processes. Due to the rapidly growing number of users, online social networks have become ideal platforms for politicians to interact with their potential voters. Exploiting the available potential of these tools to maximize their influence over voters is one of politicians' actual challenges. To step in this direction, we have analyzed the user activity in the online social network Twitter, during the 2011 Spanish Presidential electoral process, and found that such activity is correlated with the election results. We introduce a new measure to study political sentiment in Twitter, which we call the relative support. We have also characterized user behavior by analyzing the structural and dynamical patterns of the complex networks emergent from the mention and retweet networks. Our results suggest that the collective attention is driven by a very small fraction of users. Furthermore, we have analyzed the interactions taking place among politicians, observing a lack of debate. Finally, we develop a network growth model to reproduce the interactions taking place among politicians.

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.

Poincaré-Birkhoff theorem in quantum mechanics

Quantum manifestations of the dynamics around resonant tori in perturbed hamiltonian systems, dictated by the Poincaré-Birkhoff theorem, are shown to exist. They are embedded in the interactions involving states which differ in a number of quanta equal to the order of the classical resonance. Moreover, the associated classical phase space structures are mimicked in the quasiprobability density functions and their zeros.

Environmental stability of quantum chaotic ratchets

The transitory and stationary behavior of a quantum chaotic ratchet consisting of a biharmonic potential under the effect of different drivings in contact with a thermal environment is studied. For weak forcing and finite ℏ, we identify a strong dependence of the current on the structure of the chaotic region. Moreover, we have determined the robustness of the current against thermal fluctuations in the very weak coupling regime. In the case of strong forcing, the current is determined by the shape of a chaotic attractor. In both cases the temperature quickly stabilizes the ratchet, but in the latter it also destroys the asymmetry responsible for the current generation. Finally, applications to isomerization reactions are discussed.

Scars at the edge of the transition from order to chaos in the isomerizing molecular systems LiNC-LiCN and HCN-HNC, and HO2

Correlation diagrams of energy levels using ℏ as the varying parameter have proven very useful in the characterization of quantum vibrational dynamics of small polyatomic molecules, specially in relation to the transition from order to chaos. In this paper, we present calculations of such correlation diagrams for three molecular systems with very different characteristics, which can be traced down to the topology of their potential energy surfaces. By studying the broad avoided crossings existing in the diagrams, we show that the transition from regular to irregular states always corresponds to a frontier formed by scarred states.

Chaos in the classical mechanics of bound and quasi-bound HX-4He complexes with X = F, Cl, Br, CN

The classical dynamics of weakly bound floppy van der Waals complexes have been extensively studied in the past except for the weakest of all, i.e., those involving He atoms. These complexes are of considerable current interest in light of recent experimental work focussed on the study of molecules trapped in small droplets of the quantum solvent (4)He. Despite a number of quantum investigations, details on the dynamics of how quantum solvation occurs remain unclear. In this paper, the classical rotational dynamics of a series of van der Waals complexes, HX-(4)He with X = F, Cl, Br, CN, are studied. In all cases, the ground state dynamics are found to be almost entirely chaotic, in sharp contrast to other floppy complexes, such as HCl-Ar, for which chaos sets in only at relatively high energies. The consequences of this result for quantum solvation are discussed. We also investigate rotationally excited states with J = 1 which, except for HCN-(4)He, are actually resonances that decay by rotational pre-dissociation.

Solvent-induced acceleration of the rate of activation of a molecular reaction

An increase in the rates of activated processes with the coupling to the solvent has long been predicted through the phenomenological Langevin equation in the weak coupling regime. However, its direct observation in particle-based models has been elusive because the coupling typically places the processes in the spacial-diffusion limited regime wherein rates decrease with increasing friction. In this work, the forward and backward reaction rates of the LiNC<==>LiCN isomerization reaction in a bath of argon atoms at various densities have been calculated directly using molecular dynamics trajectories. The so-called Kramers turnover in the rate with microscopic friction is clearly visible, thus providing direct and unambiguous evidence for the energy-diffusion regime in which rates increase with friction.

Scarring by homoclinic and heteroclinic orbits

In addition to the well-known scarring effect of periodic orbits, we show here that homoclinic and heteroclinic orbits, which are cornerstones in the theory of classical chaos, also scar eigenfunctions of classically chaotic systems when associated closed circuits in phase space are properly quantized, thus introducing strong quantum correlations. The corresponding quantization rules are also established. This opens the door for developing computationally tractable methods to calculate eigenstates of chaotic systems.

Irreversibility with quantum trajectories

Irreversibility is an important issue for many quantum processes. Loschmidt echoes, originally introduced as a way to gauge sensitivity to perturbations in quantum mechanics, have turned out to be a useful tool for its investigation. Following the philosophy supporting this idea, and using quantum trajectories as defined in the causal interpretation of quantum mechanics due to Bohm, we introduce in this paper a more informative alternative measure for irreversibility. The method is applied to the Bunimovich stadium billiard, a paradigmatic example of chaotic system, that constitutes an excellent model for mesoscopic devices.

The onset of chaos in the vibrational dynamics of LiNC∕LiCN

Recent advances in vibrational spectroscopy have greatly enhanced the possibilities of research of highly excited states in molecular systems of moderate size. At sufficiently high level of excitation the correspondence principle holds, and classical mechanical arguments constitute a useful interpretative tool. The corresponding dynamics often become very complex specially in systems with floppy degrees of freedom, and periodic motion plays an important role for its understanding. In this paper, we present a computational procedure to systematically calculate periodic orbits of LiNCLiCN with a given symmetry, that has the additional advantage of providing a useful insight into the onset of chaos in this system.

Vibrational dynamics of the floppy LiNC∕LiCN molecular system

Modern spectroscopical techniques allow the efficient experimental investigation of highly excited vibrational states in molecular systems. On the theoretical side, powerful computational methods have also been developed for the calculation of the corresponding energy levels and wave functions, and their interpretation. In this paper we use a combination of two such methods, namely, the distribution of zeroes in the Husimi function and energy-level correlation diagrams, to discuss a classification scheme, for the lowest hundred vibrational levels of the LiNC/LiCN floppy molecular system, based on their dynamical characteristics.

Homoclinic motions in the vibrational spectra of floppy systems: The LiCN molecule

Recent experimental and theoretical methods allowed the efficient investigation of highly excited rovibrational states of molecular systems. At these levels of excitation the correspondence principle holds, and then classical mechanics can provide intuitive views of the involved processes. In this respect, we have recently shown that for completely hyperbolic systems, homoclinic motions, which are known to organize the classical chaotic region in Hamiltonian systems, imprint a clear signature in the corresponding highly excited quantum spectra. In this Communication we show that this result also holds in mixed systems, by considering an application to the floppy LiNCLiCN molecular system.

Signatures of homoclinic motion in quantum chaos

Homoclinic motion plays a key role in the organization of classical chaos in Hamiltonian systems. In this Letter, we show that it also imprints a clear signature in the corresponding quantum spectra. By numerically studying the fluctuations of the widths of wave functions localized along periodic orbits we reveal the existence of an oscillatory behavior that is explained solely in terms of the primary homoclinic motion. Furthermore, our results indicate that it survives the semiclassical limit.

Classical invariants and the quantization of chaotic systems

Due to their exponential proliferation, long periodic orbits constitute a serious drawback in Gutzwiller's theory of chaotic systems. Therefore, it would be desirable that other classical invariants, not suffering from the same problem, could be used in alternative semiclassical quantization schemes. In this Rapid Communication, we demonstrate how a suitable dynamical analysis of chaotic quantum spectra unveils the role played, in this respect, by classical invariant areas related to the stable and unstable manifolds of short periodic orbits.

Topology of the distribution of zeros of the Husimi function in the LiNC/LiCN molecular system

Phase space representations of quantum mechanics constitute useful tools to study vibrations in molecular systems. Among all possibilities, the Husimi function or coherent state representation is very widely used, its maxima indicating which regions of phase space are relevant in the dynamics of the system. The corresponding zeros are also a good indicator to investigate the characteristics of the eigenstates, and it has been shown how the corresponding distributions can discriminate between regular, irregular, and scarred wave functions. In this paper, we discuss how this result can be understood in terms of the overlap between coherent states and system eigenfunctions.

Different time scales in wave function intensity statistics

Unstable periodic orbits scar wave functions in chaotic systems. The local short term dynamics also influences the associated spectra that follow the otherwise universal Porter-Thomas intensity distribution. We show here how this deviation extends to other longer periodic orbits sharing some common dynamical characteristics. This indicates that the quantum mechanics of the system can be described quite simply with few orbits, up to the resolution associated with the corresponding lengths.