Prediction of aqueous solubility of a strongly soluble solute from molecular simulation.

The prediction of solubilities of compounds by means of molecular simulation has been receiving increasing attention due to the key role played by solubility in countless applications. We have predicted the aqueous solubility of urea at 300 K from chemical potential calculations for two urea model combinations: Özpinar/TIP3P and Hölzl/(TIP4P/2005). The methodology assumes that the intramolecular contribution of the urea molecule to the chemical potentials is identical in the crystal and in solution and, hence, cancels out. In parallel to the chemical potential calculations, we also performed direct coexistence simulations of a urea crystal slab in contact with urea-water solutions with the aim to identify upper and lower bounds to the solubility value using an independent route. The chemical potential approach yielded similar solubilities for both urea models, despite the actual chemical potential values showing a significant dependence on the force field. The predicted solubilities for the two models were 0.013-0.018 (Özpinar) and 0.008-0.012 (Hölzl) mole fraction, which are an order of magnitude lower than the experimental solubility that lies in a range of 0.125-0.216 mole fraction. The direct coexistence solubility bounds were relatively wide and did not encompass the chemical potential based solubilities, although the latter were close to the lower bound values.

Interaction of Water and Oxygen Molecules with Phosphorene: An Ab Initio Study.

Phosphorene, the 2D form of black phosphorus, has recently attracted interest for optoelectronic and tribological applications. However, its promising properties are affected by the strong tendency of the layers to oxidize in ambient conditions. A significant effort has been made to identify the role of oxygen and water in the oxidation process. In this work, we introduce a first-principles study of the phosphorene phase diagram and provide a quantitative estimate of the interaction of pristine and fully oxidized phosphorene layers with oxygen and water molecules. Specifically, we study oxidized layers with oxygen coverages of 25% and 50% that keep the typical anisotropic structure of the layers. We found that hydroxilated and hydrogenated phosphorene layers are both energetically unfavorable, leading to structural distortions. We also studied the water physisorption on both pristine and oxidized layers, finding that the adsorption energy gain doubled on the oxidized layers, whereas dissociative chemisorption was always energetically unfavorable. At the same time, further oxidation (i.e., the dissociative chemisorption of O2) was always favorable, even on oxidized layers. Ab initio molecular dynamics simulations of water intercalated between sliding phosphorene layers showed that even under harsh tribological conditions water dissociation was not activated, thus further strengthening the results obtained from our static calculations. Overall, our results provide a quantitative description of the interaction of phosphorene with chemical species that are commonly found in ambient conditions at different concentrations. The phase diagram that we introduced confirms the tendency of phosphorene layers to fully oxidize due to the presence of O2, resulting in a material with improved hydrophilicity, a piece of information that is relevant for the application of phosphorene, e.g., as a solid lubricant. At the same time, the structural deformations found for the H- and OH- terminated layers undermine their electrical, mechanical, and tribological anisotropic properties and, therefore, the usage of phosphorene.

Thermal Diffusion in Binary Mixtures: Transient Behavior and Transport Coefficients from Equilibrium and Nonequilibrium Molecular Dynamics.

Equilibrium and nonequilibrium molecular dynamics simulations are combined to compute the full set of coefficients that appear in the phenomenological equations describing thermal transport in a binary mixture subject to a constant thermal gradient. The Dynamical Non-Equilibrium Molecular Dynamics approach (D-NEMD) is employed to obtain the microscopic time evolution of the density and temperature fields, together with that of the mass and energy fluxes. D-NEMD enables one to study not only the steady state, but also the evolution of the fields during the transient that follows the onset of the thermal gradient, up to the establishment of the steady state. This makes it possible to ensure that the system has indeed reached a stationary condition, and to analyze the transient mechanisms and time scales of the mass and energy transport. A local time averaging procedure is applied to each trajectory contributing to the calculation to improve the signal-to-noise ratio in the estimation of the fluxes and to obtain a clear signal with the, relatively limited, statistics available.

Size-dependent commensurability and its possible role in determining the frictional behavior of adsorbed systems.

Recent nanofriction experiments of xenon on graphene revealed that the slip onset can be induced by increasing the adsorbate coverage above a critical value, which depends on temperature. Moreover, the xenon slippage on gold is much higher than on graphene in spite of the same physical nature of the interactions. To shed light on these intriguing results we have performed molecular dynamics simulations relying on ab initio derived potentials. By monitoring the interfacial structure factor as a function of coverage and temperature, we show that the key mechanism to interpret the observed frictional phenomena is the size-dependence of the island commensurability. The latter quantity is deeply affected also by the lattice misfit, which explains the different frictional behavior of Xe on graphene and gold.

Transient behavior of a model fluid under applied shear.

We study the transient behavior of a model fluid composed by soft repulsive spheres subjected to a planar uniform shear. To this aim, we use a dynamical non-equilibrium molecular dynamics method originally developed by Ciccotti and Jacucci [Phys. Rev. Lett. 35, 789 (1975)] and recently applied to the study of the transient regimes in various fluid systems. We show that the dynamical method allows one to study the transient behavior of the viscous time-dependent response over a wide range of applied shear rates, provided that a temperature control is enforced on the system. In this study, we adopt in particular the configurational thermostat of Braga and Travis [J. Chem. Phys. 123, 134101 (2005)]. The initial behavior of the dynamical response to a θ-like perturbation is characterized by a rapid increase, culminating in a pronounced peak, later relaxing to a plateau value. The latter positively reproduces the values of the viscosity observed in standard steady-state non-equilibrium molecular dynamics.

Onset of frictional slip by domain nucleation in adsorbed monolayers.

It has been known for centuries that a body in contact with a substrate will start to slide when the lateral force exceeds the static friction force. Yet the microscopic mechanisms ruling the crossover from static to dynamic friction are still the object of active research. Here, we analyze the onset of slip of a xenon (Xe) monolayer sliding on a copper (Cu) substrate. We consider thermal-activated creep under a small external lateral force, and observe that slip proceeds by the nucleation and growth of domains in the commensurate interface between the film and the substrate. We measure the activation energy for the nucleation process considering its dependence on the external force, the substrate corrugation, and particle interactions in the film. To understand the results, we use the classical theory of nucleation and compute analytically the activation energy which turns out to be in excellent agreement with numerical results. We discuss the relevance of our results to understand experiments on the sliding of adsorbed monolayers.

High-throughput screening of the static friction and ideal cleavage strength of solid interfaces.

We present a comprehensive ab initio, high-throughput study of the frictional and cleavage strengths of interfaces of elemental crystals with different orientations. It is based on the detailed analysis of the adhesion energy as a function of lateral, γ(x, y), and perpendicular displacements, γ(z), with respect to the considered interface plane. We use the large amount of computed data to derive fundamental insight into the relation of the ideal strength of an interface plane with its adhesion. Moreover, the ratio between the frictional and cleavage strengths is provided as good indicator for the material failure mode - dislocation propagation versus crack nucleation. All raw and curated data are made available to be used as input parameters for continuum mechanic models, benchmarks, or further analysis.

Probabilistic Derivation of Spatiotemporal Correlation Functions in the Hydrodynamic Limit.

In this paper, we use probability theory to prove in suitable conditions the equivalence of equilibrium time correlation functions of microscopic density fields with the time correlation functions of local macroscopic density fields evolved by hydrodynamics in (approximate) phenomenological continuum theories of matter. We further discuss a useful and rigorous numerical algorithm, derived from this framework, to compute macroscopic space- and time-dependent behaviors (such as the hydrodynamical one) via molecular dynamics simulations.