BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stephen J. Cox

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.stat-mechcond-mat.softcond-mat.dis-nnphysics.chem-ph

Catalog footprint

What is connected

12works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

12 published item(s)

preprint2023arXiv

Classical quantum friction at water-carbon interfaces

Friction at water-carbon interfaces remains a major puzzle with theories and simulations unable to explain experimental trends in nanoscale waterflow. A recent theoretical framework -- quantum friction (QF)-- proposes to resolve these experimental observations by considering nonadiabatic coupling between dielectric fluctuations in water and graphitic surfaces. Here, using a classical model that enables fine-tuning of the solid's dielectric spectrum, we provide evidence from simulations in general support of QF. In particular, as features in the solid's dielectric spectrum begin to overlap with water's librational and Debye modes, we find an increase in friction in line with that proposed by QF. At the microscopic level, we find that this contribution to friction manifests more distinctly in the dynamics of the solid's charge density than that of water. Our findings suggest that experimental signatures of QF may be more pronounced in the solid's response rather than liquid water's.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A theory for the stabilization of polar crystal surfaces by a liquid environment

Polar crystal surfaces play an important role in the functionality of many materials, and have been studied extensively over many decades. In this article, a theoretical framework is presented that extends existing theories by placing the surrounding solution environment on an equal footing with the crystal itself; this is advantageous, e.g., when considering processes such as crystal growth from solution. By considering the polar crystal as a stack of parallel plate capacitors immersed in a solution environment, the equilibrium adsorbed surface charge density is derived by minimizing the free energy of the system. In analogy to the well-known diverging surface energy of a polar crystal surface at zero temperature, for a crystal in solution it is shown that the "polar catastrophe" manifests as a diverging free energy cost to perturb the system from equilibrium. Going further than existing theories, the present formulation predicts that fluctuations in the adsorbed surface charge density become increasingly suppressed with increasing crystal thickness. We also show how, in the slab geometry often employed in both theoretical and computational studies of interfaces, an electric displacement field emerges as an electrostatic boundary condition, the origins of which are rooted in the slab geometry itself, rather than the use of periodic boundary conditions. This aspect of the work provides a firmer theoretical basis for the recent observation that standard "slab corrections" fail to correctly describe, even qualitatively, polar crystal surfaces in solution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation?

In principle, the answer to the posed titular question is undoubtedly 'yes.' But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the 'cut-and-shift' version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice I$_{\rm h}$. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dielectric response of thin water films: A thermodynamic perspective

The surface of a polar liquid presents a special environment for the solvation and organization of charged solutes, which differ from bulk behaviors in important ways. These differences have motivated many attempts to understand electrostatic response at aqueous interfaces in terms of a spatially varying dielectric permittivity, typically concluding that the dielectric constant of interfacial water is significantly lower than in the bulk liquid. Such analyses, however, are complicated by the potentially nonlocal nature of dielectric response over the short length scales of interfacial heterogeneity. Here we circumvent this problem for thin water films by adopting a thermodynamic approach. Using molecular simulations, we calculate the solvent's contribution to the reversible work of charging a parallel plate capacitor. We find good agreement with a simple dielectric continuum model that assumes bulk dielectric permittivity all the way up to the liquid's boundary, even for very thin ($\sim 1$ nm) films. This comparison requires careful attention to the placement of dielectric boundaries between liquid and vapor, which also resolves apparent discrepancies with dielectric imaging experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How do interfaces alter the dynamics of supercooled water?

The structure of liquid water in the proximity of an interface can deviate significantly from that of bulk water, with surface-induced structural perturbations typically converging to bulk values at about ~1 nm from the interface. While these structural changes are well established it is, in contrast, less clear how an interface perturbs the dynamics of water molecules within the liquid. Here, through an extensive set of molecular dynamics simulations of supercooled bulk and interfacial water films and nano-droplets, we observe the formation of persistent, spatially extended dynamical domains in which the average mobility varies as a function of the distance from the interface. This is in stark contrast with the dynamical heterogeneity observed in bulk water, where these domains average out spatially over time. We also find that the dynamical response of water to an interface depends critically on the nature of the interface and on the choice of interface definition. Overall these results reveal a richness in the dynamics of interfacial water that opens up the prospect of tuning the dynamical response of water through specific modifications of the interface structure or confining material.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite field formalism for bulk electrolyte solutions

The manner in which electrolyte solutions respond to electric fields is crucial to understanding the behavior of these systems both at, and away from, equilibrium. The present formulation of linear response theory for such systems is inconsistent with common molecular dynamics (MD) implementations. Using the finite field formalism, suitably adapted for finite temperature MD, we investigate the response of bulk aqueous NaCl solutions to both finite Maxwell ($\mathbf{E}$) and electric displacement ($\mathbf{D}$) fields. The constant $\mathbf{E}$ Hamiltonian allows us to derive the linear response relation for the ionic conductivity in a simple manner that is consistent with the forces used in conventional MD simulations. Simulations of a simple point charge model of an electrolyte solution at constant $\mathbf{E}$ yield conductivities at infinite dilution within 15% of experimental values. The finite field approach also allows us to measure the solvent's dielectric constant from its polarization response, which is seen to decrease with increasing ionic strength. Comparison of the dielectric constant measured from polarization response versus polarization fluctuations enables direct evaluation of the dynamic contribution to this dielectric decrement, which we find to be small but not insignificant. Using the constant $\mathbf{D}$ formulation, we also rederive the Stillinger-Lovett conditions, which place strict constraints on the coupling between solvent and ionic polarization fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation

