1
|
Luo S, Misra RP, Blankschtein D. Water Electric Field Induced Modulation of the Wetting of Hexagonal Boron Nitride: Insights from Multiscale Modeling of Many-Body Polarization. ACS NANO 2024; 18:1629-1646. [PMID: 38169482 DOI: 10.1021/acsnano.3c09811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Understanding the behavior of water contacting two-dimensional materials, such as hexagonal boron nitride (hBN), is important in practical applications, including seawater desalination and energy harvesting. Water, being a polar solvent, can strongly polarize the hBN surface via the electric fields that it generates. However, there is a lack of molecular-level understanding about the role of polarization effects at the hBN/water interface, including its effect on the wetting properties of water. In this study, we develop a theoretical framework that introduces an all-atomistic polarizable force field to accurately model the interactions of water molecules with hBN surfaces. The force field is then utilized to self-consistently describe the water-induced polarization of hBN using the classical Drude oscillator model, including predicting the hBN-water binding energies which are found to be in excellent agreement with diffusion Monte Carlo (DMC) predictions. By carrying out molecular dynamics (MD) simulations, we demonstrate that the polarizable force field yields a water contact angle on multilayered hBN which is in close agreement with the recent experimentally reported values. Conversely, an implicit modeling of the hBN-water polarization energy utilizing a Lennard-Jones (LJ) potential, a commonly utilized approximation in previous MD simulation studies, leads to a considerably lower water contact angle. This difference in the predicted contact angles is attributed to the significant energy-entropy compensation resulting from the incorporation of polarization effects at the hBN-water interface. Our work highlights the importance of self-consistently modeling the hBN-water polarization energy and offers insights into the wetting-related interfacial phenomena of water on polarizable materials.
Collapse
Affiliation(s)
- Shuang Luo
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Rahul Prasanna Misra
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Daniel Blankschtein
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
2
|
Novotný M, Dubecký M, Karlický F. Toward accurate modeling of structure and energetics of bulk hexagonal boron nitride. J Comput Chem 2024; 45:115-121. [PMID: 37737623 DOI: 10.1002/jcc.27222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Revised: 08/17/2023] [Accepted: 08/25/2023] [Indexed: 09/23/2023]
Abstract
Materials that exhibit both strong covalent and weak van der Waals interactions pose a considerable challenge to many computational methods, such as DFT. This makes assessing the accuracy of calculated properties, such as exfoliation energies in layered materials like hexagonal boron nitride (h-BN) problematic, when experimental data are not available. In this paper, we investigate the accuracy of equilibrium lattice constants and exfoliation energy calculation for various DFT-based computational approaches in bulk h-BN. We contrast these results with available experiments and reference fixed-node diffusion quantum Monte Carlo (QMC) results. From our reference QMC calculation, we obtained an exfoliation energy of - 33 ± 2 meV/atom (-0.38 ± 0.02 J/m2 ).
Collapse
Affiliation(s)
- Michal Novotný
- Department of Physics, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
| | - Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
- ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, Trnava, Slovakia
| | - František Karlický
- Department of Physics, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
| |
Collapse
|
3
|
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
|
4
|
Rothe K, Néel N, Bocquet ML, Kröger J. Tracking the Interaction between a CO-Functionalized Probe and Two Ag-Phthalocyanine Conformers by Local Vertical Force Spectroscopy. J Phys Chem A 2022; 126:6890-6897. [PMID: 36154143 DOI: 10.1021/acs.jpca.2c04760] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Intentionally terminating scanning probes with a single atom or molecule belongs to a rapidly growing field in the quantum chemistry and physics at surfaces. However, the detailed understanding of the coupling between the probe and adsorbate is in its infancy. Here, an atomic force microscopy probe functionalized with a single CO molecule is approached with picometer control to two conformational isomers of Ag-phthalocyanine adsorbed on Ag(111). The isomer with the central Ag atom pointing to CO exhibits a complex evolution of the distance-dependent interaction, while the conformer with Ag bonded to the metal surface gives rise to a Lennard-Jones behavior. By virtue of spatially resolved force spectroscopy and the comparison with results obtained from microscope probes terminated with a single Ag atom, the mutual coupling of the protruding O atom of the tip and the Ag atom of the phthalocyanine molecule is identified as the cause for the unconventional variation of the force. Simulations of the entire junction within density functional theory unveil the presence of ample relaxations in the case of one conformer, which represents a rationale for the peculiar vertical-distance evolution of the interaction. The simulations highlight the role of physisorption, chemisorption, and unexpected junction distortions at the verge of bond formation in the interpretation of the distance-dependent force between two molecules.
Collapse
Affiliation(s)
- Karl Rothe
- Institut für Physik, Technische Universität Ilmenau, D-98693 Ilmenau, Germany
| | - Nicolas Néel
- Institut für Physik, Technische Universität Ilmenau, D-98693 Ilmenau, Germany
| | - Marie-Laure Bocquet
- Laboratoire de Physique de l'École Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris-Diderot, Sorbonne Paris Cité, F-75005 Paris, France
| | - Jörg Kröger
- Institut für Physik, Technische Universität Ilmenau, D-98693 Ilmenau, Germany
| |
Collapse
|
5
|
Kumar Verma A, Govind Rajan A. Surface Roughness Explains the Observed Water Contact Angle and Slip Length on 2D Hexagonal Boron Nitride. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:9210-9220. [PMID: 35866875 DOI: 10.1021/acs.langmuir.2c00972] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Hexagonal boron nitride (hBN) is a two-dimensional (2D) material that is currently being explored in a number of applications, such as atomically thin coatings, water desalination, and biological sensors. In many of these applications, the hBN surface comes into intimate contact with water. In this work, we investigate the wetting and frictional behavior of realistic 2D hBN surfaces with atomic-scale defects and roughness. We combine density functional theory calculations of the charge distribution inside hBN with free energy calculations using molecular dynamics simulations of the hBN-water interface. We find that the presence of surface roughness, but not that of vacancy defects, leads to remarkable agreement with the experimentally observed water contact angle of 66° on freshly synthesized, uncontaminated hBN. Not only that, the inclusion of surface roughness predicts with exceptional accuracy the experimental water slip length of ∼1 nm on hBN. Our results underscore the importance of considering realistic models of 2D materials with surface roughness while modeling nanomaterial-water interfaces in molecular simulations.
Collapse
Affiliation(s)
- Ashutosh Kumar Verma
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| | - Ananth Govind Rajan
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, Karnataka 560012, India
| |
Collapse
|
6
|
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
|
7
|
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
|
8
|
Bai C, Yang Z, Zhang J, Zhang B, Yu Y, Zhang J. Friction Behavior and Structural Evolution of Hexagonal Boron Nitride: A Relation to Environmental Molecules Containing -OH Functional Group. ACS APPLIED MATERIALS & INTERFACES 2022; 14:19043-19055. [PMID: 35416641 DOI: 10.1021/acsami.2c02450] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Sliding contact experiments and first-principles calculations were performed to elucidate the roles of environmental molecules containing -OH functional groups on the friction behavior and structural evolution of hexagonal boron nitride (h-BN). A significant decrease in the friction coefficient (COF) is established by the physisorption and dissociative adsorption of molecules containing -OH functional groups on h-BN, compared with that in a H2 or N2 atmosphere. A key finding is the existence of two friction mechanisms to reconstruct the sliding interface for h-BN crystallites in humid air and carbon contaminant (CH3OH and C2H5OH) atmospheres, which is verified by the friction behavior and morphologies of the wear track. There is a running-in period in the friction process to induce the formation of defects in h-BN in humid air, which facilitates dissociative adsorption of water molecules on h-BN. The formation of nanostructured water at defect sites will promote lamellar slip of h-BN crystal materials for friction reduction. In carbon contaminant environments, both molecules exhibit strong adsorption on the h-BN surface regardless of the presence of defects, thereby weakening the structural damage rate and enhancing the bearing capacity. C2H5OH molecules are more likely to dissociate and bind onto defect sites, endowing h-BN with high in-plane stress to form a coiled structure. h-BN in the annular or tubular form would exhibit a self-protective effect, facilitate incommensurate contact, and reduce the contact area to enhance lubrication. Our work may establish the fundamental basis for future applications of h-BN in new energy vehicles.
Collapse
Affiliation(s)
- Changning Bai
- Key Laboratory of Science and Technology on Wear and Protection of Materials, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou 730000, China
- Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zaixiu Yang
- Key Laboratory of Science and Technology on Wear and Protection of Materials, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou 730000, China
| | - Jing Zhang
- International School for Optoelectronic Engineering, Qilu University of Technology (Shandong Academy of Sciences), Jinan 250353, China
| | - Bin Zhang
- Key Laboratory of Science and Technology on Wear and Protection of Materials, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou 730000, China
| | - Yuanlie Yu
- Key Laboratory of Science and Technology on Wear and Protection of Materials, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou 730000, China
- Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Junyan Zhang
- Key Laboratory of Science and Technology on Wear and Protection of Materials, Lanzhou Institute of Chemical Physics, Chinese Academy of Sciences, Lanzhou 730000, China
- Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
9
|
Tsuji Y, Yoshizawa K. Competition between Hydrogen Bonding and Dispersion Force in Water Adsorption and Epoxy Adhesion to Boron Nitride: From the Flat to the Curved. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:11351-11364. [PMID: 34519515 DOI: 10.1021/acs.langmuir.1c01935] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Hexagonal boron nitride (h-BN) is a material with excellent thermal conductivity and electrical insulation, used as an additive to various matrices. To increase the affinity of h-BN to them, hydrogen bonds should be formed at the interface. In reality, however, they are not formed; the N atoms are not capable of accepting hydrogen bonds due to the delocalization of their lone pair electrons over the B-N π bonds. To make it form hydrogen bonds, one may need to break the planarity of h-BN so that the orbital overlap in the B-N π bonds can be reduced. This idea is verified with first-principles calculations on the adsorption of a water molecule on hypothetical h-BN surfaces, the planarity of which is broken. One can do it in silico but not in vitro. BN nanotubes (BNNTs) are considered as a more realistic BN surface with nonplanarity. The hydrogen bond is shown to become stronger as the curvature of the tube increases. On the contrary, the strength of the dispersion force acting at the interface becomes weaker. In water adsorption, these two interactions are in competition with each other. However, in epoxy adhesion, the interaction due to dispersion forces is overwhelmingly stronger than that due to hydrogen bonding. The smaller the curvature of the surface, the smaller the distance between more atoms at the interface; thus, the interaction due to dispersion forces maximized.
Collapse
Affiliation(s)
- Yuta Tsuji
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan
| | - Kazunari Yoshizawa
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan
| |
Collapse
|
10
|
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
|
11
|
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
|
12
|
Needs RJ, Towler MD, Drummond ND, López Ríos P, Trail JR. Variational and diffusion quantum Monte Carlo calculations with the CASINO code. J Chem Phys 2020; 152:154106. [DOI: 10.1063/1.5144288] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- R. J. Needs
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - M. D. Towler
- University College London, London WC1E 6BT, United Kingdom
| | - N. D. Drummond
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
| | - P. López Ríos
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - J. R. Trail
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
13
|
Klimeš J, Tew DP. Efficient and accurate description of adsorption in zeolites. J Chem Phys 2019; 151:234108. [PMID: 31864262 DOI: 10.1063/1.5123425] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.
Collapse
Affiliation(s)
- Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - David P Tew
- Max-Planck-Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
14
|
Budi A, Walsh TR. A Bespoke Force Field To Describe Biomolecule Adsorption at the Aqueous Boron Nitride Interface. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:16234-16243. [PMID: 31714785 DOI: 10.1021/acs.langmuir.9b03121] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Reliable manipulation of the interface between 2D nanomaterials and biomolecules represents a current frontier in nanoscience. The ability to resolve the molecular-level structures of these biointerfaces would provide a fundamental data set that is needed to enable systematic and knowledge-based progress in this area. These structures are challenging to obtain via experiment alone, and molecular simulations offer a complementary approach to address this problem. Compared with graphene, the interface between hexagonal boron nitride (h-BN) and biomolecules is relatively understudied at present. While several force fields are currently available for modeling the h-BN/water interface, there is a lack of a suitable force field that can describe the interactions between h-BN, liquid water, and biomolecules. Here, we use density functional theory calculations to create a force field, BoNi-CHARMM, to describe biomolecular interactions at the aqueous h-BN interface. Verifying our force field presents an additional challenge, given the scarcity of available experimental data for these interfaces. We test our force field against experimental evidence regarding the water/surface contact angle and confirm that the force field provides experimentally consistent values. We also present preliminary data regarding predictions of the free energy of adsorption of a selection of amino acids at the aqueous h-BN interface, revealing arginine and tryptophan to be among the strongest binders. This force field provides an opportunity to initiate a systematic progression in our current understanding of how to capture the intermolecular interactions at the h-BN biointerface.
Collapse
Affiliation(s)
- Akin Budi
- Institute for Frontier Materials , Deakin University , 75 Pigdon's Rd. , Geelong , Victoria 3216 , Australia
| | - Tiffany R Walsh
- Institute for Frontier Materials , Deakin University , 75 Pigdon's Rd. , Geelong , Victoria 3216 , Australia
| |
Collapse
|
15
|
Zen A, Brandenburg JG, Michaelides A, Alfè D. A new scheme for fixed node diffusion quantum Monte Carlo with pseudopotentials: Improving reproducibility and reducing the trial-wave-function bias. J Chem Phys 2019; 151:134105. [DOI: 10.1063/1.5119729] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrea Zen
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
| | - Jan Gerit Brandenburg
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, Im Neuenheimer Feld 205A, 69120 Heidelberg, Germany
| | - Angelos Michaelides
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Dario Alfè
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| |
Collapse
|
16
|
Dubecký M, Jurečka P, Mitas L, Ditte M, Fanta R. Toward Accurate Hydrogen Bonds by Scalable Quantum Monte Carlo. J Chem Theory Comput 2019; 15:3552-3557. [DOI: 10.1021/acs.jctc.9b00096] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Matúš Dubecký
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
- ATRI, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovakia
| | - Petr Jurečka
- Department of Physical Chemistry, Palackỳ University Olomouc, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Matej Ditte
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Roman Fanta
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| |
Collapse
|
17
|
Govind Rajan A, Strano MS, Blankschtein D. Liquids with Lower Wettability Can Exhibit Higher Friction on Hexagonal Boron Nitride: The Intriguing Role of Solid-Liquid Electrostatic Interactions. NANO LETTERS 2019; 19:1539-1551. [PMID: 30694070 DOI: 10.1021/acs.nanolett.8b04335] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
We investigate the wetting and frictional behavior of polar (water and ethylene glycol) and nonpolar (diiodomethane) liquids on the basal plane of hexagonal boron nitride (hBN) using molecular dynamics simulations. Our results for the wettability of water on the hBN basal plane (contact angle 81°) are in qualitative agreement with the experimentally deduced mild hydrophilicity of the hBN basal plane (contact angle 66°). We find that water exhibits the lowest wettability, as quantified by the highest contact angle, but the highest friction coefficient of (1.9 ± 0.4) × 105 N-s/m3 on the hBN basal plane among the three liquids considered. This intriguing finding is explained in terms of the competition between dispersion and electrostatic interactions operating between the hBN basal plane and the three liquids. We find that electrostatic interactions do not affect the wetting behavior appreciably, as quantified by a less than 3° change in the respective contact angles of the three liquids considered. On the other hand, electrostatic interactions are found to increase the friction coefficients of the three liquids in contact with hBN to different extents, indicating that despite the increased friction of water on hBN, relative to that on graphene, nonpolar liquids may exhibit similar friction coefficients on hBN and graphene. Our findings reveal that the increase in the friction coefficient, upon incorporation of solid-liquid electrostatic interactions, is brought about by a greater increase in the solid-liquid mean-squared total lateral force, as compared to a smaller reduction in the decorrelation time of the solid-liquid force.
Collapse
Affiliation(s)
- Ananth Govind Rajan
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Michael S Strano
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Daniel Blankschtein
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
18
|
First-principles modeling of water permeation through periodically porous graphene derivatives. J Colloid Interface Sci 2019; 538:367-376. [DOI: 10.1016/j.jcis.2018.11.106] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Revised: 11/27/2018] [Accepted: 11/28/2018] [Indexed: 01/24/2023]
|
19
|
Al-Hamdani YS, Tkatchenko A. Understanding non-covalent interactions in larger molecular complexes from first principles. J Chem Phys 2019; 150:010901. [PMID: 30621423 PMCID: PMC6910608 DOI: 10.1063/1.5075487] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Accepted: 12/05/2018] [Indexed: 01/02/2023] Open
Abstract
Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
20
|
Schwarz M, Duncan DA, Garnica M, Ducke J, Deimel PS, Thakur PK, Lee TL, Allegretti F, Auwärter W. Quantitative determination of a model organic/insulator/metal interface structure. NANOSCALE 2018; 10:21971-21977. [PMID: 30444513 PMCID: PMC6289171 DOI: 10.1039/c8nr06387g] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/07/2018] [Accepted: 09/28/2018] [Indexed: 05/22/2023]
Abstract
By combining X-ray photoelectron spectroscopy, X-ray standing waves and scanning tunneling microscopy, we investigate the geometric and electronic structure of a prototypical organic/insulator/metal interface, namely cobalt porphine on monolayer hexagonal boron nitride (h-BN) on Cu(111). Specifically, we determine the adsorption height of the organic molecule and show that the original planar molecular conformation is preserved in contrast to the adsorption on Cu(111). In addition, we highlight the electronic decoupling provided by the h-BN spacer layer and find that the h-BN-metal separation is not significantly modified by the molecular adsorption. Finally, we find indication of a temperature dependence of the adsorption height, which might be a signature of strongly-anisotropic thermal vibrations of the weakly bonded molecules.
Collapse
Affiliation(s)
- Martin Schwarz
- Physics Department
, Technical University of Munich
,
85748 Garching
, Germany
.
;
| | - David A. Duncan
- Diamond Light Source
, Harwell Science and Innovation Campus
,
Didcot OX11 0DE
, UK
| | - Manuela Garnica
- Physics Department
, Technical University of Munich
,
85748 Garching
, Germany
.
;
| | - Jacob Ducke
- Physics Department
, Technical University of Munich
,
85748 Garching
, Germany
.
;
| | - Peter S. Deimel
- Physics Department
, Technical University of Munich
,
85748 Garching
, Germany
.
;
| | - Pardeep K. Thakur
- Diamond Light Source
, Harwell Science and Innovation Campus
,
Didcot OX11 0DE
, UK
| | - Tien-Lin Lee
- Diamond Light Source
, Harwell Science and Innovation Campus
,
Didcot OX11 0DE
, UK
| | - Francesco Allegretti
- Physics Department
, Technical University of Munich
,
85748 Garching
, Germany
.
;
| | - Willi Auwärter
- Physics Department
, Technical University of Munich
,
85748 Garching
, Germany
.
;
| |
Collapse
|
21
|
Liu J, Pantelides ST. Electrowetting on 2D dielectrics: a quantum molecular dynamics investigation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:375001. [PMID: 30079895 DOI: 10.1088/1361-648x/aad838] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Electrowetting on dielectrics (EWOD) is widely used to manipulate the spreading of a conductive liquid on a dielectric surface by applying an electric field. 2D hydrophobic dielectrics are promising candidates for EWOD applications. In this study, extensive quantum molecular dynamics (MD) simulations are performed to investigate the electrowetting behavior of salty water on hexagonal boron nitride (h-BN) monolayer. The proximal adsorption of salt ions and the associated realignment of the dipole moments of interfacial water with the applied electric field are found to be the physical origin of the electrowetting behavior. At low salt concentration and low electric fields, the proximal adsorption and the realignment follow the applied electric field, and the cosine of the water contact angle (WCA) follows a quadratic dependence on the applied electric field. At high salt concentration and high electric fields, the proximal adsorption saturates, which restricts further realignment and causes a saturation of the WCA. This case study provides physical insights into the much debated mechanism that underlies the contact angle saturation (CAS) found in macroscopic electrowetting phenomena and also provides an avenue for further studies of electrowetting at the atomic scale.
Collapse
Affiliation(s)
- Jian Liu
- Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, United States of America
| | | |
Collapse
|
22
|
Al-Hamdani YS, Michaelides A, von Lilienfeld OA. Exploring dissociative water adsorption on isoelectronically BN doped graphene using alchemical derivatives. J Chem Phys 2018; 147:164113. [PMID: 29096500 DOI: 10.1063/1.4986314] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The design and production of novel 2-dimensional materials have seen great progress in the last decade, prompting further exploration of the chemistry of such materials. Doping and hydrogenating graphene are an experimentally realised method of changing its surface chemistry, but there is still a great deal to be understood on how doping impacts on the adsorption of molecules. Developing this understanding is key to unlocking the potential applications of these materials. High throughput screening methods can provide particularly effective ways to explore vast chemical compositions of materials. Here, alchemical derivatives are used as a method to screen the dissociative adsorption energy of water molecules on various BN doped topologies of hydrogenated graphene. The predictions from alchemical derivatives are assessed by comparison to density functional theory. This screening method is found to predict dissociative adsorption energies that span a range of more than 2 eV, with a mean absolute error <0.1 eV. In addition, we show that the quality of such predictions can be readily assessed by examination of the Kohn-Sham highest occupied molecular orbital in the initial states. In this way, the root mean square error in the dissociative adsorption energies of water is reduced by almost an order of magnitude (down to ∼0.02 eV) after filtering out poor predictions. The findings point the way towards a reliable use of first order alchemical derivatives for efficient screening procedures.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
23
|
Govind Rajan A, Strano MS, Blankschtein D. Ab Initio Molecular Dynamics and Lattice Dynamics-Based Force Field for Modeling Hexagonal Boron Nitride in Mechanical and Interfacial Applications. J Phys Chem Lett 2018. [PMID: 29528646 DOI: 10.1021/acs.jpclett.7b03443] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Hexagonal boron nitride (hBN) is an up-and-coming two-dimensional material, with applications in electronic devices, tribology, and separation membranes. Herein, we utilize density-functional-theory-based ab initio molecular dynamics (MD) simulations and lattice dynamics calculations to develop a classical force field (FF) for modeling hBN. The FF predicts the crystal structure, elastic constants, and phonon dispersion relation of hBN with good accuracy and exhibits remarkable agreement with the interlayer binding energy predicted by random phase approximation calculations. We demonstrate the importance of including Coulombic interactions but excluding 1-4 intrasheet interactions to obtain the correct phonon dispersion relation. We find that improper dihedrals do not modify the bulk mechanical properties and the extent of thermal vibrations in hBN, although they impact its flexural rigidity. Combining the FF with the accurate TIP4P/Ice water model yields excellent agreement with interaction energies predicted by quantum Monte Carlo calculations. Our FF should enable an accurate description of hBN interfaces in classical MD simulations.
Collapse
Affiliation(s)
- Ananth Govind Rajan
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Michael S Strano
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Daniel Blankschtein
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
24
|
Kayal A, Chandra A. Orientational order and dynamics of interfacial water near a hexagonal boron-nitride sheet: An ab initio molecular dynamics study. J Chem Phys 2017; 147:164704. [DOI: 10.1063/1.4991594] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Abhijit Kayal
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
25
|
Al-Hamdani YS, Rossi M, Alfè D, Tsatsoulis T, Ramberger B, Brandenburg JG, Zen A, Kresse G, Grüneis A, Tkatchenko A, Michaelides A. Properties of the water to boron nitride interaction: From zero to two dimensions with benchmark accuracy. J Chem Phys 2017; 147:044710. [DOI: 10.1063/1.4985878] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yasmine S. Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Mariana Rossi
- Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, 14195 Berlin, Germany
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Theodoros Tsatsoulis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Wien, Austria
| | - Jan Gerit Brandenburg
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AH, United Kingdom
| | - Andrea Zen
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Wien, Austria
| | - Andreas Grüneis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
26
|
Dubecký M. Noncovalent Interactions by Fixed-Node Diffusion Monte Carlo: Convergence of Nodes and Energy Differences vs Gaussian Basis-Set Size. J Chem Theory Comput 2017; 13:3626-3635. [DOI: 10.1021/acs.jctc.7b00537] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701
03 Ostrava, Czech Republic
- ATRI, Faculty of Materials
Science and Technology, Slovak University of Technology, Paulínska
16, 917 24 Trnava, Slovakia
| |
Collapse
|
27
|
Vlaisavljevich B, Huck J, Hulvey Z, Lee K, Mason JA, Neaton JB, Long JR, Brown CM, Alfè D, Michaelides A, Smit B. Performance of van der Waals Corrected Functionals for Guest Adsorption in the M 2(dobdc) Metal-Organic Frameworks. J Phys Chem A 2017; 121:4139-4151. [PMID: 28436661 PMCID: PMC5529028 DOI: 10.1021/acs.jpca.7b00076] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Revised: 04/21/2017] [Indexed: 11/29/2022]
Abstract
Small-molecule binding in metal-organic frameworks (MOFs) can be accurately studied both experimentally and computationally, provided the proper tools are employed. Herein, we compare and contrast properties associated with guest binding by means of density functional theory (DFT) calculations using nine different functionals for the M2(dobdc) (dobdc4- = 2,5-dioxido,1,4-benzenedicarboxylate) series, where M = Mg, Mn, Fe, Co, Ni, Cu, and Zn. Additionally, we perform Quantum Monte Carlo (QMC) calculations for one system to determine if this method can be used to assess the performance of DFT. We also make comparisons with previously published experimental results for carbon dioxide and water and present new methane neutron powder diffraction (NPD) data for further comparison. All of the functionals are able to predict the experimental variation in the binding energy from one metal to the next; however, the interpretation of the performance of the functionals depends on which value is taken as the reference. On the one hand, if we compare against experimental values, we would conclude that the optB86b-vdW and optB88-vdW functionals systematically overestimate the binding strength, while the second generation of van der Waals (vdW) nonlocal functionals (vdw-DF2 and rev-vdW-DF2) correct for this providing a good description of binding energies. On the other hand, if the QMC calculation is taken as the reference then all of the nonlocal functionals yield results that fall just outside the error of the higher-level calculation. The empirically corrected vdW functionals are in reasonable agreement with experimental heat of adsorptions but under bind when compared with QMC, while Perdew-Burke-Ernzerhof fails by more than 20 kJ/mol regardless of which reference is employed. All of the functionals, with the exception of vdW-DF2, predict reasonable framework and guest binding geometries when compared with NPD measurements. The newest of the functionals considered, rev-vdW-DF2, should be used in place of vdW-DF2, as it yields improved bond distances with similar quality binding energies.
Collapse
Affiliation(s)
| | | | - Zeric Hulvey
- Center for Neutron Research, National Institute of Standards and Technology , Gaithersburg, Maryland 20899, United States
- Department of Materials Science and Engineering, University of Maryland , College Park, Maryland 20742 United States
| | | | | | - Jeffrey B Neaton
- Kavli Energy Nanosciences Institute at Berkeley , Berkeley, California 94720, United States
| | | | - Craig M Brown
- Center for Neutron Research, National Institute of Standards and Technology , Gaithersburg, Maryland 20899, United States
- Department of Chemical and Biomolecular Engineering, University of Delaware , Newark, Delaware 19716, United States
| | - Dario Alfè
- Department of Earth Sciences, Thomas Young Centre and London Centre for Nanotechnology, University College London , Gower Street, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology and Department of Physics and Astronomy, Thomas Young Centre, University College London , Gower Street, London WC1E 6BT, United Kingdom
| | - Berend Smit
- Institut des Sciences et Ingénierie Chimiques, Valais, Ecole Polytechnique Fédérale de Lausanne , Rue de l'Industrie 17, CH-1951 Sion, Switzerland
| |
Collapse
|
28
|
Tsatsoulis T, Hummel F, Usvyat D, Schütz M, Booth GH, Binnie SS, Gillan MJ, Alfè D, Michaelides A, Grüneis A. A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface. J Chem Phys 2017; 146:204108. [PMID: 28571392 PMCID: PMC5446292 DOI: 10.1063/1.4984048] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2017] [Accepted: 05/11/2017] [Indexed: 12/03/2022] Open
Abstract
We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.
Collapse
Affiliation(s)
- Theodoros Tsatsoulis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Felix Hummel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| | - Denis Usvyat
- Institut für Chemie, Humboldt-Universität zu Berlin, Brook-Taylor-St. 2, D-12489 Berlin, Germany
| | - Martin Schütz
- Institut für Chemie, Humboldt-Universität zu Berlin, Brook-Taylor-St. 2, D-12489 Berlin, Germany
| | - George H Booth
- Department of Physics, King's College London, Strand, London WC2R 2LS, United Kingdom
| | - Simon S Binnie
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Michael J Gillan
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Dario Alfè
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, University College London, Gordon St., London WC1H 0AH, United Kingdom
| | - Andreas Grüneis
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, D-70569 Stuttgart, Germany
| |
Collapse
|
29
|
Hummel F, Tsatsoulis T, Grüneis A. Low rank factorization of the Coulomb integrals for periodic coupled cluster theory. J Chem Phys 2017; 146:124105. [DOI: 10.1063/1.4977994] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Ramberger B, Schäfer T, Kresse G. Analytic Interatomic Forces in the Random Phase Approximation. PHYSICAL REVIEW LETTERS 2017; 118:106403. [PMID: 28339256 DOI: 10.1103/physrevlett.118.106403] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/04/2016] [Indexed: 06/06/2023]
Abstract
We discuss that in the random phase approximation (RPA) the first derivative of the energy with respect to the Green's function is the self-energy in the GW approximation. This relationship allows us to derive compact equations for the RPA interatomic forces. We also show that position dependent overlap operators are elegantly incorporated in the present framework. The RPA force equations have been implemented in the projector augmented wave formalism, and we present illustrative applications, including ab initio molecular dynamics simulations, the calculation of phonon dispersion relations for diamond and graphite, as well as structural relaxations for water on boron nitride. The present derivation establishes a concise framework for forces within perturbative approaches and is also applicable to more involved approximations for the correlation energy.
Collapse
Affiliation(s)
- Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Vienna, Austria
| | - Tobias Schäfer
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Vienna, Austria
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090 Vienna, Austria
| |
Collapse
|
31
|
Al-Hamdani YS, Alfè D, Michaelides A. How strongly do hydrogen and water molecules stick to carbon nanomaterials? J Chem Phys 2017. [DOI: 10.1063/1.4977180] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yasmine S. Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
32
|
Dubecký M. Bias cancellation in one-determinant fixed-node diffusion Monte Carlo: Insights from fermionic occupation numbers. Phys Rev E 2017; 95:033308. [PMID: 28415179 DOI: 10.1103/physreve.95.033308] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Indexed: 06/07/2023]
Abstract
The accuracy of the fixed-node diffusion Monte Carlo (FNDMC) depends on the node location of the supplied trial state Ψ_{T}. The practical FNDMC approaches available for large systems rely on compact yet effective Ψ_{T}, most often containing an explicitly correlated single Slater determinant (SD). However, SD nodes may be better suited to one system than to another, which may possibly lead to inaccurate FNDMC energy differences. It remains a challenge how to estimate nonequivalence or appropriateness of SDs. Here we use the differences of a measure based on the Euclidean distance between the natural orbital occupation number (NOON) vector of the SD and the exact solution in the NOON vector space, which can be viewed as a measure of SD nonequivalence and as a qualitative measure of the expected degree of nondynamic-correlation-related bias in FNDMC energy differences. This is explored on a set of small noncovalent complexes and covalent bond breaking of Si_{2} vs N_{2}. It turns out that NOON-based measures well reflect the magnitude and sign of the bias present in the data available, thus providing insights into the nature of bias cancellation in SD FNDMC energy differences.
Collapse
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic and ATRI, Faculty of Materials Science and Technology, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovakia
| |
Collapse
|
33
|
Competition of van der Waals and chemical forces on gold–sulfur surfaces and nanoparticles. Nat Rev Chem 2017. [DOI: 10.1038/s41570-017-0017] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
34
|
Wu Y, Wagner LK, Aluru NR. Hexagonal boron nitride and water interaction parameters. J Chem Phys 2017; 144:164118. [PMID: 27131542 DOI: 10.1063/1.4947094] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
The study of hexagonal boron nitride (hBN) in microfluidic and nanofluidic applications at the atomic level requires accurate force field parameters to describe the water-hBN interaction. In this work, we begin with benchmark quality first principles quantum Monte Carlo calculations on the interaction energy between water and hBN, which are used to validate random phase approximation (RPA) calculations. We then proceed with RPA to derive force field parameters, which are used to simulate water contact angle on bulk hBN, attaining a value within the experimental uncertainties. This paper demonstrates that end-to-end multiscale modeling, starting at detailed many-body quantum mechanics and ending with macroscopic properties, with the approximations controlled along the way, is feasible for these systems.
Collapse
Affiliation(s)
- Yanbin Wu
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Lucas K Wagner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801-3080, USA
| | - Narayana R Aluru
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
35
|
Zen A, Roch LM, Cox SJ, Hu XL, Sorella S, Alfè D, Michaelides A. Toward Accurate Adsorption Energetics on Clay Surfaces. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2016; 120:26402-26413. [PMID: 27917256 PMCID: PMC5126707 DOI: 10.1021/acs.jpcc.6b09559] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Revised: 10/31/2016] [Indexed: 05/31/2023]
Abstract
Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.
Collapse
Affiliation(s)
- Andrea Zen
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Loïc M. Roch
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Stephen J. Cox
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Xiao Liang Hu
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Sandro Sorella
- SISSA−International
School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- INFM
Democritos National Simulation Center, 34151 Trieste, Italy
| | - Dario Alfè
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United
Kingdom
| | - Angelos Michaelides
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
36
|
Al-Hamdani YS, Alfè D, von Lilienfeld OA, Michaelides A. Tuning dissociation using isoelectronically doped graphene and hexagonal boron nitride: Water and other small molecules. J Chem Phys 2016; 144:154706. [DOI: 10.1063/1.4945783] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yasmine S. Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Dario Alfè
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - O. Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17–19 Gordon Street, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
37
|
Affiliation(s)
- Matúš Dubecký
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
38
|
Joly L, Tocci G, Merabia S, Michaelides A. Strong Coupling between Nanofluidic Transport and Interfacial Chemistry: How Defect Reactivity Controls Liquid-Solid Friction through Hydrogen Bonding. J Phys Chem Lett 2016; 7:1381-1386. [PMID: 27012818 DOI: 10.1021/acs.jpclett.6b00280] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Defects are inevitably present in nanofluidic systems, yet the role they play in nanofluidic transport remains poorly understood. Here, we report ab initio molecular dynamics (AIMD) simulations of the friction of liquid water on defective graphene and boron nitride sheets. We show that water dissociates at certain defects and that these "reactive" defects lead to much larger friction than the "nonreactive" defects at which water molecules remain intact. Furthermore, we find that friction is extremely sensitive to the chemical structure of reactive defects and to the number of hydrogen bonds they can partake in with the liquid. Finally, we discuss how the insight obtained from AIMD can be used to quantify the influence of defects on friction in nanofluidic devices for water treatment and sustainable energy harvesting. Overall, we provide new insight into the role of interfacial chemistry on nanofluidic transport in real, defective systems.
Collapse
Affiliation(s)
- Laurent Joly
- Institut Lumière Matière, UMR5306 Université Lyon 1-CNRS, Université de Lyon, 69622 Villeurbanne, France
| | - Gabriele Tocci
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Physics and Astronomy, University College London , London WC1H 0AJ, United Kingdom
| | - Samy Merabia
- Institut Lumière Matière, UMR5306 Université Lyon 1-CNRS, Université de Lyon, 69622 Villeurbanne, France
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Physics and Astronomy, University College London , London WC1H 0AJ, United Kingdom
| |
Collapse
|
39
|
Abstract
David Craig (1919–2015) left us with a lasting legacy concerning basic understanding of chemical spectroscopy and bonding. This is expressed in terms of some of the recent achievements of my own research career, with a focus on integration of Craig’s theories with those of Noel Hush to solve fundamental problems in photosynthesis, molecular electronics (particularly in regard to the molecules synthesized by Maxwell Crossley), and self-assembled monolayer structure and function. Reviewed in particular is the relation of Craig’s legacy to: the 50-year struggle to assign the visible absorption spectrum of arguably the world’s most significant chromophore, chlorophyll; general theories for chemical bonding and structure extending Hush’s adiabatic theory of electron-transfer processes; inelastic electron-tunnelling spectroscopy (IETS); chemical quantum entanglement and the Penrose–Hameroff model for quantum consciousness; synthetic design strategies for NMR quantum computing; Gibbs free-energy measurements and calculations for formation and polymorphism of organic self-assembled monolayers on graphite surfaces from organic solution; and understanding the basic chemical processes involved in the formation of gold surfaces and nanoparticles protected by sulfur-bound ligands, ligands whose form is that of Au0-thiyl rather than its commonly believed AuI-thiolate tautomer.
Collapse
|
40
|
Fitzner M, Sosso GC, Cox SJ, Michaelides A. The Many Faces of Heterogeneous Ice Nucleation: Interplay Between Surface Morphology and Hydrophobicity. J Am Chem Soc 2015; 137:13658-69. [PMID: 26434775 DOI: 10.1021/jacs.5b08748] [Citation(s) in RCA: 105] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
What makes a material a good ice nucleating agent? Despite the importance of heterogeneous ice nucleation to a variety of fields, from cloud science to microbiology, major gaps in our understanding of this ubiquitous process still prevent us from answering this question. In this work, we have examined the ability of generic crystalline substrates to promote ice nucleation as a function of the hydrophobicity and the morphology of the surface. Nucleation rates have been obtained by brute-force molecular dynamics simulations of coarse-grained water on top of different surfaces of a model fcc crystal, varying the water-surface interaction and the surface lattice parameter. It turns out that the lattice mismatch of the surface with respect to ice, customarily regarded as the most important requirement for a good ice nucleating agent, is at most desirable but not a requirement. On the other hand, the balance between the morphology of the surface and its hydrophobicity can significantly alter the ice nucleation rate and can also lead to the formation of up to three different faces of ice on the same substrate. We have pinpointed three circumstances where heterogeneous ice nucleation can be promoted by the crystalline surface: (i) the formation of a water overlayer that acts as an in-plane template; (ii) the emergence of a contact layer buckled in an ice-like manner; and (iii) nucleation on compact surfaces with very high interaction strength. We hope that this extensive systematic study will foster future experimental work aimed at testing the physiochemical understanding presented herein.
Collapse
Affiliation(s)
- Martin Fitzner
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Gabriele C Sosso
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Stephen J Cox
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, Department of Chemistry and Thomas Young Centre, University College London , 17-19 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
41
|
Michaelides A, Martinez TJ, Alavi A, Kresse G, Manby FR. Preface: Special Topic Section on Advanced Electronic Structure Methods for Solids and Surfaces. J Chem Phys 2015; 143:102601. [DOI: 10.1063/1.4930182] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
42
|
Wu Y, Wagner LK, Aluru NR. The interaction between hexagonal boron nitride and water from first principles. J Chem Phys 2015; 142:234702. [DOI: 10.1063/1.4922491] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Yanbin Wu
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Lucas K. Wagner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801-3080, USA
| | - Narayana R. Aluru
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|