How strongly do hydrogen and water molecules stick to carbon nanomaterials?

Unraveling H2 chemisorption and physisorption on metal decorated graphene using quantum Monte Carlo

Molecular hydrogen has the potential to significantly reduce the use of carbon dioxide emitting energy processes. However, hydrogen gas storage is a major bottleneck for its large-scale use as current storage methods are energy intensive. Among different storage methods, physisorbing molecular hydrogen at ambient pressure and temperatures is a promising alternative-particularly in light of the advancements in tunable lightweight nanomaterials and high throughput screening methods. Nonetheless, understanding hydrogen adsorption in well-defined nanomaterials remains experimentally challenging and reference information is scarce despite the proliferation of works predicting hydrogen adsorption. We focus on Li, Na, Ca, and K, decorated graphene sheets as substrates for molecular hydrogen adsorption, and compute the most accurate adsorption energies available to date using quantum diffusion Monte Carlo (DMC). Building on our previous insights at the density functional theory (DFT) level, we find that a weak covalent chemisorption of molecular hydrogen, known as Kubas interaction, is feasible on Ca decorated graphene according to DMC, in agreement with DFT. This finding is in contrast to previous DMC predictions of the 4H2/Ca+ gas cluster (without graphene) where chemisorption is not favored. However, we find that the adsorption energy of hydrogen on metal decorated graphene according to a widely used DFT method is not fully consistent with DMC. The reference adsorption energies reported herein can be used to find better work-horse methods for application in large-scale modeling of hydrogen adsorption. Furthermore, the implications of this work affect strategies for finding suitable hydrogen storage materials and high-throughput methods.

Ab Initio Molecular Dynamics Investigation of the Solvation States of Hydrated Ions in Confined Water

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.

Tuning humidity sensing properties via grafting fluorine and nitrogen-containing species on single-walled carbon nanotubes

The effect of humidity on the electrical conductivity of single-walled carbon nanotube (SWCNT) films depends on both the conductivity of individual nanotubes and the electrical contacts between them. Here, we study these factors by comparing the sensor response of nanotubes with fluorine- and nitrogen-containing groups attached to the sidewalls. Experiments carried out in a wide range of relative humidity (RH) at room and elevated temperatures showed that the conductivity of non-functionalized SWCNTs and contacts between them decreases upon the adsorption of water molecules. Covalent fluorination reduces the conductivity of SWCNTs and significantly increases the sensitivity of the film to low concentrations of water vapor. The response at high RH decreases due to the large number of water molecules adsorbed on the conductive regions of the nanotubes. As a result of substitutional reactions of fluorinated SWCNTs with dimethylformamide and ethylenediamine, nitrogen-containing groups are added, the amount of which, however, is much less than the amount of fluorine. This modification of the SWCNTs improves intertube contacts in the film and increases the surface area for water adsorption. Our results show that an increase in the number of functional groups on the SWCNT surface enhances the sensitivity of the sensor to low water concentrations and worsens the response at high RH. SWCNTs modified with ethylenediamine have the highest sensitivity over the entire range of RH.

Towards an accurate description of one-dimensional pnictogen allotropes in nano-confinements

One-dimensional (1D) confined pnictogen shows a diverse range of allotropes and potential applications in electronic devices and the chemical industry. Here, we report a theoretical study aimed at an accurate assessment of the thermodynamic stability of pnictogen structures under nano-meter confinements. We develop a cylindrical potential for pnictogen, which can be integrated with density functional theory to model a confined system towards achieving ab initio accuracy. We discuss in detail the performance of confining potentials and provide insights into the understanding of complex interactions between confined pnictogen and carbon nanotubes. We reassess the thermodynamic stability of 1D pnictogens in carbon nanotubes, explaining the diverse features of confined pnictogens in recent experimental and theoretical studies.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

Multi-scale simulation of the adsorption of lithium ion on graphite surface: From quantum Monte Carlo to molecular density functional theory

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic density functional theory and molecular dynamics. Here, we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference quantum Monte Carlo calculations that allow us to benchmark various electronic density functional theory functionals. We then fit an accurate carbon–lithium pair potential, which is used in molecular density functional theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in aqueous solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer.

The interaction between atomic-scale pores and particles

