Source author record

Pablo G. Debenedetti

Pablo G. Debenedetti 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

18works
5topics
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

18 published item(s)

preprint2022arXiv

Formation of Dissipative Structures in Microscopic Models of Mixtures with Species Interconversion

The separation of substances into different phases is ubiquitous in nature and important scientifically and technologically. This phenomenon may become drastically different if the species involved, whether molecules or supramolecular assemblies, interconvert. In the presence of an external force large enough to overcome energetic differences between the interconvertible species (forced interconversion), the two alternative species will be present in equal amounts, and the striking phenomenon of steady-state, restricted phase separation into mesoscales is observed. Such microphase separation is one of the simplest examples of dissipative structures in condensed matter. In this work, we investigate the formation of such mesoscale steady-state structures through Monte Carlo and Molecular Dynamics simulations of three physically distinct microscopic models of binary mixtures that exhibit both equilibrium (natural) interconversion and a nonequilibrium source of forced interconversion. We show that this source can be introduced through an internal imbalance of intermolecular forces or an external flux of energy that promotes molecular interconversion, possible manifestations of which could include the internal nonequilibrium environment of living cells or a flux of photons. The main trends and observations from the simulations are well captured by a non-equilibrium thermodynamic theory of phase transitions affected by interconversion. We show how a nonequilibrium bicontinuous microemulsion or a spatially modulated state may be generated depending on the interplay between diffusion, natural interconversion, and forced interconversion.

preprint2022arXiv

Homogeneous ice nucleation in an ab initio machine learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine learning model trained on density-functional theory (DFT) energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation) which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

preprint2016arXiv

An Equation of State For The TIP4P/2005 Model of Water Including Negative Pressures

One of the most promising frameworks for understanding the anomalies of cold and especially supercooled water is that of two-structure thermodynamics, in which water is viewed as a non-ideal mixture of two interconvertible local structures. The non-ideality of this mixture may give rise, at very low temperatures, to a liquid-liquid phase transition (LLPT) and a liquid-liquid critical point (LLCP), at which thermodynamic response functions diverge. Various versions of the "two-structure equation of state" (TSEOS) based on this concept have shown remarkable agreement with both experimental data in real water and simulation results. However, recent experiments probing supercooled water at negative pressures reveal the inadequacy of extrapolations of equations of state developed for positive pressures, and have begun to shed additional light on the source of the anomalies of supercooled water. We have analyzed simulation results for the TIP4P/2005 model over a broad range of positive and negative pressures from ambient temperature to deep supercooling. We find that by explicitly incorporating a liquid-vapor spinodal into a two-structure equation of state, we are able to match the simulation data in TIP4P/2005 with striking accuracy. In particular, our equation of state reproduces the observed lines of minima and maxima in the density, isothermal compressibility, and isobaric heat capacity. Contrary to scenarios in which a retracing spinodal accounts for the thermodynamic anomalies of water, we find that the liquid-vapor spinodal in this model continues monotonically to lower pressures as temperature is decreased, influencing but not giving rise to the locus of density maxima and other thermodynamic anomalies. We explain the behavior of TIP4P/2005 in terms of two phenomena: the competition between two local structures and a monotonic liquid-vapor spinodal.

preprint2016arXiv

Concentration Fluctuations and Capacitive Response in Dense Ionic Solutions

We use molecular dynamics simulations in a constant potential ensemble to study the effects of solution composition on the electrochemical response of a double layer capacitor. We find that the capacitance first increases with ion concentration following its expected ideal solution behavior, but decreases upon approaching a pure ionic liquid in agreement with recent experimental observations. The non-monotonic behavior of the capacitance as a function of ion concentration results from the competition between the independent motion of solvated ions in the dilute regime and solvation fluctuations in the concentrated regime. Mirroring the capacitance, we find that the characteristic decay length of charge density correlations away from the electrode is also non-monotonic. The correlation length first decreases with ion concentration as a result of better electrostatic screening but increases with ion concentration as a result of enhanced steric interactions. When charge fluctuations induced by correlated ion-solvent fluctuations are large relative to those induced by the pure ionic liquid, such capacitive behavior is expected to be generic.

preprint2016arXiv

Density and bond-orientational relaxations in supercooled water

