Ultrafast fragmentation of highly-excited doubly-ionized deoxyribose: role of the liquid water environment.

Ab initio molecular dynamics simulations are used to investigate the fragmentation dynamics following the double ionization of 2-deoxy-D-ribose (DR), a major component in the DNA chain. Different ionization scenarios are considered to provide a complete picture. First focusing on isolated DR2+, fragmentation patterns are determined for the ground electronic state, adding randomly distributed excitation energy to the nuclei. These patterns differ for the two isomers studied. To compare thermal and electronic excitation effects, Ehrenfest dynamics are also performed, allowing to remove the two electrons from selected molecular orbitals. Two intermediate-energy orbitals, localized on the carbon chain, were selected. The dissociation pattern corresponds to the most frequent pattern obtained when adding thermal excitation. On the contrary, targeting the four deepest orbitals, localized on the oxygen atoms, leads to selective ultrafast C-O and/or O-H bond dissociation. To probe the role of environment, a system consisting of a DR molecule embedded in liquid water is then studied. The two electrons are removed from either the DR or the water molecules directly linked to the sugar through hydrogen bonds. Although the dynamics onset is similar to that of isolated DR when removing the same deep orbitals localized on the sugar oxygen atoms, the subsequent fragmentation patterns differ. Sugar damage also occurs following the Coulomb explosion of neighboring H2O2+ molecules due to interaction with the emitted O or H atoms.

Time-Resolved Graphs of Polymorphic Cycles for H-Bonded Network Identification in Flexible Biomolecules.

A novel approach based on a coarse-grained representation of topological graphs is proposed for the automatic analysis of molecular dynamics (MD) trajectories of hydrogen-bonded (H-Bonded) flexible biomolecules. Herein, our approach models an H-Bonded biomolecule by its H-Bonded cycles and its graph of cycles in which the vertices and links represent the intersections between these cycles. We propose a methodology in which each identified conformer/isomer from the MD is represented by a well-chosen set of H-Bonded cycles called a minimum cycle basis. The key component is the "polycycles" that distinguish the cycles that play the same polymorphic role in the molecule from the ones that lead to an actual conformational change of the molecule. The relevance of our proposed method is evaluated on MD trajectories of gas-phase biomolecules, for which the covalent bonds are unchanged over time and only the hydrogen bonds change over time. The polygraphs and their time evolution are shown to reveal the dynamicity of the metastructure(s) of the H-Bonded biomolecules while providing polymorphic information on the cycles. Such information on the dynamics and changes in the H-bond network, as some cycles change identity while retaining the same role in the overall structure, is not easily captured at the atomic level of representation. Such information can instead be captured by polymorphic cycles.

HiREX: High-Throughput Reactivity Exploration for Extended Databases of Transition-Metal Catalysts.

A method is introduced for the automated analysis of reactivity exploration for extended in silico databases of transition-metal catalysts. The proposed workflow is designed to tackle two key challenges for bias-free mechanistic explorations on large databases of catalysts: (1) automated exploration of the chemical space around each catalyst with unique structural and chemical features and (2) automated analysis of the resulting large chemical data sets. To address these challenges, we have extended the application of our previously developed ReNeGate method for bias-free reactivity exploration and implemented an automated analysis procedure to identify the classes of reactivity patterns within specific catalyst groups. Our procedure applied to an extended series of representative Mn(I) pincer complexes revealed correlations between structural and reactive features, pointing to new channels for catalyst transformation under the reaction conditions. Such an automated high-throughput virtual screening of systematically generated hypothetical catalyst data sets opens new opportunities for the design of high-performance catalysts as well as an accelerated method for expert bias-free high-throughput in silico reactivity exploration.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment.

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Salt Bridge Structure of Microhydrated Arginine Kinetically Trapped in the Gas Phase by Evaporative Cooling.

Amino acids and peptides generally exhibit zwitterionic forms with salt bridge (SB) structures in solution but charge-solvated (CS) motifs in the gas phase. Here, we report a study of non-covalent complexes of the protonated amino acid arginine, ArgH+(H2O)n (n = 1-5), produced in the gas phase from an aqueous solution with a controlled number of retained water molecules. These complexes were probed by cold ion spectroscopy and treated by quantum chemistry. The spectroscopic changes induced upon gradual dehydration of arginine were assigned by structural calculations to the transition from SB to CS geometries. SB conformers appear to be present for the complexes with as few as 3 retained water molecules, although energetically CS structures should become prevailing already for ArgH+ with 7-8 water molecules. We attribute the revealed kinetic trapping of arginine in native-like zwitterionic forms to evaporative cooling of the hydrated complexes to as low as below 200 K.

