Understanding the molecular-level behavior of ionic liquids (ILs) at IL-solid interfaces is of fundamental importance with respect to their application in, for example, electrochemical systems and electronic devices. Using a model system, consisting of an imidazolium-based IL ([C2Mim][NTf2]) in contact with a sapphire substrate, we have approached this problem using a complementary combination of high-resolution X-ray reflectivity measurements and atomistic molecular dynamics (MD) simulations. Our strategy enabled us to compare experimental and theoretically calculated reflectivities in a direct manner, thereby critically assessing the applicability of several force-field variants. On the other hand, using the best-matching MD description, we are able to describe the nature of the model IL-solid interface in appreciable detail. More specifically, we find that characteristic interactions between the surface hydroxyl groups and donor and acceptor sites on the IL constituents have a dominant role in inducing a multidimensional layering profile of the cations and anions.

No abstract

The active sites of the (6-4) photolyases contain two conserved histidine residues, which, in the Drosophila melanogaster enzyme, correspond to His365 and His369. While there are nine combinations in which the three possible protonation states of the two histidines (with protons on Nδ (HID), Nε (HIE), or both Nδ and Nε (HIP)) can be paired, there is presently no consensus as to which of these states is present, let alone mechanistically relevant. EPR hyperfine couplings for selected protons of the FADH(•) radical have previously been used to address this issue. Our QM/MM calculations show, however, that the experimental couplings are equally well reproduced by each of the nine combinations. Since the EPR results seemingly cannot be used to unequivocally assign the protonation states, the pKa values of the two histidines were calculated using the popular PROPKA, H++, and APBS approaches, in various environments and for several lesions. These techniques consistently indicate that, at pH = 7, both His365 and His369 should be neutral, although His369 is found to be more prone to becoming protonated. In a comparative approach, a series of molecular dynamics simulations was performed with all nine combinations, employing various reference crystal structures and different oxidation states of the FAD cofactor. The overall result of this approach is in agreement with our pKa results. Consequently, although the introduction of the reduced cofactor results in an increased stability for selected protonated states, particularly the His365═HID and His369═HIP combination, the neutral combination His365═HID and His365═HIE stands out as the most relevant state for the activity of the enzyme.

HypothesisA prototypical modelling approach is required for a full characterisation of the static and equilibrium dynamical properties of confined ionic liquids (ILs), in order to gain predictive power of properties that are difficult to extract from experiments. Such a protocol needs to be constructed by benchmarking molecular dynamics simulations against available experiments. SimulationsWe perform an in-depth study of [C 2 Mim] + [NTf 2 ]in bulk, at the vacuum and at hydroxylated alumina surface. Using the charge methods CHelpG, RESP-HF and RESP-B3LYP with charge scaling factors 1.0, 0.9 and 0.85, we search for an optimum non-polarizable force field by benchmarking against self-diffusion coefficients, surface tension, X-ray reflectivity data, and structural data. FindingsBenchmarking, which relies on establishing the significance of an appropriate size of the model systems and the length of the simulations, yields RESP-HF/0.9 as the best suited force field for this IL overall. A complete and accurate characterisation of the spatially-dependent internal configurational space and orientation of IL molecules relative to the solid and vacuum interfaces is obtained. Furthermore, the density and mobility of IL ions in the plane parallel and normal to the interfaces is evaluated and the correlation between the stratification and dynamics in the interfacial layers is detectable deep into the films.