Recent computational studies have reported evidence of a metastable liquid-liquid phase transition (LLPT) in molecular models of water under deeply supercooled conditions. A competing hypothesis suggests, however, that non-equilibrium artifacts associated with coarsening of the stable crystal phase have been mistaken for an LLPT in these models. Such artifacts are posited to arise due to a separation of time scales in which density fluctuations in the supercooled liquid relax orders of magnitude faster than those associated with bond-orientational order. Here, we use molecular simulation to investigate the relaxation of density and bond-orientational fluctuations in three molecular models of water (ST2, TIP5P and TIP4P/2005) in the vicinity of their reported LLPT. For each model, we find that density is the slowly relaxing variable under such conditions. We also observe similar behavior in the coarse-grained mW model of water. Our findings therefore challenge the key physical assumption underlying the competing hypothesis.e find that density relaxes significantly faster than bond-orientational order, as incorrectly predicted by this competing hypothesis.

preprint2016arXiv

Two-State Thermodynamics and the Possibility of a Liquid-Liquid Phase Transition in Supercooled TIP4P/2005 Water

Water shows intriguing thermodynamic and dynamic anomalies in the supercooled liquid state. One possible explanation of the origin of these anomalies lies in the existence of a metastable liquid-liquid phase transition (LLPT) between two (high and low density) forms of water. While the anomalies are observed in experiments on bulk and confined water and by computer simulation studies of water-like models, the existence of a LLPT in water is still debated. Unambiguous experimental proof of the existence of a LLPT in bulk supercooled water is hampered by fast ice nucleation which is a precursor of the hypothesized LLPT. Moreover, the hypothesized LLPT, being metastable, in principle cannot exist in the thermodynamic limit (infinite size, infinite time). Therefore, computer simulations of water models are crucial for exploring the possibility of the metastable LLPT and the nature of the anomalies. In this work, we present new simulation results in the NVT ensemble for one of the most accurate classical molecular models of water, TIP4P/2005. To describe the computed properties and explore the possibility of a LLPT we have applied two-structure thermodynamics, viewing water as a non-ideal mixture of two interconvertible local structures ("states"). The results suggest the presence of a liquid-liquid critical point and a LLPT in this model for the simulated length and time scales. We have compared the behavior of TIP4P/2005 with other popular water-like models, namely mW and ST2, and with real water, all of which are well described by two-state thermodynamics. In view of the current debate involving different studies of TIP4P/2005, we discuss consequences of metastability and finite size in observing the liquid-liquid separation. We also address the relationship between the phenomenological order parameter of two-structure thermodynamics and the microscopic nature of the low-density structure.

preprint2015arXiv

Computational investigation of structure, dynamics and nucleation kinetics of a family of modified Stillinger-Weber model fluids in bulk and free-standing thin films

In recent decades, computer simulations have found increasingly widespread use as powerful tools of studying phase transitions in wide variety of systems. In the particular and very important case of aqueous systems, the commonly used force-fields tend to offer quite different predictions with respect to a wide range of thermodynamic and kinetic properties, including the ease of ice nucleation, the propensity to freeze at a vapor-liquid interface, and the existence of a liquid-liquid phase transition. It is thus of fundamental and practical interest to understand how different features of a given water model affect its thermodynamic and kinetic properties. In this work, we use the forward-flux sampling technique to study the crystallization kinetics of a family of modified Stillinger-Weber (SW) potentials with energy ($ε$) and length ($σ$) scales taken from the monoatomic water (mW) model, but with different tetrahedrality parameters ($λ$). By increasing $λ$ from 21 to 24, we observe the nucleation rate increases by 48 orders of magnitude at a supercooling of $ζ=T/T_m=0.845$. Using classical nucleation theory, we are able to demonstrate that this change can largely be accounted for by the increase in $|Δμ|$, the thermodynamic driving force. We also perform rate calculations in freestanding thin films of the supercooled liquid, and observe a crossover from a surface-enhanced crystallization at $λ= 21$ to a bulk-dominated crystallization for $λ\ge22$.

preprint2015arXiv

Direct Calculation of Ice Homogeneous Nucleation Rate for a Molecular Model of Water

Ice formation is ubiquitous in nature, with important consequences in a variety of environments, including biological cells, soil, aircraft, transportation infrastructure and atmospheric clouds. However, its intrinsic kinetics and microscopic mechanism are difficult to discern with current experiments. Molecular simulations of ice nucleation are also challenging, and direct rate calculations have only been performed for coarse-grained models of wate. For molecular models, only indirect estimates have been obtained, e.g. by assuming the validity of classical nucleation theory. We use a path sampling approach to perform the first direct rate calculation of homogeneous nucleation of ice in a molecular model of water. We use TIP4P/Ice, the most accurate among existing molecular models for studying ice polymorphs. By using a novel topological approach to distinguish different polymorphs, we are able to identify a freezing mechanism that involves a competition between cubic and hexagonal ice in the early stages of nucleation. In this competition, the cubic polymorph takes over since the addition of new topological structural motifs consistent with cubic ice leads to the formation of more compact crystallites. This is not true for topological hexagonal motifs, which give rise to elongated crystallites that are not able to grow. This leads to transition states that are rich in cubic ice, and not the thermodynamically stable hexagonal polymorph. This mechanism provides a molecular explanation to the earlier experimental and computational observations of the preference for cubic ice in the literature.

