Source author record

Richard K. Bowles

Richard K. Bowles 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

15works
4topics
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

15 published item(s)

preprint2023arXiv

Is there a Granular Potential?

Granular materials, such as sand or grain, exhibit many structural and dynamic characteristics similar to those observed in molecular systems, despite temperature playing no role in their properties. This has led to an effort to develop a statistical mechanics for granular materials that has focused on establishing an equivalent to the microcanonical ensemble and a temperature-like thermodynamic variable. Here, we expand on these ideas by introducing a granular potential into the Edwards ensemble, as an analogue to the chemical potential, and explore its properties using a simple model of a granular system. A simple kinetic Monte Carlo simulation of the model shows the effect of mass transport leading to equilibrium and how this is connected to the redistribution of volume in the system. An exact analytical treatment of the model shows that the compactivity and the ratio of the granular potential to the compactivity determine the equilibrium between two open systems that are able to exchange volume and particles, and that mass moves from high to low values of this ratio. Analysis of the granular potential shows that adding a particle to the system increases the entropy at high compactivity, but decreases the entropy at low compactivity. Finally, we demonstrate the use of a small system thermodynamics method for the measurement of granular potential differences.

preprint2021arXiv

The Influence on Crystal Nucleation of an Order-Disorder Transition among the Subcritical Clusters

Studies of nucleation generally focus on the properties of the critical cluster, but the presence of defects within the crystal lattice means that the population of nuclei necessarily evolve through a distribution of pre-critical clusters with varying degrees of structural disorder on their way to forming a growing stable crystal. To investigate the role pre-critical clusters play in nucleation, we develop a simple thermodynamic model for crystal nucleation in terms of cluster size and the degree of cluster order that allows us to alter the work of forming the pre-critical clusters without effecting the properties of the critical cluster. The steady state and transient nucleation behaviour of the system are then studied numerically, for different microscopic ordering kinetics. We find that the models exhibits a generic order-disorder transition in the pre-critical clusters. Independent of the type of ordering kinetics, increasing the accessibility of disordered pre-critical clusters decreases both the steady state nucleation rate and the nucleation lag time. Furthermore, the interplay between the free energy surface and the microscopic ordering kinetics leads to three distinct nucleation pathways.

preprint2020arXiv

Mapping diffusivity of narrow channels into one-dimension

The diffusion of particles trapped in long narrow channels occurs predominantly in one dimension. Here, molecular dynamics simulation is used to study the inertial dynamics of two-dimensional hard disks, confined to long, narrow, structureless channels with hard walls in the no-passing regime. We show that the diffusion coefficient obtained from the mean squared displacement can be mapped onto the exact results for the diffusion of the strictly one-dimensional hard rod system through an effective occupied volume fraction obtained from either the equation of state or a geometric projection of the particle interaction diameters.

preprint2016arXiv

Competitive Heterogeneous Nucleation onto a Microscopic Impurity in a Potts Model

Many metastable systems can nucleate to multiple competing stable or intermediate metastable states. In this work, a Potts model, subject to external fields, is used to study the competitive nucleation of two phases attempting to grow on a microscopic impurity. Monte Carlo simulations are used to calculate the free energy surfaces for the system under different conditions, where the relative stability of the phases is adjusted by changing the interaction parameters, and the nucleation rates obtained using multicomponent transition state theory are compared with the rates measured using the survival probability method. We find that the two methods predict similar nucleation rates when the free energy barrier used in the transition state theory is defined as the work required to form a critical embryo from the metastable phase. An analysis of the free energy surfaces also reveals that the competition between the nucleating phases leads to an effective drying of the impurity which slows down the nucleation rate compared to the single phase case.

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.

preprint2014arXiv

An orientationally ordered helical fluid phase in a quasi-one-dimensional system of confined hard spheres

We use a series of molecular dynamics simulations, and analytical theory, to demonstrate that a system of hard spheres confined to a narrow cylindrical channel exhibits a continuous phase transition from an isotropic fluid at low densities, to an orientationally ordered, but translationally disordered, helical fluid at high densities. The ordered fluid phase contains small sections of helix separated by topological defects that change the direction of the twist, altering the local chirailty. The defects break up the translational order, but the fluid develops long range orientational order. An analysis of the particle packings show that the length separation between defects controls the geometrical properties of the helical sections, including the orientation, and that pairs of defects experience a weak, but long range attraction resulting from entropic free volumes effects. These collective long range interactions overcome the restrictions on quasi-one-dimensional transitions, even though the particle-particle interaction is short ranged.

preprint2013arXiv

A Transition State Theory for Calculating Hopping Times and Diffusion in Highly Confined Fluids

