1
|
Sundararaman R, Vigil-Fowler D, Schwarz K. Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface. Chem Rev 2022; 122:10651-10674. [PMID: 35522135 PMCID: PMC10127457 DOI: 10.1021/acs.chemrev.1c00800] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.
Collapse
Affiliation(s)
- Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th Street, Troy, New York 12180, United States
| | - Derek Vigil-Fowler
- Materials, Chemical, and Computational Science Directorate, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
2
|
Piaggi PM, Car R. Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1916634] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Pablo M. Piaggi
- Department of Chemistry, Princeton University, Princeton, NJ, USA
| | - Roberto Car
- Department of Chemistry and Department of Physics, Princeton University, Princeton, NJ, USA
| |
Collapse
|
3
|
Proton strings and rings in atypical nucleation of ferroelectricity in ice. Proc Natl Acad Sci U S A 2021; 118:2018837118. [PMID: 33443186 DOI: 10.1073/pnas.2018837118] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Ordinary ice has a proton-disordered phase which is kinetically metastable, unable to reach, spontaneously, the ferroelectric (FE) ground state at low temperature where a residual Pauling entropy persists. Upon light doping with KOH at low temperature, the transition to FE ice takes place, but its microscopic mechanism still needs clarification. We introduce a lattice model based on dipolar interactions plus a competing, frustrating term that enforces the ice rule (IR). In the absence of IR-breaking defects, standard Monte Carlo (MC) simulation leaves this ice model stuck in a state of disordered proton ring configurations with the correct Pauling entropy. A replica exchange accelerated MC sampling strategy succeeds, without open path moves, interfaces, or off-lattice configurations, in equilibrating this defect-free ice, reaching its low-temperature FE order through a well-defined first-order phase transition. When proton vacancies mimicking the KOH impurities are planted into the IR-conserving lattice, they enable standard MC simulation to work, revealing the kinetics of evolution of ice from proton disorder to partial FE order below the transition temperature. Replacing ordinary nucleation, each impurity opens up a proton ring generating a linear string, an actual FE hydrogen bond wire that expands with time. Reminiscent of those described for spin ice, these impurity-induced strings are proposed to exist in doped water ice too, where IRs are even stronger. The emerging mechanism yields a dependence of the long-time FE order fraction upon dopant concentration, and upon quenching temperature, that compares favorably with that known in real-life KOH doped ice.
Collapse
|
4
|
Xu J, Liu J, Liu J, Hu W, He X, Li J. Phase Transition of Ice at High Pressures and Low Temperatures. Molecules 2020; 25:molecules25030486. [PMID: 31979295 PMCID: PMC7037513 DOI: 10.3390/molecules25030486] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Revised: 01/08/2020] [Accepted: 01/21/2020] [Indexed: 11/21/2022] Open
Abstract
The behavior of ice under extreme conditions undergoes the change of intermolecular binding patterns and leads to the structural phase transitions, which are needed for modeling the convection and internal structure of the giant planets and moons of the solar system as well as H2O-rich exoplanets. Such extreme conditions limit the structural explorations in laboratory but open a door for the theoretical study. The ice phases IX and XIII are located in the high pressure and low temperature region of the phase diagram. However, to the best of our knowledge, the phase transition boundary between these two phases is still not clear. In this work, based on the second-order Møller–Plesset perturbation (MP2) theory, we theoretically investigate the ice phases IX and XIII and predict their structures, vibrational spectra and Gibbs free energies at various extreme conditions, and for the first time confirm that the phase transition from ice IX to XIII can occur around 0.30 GPa and 154 K. The proposed work, taking into account the many-body electrostatic effect and the dispersion interactions from the first principles, opens up the possibility of completing the ice phase diagram and provides an efficient method to explore new phases of molecular crystals.
Collapse
Affiliation(s)
- Jinjin Xu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Jinfeng Liu
- State Key Laboratory of Natural Medicines, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China;
| | - Jinyun Liu
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu 241000, China
- Correspondence: (J.L.); (X.H.); (J.L.)
| | - Wenxin Hu
- The Computer Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China;
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Correspondence: (J.L.); (X.H.); (J.L.)
| | - Jinjin Li
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu 241000, China
- Correspondence: (J.L.); (X.H.); (J.L.)
| |
Collapse
|
5
|
Shevkunov SV. Mechanism of Electric Polarization of Water Contact Layer at Its Interface with the Ion Crystal Surface. RUSS J ELECTROCHEM+ 2019. [DOI: 10.1134/s1023193519020083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
6
|
Goel H, Ling S, Ellis BN, Taconi A, Slater B, Rai N. Predicting vapor liquid equilibria using density functional theory: A case study of argon. J Chem Phys 2018; 148:224501. [PMID: 29907054 DOI: 10.1063/1.5025726] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase.
Collapse
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Sanliang Ling
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Breanna Nicole Ellis
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Anna Taconi
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| |
Collapse
|
7
|
Popov I, Lunev I, Khamzin A, Greenbaum Gutina A, Gusev Y, Feldman Y. The low-temperature dynamic crossover in the dielectric relaxation of ice I h. Phys Chem Chem Phys 2018; 19:28610-28620. [PMID: 29048435 DOI: 10.1039/c7cp05731h] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Based on the idea of defect migration as the principal mechanism in the dielectric relaxation of ice Ih, the concept of low-temperature dynamic crossover was proposed. It is known that at high temperatures, the diffusion of Bjerrum and ionic defects is high and their movement may be considered to be independent. Simple switching between these two mechanisms leads to a dynamic crossover at ∼235 K. By introducing coupling between the Bjerrum and ionic defects, it is possible to describe the smooth bend in the relaxation time at low temperatures in ice Ih. However, because the mobility of Bjerrum orientation defects slows down at low temperatures, they may create blockages for proton hopping. The trapping of ionic defects by L-D defects for a long period of time leads to an increase in the relaxation time and causes a low-temperature crossover. This model was validated by experimental dielectric measurements using various temperature protocols.
Collapse
Affiliation(s)
- Ivan Popov
- The Hebrew University of Jerusalem, Department of Applied Physics, Edmond J. Safra Campus, Givat Ram, Jerusalem, 9190401, Israel.
| | | | | | | | | | | |
Collapse
|
8
|
Liu Y, Ojamäe L. Raman and IR Spectra of Ice Ih and Ice XI with an Assessment of DFT Methods. J Phys Chem B 2016; 120:11043-11051. [PMID: 27690444 DOI: 10.1021/acs.jpcb.6b07001] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
IR and Raman spectroscopic technology can be directly used to identify the occurrence of ferroelectric ice XI in laboratory or extraterrestrial settings. The performance of 16 different DFT methods applied on the ice Ih, VIII, IX, and XI crystal phases is evaluated. Based on a selected DFT computational scheme, the IR and Raman spectra of ice Ih and XI are derived and compared. When the spectra, both IR and Raman, of ice Ih and ice XI are compared, the librational vibrations are found to be the most affected by the proton ordering. The spectroscopic fingerprint of ice XI can be used to distinguish ferroelectric ice XI from ice Ih in the universe. Furthermore, the existence of only one kind of H-bond in ice Ih is demonstrated from the overlapping subspectra for different types of H-bonded pair configurations in 16 isomers of ice Ih, which provides an illustration to the historic debate on whether one or two kinds of H-bonds existed in ice.
Collapse
Affiliation(s)
- Yuan Liu
- Department of Chemistry, IFM, Linköping University , SE-581 83 Linköping, Sweden
| | - Lars Ojamäe
- Department of Chemistry, IFM, Linköping University , SE-581 83 Linköping, Sweden
| |
Collapse
|
9
|
Pan D, Galli G. The fate of carbon dioxide in water-rich fluids under extreme conditions. SCIENCE ADVANCES 2016; 2:e1601278. [PMID: 27757424 PMCID: PMC5061492 DOI: 10.1126/sciadv.1601278] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/07/2016] [Accepted: 09/02/2016] [Indexed: 06/06/2023]
Abstract
Investigating the fate of dissolved carbon dioxide under extreme conditions is critical to understanding the deep carbon cycle in Earth, a process that ultimately influences global climate change. We used first-principles molecular dynamics simulations to study carbonates and carbon dioxide dissolved in water at pressures (P) and temperatures (T) approximating the conditions of Earth's upper mantle. Contrary to popular geochemical models assuming that molecular CO2(aq) is the major carbon species present in water under deep Earth conditions, we found that at 11 GPa and 1000 K, carbon exists almost entirely in the forms of solvated carbonate ([Formula: see text]) and bicarbonate ([Formula: see text]) ions and that even carbonic acid [H2CO3(aq)] is more abundant than CO2(aq). Furthermore, our simulations revealed that ion pairing between Na+ and [Formula: see text]/[Formula: see text] is greatly affected by P-T conditions, decreasing with increasing pressure at 800 to 1000 K. Our results suggest that in Earth's upper mantle, water-rich geofluids transport a majority of carbon in the form of rapidly interconverting [Formula: see text] and [Formula: see text] ions, not solvated CO2(aq) molecules.
Collapse
Affiliation(s)
- Ding Pan
- The Institute for Molecular Engineering, the University of Chicago, Chicago, IL 60637, USA
| | - Giulia Galli
- The Institute for Molecular Engineering, the University of Chicago, Chicago, IL 60637, USA
- Argonne National Laboratory, Argonne, IL 60439, USA
| |
Collapse
|
10
|
Zhang C, Hutter J, Sprik M. Computing the Kirkwood g-Factor by Combining Constant Maxwell Electric Field and Electric Displacement Simulations: Application to the Dielectric Constant of Liquid Water. J Phys Chem Lett 2016; 7:2696-2701. [PMID: 27352038 DOI: 10.1021/acs.jpclett.6b01127] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
In his classic 1939 paper, Kirkwood linked the macroscopic dielectric constant of polar liquids to the local orientational order as measured by the g-factor (later named after him) and suggested that the corresponding dielectric constant at short-range is effectively equal to the macroscopic value just after "a distance of molecular magnitude" [ Kirkwood, J. Chem. Phys., 1939, 7, 911 ]. Here, we show a simple approach to extract the short-ranged Kirkwood g-factor from molecular dynamics (MD) simulation by superposing the outcomes of constant electric field E and constant electric displacement D simulations [ Zhang and Sprik, Phys. Rev. B: Condens. Matter Mater. Phys., 2016, 93, 144201 ]. Rather than from the notoriously slow fluctuations of the dipole moment of the full MD cell, the dielectric constant can now be estimated from dipole fluctuations at short-range, accelerating the convergence. Exploiting this feature, we computed the bulk dielectric constant of liquid water modeled in the generalized gradient approximation (PBE) to density functional theory and found it to be at least 40% larger than the experimental value.
Collapse
Affiliation(s)
- Chao Zhang
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Jürg Hutter
- Institut für Chemie, Universität Zürich , Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Michiel Sprik
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
11
|
Gillan MJ, Alfè D, Michaelides A. Perspective: How good is DFT for water? J Chem Phys 2016; 144:130901. [DOI: 10.1063/1.4944633] [Citation(s) in RCA: 478] [Impact Index Per Article: 59.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Michael J. Gillan
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Dario Alfè
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
12
|
Engel EA, Monserrat B, Needs RJ. Vibrational renormalisation of the electronic band gap in hexagonal and cubic ice. J Chem Phys 2015; 143:244708. [DOI: 10.1063/1.4938029] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edgar A. Engel
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Bartomeu Monserrat
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
- Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019, USA
| | - Richard J. Needs
- TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
13
|
Del Ben M, Hutter J, VandeVondele J. Probing the structural and dynamical properties of liquid water with models including non-local electron correlation. J Chem Phys 2015; 143:054506. [PMID: 26254660 DOI: 10.1063/1.4927325] [Citation(s) in RCA: 75] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaric-isothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance.
Collapse
Affiliation(s)
- Mauro Del Ben
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Joost VandeVondele
- Department of Materials, ETH Zurich, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland
| |
Collapse
|
14
|
Shi L, Ni Y, Drews SEP, Skinner JL. Dielectric constant and low-frequency infrared spectra for liquid water and ice Ih within the E3B model. J Chem Phys 2015; 141:084508. [PMID: 25173022 DOI: 10.1063/1.4893792] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Two intrinsic difficulties in modeling condensed-phase water with conventional rigid non-polarizable water models are: reproducing the static dielectric constants for liquid water and ice Ih, and generating the peak at about 200 cm(-1) in the low-frequency infrared spectrum for liquid water. The primary physical reason for these failures is believed to be the missing polarization effect in these models, and consequently various sophisticated polarizable water models have been developed. However, in this work we pursue a different strategy and propose a simple empirical scheme to include the polarization effect only on the dipole surface (without modifying a model's intermolecular interaction potential). We implement this strategy for our explicit three-body (E3B) model. Our calculated static dielectric constants and low-frequency infrared spectra are in good agreement with experiment for both liquid water and ice Ih over wide temperature ranges, albeit with one fitting parameter for each phase. The success of our modeling also suggests that thermal fluctuations about local minima and the energy differences between different proton-disordered configurations play minor roles in the static dielectric constant of ice Ih. Our analysis shows that the polarization effect is important in resolving the two difficulties mentioned above and sheds some light on the origin of several features in the low-frequency infrared spectra for liquid water and ice Ih.
Collapse
Affiliation(s)
- L Shi
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | - Y Ni
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | - S E P Drews
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | - J L Skinner
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| |
Collapse
|
15
|
Gług M, Boczar M, Boda Ł, Wójcik MJ. Analysis of librational modes of ice XI studied by Car–Parrinello molecular dynamics. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.08.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
16
|
|
17
|
Abstract
We discuss the role of nuclear quantum fluctuations in ice Ih, focusing on the hydrogen-bond (HB) structure and the molecular dipole-moment distribution.
Collapse
Affiliation(s)
| | - Maurice de Koning
- Instituto de Física ‘Gleb Wataghin’
- Universidade Estadual de Campinas
- Campinas-SP
- Brazil
| |
Collapse
|
18
|
Del Ben M, VandeVondele J, Slater B. Periodic MP2, RPA, and Boundary Condition Assessment of Hydrogen Ordering in Ice XV. J Phys Chem Lett 2014; 5:4122-8. [PMID: 26278943 DOI: 10.1021/jz501985w] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Ice XV is the hydrogen-ordered form of the ice VI phase whose structure was predicted to be Cc and ferroelectric using periodic DFT approaches. However, neutron diffraction and Raman spectroscopy data show the structure to have P1̅ symmetry and to be antiferroelectric. Recent work1 using fragment-based MP2 and CCSD(T) approaches predicts the experimental structure as the ground state. We have reconsidered this problem using fully periodic MP2 and RPA approaches and find that the ferroelectric Cc structure is the lowest energy configuration. However, ubiquitously employed tinfoil boundary conditions stabilize polar structures. We suggest that ferroelectric Cc crystals can grow within a paraelectric ice VI matrix but may become unstable once a fraction of the matrix has become hydrogen-ordered. The reduction in dielectric constant causes P1̅ and other structures with small polarization to become favored, providing a possible resolution between observation and theoretical predictions.
Collapse
Affiliation(s)
- Mauro Del Ben
- †Department of Chemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Joost VandeVondele
- ‡Department of Materials, ETH Zürich, Wolfgang-Pauli-Strasse 27, CH-8093 Zürich, Switzerland
| | - Ben Slater
- §Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
19
|
Cassone G, Giaquinta PV, Saija F, Saitta AM. Effect of Electric Field Orientation on the Mechanical and Electrical Properties of Water Ices: An Ab-initio Study. J Phys Chem B 2014; 118:12717-24. [DOI: 10.1021/jp507376v] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Giuseppe Cassone
- Dipartimento
di Fisica e di Scienze della Terra, Università degli Studi di Messina, Contrada Papardo, 98166 Messina, Italy
- CNR-IPCF, Viale Ferdinando Stagno d’Alcontres
37, 98158 Messina, Italy
- UMR 7590, IMPMC,
UPMC Univ Paris 06, Sorbonne Universités, F-75005 Paris, France
- UMR 7590, IMPMC, CNRS, F-75005 Paris, France
| | - Paolo V. Giaquinta
- Dipartimento
di Fisica e di Scienze della Terra, Università degli Studi di Messina, Contrada Papardo, 98166 Messina, Italy
| | - Franz Saija
- CNR-IPCF, Viale Ferdinando Stagno d’Alcontres
37, 98158 Messina, Italy
| | - A. Marco Saitta
- UMR 7590, IMPMC,
UPMC Univ Paris 06, Sorbonne Universités, F-75005 Paris, France
- UMR 7590, IMPMC, CNRS, F-75005 Paris, France
| |
Collapse
|
20
|
Gonzalez MA, Sanz E, McBride C, Abascal JLF, Vega C, Valeriani C. Nucleation free-energy barriers with Hybrid Monte-Carlo/Umbrella Sampling. Phys Chem Chem Phys 2014; 16:24913-9. [DOI: 10.1039/c4cp02817a] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
21
|
Alfè D, Bartók AP, Csányi G, Gillan MJ. Analyzing the errors of DFT approximations for compressed water systems. J Chem Phys 2014; 141:014104. [DOI: 10.1063/1.4885440] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D. Alfè
- Department of Earth Sciences, UCL, London WC1E 6BT, United Kingdom
- London Centre for Nanotechnology, UCL, London WC1H 0AH, United Kingdom
- Thomas Young Centre, UCL, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, UCL, London WC1E 6BT, United Kingdom
| | - A. P. Bartók
- Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - G. Csányi
- Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - M. J. Gillan
- London Centre for Nanotechnology, UCL, London WC1H 0AH, United Kingdom
- Thomas Young Centre, UCL, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, UCL, London WC1E 6BT, United Kingdom
| |
Collapse
|
22
|
Pan D, Wan Q, Galli G. The refractive index and electronic gap of water and ice increase with increasing pressure. Nat Commun 2014; 5:3919. [PMID: 24861665 PMCID: PMC4050267 DOI: 10.1038/ncomms4919] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2014] [Accepted: 04/16/2014] [Indexed: 11/08/2022] Open
Abstract
Determining the electronic and dielectric properties of water at high pressure and temperature is an essential prerequisite to understand the physical and chemical properties of aqueous environments under supercritical conditions, for example, in the Earth interior. However, optical measurements of compressed ice and water remain challenging, and it has been common practice to assume that their band gap is inversely correlated with the measured refractive index, consistent with observations reported for hundreds of materials. Here we report ab initio molecular dynamics and electronic structure calculations showing that both the refractive index and the electronic gap of water and ice increase with increasing pressure, at least up to 30 GPa. Subtle electronic effects, related to the nature of interband transitions and band edge localization under pressure, are responsible for this apparently anomalous behaviour.
Collapse
Affiliation(s)
- Ding Pan
- Department of Chemistry, University of California, Davis, California 95616, USA
| | - Quan Wan
- Department of Chemistry, University of California, Davis, California 95616, USA
- The Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Giulia Galli
- The Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|