Algorithmic Graph Theory, Reinforcement Learning and Game Theory in MD Simulations: From 3D Structures to Topological 2D-Molecular Graphs (2D-MolGraphs) and Vice Versa.

This paper reviews graph-theory-based methods that were recently developed in our group for post-processing molecular dynamics trajectories. We show that the use of algorithmic graph theory not only provides a direct and fast methodology to identify conformers sampled over time but also allows to follow the interconversions between the conformers through graphs of transitions in time. Examples of gas phase molecules and inhomogeneous aqueous solid interfaces are presented to demonstrate the power of topological 2D graphs and their versatility for post-processing molecular dynamics trajectories. An even more complex challenge is to predict 3D structures from topological 2D graphs. Our first attempts to tackle such a challenge are presented with the development of game theory and reinforcement learning methods for predicting the 3D structure of a gas-phase peptide.

On the Trail of Molecular Hydrophilicity and Hydrophobicity at Aqueous Interfaces.

Uncovering microscopic hydrophilicity and hydrophobicity at heterogeneous aqueous interfaces is essential as it dictates physico/chemical properties such as wetting, the electrical double layer, and reactivity. Several molecular and spectroscopic descriptors were proposed, but a major limitation is the lack of connections between them. Here, we combine density functional theory-based MD simulations (DFT-MD) and SFG spectroscopy to explore how interfacial water responds in contact with self-assembled monolayers (SAM) of tunable hydrophilicity. We introduce a microscopic metric to track the transition from hydrophobic to hydrophilic interfaces. This metric combines the H/V descriptor, a structural descriptor based on the preferential orientation within the water network in the topmost binding interfacial layer (BIL) and spectroscopic fingerprints of H-bonded and dangling OH groups of water carried by BIL-resolved SFG spectra. This metric builds a bridge between molecular descriptors of hydrophilicity/hydrophobicity and spectroscopically measured quantities and provides a recipe to quantitatively or qualitatively interpret experimental SFG signals.

ReNeGate: A Reaction Network Graph-Theoretical Tool for Automated Mechanistic Studies in Computational Homogeneous Catalysis.

Exploration of the chemical reaction space of chemical transformations in multicomponent mixtures is one of the main challenges in contemporary computational chemistry. To remove expert bias from mechanistic studies and to discover new chemistries, an automated graph-theoretical methodology is proposed, which puts forward a network formalism of homogeneous catalysis reactions and utilizes a network analysis tool for mechanistic studies. The method can be used for analyzing trajectories with single and multiple catalytic species and can provide unique conformers of catalysts including multinuclear catalyst clusters along with other catalytic mixture components. The presented three-step approach has the integrated ability to handle multicomponent catalytic systems of arbitrary complexity (mixtures of reactants, catalyst precursors, ligands, additives, and solvents). It is not limited to predefined chemical rules, does not require prealignment of reaction mixture components consistent with a reaction coordinate, and is not agnostic to the chemical nature of transformations. Conformer exploration, reactive event identification, and reaction network analysis are the main steps taken for identifying the pathways in catalytic systems given the starting precatalytic reaction mixture as the input. Such a methodology allows us to efficiently explore catalytic systems in realistic conditions for either previously observed or completely unknown reactive events in the context of a network representing different intermediates. Our workflow for the catalytic reaction space exploration exclusively focuses on the identification of thermodynamically feasible conversion channels, representative of the (secondary) catalyst deactivation or inhibition paths, which are usually most difficult to anticipate based solely on expert chemical knowledge. Thus, the expert bias is sought to be removed at all steps, and the chemical intuition is limited to the choice of the thermodynamic constraint imposed by the applicable experimental conditions in terms of threshold energy values for allowed transformations. The capabilities of the proposed methodology have been tested by exploring the reactivity of Mn complexes relevant for catalytic hydrogenation chemistry to verify previously postulated activation mechanisms and unravel unexpected reaction channels relevant to rare deactivation events.

