Characterizing Microheterogeneity in Liquid Mixtures via Local Density Fluctuations.

We present a novel approach to characterize and quantify microheterogeneity and microphase separation in computer simulations of complex liquid mixtures. Our post-processing method is based on local density fluctuations of the different constituents in sampling spheres of varying size. It can be easily applied to both molecular dynamics (MD) and Monte Carlo (MC) simulations, including periodic boundary conditions. Multidimensional correlation of the density distributions yields a clear picture of the domain formation due to the subtle balance of different interactions. We apply our approach to the example of force field molecular dynamics simulations of imidazolium-based ionic liquids with different side chain lengths at different temperatures, namely 1-ethyl-3-methylimidazolium chloride, 1-hexyl-3-methylimidazolium chloride, and 1-decyl-3-methylimidazolium chloride, which are known to form distinct liquid domains. We put the results into the context of existing microheterogeneity analyses and demonstrate the advantages and sensitivity of our novel method. Furthermore, we show how to estimate the configuration entropy from our analysis, and we investigate voids in the system. The analysis has been implemented into our program package TRAVIS and is thus available as free software.

BILFF: All-Atom Force Field for Modeling Triazolium- and Benzoate-Based Ionic Liquids.

We present an extension of our previously developed all-atom force field BILFF (Bio-polymers in Ionic Liquids Force Field) to three different ionic liquids: 1-ethyl-3-methyl-1,2,3-triazolium acetate ([EMTr][OAc]), 1-ethyl-3-methyl-1,2,3-triazolium benzoate ([EMTr][OBz]), and 1-ethyl-3-methylimidazolium benzoate ([EMIm][OBz]). These ionic liquids are of practical importance as they have the ability to dissolve significant amounts of cellulose even at room temperature. Our force field is optimized to accurately reproduce the strong hydrogen bonding in the system with nearly quantum chemical accuracy. A very good agreement between the microstructure of the quantum chemical simulations over a wide temperature range and experimental density data with the results of BILFF were observed. Non-trivial effects, such as the solvation shell structure and π-π stacking of the cations, are also accurately reproduced. Our force field enables accurate simulations of larger systems, such as solvated cellulose in different (aqueous) ionic liquids, and is the first to present the optimized parameters for mixtures of these solvents and water.

Ultrafast Proton Transfer Pathways Mediated by Amphoteric Imidazole.

Imidazole, being an amphoteric molecule, can act both as an acid and as a base. This property enables imidazole, as an essential building block, to effectively facilitate proton transport in high-temperature proton exchange membrane fuel cells and in proton channel transmembrane proteins, enabling those systems to exhibit high energy conversion yields and optimal biological function. We explore the amphoteric properties of imidazole by following the proton transfer exchange reaction dynamics with the bifunctional photoacid 7-hydroxyquinoline (7HQ). We show with ultrafast ultraviolet-mid-infrared pump-probe spectroscopy how for imidazole, in contrast to expectations based on textbook knowledge of acid-base reactivity, the preferential reaction pathway is that of an initial proton transfer from 7HQ to imidazole, and only at a later stage a transfer from imidazole to 7HQ, completing the 7HQ tautomerization reaction. An assessment of the molecular distribution functions and first-principles calculations of proton transfer reaction barriers reveal the underlying reasons for our observations.

A force field for bio-polymers in ionic liquids (BILFF) - part 2: cellulose in [EMIm][OAc]/water mixtures.

We present the extension of our force field BILFF (Bio-Polymers in Ionic Liquids Force Field) to the bio-polymer cellulose. We already published BILFF parameters for mixtures of ionic liquid 1-ethyl-3-methylimidazolium acetate ([EMIm][OAc]) with water. Our all-atom force field focuses on a quantitative reproduction of the hydrogen bonds in the complex mixture of cellulose, [EMIm]+, [OAc]- and water when compared to reference ab initio molecular dynamics (AIMD) simulations. To enhance the sampling, 50 individual AIMD simulations starting from different initial configurations were performed for cellulose in solvent instead of one long simulation, and the resulting averages were used for force field optimization. All cellulose force field parameters were iteratively adjusted starting from the literature force field of W. Damm et al. We were able to obtain a very good agreement with respect to both the microstructure of the reference AIMD simulations and experimental results such as the system density (even at higher temperatures) and the crystal structure. Our new force field allows performing very long simulations of large systems containing cellulose solvated in (aqueous) [EMIm][OAc] with almost ab initio accuracy.

