Anisotropic Interfacial Force Field for Interfaces of Water with Hexagonal Boron Nitride

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.

Diameter-Optimum Spreading for the Impinging of Water Nanodroplets on Solid Surfaces

The impinging of water nanodroplets on solid surfaces is crucial to many nanotechnologies. Through large-scale molecular dynamics simulations, the size effect on the spreading of water nanodroplets after impinging on hydrophilic, graphite, and hydrophobic surfaces under low impinging velocities has been systematically studied. The spreading rates of nanodroplets first increase and then decrease and gradually become constant with the increase of nanodroplet diameter. The nanodroplets with the diameters of 17-19 nm possess the highest spreading rates because of the combined effect of the strongest interfacial interaction and the strongest surface interaction within water molecules. The highest water molecule densities, hydrogen bond numbers, and dielectric constants of interface and surface layers mainly contribute to the lowest interface work of adhesion and surface tension values at optimal diameters. These results unveil the nonmonotonic characteristics of spreading velocity, interface work of adhesion and surface tension with nanodroplet diameter for nanodroplets on solid surfaces.

Liquid-Solid Slip on Charged Walls: The Dramatic Impact of Charge Distribution

Nanofluidic systems show great promise for applications in energy conversion, where their performance can be enhanced by nanoscale liquid-solid slip. However, efficiency is also controlled by surface charge, which is known to reduce slip. Combining molecular dynamics simulations and analytical developments, we show the dramatic impact of surface charge distribution on the slip-charge coupling. Homogeneously charged graphene exhibits a very favorable slip-charge relation (rationalized with a new theoretical model correcting some weaknesses of the existing ones), leading to giant electrokinetic energy conversion. In contrast, slip is strongly affected on heterogeneously charged surfaces, due to the viscous drag induced by counterions trapped on the surface. In that case slip should depend on the detailed physical chemistry of the interface controlling the fraction of bound ions. Our numerical results and theoretical models provide new fundamental insight into the molecular mechanisms of liquid-solid slip, and practical guidelines for searching new functional interfaces with optimal energy conversion properties, e.g., for blue energy or waste heat harvesting.

Replica exchange MD simulations of two-dimensional water in graphene nanocapillaries: rhombic versus square structures, proton ordering, and phase transitions

The hydrogen bond patterns, proton ordering, and phase transitions of monolayer ice in two-dimensional hydrophobic confinement are fundamentally different from those found for bulk ice. To investigate the behavior of quasi-2D ice, we perform molecular dynamics simulations of water confined between fixed graphene plates at a distance of 0.65 nm. While experimental results are still limited and theoretical investigations are often based on a single, often empirically based force field model, this work presents a systematic study modeling the water-graphene interaction by effective Lennard-Jones potentials previously derived from high-level ab initio CCSD(T) calculations of water adsorbed on graphene [Phys. Chem. Chem. Phys., 2013, 15, 4995]. For the water-water interaction different water force fields, i.e. SPCE, TIP3P, TIP4P, TIP4P/ICE, and TIP5P, are used. The water occupancy of the graphene capillary at a pressure of 1000 MPa is determined to be between 13.5 and 13.9 water molecules per square nanometer, depending on the choice of the water force field. Based on these densities, we explore the structure and dynamics of quasi-2D water for temperatures ranging from 200 K to about 600 K for each of the five force fields. To ensure complete sampling of the configurational space and to overcome the barriers separating metastable structures, these simulations are based on the replica exchange molecular dynamics technique. We report different tetragonal hydrogen bond patterns, which are classified as nearly square or as rhombic. While many of these arrangements are essentially flat, in some cases puckered arrangements are found, too. Also the proton ordering of the quasi-2D water structures is considered, allowing us to identify them as ferroelectric, ferrielectric or antiferroelectric. For temperatures between 200 K and 400 K we find several second-order phase transitions from one ice structure to another, changing in many cases both the arrangements of the oxygen atoms and the proton ordering. For temperatures between 400 K and 600 K there are melting-like transitions from a monolayer of ice to a monolayer of liquid water. These first-order phase transitions have a latent heat between 3.4 and 4.0 kJ mol^-1. Both the values of the transition temperatures and of the latent heats display considerable model dependence for the five different water models investigated here.

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.

