BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

C. Vega

C. Vega contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.softcond-mat.stat-mechphysics.chem-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Free energy calculations for atomic solids through the Einstein crystal/molecule methodology using GROMACS and LAMMPS

In this work the free energy of solid phases is computed for the Lennard-Jones potential and for a model of NaCl. The free energy is evaluated through the Einstein crystal/molecule methodologies using the Molecular Dynamics programs: GROMACS and LAMMPS. The obtained results are compared with the results obtained from Monte Carlo. Good agreement between the different programs and methodologies was found. The procedure to perform the free energy calculations for the solid phase in the Molecular Dynamic programs is described. Since these programs allow to study any continuous intermolecular potential (when given in a tabulated form) this work shows that for isotropic potentials (describing for instance atomic solids or colloidal particles) free energy calculations can be performed on a routinely basis using GROMACS and/or LAMMPS.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Phase diagram of model anisotropic particles with octahedral symmetry

We computed the phase diagram for a system of model anisotropic particles with six attractive patches in an octahedral arrangement. We chose to study this model for a relatively narrow value of the patch width where the lowest-energy configuration of the system is a simple cubic crystal. At this value of the patch width, there is no stable vapour-liquid phase separation, and there are three other crystalline phases in addition to the simple cubic crystal that is most stable at low pressure. Firstly, at moderate pressures, it is more favourable to form a body-centred cubic crystal, which can be viewed as two interpenetrating, and almost non-interacting, simple cubic lattices.Secondly, at high pressures and low temperatures, an orientationally ordered face-centred cubic structure becomes favourable. Finally, at high temperatures a face-centred cubic plastic crystal is the most stable solid phase.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Equation of state and critical behavior of polymer models: A quantitative comparison between Wertheim's thermodynamic perturbation theory and computer simulations

We present an application of Wertheim's Thermodynamic Perturbation Theory (TPT1) to a simple coarse grained model made of flexibly bonded Lennard-Jones monomers. We use both the Reference Hyper-Netted-Chain (RHNC) and Mean Spherical approximation (MSA) integral equation theories to describe the properties of the reference fluid. The equation of state, the density dependence of the excess chemical potential, and the critical points of the liquid--vapor transition are compared with simulation results and good agreement is found. The RHNC version is somewhat more accurate, while the MSA version has the advantage of being almost analytic. We analyze the scaling behavior of the critical point of chain fluids according to TPT1 and find it to reproduce the mean field exponents: The critical monomer density is predicted to vanish as $n^{-1/2}$ upon increasing the chain length $n$ while the critical temperature is predicted to reach an asymptotic finite temperature that is attained as $n^{-1/2}$. The predicted asymptotic finite critical temperature obtained from the RHNC and MSA versions of TPT1 is found to be in good agreement with the $Θ$ point of our polymer model as obtained from the temperature dependence of the single chain conformations.

0 citations0 reviewsNo code linkedOpen access