Unexpectedly High Capacitance of the Metal Nanoparticle/Water Interface: Molecular-Level Insights into the Electrical Double Layer.

The electrical double-layer plays a key role in important interfacial electrochemical processes from catalysis to energy storage and corrosion. Therefore, understanding its structure is crucial for the progress of sustainable technologies. We extract new physico-chemical information on the capacitance and structure of the electrical double-layer of platinum and gold nanoparticles at the molecular level, employing single nanoparticle electrochemistry. The charge storage ability of the solid/liquid interface is larger by one order-of-magnitude than predicted by the traditional mean-field models of the double-layer such as the Gouy-Chapman-Stern model. Performing molecular dynamics simulations, we investigate the possible relationship between the measured high capacitance and adsorption strength of the water adlayer formed at the metal surface. These insights may launch the active tuning of solid-solvent and solvent-solvent interactions as an innovative design strategy to transform energy technologies towards superior performance and sustainability.

Stripping away ion hydration shells in electrical double-layer formation: Water networks matter.

The double layer at the solid/electrolyte interface is a key concept in electrochemistry. Here, we present an experimental study combined with simulations, which provides a molecular picture of the double-layer formation under applied voltage. By THz spectroscopy we are able to follow the stripping away of the cation/anion hydration shells for an NaCl electrolyte at the Au surface when decreasing/increasing the bias potential. While Na+ is attracted toward the electrode at the smallest applied negative potentials, stripping of the Cl- hydration shell is observed only at higher potential values. These phenomena are directly measured by THz spectroscopy with ultrabright synchrotron light as a source and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations.

Some opinions on MD-based vibrational spectroscopy of gas phase molecules and their assembly: An overview of what has been achieved and where to go.

We hereby review molecular dynamics simulations for anharmonic gas phase spectroscopy and provide some of our opinions of where the field is heading. With these new directions, the theoretical IR/Raman spectroscopy of large (bio)-molecular systems will be more easily achievable over longer time-scale MD trajectories for an increase in accuracy of the MD-IR and MD-Raman calculated spectra. With the new directions presented here, the high throughput 'decoding' of experimental IR/Raman spectra into 3D-structures should thus be possible, hence advancing e.g. the field of MS-IR for structural characterization by spectroscopy. We also review the assignment of vibrational spectra in terms of anharmonic molecular modes from the MD trajectories, and especially introduce our recent developments based on Graph Theory algorithms. Graph Theory algorithmic is also introduced in this review for the identification of the molecular 3D-structures sampled over MD trajectories.

Molecular Fingerprints of Hydrophobicity at Aqueous Interfaces from Theory and Vibrational Spectroscopies.

Hydrophobicity/hydrophilicity of aqueous interfaces at the molecular level results from a subtle balance in the water-water and water-surface interactions. This is characterized here via density functional theory-molecular dynamics (DFT-MD) coupled with vibrational sum frequency generation (SFG) and THz-IR absorption spectroscopies. We show that water at the interface with a series of weakly interacting materials is organized into a two-dimensional hydrogen-bonded network (2D-HB-network), which is also found above some macroscopically hydrophilic silica and alumina surfaces. These results are rationalized through a descriptor that measures the number of "vertical" and "horizontal" hydrogen bonds formed by interfacial water, quantifying the competition between water-surface and water-water interactions. The 2D-HB-network is directly revealed by THz-IR absorption spectroscopy, while the competition of water-water and water-surface interactions is quantified from SFG markers. The combination of SFG and THz-IR spectroscopies is thus found to be a compelling tool to characterize the finest details of molecular hydrophobicity at aqueous interfaces.

Ions Tune Interfacial Water Structure and Modulate Hydrophobic Interactions at Silica Surfaces.

