The water-benzene interaction: Insight from electronic structure theories

Tunable Surface Wettability via Terahertz Electrowave Controlled Vicinal Subnanoscale Water Layer

Achieving timely, reversible, and long-range remote tunability over surface wettability is highly demanded across diverse fields, including nanofluidic systems, drug delivery, and heterogeneous catalysis. Herein, using molecular dynamic simulations, we show, for the first time, a theoretical design of electrowetting to achieve remotely controllable surface wettability via using a terahertz wave. The key idea driving the design is the unique terahertz collective vibration identified in the vicinal subnanoscale water layer, which is absent in bulk water, enabling efficient energy transfer from the terahertz wave to the rotational motion of the vicinal subnanoscale water layer. Consequently, a frequency-specific alternating terahertz electric field near the critical strength can significantly affect the local hydrogen-bonding network of the contact water layer on the solid surface, thereby achieving tunable surface wettability.

Accuracy of Noncovalent Interactions Involving d-Elements by the 1-Determinant Fixed-Node Diffusion Monte Carlo Method with Effective Core Potentials

A critical assessment of effective core potential (ECP)-based single-determinant (SD) fixed-node diffusion quantum Monte Carlo (FNDMC) accuracy in prototypical noncovalent closed-shell systems involving d-elements is presented. Careful analysis of biases and elimination of possible bias sources leads to two findings of practical importance for SD FNDMC in these systems. First, in some systems (HCu:HCu, HCu:CuH), SD FNDMC reveals large biases of interaction energy differences (significantly exceeding the target 2% relative error) vs a reliable coupled-cluster CCSD(T)/CBS (complete basis set) reference. Second, the leading error of SD FNDMC with ECPs was attributed to a higher nuclear charge Z of d-group (pseudo) atoms, when compared to sp elements, in line with a previously reported finding that aggregate SD FNDMC bias tends to increase in systems with higher electronic densities. Therefore, SD FNDMC should only be used with caution in systems with a large Z.

Putting together the puzzle of ion transfer in single-digit carbon nanotubes: mean-field meets ab initio

Nature employs channel proteins to selectively pass water across cell membranes, which inspires the search for bio-mimetic analogues. Carbon nanotube porins (CNTPs) are intriguing mimics of water channels, yet ion transport in CNTPs still poses questions. As an alternative to continuum models, here we present a molecular mean-field model that transparently describes ion coupling, yet unlike continuum models, computes ab initio all required thermodynamic quantities for the KCl salt and H⁺ and OH^- ions present in water. Starting from water transfer, the model considers the transfer of free ions, along with ion-pair formation as a proxy of non-mean-field ion-ion interactions. High affinity to hydroxide, suggested by experiments, making it a dominant charge carrier in CNTPs, is revealed as an exceptionally favorable transfer of KOH pairs. Nevertheless, free ions, coexisting with less mobile ion-pairs, apparently control ion transport. The model well explains the observed effects of salt concentration and pH on conductivity, transport numbers, anion permeation and its activation energies, and current rectification. The proposed approach is extendable to other sub-nanochannels and helps design novel osmotic materials and devices.

Isotope Effects on the Vaporization of Organic Compounds from an Aqueous Solution-Insight from Experiment and Computations

An isotope fractionation analysis of organic groundwater pollutants can assess the remediation at contaminated sites yet needs to consider physical processes as potentially confounding factors. This study explores the predictability of water–air partitioning isotope effects from experiments and computational predictions for benzene and trimethylamine (both H-bond acceptors) as well as chloroform (H-bond donor). A small, but significant, isotope fractionation of different direction and magnitude was measured with ε = −0.12‰ ± 0.07‰ (benzene), ε_C = 0.49‰ ± 0.23‰ (triethylamine), and ε_H = 1.79‰ ± 0.54‰ (chloroform) demonstrating that effects do not correlate with expected hydrogen-bond functionalities. Computations revealed that the overall isotope effect arises from contributions of different nature and extent: a weakening of intramolecular vibrations in the condensed phase plus additional vibrational modes from a complexation with surrounding water molecules. Subtle changes in benzene contrast with a stronger coupling between intra- and intermolecular modes in the chloroform–water system and a very local vibrational response with few atoms involved in a specific mode of triethylamine. An energy decomposition analysis revealed that each system was affected differently by electrostatics and dispersion, where dispersion was dominant for benzene and electrostatics dominated for chloroform and triethylamine. Interestingly, overall stabilization patterns in all studied systems originated from contributions of dispersion rather than other energy terms.

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

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.

