Search

The common geometrical (symplectic) structures of classical mechanics, quantum mechanics, and classical thermodynamics are unveiled with three pictures. These cardinal theories, mainly at the non-relativistic approximation, are the cornerstones for studying chemical dynamics and chemical kinetics. Working in extended phase spaces, we show that the physical states of integrable dynamical systems are depicted by Lagrangian submanifolds embedded in phase space. Observable quantities are calculated by properly transforming the extended phase space onto a reduced space, and trajectories are integrated by solving Hamilton's equations of motion. After defining a Riemannian metric, we can also estimate the length between two states. Local constants of motion are investigated by integrating Jacobi fields and solving the variational linear equations. Diagonalizing the symplectic fundamental matrix, eigenvalues equal to one reveal the number of constants of motion. For conservative systems, geometrical quantum mechanics has proved that solving the Schrödinger equation in extended Hilbert space, which incorporates the quantum phase, is equivalent to solving Hamilton's equations in the projective Hilbert space. In classical thermodynamics, we take entropy and energy as canonical variables to construct the extended phase space and to represent the Lagrangian submanifold. Hamilton's and variational equations are written and solved in the same fashion as in classical mechanics. Solvers based on high-order finite differences for numerically solving Hamilton's, variational, and Schrödinger equations are described. Employing the Hénon-Heiles two-dimensional nonlinear model, representative results for time-dependent, quantum, and dissipative macroscopic systems are shown to illustrate concepts and methods. High-order finite-difference algorithms, despite their accuracy in low-dimensional systems, require substantial computer resources when they are applied to systems with many degrees of freedom, such as polyatomic molecules. We discuss recent research progress in employing Hamiltonian neural networks for solving Hamilton's equations. It turns out that Hamiltonian geometry, shared with all physical theories, yields the necessary and sufficient conditions for the mutual assistance of humans and machines in deep-learning processes.

Classical Hamiltonian trajectories are initiated at random points in phase space on a fixed energy shell of a model two degrees of freedom potential, consisting of two interacting minima in an otherwise flat energy plane of infinite extent. Below the energy of the plane, the dynamics are demonstrably chaotic. However, most of the work in this paper involves trajectories at a fixed energy that is 1% above that of the plane, in which regime the dynamics exhibit behavior characteristic of chaotic scattering. The trajectories are analyzed without reference to the potential, as if they had been generated in a typical direct molecular dynamics simulation. The questions addressed are whether one can recover useful information about the structures controlling the dynamics in phase space from the trajectory data alone, and whether, despite the at least partially chaotic nature of the dynamics, one can make statistically meaningful predictions of trajectory outcomes from initial conditions. It is found that key unstable periodic orbits, which can be identified on the analytical potential, appear by simple classification of the trajectories, and that the specific roles of these periodic orbits in controlling the dynamics are also readily discerned from the trajectory data alone. Two different approaches to predicting trajectory outcomes from initial conditions are evaluated, and it is shown that the more successful of them has ∼90% success. The results are compared with those from a simple neural network, which has higher predictive success (97%) but requires the information obtained from the "by-hand" analysis to achieve that level. Finally, the dynamics, which occur partly on the very flat region of the potential, show characteristics of the much-studied phenomenon called "roaming." On this potential, it is found that roaming trajectories are effectively "failed" periodic orbits and that angular momentum can be identified as a key controlling factor, despite the fact that it is not a strictly conserved quantity. It is also noteworthy that roaming on this potential occurs in the absence of a "roaming saddle," which has previously been hypothesized to be a necessary feature for roaming to occur.

Despite using potential energy surfaces, multivariable functions on molecular configuration space, to comprehend chemical dynamics for decades, the real happenings in molecules occur in phase space, in which the states of a classical dynamical system are completely determined by the coordinates and their conjugate momenta. Theoretical and numerical results are presented, employing alanine dipeptide as a model system, to support the view that geometrical structures in phase space dictate the dynamics of molecules, the fingerprints of which are traced by following the Hamiltonian flow above saddles. By properly selecting initial conditions in alanine dipeptide, we have found internally free rotor trajectories the existence of which can only be justified in a phase space perspective.This article is part of the theme issue 'Modern theoretical chemistry'., (© 2018 The Author(s).)

