Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory.

Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective we celebrate his work at these frontiers.

Molecular Polaritonics: Chemical Dynamics Under Strong Light-Matter Coupling.

Chemical manifestations of strong light-matter coupling have recently been a subject of intense experimental and theoretical studies. Here we review the present status of this field. Section 1 is an introduction to molecular polaritonics and to collective response aspects of light-matter interactions. Section 2 provides an overview of the key experimental observations of these effects, while Section 3 describes our current theoretical understanding of the effect of strong light-matter coupling on chemical dynamics. A brief outline of applications to energy conversion processes is given in Section 4. Pending technical issues in the construction of theoretical approaches are briefly described in Section 5. Finally, the summary in Section 6 outlines the paths ahead in this exciting endeavor.

Capturing Atom-Specific Electronic Structural Dynamics of Transition-Metal Complexes with Ultrafast Soft X-Ray Spectroscopy.

The atomic specificity of X-ray spectroscopies provides a distinct perspective on molecular electronic structure. For 3d metal coordination and organometallic complexes, the combination of metal- and ligand-specific X-ray spectroscopies directly interrogates metal-ligand covalency-the hybridization of metal and ligand electronic states. Resonant inelastic X-ray scattering (RIXS), the X-ray analog of resonance Raman scattering, provides access to all classes of valence excited states in transition-metal complexes, making it a particularly powerful means of characterizing the valence electronic structure of 3d metal complexes. Recent advances in X-ray free-electron laser sources have enabled RIXS to be extended to the ultrafast time domain. We review RIXS studies of two archetypical photochemical processes: charge-transfer excitation in ferricyanide and ligand photodissociation in iron pentacarbonyl. These studies demonstratefemtosecond-resolution RIXS can directly characterize the time-evolving electronic structure, including the evolution of the metal-ligand covalency.

Chemical Reaction Networks' Programming for Solving Equations.

The computational ability of the chemical reaction networks (CRNs) using DNA as the substrate has been verified previously. To solve more complex computational problems and perform the computational steps as expected, the practical design of the basic modules of calculation and the steps in the reactions have become the basic requirements for biomolecular computing. This paper presents a method for solving nonlinear equations in the CRNs with DNA as the substrate. We used the basic calculation module of the CRNs with a gateless structure to design discrete and analog algorithms and realized the nonlinear equations that could not be solved in the previous work, such as exponential, logarithmic, and simple triangle equations. The solution of the equation uses the transformation method, Taylor expansion, and Newton iteration method, and the simulation verified this through examples. We used and improved the basic calculation module of the CRN++ programming language, optimized the error in the basic module, and analyzed the error's variation over time.

Quantum Chaos in the Dynamics of Molecules.

Quantum chaos is reviewed from the viewpoint of "what is molecule?", particularly placing emphasis on their dynamics. Molecules are composed of heavy nuclei and light electrons, and thereby the very basic molecular theory due to Born and Oppenheimer gives a view that quantum electronic states provide potential functions working on nuclei, which in turn are often treated classically or semiclassically. Therefore, the classic study of chaos in molecular science began with those nuclear dynamics particularly about the vibrational energy randomization within a molecule. Statistical laws in probabilities and rates of chemical reactions even for small molecules of several atoms are among the chemical phenomena requiring the notion of chaos. Particularly the dynamics behind unimolecular decomposition are referred to as Intra-molecular Vibrational energy Redistribution (IVR). Semiclassical mechanics is also one of the main research fields of quantum chaos. We herein demonstrate chaos that appears only in semiclassical and full quantum dynamics. A fundamental phenomenon possibly giving birth to quantum chaos is "bifurcation and merging" of quantum wavepackets, rather than "stretching and folding" of the baker's transformation and the horseshoe map as a geometrical foundation of classical chaos. Such wavepacket bifurcation and merging are indeed experimentally measurable as we showed before in the series of studies on real-time probing of nonadiabatic chemical reactions. After tracking these aspects of molecular chaos, we will explore quantum chaos found in nonadiabatic electron wavepacket dynamics, which emerges in the realm far beyond the Born-Oppenheimer paradigm. In this class of chaos, we propose a notion of Intra-molecular Nonadiabatic Electronic Energy Redistribution (INEER), which is a consequence of the chaotic fluxes of electrons and energy within a molecule.

Roaming Reactions and Dynamics in the van der Waals Region.

Roaming reactions were first clearly identified in photodissociation of formaldehyde 15 years ago, and roaming dynamics are now recognized as a universal aspect of chemical reactivity. These reactions typically involve frustrated near-dissociation of a quasibound system to radical fragments, followed by reorientation at long range and intramolecular abstraction. The consequences can be unexpected formation of molecular products, depletion of the radical pool in chemical systems, and formation of products with unusual internal state distributions. In this review, I examine some current aspects of roaming reactions with an emphasis on experimental results, focusing on possible quantum effects in roaming and roaming dynamics in bimolecular systems. These considerations lead to a more inclusive definition of roaming reactions as those for which key dynamics take place at long range.