Dynamics of water in conical solvation shells around a benzene solute under different thermodynamic conditions

Water molecules in different parts of the anisotropic hydration shell of an aromatic molecule experience different interactions. In the present study, we investigate the anisotropic dynamics of water molecules in different non-overlapping conical shells around a benzene solute at sub- and supercritical conditions by means of molecular dynamics simulations using both non-polarizable and polarizable models. In addition to the dynamical properties, the effects of polarizability on the hydration structure of benzene at varying thermodynamic conditions are also investigated in the current study. The presence of πH-bonding in the solvation shell is found to be an important factor that influences the anisotropic dynamics of the benzene hydration shell. The water molecules located axial to the benzene plane are found to be maximally influenced by the πH-bonding. The extent of πH-bonding is found to be somewhat reduced on inclusion of polarizability. The πH-bonded water molecules are found to reorient through large-amplitude angular jumps where the jump-angle amplitude increases at higher temperatures and lower densities. For both non-polarizable and polarizable models, it is found that the water molecules in the axial conical shells possess faster orientational and hydrogen bond dynamics compared to those in the equatorial plane. Water molecules in the axial conical shells are also found to diffuse at a faster rate than bulk molecules due to the relatively weaker benzene-water πH-bonding interactions in the axial region of the hydration shell. The residence dynamics of water molecules in different conical solvation shells around the solute is also investigated in the current study.

Bias cancellation in one-determinant fixed-node diffusion Monte Carlo: Insights from fermionic occupation numbers

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.

Behavior of interactions between hydrogen chalcogenides and an anthracene π-system elucidated by QTAIM dual functional analysis with QC calculations

The nature of the interactions between chalcogenides and the anthracene p-system, EH₂-*-p(C₁₄H₁₀), is predicted to be close to that of EH₂-*-p(C₁₀H₈), although the partial structures around the central rings can be found in EH₂-*-p(C₆H₆).

Does borazine-water behave like benzene-water? A matrix isolation infrared and ab initio study

Borazine is isoelectronic with benzene and is popularly referred to as inorganic benzene. The study of non-covalent interactions with borazine and comparison with its organic counterpart promises to show interesting similarities and differences. The motivation of the present study of the borazine-water interaction, for the first time, stems from such interesting possibilities. Hydrogen-bonded complexes of borazine and water were studied using matrix isolation infrared spectroscopy and quantum chemical calculations. Computations were performed at M06-2X and MP2 levels of theory using 6-311++G(d,p) and aug-cc-pVDZ basis sets. At both the levels of theory, the complex involving an N-H⋯O interaction, where the N-H of borazine serves as the proton donor to the oxygen of water was found to be the global minimum, in contrast to the benzene-water system, which showed an H-π interaction. The experimentally observed infrared spectra of the complexes corroborated well with our computations for the complex corresponding to the global minimum. In addition to the global minimum, our computations also located two local minima on the borazine-water potential energy surface. Of the two local minima, one corresponded to a structure where the water was the proton donor to the nitrogen of borazine, approaching the borazine ring from above the plane of the ring; a structure that resembled the global minimum in the benzene-water H-π complex. The second local minimum corresponded to an interaction of the oxygen of water with the boron of borazine, which can be termed as the boron bond. Clearly the borazine-water system presents a richer landscape than the benzene-water system.

Dynamic and static behavior of the H...π and E...π interactions in EH₂ adducts of benzene π-system (E = O, S, Se and Te), elucidated by QTAIM dual functional analysis

