BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Eduardo Sanz

Eduardo Sanz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softphysics.chem-phcond-mat.mtrl-scicond-mat.stat-mechcond-mat.mes-hallphysics.ao-ph

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

14 published item(s)

preprint2026arXiv

Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation

Heterogeneous nucleation is the main path to ice formation on Earth. The ice nucleating ability of a certain substrate is mainly determined by both molecular interactions and the structural mismatch between the ice and the substrate lattices. We focus on the latter factor using molecular simulations of the mW model. Quantifying the effect of structural mismatch alone is challenging due to its coupling with molecular interactions. To disentangle both factors, we use a substrate composed of water molecules in such a way that any variation on the nucleation temperature can be exclusively ascribed to the structural mismatch. We find that a one per cent increase of structural mismatch leads to a decrease of approximately 4 K in the nucleation temperature. We also analyse the effect of the orientation of the substrate with respect to the liquid. The three main ice orientations (basal, primary prism and secondary prism) have a similar ice nucleating ability. We finally asses the effect of lattice flexibility by comparing substrates where molecules are immobile with others where a certain freedom to fluctuate around the lattice positions is allowed. Interestingly, we find that the latter type of substrate is more efficient in nucleating ice because it can adapt its structure to that of ice.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Heterogeneous ice nucleation on model substrates

Ice nucleation is greatly important in areas as diverse as climate change, cryobiology, geology or food industry. Predicting the ability of a substrate to induce the nucleation of ice from supercooled water is a difficult problem. Here, we use molecular simulations to analyse how the ice nucleating ability is affected by the substrate lattice structure and orientation. We focus on different model lattices: simple cubic, body centred cubic and face centred cubic, and assess their ability to induce ice nucleation by calculating nucleation rates. Several orientations are studied for the case of the face centred cubic lattice. Curiously, a hexagonal symmetry does not guarantee a better ice nucleating ability. By comparing the body centred cubic and the cubic lattices we determined that there is a significant role of the underlying crystal plane(s) on ice nucleation. The structure of the liquid layer adjacent to the substrate reveals that more efficient nucleants induce a more structured liquid. The most efficient substrates present a strong sensitivity of their ice nucleating ability to the lattice parameters. Introducing a novel methodological approach, we use Classical Nucleation Theory to estimate the contact angle of the ice nucleus on the studied substrates from the calculated nucleation rates. The method also provides the nucleation free energy barrier height, the kinetic pre-factor and the critical cluster size. The latter is in agreement with the nucleus size obtained through a microscopic analysis of the nucleation trajectories, which supports the validity of Classical Nucleation Theory down to small critical clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Homogeneous Ice Nucleation Rate in Water Droplets

To predict the radiative forcing of clouds it is necessary to know the rate with which ice homogeneously nucleates in supercooled water. Such rate is often measured in drops to avoid the presence of impurities. At large supercooling small (nanoscopic) drops must be used to prevent simultaneous nucleation events. The pressure inside such drops is larger than the atmospheric one by virtue of the Laplace equation. In this work, we take into account such pressure raise in order to predict the nucleation rate in droplets using the TIP4P/Ice water model. We start from a recent estimate of the maximum drop size that can be used at each supercooling avoiding simultaneous nucleation events [Espinosa et al. J. Chem. Phys., 2016]. We then evaluate the pressure inside the drops with the Laplace equation. Finally, we obtain the rate as a function of the supercooling by interpolating our previous results for 1 and 2000 bar [Espinosa et al. Phys. Rev. Lett. 2016] using the Classical Nucleation Theory expression for the rate. This requires, in turn, interpolating the ice-water interfacial free energy and chemical potential difference. The TIP4P/Ice rate curve thus obtained is in good agreement with most droplet-based experiments. In particular, we find a good agreement with measurements performed using nanoscopic drops, that are currently under debate. The successful comparison between model and experiments suggests that TIP4P/Ice is a reliable model to study the water-to-ice transition and that Classical Nucleation Theory is a good framework to understand it.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Seeding Approach to Bubble Nucleation in Superheated Lennard Jones Fluids

We investigate vapor homogeneous nucleation in a superheated Lennard Jones liquid with computer simulations. Special simulation techniques are required to address this study since the nucleation of a critical vapor bubble --one that has equal chance to grow or shrink-- in a moderately superheated liquid is a rare event. We use the Seeding method, that combines Classical Nucleation Theory with computer simulations of a liquid containing a vapor bubble to provide bubble nucleation rates in a wide temperature range. Seeding has been successfully applied to investigate the nucleation of crystals in supercooled fluids and here we apply it for the fist time --to the best of our knowledge-- to the liquid-to-vapor transition. We find that Seeding provides nucleation rates that are consistent with independent calculations not based on the assumptions of Classical Nucleation Theory. Different criteria to determine the radius of the critical bubble give different rate values. The accuracy of each criterion depends of the degree of superheating. Moreover, Seeding simulations show that the surface tension depends on pressure for a given temperature. Therefore, using Classical Nucleation Theory with the coexistence surface tension does not provide good estimates of the nucleation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Seeding Approach to nucleation in the NVT ensemble: the case of bubble cavitation in overstretched Lennard Jones fluids

