1
|
Dufils T, Schran C, Chen J, Geim AK, Fumagalli L, Michaelides A. Origin of dielectric polarization suppression in confined water from first principles. Chem Sci 2024; 15:516-527. [PMID: 38179530 PMCID: PMC10763014 DOI: 10.1039/d3sc04740g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Accepted: 11/23/2023] [Indexed: 01/06/2024] Open
Abstract
It has long been known that the dielectric constant of confined water should be different from that in bulk. Recent experiments have shown that it is vanishingly small, however the origin of the phenomenon remains unclear. Here we used ab initio molecular dynamics simulations (AIMD) and AIMD-trained machine-learning potentials to understand water's structure and electronic properties underpinning this effect. For the graphene and hexagonal boron-nitride substrates considered, we find that it originates in the spontaneous anti-parallel alignment of the water dipoles in the first two water layers near the solid interface. The interfacial layers exhibit net ferroelectric ordering, resulting in an overall anti-ferroelectric arrangement of confined water. Together with constrained hydrogen-bonding orientations, this leads to much reduced out-of-plane polarization. Furthermore, we directly contrast AIMD and simple classical force-field simulations, revealing important differences. This work offers insight into a property of water that is critical in modulating surface forces, the electric-double-layer formation and molecular solvation, and shows a way to compute it.
Collapse
Affiliation(s)
- T Dufils
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - C Schran
- Cavendish Laboratory, Department of Physics, University of Cambridge Cambridge CB3 0HE UK
- Lennard-Jones Centre, University of Cambridge Trinity Ln Cambridge CB2 1TN UK
| | - J Chen
- School of Physics, Peking University Beijing 100871 China
| | - A K Geim
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - L Fumagalli
- Department of Physics and Astronomy, University of Manchester Manchester M13 9PL UK
- National Graphene Institute, University of Manchester Manchester M13 9PL UK
| | - A Michaelides
- Lennard-Jones Centre, University of Cambridge Trinity Ln Cambridge CB2 1TN UK
- Yusuf Hamied Department of Chemistry, University of Cambridge Lensfield Road Cambridge CB2 1EW UK
| |
Collapse
|
2
|
Ryan MJ, Yang N, Kwac K, Wilhelm KB, Chi BK, Weix DJ, Cho M, Zanni MT. The hydrogen-bonding dynamics of water to a nitrile-functionalized electrode is modulated by voltage according to ultrafast 2D IR spectroscopy. Proc Natl Acad Sci U S A 2023; 120:e2314998120. [PMID: 38127983 PMCID: PMC10756189 DOI: 10.1073/pnas.2314998120] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Accepted: 11/15/2023] [Indexed: 12/23/2023] Open
Abstract
We report the hydrogen-bonding dynamics of water to a nitrile-functionalized and plasmonic electrode surface as a function of applied voltage. The surface-enhanced two-dimensional infrared spectra exhibit hydrogen-bonded and non-hydrogen-bonded nitrile features in similar proportions, plus cross peaks between the two. Isotopic dilution experiments show that the cross peaks arise predominantly from chemical exchange between hydrogen-bonded and non-hydrogen-bonded nitriles. The chemical exchange rate depends upon voltage, with the hydrogen bond of the water to the nitriles breaking 2 to 3 times slower (>63 vs. 25 ps) under a positive as compared to a negative potential. Spectral diffusion created by hydrogen-bond fluctuations occurs on a ~1 ps timescale and is moderately potential-dependent. Timescales from molecular dynamics simulations agree qualitatively with the experiment and show that a negative voltage causes a small net displacement of water away from the surface. These results show that the voltage applied to an electrode can alter the timescales of solvent motion at its interface, which has implications for electrochemically driven reactions.
Collapse
Affiliation(s)
- Matthew J. Ryan
- Department of Chemistry, University of Wisconsin-Madison, Madison, WI53706
| | - Nan Yang
- Department of Chemistry, University of Wisconsin-Madison, Madison, WI53706
| | - Kijeong Kwac
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul02841, Republic of Korea
- Department of Chemistry, Korea University, Seoul02841, Republic of Korea
| | - Kiera B. Wilhelm
- Department of Chemistry, University of Wisconsin-Madison, Madison, WI53706
| | - Benjamin K. Chi
- Department of Chemistry, University of Wisconsin-Madison, Madison, WI53706
| | - Daniel J. Weix
- Department of Chemistry, University of Wisconsin-Madison, Madison, WI53706
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul02841, Republic of Korea
- Department of Chemistry, Korea University, Seoul02841, Republic of Korea
| | - Martin T. Zanni
- Department of Chemistry, University of Wisconsin-Madison, Madison, WI53706
| |
Collapse
|
3
|
Foroutan M, Boudaghi A, Alibalazadeh M. Fullerenes containing water molecules: a study of reactive molecular dynamics simulations. Phys Chem Chem Phys 2023; 25:32493-32502. [PMID: 37997178 DOI: 10.1039/d3cp04420c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2023]
Abstract
A different technique was used to investigate fullerenes encapsulating a polar guest species. By reactive molecular dynamics simulations, three types of fullerenes were investigated on a gold surface: an empty C60, a single H2O molecule inside C60 (H2O@C60), and two water molecules inside C60 ((H2O)2@C60). Our findings revealed that despite the free movement of all fullerenes on gold surfaces, confined H2O molecules within the fullerenes result in a distinct pattern of motion in these systems. The (H2O)2@C60 complex had the highest displacement and average velocity, while C60 had the lowest displacement and average velocity. The symmetry of molecules and the polarity of water seem to be crucial in these cases. ReaxFF simulations showed that water molecules in an H2O molecule, H2O@C60, and (H2O)2@C60 have dipole moments of 1.76, 0.42, and 0.47 D, respectively. A combination of the non-polar C60 and polar water demonstrated a significant reduction in the dipole moment of H2O molecules due to encapsulation. The dipole moments of water molecules agreed with those in other studies, which can be useful in the development of biocompatible and high-efficiency nanocars.
Collapse
Affiliation(s)
- Masumeh Foroutan
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran, Iran.
| | - Ahmad Boudaghi
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran, Iran.
| | - Mahtab Alibalazadeh
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran, Iran.
| |
Collapse
|
4
|
Devlin SW, Jamnuch S, Xu Q, Chen AA, Qian J, Pascal TA, Saykally RJ. Agglomeration Drives the Reversed Fractionation of Aqueous Carbonate and Bicarbonate at the Air-Water Interface. J Am Chem Soc 2023; 145:22384-22393. [PMID: 37774115 DOI: 10.1021/jacs.3c05093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/01/2023]
Abstract
In the course of our investigations of the adsorption of ions to the air-water interface, we previously reported the surprising result that doubly charged carbonate anions exhibit a stronger surface affinity than singly charged bicarbonate anions. In contrast to monovalent, weakly hydrated anions, which generally show enhanced concentrations in the interfacial region, multivalent (and strongly hydrated) anions are expected to show a much weaker surface propensity. In the present work, we use resonantly enhanced deep-UV second-harmonic generation spectroscopy to measure the Gibbs free energy of adsorption of both carbonate (CO32-) and bicarbonate (HCO3-) anions to the air-water interface. Contrasting the predictions of classical electrostatic theory and in support of our previous findings from X-ray photoelectron spectroscopy, we find that carbonate anions do indeed exhibit much stronger surface affinity than do the bicarbonate anions. Extensive computer simulations reveal that strong ion pairing of CO32- with the Na+ countercation in the interfacial region results in the formation of near-neutral agglomerate clusters, consistent with a theory of interfacial ion adsorption based on hydration free energy and capillary waves. Simulated X-ray photoelectron spectra predict a 1 eV shift in the carbonate spectra compared to that of bicarbonate, further confirming our experiments. These findings not only advance our fundamental understanding of ion adsorption chemistry but also impact important practical processes such as ocean acidification, sea-spray aerosol chemistry, and mammalian respiration physiology.
Collapse
Affiliation(s)
- Shane W Devlin
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, California 94720, United States
| | - Sasawat Jamnuch
- ATLAS Materials Science Laboratory, Department of Nano Engineering and Chemical Engineering, University of California, San Diego, La Jolla, California 92023, United States
| | - Qiang Xu
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, California 94720, United States
| | - Amanda A Chen
- ATLAS Materials Science Laboratory, Department of Nano Engineering and Chemical Engineering, University of California, San Diego, La Jolla, California 92023, United States
| | - Jin Qian
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, California 94720, United States
| | - Tod A Pascal
- ATLAS Materials Science Laboratory, Department of Nano Engineering and Chemical Engineering, University of California, San Diego, La Jolla, California 92023, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92023, United States
- Sustainable Power and Energy Center, University of California San Diego, La Jolla, California 92023, United States
| | - Richard J Saykally
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, California 94720, United States
| |
Collapse
|
5
|
Dinpajooh M, Matyushov DV. Interface Dielectric Constant of Water at the Surface of a Spherical Solute. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
|
6
|
Davey SB, Cameron AP, Latham KG, Donne SW. Combined step potential electrochemical spectroscopy and electrochemical impedance spectroscopy analysis of the glassy carbon electrode in an aqueous electrolyte. Electrochim Acta 2021. [DOI: 10.1016/j.electacta.2021.139220] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
7
|
Affiliation(s)
- Franz M Geiger
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60660, United States
| |
Collapse
|
8
|
Abstract
It has long been anticipated that dielectric constants of polar liquids are reduced in the interfacial layer. Recent experiments and computer simulations support these expectations. A strong reduction of the dielectric constant is found in the direction perpendicular to a planar substrate, while the parallel response is bulk-like. This Perspective highlights recent theoretical calculations and simulations with an eye on relating them to properties observable in the laboratory. The average interface dielectric constant computed from simulations connects to thin films experiments, but this cannot be extended to screening of charges. In contrast to dielectric theories where a single dielectric constant gauges both the polarization energy and screening, these two signatures of dielectric polarization diverge on the molecular scale. The reduction of the dielectric constant of water in thin films is currently viewed as a combined effect of geometric confinement imposed by the substrate and the reconstruction of water hydrogen bonds in the surface layer.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- School of Molecular Sciences and Department of Physics, Arizona State University, PO Box 871504, Tempe, Arizona 85287-1504, United States
| |
Collapse
|
9
|
Hagopian A, Falcone A, Ben Yahia M, Filhol JS. Ab initiomodelling of interfacial electrochemical properties: beyond implicit solvation limitations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:304001. [PMID: 34108293 DOI: 10.1088/1361-648x/ac0207] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Accepted: 05/17/2021] [Indexed: 06/12/2023]
Abstract
First-principles calculations are an important tool to investigate the complex processes occurring at solid/liquid interfaces which are at the heart of modern technologies. Currently, capturing the whole electrochemical environment at an interface, including the applied potential and solvation, still remains challenging as it necessitates to couple different approaches whose interactions are not fully understood. In this work, a grand canonical density functional theory approach is coupled with solvation models to investigate the electrochemical interfaces under applied potential. We show that a parametrized polarizable continuum model (PCM) which represent solvation in a mean field approach by a continuous polarizable media, possesses catastrophic limitations for the modelling of ionic and charged interfaces. We reveal the origin of PCM instabilities under chemical or electrochemical strong oxidation to be the consequence of a phase transition in the surface Li electronic structure. Thus, PCM undergoes an unphysical response to this phase transition by penetrating within the atomic radius of surface Li atoms. To recover a physical response, an explicit first solvation shell has to be included in addition to the PCM in order to properly describe the electrochemistry of the interface. The Fukui functions show that the first solvation shell becomes involved in the redox process as solvent electron doublet is transferred to the acidic Li+. If another explicit solvent layer is added, the interface electrochemical properties become independent of the PCM parameters: in particular capacitance can then be computed from a parameter-free electrochemical approach. This is an important conclusion as the experimental electrochemical capacitance are not easily found and thus the parametrization of the PCM for electrochemical interface can be difficult. This approach can easily be applied to investigate electrochemical properties at the atomic scale and generalized to any electrochemical device for which interfaces play a crucial role.
Collapse
Affiliation(s)
- Arthur Hagopian
- ICGM, Université de Montpellier, CNRS, ENSCM, Montpellier, France
- RS2E French network on Electrochemical Energy Storage, FR5439, Amiens, France
| | - Aurélie Falcone
- ICGM, Université de Montpellier, CNRS, ENSCM, Montpellier, France
- RS2E French network on Electrochemical Energy Storage, FR5439, Amiens, France
| | - Mouna Ben Yahia
- ICGM, Université de Montpellier, CNRS, ENSCM, Montpellier, France
- RS2E French network on Electrochemical Energy Storage, FR5439, Amiens, France
| | - Jean-Sébastien Filhol
- ICGM, Université de Montpellier, CNRS, ENSCM, Montpellier, France
- RS2E French network on Electrochemical Energy Storage, FR5439, Amiens, France
| |
Collapse
|
10
|
Abstract
The dielectric response of a polar solvent to an ion is analyzed in terms of the bound charge, the net charge that accumulates near the ion as a consequence of the inhomogeneous polarization of the surrounding solvent. We demonstrate that the total bound charge arising in a full molecular treatment is identical to the total bound charge from standard continuum theory. In continuum theory, the bound charge resides in an infinitely thin layer, while in a molecular description the bound charge is spread over a region of finite width. Near simple atomic ions, the width of the bound charge distribution is roughly 1.3 nm. By simulating a sequence of ion charges from 0.1 to 2 e, where e is the magnitude of the electron charge, we analyze the applicability of linear response theory, which has been used by several authors. With increasing charge, the nonlinear response extends to an increasing distance from the ion. However, outside the region containing bound charge, the response is linear and in accord with continuum theory. Previous attempts to assign a dielectric constant to a solvent in the interfacial region are analyzed.
Collapse
Affiliation(s)
- Min-Sang Lee
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, United States
| | - Sherwin J Singer
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
11
|
Qi C, Zhu Z, Wang C, Zheng Y. Anomalously Low Dielectric Constant of Ordered Interfacial Water. J Phys Chem Lett 2021; 12:931-937. [PMID: 33439661 DOI: 10.1021/acs.jpclett.0c03299] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Despite considerable effort, the dielectric constant of interfacial water at solid surfaces is still not fully understood, thus hindering our understanding of the ubiquitous physical interactions in many materials and biological surfaces. In this study, we used molecular dynamics simulations to show that the parallel dielectric constant at the solid/water interface depends on solid-water interactions as well as the interfacial water structure on various solid crystal faces. In particular, ordered water structures can lead to a significant reduction (∼44%) in the parallel dielectric constant at the solid/water interface compared with that of bulk water. This sharp decrease in the parallel dielectric constant can be attributed to the specific antiferroelectric ordered structure of interfacial water molecules, which significantly suppresses the amplitude of the dipolar fluctuation associated with both the number of hydrogen bonds and the degree of order of interfacial water.
Collapse
Affiliation(s)
- Chonghai Qi
- School of Physics, Shandong University, Jinan 250100, China
- Key Laboratory of Interfacial Physics and Technology, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China
| | - Zhi Zhu
- School of Optical-Electrical Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
| | - Chunlei Wang
- Key Laboratory of Interfacial Physics and Technology, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China
- Zhangjiang Lab Interdisplinary Research Center, Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China
| | - Yujun Zheng
- School of Physics, Shandong University, Jinan 250100, China
| |
Collapse
|
12
|
Affiliation(s)
- Tuhin Samanta
- Department of Physics and School of Molecular Sciences, Arizona State University, P.O. Box 871604, Tempe, Arizona 85287-1604, USA
| | - Dmitry V. Matyushov
- Department of Physics and School of Molecular Sciences, Arizona State University, P.O. Box 871604, Tempe, Arizona 85287-1604, USA
| |
Collapse
|
13
|
Matyushov DV. Electrostatic solvation and mobility in uniform and non-uniform electric fields: From simple ions to proteins. BIOMICROFLUIDICS 2019; 13:064106. [PMID: 31737155 PMCID: PMC6837943 DOI: 10.1063/1.5124390] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Accepted: 10/23/2019] [Indexed: 06/10/2023]
Abstract
A number of observations related to interfacial electrostatics of polar liquids question the traditional assumption of dielectric theories that bulk dielectric properties can be continuously extended to the dividing surface separating the solute from the solvent. The deficiency of this approximation can be remedied by introducing local interface susceptibilities and the interface dielectric constant. Asymmetries of ionic hydration thermodynamics and of the mobility between cations and anions can be related to different propensities of the water molecules to orient their dipole toward and outward from solutes of opposite charges. This electrostatic asymmetry is reflected in different interface dielectric constants for cations and anions. The interface of water with neutral solutes is spontaneously polarized due to preferential water orientations in the interface. This phenomenon is responsible for a nonzero cavity potential directly related to a nonzero surface charge. This connection predicts that particles allowing a nonzero cavity potential must show mobility in an external electric field even if the net charge of the particle is zero. The theory predicts that a positive cavity potential and a positive surface charge translate to an effectively negative solute charge reported by mobility measurements. Passing of the cavity potential through a minimum found in simulations might be the origin of the maximum of mobility vs the ionic size observed experimentally. Finally, mobility of proteins in the field gradient (dielectrophoresis) is many orders of magnitude greater than predicted by the traditionally used Clausius-Mossotti equation. Two reasons contribute to this disagreement: (i) a failure of Maxwell's electrostatics to describe the cavity-field susceptibility and (ii) the neglect of the protein permanent dipole by the Clausius-Mossotti equation. An analytical relation between the dielectrophoretic susceptibility and dielectric spectroscopy of solutions provides direct access to this parameter, confirming the failure of the Clausius-Mossotti equation in application to protein dielectrophresis.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, Arizona 85287-1504, USA
| |
Collapse
|
14
|
Seyedi S, Martin DR, Matyushov DV. Screening of Coulomb interactions in liquid dielectrics. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:325101. [PMID: 31042681 DOI: 10.1088/1361-648x/ab1e6f] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The interaction of charges in dielectric materials is screened by the dielectric constant of the bulk dielectric. In dielectric theories, screening is assigned to surface charges appearing from preferential orientations of dipoles along the local field in the interface. For liquid dielectrics, such interfacial orientations are affected by the interfacial structure characterized by a separate interfacial dipolar susceptibility. We argue that dielectric properties of polar liquids should be characterized by two distinct susceptibilities responsible for local response (solvation) and long-range response (dielectric screening). A microscopic model for dipolar screening in polar liquids is developed here. It shows that the standard bulk dielectric constant is responsible for screening at large distances. The potential of mean force between ions in polar liquids becomes oscillatory at short distances. Oscillations arise from the coupling of collective longitudinal excitations (dipolarons) of the polar liquid with its interfacial density profile. Interfacial electrostatics demonstrates a crossover beyond the solute radius of ∼1 nm from the scaling laws roughly consistent with continuum expectations for small solutes to new scaling trends characterizing the much softer nano-scale solute-solvent interface. This crossover also marks the transition to a continuum-type electrostatic screening between ions, when short-distance oscillation of the potential of mean force become strongly damped. Screening oscillations are enhanced for more structured interfaces.
Collapse
Affiliation(s)
- Salman Seyedi
- Department of Physics, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, United States of America
| | | | | |
Collapse
|