Combined Crossed Molecular Beams and Ab Initio Study of the Bimolecular Reaction of Ground State Atomic Silicon (Si; 3 P) with Germane (GeH4 ; X1 A1 ).

The chemical dynamics of the elementary reaction of ground state atomic silicon (Si; 3 P) with germane (GeH4 ; X1 A1 ) were unraveled in the gas phase under single collision condition at a collision energy of 11.8±0.3 kJ mol-1 exploiting the crossed molecular beams technique contemplated with electronic structure calculations. The reaction follows indirect scattering dynamics and is initiated through an initial barrierless insertion of the silicon atom into one of the four chemically equivalent germanium-hydrogen bonds forming a triplet collision complex (HSiGeH3 ; 3 i1). This intermediate underwent facile intersystem crossing (ISC) to the singlet surface (HSiGeH3 ; 1 i1). The latter isomerized via at least three hydrogen atom migrations involving exotic, hydrogen bridged reaction intermediates eventually leading to the H3 SiGeH isomer i5. This intermediate could undergo unimolecular decomposition yielding the dibridged butterfly-structured isomer 1 p1 (Si(µ-H2 )Ge) plus molecular hydrogen through a tight exit transition state. Alternatively, up to two subsequent hydrogen shifts to i6 and i7, followed by fragmentation of each of these intermediates, could also form 1 p1 (Si(µ-H2 )Ge) along with molecular hydrogen. The overall non-adiabatic reaction dynamics provide evidence on the existence of exotic dinuclear hydrides of main group XIV elements, whose carbon analog structures do not exist.

Insights from chemical systems into Turing-type morphogenesis.

In 1952, Alan Turing proposed a theory showing how morphogenesis could occur from a simple two morphogen reaction-diffusion system [Turing, A. M. (1952) Phil. Trans. R. Soc. Lond. A 237, 37-72. (doi:10.1098/rstb.1952.0012)]. While the model is simple, it has found diverse applications in fields such as biology, ecology, behavioural science, mathematics and chemistry. Chemistry in particular has made significant contributions to the study of Turing-type morphogenesis, providing multiple reproducible experimental methods to both predict and study new behaviours and dynamics generated in reaction-diffusion systems. In this review, we highlight the historical role chemistry has played in the study of the Turing mechanism, summarize the numerous insights chemical systems have yielded into both the dynamics and the morphological behaviour of Turing patterns, and suggest future directions for chemical studies into Turing-type morphogenesis. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

The Barrier to Proton Transfer in the Dimer of Formic Acid: A Pure Rotational Study.

The rotational spectra of three C-deuterated isotopologues of the dimer of formic acid have been measured, thanks to the small dipole moment induced by asymmetric H→D substitution(s). For the DCOOH-HCOOH species, the concerted double proton transfer of the two hydroxy hydrogen atoms takes place between two equivalent minima and generates a tunneling splitting of 331.2(6) MHz. This splitting can be reproduced by a 3D model with a barrier of 2559 cm-1 (30.6 kJ mol-1 ) as obtained from theoretical calculations.

The Second Law of Thermodynamics as a Force Law.

The second law of thermodynamics states the increase of entropy, Δ S > 0 , for real processes from state A to state B at constant energy from chemistry over biological life and engines to cosmic events. The connection of entropy to information, phase-space, and heat is helpful but does not immediately convince observers of the validity and basis of the second law. This gave grounds for finding a rigorous, but more easily acceptable reformulation. Here, we show using statistical mechanics that this principle is equivalent to a force law ⟨ ⟨ f ⟩ ⟩ > 0 in systems where mass centers and forces can be identified. The sign of this net force--the average mean force along a path from A to B--determines the direction of the process. The force law applies to a wide range of processes from machines to chemical reactions. The explanation of irreversibility by a driving force appears more plausible than the traditional formulation as it emphasizes the cause instead of the effect of motions.

Characterization of Protonated Model Disaccharides from Tandem Mass Spectrometry and Chemical Dynamics Simulations.

The fragmentation mechanisms of prototypical disaccharides have been studied herein by coupling tandem mass spectrometry (MS) with collisional chemical dynamics simulations. These calculations were performed by explicitly considering the collisions between the protonated sugar and the neutral target gas, which led to an ensemble of trajectories for each system, from which it was possible to obtain reaction products and mechanisms without pre-imposing them. The ß-aminoethyl and aminopropyl derivatives of cellobiose, maltose, and gentiobiose were studied to observe differences in both the stereochemistry and the location of the glycosidic linkage. Chemical dynamics simulations of MS/MS and MS/MS/MS were used to suggest some primary and secondary fragmentation mechanisms for some experimentally observed product ions. These simulations provided some new insights into the fundamentals of the unimolecular dissociation of protonated sugars under collisional induced dissociation conditions.