Cluster Analysis in Liquids: A Novel Tool in TRAVIS.

We present a novel cluster analysis implemented in our open-source software TRAVIS and its application to realistic and complex chemical systems. The underlying algorithm is exclusively based on atom distances. Using a two-dimensional model system, we first introduce different cluster analysis functions and their application to single snapshots and trajectories including periodicity and temporal propagation. Using molecular dynamics simulations of pure water with varying system size, we show that our cluster analysis is size-independent. Furthermore, we observe a similar clustering behavior of pure water in classical and ab initio molecular dynamics simulations, showing that our cluster analysis is universal. In order to emphasize the application to more complex systems and mixtures, we additionally apply the cluster analysis to ab initio molecular dynamics simulations of the [C2C1Im][OAc] ionic liquid and its mixture with water. Using that, we show that our cluster analysis is able to analyze the clustering of the individual components in a mixture as well as the clustering of the ionic liquid with water.

Conformer Weighting and Differently Sized Cluster Weighting for Nicotine and Its Phosphorus Derivatives.

Weighting methods applied to systems with many conformers have been broadly employed to calculate thermodynamic properties, structural characteristics, and populations. To better understand and test the sensitivity of conventional weighting methods, the conformational distributions of nicotine and its phosphorus-substituted derivatives are investigated. The weighting schemes used for this are all based on Boltzmann statistics. Classical Boltzmann factors based on the electronic energy and the Gibbs free energy are calculated at different quantum chemical levels of theory and compared to cluster weights obtained by the quantum cluster equilibrium method. Furthermore, the influence of the modified rigid-rotor-harmonic-oscillator (mRRHO) approximation on the cluster weights is investigated. The substitution of the nitrogen atom in the methylpyrrolidine ring by a phosphorus atom results in more monomer conformers and clusters being populated. The conformational distribution of the monomers remained stable at different levels of theory and weighting methods. However, going to dimers and trimers, we observe a significant influence of the level of theory, weighting method, and mRRHO cutoff on the populations of these clusters. We show that mRRHO cutoff values of 50 and 100 cm-1 yield similar results, which is why 50 cm-1 is recommended as a robust choice. Furthermore, we observe that the global minimum for ΔE0 and ΔG varies in a few cases and that the global minimum is not always the dominantly occupied structure.

Raman Optical Activity of N-Acetyl-L-Cysteine in Water and in Methanol: The "Clusters-in-a-Liquid" Model and ab Initio Molecular Dynamics Simulations.

Raman and Raman Optical Activity (ROA) spectra of N-acetyl-L-cysteine (NALC), a flexible chiral molecule, were measured in water and in methanol to evaluate the solvent effects. Two different solvation approaches, that is, the DFT based "clusters-in-a-liquid" solvent model and the ab initio molecular dynamics (AIMD) simulations, were applied to simulate the Raman and ROA spectra. Systematic conformational searches were carried out using a recently developed conformational searching tool, CREST, with the inclusion of polarizable continuum model of water and of methanol. The CREST candidates of NALC and the NALC-solvent complexes were re-optimized and their Raman and ROA simulations were done at the B3LYP-D3BJ/def2-TZVP and the B3LYP-aug-cc-pVDZ//cc-pVTZ levels. Also, AIMD simulations, which includes some anharmonic effects and all intermolecular interactions in solution, were performed. By empirically weighting the computed Raman and ROA spectra of each conformer, good agreements with the experimental data were achieved with both approaches, while AIMD offered some improvements in the carbonyl and in the low wavenumber regions over the static DFT approach. The pros and cons of these two different approaches for accounting the solvent effects on Raman and ROA of this flexible chiral system will also be discussed.

Transferable Deep Learning Potential Reveals Intermediate-Range Ordering Effects in LiF-NaF-ZrF4 Molten Salt.