We report the observation of recoil inversion of the CO (v=0, J_{CO}=66) state in the UV dissociation of lab-frame oriented carbonyl sulfide (OCS). This state is ejected in the opposite direction with respect to all other (>30) states and in absence of any OCS rotation, thus resulting in spatial filtering of this particular high-J rovibrational state. This inversion is caused by resonances occurring in shallow local minima of the molecular potential, which bring the sulfur closer to the oxygen than the carbon atom, and is a striking example where such subtleties severely modify the photofragment trajectories. The resonant behavior is observed only in the photofragment trajectories and not in their population, showing that stereodynamic measurements from oriented molecules offer an indispensable probe for exploring energy landscapes.

In this review we discuss the recently described roaming mechanism for chemical reactions from the point of view of nonlinear dynamical systems in phase space. The recognition of the roaming phenomenon shows the need for further developments in our fundamental understanding of basic reaction dynamics, as is made clear by considering some questions that cut across most studies of roaming: Is the dynamics statistical? Can transition state theory be applied to estimate roaming reaction rates? What role do saddle points on the potential energy surface play in explaining the behavior of roaming trajectories? How do we construct a dividing surface that is appropriate for describing the transformation from reactants to products for roaming trajectories? How should we define the roaming region? We show that the phase space perspective on reaction dynamics provides the setting in which these questions can be properly framed and answered. We illustrate these ideas by considering photodissociation of formaldehyde. The phase-space formulation allows an unambiguous description of all possible reactive events, which also allows us to uncover the phase space mechanism that explains which trajectories roam, as opposed to evolving toward a different reactive event.

The roaming mechanism in the reaction H + MgH →Mg + HH is investigated by classical and quantum dynamics employing an accurate ab initio three-dimensional ground electronic state potential energy surface. The reaction dynamics are explored by running trajectories initialized on a four-dimensional dividing surface anchored on three-dimensional normally hyperbolic invariant manifold associated with a family of unstable orbiting periodic orbits in the entrance channel of the reaction (H + MgH). By locating periodic orbits localized in the HMgH well or involving H orbiting around the MgH diatom, and following their continuation with the total energy, regions in phase space where reactive or nonreactive trajectories may be trapped are found. In this way roaming reaction pathways are deduced in phase space. Patterns similar to periodic orbits projected into configuration space are found for the quantum bound and resonance eigenstates. Roaming is attributed to the capture of the trajectories in the neighborhood of certain periodic orbits. The complex forming trajectories in the HMgH well can either return to the radical channel or "roam" to the MgHH minimum from where the molecule may react.

A trajectory method of calculating tunneling probabilities from phase integrals along straight line tunneling paths, originally suggested by Makri and Miller [J. Chem. Phys. 91, 4026 (1989)] and recently implemented by Truhlar and co-workers [Chem. Sci. 5, 2091 (2014)], is tested for one- and two-dimensional ab initio based potentials describing hydrogen dissociation in the (1)B1 excited electronic state of pyrrole. The primary observables are the tunneling rates in a progression of bending vibrational states lying below the dissociation barrier and their isotope dependences. Several initial ensembles of classical trajectories have been considered, corresponding to the quasiclassical and the quantum mechanical samplings of the initial conditions. It is found that the sampling based on the fixed energy Wigner density gives the best agreement with the quantum mechanical dissociation rates.

We examine the phase space structures that govern reaction dynamics in the absence of critical points on the potential energy surface. We show that in the vicinity of hyperbolic invariant tori, it is possible to define phase space dividing surfaces that are analogous to the dividing surfaces governing transition from reactants to products near a critical point of the potential energy surface. We investigate the problem of capture of an atom by a diatomic molecule and show that a normally hyperbolic invariant manifold exists at large atom-diatom distances, away from any critical points on the potential. This normally hyperbolic invariant manifold is the anchor for the construction of a dividing surface in phase space, which defines the outer or loose transition state governing capture dynamics. We present an algorithm for sampling an approximate capture dividing surface, and apply our methods to the recombination of the ozone molecule. We treat both 2 and 3 degrees of freedom models with zero total angular momentum. We have located the normally hyperbolic invariant manifold from which the orbiting (outer) transition state is constructed. This forms the basis for our analysis of trajectories for ozone in general, but with particular emphasis on the roaming trajectories.

