Researcher profile

Peter H. Poole

Peter H. Poole contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

cond-mat.stat-mech cond-mat.mtrl-sci cond-mat.soft cond-mat.dis-nn physics.chem-ph

Trust snapshot

Quick read

Trust 21 - Emerging

15works

0followers

5topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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.

preprint2016arXiv

Potential Energy Landscape of the Apparent First-Order Phase Transition between Low-Density and High-Density Amorphous Ice

The potential energy landscape (PEL) formalism is a valuable approach within statistical mechanics for describing supercooled liquids and glasses. Here we use the PEL formalism and computer simulations to study the pressure-induced transformations between low-density amorphous ice (LDA) and high-density amorphous ice (HDA) at different temperatures. We employ the ST2 water model for which the LDA-HDA transformations are remarkably sharp, similar to what is observed in experiments, and reminiscent of a first-order phase transition. Our results are consistent with the view that LDA and HDA configurations are associated with two distinct regions (megabasins) of the PEL that are separated by a potential energy barrier. At higher temperature, we find that low-density liquid (LDL) configurations are located in the same megabasin as LDA, and that high-density liquid (HDL) configurations are located in the same megabasin as HDA. We show that the pressure-induced LDL-HDL and LDA-HDA transformations occur along paths that interconnect these two megabasins, but that the path followed by the liquid is different than the path followed by the amorphous solid. At higher pressure, we also study the liquid-to-ice-VII first-order phase transition, and find that the behavior of the PEL properties across this transition are qualitatively similar to the changes found during the LDA-HDA transformation. This similarity supports the interpretation that the LDA-HDA transformation is a first-order-like phase transition between out-of-equilibrium states.

preprint2015arXiv

Free energy of formation of small ice nuclei near the Widom line in simulations of supercooled water

The ST2 interaction potential has been used in a large number of simulation studies to explore the possibility of a liquid-liquid phase transition (LLPT) in supercooled water. Using umbrella sampling Monte Carlo simulations of ST2 water, we evaluate the free energy of formation of small ice nuclei in the supercooled liquid in the vicinity of the Widom line, the region above the critical temperature of the LLPT where a number of thermodynamic anomalies occur. Our results show that in this region there is a substantial free energy cost for the formation of small ice nuclei, demonstrating that the thermodynamic anomalies associated with the Widom line in ST2 water occur in a well-defined metastable liquid phase. On passing through the Widom line, we identify changes in the free energy to form small ice nuclei that illustrate how the thermodynamic anomalies associated with the LLPT may influence the ice nucleation process.

preprint2015arXiv

Phase diagram of the ST2 model of water

We evaluate the free energy of the fluid and crystal phases for the ST2 potential [F.H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)] with reaction field corrections for the long-range interactions. We estimate the phase coexistence boundaries in the temperature-pressure plane, as well as the gas-liquid critical point and gas-liquid coexistence conditions. Our study frames the location of the previously identified liquid-liquid critical point relative to the crystalline phase boundaries, and opens the way for exploring crystal nucleation in a model where the metastable liquid-liquid critical point is computationally accessible.

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.

preprint2013arXiv

Heterogeneous Nucleation in the Low Barrier Regime

In simulations of the 2D Ising model, we examine heterogeneous nucleation induced by a small impurity consisting of a line of $l$ fixed spins. As $l$ increases, we identify a limit of stability beyond which the metastable phase is not defined. We evaluate the free energy barrier for nucleation of the stable phase, and show that contrary to expectation, the barrier does not vanish on approach to the limit of stability. We also demonstrate that our values for the height of the barrier yield predictions for the nucleation time (from transition state theory) and the size of the critical cluster (from the nucleation theorem) that are in excellent agreement with direct measurements, even near the limit of stability.

preprint2012arXiv

Free energy surface of ST2 water near the liquid-liquid phase transition

We carry out umbrella sampling Monte Carlo simulations to evaluate the free energy surface of the ST2 model of water as a function two order parameters, the density and a bond-orientational order parameter. We approximate the long-range electrostatic interactions of the ST2 model using the reaction-field method. We focus on state points in the vicinity of the liquid-liquid critical point proposed for this model in earlier work. At temperatures below the predicted critical temperature we find two basins in the free energy surface, both of which have liquid-like bond orientational order, but differing in density. The pressure and temperature dependence of the shape of the free energy surface is consistent with the assignment of these two basins to the distinct low density and high density liquid phases previously predicted to occur in ST2 water.