LiF-NaF-ZrF4 multicomponent molten salts are promising candidate coolants for advanced clean energy systems owing to their desirable thermophysical and transport properties. However, the complex structures enabling these properties, and their dependence on composition, is scarcely quantified due to limitations in simulating and interpreting experimental spectra of highly disordered, intermediate-ranged structures. Specifically, size-limited ab initio simulations and accuracy-limited classical models used in the past are unable to capture a wide range of fluctuating motifs found in the extended heterogeneous structures of liquid salt. This greatly inhibits our ability to design tailored compositions and materials. Here, accurate, efficient, and transferable machine learning potentials are used to predict structures far beyond the first coordination shell in LiF-NaF-ZrF4. Neural networks trained at only eutectic compositions with 29% and 37% ZrF4 are shown to accurately simulate a wide range of compositions (11-40% ZrF4) with dramatically different coordination chemistries, while showing a remarkable agreement with theoretical and experimental Raman spectra. The theoretical Raman calculations further uncovered the previously unseen shift and flattening of bending band at ∼250 cm-1 which validated the simulated extended-range structures as observed in compositions with higher than 29% ZrF4 content. In such cases, machine learning-based simulations capable of accessing larger time and length scales (beyond 17 Å) were critical for accurately predicting both structure and ionic diffusivities.

Characterization of the Diagnostic Performance of a Novel COVID-19 PETIA in Comparison to Four Routine N-, S- and RBD-Antigen Based Immunoassays.

In 2019, a novel coronavirus emerged in Wuhan in the province of Hubei, China. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) quickly spread across the globe, causing the neoteric COVID-19 pandemic. SARS-CoV-2 is commonly transmitted by droplet infection and aerosols when coughing or sneezing, as well as high-risk exposures to infected individuals by face-to-face contact without protective gear. To date, a broad variety of techniques have emerged to assess and quantify the specific antibody response of a patient towards a SARS-CoV-2 infection. Here, we report the first comprehensive comparison of five different assay systems: Enzyme-Linked Immunosorbent Assay (ELISA), Chemiluminescence Immunoassay (CLIA), Electro-Chemiluminescence Immunoassay (ECLIA), and a new Particle-Enhanced Turbidimetric Immunoassay (PETIA) for SARS-CoV-2. Furthermore, we also evaluated the suitability of N-, S1- and RBD-antigens for quantifying the SARS-CoV-2 specific immune response. Linearity and precision, overall sensitivity and specificity of the assays, stability of samples, and cross-reactivity of general viral responses, as well as common coronaviruses, were assessed. Moreover, the reactivity of all tests to seroconversion and different sample matrices was quantified. All five assays showed good overall agreement, with 76% and 87% similarity for negative and positive samples, respectively. In conclusion, all evaluated methods showed a high consistency of results and suitability for the robust quantification of the SARS-CoV-2-derived immune response.

Liquid structure and dynamics in the choline acetate:urea 1:2 deep eutectic solvent.

We report on the thermodynamic, structural, and dynamic properties of a recently proposed deep eutectic solvent, formed by choline acetate (ChAc) and urea (U) at the stoichiometric ratio 1:2, hereinafter indicated as ChAc:U. Although the crystalline phase melts at 36-38 °C depending on the heating rate, ChAc:U can be easily supercooled at sub-ambient conditions, thus maintaining at the liquid state, with a glass-liquid transition at about -50 °C. Synchrotron high energy x-ray scattering experiments provide the experimental data for supporting a reverse Monte Carlo analysis to extract structural information at the atomistic level. This exploration of the liquid structure of ChAc:U reveals the major role played by hydrogen bonding in determining interspecies correlations: both acetate and urea are strong hydrogen bond acceptor sites, while both choline hydroxyl and urea act as HB donors. All ChAc:U moieties are involved in mutual interactions, with acetate and urea strongly interacting through hydrogen bonding, while choline being mostly involved in van der Waals mediated interactions. Such a structural situation is mirrored by the dynamic evidences obtained by means of 1H nuclear magnetic resonance techniques, which show how urea and acetate species experience higher translational activation energy than choline, fingerprinting their stronger commitments into the extended hydrogen bonding network established in ChAc:U.

Unraveling Local Structure of Molten Salts via X-ray Scattering, Raman Spectroscopy, and Ab Initio Molecular Dynamics.