Dynamic and static behavior of the interactions in the EH2 adducts of a benzene π-system (E = O, S, Se and Te) is elucidated by applying QTAIM-DFA (QTAIM dual functional analysis). Two types of H-*-π and E-*-π interactions are detected in the adducts, where the asterisk (*) emphasizes the existence of the bond critical point (BCP) on the interaction in question. Total electron energy densities Hb(rc) are plotted versus Hb(rc) -Vb(rc)/2 [=(ℏ(2)/8m)∇(2)ρb(rc)] at BCPs in QTAIM-DFA, where Vb(rc) are the potential energy densities at BCPs. Data from the fully optimized structures are analyzed by polar (R, θ) coordinate representation. Each plot for an interaction, containing data from the perturbed structures with those of the fully optimized one, shows a specific curve, which provides important information. The plot is expressed by (θp, κp): θp corresponds to the tangent line for the plot and κp is the curvature. θ and θp are measured from the y-axis and y-direction, respectively. Moreover, (R, θ) corresponds to the static nature, (θp, κp) represents the dynamic nature of interactions. While θ classifies the interaction in question, θp characterizes it. Both values are less than 90° for all H-*-π and E-*-π interactions examined in this study; therefore, they are all classified by the pure closed-shell interactions and predicted to have the character of vdW nature. However, it is suggested that E-*-π has the nature of the stronger interaction than the case of H-*-π for dynamic behavior in the same species evaluated at the MP2 and M06-2X levels. The nature of the interactions is well analyzed and specified by applying QTAIM-DFA.

Modelling of graphene functionalization

This perspective describes the available theoretical methods and models for simulating graphene functionalization based on quantum and classical mechanics.

Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations

The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.

Water on BN doped benzene: a hard test for exchange-correlation functionals and the impact of exact exchange on weak binding

Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces, whereas the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B3N3H6) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy.

Communication: Water on hexagonal boron nitride from diffusion Monte Carlo

Despite a recent flurry of experimental and simulation studies, an accurate estimate of the interaction strength of water molecules with hexagonal boron nitride is lacking. Here, we report quantum Monte Carlo results for the adsorption of a water monomer on a periodic hexagonal boron nitride sheet, which yield a water monomer interaction energy of -84 ± 5 meV. We use the results to evaluate the performance of several widely used density functional theory (DFT) exchange correlation functionals and find that they all deviate substantially. Differences in interaction energies between different adsorption sites are however better reproduced by DFT.

On the room-temperature phase diagram of high pressure hydrogen: an ab initio molecular dynamics perspective and a diffusion Monte Carlo study

The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results.

Benchmarking the performance of density functional theory and point charge force fields in their description of sI methane hydrate against diffusion Monte Carlo

High quality reference data from diffusion Monte Carlo calculations are presented for bulk sI methane hydrate, a complex crystal exhibiting both hydrogen-bond and dispersion dominated interactions. The performance of some commonly used exchange-correlation functionals and all-atom point charge force fields is evaluated. Our results show that none of the exchange-correlation functionals tested are sufficient to describe both the energetics and the structure of methane hydrate accurately, while the point charge force fields perform badly in their description of the cohesive energy but fair well for the dissociation energetics. By comparing to ice Ih, we show that a good prediction of the volume and cohesive energies for the hydrate relies primarily on an accurate description of the hydrogen bonded water framework, but that to correctly predict stability of the hydrate with respect to dissociation to ice Ih and methane gas, accuracy in the water-methane interaction is also required. Our results highlight the difficulty that density functional theory faces in describing both the hydrogen bonded water framework and the dispersion bound methane.

Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur

Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.

Intermolecular π/π and H/π interactions in dimers researched by different computational methods

The analysis of π/π and H /π interactions in complexes are a challenging aspect of theoretical research. Due to the different approximations of different levels of theory, results tend to be inconsistent. We compared the reliabilities of HF, SVWN, M06L, PW91, BLYP, B3LYP, BHandHLYP, B97D, MP2, and DFTB-D approaches in researching π/π and H /π interactions by calculating the binding energies of five benzene-containing dimers. The effects of 6-31+G**, 6-311++G** and 6-311++G(2df,2p) basis sets on the results were analyzed too. We found that the DFTB-D and B97D methods combined with the 6-311++G** basis set perform well for dimers that contain π/π and H /π interactions. With high efficiency and satisfactory precision, DFTB-D is helpful for the calculation of complexes containing π/π and H /π stacking. We further calculated the structures and properties of phenylalanine-containing dimers using the DFTB-D and B97D methods. The properties of low energy conformers such as rotational constants, dipole moments and molecular orbitals were also analyzed. These data should be helpful for research into systems that contain π/π and H /π stacking.