Attosecond vacuum UV coherent control of molecular dynamics.

High harmonic light sources make it possible to access attosecond timescales, thus opening up the prospect of manipulating electronic wave packets for steering molecular dynamics. However, two decades after the birth of attosecond physics, the concept of attosecond chemistry has not yet been realized; this is because excitation and manipulation of molecular orbitals requires precisely controlled attosecond waveforms in the deep UV, which have not yet been synthesized. Here, we present a unique approach using attosecond vacuum UV pulse-trains to coherently excite and control the outcome of a simple chemical reaction in a deuterium molecule in a non-Born-Oppenheimer regime. By controlling the interfering pathways of electron wave packets in the excited neutral and singly ionized molecule, we unambiguously show that we can switch the excited electronic state on attosecond timescales, coherently guide the nuclear wave packets to dictate the way a neutral molecule vibrates, and steer and manipulate the ionization and dissociation channels. Furthermore, through advanced theory, we succeed in rigorously modeling multiscale electron and nuclear quantum control in a molecule. The observed richness and complexity of the dynamics, even in this very simplest of molecules, is both remarkable and daunting, and presents intriguing new possibilities for bridging the gap between attosecond physics and attochemistry.

Testing Turing's theory of morphogenesis in chemical cells.

Alan Turing, in "The Chemical Basis of Morphogenesis" [Turing AM (1952) Philos Trans R Soc Lond 237(641):37-72], described how, in circular arrays of identical biological cells, diffusion can interact with chemical reactions to generate up to six periodic spatiotemporal chemical structures. Turing proposed that one of these structures, a stationary pattern with a chemically determined wavelength, is responsible for differentiation. We quantitatively test Turing's ideas in a cellular chemical system consisting of an emulsion of aqueous droplets containing the Belousov-Zhabotinsky oscillatory chemical reactants, dispersed in oil, and demonstrate that reaction-diffusion processes lead to chemical differentiation, which drives physical morphogenesis in chemical cells. We observe five of the six structures predicted by Turing. In 2D hexagonal arrays, a seventh structure emerges, incompatible with Turing's original model, which we explain by modifying the theory to include heterogeneity.

The Borderline between Reactivity and Pre-reactivity of Binary Mixtures of Gaseous Carboxylic Acids and Alcohols.

By mixing primary and secondary alcohols with carboxylic acids just before the supersonic expansion within pulsed Fourier transform microwave experiments, only the rotational spectrum of the ester was observed. However, when formic acid was mixed with tertiary alcohols, adducts were formed and their rotational spectra could be easily measured. Quantum mechanical calculations were performed to interpret the experimental evidence.

Structural and Chemical Dynamics of Pyridinic-Nitrogen Defects in Graphene.

High density and controllable nitrogen doping in graphene is a critical issue to realize high performance graphene-based devices. In this paper, we demonstrate an efficient method to selectively produce graphitic-N and pyridinic-N defects in graphene by using the mixture plasma of ozone and nitrogen. The atomic structure, electronic structure, and dynamic behavior of these nitrogen defects are systematically studied at the atomic level by using a scanning transmission electron microscopy. The pyridinic-N exhibits higher chemical activity and tends to trap a series of transition metal atoms (Mg, Al, Ca, Ti, Cr, Mn, and Fe) as individual atoms.

An automated method to find transition states using chemical dynamics simulations.

A procedure to automatically find the transition states (TSs) of a molecular system (MS) is proposed. It has two components: high-energy chemical dynamics simulations (CDS), and an algorithm that analyzes the geometries along the trajectories to find reactive pathways. Two levels of electronic structure calculations are involved: a low level (LL) is used to integrate the trajectories and also to optimize the TSs, and a higher level (HL) is used to reoptimize the structures. The method has been tested in three MSs: formaldehyde, formic acid (FA), and vinyl cyanide (VC), using MOPAC2012 and Gaussian09 to run the LL and HL calculations, respectively. Both the efficacy and efficiency of the method are very good, with around 15 TS structures optimized every 10 trajectories, which gives a total of 7, 12, and 83 TSs for formaldehyde, FA, and VC, respectively. The use of CDS makes it a powerful tool to unveil possible nonstatistical behavior of the system under study.

Relative sensitivity of hydrodynamic, thermodynamic, and chemical processes for simulating the buoyant multiphase plume and surfacing flows of an oil and gas blowout.