In this work, we resolve a long-standing issue concerning the local structure of molten MgCl2 by employing a multimodal approach, including X-ray scattering and Raman spectroscopy, along with the theoretical modeling of the experimental spectra based on ab initio molecular dynamics (AIMD) simulations utilizing several density functional theory (DFT) methods. We demonstrate the reliability of AIMD simulations in achieving excellent agreement between the experimental and simulated spectra for MgCl2 and 50 mol % MgCl2 + 50 mol % KCl, and ZnCl2, thus allowing structural insights not directly available from experiment alone. A thorough computational analysis using five DFT methods provides a convergent view that octahedrally coordinated magnesium in pure MgCl2 upon melting preferentially coordinates with five chloride anions to form distorted square pyramidal polyhedra that are connected via corners and to a lesser degree via edges. This is contrasted with the results for ZnCl2, which does not change its tetrahedral coordination on melting. Although the five-coordinate MgCl53- complex was not considered in the early literature, together with an increasing tendency to form a tetrahedrally coordinated complex with decreasing the MgCl2 content in the mixture with alkali metal chloride systems, current work reconciles the results of most previous seemingly contradictory experimental studies.

Optimized Atomic Partial Charges and Radii Defined by Radical Voronoi Tessellation of Bulk Phase Simulations.

We present a novel method for the computation of well-defined optimized atomic partial charges and radii from the total electron density. Our method is based on a two-step radical Voronoi tessellation of the (possibly periodic) system and subsequent integration of the total electron density within each Voronoi cell. First, the total electron density is partitioned into the contributions of each molecule, and subsequently the electron density within each molecule is assigned to the individual atoms using a second set of atomic radii for the radical Voronoi tessellation. The radii are optimized on-the-fly to minimize the fluctuation (variance) of molecular and atomic charges. Therefore, our method is completely free of empirical parameters. As a by-product, two sets of optimized atomic radii are produced in each run, which take into account many specific properties of the system investigated. The application of an on-the-fly interpolation scheme reduces discretization noise in the Voronoi integration. The approach is particularly well suited for the calculation of partial charges in periodic bulk phase systems. We apply the method to five exemplary liquid phase simulations and show how the optimized charges can help to understand the interactions in the systems. Well-known effects such as reduced ion charges below unity in ionic liquid systems are correctly predicted without any tuning, empiricism, or rescaling. We show that the basis set dependence of our method is very small. Only the total electron density is evaluated, and thus, the approach can be combined with any electronic structure method that provides volumetric total electron densities-it is not limited to Hartree-Fock or density functional theory (DFT). We have implemented the method into our open-source software tool TRAVIS.

Switching between Proton Vacancy and Excess Proton Transfer Pathways in the Reaction between 7-Hydroxyquinoline and Formate.

Bifunctional or amphoteric photoacids simultaneously present donor (acidic) and acceptor (basic) properties making them useful tools to analyze proton transfer reactions. In protic solvents, the proton exchange between the acid and the base is controlled by the acidity or basicity strength and typically occurs on two different pathways known as protolysis and hydrolysis. We report here how the addition of a formate base will alter the relative importance of the possible reaction pathways of the bifunctional photoacid 7-hydroxyquinoline (7HQ), which has been recently understood to predominantly involve a hydroxide/methoxide transport mechanism between the basic proton-accepting quinoline nitrogen site toward the proton-donating OH group with a time constant of 360 ps in deuterated methanol (CD3OD). We follow the reaction dynamics by probing the IR-active marker modes of the different charged forms of photoexcited 7HQ, and of formic acid (HCOOD) in CD3OD solution. A comparison of the transient IR spectra as a function of formate concentration, and classical molecular dynamics simulations enables us to identify distinct contributions of "tight" (meaning "contact") and "loose" (i.e., "solvent-separated") 7HQ-formate reaction pairs in our data. Our results suggest that depending on the orientation of the OH group with respect to the quinoline aromatic ring system, the presence of the formate molecule in a proton relay pathway facilitates a net proton transfer from the proton-donating OH group of 7HQ-N* via the methanol/formate bridge toward the quinoline N site.