We re-examine the prototypical roaming reaction--hydrogen atom roaming in formaldehyde decomposition--from a phase space perspective. Specifically, we address the question "why do trajectories roam, rather than dissociate through the radical channel?" We describe and compute the phase space structures that define and control all possible reactive events for this reaction, as well as provide a dynamically exact description of the roaming region in phase space. Using these phase space constructs, we show that in the roaming region, there is an unstable periodic orbit whose stable and unstable manifolds define a conduit that both encompasses all roaming trajectories exiting the formaldehyde well and shepherds them toward the H2···CO well.

A model Hamiltonian for the reaction CH4(+) -> CH3(+) + H, parametrized to exhibit either early or late inner transition states, is employed to investigate the dynamical characteristics of the roaming mechanism. Tight/loose transition states and conventional/roaming reaction pathways are identified in terms of time-invariant objects in phase space. These are dividing surfaces associated with normally hyperbolic invariant manifolds (NHIMs). For systems with two degrees of freedom NHIMS are unstable periodic orbits which, in conjunction with their stable and unstable manifolds, unambiguously define the (locally) non-recrossing dividing surfaces assumed in statistical theories of reaction rates. By constructing periodic orbit continuation/bifurcation diagrams for two values of the potential function parameter corresponding to late and early transition states, respectively, and using the total energy as another parameter, we dynamically assign different regions of phase space to reactants and products as well as to conventional and roaming reaction pathways. The classical dynamics of the system are investigated by uniformly sampling trajectory initial conditions on the dividing surfaces. Trajectories are classified into four different categories: direct reactive and non-reactive trajectories, which lead to the formation of molecular and radical products respectively, and roaming reactive and non-reactive orbiting trajectories, which represent alternative pathways to form molecular and radical products. By analysing gap time distributions at several energies, we demonstrate that the phase space structure of the roaming region, which is strongly influenced by nonlinear resonances between the two degrees of freedom, results in nonexponential (nonstatistical) decay.

The diffuse vibrational bands, observed in the ultraviolet photodissociation spectrum of nitrous oxide by exciting the molecule in the first (1)A' state, have recently been attributed to resonances localized mainly in the NN stretch and bend degrees of freedom. To further investigate the origin of this localization, fundamental families of periodic orbits emanating from several stationary points of the (1)A' potential energy surface and bifurcations of them are computed. We demonstrate that center-saddle bifurcations of periodic orbits are the main mechanism for creating stable regions in phase space that can support the partial trapping of the wave packet, and thus they explain the observed spectra. A non-linear mechanical methodology, which involves the calculation of equilibria, periodic orbits, and transition states in normal form coordinates, is applied for an in detail exploration of phase space. The fingerprints of the phase space structures in the quantum world are identified by solving the time dependent Schrödinger equation and calculating autocorrelation functions. This demonstrates that different reaction channels could be controlled if exact knowledge of the phase space structure is available to guide the initial excitation of the molecule.

Cytochrome c oxidase (CcO) catalyzes the four-electron reduction of molecular oxygen to water and couples this reduction to the pumping of four protons through the protein matrix. Water molecules inside the protein are involved in the proton pumping activity as proton carriers. A highly conserved water molecule, among different CcO enzymes, lies between the heme a(3) propionates. Here, we show, by quantum mechanical/molecular mechanical (QM/MM) simulations, that this conserved water molecule can transfer its proton to propionate-A. His403 residue coordinates to the Mg site near the so-called water pool. By both QM/MM and molecular dynamics calculations, we demonstrate that the also conserved His403 residue, adjacent to the heme a(3) propionate-A, plays a role of a valve controlling the protonation state of the propionate-A/Asp399 pair. This, in turn, controls the oxidation state of the heme a(3) iron, linking in this way, the D-proton pathway to the water pool.

Theoretical predictions and assignment of highly excited vibrational states and their organization is one of the most important challenges in molecular spectroscopy. A systematic procedure to investigate such problems is locating the principal families of periodic orbits that emanate from the stationary points of the molecule and then following their evolution with the total energy. This results in constructing continuation/bifurcation diagrams that assist in locating the critical bifurcation energies and to discover new types of vibrational modes. Another parameter that may influence the dynamics of a molecule is isotopic mass substitution. In this article, we investigate the effect of symmetry breaking by isotopic mass substitution of triatomic molecules with C(2v) symmetry in classical and quantum dynamics. Sulfur dioxide and water molecules in their ground electronic state are studied by employing accurate potential energy surfaces. Continuation/bifurcation diagrams of periodic orbits are constructed by varying the energy and the mass of one oxygen atom of sulfur dioxide and one hydrogen atom of a water molecule. The transition from normal-to-local mode vibrations is studied in terms of a pitchfork to a center-saddle elementary bifurcation of periodic orbits. The results presented in this article aim to help the assignment of experimentally obtained spectra.