preprint2011arXiv

Dynamical Behavior Near a Liquid-Liquid Phase Transition in Simulations of Supercooled Water

We examine the behavior of the diffusion coefficient of the ST2 model of water over a broad region of the phase diagram via molecular dynamics simulations. The ST2 model has an accessible liquid-liquid transition between low-density and high-density phases, making the model an ideal candidate to explore the impacts of the liquid-liquid transition on dynamics. We locate characteristic dynamical loci in the phase diagram and compare them with the previously investigated thermodynamic loci. The low-density liquid phase shows a crossover from non-Arrhenius to Arrhenius behavior, signaling the onset of a crossover from fragile-to-strong behavior. We explain this crossover in terms of the asymptotic approach of the low-density liquid to a random tetrahedral network, and show that the temperature dependence of the diffusion coefficient over a wide temperature range can be simply related to the concentration of defects in the network. Our findings thus confirm that the low-density phase of ST2 water is a well-defined metastable liquid.

preprint2011arXiv

Mixture-like behavior near a liquid-liquid phase transition in simulations of supercooled water

In simulations of a water-like model (ST2) that exhibits a liquid-liquid phase transition, we test for the occurrence of a thermodynamic region in which the liquid can be modelled as a two-component mixture. We assign each molecule to one of two species based on the distance to its fifth-nearest neighbor, and evaluate the concentration of each species over a wide range of temperature and density. Our concentration data compare well with mixture-model predictions in a region between the liquid-liquid critical temperature and the temperature of maximum density. Fits of the model to the data in this region yield accurate estimates for the location of the critical point. We also show that the liquid outside the region of density anomalies is poorly modelled as a simple mixture.

preprint2010arXiv

Simulations of a lattice model of two-headed linear amphiphiles: influence of amphiphile asymmetry

Using a 2D lattice model, we conduct Monte Carlo simulations of micellar aggregation of linear-chain amphiphiles having two solvophilic head groups. In the context of this simple model, we quantify how the amphiphile architecture influences the critical micelle concentration (CMC), with a particular focus on the role of the asymmetry of the amphiphile structure. Accordingly, we study all possible arrangements of the head groups along amphiphile chains of fixed length $N=12$ and 16 molecular units. This set of idealized amphiphile architectures approximates many cases of symmetric and asymmetric gemini surfactants, double-headed surfactants and boloform surfactants. Consistent with earlier results, we find that the number of spacer units $s$ separating the heads has a significant influence on the CMC, with the CMC increasing with $s$ for $s<N/2$. In comparison, the influence of the asymmetry of the chain architecture on the CMC is much weaker, as is also found experimentally.

preprint2010arXiv

Spatial Correlation of the Dynamic Propensity in a Glass-Forming Liquid

We present computer simulation results on the dynamic propensity [as defined by Widmer- Cooper, Harrowell, and Fynewever, Phys. Rev. Lett. 93, 135701 (2004)] in a Kob-Andersen binary Lennard-Jones liquid system consisting of 8788 particles. We compute the spatial correlation function of the dynamic propensity as a function of both the reduced temperature T, and the time scale on which the particle displacements are measured. For T\leq0.6, we find that nonzero correlations occur at the largest length scale accessible in our system. We also show that a cluster-size analysis of particles with extremal values of the dynamic propensity, as well as 3D visualizations, reveal spatially correlated regions that approach the size of our system as T decreases, consistent with the behavior of the spatial correlation function. Next, we define and examine the "coordination propensity", the isoconfigurational average of the coordination number of the minority B particles around the majority A particles. We show that a significant correlation exists between the spatial fluctuations of the dynamic and coordination propensities. In addition, we find non-zero correlations of the coordination propensity occurring at the largest length scale accessible in our system for all T in the range 0.466<T<1.0. We discuss the implications of these results for understanding the length scales of dynamical heterogeneity in glass-forming liquids.

preprint2009arXiv

Crystal Nucleation in a Supercooled Liquid with Glassy Dynamics

In simulations of supercooled, high-density liquid silica we study a range of temperature T in which we find both crystal nucleation, as well as the characteristic dynamics of a glass forming liquid, including a breakdown of the Stokes-Einstein relation. We find that the liquid cannot be observed below a homogeneous nucleation limit (HNL) at which the liquid crystallizes faster than it can equilibrate. We show that the HNL would occur at lower T, and perhaps not at all, if the Stokes-Einstein relation were obeyed, and hence that glassy dynamics plays a central role in setting a crystallization limit on the liquid state in this case. We also explore the relation of the HNL to the Kauzmann temperature, and test for spinodal-like effects near the HNL.

