The ionic liquid–vacuum outer atomic surface: a low-energy ion scattering study

We have identified elements present in the ionic liquid–vacuum outer atomic surface of 23 ionic liquids using high sensitivity low-energy ion scattering (LEIS), a very surface sensitive technique. We show that the probability of cationic heteroatoms being present at the ionic liquid–vacuum outer atomic surface is very low; we detected imidazolium nitrogen for only one of the 18 imidazolium based ionic liquids investigated, no nitrogen for the two ammonium based ionic liquids and a very small amount of phosphorus for two of the three phosphonium-based ionic liquids. We determine that the anion is always present at the ionic liquid–vacuum outer atomic surface, even for very large cations containing dodecyl alkyl chains or longer; these chains dominate the ionic liquid–vacuum outer atomic surface, but are not sufficiently densely packed to completely cover the anions. We demonstrate the presence of strong hydrogen bond acceptor adsorption sites at the ionic liquid–vacuum outer atomic surface. We demonstrate that the amount of ion present at the ionic liquid–vacuum outer atomic surface can be tuned by varying the size of the other ion; larger cations (or anions) occupy more of the ionic liquid–vacuum outer atomic surface, leaving less room for anions (or cations). By identifying elements present at the ionic liquid–vacuum outer atomic surface, conclusions can be drawn on the orientations of anions nearest the vacuum. We show that for five different anions there is a most probable ion orientation, but other anion orientations also exist, demonstrating the presence of multiple anion orientations. The imidazolium cations nearest to the vacuum also show similar multi-orientation behaviour. This variety of atoms present and therefore ion orientations is expected to be central to controlling surface reactivity. In addition, our results can be used to quantitatively validate simulations of the ionic liquid–vacuum surface at a molecular level. Overall, our studies, in combination with literature data from different techniques and simulations, provide a clear picture of ionic liquid–vacuum outer atomic surfaces

Plate tectonics of virus shell assembly and reorganization in phage φ8, a distant relative of mammalian reoviruses

The hallmark of a virus is its capsid, which harbors the viral genome and is formed from protein subunits, which assemble following precise geometric rules. dsRNA viruses use an unusual protein multiplicity (120 copies) to form their closed capsids. We have determined the atomic structure of the capsid protein (P1) from the dsRNA cystovirus Φ8. In the crystal P1 forms pentamers, very similar in shape to facets of empty procapsids, suggesting an unexpected assembly pathway that proceeds via a pentameric intermediate. Unlike the elongated proteins used by dsRNA mammalian reoviruses, P1 has a compact trapezoid-like shape and a distinct arrangement in the shell, with two near-identical conformers in nonequivalent structural environments. Nevertheless, structural similarity with the analogous protein from the mammalian viruses suggests a common ancestor. The unusual shape of the molecule may facilitate dramatic capsid expansion during phage maturation, allowing P1 to switch interaction interfaces to provide capsid plasticity

Properties of the Liquid-Vapor Interface of Acetone-Water Mixtures. A Computer Simulation and ITIM Analysis Study

Molecular dynamics simulations of the liquid-vapor interface of acetone-water mixtures of different compositions, covering the entire composition range have been performed on the canonical (N, V, T) ensemble at 298 K, using a model combination that excellently describes the mixing properties of these compounds. The properties of the intrinsic liquid surfaces have been analyzed in terms of the Identification of the Truly Interfacial Molecules (ITIM) method. Thus, the composition, width, roughness, and separation of the subsurface molecular layers, as well as self-association, orientation, and dynamics of exchange with the bulk phase of the surface molecules have been analyzed in detail. Our results show that acetone molecules are strongly adsorbed at the liquid surface, and this adsorption extends to several molecular layers. Like molecules in the surface layer are found to form relatively large lateral self-associates. The effect of the vicinity of the vapor phase on a number of properties of the liquid phase vanishes beyond the first molecular layer, with the second subsurface layer already part of the bulk liquid phase in these respects. The orientational preferences of the surface molecules are governed primarily by the dipole-dipole interaction of the neighboring acetone molecules, and hydrogen bonding interaction of the neighboring acetone-water pairs. (Figure Presented). © 2015 American Chemical Society

Floating Patches of HCN at the Surface of Their Aqueous Solutions - Can They Make "HCN World" Plausible?

