1
|
Nakano K, Sorella S, Alfè D, Zen A. Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons. J Chem Theory Comput 2024. [PMID: 38788330 DOI: 10.1021/acs.jctc.4c00139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2024]
Abstract
Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.
Collapse
Affiliation(s)
- Kousuke Nakano
- Center for Basic Research on Materials, National Institute for Materials Science (NIMS), Tsukuba, Ibaraki 305-0047, Japan
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Sandro Sorella
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Dario Alfè
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, 80126 Napoli, Italy
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, U.K
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, U.K
| | - Andrea Zen
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, 80126 Napoli, Italy
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, U.K
| |
Collapse
|
2
|
Finney AR, Salvalaglio M. Properties of aqueous electrolyte solutions at carbon electrodes: effects of concentration and surface charge on solution structure, ion clustering and thermodynamics in the electric double layer. Faraday Discuss 2024; 249:334-362. [PMID: 37781909 DOI: 10.1039/d3fd00133d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/03/2023]
Abstract
Surfaces are able to control physical-chemical processes in multi-component solution systems and, as such, find application in a wide range of technological devices. Understanding the structure, dynamics and thermodynamics of non-ideal solutions at surfaces, however, is particularly challenging. Here, we use Constant Chemical Potential Molecular Dynamics (CμMD) simulations to gain insight into aqueous NaCl solutions in contact with graphite surfaces at high concentrations and under the effect of applied surface charges: conditions where mean-field theories describing interfaces cannot (typically) be reliably applied. We discover an asymmetric effect of surface charge on the electric double layer structure and resulting thermodynamic properties, which can be explained by considering the affinity of the surface for cations and anions and the cooperative adsorption of ions that occurs at higher concentrations. We characterise how the sign of the surface charge affects ion densities and water structure in the double layer and how the capacitance of the interface-a function of the electric potential drop across the double layer-is largely insensitive to the bulk solution concentration. Notably, we find that negatively charged graphite surfaces induce an increase in the size and concentration of extended liquid-like ion clusters confined to the double layer. Finally, we discuss how concentration and surface charge affect the activity coefficients of ions and water at the interface, demonstrating how electric fields in this region should be explicitly considered when characterising the thermodynamics of both solute and solvent at the solid/liquid interface.
Collapse
Affiliation(s)
- Aaron R Finney
- Thomas Young Centre and Department of Chemical Engineering, University College London, London WC1E 7JE, UK.
| | - Matteo Salvalaglio
- Thomas Young Centre and Department of Chemical Engineering, University College London, London WC1E 7JE, UK.
| |
Collapse
|
3
|
Hou R, Li C, Pan D. Raman and IR spectra of water under graphene nanoconfinement at ambient and extreme pressure-temperature conditions: a first-principles study. Faraday Discuss 2024; 249:181-194. [PMID: 37791622 DOI: 10.1039/d3fd00111c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The nanoconfinement of water can result in dramatic differences in its physical and chemical properties compared to bulk water. However, a detailed molecular-level understanding of these properties is still lacking. Vibrational spectroscopy, such as Raman and infrared, is a popular experimental tool for studying the structure and dynamics of water, and is often complemented by atomistic simulations to interpret experimental spectra, but there have been few theoretical spectroscopy studies of nanoconfined water using first-principles methods at ambient conditions, let alone under extreme pressure-temperature conditions. Here, we compute the Raman and IR spectra of water nanoconfined by graphene at ambient and extreme pressure-temperature conditions using ab initio simulations. Our results revealed alterations in the Raman stretching and low-frequency bands due to the graphene confinement. We also found spectroscopic evidence indicating that nanoconfinement considerably changes the tetrahedral hydrogen bond network, which is typically found in bulk water. Furthermore, we observed an unusual bending band in the Raman spectrum at ∼10 GPa and 1000 K, which is attributed to the unique molecular structure of confined ionic water. Additionally, we found that at ∼20 GPa and 1000 K, confined water transformed into a superionic fluid, making it challenging to identify the IR stretching band. Finally, we computed the ionic conductivity of confined water in the ionic and superionic phases. Our results highlight the efficacy of Raman and IR spectroscopy in studying the structure and dynamics of nanoconfined water in a large pressure-temperature range. Our predicted Raman and IR spectra can serve as a valuable guide for future experiments.
Collapse
Affiliation(s)
- Rui Hou
- Department of Physics, Hong Kong University of Science and Technology, Hong Kong, China.
- HKUST Shenzhen-Hong Kong Collaborative Innovation Research Institute, Shenzhen, China
| | - Chu Li
- Department of Physics, Hong Kong University of Science and Technology, Hong Kong, China.
- HKUST Shenzhen-Hong Kong Collaborative Innovation Research Institute, Shenzhen, China
| | - Ding Pan
- Department of Physics, Hong Kong University of Science and Technology, Hong Kong, China.
- Department of Chemistry, Hong Kong University of Science and Technology, Hong Kong, China
- HKUST Shenzhen-Hong Kong Collaborative Innovation Research Institute, Shenzhen, China
| |
Collapse
|
4
|
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
|
5
|
Feng Z, Lei Z, Yao Y, Liu J, Wu B, Ouyang W. Anisotropic Interfacial Force Field for Interfaces of Water with Hexagonal Boron Nitride. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:18198-18207. [PMID: 38063463 DOI: 10.1021/acs.langmuir.3c01612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
This study introduces an anisotropic interfacial potential that provides an accurate description of the van der Waals (vdW) interactions between water and hexagonal boron nitride (h-BN) at their interface. Benchmarked against the strongly constrained and appropriately normed functional, the developed force field demonstrates remarkable consistency with reference data sets, including binding energy curves and sliding potential energy surfaces for various configurations involving a water molecule adsorbed atop the h-BN surface. These findings highlight the significant improvement achieved by the developed force field in empirically describing the anisotropic vdW interactions of the water/h-BN heterointerfaces. Utilizing this anisotropic force field, molecular dynamics simulations demonstrate that atomically flat, pristine h-BN exhibits inherent hydrophobicity. However, when atomic-step surface roughness is introduced, the wettability of h-BN undergoes a significant change, leading to a hydrophilic nature. The calculated water contact angle (WCA) for the roughened h-BN surface is approximately 64°, which closely aligns with experimental WCA values ranging from 52° to 67°. These findings indicate the high probability of the presence of atomic steps on the surfaces of the experimental h-BN samples, emphasizing the need for further experimental verification. The development of the anisotropic interfacial force field for accurately describing interactions at the water/h-BN heterointerfaces is a significant advancement in accurately simulating the wettability of two-dimensional (2D) materials, offering a reliable tool for studying the dynamic and transport properties of water at these interfaces, with implications for materials science and nanotechnology.
Collapse
Affiliation(s)
- Zhicheng Feng
- Department of Engineering Mechanics, School of Civil Engineering, Wuhan University, Wuhan, Hubei 430072, China
| | - Zhangke Lei
- Department of Engineering Mechanics, School of Civil Engineering, Wuhan University, Wuhan, Hubei 430072, China
| | - Yuanpeng Yao
- Department of Engineering Mechanics, School of Civil Engineering, Wuhan University, Wuhan, Hubei 430072, China
| | - Jianxin Liu
- Department of Engineering Mechanics, School of Civil Engineering, Wuhan University, Wuhan, Hubei 430072, China
| | - Bozhao Wu
- Department of Engineering Mechanics, School of Civil Engineering, Wuhan University, Wuhan, Hubei 430072, China
| | - Wengen Ouyang
- Department of Engineering Mechanics, School of Civil Engineering, Wuhan University, Wuhan, Hubei 430072, China
- State Key Laboratory of Water Resources & Hydropower Engineering Science, Wuhan University, Wuhan, Hubei 430072, China
| |
Collapse
|
6
|
Nezval D, Bartošík M, Mach J, Švarc V, Konečný M, Piastek J, Špaček O, Šikola T. DFT study of water on graphene: Synergistic effect of multilayer p-doping. J Chem Phys 2023; 159:214710. [PMID: 38047516 DOI: 10.1063/5.0161160] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Accepted: 11/07/2023] [Indexed: 12/05/2023] Open
Abstract
Recent experiments related to a study concerning the adsorption of water on graphene have demonstrated the p-doping of graphene, although most of the ab initio calculations predict nearly zero doping. To shed more light on this problem, we have carried out van der Waals density functional theory calculations of water on graphene for both individual water molecules and continuous water layers with coverage ranging from one to eight monolayers. Furthermore, we have paid attention to the influence of the water molecule orientation toward graphene on its doping properties. In this article, we present the results of the band structure and the Bader charge analysis, showing the p-doping of graphene can be synergistically enhanced by putting 4-8 layers of an ice-like water structure on graphene having the water molecules oriented with oxygen atoms toward graphene.
Collapse
Affiliation(s)
- D Nezval
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
| | - M Bartošík
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
- Department of Physics and Materials Engineering, Faculty of Technology, Tomas Bata University in Zlín, Vavrečkova 275, 760 01 Zlín, Czech Republic
| | - J Mach
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
| | - V Švarc
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
| | - M Konečný
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
| | - J Piastek
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
| | - O Špaček
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
| | - T Šikola
- Institute of Physical Engineering, Brno University of Technology, Technická 2, 616 69 Brno, Czech Republic
- Central European Institute of Technology, Brno University of Technology, Purkyňova 656/123, 612 00 Brno, Czech Republic
| |
Collapse
|
7
|
Shi B, Zen A, Kapil V, Nagy PR, Grüneis A, Michaelides A. Many-Body Methods for Surface Chemistry Come of Age: Achieving Consensus with Experiments. J Am Chem Soc 2023; 145:25372-25381. [PMID: 37948071 PMCID: PMC10683001 DOI: 10.1021/jacs.3c09616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Revised: 10/15/2023] [Accepted: 10/17/2023] [Indexed: 11/12/2023]
Abstract
The adsorption energy of a molecule onto the surface of a material underpins a wide array of applications, spanning heterogeneous catalysis, gas storage, and many more. It is the key quantity where experimental measurements and theoretical calculations meet, with agreement being necessary for reliable predictions of chemical reaction rates and mechanisms. The prototypical molecule-surface system is CO adsorbed on MgO, but despite intense scrutiny from theory and experiment, there is still no consensus on its adsorption energy. In particular, the large cost of accurate many-body methods makes reaching converged theoretical estimates difficult, generating a wide range of values. In this work, we address this challenge, leveraging the latest advances in diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] to obtain accurate predictions for CO on MgO. These reliable theoretical estimates allow us to evaluate the inconsistencies in published temperature-programed desorption experiments, revealing that they arise from variations in employed pre-exponential factors. Utilizing this insight, we derive new experimental estimates of the (electronic) adsorption energy with a (more) precise pre-exponential factor. As a culmination of all of this effort, we are able to reach a consensus between multiple theoretical calculations and multiple experiments for the first time. In addition, we show that our recently developed cluster-based CCSD(T) approach provides a low-cost route toward achieving accurate adsorption energies. This sets the stage for affordable and reliable theoretical predictions of chemical reactions on surfaces to guide the realization of new catalysts and gas storage materials.
Collapse
Affiliation(s)
- Benjamin
X. Shi
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Andrea Zen
- Dipartimento
di Fisica Ettore Pancini, Università
di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
- Department
of Earth Sciences, University College London, Gower Street, WC1E 6BT London, U.K.
| | - Venkat Kapil
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Andreas Grüneis
- Institute
for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| |
Collapse
|
8
|
Neufeld VA, Berkelbach TC. Highly Accurate Electronic Structure of Metallic Solids from Coupled-Cluster Theory with Nonperturbative Triple Excitations. PHYSICAL REVIEW LETTERS 2023; 131:186402. [PMID: 37977636 DOI: 10.1103/physrevlett.131.186402] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Accepted: 09/28/2023] [Indexed: 11/19/2023]
Abstract
Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T))-often considered the "gold standard" of main-group quantum chemistry-is inapplicable to three-dimensional metals due to an infrared divergence, preventing its application to many important problems in materials science. We study the full, nonperturbative inclusion of triple excitations (CCSDT) and propose a new, iterative method, which we call ring-CCSDT, that resums the essential triple excitations with the same N^{7} run-time scaling as CCSD(T). CCSDT and ring-CCSDT are used to calculate the correlation energy of the uniform electron gas at metallic densities and the structural properties of solid lithium. Inclusion of connected triple excitations is shown to be essential to achieving high accuracy. We also investigate semiempirical CC methods based on spin-component scaling and the distinguishable cluster approximation and find that they enhance the accuracy of their parent ab initio methods.
Collapse
Affiliation(s)
- Verena A Neufeld
- 1Department of Chemistry, Columbia University, New York, New York 10027, USA
- 2Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | |
Collapse
|
9
|
Qian C, Zhou K. Ab Initio Molecular Dynamics Investigation of the Solvation States of Hydrated Ions in Confined Water. Inorg Chem 2023; 62:17756-17765. [PMID: 37855150 DOI: 10.1021/acs.inorgchem.3c02443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
Ionic transport in nanoscale channels with a critical size comparable to that of ions and solutes exhibits exceptional performance in water desalination, ion separation, electrocatalysts, and supercapacitors. However, the solvation states (SSs), i.e., the hydration structures and probability distribution, of hydrated ions in nanochannels differ from those in the bulk and the perspective of continuum theory. In this work, we conduct ab initio enhanced-sampling atomistic simulations to investigate the ion-specific SSs of monovalent ions (including Li+, Na+, K+, F-, Cl-, and I-) in the graphene channel with a width of 1 nm. Our findings highlight that the SSs of those ions are primarily determined by ion-water hydration, where ion-wall interactions play a minor role. The distribution of ions in layered confined water is a result of ion-specific hydration, which arises from the synergy of entropy and enthalpy. The free energy barriers for transitions between SSs are on the order of 1kBT, allowing for modulation through applying external fields or modifying surface properties. As the ion-wall interaction strengthens, as observed in vermiculite and carbides and nitrides of transition metal channels, the probability of near-wall SSs increases. These results help to improve the performance of nanofluidic devices and provide crucial insights for developing accurate force fields of molecular simulations or advanced theoretical approaches for ion dynamics in confined channels.
Collapse
Affiliation(s)
- Chen Qian
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
- Department of Chemistry, City University of Hong Kong, 83 Tat Chee Avenue, Kowloon 999077, Hong Kong, China
| | - Ke Zhou
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| |
Collapse
|
10
|
Riemelmoser S, Verdi C, Kaltak M, Kresse G. Machine Learning Density Functionals from the Random-Phase Approximation. J Chem Theory Comput 2023; 19:7287-7299. [PMID: 37800677 PMCID: PMC10601474 DOI: 10.1021/acs.jctc.3c00848] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Indexed: 10/07/2023]
Abstract
Kohn-Sham density functional theory (DFT) is the standard method for first-principles calculations in computational chemistry and materials science. More accurate theories such as the random-phase approximation (RPA) are limited in application due to their large computational cost. Here, we use machine learning to map the RPA to a pure Kohn-Sham density functional. The machine learned RPA model (ML-RPA) is a nonlocal extension of the standard gradient approximation. The density descriptors used as ingredients for the enhancement factor are nonlocal counterparts of the local density and its gradient. Rather than fitting only RPA exchange-correlation energies, we also include derivative information in the form of RPA optimized effective potentials. We train a single ML-RPA functional for diamond, its surfaces, and liquid water. The accuracy of ML-RPA for the formation energies of 28 diamond surfaces reaches that of state-of-the-art van der Waals functionals. For liquid water, however, ML-RPA cannot yet improve upon the standard gradient approximation. Overall, our work demonstrates how machine learning can extend the applicability of the RPA to larger system sizes, time scales, and chemical spaces.
Collapse
Affiliation(s)
- Stefan Riemelmoser
- Faculty
of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, A-1090 Vienna, Austria
| | - Carla Verdi
- Faculty
of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
- School
of Physics, The University of Sydney, Sydney, New South Wales 2006, Australia
- School
of Mathematics and Physics, The University
of Queensland, Brisbane, Queensland 4072, Australia
| | - Merzuk Kaltak
- VASP
Software GmbH, Sensengasse
8/12, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty
of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
- VASP
Software GmbH, Sensengasse
8/12, A-1090 Vienna, Austria
| |
Collapse
|
11
|
Ehlert C, Piras A, Gryn’ova G. CO 2 on Graphene: Benchmarking Computational Approaches to Noncovalent Interactions. ACS OMEGA 2023; 8:35768-35778. [PMID: 37810719 PMCID: PMC10551916 DOI: 10.1021/acsomega.3c03251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Accepted: 09/05/2023] [Indexed: 10/10/2023]
Abstract
Designing and optimizing graphene-based gas sensors in silico entail constructing appropriate atomistic representations for the physisorption complex of an analyte on an infinite graphene sheet, then selecting accurate yet affordable methods for geometry optimizations and energy computations. In this work, diverse density functionals (DFs), coupled cluster theory, and symmetry-adapted perturbation theory (SAPT) in conjunction with a range of finite and periodic surface models of bare and supported graphene were tested for their ability to reproduce the experimental adsorption energies of CO2 on graphene in a low-coverage regime. Periodic results are accurately reproduced by the interaction energies extrapolated from finite clusters to infinity. This simple yet powerful scheme effectively removes size dependence from the data obtained using finite models, and the latter can be treated at more sophisticated levels of theory relative to periodic systems. While for small models inexpensive DFs such as PBE-D3 afford surprisingly good agreement with the gold standard of quantum chemistry, CCSD(T), interaction energies closest to experiment are obtained by extrapolating the SAPT results and with nonlocal van der Waals functionals in the periodic setting. Finally, none of the methods and models reproduce the experimentally observed CO2 tilted adsorption geometry on the Pt(111) support, calling for either even more elaborate theoretical approaches or a revision of the experiment.
Collapse
Affiliation(s)
- Christopher Ehlert
- Heidelberg
Institute for Theoretical Studies (HITS gGmbH), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Interdisciplinary
Center for Scientific Computing (IWR), Heidelberg
University, Im Neuenheimer
Feld 368, 69120 Heidelberg, Germany
| | - Anna Piras
- Heidelberg
Institute for Theoretical Studies (HITS gGmbH), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Interdisciplinary
Center for Scientific Computing (IWR), Heidelberg
University, Im Neuenheimer
Feld 368, 69120 Heidelberg, Germany
| | - Ganna Gryn’ova
- Heidelberg
Institute for Theoretical Studies (HITS gGmbH), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Interdisciplinary
Center for Scientific Computing (IWR), Heidelberg
University, Im Neuenheimer
Feld 368, 69120 Heidelberg, Germany
| |
Collapse
|
12
|
Freyre P, St Pierre E, Rybolt T. Carbon Dioxide Capture by Adsorption in a Model Hydroxy-Modified Graphene Pore. Int J Mol Sci 2023; 24:11452. [PMID: 37511209 PMCID: PMC10380235 DOI: 10.3390/ijms241411452] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Revised: 07/09/2023] [Accepted: 07/10/2023] [Indexed: 07/30/2023] Open
Abstract
Concerns regarding the environmental impact of increasing levels of anthropogenic carbon dioxide have led to a variety of studies examining solid surfaces for their ability to trap this greenhouse gas (GHG). Atmospheric or post-combustion carbon capture requires an efficient separation of carbon dioxide and nitrogen gas. We used the molecular mechanics MM3 parameter set (previously shown to provide good estimates of molecule-surface binding energies) to calculate theoretical surface binding energies for carbon dioxide ∆E(CO2) and nitrogen ∆E(N2). For efficient separation, differentiation of these two gas-surface adsorption energies is required. Examined structures based on graphene, carbon slit width pore, and carbon nanotube gave ∆E(CO2) to ∆E(N2) ratios of 1.7, 1.8, and 1.9, respectively. To enhance the CO2 adsorption, we developed a model graphene surface pore lined with four hydroxy groups whose orientation allowed them to form hydrogen bonds with the oxygens in CO2. Both the single-layer and double-layer versions of this pore gave significant enhancement in the ability to trap CO2 preferentially to N2. The two-layer version of this pore gave ∆E(CO2) = 73 and ∆E(N2) = 6.8 kJ/mol. The one- and two-layer versions of this novel pore averaged a ∆E(CO2) to ∆E(N2) ratio of 12.
Collapse
Affiliation(s)
- Paige Freyre
- Department of Chemistry and Physics, University of Tennessee at Chattanooga, Chattanooga, TN 37403, USA
| | - Emalee St Pierre
- Department of Chemistry and Physics, University of Tennessee at Chattanooga, Chattanooga, TN 37403, USA
| | - Thomas Rybolt
- Department of Chemistry and Physics, University of Tennessee at Chattanooga, Chattanooga, TN 37403, USA
| |
Collapse
|
13
|
Tyagi R, Zen A, Voora VK. Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules. J Phys Chem A 2023. [PMID: 37406194 DOI: 10.1021/acs.jpca.3c02291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2023]
Abstract
Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C6X6)n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.
Collapse
Affiliation(s)
- Ritaj Tyagi
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Andrea Zen
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| |
Collapse
|
14
|
Tonel MZ, Abal JPK, Fagan SB, Barbosa MC. Ab initio study of water anchored in graphene pristine and vacancy-type defects. J Mol Model 2023; 29:198. [PMID: 37268861 DOI: 10.1007/s00894-023-05611-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 05/30/2023] [Indexed: 06/04/2023]
Abstract
CONTEXT In this paper, we have addressed two issues that are relevant to the interaction of water in pristine and vacant graphene through first-principles calculations based on the Density Functional Theory (DFT). The results showed that for the interaction of pristine graphene with water, the DOWN configuration (with the hydrogen atoms facing downwards) was the most stable, presenting binding energies in the order of -13.62 kJ/mol at a distance of 2.375 Å in the TOP position. We also evaluated the interaction of water with two vacancy models, removing one carbon atom (Vac-1C) and four atoms (Vac-4C). In the Vac-1C system, the most favourable system was the DOWN configuration, with binding energies ranging from -20.60 kJ/mol to -18.41 kJ/mol in the TOP and UP positions, respectively. A different behaviour was observed for the interaction of water with Vac-4C; regardless of the configuration of the water, it is always more favourable for the interaction to occur through the vacancy centre, with binding energies between -13.28 kJ/mol and -20.49 kJ/mol. Thus, the results presented open perspectives for the technological development of nanomembranes as well as providing a better understanding of the wettability effects of graphene sheets, whether pristine or with defects. METHOD We evaluated the interaction of pristine and vacant graphene with the water molecule, through calculations based on Density Functional Theory (DFT); implemented by the SIESTA program. The electronic, energetic, and structural properties were analyzed by solving self-consistent Kohn-Sham equations. In all calculations, a double ζ plus a polarized function (DZP) was used for the numerical baise set. Local Density Approximation (LDA) with the Perdew and Zunger (PZ) parameterisation along with a basis set superposition error (BSSE) correction were used to describe the exchange and correlation potential (Vxc). The water and isolated graphene structures were relaxed until the residual forces were less than 0.05 eV/Å-1 in all atomic coordinates.
Collapse
Affiliation(s)
- Mariana Zancan Tonel
- Universidade Franciscana-UFN, PPGNANO - Postgraduate Program in Nanoscience, Rua dos Andradas, 1614, ZIP, Santa Maria, RS, 97010-032, Brazil.
| | - João Pedro Kleinubing Abal
- Universidade Federal do Rio Grande do Sul- UFRGS, Institute of Physics, Av. Bento Gonçalves, 9500 - Agronomia, ZIP, Porto Alegre, RS, 91501-970, Brazil
| | - Solange Binotto Fagan
- Universidade Franciscana-UFN, PPGNANO - Postgraduate Program in Nanoscience, Rua dos Andradas, 1614, ZIP, Santa Maria, RS, 97010-032, Brazil
| | - Marcia Cristina Barbosa
- Universidade Federal do Rio Grande do Sul- UFRGS, Institute of Physics, Av. Bento Gonçalves, 9500 - Agronomia, ZIP, Porto Alegre, RS, 91501-970, Brazil
| |
Collapse
|
15
|
Lee S, Lee J, Zhai H, Tong Y, Dalzell AM, Kumar A, Helms P, Gray J, Cui ZH, Liu W, Kastoryano M, Babbush R, Preskill J, Reichman DR, Campbell ET, Valeev EF, Lin L, Chan GKL. Evaluating the evidence for exponential quantum advantage in ground-state quantum chemistry. Nat Commun 2023; 14:1952. [PMID: 37029105 PMCID: PMC10082187 DOI: 10.1038/s41467-023-37587-6] [Citation(s) in RCA: 16] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Accepted: 03/22/2023] [Indexed: 04/09/2023] Open
Abstract
Due to intense interest in the potential applications of quantum computing, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation, for generic chemical problems where heuristic quantum state preparation might be assumed to be efficient. The availability of exponential quantum advantage then centers on whether features of the physical problem that enable efficient heuristic quantum state preparation also enable efficient solution by classical heuristics. Through numerical studies of quantum state preparation and empirical complexity analysis (including the error scaling) of classical heuristics, in both ab initio and model Hamiltonian settings, we conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for ground-state quantum chemistry through polynomial speedups, it may be prudent to assume exponential speedups are not generically available for this problem.
Collapse
Affiliation(s)
- Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, NY, 10027, USA
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Yu Tong
- Department of Mathematics, University of California, Berkeley, CA, 94720, USA
| | | | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, VA, 24061, USA
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Phillip Helms
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Johnnie Gray
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Zhi-Hao Cui
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Wenyuan Liu
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Michael Kastoryano
- AWS Center for Quantum Computing, Pasadena, CA, 91125, USA
- Amazon Quantum Solutions Lab, Seattle, WA, 98170, USA
| | - Ryan Babbush
- Google Quantum AI, 340 Main Street, Venice, CA, 90291, USA
| | - John Preskill
- AWS Center for Quantum Computing, Pasadena, CA, 91125, USA
- Institute for Quantum Information and Matter, California Institute of Technology, Pasadena, CA, 91125, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, NY, 10027, USA
| | | | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, VA, 24061, USA
| | - Lin Lin
- Department of Mathematics, University of California, Berkeley, CA, 94720, USA
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| |
Collapse
|
16
|
Hydrogen Bonding to Graphene Surface: A Comparative Computational Study. Inorganica Chim Acta 2023. [DOI: 10.1016/j.ica.2023.121454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
|
17
|
Robert A, Berthoumieux H, Bocquet ML. Coupled Interactions at the Ionic Graphene-Water Interface. PHYSICAL REVIEW LETTERS 2023; 130:076201. [PMID: 36867792 DOI: 10.1103/physrevlett.130.076201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Accepted: 01/13/2023] [Indexed: 06/18/2023]
Abstract
We compute ionic free energy adsorption profiles at an aqueous graphene interface by developing a self-consistent approach. To do so, we design a microscopic model for water and put the liquid on an equal footing with the graphene described by its electronic band structure. By evaluating progressively the electronic and dipolar coupled electrostatic interactions, we show that the coupling level including mutual graphene and water screening permits one to recover remarkably the precision of extensive quantum simulations. We further derive the potential of mean force evolution of several alkali cations.
Collapse
Affiliation(s)
- Anton Robert
- PASTEUR, Département de chimie, École normale supérieure, Université PSL, CNRS, Sorbonne Université, 75005 Paris, France
| | - Hélène Berthoumieux
- Sorbonne Université, CNRS, Laboratoire de Physique Théorique de la Matière Condensée (LPTMC, UMR 7600), F-75005 Paris, France
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, Berlin 14195, Germany
| | - Marie-Laure Bocquet
- LPENS, École normale supérieure, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, F-75005 Paris, France
| |
Collapse
|
18
|
Vekeman J, García Cuesta I, Faginas-Lago N, Sánchez-Marín J, Sánchez de Merás AMJ. Development of accurate potentials for the physisorption of water on graphene. J Chem Phys 2023; 158:024104. [PMID: 36641401 DOI: 10.1063/5.0131626] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
From coupled-cluster singles and doubles model including connected triples corrections [CCSD(T)] calculations on the water dimer and B97D/CC on the water-circumcoronene complex at a large number of randomly generated conformations, interaction potentials for the physisorption of water on graphene are built, accomplishing almost sub-chemical accuracy. The force fields were constructed by decomposing the interaction into electrostatic and van der Waals contributions, the latter represented through improved Lennard-Jones potentials. Besides, a Chemistry at Harvard Macromolecular Mechanics (CHARMM)-like term was included in the water-water potential to improve the description of hydrogen bonds, and an induction term was added to model the polarization effects in the interaction between water and polyaromatic hydrocarbons (PAHs) or graphene. Two schemes with three and six point charges were considered for the interactions water-water and water-PAH, as Coulomb contributions are zero in the water-graphene system. The proposed fitted potentials reproduce the ab initio data used to build them in the whole range of distances and conformations and provide results for selected points very close to CCSD(T) benchmarks. When applied to the water-graphene system, the obtained results are in excellent agreement with p-CCSD(T), revised symmetry-adapted perturbation theory based on density functional theory monomer properties (DFT-SAPT), and diffusion Monte Carlo reference values. Furthermore, the stability of the various conformers water-PAH and water-graphene, as well as the different trends observed between these systems are rationalized in terms of the modifications of the electrostatic contribution.
Collapse
Affiliation(s)
- Jelle Vekeman
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark-Zwijnaarde 46, 9052 Zwijnaarde, Belgium
| | - Inmaculada García Cuesta
- Departamento de Química Física, Universidad de Valencia, Avda. Dr. Moliner 50, E-46100 Burjassot, Spain
| | - Noelia Faginas-Lago
- Dipartimento di Chimica, Biologia e Biotecnologie, Università Degli Studi di Perugia, Via Elce di Sotto 8, I-06123 Perugia, Italy
| | - José Sánchez-Marín
- Instituto de Ciencia Molecular, Universidad de Valencia, C/Catedrático José Beltrán 2, E-46980 Paterna, Spain
| | | |
Collapse
|
19
|
Garcia R. Interfacial Liquid Water on Graphite, Graphene, and 2D Materials. ACS NANO 2023; 17:51-69. [PMID: 36507725 PMCID: PMC10664075 DOI: 10.1021/acsnano.2c10215] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Accepted: 12/08/2022] [Indexed: 06/17/2023]
Abstract
The optical, electronic, and mechanical properties of graphite, few-layer, and two-dimensional (2D) materials have prompted a considerable number of applications. Biosensing, energy storage, and water desalination illustrate applications that require a molecular-scale understanding of the interfacial water structure on 2D materials. This review introduces the most recent experimental and theoretical advances on the structure of interfacial liquid water on graphite-like and 2D materials surfaces. On pristine conditions, atomic-scale resolution experiments revealed the existence of 1-3 hydration layers. Those layers were separated by ∼0.3 nm. The experimental data were supported by molecular dynamics simulations. However, under standard working conditions, atomic-scale resolution experiments revealed the presence of 2-3 hydrocarbon layers. Those layers were separated by ∼0.5 nm. Linear alkanes were the dominant molecular specie within the hydrocarbon layers. Paradoxically, the interface of an aged 2D material surface immersed in water does not have water molecules on its vicinity. Free-energy considerations favored the replacement of water by alkanes.
Collapse
Affiliation(s)
- Ricardo Garcia
- Instituto de Ciencia de Materiales
de Madrid, CSIC, c/Sor Juana Inés de la Cruz 3, 28049Madrid, Spain
| |
Collapse
|
20
|
Sacchi M, Tamtögl A. Water adsorption and dynamics on graphene and other 2D materials: Computational and experimental advances. ADVANCES IN PHYSICS: X 2022; 8:2134051. [PMID: 36816858 PMCID: PMC7614201 DOI: 10.1080/23746149.2022.2134051] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Revised: 09/26/2022] [Accepted: 09/27/2022] [Indexed: 06/18/2023] Open
Abstract
The interaction of water and surfaces, at molecular level, is of critical importance for understanding processes such as corrosion, friction, catalysis and mass transport. The significant literature on interactions with single crystal metal surfaces should not obscure unknowns in the unique behaviour of ice and the complex relationships between adsorption, diffusion and long-range inter-molecular interactions. Even less is known about the atomic-scale behaviour of water on novel, non-metallic interfaces, in particular on graphene and other 2D materials. In this manuscript, we review recent progress in the characterisation of water adsorption on 2D materials, with a focus on the nano-material graphene and graphitic nanostructures; materials which are of paramount importance for separation technologies, electrochemistry and catalysis, to name a few. The adsorption of water on graphene has also become one of the benchmark systems for modern computational methods, in particular dispersion-corrected density functional theory (DFT). We then review recent experimental and theoretical advances in studying the single-molecular motion of water at surfaces, with a special emphasis on scattering approaches as they allow an unparalleled window of observation to water surface motion, including diffusion, vibration and self-assembly.
Collapse
Affiliation(s)
- M. Sacchi
- Department of Chemistry, University of Surrey, Guildford GU2 7XH, UK
| | - A. Tamtögl
- Institute of Experimental Physics, Graz University of Technology, 8010 Graz, Austria
| |
Collapse
|
21
|
Nanoconfinement facilitates reactions of carbon dioxide in supercritical water. Nat Commun 2022; 13:5932. [PMID: 36209274 PMCID: PMC9547913 DOI: 10.1038/s41467-022-33696-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 09/27/2022] [Indexed: 11/08/2022] Open
Abstract
The reactions of CO2 in water under extreme pressure-temperature conditions are of great importance to the carbon storage and transport below Earth’s surface, which substantially affect the carbon budget in the atmosphere. Previous studies focus on the CO2(aq) solutions in the bulk phase, but underground aqueous solutions are often confined to the nanoscale, and nanoconfinement and solid-liquid interfaces may substantially affect chemical speciation and reaction mechanisms, which are poorly known on the molecular scale. Here, we apply extensive ab initio molecular dynamics simulations to study aqueous carbon solutions nanoconfined by graphene and stishovite (SiO2) at 10 GPa and 1000 ~ 1400 K. We find that CO2(aq) reacts more in nanoconfinement than in bulk. The stishovite-water interface makes the solutions more acidic, which shifts the chemical equilibria, and the interface chemistry also significantly affects the reaction mechanisms. Our findings suggest that CO2(aq) in deep Earth is more active than previously thought, and confining CO2 and water in nanopores may enhance the efficiency of mineral carbonation. Aqueous CO2 under nanoconfinement is of great importance to the carbon storage and transport in Earth. Here, the authors apply ab initio molecular dynamics simulations to study the effects of confinement and interfaces, and show that that CO(aq) reacts more in nanoconfinement than in bulk.
Collapse
|
22
|
The first-principles phase diagram of monolayer nanoconfined water. Nature 2022; 609:512-516. [PMID: 36104556 DOI: 10.1038/s41586-022-05036-x] [Citation(s) in RCA: 32] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Accepted: 06/28/2022] [Indexed: 11/09/2022]
Abstract
Water in nanoscale cavities is ubiquitous and of central importance to everyday phenomena in geology and biology. However, the properties of nanoscale water can be substantially different from those of bulk water, as shown, for example, by the anomalously low dielectric constant of water in nanochannels1, near frictionless water flow2 or the possible existence of a square ice phase3. Such properties suggest that nanoconfined water could be engineered for technological applications in nanofluidics4, electrolyte materials5 and water desalination6. Unfortunately, challenges in experimentally characterizing water at the nanoscale and the high cost of first-principles simulations have prevented the molecular-level understanding required to control the behaviour of water. Here we combine a range of computational approaches to enable a first-principles-level investigation of a single layer of water within a graphene-like channel. We find that monolayer water exhibits surprisingly rich and diverse phase behaviour that is highly sensitive to temperature and the van der Waals pressure acting within the nanochannel. In addition to multiple molecular phases with melting temperatures varying non-monotonically by more than 400 kelvins with pressure, we predict a hexatic phase, which is an intermediate between a solid and a liquid, and a superionic phase with a high electrical conductivity exceeding that of battery materials. Notably, this suggests that nanoconfinement could be a promising route towards superionic behaviour under easily accessible conditions.
Collapse
|
23
|
Hutchison P, Rice PS, Warburton RE, Raugei S, Hammes-Schiffer S. Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe-N-C Materials. J Am Chem Soc 2022; 144:16524-16534. [PMID: 36001092 DOI: 10.1021/jacs.2c05779] [Citation(s) in RCA: 26] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The systematic improvement of Fe-N-C materials for fuel cell applications has proven challenging, due in part to an incomplete atomistic understanding of the oxygen reduction reaction (ORR) under electrochemical conditions. Herein, a multilevel computational approach, which combines ab initio molecular dynamics simulations and constant potential density functional theory calculations, is used to assess proton-coupled electron transfer (PCET) processes and adsorption thermodynamics of key ORR intermediates. These calculations indicate that the potential-limiting step for ORR on Fe-N-C materials is the formation of the FeIII-OOH intermediate. They also show that an active site model with a water molecule axially ligated to the iron center throughout the catalytic cycle produces results that are consistent with the experimental measurements. In particular, reliable prediction of the ORR onset potential and the Fe(III/II) redox potential associated with the conversion of FeIII-OH to FeII and desorbed H2O requires an axial H2O co-adsorbed to the iron center. The observation of a five-coordinate rather than four-coordinate active site has significant implications for the thermodynamics and mechanism of ORR. These findings highlight the importance of solvent-substrate interactions and surface charge effects for understanding the PCET reaction mechanisms and transition-metal redox couples under realistic electrochemical conditions.
Collapse
Affiliation(s)
- Phillips Hutchison
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Peter S Rice
- Center for Molecular Electrocatalysis, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Robert E Warburton
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Simone Raugei
- Center for Molecular Electrocatalysis, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | |
Collapse
|
24
|
Thiemann F, Schran C, Rowe P, Müller EA, Michaelides A. Water Flow in Single-Wall Nanotubes: Oxygen Makes It Slip, Hydrogen Makes It Stick. ACS NANO 2022; 16:10775-10782. [PMID: 35726839 PMCID: PMC9331139 DOI: 10.1021/acsnano.2c02784] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
Experimental measurements have reported ultrafast and radius-dependent water transport in carbon nanotubes which are absent in boron nitride nanotubes. Despite considerable effort, the origin of this contrasting (and fascinating) behavior is not understood. Here, with the aid of machine learning-based molecular dynamics simulations that deliver first-principles accuracy, we investigate water transport in single-wall carbon and boron nitride nanotubes. Our simulations reveal a large, radius-dependent hydrodynamic slippage on both materials, with water experiencing indeed a ≈5 times lower friction on carbon surfaces compared to boron nitride. Analysis of the diffusion mechanisms across the two materials reveals that the fast water transport on carbon is governed by facile oxygen motion, whereas the higher friction on boron nitride arises from specific hydrogen-nitrogen interactions. This work not only delivers a clear reference of quantum mechanical accuracy for water flow in single-wall nanotubes but also provides detailed mechanistic insight into its radius and material dependence for future technological application.
Collapse
Affiliation(s)
- Fabian
L. Thiemann
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, Gower Street, London WC1E 6BT, United
Kingdom
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Christoph Schran
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, Gower Street, London WC1E 6BT, United
Kingdom
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Patrick Rowe
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, Gower Street, London WC1E 6BT, United
Kingdom
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Erich A. Müller
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Angelos Michaelides
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, Gower Street, London WC1E 6BT, United
Kingdom
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
25
|
Iyer GR, Rubenstein BM. Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry. J Phys Chem A 2022; 126:4636-4646. [PMID: 35820033 DOI: 10.1021/acs.jpca.2c01957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H2O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.
Collapse
Affiliation(s)
- Gopal R Iyer
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
26
|
Dumi A, Upadhyay S, Bernasconi L, Shin H, Benali A, Jordan KD. The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations. J Chem Phys 2022; 156:144702. [DOI: 10.1063/5.0085982] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.
Collapse
Affiliation(s)
- Amanda Dumi
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Shiv Upadhyay
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Leonardo Bernasconi
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
- Center for Research Computing, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Hyeondeok Shin
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
27
|
Gasparotto P, Fitzner M, Cox SJ, Sosso GC, Michaelides A. How do interfaces alter the dynamics of supercooled water? NANOSCALE 2022; 14:4254-4262. [PMID: 35244128 DOI: 10.1039/d2nr00387b] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The structure of liquid water in the proximity of an interface can deviate significantly from that of bulk water, with surface-induced structural perturbations typically converging to bulk values at about ∼1 nm from the interface. While these structural changes are well established it is, in contrast, less clear how an interface perturbs the dynamics of water molecules within the liquid. Here, through an extensive set of molecular dynamics simulations of supercooled bulk and interfacial water films and nano-droplets, we observe the formation of persistent, spatially extended dynamical domains in which the average mobility varies as a function of the distance from the interface. This is in stark contrast with the dynamical heterogeneity observed in bulk water, where these domains average out spatially over time. We also find that the dynamical response of water to an interface depends critically on the nature of the interface and on the choice of interface definition. Overall these results reveal a richness in the dynamics of interfacial water that opens up the prospect of tuning the dynamical response of water through specific modifications of the interface structure or confining material.
Collapse
Affiliation(s)
- Piero Gasparotto
- Scientific Computing Division, Paul Scherrer Institute, Villigen 5232, Switzerland.
| | - Martin Fitzner
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, London WC1E 6BT, UK
| | - Stephen James Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | - Gabriele Cesare Sosso
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| |
Collapse
|
28
|
Keller E, Tsatsoulis T, Reuter K, Margraf JT. Regularized second-order correlation methods for extended systems. J Chem Phys 2022; 156:024106. [PMID: 35032995 DOI: 10.1063/5.0078119] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Second-order Møller-Plesset perturbation theory (MP2) constitutes the simplest form of many-body wavefunction theory and often provides a good compromise between efficiency and accuracy. There are, however, well-known limitations to this approach. In particular, MP2 is known to fail or diverge for some prototypical condensed matter systems like the homogeneous electron gas (HEG) and to overestimate dispersion-driven interactions in strongly polarizable systems. In this paper, we explore how the issues of MP2 for metallic, polarizable, and strongly correlated periodic systems can be ameliorated through regularization. To this end, two regularized second-order methods (including a new, size-extensive Brillouin-Wigner approach) are applied to the HEG, the one-dimensional Hubbard model, and the graphene-water interaction. We find that regularization consistently leads to improvements over the MP2 baseline and that different regularizers are appropriate for the various systems.
Collapse
Affiliation(s)
- Elisabeth Keller
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Theodoros Tsatsoulis
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Johannes T Margraf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
29
|
Keal TW, Elena AM, Sokol AA, Stoneham K, Probert MIJ, Cucinotta CS, Willock DJ, Logsdail AJ, Zen A, Hasnip PJ, Bush IJ, Watkins M, Alfe D, Skylaris CK, Curchod BFE, Cai Q, Woodley SM. Materials and Molecular Modeling at the Exascale. Comput Sci Eng 2022. [DOI: 10.1109/mcse.2022.3141328] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
| | | | | | | | | | | | | | | | | | | | - Ian J. Bush
- STFC Rutherford Appleton Laboratory, Oxfordshire, U.K
| | | | | | | | | | - Qiong Cai
- University of Surrey, Guildford, U.K
| | | |
Collapse
|
30
|
Schäfer T, Gallo A, Irmler A, Hummel F, Grüneis A. Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: The case of water on graphitic carbon nitride. J Chem Phys 2021; 155:244103. [PMID: 34972356 DOI: 10.1063/5.0074936] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride is reported employing an embedding approach for many-electron correlation methods. To this end, a plane-wave based implementation to obtain intrinsic atomic orbitals and Wannier functions for arbitrary localization potentials is presented. In our embedding scheme, the localized occupied orbitals allow for a separate treatment of short-range and long-range correlation contributions to the adsorption energy by a fragmentation of the simulation cell. In combination with unoccupied natural orbitals, the coupled cluster ansatz with single, double, and perturbative triple particle-hole excitation operators is used to capture the correlation in local fragments centered around the adsorption process. For the long-range correlation, a seamless embedding into the random phase approximation yields rapidly convergent adsorption energies with respect to the local fragment size. Convergence of computed binding energies with respect to the virtual orbital basis set is achieved employing a number of recently developed techniques. Moreover, we discuss fragment size convergence for a range of approximate many-electron perturbation theories. The obtained benchmark results are compared to a number of density functional calculations.
Collapse
Affiliation(s)
- Tobias Schäfer
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Alejandro Gallo
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Andreas Irmler
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Felix Hummel
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Andreas Grüneis
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| |
Collapse
|
31
|
Barbosa N, Pagliai M, Sinha S, Barone V, Alfè D, Brancato G. Enhancing the Accuracy of Ab Initio Molecular Dynamics by Fine Tuning of Effective Two-Body Interactions: Acetonitrile as a Test Case. J Phys Chem A 2021; 125:10475-10484. [PMID: 34843249 DOI: 10.1021/acs.jpca.1c07576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Grimme's dispersion-corrected density functional theory (DFT-D) methods have emerged among the most practical approaches to perform accurate quantum mechanical calculations on molecular systems ranging from small clusters to microscopic and mesoscopic samples, i.e., including hundreds or thousands of molecules. Moreover, DFT-D functionals can be easily integrated into popular ab initio molecular dynamics (MD) software packages to carry out first-principles condensed-phase simulations at an affordable computational cost. Here, starting from the well-established D3 version of the dispersion-correction term, we present a simple protocol to improve the accurate description of the intermolecular interactions of molecular clusters of growing size, considering acetonitrile as a test case. Optimization of the interaction energy was performed with reference to diffusion quantum Monte Carlo calculations, successfully reaching the same inherent accuracy of the latter (statistical error of ∼0.1 kcal/mol per molecule). The refined DFT-D3 model was then used to perform ab initio MD simulations of liquid acetonitrile, again showing significant improvements toward available experimental data with respect to the default correction.
Collapse
Affiliation(s)
- Nuno Barbosa
- Scuola Normale Superiore and CSGI, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Marco Pagliai
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| | - Sourab Sinha
- Scuola Normale Superiore and CSGI, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore and CSGI, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Dario Alfè
- Department of Earth Sciences, Thomas Young Center, University College London, 5 Gower Place, WC1E 6BS London, United Kingdom.,London Centre for Nanotechnology, Thomas Young Centre, University College London, 17-19 Gordon Street, WC1H 0AH London, United Kingdom.,Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, 80126 Napoli, Italy
| | - Giuseppe Brancato
- Scuola Normale Superiore and CSGI, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| |
Collapse
|
32
|
Joly L, Meißner RH, Iannuzzi M, Tocci G. Osmotic Transport at the Aqueous Graphene and hBN Interfaces: Scaling Laws from a Unified, First-Principles Description. ACS NANO 2021; 15:15249-15258. [PMID: 34491721 DOI: 10.1021/acsnano.1c05931] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Osmotic transport in nanoconfined aqueous electrolytes provides alternative venues for water desalination and "blue energy" harvesting. The osmotic response of nanofluidic systems is controlled by the interfacial structure of water and electrolyte solutions in the so-called electrical double layer (EDL), but a molecular-level picture of the EDL is to a large extent still lacking. Particularly, the role of the electronic structure has not been considered in the description of electrolyte/surface interactions. Here, we report enhanced sampling simulations based on ab initio molecular dynamics, aiming at unravelling the free energy of prototypical ions adsorbed at the aqueous graphene and hBN interfaces, and its consequences on nanofluidic osmotic transport. Specifically, we predicted the zeta potential, the diffusio-osmotic mobility, and the diffusio-osmotic conductivity for a wide range of salt concentrations from the ab initio water and ion spatial distributions through an analytical framework based on Stokes equation and a modified Poisson-Boltzmann equation. We observed concentration-dependent scaling laws, together with dramatic differences in osmotic transport between the two interfaces, including diffusio-osmotic flow and current reversal on hBN but not on graphene. We could rationalize the results for the three osmotic responses with a simple model based on characteristic length scales for ion and water adsorption at the surface, which are quite different on graphene and on hBN. Our work provides fundamental insights into the structure and osmotic transport of aqueous electrolytes on 2D materials and explores alternative pathways for efficient water desalination and osmotic energy conversion.
Collapse
Affiliation(s)
- Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
- Institut Universitaire de France (IUF), 1 rue Descartes, 75005 Paris, France
| | - Robert H Meißner
- Hamburg University of Technology, Insitute of Polymers and Composites, Hamburg 21073, Germany
- Helmholtz-Zentrum Hereon, Institute of Surface Science, Geesthacht 21502, Germany
| | - Marcella Iannuzzi
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| |
Collapse
|
33
|
Galimberti DR, Sauer J. Chemically Accurate Vibrational Free Energies of Adsorption from Density Functional Theory Molecular Dynamics: Alkanes in Zeolites. J Chem Theory Comput 2021; 17:5849-5862. [PMID: 34459582 PMCID: PMC8444336 DOI: 10.1021/acs.jctc.1c00519] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present a methodology
to compute, at reduced computational cost,
Gibbs free energies, enthalpies, and entropies of adsorption from
molecular dynamics. We calculate vibrational partition functions from
vibrational energies, which we obtain from the vibrational density
of states by projection on the normal modes. The use of a set of well-chosen
reference structures along the trajectories accounts for the anharmonicities
of the modes. For the adsorption of methane, ethane, and propane in
the H-CHA zeolite, we limit our treatment to a set of vibrational
modes localized at the adsorption site (zeolitic OH group) and the
alkane molecule interacting with it. Only two short trajectories (1–20
ps) are required to reach convergence (<1 kJ/mol) for the thermodynamic
functions. The mean absolute deviations from the experimentally measured
values are 2.6, 2.8, and 4.7 kJ/mol for the Gibbs free energy, the
enthalpy, and the entropy term (−TΔS),
respectively. In particular, the entropy terms show a major improvement
compared to the harmonic approximation and almost reach the accuracy
of the previous use of anharmonic frequencies obtained with curvilinear
distortions of individual modes. The thermodynamic functions so obtained
follow the trend of the experimental values for methane, ethane, and
propane, and the Gibbs free energy of adsorption at experimental conditions
is correctly predicted to change from positive for methane (5.9 kJ/mol)
to negative for ethane (−4.8 kJ/mol) and propane (−7.1
kJ/mol).
Collapse
Affiliation(s)
- Daria Ruth Galimberti
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany.,Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany
| |
Collapse
|
34
|
Ho TA, Wang Y. Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations. ACS APPLIED MATERIALS & INTERFACES 2021; 13:41330-41338. [PMID: 34410713 DOI: 10.1021/acsami.1c11540] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO2 or CH4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.
Collapse
Affiliation(s)
- Tuan A Ho
- Geochemistry Department, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Yifeng Wang
- Nuclear Waste Disposal Research and Analysis Department, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| |
Collapse
|
35
|
Interactions between large molecules pose a puzzle for reference quantum mechanical methods. Nat Commun 2021; 12:3927. [PMID: 34168142 PMCID: PMC8225865 DOI: 10.1038/s41467-021-24119-3] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2020] [Accepted: 06/02/2021] [Indexed: 02/05/2023] Open
Abstract
Quantum-mechanical methods are used for understanding molecular interactions throughout the natural sciences. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are state-of-the-art trusted wavefunction methods that have been shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. We show that CCSD(T) and DMC interaction energies are not consistent for a set of polarizable supramolecules. Whilst there is agreement for some of the complexes, in a few key systems disagreements of up to 8 kcal mol-1 remain. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules.
Collapse
|
36
|
Callahan JM, Lange MF, Berkelbach TC. Dynamical correlation energy of metals in large basis sets from downfolding and composite approaches. J Chem Phys 2021; 154:211105. [PMID: 34240964 DOI: 10.1063/5.0049890] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Coupled-cluster theory with single and double excitations (CCSD) is a promising ab initio method for the electronic structure of three-dimensional metals, for which second-order perturbation theory (MP2) diverges in the thermodynamic limit. However, due to the high cost and poor convergence of CCSD with respect to basis size, applying CCSD to periodic systems often leads to large basis set errors. In a common "composite" method, MP2 is used to recover the missing dynamical correlation energy through a focal-point correction, but the inadequacy of finite-order perturbation theory for metals raises questions about this approach. Here, we describe how high-energy excitations treated by MP2 can be "downfolded" into a low-energy active space to be treated by CCSD. Comparing how the composite and downfolding approaches perform for the uniform electron gas, we find that the latter converges more quickly with respect to the basis set size. Nonetheless, the composite approach is surprisingly accurate because it removes the problematic MP2 treatment of double excitations near the Fermi surface. Using this method to estimate the CCSD correlation energy in the combined complete basis set and thermodynamic limits, we find that CCSD recovers 85%-90% of the exact correlation energy at rs = 4. We also test the composite approach with the direct random-phase approximation used in place of MP2, yielding a method that is typically (but not always) more cost effective due to the smaller number of orbitals that need to be included in the more expensive CCSD calculation.
Collapse
Affiliation(s)
- James M Callahan
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Malte F Lange
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
37
|
Kim D, Kim E, Park S, Kim S, Min BK, Yoon HJ, Kwak K, Cho M. Wettability of graphene and interfacial water structure. Chem 2021. [DOI: 10.1016/j.chempr.2021.03.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
38
|
dos Santos TC, Mancera RC, Rocha MV, da Silva AF, Furtado IO, Barreto J, Stavale F, Archanjo BS, de M. Carneiro JW, Costa LT, Ronconi CM. CO2 and H2 adsorption on 3D nitrogen-doped porous graphene: Experimental and theoretical studies. J CO2 UTIL 2021. [DOI: 10.1016/j.jcou.2021.101517] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
39
|
Pham TT, Pham TN, Chihaia V, Vu QA, Trinh TT, Pham TT, Van Thang L, Son DN. How do the doping concentrations of N and B in graphene modify the water adsorption? RSC Adv 2021; 11:19560-19568. [PMID: 35479230 PMCID: PMC9033564 DOI: 10.1039/d1ra01506k] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2021] [Accepted: 05/14/2021] [Indexed: 12/17/2022] Open
Abstract
Understanding the interaction of water and graphene is crucial for various applications such as water purification, desalination, and electrocatalysis. Experimental and theoretical studies have already investigated water adsorption on N- and B-doped graphene. However, there are no reports available that elucidate the influences of the N and B doping content in graphene on the microscopic geometrical structure and the electronic properties of the adsorbed water. Thus, this work is devoted to solving this problem using self-consistent van der Waals density functional theory calculations. The N and B doping contents of 0.0, 3.1, 6.3, and 9.4% were considered. The results showed that the binding energy of water increases almost linearly as a function of doping content at all concentrations for N-doped graphene but below 6.3% for B-doped graphene. In the linear range, the binding energy increases by approximately 30 meV for each increment of the doping ratio. Analyses of the geometric and electronic structures explained the enhancement of the water-graphene interaction with the variation in doping percentage.
Collapse
Affiliation(s)
- Thi Tan Pham
- Ho Chi Minh City University of Technology 268 Ly Thuong Kiet Street, Ward 14, District 10 Ho Chi Minh City Vietnam .,Vietnam National University Ho Chi Minh City Quarter 6, Linh Trung Ward, Thu Duc District Ho Chi Minh City Vietnam
| | - Thanh Ngoc Pham
- Ho Chi Minh City University of Technology 268 Ly Thuong Kiet Street, Ward 14, District 10 Ho Chi Minh City Vietnam .,Vietnam National University Ho Chi Minh City Quarter 6, Linh Trung Ward, Thu Duc District Ho Chi Minh City Vietnam
| | - Viorel Chihaia
- Institute of Physical Chemistry "Ilie Murgulescu" of the Romanian Academy Splaiul Independentei 202, Sector 6 060021 Bucharest Romania
| | - Quang Anh Vu
- Ho Chi Minh City University of Technology 268 Ly Thuong Kiet Street, Ward 14, District 10 Ho Chi Minh City Vietnam .,Vietnam National University Ho Chi Minh City Quarter 6, Linh Trung Ward, Thu Duc District Ho Chi Minh City Vietnam
| | - Thuat T Trinh
- Department of Civil and Environmental Engineering, Norwegian University of Science and Technology NO-7491 Trondheim Norway
| | - Trung Thanh Pham
- Namur Institute of Structured Matter (NISM), Department of Physics, University of Namur 61 Rue de Bruxelles B-5000 Namur Belgium
| | - Le Van Thang
- Ho Chi Minh City University of Technology 268 Ly Thuong Kiet Street, Ward 14, District 10 Ho Chi Minh City Vietnam .,Vietnam National University Ho Chi Minh City Quarter 6, Linh Trung Ward, Thu Duc District Ho Chi Minh City Vietnam
| | - Do Ngoc Son
- Ho Chi Minh City University of Technology 268 Ly Thuong Kiet Street, Ward 14, District 10 Ho Chi Minh City Vietnam .,Vietnam National University Ho Chi Minh City Quarter 6, Linh Trung Ward, Thu Duc District Ho Chi Minh City Vietnam
| |
Collapse
|
40
|
Water friction in nanofluidic channels made from two-dimensional crystals. Nat Commun 2021; 12:3092. [PMID: 34035239 PMCID: PMC8149694 DOI: 10.1038/s41467-021-23325-3] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2020] [Accepted: 04/20/2021] [Indexed: 11/08/2022] Open
Abstract
Membrane-based applications such as osmotic power generation, desalination and molecular separation would benefit from decreasing water friction in nanoscale channels. However, mechanisms that allow fast water flows are not fully understood yet. Here we report angstrom-scale capillaries made from atomically flat crystals and study the effect of confining walls' material on water friction. A massive difference is observed between channels made from isostructural graphite and hexagonal boron nitride, which is attributed to different electrostatic and chemical interactions at the solid-liquid interface. Using precision microgravimetry and ion streaming measurements, we evaluate the slip length, a measure of water friction, and investigate its possible links with electrical conductivity, wettability, surface charge and polarity of the confining walls. We also show that water friction can be controlled using hybrid capillaries with different slip lengths at opposing walls. The reported advances extend nanofluidics' toolkit for designing smart membranes and mimicking manifold machinery of biological channels.
Collapse
|
41
|
Abstract
AbstractGraphene as a two-dimensional material is prone to hydrocarbon contaminations, which can significantly alter its intrinsic electrical properties. Herein, we implement a facile hydrogenation-dehydrogenation strategy to remove hydrocarbon contaminations and preserve the excellent transport properties of monolayer graphene. Using electron microscopy we quantitatively characterized the improved cleanness of hydrogenated graphene compared to untreated samples. In situ spectroscopic investigations revealed that the hydrogenation treatment promoted the adsorption ofytyt water at the graphene surface, resulting in a protective layer against the re-deposition of hydrocarbon molecules. Additionally, the further dehydrogenation of hydrogenated graphene rendered a more pristine-like basal plane with improved carrier mobility compared to untreated pristine graphene. Our findings provide a practical post-growth cleaning protocol for graphene with maintained surface cleanness and lattice integrity to systematically carry a range of surface chemistry in the form of a well-performing and reproducible transistor.
Collapse
|
42
|
Goujon F, Ghoufi A, Malfreyt P. Associated molecular liquids at the graphene monolayer interface. J Chem Phys 2021; 154:104504. [PMID: 33722040 DOI: 10.1063/5.0042438] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
We report molecular simulations of the interaction between a graphene sheet and different liquids such as water, ethanol, and ethylene glycol. We describe the structural arrangements at the graphene interface in terms of density profiles, number of hydrogen bonds (HBs), and local structuration in neighboring layers close to the surface. We establish the formation of a two-dimensional HB network in the layer closest to the graphene. We also calculate the interfacial tension of liquids with a graphene monolayer and its profile along the direction normal to the graphene to rationalize and quantify the strengthening of the intermolecular interactions in the liquid due to the presence of the surface.
Collapse
Affiliation(s)
- Florent Goujon
- Université Clermont Auvergne, CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand (ICCF), F-63000 Clermont-Ferrand, France
| | - Aziz Ghoufi
- Université de Rennes, CNRS, IPR (Institut de Physique de Rennes), UMR 6251, F-35000 Rennes, France
| | - Patrice Malfreyt
- Université Clermont Auvergne, CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand (ICCF), F-63000 Clermont-Ferrand, France
| |
Collapse
|
43
|
Ramberger B, Kresse G. New insights into the 1D carbon chain through the RPA. Phys Chem Chem Phys 2021; 23:5254-5260. [PMID: 33629671 DOI: 10.1039/d0cp06607a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigated the electronic and structural properties of the infinite linear carbon chain (carbyne) using density functional theory (DFT) and the random phase approximation (RPA) to the correlation energy. The studies are performed in vacuo and for carbyne inside a carbon nano tube (CNT). In the vacuum, semi-local DFT and RPA predict bond length alternations of about 0.04 Å and 0.13 Å, respectively. The frequency of the highest optical mode at the Γ point is 1219 cm-1 and about 2000 cm-1 for DFT and the RPA. Agreement of the RPA to previous high level quantum chemistry and diffusion Monte-Carlo results is excellent. For the RPA we calculate the phonon-dispersion in the full Brillouine zone and find marked quantitative differences to DFT calculations not only at the Γ point but also throughout the entire Brillouine zone. To model carbyne inside a carbon nanotube, we considered a (10,0) CNT. Here the DFT calculations are even qualitatively sensitive to the k-points sampling. At the limes of a very dense k-points sampling, semi-local DFT predicts no bond length alternation (BLA), whereas in the RPA a sizeable BLA of 0.09 Å prevails. The reduced BLA leads to a significant red shift of the vibrational frequencies of about 350 cm-1, so that they are in good agreement with experimental estimates. Overall, the good agreement between the RPA and previously reported results from correlated wavefunction methods and experimental Raman data suggests that the RPA provides reliable results at moderate computational costs. It hence presents a useful addition to the repertoire of correlated wavefunction methods and its accuracy clearly prevails for low dimensional systems, where semi-local density functionals struggle to yield even qualitatively correct results.
Collapse
Affiliation(s)
- Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Kolingasse 14-16, 1090 Vienna, Austria.
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Kolingasse 14-16, 1090 Vienna, Austria.
| |
Collapse
|
44
|
Lau BTG, Knizia G, Berkelbach TC. Regional Embedding Enables High-Level Quantum Chemistry for Surface Science. J Phys Chem Lett 2021; 12:1104-1109. [PMID: 33475362 DOI: 10.1021/acs.jpclett.0c03274] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Compared to common density functionals, ab initio wave function methods can provide greater reliability and accuracy, which could prove useful when modeling adsorbates or defects of otherwise periodic systems. However, the breaking of translational symmetry necessitates large supercells that are often prohibitive for correlated wave function methods. As an alternative, this paper introduces the regional embedding approach, which enables correlated wave function treatments of only a target fragment of interest through small, fragment-localized orbital spaces constructed using a simple overlap criterion. Applications to the adsorption of water on lithium hydride, hexagonal boron nitride, and graphene substrates show that regional embedding combined with focal-point corrections can provide converged CCSD(T) (coupled-cluster) adsorption energies with very small fragment sizes.
Collapse
Affiliation(s)
- Bryan T G Lau
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Gerald Knizia
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Timothy C Berkelbach
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
45
|
Zhao Z, Hou T, Wu N, Jiao S, Zhou K, Yin J, Suk JW, Cui X, Zhang M, Li S, Qu Y, Xie W, Li XB, Zhao C, Fu Y, Hong RD, Guo S, Lin D, Cai W, Mai W, Luo Z, Tian Y, Lai Y, Liu Y, Colombo L, Hao Y. Polycrystalline Few-Layer Graphene as a Durable Anticorrosion Film for Copper. NANO LETTERS 2021; 21:1161-1168. [PMID: 33411539 DOI: 10.1021/acs.nanolett.0c04724] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Corrosion of metals in atmospheric environments is a worldwide problem in industry and daily life. Traditional anticorrosion methods including sacrificial anodes or protective coatings have performance limitations. Here, we report atomically thin, polycrystalline few-layer graphene (FLG) grown by chemical vapor deposition as a long-term protective coating film for copper (Cu). A six-year old, FLG-protected Cu is visually shiny and detailed material characterizations capture no sign of oxidation. The success of the durable anticorrosion film depends on the misalignment of grain boundaries between adjacent graphene layers. Theoretical calculations further found that corrosive molecules always encounter extremely high energy barrier when diffusing through the FLG layers. Therefore, the FLG is able to prevent the corrosive molecules from reaching the underlying Cu surface. This work highlights the interesting structures of polycrystalline FLG and sheds insight into the atomically thin coatings for various applications.
Collapse
Affiliation(s)
- Zhijuan Zhao
- Siyuan Laboratory, Guangzhou Key Laboratory of Vacuum Coating Technologies and New Energy Materials, Guangdong Provincial Engineering Technology Research Center of Vacuum Coating Technologies and New Energy Materials, Department of Physics, Jinan University, Guangzhou, Guangdong 510632, China
| | - Tianyu Hou
- National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Jiangsu Key Laboratory of Artificial Functional Materials, School of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210023, China
| | - Nannan Wu
- National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Jiangsu Key Laboratory of Artificial Functional Materials, School of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210023, China
| | - Shuping Jiao
- Shanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai Institute of Applied Mathematics and Mechanics, School of Mechanics, and Engineering Science, Shanghai University, Shanghai, 200444, China
| | - Ke Zhou
- State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Jun Yin
- State Key Laboratory of Mechanics and Control of Mechanical Structures and MOE Key Laboratory for Intelligent Nano Materials and Devices, College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
| | - Ji Won Suk
- School of Mechanical Engineering and SKKU Advanced Institute of Nanotechnology (SAINT), Sungkyunkwan University, Suwon, Gyeonggi-do 16419, Republic of Korea
| | - Xu Cui
- AutoX Technologies Inc., San Jose, California 95131, United States
| | - Mingfei Zhang
- National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Jiangsu Key Laboratory of Artificial Functional Materials, School of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210023, China
| | - Shaopeng Li
- National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Jiangsu Key Laboratory of Artificial Functional Materials, School of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210023, China
| | - Yan Qu
- National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Jiangsu Key Laboratory of Artificial Functional Materials, School of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210023, China
- The Sixth Element Materials Technology Co., Ltd., Changzhou 213000, China
| | - Weiguang Xie
- Siyuan Laboratory, Guangzhou Key Laboratory of Vacuum Coating Technologies and New Energy Materials, Guangdong Provincial Engineering Technology Research Center of Vacuum Coating Technologies and New Energy Materials, Department of Physics, Jinan University, Guangzhou, Guangdong 510632, China
| | - Xi-Bo Li
- Siyuan Laboratory, Guangzhou Key Laboratory of Vacuum Coating Technologies and New Energy Materials, Guangdong Provincial Engineering Technology Research Center of Vacuum Coating Technologies and New Energy Materials, Department of Physics, Jinan University, Guangzhou, Guangdong 510632, China
| | - Chuanxi Zhao
- Siyuan Laboratory, Guangzhou Key Laboratory of Vacuum Coating Technologies and New Energy Materials, Guangdong Provincial Engineering Technology Research Center of Vacuum Coating Technologies and New Energy Materials, Department of Physics, Jinan University, Guangzhou, Guangdong 510632, China
| | - Yong Fu
- Siyuan Laboratory, Guangzhou Key Laboratory of Vacuum Coating Technologies and New Energy Materials, Guangdong Provincial Engineering Technology Research Center of Vacuum Coating Technologies and New Energy Materials, Department of Physics, Jinan University, Guangzhou, Guangdong 510632, China
| | - Rong-Dun Hong
- Department of Physics and State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
| | - Shengshi Guo
- Department of Physics and State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
| | - Dingqu Lin
- Department of Physics and State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
| | - Weiwei Cai
- Department of Physics and State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
| | - Wenjie Mai
- Siyuan Laboratory, Guangzhou Key Laboratory of Vacuum Coating Technologies and New Energy Materials, Guangdong Provincial Engineering Technology Research Center of Vacuum Coating Technologies and New Energy Materials, Department of Physics, Jinan University, Guangzhou, Guangdong 510632, China
| | - Zhengtang Luo
- Department of Chemical and Biological Engineering, William Mong Institute of Nano Science and Technology, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong 999077, China
| | - Yongtao Tian
- School of Physics and Microelectronics, Zhengzhou University, Zhengzhou, Henan 450001, China
| | - Yun Lai
- National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Jiangsu Key Laboratory of Artificial Functional Materials, School of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210023, China
| | - Yuanyue Liu
- Texas Materials Institute and Department of Mechanical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Luigi Colombo
- Department of Materials Science and Engineering, The University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Yufeng Hao
- National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Jiangsu Key Laboratory of Artificial Functional Materials, School of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210023, China
- Haian Institute of New Technology, Nanjing University, Haian, 226600, China
| |
Collapse
|
46
|
Paniagua-Guerra LE, Gonzalez-Valle CU, Ramos-Alvarado B. Effects of the Interfacial Modeling Approach on Equilibrium Calculations of Slip Length for Nanoconfined Water in Carbon Slits. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:14772-14781. [PMID: 33215929 DOI: 10.1021/acs.langmuir.0c02718] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
In this investigation, equilibrium molecular dynamics simulations were conducted to assess the influence of the interface modeling approach on the calculation of hydrodynamic slip in carbon nanochannels. A Green-Kubo formalism was implemented for the calculation of the slip length in water confined by graphite layers. The nonbonded interactions between solid and liquid atoms (interface models) were modeled using parameters optimized to represent the wetting behavior and adsorption energy curves from electronic structure calculations. Conventional carbon-oxygen-only interaction models were compared against comprehensive models able to represent the molecular-orientation-dependent energy of interaction. Quasi-universal relationships built under the premise of the slip length dependence on the water-graphite affinity and characterized by macroscopic wettability were critically assessed. It was found that the wetting behavior cannot fully characterize the hydrodynamic slip because interface models that produced the same surface wettability yielded different values of the friction coefficient. Alternatively, the density depletion length, used to characterize the interfacial liquid structuring and the availability of momentum carriers (interfacial water molecules), was able to accurately represent the slip length trends independently of the interface model. These findings reassert the importance of physically sound interface models to study interfacial transport properties and the need of reliable parameters and characterization procedures to support theoretical models that seek to unveil the inconsistencies in hydrodynamic slip calculations.
Collapse
Affiliation(s)
- Luis E Paniagua-Guerra
- Department of Mechanical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - C Ulises Gonzalez-Valle
- Department of Mechanical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Bladimir Ramos-Alvarado
- Department of Mechanical Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
47
|
Computational Investigations of Dispersion Interactions between Small Molecules and Graphene-like Flakes. J Phys Chem A 2020; 124:9552-9561. [PMID: 33166136 DOI: 10.1021/acs.jpca.0c06595] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate dispersion interactions in a selection of atomic, molecular, and molecule-surface systems, comparing high-level correlated methods with empirically corrected density functional theory (DFT). We assess the efficacy of functionals commonly used for surface-based calculations, with and without the D3 correction of Grimme. We find that the inclusion of the correction is essential to get meaningful results, but there is otherwise little to distinguish between the functionals. We also present coupled-cluster quality interaction curves for H2, NO2, H2O, and Ar interacting with large carbon flakes, acting as models for graphene surfaces, using novel absolutely localized molecular orbital based methods. These calculations demonstrate that the problems with empirically corrected DFT when investigating dispersion appear to compound as the system size increases, with important implications for future computational studies of molecule-surface interactions.
Collapse
|
48
|
Guo J, Zhou L, Zen A, Michaelides A, Wu X, Wang E, Xu L, Chen J. Hydration of NH_{4}^{+} in Water: Bifurcated Hydrogen Bonding Structures and Fast Rotational Dynamics. PHYSICAL REVIEW LETTERS 2020; 125:106001. [PMID: 32955332 DOI: 10.1103/physrevlett.125.106001] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 08/04/2020] [Indexed: 05/08/2023]
Abstract
Understanding the hydration and diffusion of ions in water at the molecular level is a topic of widespread importance. The ammonium ion (NH_{4}^{+}) is an exemplar system that has received attention for decades because of its complex hydration structure and relevance in industry. Here we report a study of the hydration and the rotational diffusion of NH_{4}^{+} in water using ab initio molecular dynamics simulations and quantum Monte Carlo calculations. We find that the hydration structure of NH_{4}^{+} features bifurcated hydrogen bonds, which leads to a rotational mechanism involving the simultaneous switching of a pair of bifurcated hydrogen bonds. The proposed hydration structure and rotational mechanism are supported by existing experimental measurements, and they also help to rationalize the measured fast rotation of NH_{4}^{+} in water. This study highlights how subtle changes in the electronic structure of hydrogen bonds impacts the hydration structure, which consequently affects the dynamics of ions and molecules in hydrogen bonded systems.
Collapse
Affiliation(s)
- Jianqing Guo
- International Center for Quantum Materials, Peking University, Beijing 100871, People's Republic of China
- School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Liying Zhou
- International Center for Quantum Materials, Peking University, Beijing 100871, People's Republic of China
- School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Andrea Zen
- Department of Physics and Astronomy, Thomas Young Centre and London Centre for Nanotechnology University College London, Gower Street, London WC1E 6BT, United Kingdom
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- Department of Physics and Astronomy, Thomas Young Centre and London Centre for Nanotechnology University College London, Gower Street, London WC1E 6BT, United Kingdom
- Max Planck Institute for Solid State Research, Stuttgart 70569, Germany
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Enge Wang
- International Center for Quantum Materials, Peking University, Beijing 100871, People's Republic of China
- School of Physics, Peking University, Beijing 100871, People's Republic of China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, People's Republic of China
- Songshan Lake Materials Lab, Institute of Physics, Chinese Academy of Sciences, Guangdong 523808, People's Republic of China
- School of Physics, Liaoning University, Shenyang 110136, People's Republic of China
| | - Limei Xu
- International Center for Quantum Materials, Peking University, Beijing 100871, People's Republic of China
- School of Physics, Peking University, Beijing 100871, People's Republic of China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, People's Republic of China
| | - Ji Chen
- School of Physics, Peking University, Beijing 100871, People's Republic of China
- Max Planck Institute for Solid State Research, Stuttgart 70569, Germany
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, People's Republic of China
| |
Collapse
|
49
|
Elliott JD, Troisi A, Carbone P. A QM/MD Coupling Method to Model the Ion-Induced Polarization of Graphene. J Chem Theory Comput 2020; 16:5253-5263. [PMID: 32644791 DOI: 10.1021/acs.jctc.0c00239] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We report a new Quantum Mechanical/Molecular Dynamics (QM/MD) simulation loop to model the coupling between the electron and atom dynamics in solid/liquid interfacial systems. The method can describe simultaneously both the quantum mechanical surface polarizability emerging from the proximity to the electrolyte and the electrolyte structure and dynamics. In the current setup, Density Functional Tight Binding calculations for the electronic structure calculations of the surface are coupled with classical molecular dynamics to simulate the electrolyte solution. The reduced computational cost of the QM part makes the coupling with a classical simulation engine computationally feasible and allows simulation of large systems for hundreds of nanoseconds. We tested the method by simulating both a noncharged graphene flake and a noncharged and charged infinite graphene sheet immersed in an NaCl electrolyte solution. We found that, when no bias is applied, ions preferentially remained in solution, and only cations are mildly attracted to the surface of the graphene. This preferential adsorption of cations vs anions seems to persist also when the surface is moderately charged and rules out any substantial ions/surface charge transfer.
Collapse
Affiliation(s)
- Joshua D Elliott
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Alessandro Troisi
- Department of Chemistry, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| | - Paola Carbone
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| |
Collapse
|
50
|
Tocci G, Bilichenko M, Joly L, Iannuzzi M. Ab initio nanofluidics: disentangling the role of the energy landscape and of density correlations on liquid/solid friction. NANOSCALE 2020; 12:10994-11000. [PMID: 32426791 DOI: 10.1039/d0nr02511a] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Despite relevance to water purification and renewable energy conversion membranes, the molecular mechanisms underlying water slip are poorly understood. We disentangle the static and dynamical origin of water slippage on graphene, hBN and MoS2 by means of large-scale ab initio molecular dynamics. Accounting for the role of the electronic structure of the interface is essential to determine that water slips five and eleven times faster on graphene compared to hBN and to MoS2, respectively. Intricate changes in the water energy landscape as well as in the density correlations of the fluid provide, respectively, the main static and dynamical origin of water slippage. Surprisingly, the timescales of the density correlations are the same on graphene and hBN, whereas they are longer on MoS2 and yield a 100% slowdown in the flow of water on this material. Our results pave the way for an in silico first principles design of materials with enhanced water slip, through the modification of properties connected not only to the structure, but also to the dynamics of the interface.
Collapse
Affiliation(s)
- Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland.
| | - Maria Bilichenko
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland.
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France and Institut Universitaire de France (IUF), France
| | - Marcella Iannuzzi
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland.
| |
Collapse
|