A data-driven and topological mapping approach for the a priori prediction of stable molecular crystalline hydrates.

Predictions of the structures of stoichiometric, fractional, or nonstoichiometric hydrates of organic molecular crystals are immensely challenging due to the extensive search space of different water contents, host molecular placements throughout the crystal, and internal molecular conformations. However, the dry frameworks of these hydrates, especially for nonstoichiometric or isostructural dehydrates, can often be predicted from a standard anhydrous crystal structure prediction (CSP) protocol. Inspired by developments in the field of drug binding, we introduce an efficient data-driven and topologically aware approach for predicting organic molecular crystal hydrate structures through a mapping of water positions within the crystal structure. The method does not require a priori specification of water content and can, therefore, predict stoichiometric, fractional, and nonstoichiometric hydrate structures. This approach, which we term a mapping approach for crystal hydrates (MACH), establishes a set of rules for systematic determination of favorable positions for water insertion within predicted or experimental crystal structures based on considerations of the chemical features of local environments and void regions. The proposed approach is tested on hydrates of three pharmaceutically relevant compounds that exhibit diverse crystal packing motifs and void environments characteristic of hydrate structures. Overall, we show that our mapping approach introduces an advance in the efficient performance of hydrate CSP through generation of stable hydrate stoichiometries at low cost and should be considered an integral component for CSP workflows.

An interoperable implementation of collective-variable based enhanced sampling methods in extended phase space within the OpenMM package.

Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. In this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.

Microswimmers under the spotlight: interplay between agents with different levels of activity.

The ability of active matter to assemble into reconfigurable nonequilibrium structures has drawn considerable interest in recent years. We investigate how active fluids respond to spatial light patterns through simulations and experiments on light-activated self-propelled colloidal particles. We examine the processes of inverse templated assembly, which involves creating a region without active particles through a bright pattern, and templated assembly, which promotes the formation of dense particle regions through a dark pattern. We identify scaling relations for the characteristic times for both processes that quantify the interplay between the dimension of the applied pattern and the intrinsic properties of the active fluid. We also explore the assembly mechanism and dynamics of large clusters and show how assembly and inverse assembly can be combined to create any arbitrarily complex template. In addition to providing protocols for templated assembly via light patterning, our results demonstrate how the local packing fraction can be fine-tuned by modulation of the light intensity. The protocol so obtained exceeds the capabilities of conventional assembly strategies, in which packing fraction is dictated by thermodynamics, and opens the door to arbitrarily precise and programmable nonequilibrium assembly strategies in active matter.

Deep Eutectic Solvents: A Review of Fundamentals and Applications.

Deep eutectic solvents (DESs) are an emerging class of mixtures characterized by significant depressions in melting points compared to those of the neat constituent components. These materials are promising for applications as inexpensive "designer" solvents exhibiting a host of tunable physicochemical properties. A detailed review of the current literature reveals the lack of predictive understanding of the microscopic mechanisms that govern the structure-property relationships in this class of solvents. Complex hydrogen bonding is postulated as the root cause of their melting point depressions and physicochemical properties; to understand these hydrogen bonded networks, it is imperative to study these systems as dynamic entities using both simulations and experiments. This review emphasizes recent research efforts in order to elucidate the next steps needed to develop a fundamental framework needed for a deeper understanding of DESs. It covers recent developments in DES research, frames outstanding scientific questions, and identifies promising research thrusts aligned with the advancement of the field toward predictive models and fundamental understanding of these solvents.

An exploration of machine learning models for the determination of reaction coordinates associated with conformational transitions.

Determining collective variables (CVs) for conformational transitions is crucial to understanding their dynamics and targeting them in enhanced sampling simulations. Often, CVs are proposed based on intuition or prior knowledge of a system. However, the problem of systematically determining a proper reaction coordinate (RC) for a specific process in terms of a set of putative CVs can be achieved using committor analysis (CA). Identifying essential degrees of freedom that govern such transitions using CA remains elusive because of the high dimensionality of the conformational space. Various schemes exist to leverage the power of machine learning (ML) to extract an RC from CA. Here, we extend these studies and compare the ability of 17 different ML schemes to identify accurate RCs associated with conformational transitions. We tested these methods on an alanine dipeptide in vacuum and on a sarcosine dipeptoid in an implicit solvent. Our comparison revealed that the light gradient boosting machine method outperforms other methods. In order to extract key features from the models, we employed Shapley Additive exPlanations analysis and compared its interpretation with the "feature importance" approach. For the alanine dipeptide, our methodology identifies ϕ and θ dihedrals as essential degrees of freedom in the C7ax to C7eq transition. For the sarcosine dipeptoid system, the dihedrals ψ and ω are the most important for the cisαD to transαD transition. We further argue that analysis of the full dynamical pathway, and not just endpoint states, is essential for identifying key degrees of freedom governing transitions.