The structure and ultrafast dynamics of the electric double layer (EDL) are central to chemical reactivity and physical properties at solid/aqueous interfaces. While the Gouy-Chapman-Stern model is widely used to describe EDLs, it is solely based on the macroscopic electrostatic attraction of electrolytes for the charged surfaces. Structure and dynamics in the Stern layer are, however, more complex because of competing effects due to the localized surface charge distribution, surface-solvent-ion correlations, and the interfacial hydrogen bonding environment. Here, we report combined time-resolved vibrational sum frequency generation (TR-vSFG) spectroscopy with ab initio DFT-based molecular dynamics simulations (AIMD/DFT-MD) to get direct access to the molecular-level understanding of how ions change the structure and dynamics of the EDL. We show that innersphere adsorbed ions tune the hydrophobicity of the silica-aqueous interface by shifting the structural makeup in the Stern layer from dominant water-surface interactions to water-water interactions. This drives an initially inhomogeneous interfacial water coordination landscape observed at the neat interface toward a homogeneous, highly interconnected in-plane 2D hydrogen bonding (2D-HB) network at the ionic interface, reminiscent of the canonical, hydrophobic air-water interface. This ion-induced transformation results in a characteristic decrease of the vibrational lifetime (T1) of excited interfacial O-H stretching modes from T1 ∼ 600 fs to T1 ∼ 250 fs. Hence, we propose that the T1 determined by TR-vSFG in combination with DFT-MD simulations can be widely used for a quantitative spectroscopic probe of the ion kosmotropic/chaotropic effect at aqueous interfaces as well as of the ion-induced surface hydrophobicity.

Enhanced conductivity of water at the electrified air-water interface: a DFT-MD characterization.

DFT-based molecular dynamics simulations of the electrified air-liquid water interface are presented, where a homogeneous field is applied parallel to the surface plane. We unveil the field intensity for the onset of proton transfer and molecular dissociation; the protonic current/proton conductivity is measured as a function of the field intensity/voltage. The air-water interface is shown to exhibit a proton conductivity twice the one in the liquid water for field intensities below 0.40 V Å-1. We show that this difference arises from the very specific organization of water in the binding interfacial layer (BIL, i.e. the air-water interface region) into a 2D-HBond-network that is maintained and enforced at the electrified interface. Beyond fields of 0.40 V Å-1, water in the BIL and in the bulk liquid are aligned in the same way by the rather intense fields, hence leading to the same proton conductivity in both BIL and bulk water.

Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain.

Gas-phase, double resonance IR spectroscopy has proven to be an excellent approach to obtain structural information on peptides ranging from single amino acids to large peptides and peptide clusters. In this review, we discuss the state-of-the-art of infrared action spectroscopy of peptides in the far-IR and THz regime. An introduction to the field of far-IR spectroscopy is given, thereby highlighting the opportunities that are provided for gas-phase research on neutral peptides. Current experimental methods, including spectroscopic schemes, have been reviewed. Structural information from the experimental far-IR spectra can be obtained with the help of suitable theoretical approaches such as dynamical DFT techniques and the recently developed Graph Theory. The aim of this review is to underline how the synergy between far-IR spectroscopy and theory can provide an unprecedented picture of the structure of neutral biomolecules in the gas phase. The far-IR signatures of the discussed studies are summarized in a far-IR map, in order to gain insight into the origin of the far-IR localized and delocalized motions present in peptides and where they can be found in the electromagnetic spectrum.

Deconvolution of BIL-SFG and DL-SFG spectroscopic signals reveals order/disorder of water at the elusive aqueous silica interface.

Through the prism of the rather controversial and elusive silica/water interface, ab initio DFT-based molecular dynamics simulations of the structure and non-linear SFG spectroscopy of the interface are analysed. Following our recent work [Phys. Chem. Chem. Phys., 2018, 20, 5190-5199], we show that once the interfacial water is decomposed into BIL (Binding Interfacial Layer) and DL (Diffuse Layer) interfacial regions, the SFG signals can be deconvolved and unambiguously interpreted, and a global microscopic understanding on silica/water interfaces can be obtained. By comparing crystalline quartz/water and amorphous (fused) silica/water interfaces, the dependence of interfacial structural and spectroscopic properties on the degree of surface crystallinity is established, while by adding KCl electrolytes at the quartz/water interface, the chaotropic effect of ions on the interfacial molecular arrangement is unveiled. The evolution of structure and SFG spectra of silica/water interfaces with respect to increasing surface deprotonation, i.e., with respect to pH conditions, is also evaluated. Spectroscopic BIL-SFG markers that experimentally allow one detect the water order/disorder in the BIL as a function of surface hydroxylation and ion concentration are revealed, while the pH-induced modulations in the experimentally recorded SFG spectra are rationalized in terms of changes in both BIL and DL SFG signatures.