preprint2005arXiv

Density minimum and liquid-liquid phase transition

We present a high-resolution computer simulation study of the equation of state of ST2 water, evaluating the liquid-state properties at 2718 state points, and precisely locating the liquid-liquid critical point (LLCP) occurring in this model. We are thereby able to reveal the interconnected set of density anomalies, spinodal instabilities and response function extrema that occur in the vicinity of a LLCP for the case of a realistic, off-lattice model of a liquid with local tetrahedral order. In particular, we unambiguously identify a density minimum in the liquid state, define its relationship to other anomalies, and show that it arises due to the approach of the liquid structure to a defect-free random tetrahedral network of hydrogen bonds.

preprint2004arXiv

Phase diagram of silica from computer simulation

We evaluate the phase diagram of the ``BKS'' potential [Van Beest, Kramer and van Santen, Phys. Rev. Lett. 64, 1955 (1990)], a model of silica widely used in molecular dynamics (MD) simulations. We conduct MD simulations of the liquid, and three crystals (beta-quartz, coesite and stishovite) over wide ranges of temperature and density, and evaluate the total Gibbs free energy of each phase. The phase boundaries are determined by the intersection of these free energy surfaces. Not unexpectedly for a classical pair potential, our results reveal quantitative discrepancies between the locations of the BKS and real silica phase boundaries. At the same time, we find that the topology of the real phase diagram is reproduced, confirming that the BKS model provides a satisfactory qualitative description of a silica-like material. We also compare the phase boundaries with the locations of liquid-state thermodynamic anomalies identified in previous studies of the BKS model.

preprint2003arXiv

Free energy and configurational entropy of liquid silica: fragile-to-strong crossover and polyamorphism

Recent molecular dynamics (MD) simulations of liquid silica, using the ``BKS'' model [Van Beest, Kramer and van Santen, Phys. Rev. Lett. {\bf 64}, 1955 (1990)], have demonstrated that the liquid undergoes a dynamical crossover from super-Arrhenius, or ``fragile'' behavior, to Arrhenius, or ``strong'' behavior, as temperature $T$ is decreased. From extensive MD simulations, we show that this fragile-to-strong crossover (FSC) can be connected to changes in the properties of the potential energy landscape, or surface (PES), of the liquid. To achieve this, we use thermodynamic integration to evaluate the absolute free energy of the liquid over a wide range of density and $T$. We use this free energy data, along with the concept of ``inherent structures'' of the PES, to evaluate the absolute configurational entropy $S_c$ of the liquid. We find that the temperature dependence of the diffusion coefficient and of $S_c$ are consistent with the prediction of Adam and Gibbs, including in the region where we observe the FSC to occur. We find that the FSC is related to a change in the properties of the PES explored by the liquid, specifically an inflection in the $T$ dependence of the average inherent structure energy. In addition, we find that the high $T$ behavior of $S_c$ suggests that the liquid entropy might approach zero at finite $T$, behavior associated with the so-called Kauzmann paradox. However, we find that the change in the PES that underlies the FSC is associated with a change in the $T$ dependence of $S_c$ that elucidates how the Kauzmann paradox is avoided in this system. Finally, we also explore the relation of the observed PES changes to the recently discussed possibility that BKS silica exhibits a liquid-liquid phase transition, a behavior that has been proposed to underlie the observed polyamorphism of amorphous solid silica.

preprint2002arXiv

Bulk motion of granular matter in an agitated cylindrical bed

Experimental results are reported for the bulk motion induced in a bed of granular matter contained in a cylindrical pan with a flat bottom subjected to simultaneous vertical and horizontal vibrations. The motion in space of the moving pan is quantified. A number of distinct bulk dynamical modes are observed in which the particle bed adopts different shapes and motions. At the lowest pan excitation frequency $ω$, the bed forms a ``heap,'' and rotates about the cylinder axis. As $ω$ is increased, a more complex ``toroidal'' mode appears in which the bed takes the shape of a torus; in this mode, circulation occurs both about the cylinder axis, and also radially, with particles moving from the outer edge of the pan to the centre on the top surface of the bed, and back to the outer edge along the pan bottom. At the highest $ω$, surface modulations (``surface waves'' and ``sectors'') of the toroidal mode occur. The origin of this family of behavior in terms of the pan motion is discussed.