Monte Carlo simulation is used to study the dynamical crossover from single file diffusion to normal diffusion in fluids confined to narrow channels. We show that the long time diffusion coefficients for a series of systems involving hard and soft interaction potentials can be described in terms of a hopping time that measures the time it takes for a particle to escape the cage formed by its neighbors in the pore. Free energy barriers for the particle hopping process are calculated and used to show that transition state theory effectively describes the hopping time for all the systems studied, over a range of pore diameters. Our work suggests that the combination of hopping times and transition state theory offers a useful and general framework to describe the dynamics of these highly confined fluids.

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

A Fragile-Strong Fluid Crossover and Universal Relaxation Times in a Confined Hard Disc Fluid

We show that a system of hard discs confined to a narrow channel exhibits a fragile-strong fluid crossover located at the maximum of the isobaric heat capacity and that the relaxation times for different channel widths fall onto a single master curve when rescaled by the relaxation times and temperatures of the crossover. Calculations of the configurational entropy and the inherent structure equation of state find that the crossover is related to properties of the jamming landscape for the model but that the Adams-Gibbs relation does not predict the relaxation behavior. We also show that a facilitated dynamics description of the system, where kinetically excited regions are identified with local packing arrangements of the discs, successfully describes the fragile-strong crossover.

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

Edwards entropy and compactivity in a model of granular matter

Formulating a statistical mechanics for granular matter remains a significant challenge, in part, due to the difficulty associated with a complete characterization of the systems under study. We present a fully characterized model of a granular material consisting of $N$ two-dimensional, frictionless, hard discs, confined between hard walls, including a complete enumeration of all possible jammed structures. We show the properties of the jammed packings are independent of the distribution of defects within the system and that all the packings are isostatic. This suggests the assumption of equal probability for states of equal volume, which provides one possible way of constructing the equivalent of a microcanonical ensemble, is likely to be vaild for our model. An application of the second law of thermodynamics involving two subsystems in contact shows that the expected spontaneous equilibration of defects between the two is accompanied by an increase in entropy and that the equilibirum, obtained by entropy maximization, is characterized by the equality of compactivities. Finally, we explore the properties of the equivalent to the canonical ensemble for this system.

preprint2011arXiv

Heterogeneous condensation of the Lennard-Jones vapor onto a nanoscale seed particle

The heterogeneous condensation of a Lennard-Jones vapor onto a nanoscale seed particle is studied using molecular dynamics simulations. Measuring the nucleation rate and the height of the free energy barrier using the mean first passage time method shows that the presence of a weakly interacting seed has little effect on the work of forming very small cluster embryos but accelerates the rate by lowering the barrier for larger clusters. We suggest that this results from a competition between the energetic and entropic features of cluster formation in the bulk and at the heterogeneity. As the interaction is increased, the free energy of formation is reduced for all cluster sizes. We also develop a simple phenomenological model of film formation on a small seed that captures the general features of the nucleation process for small heterogeneities. A comparison of our simulation results with the model shows that heterogeneous classical nucleation theory provides a good estimate of the critical size of the film but significantly over-estimates the size of the barrier.

preprint2010arXiv

The role of FCC tetrahedral subunits in the phase behavior of medium sized Lennard-Jones clusters

The free energy of a 600-atom Lennard-Jones cluster is calculated as a function of surface and bulk crystallinity in order to study the structural transformations that occur in the core of medium sized clusters. Within the order parameter range studied, we find the existence of two free energy minima at temperatures near freezing. One minimum, at low values of both bulk and surface order, belongs to the liquid phase. The second minimum exhibits a highly ordered core with a disordered surface and is related to structures containing a single FCC-tetrahedral subunit, with an edge length of seven atoms (l=7), located in the particle core. At lower temperatures, a third minimum appears at intermediate values of the bulk order parameter which is shown to be related to the formation of multiple l=6 tetrahedra in the core of the cluster. We also use molecular dynamics simulations to follow a series of nucleation events and find that the clusters freeze to structures containing l=5,6,7 and 8 sized tetrahedra as well as those containing no tetrahedral units. The structural correlations between bulk and surface order with the size of the tetrahedral units in the cluster core are examined. Finally, the relationships between the formation of FCC tetrahedral subunits in the core, the phase behavior of medium sized clusters and the nucleation of noncrystalline global structures such as icosahedra and decahedra are discussed.

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.

preprint2009arXiv

The Complete Jamming Landscape of Confined Hard Discs

An exact description of the complete jamming landscape is developed for a system of hard discs of diameter $σ$, confined between two lines separated by a distance $1+\sqrt{3/4} < H/σ< 2$. By considering all possible local packing arrangements, the generalized ensemble partition function of jammed states is obtained using the transfer matrix method, which allows us to calculate the configurational entropy and the equation of state for the packings. Exploring the relationship between structural order and packing density, we find that the geometric frustration between local packing environments plays an important role in determining the density distribution of jammed states and that structural "randomness" is a non-monotonic function of packing density. Molecular dynamics simulations show that the properties of the equilibrium liquid are closely related to those of the landscape.