preprint2015arXiv

Pathways to dewetting in hydrophobic confinement

Liquid water can become metastable with respect to its vapor in hydrophobic confinement. The resulting dewetting transitions are often impeded by large kinetic barriers. According to macroscopic theory, such barriers arise from the free energy required to nucleate a critical vapor tube that spans the region between two hydrophobic surfaces - tubes with smaller radii collapse, whereas larger ones grow to dry the entire confined region. Using extensive molecular simulations of water between two nanoscopic hydrophobic surfaces, in conjunction with advanced sampling techniques, here we show that for inter-surface separations that thermodynamically favor dewetting, the barrier to dewetting does not correspond to the formation of a (classical) critical vapor tube. Instead, it corresponds to an abrupt transition from an isolated cavity adjacent to one of the confining surfaces to a gap-spanning vapor tube that is already larger than the critical vapor tube anticipated by macroscopic theory. Correspondingly, the barrier to dewetting is also smaller than the classical expectation. We show that the peculiar nature of water density fluctuations adjacent to extended hydrophobic surfaces - namely, the enhanced likelihood of observing low-density fluctuations relative to Gaussian statistics - facilitates this non-classical behavior. By stabilizing isolated cavities relative to vapor tubes, enhanced water density fluctuations thus stabilize novel pathways, which circumvent the classical barriers and offer diminished resistance to dewetting. Our results thus suggest a key role for fluctuations in speeding up the kinetics of numerous phenomena ranging from Cassie-Wenzel transitions on superhydrophobic surfaces, to hydrophobically-driven biomolecular folding and assembly.

preprint2015arXiv

Thermodynamic and Kinetic Anisotropies in Octane Thin Films

Confinement breaks the translational symmetry of materials. Such symmetry breaking can be used to obtain configurations that are not otherwise accessible in the bulk. Here, we explore the effect of substrate-liquid interactions on the induced thermodynamic and kinetic anisotropies. We consider n-octane nanofilms that are in contact with substrates with varying degrees of attraction. Complete freezing of octane nanofilms is observed at low temperatures, while at intermediate temperatures, a frozen monolayer emerges at both interfaces. By carefully inspecting the profiles of translational and orientational relaxation times, we confirm that the translational and orientational degrees of freedom are decoupled at these frozen monolayers. At sufficiently high temperatures, however, free interfaces and solid-liquid interfaces close to loose substrates undergo pre-freezing, characterized by mild peaks in several thermodynamic quantities. Two distinct dynamic regimes are observed. The dynamics is accelerated in the vicinity of loose substrates, while sticky substrates decelerate dynamics, sometimes by as much as two orders of magnitude. These two distinct dynamical regimes have been previously by us [JCP 141: 024506, 2014] for a model atomic glass-forming liquid. We also confirm the existence of two correlations proposed in the above-mentioned work in solid-liquid subsurface regions of octane films, i.e., a correlation between density and normal stress, and between atomic translational relaxation time and lateral stress. Finally, we inspect the ability of different regions of a film to explore the potential energy landscape, and observe no noticeable difference between the free surface and the bulk. This is unlike the films of model atomic glass formers that tend to sample their respective landscape more efficiently at free surfaces.

preprint2014arXiv

Response to Comment [arXiv:1407.6854] on Palmer et al., Nature, 510, 385, 2014

We respond to a Comment [arXiv:1407.6854 (2014)] on our recent Nature paper [Nature, 510, 385 (2014)]. We categorically disprove the arguments provided in arXiv:1407.6854 (2014) and thereby further substantiate the evidence we presented in our recent study, demonstrating the existence of a metastable liquid-liquid transition in a molecular model of water. We will make our code publicly available shortly along with proper user documentation that is currently under development.

preprint2014arXiv

Suppression of Sub-surface Freezing in Free-Standing Thin Films of a Coarse-grained Model of Water