Simulations are widely used to study nucleation in first order phase transitions due to the fact that they have access to the relevant length and time scales. However, simulations face the problem that nucleation is an activated process. Therefore, rare event simulation techniques are needed to promote the formation of the critical nucleus. The Seeding method, where the simulations are started with the nucleus already formed, has proven quite useful in efficiently providing estimates of the nucleation rate for a wide range of orders of magnitude. So far, Seeding has been employed in the NPT ensemble, where the nucleus either grows or redissolves. Thus, several trajectories have to be run in order to find the thermodynamic conditions that make the seeded nucleus critical. Moreover, the nucleus lifetime is short and the statistics for obtaining its properties is consequently poor. To deal with these shortcomings we extend the Seeding method to the NVT ensemble. We focus on the problem of bubble nucleation in a mestastable Lennard Jones fluid. We show that, in the NVT ensemble, it is possible to equilibrate and stabilise critical bubbles for a long time. The nucleation rate inferred from NVT-Seeding is fully consistent with that coming from NPT-Seeding. The former is quite suitable to obtain the nucleation rate along isotherms, whereas the latter is preferable if the dependence of the rate with temperature at constant pressure is required. Care should be taken with finite size effects when using NVT-Seeding. Further work is required to extend NVT seeding to other sorts of phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure and fluctuations of the premelted liquid film of ice at the triple point

In this paper we study the structure of the ice/vapor interface in the neighborhood of the triple point for the TIP4P/2005 model. We probe the fluctuations of the ice/film and film/vapor surfaces that separate the liquid film from the coexisting bulk phases at basal, primary prismatic and secondary prismatic planes. The results are interpreted using a coupled sine Gordon plus Interface Hamiltonian model. At large length-scales, the two bounding surfaces are correlated and behave as a single complex ice/vapor interface. For small length, on the contrary, the ice/film and film/vapor surfaces behave very much like independent ice/water and water/vapor interfaces. The study suggests that the basal facet of the TIP4P/2005 model is smooth, the prismatic facet is close to a roughening transition, and the secondary prismatic facet is rough. For the faceted basal face, our fluctuation analysis allows us to estimate the step free energy in good agreement with experiment. Our results allow for a quantitative characterization of the extent to which the adsorbed quasi-liquid layer behaves as water, and explains experimental observations which reveal similar activation energies for crystals grown in bulk vapor or bulk water.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Premelting-Induced Smoothening of the Ice-Vapor Interface

We perform computer simulations of the quasiliquid layer of ice formed at the ice-vapor interface close to the ice Ih-liquid-vapor triple point of water. Our study shows that the two distinct surfaces bounding the film behave at small wavelengths as atomically rough and independent ice-water and water-vapor interfaces. For long wavelengths, however, the two surfaces couple, large scale parallel fluctuations are inhibited, and the ice-vapor interface becomes smooth. Our results could help explain the complex morphology of ice crystallites.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A study of the ice-water interface using the TIP4P/2005 water model

In this work we study the ice-water interface under coexistence conditions by means of molecular simulations using the TIP4P/2005 water model. Following the methodology proposed by Hoyt and co-workers [J. J. Hoyt, M. Asta and A. Karma, Phys. Rev. Lett., 86, 5530, (2001)] we measure the interfacial free energy of ice with liquid water by analysing the spectrum of capillary fluctuations of the interface. We get an orientationally averaged interfacial free energy of 27(2) mN/m, in good agreement with a recent estimate obtained from simulation data of the size of critical clusters [E. Sanz, C. Vega, J. R. Espinosa, R. Caballero-Bernal, J. L. F. Abascal and C. Valeriani, JACS, 135, 15008, (2013)]. We also estimate the interfacial free energy of different planes and obtain 27(2), 28(2)and 28(2) mN/m for the basal, the primary prismatic and the secondary prismatic planes respectively. Finally, we inspect the structure of the interface and find that its thickness is of approximately 4-5 molecular diameters. Moreover, we find that when the basal plane is exposed to the fluid the interface alternates regions of cubic ice with regions of hexagonal ice.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Avalanches mediate crystallization in a hard-sphere glass