Cytochrome c oxidase (CcO), found in the inner mitochondrial membranes or in many bacteria, catalyzes the four-electron reduction of molecular oxygen to water. Four protons are pumped across the inner mitochondrial membrane through CcO. In this study, quantum mechanics/molecular mechanics and molecular dynamics calculations are used to probe the spectroscopic characteristics of the ferryl intermediates in the aa(3) CcO/O(2) reaction. These highly elaborate calculations, supported by several calculations on smaller model systems, demonstrate the sensitivity of vibrational frequencies on the Coulombic field of heme a(3) and their dependence on the distance of the adjacent Cu(B) to the heme a(3)-Fe atom. This distance seems to be associated with the protonation state of the heme a(3) propionate A, and we propose that it plays a crucial role on the mechanism of action of CcO. In detail, we link proton pumping activity in CcO enzyme (a) to a multiple (1:1:2) resonance among the frequencies of Fe(IV)=O bond stretching, the breathing mode of Histidine 411, and a bending mode of the His411-Fe(IV)=O species (aa(3) from Paracoccus denitrificans numbering) and (b) to Cu(B) displacement by electrostatic interactions toward the heme a(3) iron. We find that the vibrations of the His411-Fe(IV)=O unit become highly coupled depending on the protonation state of the heme a(3) ring A propionate/Asp399 pair, and we propose a mechanism for the resonance Raman enhancement of the bending mode delta(His411-Fe(IV)=O). Calculations on model systems demonstrate that the position of Cu(B) in relation to heme a(3) iron-oxo plays a crucial role in regulating that resonance. We also discuss the origin of the coupling between bending, delta(His411-Fe(IV)=O) and nu(Fe=O) stretching modes, and the role played by such vibrational coupling interactions or Cu(B) position in controlling functional properties of the enzyme, including electron/proton coupling as well as experimental spectra.

Intracavity molecular dynamics studies of photodissociated carbon monoxide from ba(3)-cytochrome c oxidase have been performed by sampling the phase space with several hundreds of trajectories each integrated up to 100 ps time interval. It is shown that the cis conformation of protonated ring-D propionate of heme a(3) and its trans conformation for the deprotonated species control the CO location by creating two distinct equilibrium states for CO confined in a cavity internal to the distal heme pocket. Thus, these cis (closed gate) and trans (open gate) conformations of heme a(3) propionate D play the role of a switch, opening or closing a gate for confining CO in a cavity internal to the heme pocket or releasing it to a bigger outer cavity. The geometry of the inner cavity and the validity of the potential function employed are further investigated by Density Functional Theory calculations for the active site, potential of mean force curves along the copper-CO bond, as well as with Quantum Mechanics/Molecular Mechanics calculations. In the light of the present study, trajectory scenarios for the dissociation of CO previously suggested from time-resolved infrared spectroscopy are re-examined.

Complexity is inherent in biological molecules not only because of the large number of atoms but also because of their nonlinear interactions responsible for chaotic behaviours, localized motions, and bifurcation phenomena. Thus, versatile spectroscopic techniques have been invented to achieve temporal and spatial resolution to minimize the uncertainties in assigning the spectra of complex molecules. Can we associate spectral lines to specific chemical bonds or species in a large molecule? Can energy stay localized in a bond for a substantial period of time to leave its spectroscopic signature? These longstanding problems are investigated by studying the resonance Raman spectra of ferryl-oxo intermediates of cytochrome c oxidase. The difference spectra of isotopically substituted ferryl oxygen ((16)O minus (18)O) in the cytochrome c oxidase recorded in several laboratories show one or two prominent positive peaks which have not been completely elucidated yet. By applying the hierarchical methods of nonlinear mechanics, and particularly the study of periodic orbits in the active site of the enzyme, in conjunction with molecular dynamics calculations of larger systems which include the embraced active site by the protein and selected protonated/deprotonated conformations of amino acids, we translate the spectral lines to molecular motions. It is demonstrated that for the active site stable periodic orbits exist for a substantial energy range. Families of periodic orbits which are associated with the vibrations of Fe(IV)=O bond mark the regions of phase space where nearby trajectories remain localized, as well as assign the spectral bands of the active site in the protein matrix. We demonstrate that proton movement adjacent to active site, which occurs during the P --> F transition, can lead to significant perturbations of the Fe(IV)=O isotopic difference vibrational spectra in cytochrome c oxidase, without a change in oxidation state of the metal sites. This finding links spectroscopic characteristics to protonation events occurring during enzymatic turnover.