Silver iodide is one of the most potent inorganic ice nucleating particles known, a feature generally attributed to the excellent lattice match between its basal Ag-(0001) and I-(000-1) surfaces, and ice. This crystal termination, however, is a type-III polar surface, and its surface energy therefore diverges with crystal size unless a polarity compensation mechanism prevails. In this simulation study, we investigate to what extent the surrounding aqueous environment is able to provide such polarity compensation. On its own, we find that pure water is unable to stabilize the AgI crystal in a physically reasonable manner, and that mobile charge carriers such as dissolved ions, are essential. In other words, proximate dissolved ions must be considered an integral part of the heterogeneous ice formation mechanism. The simulations we perform utilize recent advances in simulation methodology in which appropriate electric and electric displacement fields are imposed. A useful by-product of this study is the direct comparison to the commonly used Yeh-Berkowitz method that this enables. Here we find that naive application of the latter leads to physically unreasonable results, and greatly influences the structure of water in the contact layer. We therefore expect these results to be of general importance to those studying polar/charged surfaces in aqueous environments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insight into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that in the last few decades have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state of the art computational methods, by reviewing simulations of e.g. ice nucleation or crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insight into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that in doing so the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that by improving (i.) existing interatomic potentials; and (ii.) currently available enhanced sampling methods, the community can move towards accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular simulations of heterogeneous ice nucleation. I. Controlling ice nucleation through surface hydrophilicity

Ice formation is one of the most common and important processes on earth and almost always occurs at the surface of a material. A basic understanding of how the physicochemical properties of a material's surface affect its ability to form ice has remained elusive. Here, we use molecular dynamics simulations to directly probe heterogeneous ice nucleation at a hexagonal surface of a nanoparticle of varying hydrophilicity. Surprisingly, we find that structurally identical surfaces can both inhibit and promote ice formation and analogous to a chemical catalyst, it is found that an optimal interaction between the surface and the water exists for promoting ice nucleation. We use our microscopic understanding of the mechanism to design a modified surface in silico with enhanced ice nucleating ability.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular simulations of heterogeneous ice nucleation. II. Peeling back the layers

Coarse grained molecular dynamics simulations are presented in which the sensitivity of the ice nucleation rate to the hydrophilicity of a graphene nanoflake is investigated. We find that an optimal interaction strength for promoting ice nucleation exists, which coincides with that found previously for a face centered cubic (111) surface. We further investigate the role that the layering of interfacial water plays in heterogeneous ice nucleation and demonstrate that the extent of layering is not a good indicator of ice nucleating ability for all surfaces. Our results suggest that to be an efficient ice nucleating agent, a surface should not bind water too strongly if it is able to accommodate high coverages of water.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Benchmarking the performance of Density Functional Theory and Point Charge Force Fields in their Description of sI Methane Hydrate against Diffusion Monte Carlo

High quality reference data from diffusion Monte Carlo calculations are presented for bulk sI methane hydrate, a complex crystal exhibiting both hydrogen-bond and dispersion dominated interactions. The performance of some commonly used exchange-correlation functionals and all-atom point charge force fields is evaluated. Our results show that none of the exchange-correlation functionals tested are sufficient to describe both the energetics and the structure of methane hydrate accurately, whilst the point charge force fields perform badly in their description of the cohesive energy but fair well for the dissociation energetics. By comparing to ice Ih, we show that a good prediction of the volume and cohesive energies for the hydrate relies primarily on an accurate description of the hydrogen bonded water framework, but that to correctly predict stability of the hydrate with respect to dissociation to ice Ih and methane gas, accuracy in the water-methane interaction is also required. Our results highlight the difficulty that density functional theory faces in describing both the hydrogen bonded water framework and the dispersion bound methane.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Microscopic Features of Heterogeneous Ice Nucleation May Affect the Macroscopic Morphology of Atmospheric Ice Crystals

It is surprisingly difficult to freeze water. Almost all ice that forms under "mild" conditions (temperatures > -40 degrees Celsius) requires the presence of a nucleating agent - a solid particle that facilitates the freezing process - such as clay mineral dust, soot or bacteria. In a computer simulation, the presence of such ice nucleating agents does not necessarily alleviate the difficulties associated with forming ice on accessible timescales. Nevertheless, in this work we present results from molecular dynamics simulations in which we systematically compare homogeneous and heterogeneous ice nucleation, using the atmospherically important clay mineral kaolinite as our model ice nucleating agent. From our simulations, we do indeed find that kaolinite is an excellent ice nucleating agent but that contrary to conventional thought, non-basal faces of ice can nucleate at the basal face of kaolinite. We see that in the liquid phase, the kaolinite surface has a drastic effect on the density profile of water, with water forming a dense, tightly bound first contact layer. Monitoring the time evolution of the water density reveals that changes away from the interface may play an important role in the nucleation mechanism. The findings from this work suggest that heterogeneous ice nucleating agents may not only enhance the ice nucleation rate, but also alter the macroscopic structure of the ice crystals that form.

0 citations0 reviewsNo code linkedOpen access