Coexistence of solvated electron and benzene-centered valence anion in the negatively charged benzene-water clusters

We present a combined M06 functional calculation and ab initio molecular dynamics simulation study of an excess electron (EE) in a microhydrated aromatic complex (modeled by benzene (Bz)-water binary clusters, Bz(H(2)O)(n)). Calculated results illustrate that Bz ring and water clusters are indeed linked through the π···HO interactions in the neutral Bz(H(2)O)(n) (n = 1-8) clusters, and the size of the water cluster does not influence the nature of its interaction with the π system for the oligo-hydrated complexes. The states and the dynamics of an EE trapped in such Bz-water clusters were also determined. All of possible localized states for the EE can be roughly classified into two types: (i) single, ring-localized states (the Bz-centered valence anions) in which an EE occupies the LUMO of the complexes originating from the LUMO (π*) of the Bz ring, and the π···HO interactions are enhanced for increase of electron density of the Bz ring. In this mode, the carbon skeleton of the Bz part is significantly deformed due to increase of electron density and nonsymmetric distribution of electron density induced by the interacting H-O bonds; (ii) solvated states, in which an EE is trapped directly as a surface state by the dangling hydrogen atoms of water molecules or as a solvated state in a mixed cavity formed by Bz and water cluster. In the latter case, Bz may also participate in capturing an EE using its C-H bonds in the side edge of the aromatic ring as a part of the cavity. In general, a small water cluster is favorable to the Bz-centered valence anion state, while a large one prefers a solvated electron state. Fluctuations and rearrangement of water molecules can sufficiently modify the relative energies of the EE states to permit facile conversion from the Bz-centered to the water cluster-centered state. This indicates that aromatic Bz can be identified as a stepping stone in electron transfer and the weak π···HO interaction plays an important role as the driving force in conversion of the two states.

Studies on the structure, stability, and spectral signatures of hydride ion-water clusters

The gas-phase structure, stability, spectra, and electron density topography of H(-)W(n) clusters (where n = 1-8) have been calculated using coupled-cluster CCSD(T) and Møller-Plesset second-order perturbation (MP2) theory combined with complete basis set (CBS) approaches. The performance of various density functional theory (DFT) based methods such as B3LYP, M05-2X, M06, M06-L, and M06-2X using 6-311++G(d,p), and aug-cc-pVXZ (aVXZ, where X = D, T, and Q) basis sets has also been assessed by considering values calculated using CCSD(T)/CBS limit as reference. The performance of the functionals has been ranked based on the mean signed/unsigned error. The comparison of geometrical parameters elicits that the geometrical parameters predicted by B3LYP/aVTZ method are in good agreement with those values obtained at MP2/aVTZ level of theory. Results show that M05-2X functional outperform other functionals in predicting the energetics when compared to CCSD(T)/CBS value. On the other hand, values predicted by M06-2X, and M06 methods, are closer to those values obtained from MP2/CBS approach. It is evident from the calculations that H(-)W(n) (where n = 5-8) clusters adopt several interesting structural motifs such as pyramidal, prism, book, Clessidra, cubic, cage, and bag. The important role played by ion-water (O-H···H(-)) and water-water (O-H···O) interactions in determining the stability of the clusters has also been observed. Analysis of the results indicates that the most stable cluster is made up of minimum number of O-H···H(-) interaction in conjugation with the maximum number of O-H···O interactions. The Bader theory of atoms in molecules (AIM) and natural bond orbital (NBO) analyses has also been carried out to characterize the nature of interactions between hydride ion and water molecules. It can be observed from the vibrational spectra of H(-)W(n) clusters, the stretching frequencies involving ion-water interaction always exhibit larger redshift and intensities than that of water-water (inter solvent) interactions.