Efficient EOM-CC-based Protocol for the Calculation of Electron Affinity of Solvated Nucleobases: Uracil as a Case Study.

We present an explicit solvation protocol for the calculation of electron affinity values of the solvated nucleobases. The protocol uses a quantum mechanics/molecular mechanics (QM/MM) approach based on the newly implemented domain-based pair natural orbital EOM-CCSD (equation-of-motion coupled-cluster single-double) method. The stability of the solvated nucleobase anion is sensitive to the local distribution of the water molecules around the nucleobase, and the calculated electron affinity values converge slowly with respect to the number of snapshots and the size of the water box. The use of nonpolarizable water molecules leads to an overestimation of the electron affinity and makes the result sensitive to the size of the QM region in the QM/MM calculation. The electron affinity values, although sensitive to the size of the basis set, lead to an almost constant blue shift of the electron affinity upon the increase in the basis set. The present protocol allows for a controllable description of the various parameters affecting the electron affinity value, and the calculated adiabatic electron affinity values are in excellent agreement with experimental results.

A force field for bio-polymers in ionic liquids (BILFF) - part 1: [EMIm][OAc]/water mixtures.

We present BILFF, a novel force field for bio-polymers in ionic liquids. In the first part of our study, we introduce optimized force field parameters for mixtures of the ionic liquid (IL) 1-ethyl-3-methylimidazolium acetate ([EMIm][OAc]) with water. This imidazolium-based IL is of particular practical importance as it can dissolve significant amounts of cellulose even at room temperature. An understanding of this dissolution process via molecular dynamics simulations requires a quantitative description of the microscopic structure and the strong hydrogen bonds with a method able of simulating at least several dozen nanoseconds, which is the main aim of our novel force field. To reach this goal, we optimize the force field parameters to reproduce radial, spatial, and combined distribution functions, hydrogen bond lifetimes, diffusion coefficients, and several other quantities from reference ab initio molecular dynamics (AIMD) simulations. Non-trivial effects such as dispersion interactions between the side chains and π-π stacking of the cations are reproduced very well. We further validate the force field by comparison to experimental data such as thermal expansion coefficients, bulk modulus, and density at different temperatures, which yields good agreement and correct trends. No other force field with optimized parameters for mixtures of [EMIm][OAc] and water has been presented in the literature yet. Optimized force field parameters for cellulose and other ILs will be published in upcoming articles.

Exploring Free Energy Profiles of Enantioselective Organocatalytic Aldol Reactions under Full Solvent Influence.

We present a computational study on the enantioselectivity of organocatalytic proline-catalyzed aldol reactions between aldehydes in dimethylformamide (DMF). To explore the free energy surface of the reaction, we apply two-dimensional metadynamics on top of ab initio molecular dynamics (AIMD) simulations with explicit solvent description on the DFT level of theory. We avoid unwanted side reactions by utilizing our newly developed hybrid AIMD (HyAIMD) simulation scheme, which adds a simple force field to the AIMD simulation to prevent unwanted bond breaking and formation. Our condensed phase simulation results are able to nicely reproduce the experimental findings, including the main stereoisomer that is formed, and give a correct qualitative prediction of the change in syn:anti product ratio with different substituents. Furthermore, we give a microscopic explanation for the selectivity. We show that both the explicit description of the solvent and the inclusion of entropic effects are vital to a good outcome-metadynamics simulations in vacuum and static nudged elastic band (NEB) calculations yield significantly worse predictions when compared to the experiment. The approach described here can be applied to a plethora of other enantioselective or organocatalytic reactions, enabling us to tune the catalyst or determine the solvent with the highest stereoselectivity.

Characterization of Aqueous Lower-Polarity Solvation Shells Around Amphiphilic 2,2,6,6-Tetramethylpiperidine-1-oxyl Radicals in Water.

