Researcher profile

Stephen J. Cox

Stephen J. Cox contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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.

Published work

7 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.

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.

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.

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.

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.

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.

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.