Interruption of Hydrogen Bonding Networks of Water in Carbon Nanotubes Due to Strong Hydration Shell Formation

We present the structures of NaCl aqueous solution in carbon nanotubes with diameters of 1, 2, and 3 nm based on an analysis performed using X-ray diffraction and canonical ensemble Monte Carlo simulations. Anomalously longer nearest-neighbor distances were observed in the electrolyte for the 1-nm-diameter carbon nanotubes; in contrast, in the 2 and 3 nm carbon nanotubes, the nearest-neighbor distances were shorter than those in the bulk electrolyte. We also observed similar properties for water in carbon nanotubes, which was expected because the main component of the electrolyte was water. However, the nearest-neighbor distances of the electrolyte were longer than those of water in all of the carbon nanotubes; the difference was especially pronounced in the 2-nm-diameter carbon nanotubes. Thus, small numbers of ions affected the entire structure of the electrolyte in the nanopores of the carbon nanotubes. The formation of strong hydration shells between ions and water molecules considerably interrupted the hydrogen bonding between water molecules in the nanopores of the carbon nanotubes. The hydration shell had a diameter of approximately 1 nm, and hydration shells were thus adopted for the nanopores of the 2-nm-diameter carbon nanotubes, providing an explanation for the large difference in the nearest-neighbor distances between the electrolyte and water in these nanopores.

Molecular origin of drug release by water boiling inside carbon nanotubes from reactive molecular dynamics simulation and DFT perspectives

Owing to their nanosized hollow cylindrical structure, CNTs hold the promise to be utilized as desired materials for encapsulating molecules which demonstrate wide inferences in drug delivery. Here we evaluate the possibility of drug release from the CNTs with various types and edge chemistry by reactive MD simulation to explain the scientifically reliable relations for proposed process. It was shown that heating of CNTs (up to 750 K) cannot be used for release of incorporated drug (phenylalanine) into water and even carbonated water solvent with very low boiling temperature. This is due to the strong physisorption (π-stacking interaction) between the aromatic of encapsulated drug and CNT sidewall which causes the drug to bind the nanotube sidewall. We have further investigated the interaction nature and release mechanism of water and drug confined/released within/from the CNTs by DFT calculations and the results confirmed our MD simulation findings. The accuracy of DFT method was also validated against the experimental and theoretical values at MP2/CCSD level. Therefore, we find that boiling of water/carbonated water confined within the CNTs could not be a suitable technique for efficient drug release. Our atomistic simulations provide a well-grounded understanding for the release of drug molecules confined within CNTs.

Solubility origin at the nanoscale: enthalpic and entropic contributions in polar and nonpolar environments

Nanostructures are known to be poorly soluble, irrespective of their elemental composition, shape, electronic structure, dipole moment, hydrophobicity/hydrophilicity and the employed solvent. The methods of colloid chemistry allow for preparing suspensions - metastable systems, the stabilities of which differ greatly from one another - but not real solutions. A systematic investigation of the solubility origin at the nanoscale is hereby reported in terms of its fundamental constituents: enthalpy and entropy. Slightly different one-dimensional solutes - narrow carbon nanotubes (CNTs) of different lengths - were considered in hydrophilic (water) and hydrophobic (benzene) environments. We decompose the process of solvation into the solid → gas transition (sublimation) and the gas → liquid transition (condensation). Sublimation is a thermodynamically unfavorable process under room conditions, while the condensation transition depends on the solvent-solute interactions (enthalpic contribution). Unlike solvation of small molecules, solvation of the nanostructures results in a significant alteration of entropy. This alteration is proportional to the linear dimensions of the nanostructure. If the solvent exhibits peculiar solvent-solvent interactions (such as hydrogen bonding in water), solvation is entropically forbidden, irrespective of the solute nature and its nanoscale dimensions. In the case of the hydrophobic solvent (benzene), the condensation transition can be both enthalpically and entropically favorable. The free energy of solvation is in direct proportion to the CNT length. While highlighting principal difficulties in solvating nanostructures, this paper discusses an optimal choice of solvents for solutes exhibiting hydrophobic and hydrophilic interactions with their environments. Our results allow us to predict the solvation of an arbitrary nanostructure using its small, about 2 nm, atomistic model.

Molecular dynamics simulations of electric field induced water flow inside a carbon nanotorus: a molecular cyclotron