The liquid/vapor interface of the aqueous solutions of HCN of different concentrations has been investigated using molecular dynamics simulation and intrinsic surface analysis. Although HCN is fully miscible with water, strong interfacial adsorption of HCN is observed at the surface of its aqueous solutions, and, at the liquid surface, the HCN molecules tend to be located even at the outer edge of the surface layer. It turns out that in dilute systems the HCN concentration can be about an order of magnitude larger in the surface layer than in the bulk liquid phase. Furthermore, HCN molecules show a strong lateral self-association behavior at the liquid surface, forming thus floating HCN patches at the surface of their aqueous solutions. Moreover, HCN molecules are staying, on average, an order of magnitude longer at the liquid surface than water molecules, and this behavior is more pronounced at smaller HCN concentrations. Because of this enhanced dynamical stability, the floating HCN patches can provide excellent spots for polymerization of HCN, which can be the key step in the prebiotic synthesis of partially water-soluble adenine. All of these findings make the hypothesis of "HCN world" more plausible

The phase diagram of water at high pressures as obtained by computer simulations of the TIP4P/2005 model: the appearance of a plastic crystal phase

In this work the high pressure region of the phase diagram of water has been studied by computer simulation by using the TIP4P/2005 model of water. Free energy calculations were performed for ices VII and VIII and for the fluid phase to determine the melting curve of these ices. In addition molecular dynamics simulations were performed at high temperatures (440K) observing the spontaneous freezing of the liquid into a solid phase at pressures of about 80000 bar. The analysis of the structure obtained lead to the conclusion that a plastic crystal phase was formed. In the plastic crystal phase the oxygen atoms were arranged forming a body center cubic structure, as in ice VII, but the water molecules were able to rotate almost freely. Free energy calculations were performed for this new phase, and it was found that for TIP4P/2005 this plastic crystal phase is thermodynamically stable with respect to ices VII and VIII for temperatures higher than about 400K, although the precise value depends on the pressure. By using Gibbs Duhem simulations, all coexistence lines were determined, and the phase diagram of the TIP4P/2005 model was obtained, including ices VIII and VII and the new plastic crystal phase. The TIP4P/2005 model is able to describe qualitatively the phase diagram of water. It would be of interest to study if such a plastic crystal phase does indeed exist for real water. The nearly spherical shape of water makes possible the formation of a plastic crystal phase at high temperatures. The formation of a plastic crystal phase at high temperatures (with a bcc arrangements of oxygen atoms) is fast from a kinetic point of view occurring in about 2ns. This is in contrast to the nucleation of ice Ih which requires simulations of the order of hundreds of ns

Computational approaches to understanding reaction outcomes of organic processes in ionic liquids

This review considers how various computational methods have been applied to explain the changes in reaction outcome on moving from a molecular to an ionic liquid solvent. Initially, different conceptual approaches to modelling ionic liquids are discussed, followed by a consideration of the limitations and constraints of these approaches. A series of case studies demonstrating the utility of computational approaches to explain processes in ionic liquids are considered; some of these address the solubility of species in ionic liquids while others examine classes of reaction where the outcome in ionic liquids can be explained through the application of computational approaches. Overall, the utility of computational methods to explain, and potentially predict, the effect of ionic liquids on reaction outcome is demonstrated

On the numerical treatment of dissipative particle dynamics and related systems

We review and compare numerical methods that simultaneously control temperature while preserving the momentum, a family of particle simulation methods commonly used for the modelling of complex fluids and polymers. The class of methods considered includes dissipative particle dynamics (DPD) as well as extended stochastic-dynamics models incorporating a generalized pairwise thermostat scheme in which stochastic forces are eliminated and the coefficient of dissipation is treated as an additional auxiliary variable subject to a feedback (kinetic energy) control mechanism. In the latter case, we consider the addition of a coupling of the auxiliary variable, as in the Nos\'{e}-Hoover-Langevin (NHL) method, with stochastic dynamics to ensure ergodicity, and find that the convergence of ensemble averages is substantially improved. To this end, splitting methods are developed and studied in terms of their thermodynamic accuracy, two-point correlation functions, and convergence. In terms of computational efficiency as measured by the ratio of thermodynamic accuracy to CPU time, we report significant advantages in simulation for the pairwise NHL method compared to popular alternative schemes (up to an 80\% improvement), without degradation of convergence rate. The momentum-conserving thermostat technique described here provides a consistent hydrodynamic model in the low-friction regime, but it will also be of use in both equilibrium and nonequilibrium molecular simulation applications owing to its efficiency and simple numerical implementation