Researcher profile

E. Sanz

E. Sanz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

cond-mat.soft cond-mat.stat-mech cond-mat.mtrl-sci physics.chem-ph physics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging

7works

0followers

5topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Equivalence between condensation and boiling in a Lennard Jones fluid

Condensation and boiling are phase transitions highly relevant to industry, geology or atmospheric science. These phase transitions are initiated by the nucleation of a drop in a supersaturated vapor and of a bubble in an overstretched liquid respectively. The surface tension between both phases, liquid and vapor, is a key parameter in the development of such nucleation stage. Whereas the surface tension can be readily measured for a flat interface, there are technical and conceptual limitations to obtain it for the curved interface of the nucleus. On the technical side, it is quite difficult to observe a critical nucleus in experiments. From a conceptual point of view, the interfacial free energy depends on the choice of the dividing surface, being the surface of tension the one relevant for nucleation. We bypass the technical limitation by performing simulations of a Lennard Jones fluid where we equilibrate critical nuclei (both drops and bubbles). Regarding the conceptual hurdle, we find the relevant cluster size by searching the radius that correctly predicts nucleation rates and nucleation free energy barriers when combined with Classical Nucleation Theory. With such definition of the cluster size we find the same value of the surface tension for drops and bubbles of a given radius. Thus, condensation and boiling can be viewed as two sides of the same coin. Finally, we combine the data coming from drops and bubbles to obtain, via two different routes, estimates of the Tolman length, a parameter that allows describing the curvature dependence of the surface tension in a theoretical framework.

preprint2014arXiv

Homogeneous ice nucleation evaluated for several water models

In this work, we evaluate by means of computer simulations the rate for ice homogeneous nucleation for several water models such as TIP4P, TIP4P/2005,TIP4P/ICE, and mW (following the same procedure as in Sanz et al. [J. Am. Chem. Soc.135, 15008 (2013)]) in a broad temperature range. We estimate the ice-liquid interfacial free-energy, and conclude that for all water models γ decreases as the temperature decreases. Extrapolating our results to the melting temperature, we obtain a value of the interfacial free-energy between 25 and 32 mN/m in reasonable agreement with the reported experimental values. Moreover, we observe that the values of γ depend on the chosen water model and this is a key factor when numerically evaluating nucleation rates, given that the kinetic prefactor is quite similar for all water models with the exception of the mW (due to the absence of hydrogens). Somewhat surprisingly the estimates of the nucleation rates found in this work for TIP4P/2005 are slightly higher than those of the mW model, even though the former has explicit hydrogens. Our results suggest that it may be possible to observe in computer simulations spontaneous crystallization of TIP4P/2005 at about 60 K below the melting point.

preprint2014arXiv

The mold integration method for the calculation of the crystal-fluid interfacial free energy from simulations

The interfacial free energy between a crystal and a fluid, γ cf, is a highly relevant parameter in phenomena such as wetting or crystal nucleation and growth. Due to the difficulty of measuring γ cf experimentally, computer simulations are often used to study the crystal-fluid interface. Here, we present a novel simulation methodology for the calculation of γ cf . The methodology consists in using a mold composed of potential energy wells to induce the formation of a crystal slab in the fluid at coexistence conditions. This induction is done along a reversible pathway along which the free energy difference between the initial and the final states is obtained by means of thermodynamic integration. The structure of the mold is given by that of the crystal lattice planes, which allows to easily obtain the free energy for different crystal orientations. The method is validated by calculating γ cf for previously studied systems, namely, the hard spheres and the Lennard-Jones systems. Our results for the latter show that the method is accurate enough to deal with the anisotropy of γ cf with respect to the crystal orientation. We also calculate γ cf for a recently proposed continuous version of the hard sphere potential and obtain the same γ cf as for the pure hard sphere system. The method can be implemented both in Monte Carlo and Molecular Dynamics. In fact, we show that it can be easily used in combination with the popular Molecular Dynamics package GROMACS.

preprint2013arXiv

Homogeneous ice nucleation at moderate supercooling from molecular simulation