Using first-principles calculations for angstrom-sized pores (3-10 Å), we investigate pore-particle interaction. The translocation energy barrier (TEB) plays important role for the angstrom-scale pores created in 2D-materials such as graphene which is calculated for the translocation of rare gases (He, Ne, Ar, Xe), diatomic molecules (H₂and N₂), CO₂, and CH₄. The critical incident angle (the premeance beyond that is zero) was found to be 40°, which is different from classical model's prediction of 19-37°. The calculated TEB (Δ) and the surface diffusion energy barrier (Δ') for the particles with small kinetic diameter (He, Ne and H₂), show that the direct flow is the dominant permeation mechanism (Δ ≈ 0 and Δ' > 30 meV). For the other particles with larger kinetic diameters (Ar, Kr, N₂, CH₄and CO₂), we found that both surface diffusion and direct flow mechanisms are possible, i.e. Δ and Δ' ≠ 0. This work provides important insights into the gas permeation theory and into the design and development of gas separation and filtration devices.

Interactions between large molecules pose a puzzle for reference quantum mechanical methods

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.

Molecular Simulation of Electrode-Solution Interfaces

Many key industrial processes, from electricity production, conversion, and storage to electrocatalysis or electrochemistry in general, rely on physical mechanisms occurring at the interface between a metallic electrode and an electrolyte solution, summarized by the concept of an electric double layer, with the accumulation/depletion of electrons on the metal side and of ions on the liquid side. While electrostatic interactions play an essential role in the structure, thermodynamics, dynamics, and reactivity of electrode-electrolyte interfaces, these properties also crucially depend on the nature of the ions and solvent, as well as that of the metal itself. Such interfaces pose many challenges for modeling because they are a place where quantum chemistry meets statistical physics. In the present review, we explore the recent advances in the description and understanding of electrode-electrolyte interfaces with classical molecular simulations, with a focus on planar interfaces and solvent-based liquids, from pure solvent to water-in-salt electrolytes.

Mechanical hydrolysis imparts self-destruction of water molecules under steric confinement

Decoding behavioral aspects associated with the water molecules in confined spaces such as an interlayer space of two-dimensional nanosheets is key for the fundamental understanding of water-matter interactions and identifying unexpected phenomena of water molecules in chemistry and physics. Although numerous studies have been conducted on the behavior of water molecules in confined spaces, their reach stops at the properties of the planar ice-like formation, where van der Waals interactions are the predominant interactions and many questions on the confined space such as the possibility of electron exchange and excitation state remain unsettled. We used density functional theory and reactive molecular dynamics to reveal orbital overlap and induction bonding between water molecules and graphene sheets under much less pressure than graphene fractures. Our study demonstrates high amounts of charge being transferred between water and the graphene sheets, as the interlayer space becomes smaller. As a result, the inner face of the graphene nanosheets is functionalized with hydroxyl and epoxy functional groups while released hydrogen in the form of protons either stays still or traverses a short distance inside the confined space via the Grotthuss mechanism. We found signatures of a new hydrolysis mechanism in the water molecules, i.e. mechanical hydrolysis, presumably responsible for relieving water from extremely confined conditions. This phenomenon where water reacts under extreme confinement by disintegration rather than forming ice-like structures is observed for the first time, illustrating the prospect of treating ultrafine porous nanostructures as a driver for water splitting and material functionalization, potentially impacting the modern design of nanofilters, nanochannels, nano-capacitators, sensors, and so on.

Polycrystalline Few-Layer Graphene as a Durable Anticorrosion Film for Copper

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.

Understanding water slippage through carbon nanotubes

It is a formidable challenge to understand water slippage through carbon nanotubes (CNTs), despite its great significance in fundamental research and technology. Herein, we propose an effective scheme to describe water slippage properties by extending two friction models - the phononic friction model and Einstein's diffusion model, both relying on the potential corrugation of water slippage. Our scheme effectively captures the tube-size effect on the viscosity and slippage of water molecules through CNTs. It also identifies the experimentally reported size-dependent transition from continuum to sub-continuum flow and further reveals that this transition is likely to be determined by the hydrogen bond instead of the structural transition or entropic change. Besides, the size-dependence of slip lengths is found to be controllable by temperature. Our methods are thus expected to be a useful basis for further studies on substance transport under confinement.

Bromine polycondensation in pristine and fluorinated graphitic carbons

Despite decades of study the precise behavior of bromine in graphitic carbons remains unclear. In this report, using Raman spectroscopy, we reveal two types of bromine structure in graphitic carbon materials. Between fluorinated graphene layers with a composition close to C2F, Br2 molecules are intercalated in a form similar to liquid bromine. Bromination of pristine and low-fluorinated graphitic carbons behaves very differently with distinct Br-related Raman spectra. With the guidance of density functional theory (DFT) calculations, all Raman features are assigned to normal vibration modes of specific bromine species over graphene and fluorinated graphene. When intercalated between extended non-fluorinated sp2-hybridized carbon regions, physisorbed Br2 molecules move freely across the non-functionalized region toward the CF border. Multiple Br2 molecules then combine spontaneously into Br3-based chains, whose coupling activates otherwise Raman inactive modes. Significant charge transfer to bromine species occurs in this case. DFT calculated frequencies match precisely the experimental Br-related Raman bands observed in the intercalation carbon compounds. The fluorine-catalyzed bromine chain-formation process shown here is general and should also operate with edges and other defect species.

Liquids with Lower Wettability Can Exhibit Higher Friction on Hexagonal Boron Nitride: The Intriguing Role of Solid-Liquid Electrostatic Interactions

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) × 10⁵ N-s/m³ 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.

