Percolation transition and bimodal density distribution in hydrogen fluoride.

Hydrogen-bond networks in associating fluids can be extremely robust and characterize the topological properties of the liquid phase, as in the case of water, over its whole domain of stability and beyond. Here, we report on molecular dynamics simulations of hydrogen fluoride (HF), one of the strongest hydrogen-bonding molecules. HF has more limited connectivity than water but can still create long, dynamic chains, setting it apart from most other small molecular liquids.

Surface viscosity of liquid interfaces from Green-Kubo relations.

The precise determination of surface transport coefficients at liquid interfaces is critical to an array of processes, ranging from atmospheric chemistry to catalysis. Building on our prior results that highlighted the emergence of a greatly reduced surface viscosity in simple liquids via the dispersion relation of surface excitations [Malgaretti et al., J.

Renormalized charge and dielectric effects in colloidal interactions: a numerical solution of the nonlinear Poisson-Boltzmann equation for unknown boundary conditions.

 The Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, introduced more than 70 years ago, is a hallmark of colloidal particle modeling. For highly charged particles in the dilute regime, it is often supplemented by Alexander's prescription (Alexander et al. in J Chem Phys 80:5776, 1984) for using a renormalized charge.

Local structure of liquid/vapour interfaces approaching the critical point.

Investigating the structure of fluid interfaces at high temperatures is a particularly delicate task that requires effective ways of discriminating liquid from vapour and identifying the location of the liquid phase boundary, thereby allowing to distinguish intrinsic from capillary fluctuations. Several numerical approaches require introducing a coarse-graining length scale, often heuristically chosen to be the molecular size, to determine the location of the liquid phase boundary. Here, we propose an alternative rationale for choosing this coarse-graining length scale; we require the average position of the local liquid phase dividing surface to match its flat, macroscopic counterpart.

Surface viscosity in simple liquids.

The response of Newtonian liquids to small perturbations is usually considered to be fully described by homogeneous transport coefficients like shear and dilatational viscosity. However, the presence of strong density gradients at the liquid/vapor boundary of fluids hints at the possible existence of an inhomogeneous viscosity. Here, we show that a surface viscosity emerges from the collective dynamics of interfacial layers in molecular simulations of simple liquids.

High-density liquid (HDL) adsorption at the supercooled water/vapor interface and its possible relation to the second surface tension inflection point.

We investigate the properties of water along the liquid/vapor coexistence line in the supercooled regime down to the no-man's land. Extensive molecular dynamics simulations of the TIP4P/2005 liquid/vapor interface in the range 198-348 K allow us to locate the second surface tension inflection point with a high accuracy at 283 ± 5 K, close to the temperature of maximum density. This temperature also coincides with the appearance of a density anomaly at the interface known as the apophysis.

Accuracy and performance of the lattice Boltzmann method with 64-bit, 32-bit, and customized 16-bit number formats.

Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here we evaluate the possibility to use even FP16 and posit16 (half) precision for storing fluid populations, while still carrying arithmetic operations in FP32.

Contribution of Different Molecules and Moieties to the Surface Tension in Aqueous Surfactant Solutions. II: Role of the Size and Charge Sign of the Counterions.

Understanding the role of the counterion species in surfactant solutions is a complicated task, made harder by the fact that, experimentally, it is not possible to vary independently bulk and surface quantities. Here, we perform molecular dynamics simulations at constant surface coverage of the liquid/vapor interface of lithium, sodium, potassium, rubidium, and cesium dodecyl sulfate aqueous solutions. We investigate the effect of counterion type and charge sign on the surface tension of the solution, analyzing the contribution of different species and moieties to the lateral pressure profile.

Monolayer Structures of Supramolecular Antagonistic Salt Aggregates.

The speculated presence of monomolecular lamellae of antagonistic salts in oil-water mixtures has left several open questions besides their hypothetical existence, including their microscopic structure and stabilization mechanism. Here, we simulate the spontaneous formation of supramolecular aggregates of the antagonistic salt sodium tetraphenylborate (NaBPh) in water and 3-methylpyridine (3-MP) at the atomistic level. We show that, indeed, the lamellae are formed by a monomolecular layer of the anion, enveloped by 3-MP and hydrated sodium counterions.

Weak scaling of the contact distance between two fluctuating interfaces with system size.

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions.

Single-Particle Dynamics at the Intrinsic Surface of Aqueous Alkali Halide Solutions.

The distribution of ions in the proximity of the liquid-vapor interface of their aqueous solution has been the subject of an intense debate during the last decade. The effects of ionic polarizability have been one of its salient aspects. Much less has been said about the corresponding dynamical properties, which are substantially unexplored.

Capillary interactions between soft capsules protruding through thin fluid films.

When a suspension dries, the suspending fluid evaporates, leaving behind a dry film composed of the suspended particles. During the final stages of drying, the height of the fluid film on the substrate drops below the particle size, inducing local interface deformations that lead to strong capillary interactions among the particles. Although capillary interactions between rigid particles are well studied, much is still to be understood about the behaviour of soft particles and the role of their softness during the final stages of film drying.

Ab initio structure and thermodynamics of the RPBE-D3 water/vapor interface by neural-network molecular dynamics.

Aided by a neural network representation of the density functional theory potential energy landscape of water in the Revised Perdew-Burke-Ernzerhof approximation corrected for dispersion, we calculate several structural and thermodynamic properties of its liquid/vapor interface. The neural network speed allows us to bridge the size and time scale gaps required to sample the properties of water along its liquid/vapor coexistence line with unprecedented precision.

Surface Affinity of Alkali and Halide Ions in Their Aqueous Solution: Insight from Intrinsic Density Analysis.

The surface tension of all aqueous alkali halide solutions is higher than that of pure water. According to the Gibbs adsorption equation, this indicates a net depletion of these ions in the interfacial region. However, simulations and experiments show that large, soft ions, such as I, can accumulate at the liquid/vapor interface.

PDADMAC/PSS Oligoelectrolyte Multilayers: Internal Structure and Hydration Properties at Early Growth Stages from Atomistic Simulations.

We analyze the internal structure and hydration properties of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) oligoelectrolyte multilayers at early stages of their layer-by-layer growth process. Our study is based on large-scale molecular dynamics simulations with atomistic resolution that we presented recently [Sánchez et al., , 15, 9437], in which we produced the first four deposition cycles of a multilayer obtained by alternate exposure of a flat silica substrate to aqueous electrolyte solutions of such polymers at 0.