A Green's Function Approach for Determining Surface Induced Broadening and Shifting of Molecular Energy Levels.

Upon adsorption of a molecule onto a surface, the molecular energy levels (MELs) broaden and change their alignment. This phenomenon directly affects electron transfer across the interface and is, therefore, a fundamental observable that influences electrochemical device performance. Here, we propose a rigorous parameter-free framework, built upon the theoretical construct of Green's functions, for studying the interface between a molecule and a bulk surface and its effect on MELs. The method extends beyond the usual wide-band limit approximation, and its generality allows its use with any level of electronic structure theory. We demonstrate its ability to predict the broadening and shifting of MELs as a function of intramolecular coupling, molecule/surface coupling, and the surface density of states for a molecule with two MELs adsorbed on a one-dimensional model metal surface. The new approach could help provide guidelines for the design and experimental characterization of electrochemical devices with optimal electron transport.

Imaginary-time open-chain path-integral approach for two-state time correlation functions and applications in charge transfer.

Quantum time correlation functions (TCFs) involving two states are important for describing nonadiabatic dynamical processes such as charge transfer (CT). Based on a previous single-state method, we propose an imaginary-time open-chain path-integral (OCPI) approach for evaluating the two-state symmetrized TCFs. Expressing the forward and backward propagation on different electronic potential energy surfaces as a complex-time path integral, we then transform the path variables to average and difference variables such that the integration over the difference variables up to the second order can be performed analytically. The resulting expression for the symmetrized TCF is equivalent to sampling the open-chain configurations in an effective potential that corresponds to the average surface. Using importance sampling over the extended OCPI space via open path-integral molecular dynamics, we tested the resulting path-integral approximation by calculating the Fermi's golden rule CT rate constant within a widely used spin-boson model. Comparing with the real-time linearized semiclassical method and analytical result, we show that the imaginary-time OCPI provides an accurate two-state symmetrized TCF and rate constant in the typical turnover region. It is shown that the first bead of the open chain corresponds to physical zero-time and that the endpoint bead corresponds to final time t; oscillations of the end-to-end distance perfectly match the nuclear mode frequency. The two-state OCPI scheme is seen to capture the tested model's electronic quantum coherence and nuclear quantum effects accurately.

Imidacloprid Crystal Polymorphs for Disease Vector Control and Pollinator Protection.

Imidacloprid, the world's leading insecticide, has been approved recently for controlling infectious disease vectors; yet, in agricultural settings, it has been implicated in the frightening decline of pollinators. This argues for strategies that sharply reduce the environmental impact of imidacloprid. When used as a contact insecticide, the effectiveness of imidacloprid relies on physical contact between its crystal surfaces and insect tarsi. Herein, seven new imidacloprid crystal polymorphs are reported, adding to two known forms. Anticipating that insect uptake of imidacloprid molecules would depend on the respective free energies of crystal polymorph surfaces, measurements of insect knockdown times for the metastable crystal forms were as much as nine times faster acting than the commercial form against Aedes, Anopheles, and Culex mosquitoes as well as Drosophila (fruit flies). These results suggest that replacement of commercially available imidacloprid crystals (a.k.a. Form I) in space-spraying with any one of three new polymorphs, Forms IV, VI, IX, would suppress vector-borne disease transmission while reducing environmental exposure and harm to nontarget organisms.

Neural-Network-Based Path Collective Variables for Enhanced Sampling of Phase Transformations.

The investigation of the microscopic processes underlying structural phase transformations in solids is extremely challenging for both simulation and experiment. Atomistic simulations of solid-solid phase transitions require extensive sampling of the corresponding high-dimensional and often rugged energy landscape. Here, we propose a rigorous construction of a 1D path collective variable that is used in combination with enhanced sampling techniques for efficient exploration of the transformation mechanisms. The path collective variable is defined in a space spanned by global classifiers that are derived from local structural units. A reliable identification of the local structural environments is achieved by employing a neural-network-based classification scheme. The proposed path collective variable is generally applicable and enables the investigation of both transformation mechanisms and kinetics.

Unified Efficient Thermostat Scheme for the Canonical Ensemble with Holonomic or Isokinetic Constraints via Molecular Dynamics.

We have recently proposed a new unified theoretical scheme (the "middle" scheme) for thermostat algorithms for efficient and accurate configurational sampling of the canonical ensemble. In this paper, we extend the "middle" scheme to molecular dynamics algorithms for configurational sampling in systems subject to constraints. Holonomic constraints and isokinetic constraints are used for demonstration. Numerical examples indicate that the "middle" scheme presents a promising approach to calculate configuration-dependent thermodynamic properties and their thermal fluctuations.