Solvation of the amphiphilic nitroxide radical 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) and hydrophilic 4-oxo-2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPONE) in water and tetrahydrofuran (THF) is studied in detail. The existence of pure water shells enclosing TEMPO in an aqueous solution that leads to significantly reduced local polarity at the nitroxide moiety is shown with multifrequency electron paramagnetic resonance (EPR) spectroscopy at X- and Q-bands as well as spectral simulations. These aqueous lower-polarity solvation shells (ALPSS) offer TEMPO a local polarity that is similar to that in organic solvents like THF. Furthermore, using double electron-electron resonance spectroscopy, local enrichment and inhomogeneous distribution without direct molecular encounters of dissolved TEMPO in water are found that can be correlated with potentially attractive interactions mediated through ALPSS. However, no local enrichment of TEMPO is found in organic solvents such as THF. In contrast to TEMPO, the structurally very similar nitroxide radical TEMPONE shows no ALPSS encapsulation behavior with water molecules in aqueous solutions. Ensemble-averaging methods such as dynamic light scattering and electrospray ionization mass spectrometry substantiate the EPR spectroscopically obtained results of ALPSS-encased TEMPO and attractive interactions between them, leading to a higher local concentration. Furthermore, force field molecular dynamics simulations and metadynamics deliver support for our conclusions.

Dissolving Cellulose in 1,2,3-Triazolium- and Imidazolium-Based Ionic Liquids with Aromatic Anions.

We present 1,2,3-triazolium- and imidazolium-based ionic liquids (ILs) with aromatic anions as a new class of cellulose solvents. The two anions in our study, benzoate and salicylate, possess a lower basicity when compared to acetate and therefore should lead to a lower amount of N-heterocyclic carbenes (NHCs) in the ILs. We characterize their physicochemical properties and find that all of them are liquids at room temperature. By applying force field molecular dynamics (MD) simulations, we investigate the structure and dynamics of the liquids and find strong and long-lived hydrogen bonds, as well as significant π-π stacking between the aromatic anion and cation. Our ILs dissolve up to 8.5 wt.-% cellulose. Via NMR spectroscopy of the solution, we rule out chain degradation or derivatization, even after several weeks at elevated temperature. Based on our MD simulations, we estimate the enthalpy of solvation and derive a simple model for semi-quantitative prediction of cellulose solubility in ILs. With the help of Sankey diagrams, we illustrate the hydrogen bond network topology of the solutions, which is characterized by competing hydrogen bond donors and acceptors. The hydrogen bonds between cellulose and the anions possess average lifetimes in the nanosecond range, which is longer than found in common pure ILs.

Exploring non-equilibrium molecular dynamics of mobile protons in the solid acid CsH2PO4 at the micrometer and microsecond scale.

We explicitly compute the non-equilibrium molecular dynamics of protons in the solid acid CsH2PO4 on the micrometer length scale via a multiscale Markov model: The molecular dynamics/matrix propagation (MDM) method. Within the MDM approach, the proton dynamics information of an entire molecular dynamics simulation can be condensed into a single M × M matrix (M is the number of oxygen atoms in the simulated system). Due to this drastic reduction in the complexity, we demonstrate how to increase the length and time scales in order to enable the simulation of inhomogeneities of CsH2PO4 systems at the nanometer scale. We incorporate explicit correlation of protonation dynamics with the protonation state of the neighboring proton sites and illustrate that this modification conserves the Markov character of the MDM method. We show that atomistic features such as the mean square displacement and the diffusion coefficient of the protons can be computed quantitatively from the matrix representation. Furthermore, we demonstrate the application potential of the scheme by computing the explicit dynamics of a non-equilibrium process in an 8 µm CsH2PO4 system during 5 ms.

Dynamical matrix propagator scheme for large-scale proton dynamics simulations.

We derive a matrix formalism for the simulation of long range proton dynamics for extended systems and timescales. On the basis of an ab initio molecular dynamics simulation, we construct a Markov chain, which allows us to store the entire proton dynamics in an M × M transition matrix (where M is the number of oxygen atoms). In this article, we start from common topology features of the hydrogen bond network of good proton conductors and utilize them as constituent constraints of our dynamic model. We present a thorough mathematical derivation of our approach and verify its uniqueness and correct asymptotic behavior. We propagate the proton distribution by means of transition matrices, which contain kinetic data from both ultra-short (sub-ps) and intermediate (ps) timescales. This concept allows us to keep the most relevant features from the microscopic level while effectively reaching larger time and length scales. We demonstrate the applicability of the transition matrices for the description of proton conduction trends in proton exchange membrane materials.