Source author record

G. Malescio

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

4 published item(s)

preprint2015arXiv

On the determination of phase boundaries via thermodynamic integration across coexistence regions

Specialized Monte Carlo methods are nowadays routinely employed, in combination with thermodynamic integration (TI), to locate phase boundaries of classical many-particle systems. This is especially useful for the fluid-solid transition, where a critical point does not exist and both phases may notoriously go deeply metastable. Using the Lennard-Jones model for demonstration, we hereby investigate on the alternate possibility of tracing reasonably accurate transition lines directly by integrating the pressure equation of state computed in a canonical-ensemble simulation with local moves. The recourse to this method would become a necessity when the stable crystal structure is not known. We show that, rather counterintuitively, metastability problems can be alleviated by reducing (rather than increasing) the size of the system. In particular, the location of liquid-vapor coexistence can exactly be predicted by just TI. On the contrary, TI badly fails in the solid-liquid region, where a better assessment (to within 10\% accuracy) of the coexistence pressure can be made by following the expansion, until melting, of the defective solid which has previously emerged from the decay of the metastable liquid.

preprint2015arXiv

Shapes of a liquid droplet in a periodic box

Within the coexistence region between liquid and vapor the equilibrium pressure of a simulated fluid exhibits characteristic jumps and plateaus when plotted as a function of density at constant temperature. These features exclusively pertain to a finite-size sample in a periodic box, as they are washed out in the bulk limit. Below the critical density, at each pressure jump the shape of the liquid drop undergoes a morphological transition, changing from spherical to cylindrical to slab-like as the density is increased. We formulate a simple theory of these shape transitions, which is adapted from a calculation originally developed by Binder and coworkers [{\em J. Chem. Phys.} {\bf 120}, 5293 (2004)]. Our focus is on the pressure equation of state (rather than on the chemical potential, as in the original work) and includes an extension to elongated boxes. Predictions based on this theory well agree with extensive Monte Carlo data for the cut-and-shifted Lennard-Jones fluid. We further discuss on the thermodynamic stability of liquid drops with shapes other than the three mentioned above, like those found deep inside the liquid-vapor region in simulations starting from scratch. Our theory classifies these more elaborate shapes as metastable.

preprint2014arXiv

Phase behavior of a double-Gaussian fluid displaying water-like features

Pair potentials that are bounded at the origin provide an accurate description of the effective interaction for many systems of dissolved soft macromolecules (e.g., flexible dendrimers). Using numerical free-energy calculations, we reconstruct the equilibrium phase diagram of a system of particles interacting through a potential that brings together a Gaussian repulsion with a much weaker Gaussian attraction, close to the thermodynamic stability threshold. Compared to the purely-repulsive model, only the reentrant branch of the melting line survives, since for lower densities solidification is overridden by liquid-vapor separation. As a result, the phase diagram of the system recalls that of water up to moderate (i.e., a few tens MPa) pressures. Upon superimposing a suitable hard core on the double-Gaussian potential, a further transition to a more compact solid phase is induced at high pressure, which might be regarded as the analog of the ice I to ice III transition in water.

preprint2014arXiv

Twofold reentrant melting in a double-Gaussian fluid

Isotropic pair potentials that are bounded at the origin have been proposed from time to time as models of the effective interaction between macromolecules of interest in the chemical physics of soft matter. We present a thorough study of the phase behavior of point particles interacting through a potential which combines a bounded short-range repulsion with a much weaker attraction at moderate distances, both of Gaussian shape. Notwithstanding the fact that the attraction acts as a small perturbation of the Gaussian-core model potential, the phase diagram of the double-Gaussian model (DGM) is far richer, showing two fluid phases and four distinct solid phases in the case that we have studied. Using free-energy calculations, the various regions of confluence of three distinct phases in the DGM system have all been characterized in detail. Moreover, two distinct lines of reentrant melting are found, and for each of them a rationale is provided in terms of the elastic properties of the solid phases.