-
Uncovering low-frequency vibrations in surface-enhanced Raman of organic molecules
Authors:
Alexandra Boehmke,
Roberto A Boto,
Eoin Elliot,
Bart de Nijs,
Ruben Esteban,
Tamás Földes,
Fernando Aguilar-Galindo,
Edina Rosta,
Javier Aizpurua,
Jeremy J Baumberg
Abstract:
Accessing the terahertz (THz) spectral domain through surface-enhanced Raman spectroscopy (SERS) is challenging and opens up the study of low-frequency molecular and electronic excitations. Compared to direct THz probing of heterogenous ensembles, the extreme plasmonic confinement of visible light to deep sub-wavelength scales allows the study of hundreds or even single molecules. We show that sel…
▽ More
Accessing the terahertz (THz) spectral domain through surface-enhanced Raman spectroscopy (SERS) is challenging and opens up the study of low-frequency molecular and electronic excitations. Compared to direct THz probing of heterogenous ensembles, the extreme plasmonic confinement of visible light to deep sub-wavelength scales allows the study of hundreds or even single molecules. We show that self-assembled molecular monolayers of a set of simple aromatic thiols confined inside single-particle plasmonic nanocavities can be distinguished by their low-wavenumber spectral peaks below 200 cm-1, after removal of a bosonic inelastic contribution and an exponential background from the spectrum Developing environment-dependent density-functional-theory simulations of the metal-molecule configuration enables the assignment and classification of their THz vibrations as well as the identification of intermolecular coupling effects and of the influence of the gold surface configuration Furthermore, we show dramatically narrower THz SERS spectra from individual molecules at picocavities, which indicates the possibility to study intrinsic vibrational properties beyond inhomogeneous broadening further supporting the key role of local environment.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Kinetics-Optimized Enhanced Sampling Using Mean First Passage Times
Authors:
Tiejun Wei,
Balint Dudas,
Edina Rosta
Abstract:
Molecular dynamics simulations have become essential in many areas of atomistic modelling from drug discovery to materials science. They provide critical atomic-level insights into key dynamical events experiments cannot easily capture. However, their impact often falls short as the timescales of the important processes are inaccessible using standard molecular dynamics. Enhanced sampling methods…
▽ More
Molecular dynamics simulations have become essential in many areas of atomistic modelling from drug discovery to materials science. They provide critical atomic-level insights into key dynamical events experiments cannot easily capture. However, their impact often falls short as the timescales of the important processes are inaccessible using standard molecular dynamics. Enhanced sampling methods provided avenues to access such crucial rare events, for example key slow conformational changes of biomolecules. However, the bias in enhanced sampling simulations is rarely optimized, and even if they are, the optimization criteria is based on the thermodynamics or Hamiltonian of the system, but do not directly consider molecular kinetics. Here, we introduce a novel enhanced sampling algorithm that adaptively optimizes the bias based on the kinetics of the system for the first time. We identify the optimal bias that minimizes a key physical observable, the mean first passage time (MFPT) from a starting state to a target state. Our algorithm makes use of the relation between biased and unbiased kinetics obtained from discretized Markov state models (MSMs), as established in the dynamic histogram analysis method (DHAM). We demonstrate the applicability of the method for different 1D and 2D analytical potential-based model examples, NaCl dissociation in explicit water, and phosphate unbinding in Ras GTPase. Our algorithm has excellent performance compared with state-of-art enhanced sampling methods in terms of the timescales required to reach the final state in the benchmarking systems. Our findings provide a novel, kinetics-driven enhanced sampling strategy, signatured by a targeted approach to facilitate mapping rare events, with the potential for breakthrough applications in drug discovery and materials science.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Large Scale Training of Graph Neural Networks for Optimal Markov-Chain Partitioning Using the Kemeny Constant
Authors:
Sam Alexander Martino,
João Morado,
Chenghao Li,
Zhenghao Lu,
Edina Rosta
Abstract:
Traditional clustering algorithms often struggle to capture the complex relationships within graphs and generalise to arbitrary clustering criteria. The emergence of graph neural networks (GNNs) as a powerful framework for learning representations of graph data provides new approaches to solving the problem. Previous work has shown GNNs to be capable of proposing partitionings using a variety of c…
▽ More
Traditional clustering algorithms often struggle to capture the complex relationships within graphs and generalise to arbitrary clustering criteria. The emergence of graph neural networks (GNNs) as a powerful framework for learning representations of graph data provides new approaches to solving the problem. Previous work has shown GNNs to be capable of proposing partitionings using a variety of criteria, however, these approaches have not yet been extended to work on Markov chains or kinetic networks. These arise frequently in the study of molecular systems and are of particular interest to the biochemical modelling community. In this work, we propose several GNN-based architectures to tackle the graph partitioning problem for Markov Chains described as kinetic networks. This approach aims to minimize how much a proposed partitioning changes the Kemeny constant. We propose using an encoder-decoder architecture and show how simple GraphSAGE-based GNNs with linear layers can outperform much larger and more expressive attention-based models in this context. As a proof of concept, we first demonstrate the method's ability to cluster randomly connected graphs. We also use a linear chain architecture corresponding to a 1D free energy profile as our kinetic network. Subsequently, we demonstrate the effectiveness of our method through experiments on a data set derived from molecular dynamics. We compare the performance of our method to other partitioning techniques such as PCCA+. We explore the importance of feature and hyperparameter selection and propose a general strategy for large-scale parallel training of GNNs for discovering optimal graph partitionings.
△ Less
Submitted 25 December, 2023; v1 submitted 22 December, 2023;
originally announced December 2023.
-
Cooperative O-H$\cdotsπ$ and C-H$\cdots$O Hydrogen Bonding in Benzene-Methanol Solution: Strong Structures from Weak Interactions
Authors:
Camilla Di Mino,
Andrew G. Seel,
Adam J. Clancy,
Andrea Sella,
Thomas F. Headen,
Támas Földes,
Edina Rosta,
Neal T. Skipper
Abstract:
Weak hydrogen bonds, such as O-H$\cdotsπ$ and C-H$\cdots$O, are pivotal in a wide range of important natural and industrial processes including biochemical assembly, molecular recognition, and chemical selectivity. In this study we use neutron diffraction in conjunction with comprehensive H/D isotopic substitution to obtain a detailed spatial and orientational picture of the structure in benzene-m…
▽ More
Weak hydrogen bonds, such as O-H$\cdotsπ$ and C-H$\cdots$O, are pivotal in a wide range of important natural and industrial processes including biochemical assembly, molecular recognition, and chemical selectivity. In this study we use neutron diffraction in conjunction with comprehensive H/D isotopic substitution to obtain a detailed spatial and orientational picture of the structure in benzene-methanol solution. This system provides us with a prototypical situation where the aromatic ring can act as an hydrogen bond acceptor (via the $π$ electron density) and/or a hydrogen bond donor (via the CH groups), with the potential for cooperative effects. Our analysis places benzene at the centre of our frame-of-reference, and reveals for the first time that in solution the O-H$\cdotsπ$ interaction is highly localised and directional, the hydrogen atom being located directly above/below the ring centroid at a distance of 2.30 Å and with the hydroxyl bond axis normal to the aromatic plane. The tendency of methanol to form chain and cyclic motifs in the bulk liquid is manifest in a highly templated, symmetrical equatorial solvation structure; the methanol molecules surround the benzene so that the O-H bonds are coplanar with the aromatic ring while the oxygens interact with C-H groups through simultaneous bifurcated hydrogen bonds. By contrast, C-H$\cdotsπ$ interactions are relegated to the role of more distant spectators. The experimentally observed solvation therefore demonstrates that weak hydrogen bonding can give rise to strongly-ordered cooperative structural motifs also in the liquid phase.
△ Less
Submitted 20 December, 2022;
originally announced December 2022.
-
Molecular Vibration Explorer: an online database and toolbox for surface-enhanced frequency conversion, infrared and Raman spectroscopy
Authors:
Zsuzsanna Koczor-Benda,
Philippe Roelli,
Christophe Galland,
Edina Rosta
Abstract:
We present Molecular Vibration Explorer, a freely accessible online database and interactive tool for exploring vibrational spectra and tensorial light-vibration coupling strengths of a large collection of thiolated molecules. The `Gold' version of the database gathers the results from density functional theory calculations on 2'800 commercially available thiol compounds linked to a gold atom, wit…
▽ More
We present Molecular Vibration Explorer, a freely accessible online database and interactive tool for exploring vibrational spectra and tensorial light-vibration coupling strengths of a large collection of thiolated molecules. The `Gold' version of the database gathers the results from density functional theory calculations on 2'800 commercially available thiol compounds linked to a gold atom, with the main motivation to screen the best molecules for THz and mid-infrared to visible upconversion. Additionally, the `Thiol' version of the database contains results for 1'900 unbound thiolated compounds. They both provide access to a comprehensive set of computed spectroscopic parameters for all vibrational modes of all molecules in the database. The user can simultaneously investigate infrared absorption, Raman scattering and vibrational sum- and difference frequency generation cross sections. Molecules can be screened for various parameters in custom frequency ranges, such as large Raman cross-section under specific molecular orientation, or large orientation-averaged sum-frequency generation (SFG) efficiency. The user can select polarization vectors for the electromagnetic fields, set the orientation of the molecule and customize parameters for plotting the corresponding IR, Raman and sum-frequency spectra. We illustrate the capabilities of this tool with selected applications in the field of surface-enhanced spectroscopy.
△ Less
Submitted 21 December, 2021;
originally announced December 2021.
-
Detecting mid-infrared light by molecular frequency upconversion with dual-wavelength hybrid nanoantennas
Authors:
Angelos Xomalis,
Xuezhi Zheng,
Rohit Chikkaraddy,
Zsuzsanna Koczor-Benda,
Ermanno Miele,
Edina Rosta,
Guy A E Vandenbosch,
Alejandro Martínez,
Jeremy J Baumberg
Abstract:
Coherent interconversion of signals between optical and mechanical domains is enabled by optomechanical interactions. Extreme light-matter coupling produced by confining light to nanoscale mode volumes can then access single mid-infrared (MIR) photon sensitivity. Here we utilise the infrared absorption and Raman activity of molecular vibrations in plasmonic nanocavities to demonstrate frequency up…
▽ More
Coherent interconversion of signals between optical and mechanical domains is enabled by optomechanical interactions. Extreme light-matter coupling produced by confining light to nanoscale mode volumes can then access single mid-infrared (MIR) photon sensitivity. Here we utilise the infrared absorption and Raman activity of molecular vibrations in plasmonic nanocavities to demonstrate frequency upconversion. We convert λ~10 μm incoming light to visible via surface-enhanced Raman scattering (SERS) in doubly-resonant antennas that enhance upconversion by >10^10. We show >200% amplification of the SERS antiStokes emission when a MIR pump is tuned to a molecular vibrational frequency, obtaining lowest detectable powers ~1 μW/μm^2 at room temperature. These results have potential for low-cost and large-scale infrared detectors and spectroscopic techniques, and bring single-molecule sensing into the infrared
△ Less
Submitted 6 July, 2021;
originally announced July 2021.
-
Assessing Position-Dependent Diffusion from Biased Simulations and Markov State Model Analysis
Authors:
Francois Sicard,
Vladimir Koskin,
Alessia Annibale,
Edina Rosta
Abstract:
A variety of enhanced statistical and numerical methods are now routinely used to extract comprehensible and relevant thermodynamic information from the vast amount of complex, high-dimensional data obtained from intensive molecular simulations. The characterization of kinetic properties, such as diffusion coefficients, of molecular systems with significantly high energy barriers, on the other han…
▽ More
A variety of enhanced statistical and numerical methods are now routinely used to extract comprehensible and relevant thermodynamic information from the vast amount of complex, high-dimensional data obtained from intensive molecular simulations. The characterization of kinetic properties, such as diffusion coefficients, of molecular systems with significantly high energy barriers, on the other hand, has received less attention. Among others, Markov state models, in which the long-time statistical dynamics of a system is approximated by a Markov chain on a discrete partition of configuration space, have seen widespread use in recent years, with the aim of tackling these fundamental issues. Here we propose a general, automatic method to assess multidimensional position-dependent diffusion coefficients within the framework of Markovian stochastic processes and Kramers-Moyal expansion. We apply the formalism to one- and two-dimensional analytic potentials and data from explicit solvent molecular dynamics simulations, including the water-mediated conformations of alanine dipeptide. Importantly, the developed algortihm presents significant improvement compared to standard methods when the transport of solute across three-dimensional heterogeneous porous media is studied, for example, the prediction of membrane permeation of drug molecules.
△ Less
Submitted 13 October, 2020;
originally announced October 2020.
-
Breaking the Selection Rules of Spin-Forbidden Molecular Absorption in Plasmonic Nanocavities
Authors:
Oluwafemi S. Ojambati,
William M. Deacon,
Rohit Chikkaraddy,
Charlie Readman,
Qianqi Lin,
Zsuzsanna Koczor-Benda,
Edina Rosta,
Oren A. Scherman,
Jeremy J. Baumberg
Abstract:
Controlling absorption and emission of organic molecules is crucial for efficient light-emitting diodes, organic solar cells and single-molecule spectroscopy. Here, a new molecular absorption is activated inside a gold plasmonic nanocavity, and found to break selection rules via spin-orbit coupling. Photoluminescence excitation scans reveal absorption from a normally spin-forbidden singlet to trip…
▽ More
Controlling absorption and emission of organic molecules is crucial for efficient light-emitting diodes, organic solar cells and single-molecule spectroscopy. Here, a new molecular absorption is activated inside a gold plasmonic nanocavity, and found to break selection rules via spin-orbit coupling. Photoluminescence excitation scans reveal absorption from a normally spin-forbidden singlet to triplet state transition, while drastically enhancing the emission rate by several thousand fold. The experimental results are supported by density functional theory, revealing the manipulation of molecular absorption by nearby metallic gold atoms.
△ Less
Submitted 11 May, 2020;
originally announced May 2020.
-
Correlation Functions, Mean First Passage Times and the Kemeny Constant
Authors:
Adam Kells,
Edina Rosta,
Alessia Annibale
Abstract:
Markov processes are widely used models for investigating kinetic networks. Here we collate and present a variety of results pertaining to kinetic network models, in a unified framework. The aim is to lay out explicit links between several important quantities commonly studied in the field, including mean first passage times (MFPTs), correlation functions and the Kemeny constant, and highlight som…
▽ More
Markov processes are widely used models for investigating kinetic networks. Here we collate and present a variety of results pertaining to kinetic network models, in a unified framework. The aim is to lay out explicit links between several important quantities commonly studied in the field, including mean first passage times (MFPTs), correlation functions and the Kemeny constant, and highlight some of the subtleties which are often overlooked in the literature, while providing new insights. Results include (i) a simple physical interpretation of the Kemeny constant, (ii) a recipe to infer equilibrium distributions and rate matrices from measurements of MFPTs, potentially useful in applications, including milestoning in molecular dynamics, and (iii) a protocol to reduce the dimensionality of kinetic networks, based on specific requirements that the MFPTs in the coarse-grained system should satisfy. It is proven that this protocol coincides with the one proposed by Hummer and Szabo in [1] and it leads to a variational principle for the Kemeny constant. We hope that this study will serve as a useful reference for readers interested in theoretical aspects of kinetic networks, some of which underpin useful applications, including milestoning and coarse-graining.
△ Less
Submitted 5 November, 2019;
originally announced November 2019.
-
Special Topic: Markov Models of Molecular Kinetics
Authors:
Frank Noé,
Edina Rosta
Abstract:
The Journal of Chemical Physics (JCP) article collection on Markov Models of Molecular Kinetics (MMMK) features recent advances developing and using Markov State Models (MSMs) in atomistic molecular simulations and related applications. This editorial provides a brief overview of the state of the art in the field and relates it to the articles in this JCP collection.
The Journal of Chemical Physics (JCP) article collection on Markov Models of Molecular Kinetics (MMMK) features recent advances developing and using Markov State Models (MSMs) in atomistic molecular simulations and related applications. This editorial provides a brief overview of the state of the art in the field and relates it to the articles in this JCP collection.
△ Less
Submitted 2 November, 2019;
originally announced November 2019.
-
Identification and Analysis of Transition and Metastable Markov States
Authors:
Linda Martini,
Adam Kells,
Gerhard Hummer,
Nicolae-Viorel Buchete,
Edina Rosta
Abstract:
We present a new method that enables the identification and analysis of both transition and metastable conformational states from atomistic or coarse-grained molecular dynamics (MD) trajectories. Our algorithm is presented and studied by using both analytical and actual examples from MD simulations of the helix-forming peptide Ala5, and of a larger system, the epidermal growth factor receptor (EGF…
▽ More
We present a new method that enables the identification and analysis of both transition and metastable conformational states from atomistic or coarse-grained molecular dynamics (MD) trajectories. Our algorithm is presented and studied by using both analytical and actual examples from MD simulations of the helix-forming peptide Ala5, and of a larger system, the epidermal growth factor receptor (EGFR) protein. In all cases, our method identifies automatically the corresponding transition states and metastable conformations in an optimal way, with the input of a set of relevant coordinates, by capturing accurately the intrinsic slowest relaxation rate. Our approach provides a general and easy to implement analysis method that provides unique insight into the molecular mechanism and the rare but crucial rate limiting conformational pathways occurring in complex dynamical systems such as molecular trajectories.
△ Less
Submitted 13 May, 2016;
originally announced May 2016.
-
Statistically optimal analysis of state-discretized trajectory data from multiple thermodynamic states
Authors:
Hao Wu,
Antonia S. J. S. Mey,
Edina Rosta,
Frank Noé
Abstract:
We propose a discrete transition-based reweighting analysis method (dTRAM) for analyzing configuration-space-discretized simulation trajectories produced at different thermodynamic states (temperatures, Hamiltonians, etc.) dTRAM provides maximum-likelihood estimates of stationary quantities (probabilities, free energies, expectation values) at any thermodynamic state. In contrast to the weighted h…
▽ More
We propose a discrete transition-based reweighting analysis method (dTRAM) for analyzing configuration-space-discretized simulation trajectories produced at different thermodynamic states (temperatures, Hamiltonians, etc.) dTRAM provides maximum-likelihood estimates of stationary quantities (probabilities, free energies, expectation values) at any thermodynamic state. In contrast to the weighted histogram analysis method (WHAM), dTRAM does not require data to be sampled from global equilibrium, and can thus produce superior estimates for enhanced sampling data such as parallel/simulated tempering, replica exchange, umbrella sampling, or metadynamics. In addition, dTRAM provides optimal estimates of Markov state models (MSMs) from the discretized state-space trajectories at all thermodynamic states. Under suitable conditions, these MSMs can be used to calculate kinetic quantities (e.g. rates, timescales). In the limit of a single thermodynamic state, dTRAM estimates a maximum likelihood reversible MSM, while in the limit of uncorrelated sampling data, dTRAM is identical to WHAM. dTRAM is thus a generalization to both estimators.
△ Less
Submitted 30 November, 2014; v1 submitted 11 November, 2014;
originally announced November 2014.