Freezing in the vicinity of water-vapor interfaces is of considerable interest to a wide range of disciplines, most notably the atmospheric sciences. In this work, we use molecular dynamics and two advanced sampling techniques, forward flux sampling and umbrella sampling, to study homogeneous nucleation of ice in free-standing thin films of supercooled water. We use a coarse-grained monoatomic model of water, known as mW, and we find that in this model a vapor-liquid interface suppresses crystallization in its vicinity. This suppression occurs in the vicinity of flat interfaces where no net Laplace pressure in induced. Our free energy calculations reveal that the pre-critical crystalline nuclei that emerge near the interface are thermodynamically less stable than those that emerge in the bulk. We investigate the origin of this instability by computing the average asphericity of nuclei that form in different regions of the film, and observe that average asphericity increases closer to the interface, which is consistent with an increase in the free energy due to increased surface-to-volume ratios.

preprint2014arXiv

The Effect of Substrate on Thermodynamic and Kinetic Anisotropies in Atomic Thin Films

Glasses have a wide range of technological applications. The recent discovery of ultrastable glasses that are obtained by depositing the vapor of a glass-forming liquid onto the surface of a cold substrate has sparked renewed interest in the effects of confinements on physicochemical properties of liquids and glasses. Here we use molecular dynamics simulations to study the effect of substrate on thin films of a model glass-forming liquid, the Kob-Andersen binary Lennard-Jones system, and compute profiles of several thermodynamic and kinetic properties across the film. We observe that the substrate can induce large oscillations in profiles of thermodynamic properties such as density, composition and stress, and we establish a correlation between the oscillations in total density and the oscillations in normal stress. We also demonstrate that the kinetic properties of an atomic film can be readily tuned by changing the strength of interactions between the substrate and the liquid. Most notably, we show that a weakly attractive substrate can induce the emergence of a highly mobile region in its vicinity. In this highly mobile region, structural relaxation is several times faster than in the bulk, and the exploration of the potential energy landscape is also more efficient. In the subsurface region near a strongly attractive substrate, however, the dynamics is decelerated and the sampling of the potential energy landscape becomes less efficient than the bulk. We explain these two distinct behaviors by establishing a correlation between the oscillations in kinetic properties and the oscillations in lateral stress. Our findings offer interesting opportunities for designing better substrates for the vapor deposition process or developing alternative procedures for situations where vapor deposition is not feasible.

preprint2014arXiv

Two-state thermodynamics of the ST2 model for supercooled water

Thermodynamic properties of the ST2 model for supercooled liquid water exhibit anomalies similar to those observed in real water. A possible explanation of these anomalies is the existence of a metastable, liquid-liquid transition terminated by a critical point. This phenomenon, whose possible existence in real water is the subject of much current experimental work, has been unambiguously demonstrated for this particular model by most recent simulations. In this work, we reproduce the anomalies of two versions of the ST2 model with an equation of state describing water as a non-ideal "mixture" of two different types of local molecular order. We show that the liquid-liquid transition in the ST2 water is energy-driven. This is in contrast to another popular model, mW, in which non-ideality in mixing of two alternative local molecular orders is entropy-driven, and is not sufficiently strong to induce a liquid-liquid transition.

preprint2012arXiv

Temperature and Length Scale Dependence of Solvophobic Solvation in a Single-site Water-like Liquid

The temperature and length scale dependence of solvation properties of spherical hard solvophobic solutes is investigated in the Jagla liquid, a simple liquid that consists of particles interacting via a spherically symmetric potential combining a hard core repulsion and a longer ranged soft core interaction, yet exhibits water-like anomalies. The results are compared with equivalent calculations for a model of a typical atomic liquid, the Lennard-Jones (LJ) potential, and with predictions for hydrophobic solvation in water using the cavity equation of state and the extended simple point charge (SPC/E) model. We find that the Jagla liquid captures the qualitative thermodynamic behavior of hydrophobic hydration as a function of temperature for both small and large length scale solutes. In particular, for both the Jagla liquid and water, we observe temperature-dependent enthalpy and entropy of solvation for all solute sizes as well as a negative solvation entropy for sufficiently small solutes at low temperature. The results suggest that, compared to a simple liquid, it is the presence of a second thermally accessible repulsive energy scale, acting to increasingly favor larger separations for decreasing temperature, that is the essential characteristic of a liquid that favors low-density, open structures and models hydrophobic hydration, and that it is the presence of this second energy scale that leads to the similarity in the behavior of water and the Jagla liquid. The implications of the temperature and length scale dependence of solvation free energies in water-like liquids are explored with a simple model for the aggregation of solvophobic solutes. We show how aggregate stability depends upon the size of the aggregate and the size of its constituent solutes, and we relate this dependence to cold-induced destabilization phenomena such as the cold-induced denaturation of proteins.