A nano-flow is induced by applying gigahertz rotating electric fields (EFs) of different strengths and frequencies on a carbon nanotorus filled with water molecules, using molecular dynamics simulations.

Effects of confinement between attractive and repulsive walls on the thermodynamics of an anomalous fluid

We study via molecular-dynamics simulations the thermodynamics of an anomalous fluid confined in a slit pore with one wall structured and attractive and another unstructured and repulsive. We find that the phase diagram of the homogeneous part of the confined fluid is shifted to higher temperatures, densities, and pressures with respect to the bulk, but it can be rescaled on the bulk case. We calculate a moderate increase of mobility of the homogeneous confined fluid that we interpret as a consequence of the layering due to confinement and the collective modes due to long-range correlations. We show that, as in bulk, the confined fluid has structural, diffusion, and density anomalies that order in the waterlike hierarchy, and a liquid-liquid critical point (LLCP). The overall anomalous region moves to higher temperatures, densities, and pressure, and the LLCP displaces to higher temperature compared to bulk. Motivated by experiments, we calculate also the phase diagram not just for the homogeneous part of the confined fluid but for the entire fluid in the pore, and we show that it is shifted toward higher pressures but preserves the thermodynamics, including the LLCP. Because our model has waterlike properties, we argue that in experiments with supercooled water confined in slit pores with a width of >3 nm if hydrophilic and of >1.5 nm if hydrophobic, the existence of the LLCP could be easier to test than in bulk, where it is not directly accessible.

Spontaneous Insertion, Helix Formation, and Hydration of Polyethylene Oxide in Carbon Nanotubes

Hydration strongly affects macromolecular conformation in solution and under nanoconfinement as encountered in nature and nanomaterials. Using atomistic molecular dynamics simulations we demonstrate that polyethylene oxide spontaneously enters single wall carbon nanotubes (CNTs) from aqueous solutions and forms rodlike, helix, and wrapped chain conformations depending on the CNT diameter. We show that water organization and the stability of the polyethylene oxide hydration shell under confinement is responsible for the helix formation, which can have significant implications for nanomaterial design.

Curvature-dependent adsorption of water inside and outside armchair carbon nanotubes

The curvature dependence of the physisorption properties of a water molecule inside and outside an armchair carbon nanotube (CNT) is investigated by an incremental density-fitting local coupled cluster treatment with single and double excitations and perturbative triples (DF-LCCSD(T)) study. Our results show that a water molecule outside and inside (n, n) CNTs (n = 4, 5, 6, 7, 8, 10) is stabilized by electron correlation. The adsorption energy of water inside CNTs decreases quickly with the decrease of curvature (increase of radius) and the configuration with the oxygen pointing toward the CNT wall is the most stable one. However, when the water molecule is adsorbed outside the CNT, the adsorption energy varies only slightly with the curvature and the configuration with hydrogens pointing toward the CNT wall is the most stable one. We also use the DF-LCCSD(T) results to parameterize Lennard-Jones (LJ) force fields for the interaction of water both with the inner and outer sides of CNTs and with graphene representing the zero curvature limit. It is not possible to reproduce all DF-LCCSD(T) results for water inside and outside CNTs of different curvature by a single set of LJ parameters, but two sets have to be used instead. Each of the two resulting sets can reproduce three out of four minima of the effective potential curves reasonably well. These LJ models are then used to calculate the water adsorption energies of larger CNTs, approaching the graphene limit, thus bridging the gap between CNTs of increasing radius and flat graphene sheets.

Confinement and surface effects of aqueous solutions within charged carbon nanotubes

Size-charge effects of brines in charged nanotubes from a molecular dynamics investigation of ion hydration, water coordination, and hydrogen bonding.

Molecular dynamics simulations of proton-ordered water confined in low-diameter carbon nanotubes

Metastable ferroelectric, ferrielectric, and antiferroelectric structures of water confined in carbon nanotubes and their abrupt or continuous transitions are identified.

What are preferred water-aromatic interactions in proteins and crystal structures of small molecules?

The distribution of water molecules around aromatic rings in the protein structures and crystal structures of small molecules shows quite a small number of the strongest OH-π interactions, a larger number of parallel interactions, and the largest number of the weakest CH-O interactions.