Physisorption of Water on Graphene: Subchemical Accuracy from Many-Body Electronic Structure Methods

Wet carbon interfaces are ubiquitous in the natural world and exhibit anomalous properties, which could be exploited by emerging technologies. However, progress is limited by lack of understanding at the molecular level. Remarkably, even for the most fundamental system (a single water molecule interacting with graphene), there is no consensus on the nature of the interaction. We tackle this by performing an extensive set of complementary state-of-the-art computer simulations on some of the world's largest supercomputers. From this effort a consensus on the water-graphene interaction strength has been obtained. Our results have significant impact for the physical understanding, as they indicate that the interaction is weaker than predicted previously. They also pave the way for more accurate and reliable studies of liquid water at carbon interfaces.

Understanding non-covalent interactions in larger molecular complexes from first principles

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.

Theory and practice of modeling van der Waals interactions in electronic-structure calculations

The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.

First-Principles Study of Water Nanotubes Captured Inside Carbon/Boron Nitride Nanotubes

Water confined to nanopores such as carbon nanotubes (CNTs) exhibits different states, enabling the study of solidlike water nanotubes (WNTs) and the potential application of their properties due to confined effects. Herein, we report the interfacial interaction and particular stabilized boundaries of confined WNTs within CNTs and boron nitride nanotubes (BNNTs) using first-principles calculations. We demonstrate that the intermolecular potential of nanotube walls exerts diameter-dependent additive or subtractive van der Waals (vdW) pressure on the WNTs, altering the phase boundaries. Our results reveal that the most stable WNT@CNT is associated with a CNT diameter of 10.5 Å. By correlating the stability of WNTs with interfacial properties such as the vdW pressure and vibrational phonon modes of confined WNTs, we decode and compare various synergies in water interaction and stabilized states within the CNTs and BNNTs, including interfacial properties of WNT@BNNTs that are more significant than those of WNT@CNTs. Our results suggest that the transition of a water tube to an ice tube is strongly dependent on the diameter of the confining CNT or BNNT, providing new insights on leveraging the interfacial interaction mechanism of confined WNTs and their potential application for fabricating nanochannels and nanocapacitors.

Exploring dissociative water adsorption on isoelectronically BN doped graphene using alchemical derivatives

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.

Understanding Fast and Robust Thermo-osmotic Flows through Carbon Nanotube Membranes: Thermodynamics Meets Hydrodynamics

Following our recent theoretical prediction of the giant thermo-osmotic response of the water-graphene interface, we explore the practical implementation of waste heat harvesting with carbon-based membranes, focusing on model membranes of carbon nanotubes (CNT). To that aim, we combine molecular dynamics simulations and an analytical model considering the details of hydrodynamics in the membrane and at the tube entrances. The analytical model and the simulation results match quantitatively, highlighting the need to take into account both thermodynamics and hydrodynamics to predict thermo-osmotic flows through membranes. We show that, despite viscous entrance effects and a thermal short-circuit mechanism, CNT membranes can generate very fast thermo-osmotic flows, which can overcome the osmotic pressure of seawater. We then show that in small tubes confinement has a complex effect on the flow and can even reverse the flow direction. Beyond CNT membranes, our analytical model can guide the search for other membranes to generate fast and robust thermo-osmotic flows.

What Controls Thermo-osmosis? Molecular Simulations Show the Critical Role of Interfacial Hydrodynamics

Thermo-osmotic and related thermophoretic phenomena can be found in many situations from biology to colloid science, but the underlying molecular mechanisms remain largely unexplored. Using molecular dynamics simulations, we measure the thermo-osmosis coefficient by both mechanocaloric and thermo-osmotic routes, for different solid-liquid interfacial energies. The simulations reveal, in particular, the crucial role of nanoscale interfacial hydrodynamics. For nonwetting surfaces, thermo-osmotic transport is largely amplified by hydrodynamic slip at the interface. For wetting surfaces, the position of the hydrodynamic shear plane plays a key role in determining the amplitude and sign of the thermo-osmosis coefficient. Finally, we measure a giant thermo-osmotic response of the water-graphene interface, which we relate to the very low interfacial friction displayed by this system. These results open new perspectives for the design of efficient functional interfaces for, e.g., waste-heat harvesting.

Simulations of water nano-confined between corrugated planes

Water confined to nanoscale widths in two dimensions between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water to the imposition of a periodicity in the confinement by means of computer simulations, both using empirical potentials and from first-principles. For that we propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice and other values that would correspond to more realistic substrates. For the former, the phase diagram shows an overall rise of the melting temperature. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude, the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials (TIP4P/2005), both agree remarkably well for the main conclusions of the study.