Atomistic simulation of PDADMAC/PSS oligoelectrolyte multilayers: overall comparison of tri- and tetra-layer systems.

By employing large-scale molecular dynamics simulations of atomistically resolved oligoelectrolytes in aqueous solutions, we study in detail the first four layer-by-layer deposition cycles of an oligoelectrolyte multilayer made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS). The multilayers are grown on a silica substrate in 0.1 M NaCl electrolyte solutions and the swollen structures are then subsequently exposed to varying added salt concentration.

Phase stability of the ice XVII-based CO chiral hydrate from molecular dynamics simulations.

We computed the phase diagram of CO hydrates at high pressure (HP), from 0.3 to 20 kbar, by means of molecular dynamics simulations. The two CO hydrates known to occur in this pressure range are the cubic structure I (sI) clathrate and the HP hydrate, whose water framework is the recently discovered ice XVII.

Effect of general anesthetics on the properties of lipid membranes of various compositions.

Computer simulations of four lipid membranes of different compositions, namely neat DPPC and PSM, and equimolar DPPC-cholesterol and PSM-cholesterol mixtures, are performed in the presence and absence of the general anesthetics diethylether and sevoflurane both at 1 and 600 bar. The results are analyzed in order to identify membrane properties that are potentially related to the molecular mechanism of anesthesia, namely that change in the same way in any membrane with any anesthetics, and change oppositely with increasing pressure. We find that the lateral lipid density satisfies both criteria: it is decreased by anesthetics and increased by pressure.

Electrokinetic droplet transport from electroosmosis to electrophoresis.

Droplet transport in microfluidic channels by electrically induced flows often entails the simultaneous presence of electroosmosis and electrophoresis. Here we make use of coupled lattice-Boltzmann/molecular dynamics simulations to compute the mobility of a droplet in a microchannel under the effect of an external electric field. By varying the droplet solvation free energy of the counterions released at the channel walls, we observe the continuous transition between the electroosmotic and electrophoretic regime.

Pytim: A python package for the interfacial analysis of molecular simulations.

Pytim is a versatile python framework for the analysis of interfacial properties in molecular simulations. The code implements several algorithms for the identification of instantaneous interfaces of arbitrary shape, and analysis tools written specifically for the study of interfacial properties, such as intrinsic profiles. The code is written in the python language, and makes use of the numpy and scipy packages to deliver high computational performances.

The impact of tensorial temperature on equilibrium thermodynamics.

Thermodynamic temperature is a scalar. However, the connection with the kinetic energy tensor in statistical mechanics leaves open the possibility to define a tensorial temperature. This concept has sometimes been used to simulate isothermal conditions in out-of-equilibrium systems.

Phase and interface determination in computer simulations of liquid mixtures with high partial miscibility.

Partially miscible solutions can represent a challenge from the computer simulation standpoint, especially if the mutual solubility of the components is so large that their concentrations do not change much from one phase to another. In this case, identifying which molecules belong to which phase becomes a complicated task. Here, we propose a density-based clustering approach with self-tuning capabilities and apply it to the case of the mixture of an ionic liquid with benzene.

Nonzero Ideal Gas Contribution to the Surface Tension of Water.

Surface tension, the tendency of fluid interfaces to behave elastically and minimize their surface, is routinely calculated as the difference between the lateral and normal components of the pressure or, invoking isotropy in momentum space, of the virial tensor. Here we show that the anisotropy of the kinetic energy tensor close to a liquid-vapor interface can be responsible for a large part of its surface tension (about 15% for water, independent from temperature).

Single Particle Dynamics at the Intrinsic Surface of Various Apolar, Aprotic Dipolar, and Hydrogen Bonding Liquids As Seen from Computer Simulations.

We investigate the single molecule dynamics at the intrinsic liquid/vapor interface of five different molecular liquids (carbon tetrachloride, acetone, acetonitrile, methanol, and water). After assessing that the characteristic residence times in the surface layer are long enough for a meaningful definition of several transport properties within the layer itself, we characterize the dynamics of the individual molecules at the liquid surface by analyzing their normal and lateral mean-square displacements and lateral velocity autocorrelation functions and, in the case of the hydrogen bonding liquids (i.e.

Lateral Pressure Profile and Free Volume Properties in Phospholipid Membranes Containing Anesthetics.

The effect of four general anesthetics, namely chloroform, halothane, diethyl ether, and enflurane on the free volume fraction and lateral pressure profiles in a fully hydrated dipalmitoylphosphatidylcholime (DPPC) membrane is investigated by means of computer simulation. In order to find changes that can be related to the molecular mechanism of anesthesia as well as its pressure reversal, the simulations are performed both at atmospheric and high (1000 bar) pressures. The obtained results show that the additional free volume occurring in the membrane is localized around the anesthetic molecules themselves.