Deepwater hydrocarbon releases experience complex chemical and physical processes. To assess simplifications of these processes on model predictions, we present a sensitivity analysis using simulations for the Deepwater Horizon oil spill. We compare the buoyant multiphase plume metrics (trap height, rise time etc), the hydrocarbon mass flowrates at the near-field plume termination and their mass fractions dissolved in the water column and reaching the water surface. The baseline simulation utilizes a 19-component hydrocarbon model, live-fluid state equations, hydrate dynamics, and heat and mass transfer. Other simulations turn-off each of these processes, with the simplest one using inert oil and methane gas. Plume metrics are the least sensitive to the modeled processes and can be matched by adjusting the release buoyancy flux. The mass flowrate metrics are more sensitive. Both liquid- and gas-phase mass transfer should be modeled for accurate tracking of soluble components (e.g. C1 - C7 hydrocarbons) in the environment.

Phase Coexistence and Structural Dynamics of Redox Metal Catalysts Revealed by Operando TEM.

Metal catalysts play an important role in industrial redox reactions. Although extensively studied, the state of these catalysts under operating conditions is largely unknown, and assignments of active sites remain speculative. Herein, an operando transmission electron microscopy study is presented, which interrelates the structural dynamics of redox metal catalysts to their activity. Using hydrogen oxidation on copper as an elementary redox reaction, it is revealed how the interaction between metal and the surrounding gas phase induces complex structural transformations and drives the system from a thermodynamic equilibrium toward a state controlled by the chemical dynamics. Direct imaging combined with the simultaneous detection of catalytic activity provides unparalleled structure-activity insights that identify distinct mechanisms for water formation and reveal the means by which the system self-adjusts to changes of the gas-phase chemical potential. Density functional theory calculations show that surface phase transitions are driven by chemical dynamics even when the system is far from a thermodynamic phase boundary. In a bottom-up approach, the dynamic behavior observed here for an elementary reaction is finally extended to more relevant redox reactions and other metal catalysts, which underlines the importance of chemical dynamics for the formation and constant re-generation of transient active sites during catalysis.

Unimolecular Fragmentation Properties of Thermometer Ions from Chemical Dynamics Simulations.

Thermometer ions are widely used to calibrate the internal energy of the ions produced by electrospray ionization in mass spectrometry. Typically, benzylpyridinium ions with different substituents are used. More recently, benzhydrylpyridinium ions were proposed for their lower bond dissociation energies. Direct dynamics simulations using M06-2X/6-31G(d), DFTB, and PM6-D3 are performed to characterize the activation energies of two representative systems: para-methylbenzylpyridinium ion (p-Me-BnPy+) and methyl,methylbenzhydrylpyridinium ion (Me,Me-BhPy+). Simulation results are used to calculate rate constants for the two systems. These rate constants and their uncertainties are used to find the Arrhenius activation energies and RRK fitted threshold energies which give reasonable agreement with calculated bond dissociation energies at the same level of theory. There is only one fragmentation mechanism observed for both systems, which involves C-N bond dissociation via a loose transition state, to generate either benzylium or benzhydrylium ion and a neutral pyridine molecule. For p-Me-BnPy+ using DFTB and PM6-D3 the formation of tropylium ion, from rearrangement of benzylium ion, was observed but only at higher excitation energies and for longer simulation times. These observations suggest that there is no competition between reaction pathways that could affect the reliability of internal energy calibrations. Finally, we suggest using DFTB with a modified-Arrhenius model in future studies.

Role of Chemical Dynamics Simulations in Mass Spectrometry Studies of Collision-Induced Dissociation and Collisions of Biological Ions with Organic Surfaces.

In this article, a perspective is given of chemical dynamics simulations of collisions of biological ions with surfaces and of collision-induced dissociation (CID) of ions. The simulations provide an atomic-level understanding of the collisions and, overall, are in quite good agreement with experiment. An integral component of ion/surface collisions is energy transfer to the internal degrees of freedom of both the ion and the surface. The simulations reveal how this energy transfer depends on the collision energy, incident angle, biological ion, and surface. With energy transfer to the ion's vibration fragmentation may occur, i.e. surface-induced dissociation (SID), and the simulations discovered a new fragmentation mechanism, called shattering, for which the ion fragments as it collides with the surface. The simulations also provide insight into the atomistic dynamics of soft-landing and reactive-landing of ions on surfaces. The CID simulations compared activation by multiple "soft" collisions, resulting in random excitation, versus high energy single collisions and nonrandom excitation. These two activation methods may result in different fragment ions. Simulations provide fragmentation products in agreement with experiments and, hence, can provide additional information regarding the reaction mechanisms taking place in experiment. Such studies paved the way on using simulations as an independent and predictive tool in increasing fundamental understanding of CID and related processes.