Among all the freezing transitions, that of water into ice is probably the most relevant to biology, physics, geology or atmospheric science. In this work we investigate homogeneous ice nucleation by means of computer simulations. We evaluate the size of the critical cluster and the nucleation rate for temperatures ranging between 15K and 35K below melting. We use the TIP4P/2005 and the TIP4P/Ice water models. Both give similar results when compared at the same temperature difference with the model's melting temperature. The size of the critical cluster varies from $\sim$8000 molecules (radius$ = 4$nm) at 15K below melting to $\sim$600 molecules (radius$ = 1.7$nm) at 35K below melting. We use Classical Nucleation Theory (CNT) to estimate the ice-water interfacial free energy and the nucleation free energy barrier. We obtain an interfacial free energy of 29(3)mN/m from an extrapolation of our results to the melting temperature. This value is in good agreement both with experimental measurements and with previous estimates from computer simulations of TIP4P-like models. Moreover, we obtain estimates of the nucleation rate from simulations of the critical cluster at the barrier top. The values we get for both models agree within statistical error with experimental measurements. At temperatures higher than 20K below melting we get nucleation rates slower than the appearance of a critical cluster in all the water of the hydrosphere in the age of the universe. Therefore, our simulations predict that water freezing above this temperature must necessarily be heterogeneous.

preprint2012arXiv

Calculation of the melting point of alkali halides by means of computer simulations

In this manuscript we study the liquid-solid coexistence of NaCl-type alkali halides, described by interaction potentials such as Tosi-Fumi (TF), Smith-Dang (SD) and Joung-Cheatham (JC), and compute their melting temperature (Tm) at 1 bar via three independent routes: 1) liquid/solid direct coexistence, 2) free-energy calculations and 3) Hamiltonian Gibbs-Duhem integration. The melting points obtained by the three routes are consistent with each other. The calculated Tm of the Tosi-Fumi model of NaCl is in good agreement with the experimental value as well as with other numerical calculations. However, the other two models considered for NaCl, SD and JC, overestimate the melting temperature of NaCl by more than 200 K. We have also computed the melting temperature of other alkali halides using the Tosi-Fumi interaction potential and observed that the predictions are not always as close to the experimental values as they are for NaCl. It seems that there is still room for improvement in the area of force-fields for alkaline halides, given that so far most models are still unable to describe a simple yet important property such as the melting point.

preprint2011arXiv

Irreducible finite-size effects in surface free energies from crystal-nucleation data

In this Letter we report a simulation study in which we compare the solid-liquid interfacial free energy of NaCl at coexistence, with the value that follows from the height of the homogeneous nucleation barrier. We find that the two estimates differ by more than 100%. Similar, although smaller discrepancies are found for crystals of hard-sphere colloids and of Lennard-Jones ("argon") particles. We consider a variety of possible causes for this discrepancy and are forced to conclude that it is due to a finite-size effect that cannot be corrected for by any simple thermodynamic procedure. Importantly, we find that the surface free energies that follow from real nucleation experiments should be subject to large finite size effects. Taking this in to account, we obtain quantitative agreement between the simulation data and the surface free energy of NaCl that follows from nucleation experiments. Our finding suggests that most published solid-liquid surface free energies derived from nucleation experiments will have to be revised.

preprint2009arXiv

Hard Spheres: Crystallization and Glass Formation

Motivated by old experiments on colloidal suspensions, we report molecular dynamics simulations of assemblies of hard spheres, addressing crystallization and glass formation. The simulations cover wide ranges of polydispersity s (standard deviation of the particle size distribution divided by its mean) and particle concentration. No crystallization is observed for s > 0.07. For 0.02 < s < 0.07, we find that increasing the polydispersity at a given concentration slows down crystal nucleation. The main effect here is that polydispersity reduces the supersaturation since it tends to stabilise the fluid but to destabilise the crystal. At a given polydispersity (< 0.07) we find three regimes of nucleation: standard nucleation and growth at concentrations in and slightly above the coexistence region; "spinodal nucleation", where the free energy barrier to nucleation appears to be negligible, at intermediate concentrations; and, at the highest concentrations, a new mechanism, still to be fully understood, which only requires small re-arrangement of the particle positions. The cross-over between the second and third regimes occurs at a concentration, around 58% by volume, where the colloid experiments show a marked change in the nature of the crystals formed and the particle dynamics indicate an "ideal" glass transition.