Papers by Santi Prestipino
Entropy, 2021
Thermodynamic phases are the most prominent manifestation of emergent behavior [...]
Entropy, 2022
We reconsider model II of Orban et al. (J. Chem. Phys. 1968, 49, 1778–1783), a two-dimensional la... more We reconsider model II of Orban et al. (J. Chem. Phys. 1968, 49, 1778–1783), a two-dimensional lattice-gas system featuring a crystalline phase and two distinct fluid phases (liquid and vapor). In this system, a particle prevents other particles from occupying sites up to third neighbors on the square lattice, while attracting (with decreasing strength) particles sitting at fourth- or fifth-neighbor sites. To make the model more realistic, we assume a finite repulsion at third-neighbor distance, with the result that a second crystalline phase appears at higher pressures. However, the similarity with real-world substances is only partial: Upon closer inspection, the alleged liquid–vapor transition turns out to be a continuous (albeit sharp) crossover, even near the putative triple point. Closer to the standard picture is instead the freezing transition, as we show by computing the free-energy barrier relative to crystal nucleation from the “liquid”.
We focus on the Gibbs free energy Δ G for nucleating a droplet of the stable phase (e.g. solid) i... more We focus on the Gibbs free energy Δ G for nucleating a droplet of the stable phase (e.g. solid) inside the metastable parent phase (e.g. liquid), close to the first-order transition temperature. This quantity is central to the theory of homogeneous nucleation, since it superintends the nucleation rate. We recently introduced a field theory describing the dependence of Δ G on the droplet volume V, taking into account besides the microscopic fuzziness of the droplet-parent interface, also small fluctuations around the spherical shape whose effect, assuming isotropy, was found to be a characteristic logarithmic term. Here we extend this theory, introducing the effect of anisotropy in the surface tension, and show that in the limit of strong anisotropy Δ G(V) once more develops a term logarithmic on V, now with a prefactor of opposite sign with respect to the isotropic case. Based on this result, we argue that the geometrical shape that large solid nuclei mostly prefer could be inferred...
As first shown by H. S. Green in 1952, the entropy of a classical fluid of identical particles ca... more As first shown by H. S. Green in 1952, the entropy of a classical fluid of identical particles can be written as a sum of many-particle contributions, each of them being a distinctive functional of all spatial distribution functions up to a given order. By revisiting the combinatorial derivation of the entropy formula, we argue that a similar correlation expansion holds for the entropy of a crystalline system. We discuss how one- and two-body entropies scale with the size of the crystal, and provide fresh numerical data to check the expectation, grounded on theoretical arguments, that both entropies are extensive quantities.
The properties of a macroscopic assembly of weakly-repulsive bosons at zero temperature are well ... more The properties of a macroscopic assembly of weakly-repulsive bosons at zero temperature are well described by Gross-Pitaevskii mean-field theory. According to this formalism the system exhibits a quantum transition from superfluid to cluster supersolid as a function of pressure. We develop a thermodynamically rigorous treatment of the different phases of the system by adopting a variational formulation of the condensate wave function --- represented as a sum of Gaussians --- that is amenable to exact manipulations. Not only is this description quantitatively accurate, but it is also capable to predict the order (and sometimes even the location) of the transition. We consider a number of crystal structures in two and three dimensions and determine the phase diagram. Depending on the lattice, the transition from fluid to solid can be first-order or continuous, a lower coordination entailing a milder transition. In two dimensions, crystallization would occur at the same pressure on thr...
I employ the van der Waals theory of Baus and coworkers to analyze the fast, adiabatic decay of a... more I employ the van der Waals theory of Baus and coworkers to analyze the fast, adiabatic decay of a supercooled liquid in a closed vessel with which the solidification process usually starts. By imposing a further constraint on either the system volume or pressure, I use the maximum-entropy method to quantify the fraction of liquid that is transformed into solid as a function of undercooling and of the amount of a foreign gas that could possibly be also present in the test tube. Upon looking at the implications of thermal and mechanical insulation for the energy cost of forming a solid droplet within the liquid, I identify one situation where the onset of solidification inevitably occurs near the wall in contact with the bath.
The phase behavior of stabilized dispersions of macromolecules is most easily described in terms ... more The phase behavior of stabilized dispersions of macromolecules is most easily described in terms of the effective interaction between the centers of mass of solute particles. For molecules like polymer chains, dendrimers, etc., the effective pair potential is finite at the origin, allowing "particles" to freely interpenetrate each other. Using a double-Gaussian model (DGM) for demonstration, we studied the behavior of the system as a function of the attraction strength η. Above a critical strength η_ c, the infinite-size system is Ruelle-unstable, in that it collapses to a cluster of finite volume. As η_ c is approached from below, the liquid-vapor region exhibits an anomalous widening at low temperature, and the liquid density apparently diverges at the stability threshold. Above η_ c, the thermodynamic plane is divided in two regions, differing in the value of the average waiting time for collapse, being finite and small on one side of the boundary line, while large or e...
We investigated numerically the high-temperature/high-pressure phase diagram of Xenon as modelled... more We investigated numerically the high-temperature/high-pressure phase diagram of Xenon as modelled through the exp-6 interaction potential, which is thought to provide a reliable description of the thermal behaviour of rare gases under extreme conditions. We performed a series of extensive NVT Monte Carlo simulations which, in conjunction with exact computation of the solid free energy by the Frenkel-Ladd method, allowed us to precisely locate the freezing and the melting thresholds at each temperature. We find that, under isothermal compression, the exp-6 fluid freezes directly into a FCC solid; however, above 4500 K, an intermediate BCC phase becomes stable in a narrow range of pressures. The chemical potential of the HCP phase does never significantly differ from that of the FCC solid of equal T and P, though the former is found to slightly overcome the latter. We discuss our results in the light of previous numerical studies of the same model system and of the experimental data a...
In standard nucleation theory, the nucleation process is characterized by computing ΔΩ(V), the re... more In standard nucleation theory, the nucleation process is characterized by computing ΔΩ(V), the reversible work required to form a cluster of volume V of the stable phase inside the metastable mother phase. However, other quantities besides the volume could play a role in the free energy of cluster formation, and this will in turn affect the nucleation barrier and the shape of the nucleus. Here we exploit our recently introduced mesoscopic theory of nucleation to compute the free energy cost of a nearly-spherical cluster of volume V and a fluctuating surface area A, whereby the maximum of ΔΩ(V) is replaced by a saddle point in ΔΩ(V,A). Compared to the simpler theory based on volume only, the barrier height of ΔΩ(V,A) at the transition state is systematically larger by a few k_BT. More importantly, we show that, depending on the physical situation, the most probable shape of the nucleus may be highly non spherical, even when the surface tension and stiffness of the model are isotropic...
We present a Monte Carlo simulation study of the phase behavior of two-dimensional classical part... more We present a Monte Carlo simulation study of the phase behavior of two-dimensional classical particles repelling each other through an isotropic Gaussian potential. As in the analogous three-dimensional case, a reentrant-melting transition occurs upon compression for not too high temperatures, along with a spectrum of water-like anomalies in the fluid phase. However, in two dimensions melting is a continuous two-stage transition, with an intermediate hexatic phase which becomes increasingly more definite as pressure grows. All available evidence supports the Kosterlitz-Thouless-Halperin-Nelson-Young scenario for this melting transition. We expect that such a phenomenology can be checked in confined monolayers of charge-stabilized colloids with a softened core.
Pair potentials that are bounded at the origin provide an accurate description of the effective i... more 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.
Phase transitions in one-dimensional classical fluids are usually ruled out by making appeal to v... more Phase transitions in one-dimensional classical fluids are usually ruled out by making appeal to van Hove's theorem. A way to circumvent the conclusions of the theorem is to consider an interparticle potential that is everywhere bounded. Such is the case of, e.g., the generalized exponential model of index 4 (GEM-4 potential), which in three dimensions gives a reasonable description of the effective repulsion between flexible dendrimers in a solution. An extensive Monte Carlo simulation of the one-dimensional GEM-4 model [S. Prestipino, Phys. Rev. E 90, 042306 (2014)] has recently provided evidence of an infinite sequence of low-temperature cluster phases, however also suggesting that upon pushing the simulation forward what seemed a true transition may eventually prove to be only a sharp crossover. We hereby investigate this problem theoretically, by three different and increasingly sophisticated approaches ( i.e., a mean-field theory, the transfer matrix of a lattice model of c...
Specialized Monte Carlo methods are nowadays routinely employed, in combination with thermodynami... more 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...
Effective pair interactions with a soft-repulsive component are a well-known feature of polymer s... more Effective pair interactions with a soft-repulsive component are a well-known feature of polymer solutions and colloidal suspensions, but they also provide a key to interpret the high-pressure behaviour of simple elements. We have computed the zero-temperature phase diagram of four different model potentials with various degrees of core softness. Among the reviewed crystal structures, there are also a number of non-Bravais lattices, chosen among those observed in real systems. Some of these crystals are indeed found to be stable for the selected potentials. We recognize an apparently universal trend for unbounded potentials, going from high- to low-coordinated crystal phases and back upon increasing the pressure. Conversely, a bounded repulsion may lead to intermittent appearance of compact structures with compression and no eventual settling down in a specific phase. In both cases, the fluid phase repeatedly reenters at intermediate pressures, as suggested by a cell-theory treatment...
I present a generalization of the Ehrenfest urn model that is aimed at simulating the approach to... more I present a generalization of the Ehrenfest urn model that is aimed at simulating the approach to equilibrium in a dilute gas. The present model differs from the original one in two respects: 1) the two boxes have different volumes and are divided into identical cells with either multiple or single occupancy; 2) particles, which carry also a velocity vector, are subjected to random, but elastic, collisions, both mutual and against the container walls. I show, both analytically and numerically, that the number and energy of particles in a given urn evolve eventually to an equilibrium probability density W which, depending on cell occupancy, is binomial or hypergeometric in the particle number and beta-like in the energy. Moreover, the Boltzmann entropy W takes precisely the same form as the thermodynamic entropy of an ideal gas. This exercise can be useful for pedagogical purposes in that it provides, although in an extremely simplified case, a probabilistic justification for the max...
The approach of an ideal gas to equilibrium is simulated through a generalization of the Ehrenfes... more The approach of an ideal gas to equilibrium is simulated through a generalization of the Ehrenfest ball-and-box model. In the present model, the interior of each box is discretized, i.e., balls/particles live in cells whose occupation can be either multiple or single. Moreover, particles occasionally undergo random, but elastic, collisions between each other and against the container walls. I show, both analitically and numerically, that the number and energy of particles in a given box eventually evolve to an equilibrium distribution W which, depending on cell occupations, is binomial or hypergeometric in the particle number and beta-like in the energy. Furthermore, the long-run probability density of particle velocities is Maxwellian, whereas the Boltzmann entropy W exactly reproduces the ideal-gas entropy. Besides its own interest, this exercise is also relevant for pedagogical purposes since it provides, although in a simple case, an explicit probabilistic foundation for the erg...
An ongoing problem in the study of a classical many-body system is the characterization of its eq... more An ongoing problem in the study of a classical many-body system is the characterization of its equilibrium behaviour by theory or numerical simulation. For purely repulsive particles, locating the melting line in the pressure-temperature plane can be especially hard if the interparticle potential has a softened core or contains some adjustable parameters. A method is hereby presented that yields reliable melting-curve topologies with negligible computational effort. It is obtained by combining the Lindemann melting criterion with a description of the solid phase as an elastic continuum. A number of examples are given in order to illustrate the scope of the method and possible shortcomings. For a two-body repulsion of Gaussian shape, the outcome of the present approach compares favourably with the more accurate but also more computationally demanding self-consistent harmonic approximation.
We reconsider the applicability of classical nucleation theory (CNT) to the calculation of the fr... more We reconsider the applicability of classical nucleation theory (CNT) to the calculation of the free energy of solid cluster formation in a liquid and its use to the evaluation of interface free energies from nucleation barriers. Using two different freezing transitions (hard spheres and NaCl) as test cases, we first observe that the interface-free-energy estimates based on CNT are generally in error. As successive refinements of nucleation-barrier theory, we consider corrections due to a non-sharp solid-liquid interface and to a non-spherical cluster shape. Extensive calculations for the Ising model show that corrections due to a non-sharp and thermally fluctuating interface account for the barrier shape with excellent accuracy. The experimental solid nucleation rates that are measured in colloids are better accounted for by these non-CNT terms, whose effect appears to be crucial in the interpretation of data and in the extraction of the interface tension from them.
Building structures with hierarchical order through the self-assembly of smaller blocks is not on... more Building structures with hierarchical order through the self-assembly of smaller blocks is not only a prerogative of nature, but also a strategy to design artificial materials with tailored functions. We explore in simulation the spontaneous assembly of colloidal particles into extended structures, using spheres and size-asymmetric dimers as solute particles, while treating the solvent implicitly. Besides rigid cores for all particles, we assume an effective short-range attraction between spheres and small monomers to promote, through elementary rules, dimer-mediated aggregation of spheres. Starting from a completely disordered configuration, we follow the evolution of the system at low temperature and density, as a function of the relative concentration of the two species. When spheres and large monomers are of same size, we observe the onset of elongated aggregates of spheres, either disconnected or cross-linked, and a crystalline bilayer. As spheres grow bigger, the self-assembli...
Two-dimensional crystals of classical particles are very peculiar in that melting may occur in tw... more Two-dimensional crystals of classical particles are very peculiar in that melting may occur in two steps, in a continuous fashion, via an intermediate hexatic fluid phase exhibiting quasi-long-range orientational order. On the other hand, three-dimensional spheres repelling each other through a fast-decaying bounded potential of generalized-exponential shape (GEM4 potential) can undergo freezing into cluster crystals, allowing for more that one particle per lattice site. We hereby study the combined effect of low spatial dimensionality and extreme potential softness, by investigating the phase behavior of the two-dimensional (2D) GEM4 system. Using a combination of density-functional theory and numerical free-energy calculations, we show that the 2D GEM4 system displays one ordinary and several cluster triangular-crystal phases, and that only the ordinary crystal first melts into a hexatic phase. Upon heating, the difference between the various cluster crystals fades away, eventuall...
Uploads
Papers by Santi Prestipino