1
|
Pireddu G, Fairchild CJ, Niblett SP, Cox SJ, Rotenberg B. Impedance of nanocapacitors from molecular simulations to understand the dynamics of confined electrolytes. Proc Natl Acad Sci U S A 2024; 121:e2318157121. [PMID: 38662549 PMCID: PMC11067016 DOI: 10.1073/pnas.2318157121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2023] [Accepted: 04/01/2024] [Indexed: 05/05/2024] Open
Abstract
Nanoelectrochemical devices have become a promising candidate technology across various applications, including sensing and energy storage, and provide new platforms for studying fundamental properties of electrode/electrolyte interfaces. In this work, we employ constant-potential molecular dynamics simulations to investigate the impedance of gold-aqueous electrolyte nanocapacitors, exploiting a recently introduced fluctuation-dissipation relation. In particular, we relate the frequency-dependent impedance of these nanocapacitors to the complex conductivity of the bulk electrolyte in different regimes, and use this connection to design simple but accurate equivalent circuit models. We show that the electrode/electrolyte interfacial contribution is essentially capacitive and that the electrolyte response is bulk-like even when the interelectrode distance is only a few nanometers, provided that the latter is sufficiently large compared to the Debye screening length. We extensively compare our simulation results with spectroscopy experiments and predictions from analytical theories. In contrast to experiments, direct access in simulations to the ionic and solvent contributions to the polarization allows us to highlight their significant and persistent anticorrelation and to investigate the microscopic origin of the timescales observed in the impedance spectrum. This work opens avenues for the molecular interpretation of impedance measurements, and offers valuable contributions for future developments of accurate coarse-grained representations of confined electrolytes.
Collapse
Affiliation(s)
- Giovanni Pireddu
- Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, CNRS, Sorbonne Université, Physicochimie des Électrolytes et Nanosystèmes Interfaciaux (PHENIX), CNRS, Sorbonne Université, ParisF-75005, France
| | - Connie J. Fairchild
- Yusuf Hamied Department of Chemistry, University of Cambridge, CambridgeCB2 1EW, United Kingdom
| | - Samuel P. Niblett
- Yusuf Hamied Department of Chemistry, University of Cambridge, CambridgeCB2 1EW, United Kingdom
| | - Stephen J. Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, CambridgeCB2 1EW, United Kingdom
| | - Benjamin Rotenberg
- Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, CNRS, Sorbonne Université, Physicochimie des Électrolytes et Nanosystèmes Interfaciaux (PHENIX), CNRS, Sorbonne Université, ParisF-75005, France
- Réseau sur le Stockage Electrochimique de l’Energie, Fédération de Recherche CNRS 3459, Amiens Cedex80039, France
| |
Collapse
|
2
|
Advincula XR, Backus EHG, Bonn M, Cox SJ, Diebold U, Fellows A, Finney AR, Goel G, Hedley J, Jiang Y, Jin D, Kapil V, Kavokine N, Klein J, Laage D, Mohandas N, Morgenstern K, Mukherjee T, Olvera de la Cruz M, Orlikowska-Rzeznik H, Perkin S, Piaggi PM, Rodellar CG, Ryan P, Sayer T, Seyffertitz M, Shepelenko M, Sosso GC, Thämer M, Vilangottunjalil A, Walker-Gibbons R, Wang Y, Willard AP, Zhang P. Electrified/charged aqueous interfaces: general discussion. Faraday Discuss 2024; 249:381-407. [PMID: 38170868 DOI: 10.1039/d3fd90065g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
|
3
|
Advincula XR, Blow KE, Bonn M, Bui AT, Cheng Y, Cox SJ, Della Pia F, Diebold U, Fumagalli L, Goel G, Hayton JA, Jiang Y, Kapil V, Kavokine N, Koga K, Laage D, Lahav M, Miao S, Michaelides A, Montero de Hijes P, Morgenstern K, Mukherjee T, O'Neill N, Pan D, Piaggi PM, Rempe SLB, Salvalaglio M, Salzmann CG, Sayer T, Shepelenko M, Sosso GC, Wang S, Webber B, Willard AP, Yao Y. Dynamics and nano-rheology of interfacial water: general discussion. Faraday Discuss 2024; 249:243-266. [PMID: 38174388 DOI: 10.1039/d3fd90064a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
|
4
|
Advincula XR, Backus EHG, Bartels-Rausch T, Benaglia S, Ben Ari G, Blow KE, Bonn M, Bui AT, Cox SJ, Della Pia F, Diebold U, Finney AR, Franceschi G, Fumagalli L, Goel G, Hayton JA, Holdship C, Jiang Y, Jin D, Kapil V, Kavokine N, Koga K, Laage D, Lahav M, Miao S, Michaelides A, Mohandas N, Morgenstern K, Mukherjee T, Nagata Y, Olvera de la Cruz M, Pan D, Piaggi PM, Rempe SLB, Ryan P, Salzmann CG, Sayer T, Saykally RJ, Shepelenko M, Sosso GC, Whale TF, White JJ, Willard AP, Zhang P. Ice interfaces: general discussion. Faraday Discuss 2024; 249:133-161. [PMID: 38174608 DOI: 10.1039/d3fd90063k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
|
5
|
Backus EHG, Ben Ari G, Benaglia S, Bonn M, Bui AT, Cox SJ, Della Pia F, Fraxedas J, Goel G, Jiang Y, Jin D, Koga K, Laage D, Miao S, Michaelides A, Morgenstern K, Mukherjee T, Nagata Y, Naito H, Nir O, Olvera de la Cruz M, Orlikowska-Rzeznik H, Pan D, Rempe SLB, Salzmann CG, Taira A, Vilangottunjalil A, Wang S, Willard AP, Yao Y, Yu J. Soft matter-water interface: general discussion. Faraday Discuss 2024; 249:485-520. [PMID: 38193511 DOI: 10.1039/d3fd90066e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
|
6
|
Hayton JA, Davies MB, Whale TF, Michaelides A, Cox SJ. The limit of macroscopic homogeneous ice nucleation at the nanoscale. Faraday Discuss 2024; 249:210-228. [PMID: 37791990 DOI: 10.1039/d3fd00099k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
Nucleation in small volumes of water has garnered renewed interest due to the relevance of pore condensation and freezing under conditions of low partial pressures of water, such as in the upper troposphere. Molecular simulations can in principle provide insight on this process at the molecular scale that is challenging to achieve experimentally. However, there are discrepancies in the literature as to whether the rate in confined systems is enhanced or suppressed relative to bulk water at the same temperature and pressure. In this study, we investigate the extent to which the size of the critical nucleus and the rate at which it grows in thin films of water are affected by the thickness of the film. Our results suggest that nucleation remains bulk-like in films that are barely large enough accommodate a critical nucleus. This conclusion seems robust to the presence of physical confining boundaries. We also discuss the difficulties in unambiguously determining homogeneous nucleation rates in nanoscale systems, owing to the challenges in defining the volume. Our results suggest any impact on a film's thickness on the rate is largely inconsequential for present day experiments.
Collapse
Affiliation(s)
- John A Hayton
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | - Michael B Davies
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK
| | - Thomas F Whale
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK
- School of Earth and Environment, University of Leeds, Leeds, UK
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| |
Collapse
|
7
|
Bowman GR, Cox SJ, Dellago C, DuBay KH, Eaves JD, Fletcher DA, Frechette LB, Grünwald M, Klymko K, Ku J, Omar A, Rabani E, Reichman DR, Rogers JR, Rosnik AM, Rotskoff GM, Schneider AR, Schwierz N, Sivak DA, Vaikuntanathan S, Whitelam S, Widmer-Cooper A. Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory. Annu Rev Phys Chem 2023; 74:1-27. [PMID: 36719975 DOI: 10.1146/annurev-physchem-101422-030127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective we celebrate his work at these frontiers. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 74 is April 2023. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Gregory R Bowman
- Bioengineering, Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | | | - Kateri H DuBay
- Department of Chemistry, University of Virginia, Charlottesville, Virginia, USA;
| | - Joel D Eaves
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado, USA
| | - Daniel A Fletcher
- Department of Bioengineering and Biophysics Program, University of California, Berkeley, California, USA.,Division of Biological Systems and Engineering, Lawrence Berkeley National Laboratory, Berkeley, California, USA.,Chan Zuckerberg Biohub, San Francisco, California, USA
| | - Layne B Frechette
- Martin A. Fisher School of Physics, Brandeis University, Waltham, Massachusetts, USA;
| | - Michael Grünwald
- Department of Chemistry, University of Utah, Salt Lake City, Utah, USA
| | - Katherine Klymko
- National Energy Research Scientific Computing Center, Lawrence Berkeley National Laboratory, Berkeley, California, USA
| | - JiYeon Ku
- R&D Center, Eloi Materials (EML) Co., Ltd, Suwon, Republic of Korea
| | - Ahmad Omar
- Department of Materials Science and Engineering, University of California, Berkeley, California, USA
| | - Eran Rabani
- Department of Chemistry, University of California, Berkeley, California, USA.,Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA.,The Raymond and Beverly Sackler Center of Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv, Israel
| | | | - Julia R Rogers
- Department of Systems Biology, Columbia University, New York, NY, USA;
| | | | - Grant M Rotskoff
- Department of Chemistry, Stanford University, Stanford, California, USA;
| | | | - Nadine Schwierz
- Institute of Physics, University of Augsburg, Augsburg, Germany
| | - David A Sivak
- Department of Physics, Simon Fraser University, Burnaby, British Columbia, Canada;
| | | | - Stephen Whitelam
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California, USA;
| | - Asaph Widmer-Cooper
- School of Chemistry, The University of Sydney, Sydney, New South Wales, Australia
| |
Collapse
|
8
|
Bui AT, Thiemann FL, Michaelides A, Cox SJ. Classical Quantum Friction at Water-Carbon Interfaces. Nano Lett 2023; 23:580-587. [PMID: 36626824 PMCID: PMC9881168 DOI: 10.1021/acs.nanolett.2c04187] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 01/04/2023] [Indexed: 05/20/2023]
Abstract
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.
Collapse
Affiliation(s)
- Anna T. Bui
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CambridgeCB2 1EW, United Kingdom
| | - Fabian L. Thiemann
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CambridgeCB2 1EW, United Kingdom
- Thomas
Young Centre, London Centre for Nanotechnology, and Department of
Physics and Astronomy, University College
London, Gower Street, LondonWC1E 6BT, United Kingdom
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, LondonSW7 2AZ, United Kingdom
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CambridgeCB2 1EW, United Kingdom
| | - Stephen J. Cox
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CambridgeCB2 1EW, United Kingdom
| |
Collapse
|
9
|
Abstract
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.
Collapse
Affiliation(s)
- Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
10
|
Anderson MW, Bennett M, Cedeno R, Cölfen H, Cox SJ, Cruz-Cabeza AJ, De Yoreo JJ, Drummond-Brydson R, Dudek MK, Fichthorn KA, Finney AR, Ford I, Galloway JM, Gebauer D, Grossier R, Harding JH, Hare A, Horváth D, Hunter L, Kim J, Kimura Y, Kirschhock CEA, Kiselev AA, Kras W, Kuttner C, Lee AY, Liao Z, Maini L, Nilsson Lill SO, Pellens N, Price SL, Rietveld IB, Rimer JD, Roberts KJ, Rogal J, Salvalaglio M, Sandei I, Schuszter G, Sefcik J, Sun W, Ter Horst JH, Ukrainczyk M, Van Driessche AES, Veesler S, Vekilov PG, Verma V, Whale T, Wheatcroft HP, Zeglinski J. Understanding crystal nucleation mechanisms: where do we stand? General discussion. Faraday Discuss 2022; 235:219-272. [PMID: 35789238 DOI: 10.1039/d2fd90021a] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
11
|
Abstract
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 Ih. 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.
Collapse
Affiliation(s)
- Dominic Atherton
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
12
|
Abstract
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 (∼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. Free energy calculations from molecular simulations reveal that water's interfacial dielectric response is well-described by bulk properties.![]()
Collapse
Affiliation(s)
- Stephen J. Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| | - Phillip L. Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
- Department of Chemistry, University of California, Berkeley, CA 94720, USA
| |
Collapse
|
13
|
Newman JL, Phillips JS, Cox SJ. Detecting positional vertigo using an ensemble of 2D convolutional neural networks. Biomed Signal Process Control 2021; 68:102708. [PMID: 34276807 PMCID: PMC8261823 DOI: 10.1016/j.bspc.2021.102708] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Revised: 04/04/2021] [Accepted: 04/28/2021] [Indexed: 11/18/2022]
Abstract
We trained Deep Neural Networks to detect attacks of motion provoked dizziness. 2D Convolutional Deep Neural Networks outperform 1D network architectures. Best results were provided by input features combining eye- and head-movement. An ensemble of five networks outperformed each individual network alone.
The aim of the work presented here was to develop a system that can automatically identify attacks of dizziness occurring in patients suffering from positional vertigo, which occurs when sufferers move their head into certain positions. We used our novel medical device, CAVA, to record eye- and head-movement data continually for up to 30 days in patients diagnosed with a disorder called Benign Paroxysmal Positional Vertigo. Building upon our previous work, we describe a novel ensemble of five 2D Convolutional Neural Networks, using composite recognition features, including eye-movement data and three-channel accelerometer data. We achieve an F1 score of 0.63 across an 11-fold cross-fold validation experiment, demonstrating that the system can detect a few seconds of motion provoked dizziness from within over a 100 h of normal eye-movement data. We show that the system outperforms our previous 1D Neural Network approach, and that our ensemble classifier is superior to each of the individual networks it contains. We also demonstrate that our composite recognition features provide improved performance over results obtained using the individual data sources independently.
Collapse
Affiliation(s)
- Jacob L. Newman
- The School of Computing Sciences, University of East Anglia, Norwich NR4 7TJ, United Kingdom
- Corresponding author.
| | - John S. Phillips
- The Department of Ear, Nose, and Throat Surgery, Norfolk & Norwich University Hospitals NHS Foundation Trust, Norwich NR4 7UY, United Kingdom
| | - Stephen J. Cox
- The School of Computing Sciences, University of East Anglia, Norwich NR4 7TJ, United Kingdom
| |
Collapse
|
14
|
Cox SJ, Mandadapu KK, Geissler PL. Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water. J Chem Phys 2021; 154:244502. [PMID: 34241373 DOI: 10.1063/5.0051399] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. We propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.
Collapse
Affiliation(s)
- Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Kranthi K Mandadapu
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Phillip L Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
15
|
Abstract
Vertigo is a type of dizziness characterised by the subjective feeling of movement despite being stationary. One in four individuals in the community experience symptoms of dizziness at any given time, and it can be challenging for clinicians to diagnose the underlying cause. When dizziness is the result of a malfunction in the inner-ear, the eyes flicker and this is called nystagmus. In this article we describe the first use of Deep Neural Network architectures applied to detecting nystagmus. The data used in these experiments was gathered during a clinical investigation of a novel medical device for recording head and eye movements. We describe methods for training networks using very limited amounts of training data, with an average of 11 mins of nystagmus across four subjects, and less than 24 hours of data in total, per subject. Our methods work by replicating and modifying existing samples to generate new data. In a cross-fold validation experiment, we achieve an average F1 score of 0.59 (SD = 0.24) across all four folds, showing that the methods employed are capable of identifying periods of nystagmus with a modest degree of accuracy. Notably, we were also able to identify periods of pathological nystagmus produced by a patient during an acute attack of Ménière's Disease, despite training the network on nystagmus that was induced by different means.
Collapse
|
16
|
Newman JL, Phillips JS, Cox SJ. Reconstructing animated eye movements from electrooculography data to aid the diagnosis of vestibular disorders. Int J Audiol 2021; 61:78-83. [PMID: 33825603 DOI: 10.1080/14992027.2021.1883196] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
OBJECTIVE To develop a method of visualising electrooculography data to improve the interpretability of nystagmus eye-movements captured using the Continuous Ambulatory Vestibular Assessment (CAVA®) device. DESIGN We are currently undertaking a clinical investigation to evaluate the capabilities of the CAVA® device to detect periods of pathological nystagmus. The work presented here was undertaken using unblinded data obtained from the preliminary phase of this investigation. STUDY SAMPLE One patient with Ménière's disease and one with Benign Paroxysmal Positional Vertigo. RESULTS Using the electrooculography data captured by the CAVA® device, we reconstructed 2D animations of patients' eye movements during attacks of vertigo. We were able to reanimate nystagmus produced as a consequence of two conditions. Concurrent video footage showed that the animations were visually very similar to the patient's actual eye-movements, excepting torsional eye-movements. CONCLUSIONS The reconstructed animations provide an alternative presentation modality, enabling clinicians to largely interpret electrooculography data as if they were present during a vertigo attack. We were able to recreate nystagmus from attacks experienced in the community rather than a clinical setting. This information provides an objective record of a patient's nystagmus and could be used to complement a full neurotologic history when considering diagnosis and treatment options.
Collapse
Affiliation(s)
| | - John S Phillips
- Norfolk & Norwich University Hospitals NHS Foundation Trust, Norwich, UK
| | | |
Collapse
|
17
|
Wang L, Chen J, Cox SJ, Liu L, Sosso GC, Li N, Gao P, Michaelides A, Wang E, Bai X. Microscopic Kinetics Pathway of Salt Crystallization in Graphene Nanocapillaries. Phys Rev Lett 2021; 126:136001. [PMID: 33861106 DOI: 10.1103/physrevlett.126.136001] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2020] [Revised: 02/08/2021] [Accepted: 02/26/2021] [Indexed: 06/12/2023]
Abstract
The fundamental understanding of crystallization, in terms of microscopic kinetic and thermodynamic details, remains a key challenge in the physical sciences. Here, by using in situ graphene liquid cell transmission electron microscopy, we reveal the atomistic mechanism of NaCl crystallization from solutions confined within graphene cells. We find that rock salt NaCl forms with a peculiar hexagonal morphology. We also see the emergence of a transitory graphitelike phase, which may act as an intermediate in a two-step pathway. With the aid of density functional theory calculations, we propose that these observations result from a delicate balance between the substrate-solute interaction and thermodynamics under confinement. Our results highlight the impact of confinement on both the kinetics and thermodynamics of crystallization, offering new insights into heterogeneous crystallization theory and a potential avenue for materials design.
Collapse
Affiliation(s)
- Lifen Wang
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Laboratory for Materials Science, Dongguan 523000, China
| | - Ji Chen
- School of Physics and the Collaborative Innovation Center of Quantum Matters, Peking University, Beijing 100871, China
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Lei Liu
- School of Materials Science and Engineering, Peking University, Beijing 100871, China
| | - Gabriele C Sosso
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Ning Li
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
- Electron Microscopy Laboratory, School of Physics, Peking University, Beijing 100871, China
| | - Peng Gao
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
- Electron Microscopy Laboratory, School of Physics, Peking University, Beijing 100871, China
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Department of Physics and Astronomy, and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
- London Centre for Nanotechnology, University College London, London WC1H 0AH, United Kingdom
| | - Enge Wang
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Laboratory for Materials Science, Dongguan 523000, China
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
- School of Physics, Liaoning University, Shenyang 110036, China
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing 100871, China
| | - Xuedong Bai
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Laboratory for Materials Science, Dongguan 523000, China
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
18
|
Phillips JS, Newman JL, Cox SJ. Clinical Techniques and Technology: Vestibular Telemetry. Otolaryngol Head Neck Surg 2021; 165:751-753. [PMID: 33650899 PMCID: PMC8564261 DOI: 10.1177/0194599821993411] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
When a patient presents to a clinician with dizziness, it can be difficult for the
patient to describe their symptoms in a clear manner, and clinical examination often
yields entirely normal results. Ideally, it would be favorable to measure key
physiological parameters during their episodes of dizziness. From a clinical perspective,
this would allow a more timely and more accurate diagnosis. From a research perspective,
it would allow a greater understanding of how the vestibular system malfunctions as a
consequence of vestibular disease. The authors of this report have been funded by the UK
Medical Research Council to develop and test a novel technology to measure, record, and
analyze key physiological parameters provided by the dizzy individual during an episode of
dizziness while active in the community. We provide the context to evolving work in this
field, the outcome of preliminary studies, and a consideration of future
opportunities.
Collapse
Affiliation(s)
- John S Phillips
- Norfolk & Norwich University Hospitals NHS Foundation Trust, Norwich, UK
| | | | | |
Collapse
|
19
|
Phillips JS, Newman JL, FitzGerald JE, Cox SJ. Implications of vestibular telemetry for the management of Ménière’s Disease—Our experience with three adults. Clin Otolaryngol 2021. [PMCID: PMC7983932 DOI: 10.1111/coa.13676] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- John S. Phillips
- Norfolk & Norwich University Hospitals NHS Foundation Trust Norwich UK
| | | | | | | |
Collapse
|
20
|
Abstract
BACKGROUND Foam sclerotherapy is the process of using an aqueous foam to deliver surfactant to a varicose vein to damage vein wall endothelial cells, causing the vein to spasm, collapse and ultimately be re-absorbed into the body. Aqueous foams are complex fluids that can exhibit a significant yield stress and high effective viscosity which depend on their composition, particularly the bubble size and liquid fraction. OBJECTIVE To characterise the properties of foams used for varicose vein sclerotherapy and determine their effectiveness in the displacement of blood during sclerotherapy. METHODS Foams are modelled as yield stress fluids and their flow profiles in a model vein are predicted. Values of the yield stress are determined from experimental data for three different foams using the Sauter mean of the bubble size distribution. Along with the measured liquid fraction of the foams, this information is collected into a Bingham number which entirely characterises the process of sclerotherapy. RESULTS Polydispersity in bubble size has a strong effect on the yield stress of a foam and the Sauter mean of the size distribution better captures the effects of a few large bubbles. Reducing the polydispersity increases the yield stress, and a higher yield stress results in a larger plug region moving along the vein, which is more effective in displacing blood. The width of the plug region is proportional to the Bingham number, which also has a quadratic dependence on the liquid fraction of the foam. Assuming typical values for the rate of injection of a foam, we predict that for a vein of diameter 5 mm, the most effective foams have low liquid fraction, a narrow size distribution, and a Bingham number B ≈ 4.5. CONCLUSIONS The Sauter mean radius provides the most appropriate measure of the bubble size for sclerotherapy and the Bingham number then provides a simple measure of the efficacy of foam sclerotherapy in a vein of a given size, and explains the ability of different foams to remove varicose veins. Foams containing small bubbles, with a narrow size distribution, and a low liquid fraction are beneficial for sclerotherapy.
Collapse
Affiliation(s)
- T G Roberts
- Department of Mathematics, Aberystwyth University, UK
| | - S J Cox
- Department of Mathematics, Aberystwyth University, UK
| | - A L Lewis
- Biocompatibles UK Ltd, a Boston Scientific Company, Lakeview, Camberley, Surrey, UK
| | - S A Jones
- Biocompatibles UK Ltd, a Boston Scientific Company, Lakeview, Camberley, Surrey, UK
| |
Collapse
|
21
|
Abstract
Attaining accurate average structural properties in a molecular simulation should be considered a prerequisite if one aims to elicit meaningful insights into a system's behavior. For charged surfaces in contact with an electrolyte solution, an obvious example is the density profile of ions along the direction normal to the surface. Here, we demonstrate that, in the slab geometry typically used in simulations, imposing an electric displacement field D determines the integrated surface charge density of adsorbed ions at charged interfaces. This allows us to obtain macroscopic surface charge densities irrespective of the slab thickness used in our simulations. We also show that the commonly used Yeh-Berkowitz method and the "mirrored slab" geometry both impose vanishing integrated surface charge densities. We present results both for relatively simple rocksalt (1 1 1) interfaces and the more complex case of kaolinite's basal faces in contact with an aqueous electrolyte solution.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stephen J Cox
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
22
|
Cox SJ, Thorpe DG, Shaffer PR, Geissler PL. Assessing long-range contributions to the charge asymmetry of ion adsorption at the air-water interface. Chem Sci 2020; 11:11791-11800. [PMID: 34094413 PMCID: PMC8162909 DOI: 10.1039/d0sc01947j] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Anions generally associate more favorably with the air–water interface than cations. In addition to solute size and polarizability, the intrinsic structure of the unperturbed interface has been discussed as an important contributor to this bias. Here we assess quantitatively the role that intrinsic charge asymmetry of water's surface plays in ion adsorption, using computer simulations to compare model solutes of various size and charge. In doing so, we also evaluate the degree to which linear response theory for solvent polarization is a reasonable approach for comparing the thermodynamics of bulk and interfacial ion solvation. Consistent with previous works on bulk ion solvation, we find that the average electrostatic potential at the center of a neutral, sub-nanometer solute at the air–water interface depends sensitively on its radius, and that this potential changes quite nonlinearly as the solute's charge is introduced. The nonlinear response closely resembles that of the bulk. As a result, the net nonlinearity of ion adsorption is weaker than in bulk, but still substantial, comparable to the apparent magnitude of macroscopically nonlocal contributions from the undisturbed interface. For the simple-point-charge model of water we study, these results argue distinctly against rationalizing ion adsorption in terms of surface potentials inherent to molecular structure of the liquid's boundary. Cations and anions have different affinities for the air-water interface. The intrinsic orientation of surface molecules suggests such an asymmetry, but the bias is dominated by solvent response that is spatially local and significantly nonlinear.![]()
Collapse
Affiliation(s)
- Stephen J Cox
- Department of Chemistry, University of Cambridge Lensfield Road Cambridge CB2 1EW UK
| | - Dayton G Thorpe
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA.,Department of Physics, University of California Berkeley CA 94720 USA
| | - Patrick R Shaffer
- Department of Chemistry, University of California Berkeley CA 94720 USA
| | - Phillip L Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA.,Department of Chemistry, University of California Berkeley CA 94720 USA
| |
Collapse
|
23
|
Jia M, Zhang C, Cox SJ, Sprik M, Cheng J. Computing Surface Acidity Constants of Proton Hopping Groups from Density Functional Theory-Based Molecular Dynamics: Application to the SnO 2(110)/H 2O Interface. J Chem Theory Comput 2020; 16:6520-6527. [PMID: 32794753 DOI: 10.1021/acs.jctc.0c00021] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Proton transfer at metal oxide/water interfaces plays an important role in electrochemistry, geochemistry, and environmental science. The key thermodynamic quantity to characterize this process is the surface acidity constant. An ab initio method that combines density functional theory-based molecular dynamics (DFTMD) and free energy perturbation theory has been established for computing surface acidity constants. However, it involves a reversible proton insertion procedure in which frequent proton hopping, e.g., for strong bases and some oxide surfaces (e.g., SnO2), can cause instability issues in electronic structure calculation. In the original implementation, harmonic restraining potentials are imposed on all O-H bonds (denoted by the VrH scheme) to prevent proton hopping and thus may not be applicable for systems involving spontaneous proton hopping. In this work, we introduce an improved restraining scheme with a repulsive potential Vrep to compute the surface acidities of systems in which proton hopping is spontaneous and fast. In this Vrep scheme, a Buckingham-type repulsive potential Vrep is applied between the deprotonation site and all other protons in DFTMD simulations. We first verify the Vrep scheme by calculating the pKa values of H2O and aqueous HS- solution (i.e., strong conjugate bases) and then apply it to the SnO2(110)/H2O interface. It is found that the Vrep scheme leads to a prediction of the point of zero charge (PZC) of 4.6, which agrees well with experiment. The intrinsic individual pKa values of the terminal five-coordinated Sn site (Sn5cOH2) and bridge oxygen site (Sn2ObrH+) are 4.4 and 4.7, respectively, both being almost the same as the PZC. The similarity of the two pKa values indicates that dissociation of terminal water has almost zero free energy at this proton hopping interface (i.e., partial water dissociation), as expected from the acid-base equilibrium on SnO2.
Collapse
Affiliation(s)
- Mei Jia
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, P. R. China
| | - Chao Zhang
- Department of Chemistry-Ångström Laboratory, Uppsala University, Lägerhyddsvägen 1, PO Box 538, Uppsala 75121, Sweden
| | - Stephen J Cox
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Michiel Sprik
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Jun Cheng
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, P. R. China
| |
Collapse
|
24
|
Abstract
OBJECTIVE We are currently undertaking a clinical investigation to evaluate the diagnostic capability of a system for detecting periods of pathological dizziness. This article presents an analysis of the data captured during an acute attack of Ménière's disease. DESIGN The Continuous Ambulatory Vestibular Assessment (CAVA) device is worn by patients in the community, and continuously records eye and head movement data (vestibular telemetry). STUDY SAMPLE A 53-year-old lady with a fifteen-year history of left-sided unilateral Ménière's disease. RESULTS The patient wore the device nearly continuously for thirty days. The data revealed a three-hour long attack of vertigo consisting of four separate phases of nystagmus. The duration, beat-direction and slow phase velocity of the nystagmus evolved through time. The first phase contained isolated nystagmus beats which preceded the patient's record of the vertigo attack onset but coincided with anticipation of an impending vertigo attack. CONCLUSIONS CAVA provides a unique insight into the physiological parameters present during episodes of dizziness. Here, it has provided the first full example of an acute Ménière's attack, including a period of prodrome. These findings have implications for the prediction of vertigo attack onset, for the diagnosis of Ménière's disease and other diseases resulting in dizziness.
Collapse
Affiliation(s)
- John S Phillips
- Norfolk & Norwich University Hospitals NHS Foundation Trust, Norwich, UK
| | | | | | - John FitzGerald
- Norfolk & Norwich University Hospitals NHS Foundation Trust, Norwich, UK
| |
Collapse
|
25
|
Newman JL, Phillips JS, Cox SJ, FitzGerald J, Bath A. Automatic nystagmus detection and quantification in long-term continuous eye-movement data. Comput Biol Med 2019; 114:103448. [DOI: 10.1016/j.compbiomed.2019.103448] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2019] [Revised: 09/13/2019] [Accepted: 09/13/2019] [Indexed: 01/18/2023]
|
26
|
Affiliation(s)
- Stephen J. Cox
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Michiel Sprik
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
27
|
Zhang T, Stackhouse PW, Cox SJ, Mikovitz JC, Long CN. Clear-Sky Shortwave Downward Flux at the Earth's Surface: Ground-Based Data vs. Satellite-Based Data. J Quant Spectrosc Radiat Transf 2019; 224:247-260. [PMID: 33505085 PMCID: PMC7837425 DOI: 10.1016/j.jqsrt.2018.11.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The radiative flux data and other meteorological data in the BSRN archive start from 1992, but the RadFlux data, the clear-sky radiative fluxes at the BSRN sites derived through regression analyses of actually observed clear-sky fluxes, did not come into existence until the early 2000s, and at first, they were limited to the 7 NOAA SURFRAD and 4 DOE ARM sites, a subset of the BSRN sites. Recently, the RadFlux algorithm was applied more extensively to the BSRN sites for the production of clear-sky ground-based fluxes. At the time of this writing, there are 7119 site-months of clear-sky fluxes at 42 BSRN sites spanning the time from 1992 to late 2017. These data provide an unprecedented opportunity to validate the satellite-based clear-sky fluxes. In this paper, the GEWEX SRB GSW(V3.0) shortwave downward fluxes spanning 24.5 years from 1983-07 to 2007-12, the CERES SYN1deg(Ed4A) and EBAF(Ed4.0) shortwave fluxes spanning 2000-03 to mid-2017 are compared with their RadFlux counterparts on the hourly, 3-hourly, daily and monthly time scales. All the three datasets show reasonable agreement with their ground-based counterparts. Comparison of the satellite-based surface shortwave clear-sky radiative fluxes to the BSRN RadFlux analysis shows negative biases. Further analysis shows that the satellite-based atmosphere contains greater aerosol optical paths as well as more precipitable water than RadFlux analysis estimates.
Collapse
Affiliation(s)
- Taiping Zhang
- SSAI/NASA Langley Research Center, One Enterprise Parkway, Suite 200, Hampton, VA 23666-5845
| | | | - Stephen J. Cox
- SSAI/NASA Langley Research Center, One Enterprise Parkway, Suite 200, Hampton, VA 23666-5845
| | - J. Colleen Mikovitz
- SSAI/NASA Langley Research Center, One Enterprise Parkway, Suite 200, Hampton, VA 23666-5845
| | - Charles N. Long
- NOAA ESRL Global Monitoring Division, 325 Broadway R/GMD, Boulder, CO 80305-3328
| |
Collapse
|
28
|
Sayer T, Cox SJ. Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation. Phys Chem Chem Phys 2019; 21:14546-14555. [DOI: 10.1039/c9cp02193k] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
AgI is a potent inorganic ice nucleating particle, a feature often attributed to the lattice match between its {0001} surfaces and ice. Dissolved ions are found to be essential to the stability of these polar surfaces, and crucial to ice formation.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry
- University of Cambridge
- Cambridge CB2 1EW
- UK
| | - Stephen J. Cox
- Department of Chemistry
- University of Cambridge
- Cambridge CB2 1EW
- UK
| |
Collapse
|
29
|
Zaccagnino F, Audebert A, Cox SJ. Simulation of surfactant transport during the rheological relaxation of two-dimensional dry foams. Phys Rev E 2018; 98:022801. [PMID: 30253619 DOI: 10.1103/physreve.98.022801] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Indexed: 11/07/2022]
Abstract
We describe a numerical model to predict the rheology of two-dimensional dry foams. The model accurately describes soap film curvature and viscous friction with the walls, and includes the transport of surfactant within the films and across the vertices where films meet. It accommodates the changes in foam topology that occur when a foam flows and, in particular, accurately represents the relaxation of the foam following a topological change. The model is validated against experimental data, allowing the prediction of elastic and viscous parameters associated with different surfactant solutions.
Collapse
Affiliation(s)
- F Zaccagnino
- Department of Mathematics, Aberystwyth University, Aberystwyth, Ceredigion, SY23 3BZ, United Kingdom
| | - A Audebert
- STLO, UMR1253, INRA, Agrocampus Ouest, 35000 Rennes, France
| | - S J Cox
- Department of Mathematics, Aberystwyth University, Aberystwyth, Ceredigion, SY23 3BZ, United Kingdom
| |
Collapse
|
30
|
Cox SJ, Kraynik AM, Weaire D, Hutzler S. Ideal wet two-dimensional foams and emulsions with finite contact angle. Soft Matter 2018; 14:5922-5929. [PMID: 29972190 DOI: 10.1039/c8sm00739j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We present simulations that show that the equilibrium structure of an ideal two-dimensional foam with a finite contact angle develops an inhomogeneity for high liquid fraction φ. In liquid-liquid emulsions this inhomogeneity is known as flocculation. In the case of an ordered foam this requires a perturbation, but in a disordered foam inhomogeneity grows steadily and spontaneously with φ, as demonstrated in our simulations performed with the Surface Evolver.
Collapse
Affiliation(s)
- S J Cox
- Department of Mathematics, Aberystwyth University, Aberystwyth SY23 3BZ, UK.
| | | | | | | |
Collapse
|
31
|
Affiliation(s)
- Stephen J. Cox
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Phillip L. Geissler
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
32
|
Cox SJ, Taylor DJF, Youngs TGA, Soper AK, Totton TS, Chapman RG, Arjmandi M, Hodges MG, Skipper NT, Michaelides A. Formation of Methane Hydrate in the Presence of Natural and Synthetic Nanoparticles. J Am Chem Soc 2018; 140:3277-3284. [PMID: 29401390 PMCID: PMC5860788 DOI: 10.1021/jacs.7b12050] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Natural gas hydrates occur widely
on the ocean-bed and in permafrost
regions, and have potential as an untapped energy resource. Their
formation and growth, however, poses major problems for the energy
sector due to their tendency to block oil and gas pipelines, whereas
their melting is viewed as a potential contributor to climate change.
Although recent advances have been made in understanding bulk methane
hydrate formation, the effect of impurity particles, which are always
present under conditions relevant to industry and the environment,
remains an open question. Here we present results from neutron scattering
experiments and molecular dynamics simulations that show that the
formation of methane hydrate is insensitive to the addition of a wide
range of impurity particles. Our analysis shows that this is due to
the different chemical natures of methane and water, with methane
generally excluded from the volume surrounding the nanoparticles.
This has important consequences for our understanding of the mechanism
of hydrate nucleation and the design of new inhibitor molecules.
Collapse
Affiliation(s)
- Stephen J Cox
- Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom.,Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom
| | - Diana J F Taylor
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom.,Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
| | - Tristan G A Youngs
- ISIS Facility , STFC Rutherford Appleton Laboratory , Harwell Oxford , Didcot OX11 0QX , United Kingdom
| | - Alan K Soper
- ISIS Facility , STFC Rutherford Appleton Laboratory , Harwell Oxford , Didcot OX11 0QX , United Kingdom
| | - Tim S Totton
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Richard G Chapman
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Mosayyeb Arjmandi
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Michael G Hodges
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Neal T Skipper
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom.,Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom.,Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
| |
Collapse
|
33
|
Abstract
We investigate the equilibrium properties of a single area-minimizing bubble trapped between two narrowly separated parallel curved plates. We begin with the case of a bubble trapped between concentric spherical plates. We develop a model which shows that the surface energy of the bubble is lower when confined between spherical plates than between flat plates. We confirm our findings by comparing against Surface Evolver simulations. We then derive a simple model for a bubble between arbitrarily curved parallel plates. The energy is found to be higher when the local Gaussian curvature of the plates is negative and lower when the curvature is positive. To check the validity of the model, we consider a bubble trapped between concentric tori. In the toroidal case, we find that the sensitivity of the bubble's energy to the local curvature acts as a geometric potential capable of driving bubbles from regions with negative to positive curvature.
Collapse
Affiliation(s)
- A Mughal
- Department of Mathematics, Aberystwyth University, Aberystwyth, Ceredigion SY23 3BZ, UK.,Theoretische Physik, Friedrich-Alexander-Universität Erlangen-Nürnberg - Staudtstr. 7, 91058 Erlangen, Germany
| | - S J Cox
- Department of Mathematics, Aberystwyth University, Aberystwyth, Ceredigion SY23 3BZ, UK
| | - G E Schröder-Turk
- Theoretische Physik, Friedrich-Alexander-Universität Erlangen-Nürnberg - Staudtstr. 7, 91058 Erlangen, Germany.,Murdoch University, School of Engineering and IT, Murdoch, Perth, WA6162, Australia
| |
Collapse
|
34
|
Abstract
We present a numerical study of quasiperiodic foams, in which the bubbles are generated as duals of quasiperiodic Frank-Kasper phases. These foams are investigated as potential candidates to the celebrated Kelvin problem for the partition of three-dimensional space with equal volume bubbles and minimal surface area. Interestingly, one of the computed structures falls close to (but still slightly above) the best known Weaire-Phelan periodic candidate. In addition we find a correlation between the normalized bubble surface area and the root mean squared deviation of the number of faces, giving an additional clue to understanding the main geometrical ingredients driving the Kelvin problem.
Collapse
Affiliation(s)
- S J Cox
- Department of Mathematics, Aberystwyth University, SY23 3BZ, United Kingdom
| | | | | | | |
Collapse
|
35
|
Kiselev A, Bachmann F, Pedevilla P, Cox SJ, Michaelides A, Gerthsen D, Leisner T. Active sites in heterogeneous ice nucleation—the example of K-rich feldspars. Science 2016; 355:367-371. [DOI: 10.1126/science.aai8034] [Citation(s) in RCA: 180] [Impact Index Per Article: 22.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2016] [Accepted: 11/23/2016] [Indexed: 11/02/2022]
Abstract
Ice formation on aerosol particles is a process of crucial importance to Earth’s climate and the environmental sciences, but it is not understood at the molecular level. This is partly because the nature of active sites, local surface features where ice growth commences, is still unclear. Here we report direct electron-microscopic observations of deposition growth of aligned ice crystals on feldspar, an atmospherically important component of mineral dust. Our molecular-scale computer simulations indicate that this alignment arises from the preferential nucleation of prismatic crystal planes of ice on high-energy (100) surface planes of feldspar. The microscopic patches of (100) surface, exposed at surface defects such as steps, cracks, and cavities, are thought to be responsible for the high ice nucleation efficacy of potassium (K)–feldspar particles.
Collapse
|
36
|
Zen A, Roch LM, Cox SJ, Hu XL, Sorella S, Alfè D, Michaelides A. Toward Accurate Adsorption Energetics on Clay Surfaces. J Phys Chem C Nanomater Interfaces 2016; 120:26402-26413. [PMID: 27917256 PMCID: PMC5126707 DOI: 10.1021/acs.jpcc.6b09559] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Revised: 10/31/2016] [Indexed: 05/31/2023]
Abstract
Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.
Collapse
Affiliation(s)
- Andrea Zen
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Loïc M. Roch
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Stephen J. Cox
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Xiao Liang Hu
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Sandro Sorella
- SISSA−International
School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- INFM
Democritos National Simulation Center, 34151 Trieste, Italy
| | - Dario Alfè
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United
Kingdom
| | - Angelos Michaelides
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
37
|
Affiliation(s)
| | - S J Cox
- St Thomas' Hospital, London SEI 7EH
| |
Collapse
|
38
|
Sosso G, Chen J, Cox SJ, Fitzner M, Pedevilla P, Zen A, Michaelides A. Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations. Chem Rev 2016; 116:7078-116. [PMID: 27228560 PMCID: PMC4919765 DOI: 10.1021/acs.chemrev.5b00744] [Citation(s) in RCA: 363] [Impact Index Per Article: 45.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2015] [Indexed: 11/28/2022]
Abstract
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 insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past 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 such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, 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 toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.
Collapse
Affiliation(s)
- Gabriele
C. Sosso
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Ji Chen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | | | - Martin Fitzner
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Philipp Pedevilla
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Andrea Zen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Angelos Michaelides
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| |
Collapse
|
39
|
Pedevilla P, Cox SJ, Slater B, Michaelides A. Can Ice-Like Structures Form on Non-Ice-Like Substrates? The Example of the K-feldspar Microcline. J Phys Chem C Nanomater Interfaces 2016; 120:6704-6713. [PMID: 27917255 PMCID: PMC5127609 DOI: 10.1021/acs.jpcc.6b01155] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 03/02/2016] [Indexed: 05/12/2023]
Abstract
Feldspar minerals are the most common rock formers in Earth's crust. As such they play an important role in subjects ranging from geology to climate science. An atomistic understanding of the feldspar structure and its interaction with water is therefore desirable, not least because feldspar has been shown to dominate ice nucleation by mineral dusts in Earth's atmosphere. The complexity of the ice/feldspar interface arising from the numerous chemical motifs expressed on the surface makes it a challenging system. Here we report a comprehensive study of this challenging system with ab initio density functional theory calculations. We show that the distribution of Al atoms, which is crucial for the dissolution kinetics of tectosilicate minerals, differs significantly between the bulk environment and on the surface. Furthermore, we demonstrate that water does not form ice-like overlayers in the contact layer on the most easily cleaved (001) surface of K-feldspar. We do, however, identify contact layer structures of water that induce ice-like ordering in the second overlayer. This suggests that even substrates without an apparent match with the ice structure may still act as excellent ice nucleating agents.
Collapse
Affiliation(s)
- Philipp Pedevilla
- Thomas
Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- London
Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Stephen J. Cox
- Thomas
Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- London
Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Ben Slater
- Thomas
Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- London
Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
- Thomas
Young Centre and Department of Physics and Astronomy, University College London, Gower Street, London WC1E
6BT, United Kingdom
- E-mail: . Phone: +44 207 679
0647
| |
Collapse
|
40
|
Fitzner M, Sosso GC, Cox SJ, Michaelides A. The Many Faces of Heterogeneous Ice Nucleation: Interplay Between Surface Morphology and Hydrophobicity. J Am Chem Soc 2015; 137:13658-69. [PMID: 26434775 DOI: 10.1021/jacs.5b08748] [Citation(s) in RCA: 105] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
What makes a material a good ice nucleating agent? Despite the importance of heterogeneous ice nucleation to a variety of fields, from cloud science to microbiology, major gaps in our understanding of this ubiquitous process still prevent us from answering this question. In this work, we have examined the ability of generic crystalline substrates to promote ice nucleation as a function of the hydrophobicity and the morphology of the surface. Nucleation rates have been obtained by brute-force molecular dynamics simulations of coarse-grained water on top of different surfaces of a model fcc crystal, varying the water-surface interaction and the surface lattice parameter. It turns out that the lattice mismatch of the surface with respect to ice, customarily regarded as the most important requirement for a good ice nucleating agent, is at most desirable but not a requirement. On the other hand, the balance between the morphology of the surface and its hydrophobicity can significantly alter the ice nucleation rate and can also lead to the formation of up to three different faces of ice on the same substrate. We have pinpointed three circumstances where heterogeneous ice nucleation can be promoted by the crystalline surface: (i) the formation of a water overlayer that acts as an in-plane template; (ii) the emergence of a contact layer buckled in an ice-like manner; and (iii) nucleation on compact surfaces with very high interaction strength. We hope that this extensive systematic study will foster future experimental work aimed at testing the physiochemical understanding presented herein.
Collapse
Affiliation(s)
- Martin Fitzner
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Gabriele C Sosso
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Stephen J Cox
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
41
|
Cox SJ, Kathmann SM, Slater B, Michaelides A. Molecular simulations of heterogeneous ice nucleation. II. Peeling back the layers. J Chem Phys 2015; 142:184705. [PMID: 25978903 DOI: 10.1063/1.4919715] [Citation(s) in RCA: 67] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Stephen J Cox
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Ben Slater
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
42
|
Cox SJ, Kathmann SM, Slater B, Michaelides A. Molecular simulations of heterogeneous ice nucleation. I. Controlling ice nucleation through surface hydrophilicity. J Chem Phys 2015; 142:184704. [PMID: 25978902 DOI: 10.1063/1.4919714] [Citation(s) in RCA: 104] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Stephen J Cox
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Ben Slater
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
43
|
Cox SJ, Towler MD, Alfè D, Michaelides A. Benchmarking the performance of density functional theory and point charge force fields in their description of sI methane hydrate against diffusion Monte Carlo. J Chem Phys 2015; 140:174703. [PMID: 24811651 DOI: 10.1063/1.4871873] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
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, while 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.
Collapse
Affiliation(s)
- Stephen J Cox
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Michael D Towler
- Department of Earth Sciences, University College London Gower Street, London WC1E 6BT, United Kingdom
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| |
Collapse
|
44
|
Cox SJ, Raza Z, Kathmann SM, Slater B, Michaelides A. The microscopic features of heterogeneous ice nucleation may affect the macroscopic morphology of atmospheric ice crystals. Faraday Discuss 2015; 167:389-403. [PMID: 24640502 DOI: 10.1039/c3fd00059a] [Citation(s) in RCA: 70] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
It is surprisingly difficult to freeze water. Almost all ice that forms under "mild" conditions (temperatures > -40 degrees C) 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.
Collapse
|
45
|
Kaiping GA, Jacobs GS, Cox SJ, Sluckin TJ. Nonequivalence of updating rules in evolutionary games under high mutation rates. Phys Rev E Stat Nonlin Soft Matter Phys 2014; 90:042726. [PMID: 25375542 DOI: 10.1103/physreve.90.042726] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/10/2014] [Indexed: 06/04/2023]
Abstract
Moran processes are often used to model selection in evolutionary simulations. The updating rule in Moran processes is a birth-death process, i. e., selection according to fitness of an individual to give birth, followed by the death of a random individual. For well-mixed populations with only two strategies this updating rule is known to be equivalent to selecting unfit individuals for death and then selecting randomly for procreation (biased death-birth process). It is, however, known that this equivalence does not hold when considering structured populations. Here we study whether changing the updating rule can also have an effect in well-mixed populations in the presence of more than two strategies and high mutation rates. We find, using three models from different areas of evolutionary simulation, that the choice of updating rule can change model results. We show, e. g., that going from the birth-death process to the death-birth process can change a public goods game with punishment from containing mostly defectors to having a majority of cooperative strategies. From the examples given we derive guidelines indicating when the choice of the updating rule can be expected to have an impact on the results of the model.
Collapse
Affiliation(s)
- G A Kaiping
- Computational Engineering and Design, University of Southampton, Southampton, United Kingdom
| | - G S Jacobs
- Applied Mathematics, University of Southampton, Southampton, United Kingdom
| | - S J Cox
- Computational Engineering and Design, University of Southampton, Southampton, United Kingdom
| | - T J Sluckin
- Applied Mathematics, University of Southampton, Southampton, United Kingdom
| |
Collapse
|
46
|
Ishihara S, Sugimura K, Cox SJ, Bonnet I, Bellaïche Y, Graner F. Comparative study of non-invasive force and stress inference methods in tissue. Eur Phys J E Soft Matter 2013; 36:9859. [PMID: 23615875 DOI: 10.1140/epje/i2013-13045-8] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2012] [Revised: 02/01/2013] [Accepted: 03/19/2013] [Indexed: 06/02/2023]
Abstract
In the course of animal development, the shape of tissue emerges in part from mechanical and biochemical interactions between cells. Measuring stress in tissue is essential for studying morphogenesis and its physical constraints. For that purpose, a possible new approach is force inference (up to a single prefactor) from cell shapes and connectivity. It is non-invasive and can provide space-time maps of stress in a whole tissue, unlike existing methods. To validate this approach, three force-inference methods, which differ in their approach of treating indefiniteness in an inverse problem between cell shapes and forces, were compared. Tests using two artificial and two experimental data sets consistently indicate that our Bayesian force inference, by which cell-junction tensions and cell pressures are simultaneously estimated, performs best in terms of accuracy and robustness. Moreover, by measuring the stress anisotropy and relaxation, we cross-validated the force inference and the global annular ablation of tissue, each of which relies on different prefactors. A practical choice of force-inference methods in different systems of interest is discussed.
Collapse
Affiliation(s)
- S Ishihara
- Graduate School of Arts and Sciences, The University of Tokyo, Komaba 3-8-1, Meguro-ku, Tokyo, Japan.
| | | | | | | | | | | |
Collapse
|
47
|
Holden JD, Cox SJ, Hargreaves SA. Revalidation: helping leaders to promote insight and avoid isolation. Educ Prim Care 2013; 24:19-21. [PMID: 23356758 DOI: 10.1080/14739879.2013.11493450] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
48
|
Cox SJ. Developing a "Research Test Bed" to introduce innovative Emission Testing Technology to improve New Zealand's Vehicle Emission Standards. ACTA ACUST UNITED AC 2012. [DOI: 10.1088/1742-6596/364/1/012106] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
49
|
|
50
|
Cox SJ, Kathmann SM, Purton JA, Gillan MJ, Michaelides A. Non-hexagonal ice at hexagonal surfaces: the role of lattice mismatch. Phys Chem Chem Phys 2012; 14:7944-9. [DOI: 10.1039/c2cp23438f] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|