The complete vibrational spectrum of the HO2(X(2)A' ') radical, up to the H + O2 dissociation limit, has been determined quantum mechanically on an accurate potential energy surface (PES), based on approximately 15000 ab initio points at the icMRCI+Q/aug-cc-pVQZ level of theory. The vibrational states are found to be assignable at low energies but become more irregular as the energy approaches the dissociation limit. However, even at very high energies, regularity still exists, in sharp contrast to earlier results based on the double many-body expansion (DMBE) IV potential. Several Fermi resonances have been identified, and the spectrum is fit with a spectroscopic Hamiltonian. In addition, the vibrational dynamics is analyzed using a periodic orbit approach.

Small and large molecules may localize their energy in specific bonds or generally in vibrational modes for extended periods of time, an effect which may have dramatic consequences in reaction dynamics. Periodic orbits offer the means to identify phase space regions with localized motions. The author demonstrate that techniques to locate periodic orbits developed for small molecules can be applied to large molecules such as alanine dipeptide. The widely used empirical force fields are employed and principal families of periodic orbits associated with local-type motions and emanated from the lowest energy minima and saddle points are investigated. Continuation of these families at high energies unravels the stable and unstable regions of phase space as well as elementary bifurcations such as saddle nodes.

The interaction of the glycine radical on the side walls of both armchair and zigzag single walled carbon nanotubes is investigated by density functional theory. It is found that the interaction potential of the N-centered glycine radical with the tubes has a minimum of 16.9 (armchair) and 20.2 (zigzag) kcal/mol with respect to the dissociation products. In contrast, the C-centered radical, which is 22.7 kcal/mol lower in energy than the N-centered radical, does not form stable complexes with both types of carbon nanotubes.

Low-lying resonances of the CH2(a 1A1) system (J=0) in an accurate ab initio potential energy surface are studied using a filter-diagonalization method. The width of these resonances fluctuates by more than two orders of magnitude and on average increases with the energy. Analysis of the resonance states concludes that the unimolecular decay of the excited molecular system near the dissociation threshold is neither mode specific nor statistical state specific. This is apparently due to remnant regularity embedded in the largely chaotic classical phase space, as evidenced by periodic orbit analysis. As a result, the Rice-Ramsperger-Kassel-Marcus and statistical adiabatic channel models overestimate the average unimolecular decay rate. The implications of the resonances for the bimolecular C(1D)+H2 reaction are also discussed.

The Huggins band of ozone is investigated by means of dynamics calculations using a new (diabatic) potential energy surface for the 3 (1)A'(1B2) state. The good overall agreement of the calculated spectrum of vibrational energies and intensities with the experimental spectrum, especially at low to intermediate excitation energies, is considered as evidence that the Huggins band is due to the two C(s) potential wells of the 1B2 state rather than the single C2v well of the 2 (1)A'(1A1) state. The vibrational assignment of the "cold bands," based on the nodal structure of wave functions, on the whole supports the most recent experimental assignment [J. Chem. Phys. 115, 9311 (2001)]. The quantum mechanical spectrum is analyzed in terms of classical periodic orbits and the structure of the classical phase space., ((c) 2004 American Institute of Physics)

The Huggins band of ozone is investigated by means of exact dynamics calculations using a new (diabatic) potential energy surface for the (1)B(2) state. The remarkable agreement with the measured spectrum strongly suggests that the Huggins band is due to the two C(s) potential wells of the (1)B(2) state. The vibrational assignment, based on the nodal structure of wave functions, supports the most recent experimental assignment., ((c) 2004 American Institute of Physics.)