By molecular-dynamics simulations, we have studied the devitrification (or crystallization) of aged hard-sphere glasses. First, we find that the dynamics of the particles are intermittent: Quiescent periods, when the particles simply "rattle" in their nearest-neighbor cages, are interrupted by abrupt "avalanches," where a subset of particles undergo large rearrangements. Second, we find that crystallization is associated with these avalanches but that the connection is not straightforward. The amount of crystal in the system increases during an avalanche, but most of the particles that become crystalline are different from those involved in the avalanche. Third, the occurrence of the avalanches is a largely stochastic process. Randomizing the velocities of the particles at any time during the simulation leads to a different subsequent series of avalanches. The spatial distribution of avalanching particles appears random, although correlations are found among avalanche initiation events. By contrast, we find that crystallization tends to take place in regions that already show incipient local order.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computer simulation study of surface wave dynamics at the crystal--melt interface

We study, by means of computer simulations, the crystal-melt interface of three different systems: hard-spheres, Lennard Jones and the TIP4P/2005 water model. In particular, we focus on the dynamics of surface waves. We observe that the processes involved in the relaxation of surface waves are characterized by distinct time scales: a slow one related to the continuous recrystallization and melting, that is governed by capillary forces; and a fast one which we suggest to be due to a combination of processes that quickly cause small perturbations to the shape of the interface (like e. g. Rayleigh waves, subdiffusion, or attachment/detachment of particles to/from the crystal). The relaxation of surface waves becomes dominated by the slow process as the wavelength increases. Moreover, we see that the slow relaxation is not influenced by the details of the microscopic dynamics. In a time scale characteristic for the diffusion of the liquid phase, the relaxation dynamics of the crystal-melt interface of water is around one order of magnitude slower than that of Lennard Jones or hard spheres, which we ascribe to the presence of orientational degrees of freedom in the water molecule. Finally, we estimate the rate of crystal growth from our analysis of the capillary wave dynamics and compare it with previous simulation studies and with experiments for the case of water.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From compact to fractal crystalline clusters in concentrated systems of monodisperse hard spheres

We address the crystallization of monodisperse hard spheres in terms of the properties of finite- size crystalline clusters. By means of large scale event-driven Molecular Dynamics simulations, we study systems at different packing fractions ϕ ranging from weakly supersaturated state points to glassy ones, covering different nucleation regimes. We find that such regimes also result in different properties of the crystalline clusters: compact clusters are formed in the classical-nucleation-theory regime (ϕ \leq 0.54), while a crossover to fractal, ramified clusters is encountered upon increasing packing fraction (ϕ \geq 0.56), where nucleation is more spinodal-like. We draw an analogy between macroscopic crystallization of our clusters and percolation of attractive systems to provide ideas on how the packing fraction influences the final structure of the macroscopic crystals. In our previous work (Phys. Rev. Lett., 106, 215701, 2011), we have demonstrated how crystallization from a glass (at ϕ > 0.58) happens via a gradual (many-step) mechanism: in this paper we show how the mechanism of gradual growth seems to hold also in super-saturated systems just above freezing showing that static properties of clusters are not much affected by dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Crystallization Mechanism of Hard Sphere Glasses

In supercooled liquids, vitrification generally suppresses crystallization. Yet some glasses can still crystallize despite the arrest of diffusive motion. This ill-understood process may limit the stability of glasses, but its microscopic mechanism is not yet known. Here we present extensive computer simulations addressing the crystallization of monodisperse hard-sphere glasses at constant volume (as in a colloid experiment). Multiple crystalline patches appear without particles having to diffuse more than one diameter. As these patches grow, the mobility in neighbouring areas is enhanced, creating dynamic heterogeneity with positive feedback. The future crystallization pattern cannot be predicted from the coordinates alone: crystallization proceeds by a sequence of stochastic micro-nucleation events, correlated in space by emergent dynamic heterogeneity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Colloidal gels assembled via a temporary interfacial scaffold

The liquid-liquid phase separation of a binary solvent can be arrested by colloidal particles trapped at the interface [K. Stratford et al., Science 309, 5744 (2005)]. We show experimentally that the colloidal network so formed can remain stable after fully remixing the liquids, creating a new type of gel in which colloids in a single-phase solvent have locally planar coordination. We argue that this structure is likely maintained by primary-minimum DLVO bonding of our charged colloids, created under strong compression by capillary forces. We present simulation evidence that the combination of a short-ranged attraction with a repulsive barrier can strongly stabilize such locally planar gels.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Evidence for out-of-equilibrium crystal nucleation in suspensions of oppositely charged colloids

We report a numerical study of the rate of crystal nucleation in a binary suspension of oppositely charged colloids. Two different crystal structures compete in the thermodynamic conditions under study. We find that the crystal phase that nucleates is metastable and, more surprisingly, its nucleation free energy barrier is not the lowest one. This implies that, during nucleation, there is insufficient time for sub-critical nuclei to relax to their lowest free-energy structure. Such behavior is in direct contradiction with the common assumption that the phase that crystallizes most readily is the one with the lowest free-energy barrier for nucleation. The phenomenon that we describe should be relevant for crystallization experiments where competing solid structures are not connected by an easy transformation.

0 citations0 reviewsNo code linkedOpen access