Effects of the bond polarity on the structural and dynamical properties of silica-like liquids

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

Molecular dynamics study of water channels in natural and synthetic amyloid-β fibrils

We compared all-atom explicit solvent molecular dynamics simulations of three types of Aβ(1-40) fibrils: brain-seeded fibrils (2M4J, with a threefold axial symmetry) and the other two, all-synthetic fibril polymorphs (2LMN and 2LMP, made under different fibrillization conditions). Fibril models were constructed using either a finite or an infinite number of layers made using periodic images. These studies yielded four conclusions. First, finite fibrils tend to unravel in a manner reminiscent of fibril dissolution, while infinite fibrils were more stable during simulations. Second, salt bridges in these fibrils remained stable in those fibrils that contained them initially, and those without salt bridges did not develop them over the time course of the simulations. Third, all fibrils tended to develop a "stagger" or register shift of β-strands along the fibril axis. Fourth and most importantly, the brain-seeded, 2M4J, infinite fibrils allowed bidirectional transport of water in and out of the central longitudinal core of the fibril by rapidly developing gaps at the fibril vertices. 2LMP fibrils also showed this behavior, although to a lesser extent. The diffusion of water molecules in the fibril core region involved two dynamical states: a localized state and directed diffusion in the presence of obstacles. These observations provided support for the hypothesis that Aβ fibrils could act as nanotubes. At least some Aβ oligomers resembled fibrils structurally in having parallel, in-register β-sheets and a sheet-turn-sheet motif. Thus, our findings could have implications for Aβ cytotoxicity, which may occur through the ability of oligomers to form abnormal water and ion channels in cell membranes.

MD simulations of charged binary mixtures reveal a generic relation between high- and low-temperature behavior

Experimental studies of the glassy slowdown in molecular liquids indicate that the high-temperature activation energy E_∞ of glass-forming liquids is directly related to their glass transition temperature T_g. To further investigate such a possible relation between high- and low-temperature dynamics in glass-forming liquids, we analyze the glassy dynamics of binary mixtures using molecular dynamics simulations. We consider a binary mixture of charged Lennard-Jones particles and vary the partial charges of the particles and, thus, the high-temperature activation energy and the glass transition temperature of the system. Based on previous results, we introduce a phenomenological model describing relaxation times over the whole temperature regime from high temperatures to temperatures well inside the supercooled regime. By investigating the dynamics of both particle species on molecular and diffusive length scales along isochoric and isobaric pathways, we find a quadratic charge dependence of both E_∞ and T_g, resulting in an approximately constant ratio of both quantities independent of the underlying observable, the thermodynamic ensemble, and the particle species, and this result is robust against the actual definition of T_g. This generic relation between the activation energy and the glass transition temperature indicates that high-temperature dynamics and the glassy slowdown are related phenomena, and the knowledge of E_∞ may allow us to approximately predict T_g.

Effects of confinement on supercooled tetrahedral liquids

We use molecular dynamics simulations to ascertain the effects of geometrical restriction on glass-forming tetrahedral liquids. Striving for a broad approach, we study families of waterlike and silicalike liquids, for which we systematically scale the partial charges and, hence, the relevance of the tetrahedral networks. The confined liquids and the confining matrices consist of the same type of particles to avoid disruptive interactions and distorted structures at the interfaces. Spatially resolved analyses show that these neutral confinements still impose static mobility gradients and density correlations on the liquids. We quantify the increasing degree and range of the altered properties upon cooling. For both families of models, common relations describe the confinement effects of all systems with tetrahedral order, while deviations occur for systems with lower polarities and different structures. The observations are rationalized by considering the fact that a pinned wall imprints a static energy landscape to a neighboring liquid. We explore the properties of this landscape based on changes in vibrational motion and structural relaxation and find that typical barrier heights amount to two to three times the activation energy of bulk dynamics. Combining the present and previous results, we predict the evolution of confinement effects down to the glass transition temperature for liquids without fragile-to-strong crossover. In addition, it is found for silicalike liquids that the temperature dependence of dynamic and static correlation lengths from confinement studies is not affected when cooling through fragile-to-strong transitions of the bulk materials, casting doubt on the relevance of these length scales for the glassy slowdown.

Molecular Dynamics Simulations of Water, Silica, and Aqueous Mixtures in Bulk and Confinement

Aqueous systems are omnipresent in nature and technology. They show complex behaviors, which often originate in the existence of hydrogen-bond networks. Prominent examples are the anomalies of water and the non-ideal behaviors of aqueous solutions. The phenomenology becomes even richer when aqueous liquids are subject to confinement. To this day, many properties of water and its mixtures, in particular, under confinement, are not understood. In recent years, molecular dynamics simulations developed into a powerful tool to improve our knowledge in this field. Here, our simulation results for water and aqueous mixtures in the bulk and in various confinements are reviewed and some new simulation data are added to improve our knowledge about the role of interfaces. Moreover, findings for water are compared with results for silica, exploiting that both systems form tetrahedral networks.