preprint2011arXiv

Evaporation Rate of Water in Hydrophobic Confinement

The drying of hydrophobic cavities is believed to play an important role in biophysical phenomena such as the folding of globular proteins, the opening and closing of ligand-gated ion channels, and ligand binding to hydrophobic pockets. We use forward flux sampling, a molecular simulation technique, to compute the rate of capillary evaporation of water confined between two hydrophobic surfaces separated by nanoscopic gaps, as a function of gap, surface size and temperature. Over the range of conditions investigated (gaps between 9 and 14 Å and surface areas between 1 and 9 nm^2) the free energy barrier to evaporation scales linearly with the gap between hydrophobic surfaces, suggesting that line tension makes the predominant contribution to the free energy barrier. The exponential dependence of the evaporation rate on the gap between confining surfaces causes a ten order-of-magnitude decrease in the rate when the gap increases from 9 to 14 Å. The computed free energy barriers are of the order of 50kT, and are predominantly enthalpic. Evaporation rates per unit area are found to be two orders of magnitude faster in confinement by the larger (9 nm^2) than by the smaller (1nm^2) surfaces considered here, at otherwise identical conditions. We show that this is a direct consequence of the dependence of hydrophobic hydration on the size of solvated objects. For sufficiently large surfaces, the critical nucleus for the evaporation process is a gap-spanning cylindrical vapor tube.

preprint2011arXiv

Modeling Amphiphilic Solutes in a Jagla Solvent

Methanol is an amphiphilic solute whose aqueous solutions exhibit distinctive physical properties. The volume change upon mixing, for example, is negative across the entire composition range, indicating strong association. We explore the corresponding behavior of a Jagla solvent, which has been previously shown to exhibit many of the anomalous properties of water. We consider two models of an amphiphilic solute: (i) a "dimer" model, which consists of one hydrophobic hard sphere linked to a Jagla particle with a permanent bond, and (ii) a "monomer" model, which is a limiting case of the dimer, formed by concentrically overlapping a hard sphere and a Jagla particle. Using discrete molecular dynamics, we calculate the thermodynamic properties of the resulting solutions. We systematically vary the set of parameters of the dimer and monomer models and find that one can readily reproduce the experimental behavior of the excess volume of the methanolwater system as a function of methanol volume fraction. We compare the pressure and temperature dependence of the excess volume and the excess enthalpy of both models with experimental data on methanol-water solutions and find qualitative agreement in most cases. We also investigate the solute effect on the temperature of maximum density and find that the effect of concentration isorders of magnitude stronger than measured experimentally.

preprint2006arXiv

Effect of pressure on the phase behavior and structure of water confined between nanoscale hydrophobic and hydrophilic plates

We perform systematic molecular dynamics simulations of water confined between two nanoscale plates at T=300 K. We investigate the effect of pressure (-0.15 GPa <= P <= 0.2 GPa) and plate separation (0.4 nm <= d <= 1.6 nm) on the phase behavior of water when the plates are either hydrophobic or hydrophilic. When water is confined between hydrophobic plates, capillary evaporation occurs between the plates at low enough P. The threshold value of d at which this transition occurs decreases with P (e.g., 1.6 nm at P ~ -0.05 GPa; 0.5 nm at P \~ 0.1 GPa), until, at high P, no capillary evaporation occurs. For d ~ 0.6 nm and P >= 0.1 GPa, the system crystallizes into a bilayer ice. A P-d phase diagram showing the vapor, liquid, and bilayer ice phases is proposed. When water is confined by hydrophilic (hydroxylated silica) plates, it remains in the liquid phase at all P and d studied. Interestingly, we observe, for this case, that even at the P at which bulk water cavitates, the confined water remains in the liquid state. We also study systematically the state of hydration at different P for both kinds of plates. For the range of conditions studied here, we find that in the presence of hydrophobic plates the effect of P is to enhance water structure and to push water molecules toward the plates. The average orientation of water molecules next to the hydrophobic plates does not change upon pressurization. In contrast, in the presence of hydrophilic plates, water structure is insensitive to P. Hence, our results suggest that upon pressurization, hydrophobic plates behave as `soft' surfaces (in the sense of accommodating pressure-dependent changes in water structure) while hydrophilic walls behave as `hard' surfaces.