Dynamics of confined water inside carbon nanotubes based on studying tetrahedral order parameters

Water dynamics inside hydrophobic confinement, such as carbon nanotubes (CNTs), has garnered significant attention, focusing on water diffusion. However, a crucial aspect remains unexplored - the influence of confinement size on water ordering and intrinsic hydrogen bond dynamics. To address this gap, we conducted extensive molecular dynamics simulations to investigate local ordering and intrinsic hydrogen bond dynamics of water molecules within CNTs of various sizes (length:20 nm, diameters: 1.0 nm to 5.0 nm) over a wide range of temperatures (260K, 280K, 300K, and 320K). A striking observation emerged: in smaller CNTs, water molecules adopt an icy structure near tube walls while maintaining liquid state towards the center. Notably, water behavior within a 2.0 nm CNT stands out as an anomaly, distinct from other CNT sizes considered in this study. This anomaly was explained through the formation of water layers inside CNTs. The hydrogen bond correlation function of water within CNTs decayed more slowly than bulk water, with an increasing rate as CNT diameter increased. In smaller CNTs, water molecules hold onto their hydrogen bond longer than larger ones. Interestingly, in larger CNTs, the innermost layer's hydrogen bond lasts a shorter time compared to the other layers, and this changes with temperature.

Predicting Quantum-Mechanical Partial Charges in Arbitrarily Long Boron Nitride Nanotubes to Accurately Simulate Nanoscale Water Transport

Single-walled boron nitride nanotubes (BNNTs) have been explored for various applications, ranging from water desalination to osmotic power harvesting. However, no simulation work so far has modeled the changes in the partial charge distribution when a flat sheet is rolled into a tube, hindering the ability to perform accurate molecular dynamics (MD) simulations of water flow through BNNTs. To address this knowledge gap, we employ electronic density functional theory (DFT) calculations to precisely estimate quantum-mechanically derived partial charges on boron (B) and nitrogen (N) atoms in BNNTs of varying lengths and diameters. We observe a spatially varying charge distribution inside both armchair and zigzag nanotubes of finite lengths. Performing DFT calculations for longer BNNTs is computationally intractable, even with state-of-the-art computing resources. To solve this issue, we devise a charge assignment scheme to predict partial charges for longer BNNTs using DFT data for shorter nanotubes, thus overcoming the need to perform more expensive DFT calculations. We show that these charges reproduce the electrostatic potential predicted from first-principles simulations. Subsequently, we carried out MD simulations to predict the effect of the charge distribution inside BNNTs on water flow enhancement via them. We find that using uniform charges leads to an underprediction in flow enhancement, as compared to using quantum-mechanical charges for both armchair and zigzag BNNTs. We also incorporate atomic vibrations into our simulations and show that these vibrations lead to a reduction in the water flow through aperiodic BNNTs. Our work demonstrates the requirement of a quantum-mechanical charge assignment scheme for BNNTs and evolves a framework to assign charges to nanotubes of arbitrary length, thus allowing realistic MD simulations of long BNNTs using accurate partial charges.

Effects of polarizability and charge transfer on water dynamics and the underlying activation energies

A large number of force fields have been proposed for describing the behavior of liquid water within classical atomistic simulations, particularly molecular dynamics. In the past two decades, models that incorporate molecular polarizability and even charge transfer have become more prevalent, in attempts to develop more accurate descriptions. These are frequently parameterized to reproduce the measured thermodynamics, phase behavior, and structure of water. On the other hand, the dynamics of water is rarely considered in the construction of these models, despite its importance in their ultimate applications. In this paper, we explore the structure and dynamics of polarizable and charge-transfer water models, with a focus on timescales that directly or indirectly relate to hydrogen bond (H-bond) making and breaking. Moreover, we use the recently developed fluctuation theory for dynamics to determine the temperature dependence of these properties to shed light on the driving forces. This approach provides key insight into the timescale activation energies through a rigorous decomposition into contributions from the different interactions, including polarization and charge transfer. The results show that charge transfer effects have a negligible effect on the activation energies. Furthermore, the same tension between electrostatic and van der Waals interactions that is found in fixed-charge water models also governs the behavior of polarizable models. The models are found to involve significant energy-entropy compensation, pointing to the importance of developing water models that accurately describe the temperature dependence of water structure and dynamics.

Hydrogen Bond Dynamics and Phase Transitions of Water inside Carbon Nanotubes

Water dynamics in nanochannels are altered by confinement, particularly in small carbon nanotubes (CNTs). However, the mechanisms behind these effects remain unclear. To address these issues, we carried out extensive molecular dynamics (MD) simulations to investigate the structure and dynamics of water inside CNTs of different sizes (length of 20 nm and diameters vary from 0.8 nm to 5.0 nm) at different temperatures (from 200 K to 420 K). The radial density profile of water inside CNTs shows a single peak near the CNT walls for small nanotubes. For CNTs with larger sizes, water molecules are arranged into coaxial tubular sheets, the number of which increases with the CNT size. Subdiffusive behavior is observed for ultranarrow CNTs with diameters of 0.8 nm and 1 nm. As the size of CNTs increases, Fickian diffusion becomes evident. The hydrogen bond correlation function of water inside CNT decays slower than in bulk water, and the decay rate decreases as we increase the diameter of the CNTs. In large CNTs, the hydrogen bond lifetime of the innermost layer is shorter than the other layers and depends on temperature. Additional analysis of our results reveals that water molecules along the CNT axis show a non-Arrhenius to Arrhenius diffusion crossover. In general, the diffusion transition temperature is higher than that of bulk water, but it depends on the size of the CNT.

Ab Initio Molecular Dynamics Simulation of Water Transport through Short Carbon Nanotubes

Water transport through short single-walled (6, 6) carbon nanotubes (CNTs) was investigated with ab initio molecular dynamics (AIMD) simulation at different temperatures. The water molecules under extreme confinement present a one-dimensional jagged pattern owing to hydrogen bonding, with the near-perfect alignment of the dipole orientations. CNTs ending with dangling bonds can promote water dissociation near the entrance and the occurrence of dipole flipping along the water wire at high temperatures, accompanied by the formation of D defects and L defects in the hydrogen-bond network. In contrast, dissociation of water molecules rarely takes place if the dangling bonds at the ends of the CNTs are terminated with H atoms. Angular jumps of water molecules are commonplace inside the narrow CNTs, implying a low-energy barrier for hydrogen-bond exchange among water molecules in narrow CNTs. The simulation results demonstrate the high activity of dangling bonds at the ends of short CNTs, accompanying passivation processes and their profound impact on water structure and transport, which is important for diverse technological applications.

Dielectric Profile and Electromelting of a Monolayer of Water Confined in Graphene Slit Pore

A monolayer of water confined between two parallel graphene sheets exists in many different phases and exhibits fascinating dielectric properties that have been studied in experiments. In this work, we use molecular dynamics simulations to study how the dielectric properties of a confined monolayer of water is affected by its structure. We consider six of the popular nonpolarizable water models-SPC/E, SPC/Fw, TIP3P, TIP3P_M (modified), TIP4P-2005, and TIP4P-2005f-and find that the in-plane structure of the water molecules at ambient temperature and pressure is strongly dependent on the water model: all the 3-point water models considered here show square ice formation, whereas no such structural ordering is observed for the 4-point water models. This allows us to investigate the role of the in-plane structure of the water monolayer on its dielectric profile. Our simulations show an anomalous perpendicular dielectric constant compared to the bulk, and the models that do not exhibit ice formation show very different dielectric response along the channel width compared to models that exhibit square ice formation. We also demonstrate the occurrence of electromelting of the in-plane ordered water under the application of a perpendicular electric field and find that the critical field for electromelting strongly depends on the water model. Together, we have shown the dependence of confined water properties on the different water structures that it may take when sandwiched between bilayer graphene. These remarkable properties of confined water can be exploited in various nanofluidic devices, artificial ion channels, and molecular sieving.

Transport of hydrated nitrate and nitrite ions through graphene nanopores in aqueous medium

Nitrate ( NO 3 - ) and nitrite ( NO 2 - ) ions are naturally occurring inorganic ions that are part of the nitrogen cycle. High doses of these ions in drinking water impose a potential risk to public health. In this work, molecular dynamics simulations are carried out to study the passage of nitrate and nitrite ions from water through graphene nanosheets (GNS) with hydrogen-functionalized narrow pores in presence of an external electric field. The passage of ions through the pores is investigated through calculations of ion flux, and the results are analyzed through calculations of various structural and thermodynamic properties such as the density of ions and water, ion-water radial distribution functions, two-dimensional density distribution functions, and the potentials of mean force of the ions. Current simulations show that the nitrite ions can pass more in numbers than the nitrate ions in a given time through GNS hydrogen-functionalized pore of different geometry. It is found that the nitrite ions can permeate faster than the nitrate ions despite the former having higher hydration energy in the bulk. This can be explained in terms of the competition between the number density of the ions along the pore axis and the free energy barrier calculated from the potential of mean force. Also, an externally applied electric field is found to be important for faster permeation of the nitrite over the nitrate ions. The current study suggests that graphene nanosheets with carefully created pores can be effective in achieving selective passage of ions from aqueous solutions.

Enhanced wettability of long narrow carbon nanotubes in a double-walled hetero-structure: unraveling the effects of a boron nitride nanotube as the exterior

Studying the structure and dynamics of nano-confined water inside carbon nanotubes has consistently attracted the wide-spread interest of researchers. In the present work, molecular dynamics simulations indicated internal nonwetting behavior for the central region of the long and narrow single-wall carbon nanotube (5,5) (SWNT) and showed that continuous single-file water molecules are not formed through it. Unlike the SWNT, by adding boron nitride nanotubes (6,6) as an outer wall to the SWNT, a continuously long single-file water chain is formed through the double-walled carbon and boron nitride hetero-nanotube (DWHNT) and thorough internal wetting of the DWHNT is observed. The position and the number of water molecules, electrostatic potential heatmap of the nanotube's wall, free energy profile of nano-confined water, and number of hydrogen bonds between them confirmed the aforementioned results and complete internal wetting of the DWHNT. After using the boron nitride nanotube (6,6) as the outer wall, an homogeneous electrostatic potential distribution in the DWHNT and increase in the hydrophilic characteristics of the nano-channel wall are observed, bringing about gradual trapping of more water molecules through it. Finally, water molecules occupied the central region of the DWHNT and a thorough single-file water chain is formed inside the nano-channel. Water dipole orientation inside the DWHNT and their radial distribution function asserted the occurrence of the liquid-solid quasi-phase transition of single-file water molecules confined inside the long and narrow carbon nanotube (5,5) under ambient conditions.

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.

Structure and dynamic properties of stretched water in graphene nanochannels by molecular dynamics simulation: effects of stretching extent

Water confined in nanochannels can be stretched with the variation of external pressure, leading to the more disordered microstructure and higher diffusion coefficient than bulk water.

VP, Sathian SP. Water flow in carbon nanotubes: the role of tube chirality

Flow rate of water in CNTs of different types.

Water desalination using graphene nanopores: influence of the water models used in simulations

Water desalination using graphene nanopores was studied using different water models. The water permeation was found to be influenced by the bulk transport properties and the hydrogen-bond dynamics of the simulated water.

Molecular dynamics simulation strategies for designing carbon-nanotube-based targeted drug delivery

The carbon nanotube (CNT)-based target-specific delivery of drugs, or other molecular cargo, has emerged as one of the most promising biomedical applications of nanotechnology. To achieve efficient CNT-based drug delivery, the interactions between the drug, CNT and biomolecular target need to be properly optimized. Recent advances in the computer-aided molecular design tools, in particular molecular dynamics (MD) simulation studies, offer an appropriate low-cost approach for such optimization. This review highlights the various potential MD approaches for the simulation of CNT interactions with cell membranes while emphasizing various methods of cellular internalization and toxicities of CNTs to build new strategies for designing rational CNT-based targeted drug delivery to circumvent the limitations associated with the various clinically available nonspecific therapeutic agents.

DNA-Assisted Dispersion of Carbon Nanotubes and Comparison with Other Dispersing Agents

Separation and sorting of pristine carbon nanotubes (CNTs) from bundle geometry is a very challenging task due to the insoluble and nondispersive nature of CNTs in aqueous medium. Recently, many studies have been performed to address this problem using various organic and inorganic solutions, surfactant molecules, and biomolecules as dispersing agents. Recent experimental studies have reported the DNA to be highly efficient in dispersing CNTs from bundle geometry. However, there is no microscopic study and also quantitative estimation of the dispersion efficiency of the DNA. Using all-atom molecular dynamics simulation, we study the structure and stability of single-stranded DNA (ssDNA)-single-walled carbon nanotube (SWNT) (6,5) complex. To quantify the dispersion efficiency of various DNA sequences, we perform potential of mean forces (PMF) calculation between two bare SWNTs as well ssDNA-wrapped CNTs for different base sequences. From the PMF calculation, we find the PMF between two bare (6,5) SWNTs to be approximately -29 kcal/mol. For the ssDNA-wrapped SWNTs, the PMF reduces significantly and becomes repulsive. In the presence of ssDNA of different polynucleotide bases (A, T, G, and C), we present a microscopic picture of the ssDNA-SWNT (6,5) complex and also a quantitative estimate of the interaction strength between nanotubes from PMF calculation. From PMF, we show the sequence of dispersion efficiency for four different nucleic bases to be T > A > C > G. We have also presented a comparison of the dispersion efficiencies of ssDNA, flavin mononucleotide surfactant, and poly(amidoamine) (PAMAM) dendrimer by comparing their respective PMF values.

Confined Water: Structure, Dynamics, and Thermodynamics

Understanding the properties of strongly confined water is important for a variety of applications such as fast flow and desalination devices, voltage generation, flow sensing, and nanofluidics. Confined water also plays an important role in many biological processes such as flow through ion channels. Water in the bulk exhibits many unusual properties that arise primarily from the presence of a network of hydrogen bonds. Strong confinement in structures such as carbon nanotubes (CNTs) substantially modifies the structural, thermodynamic, and dynamic (both translational and orientational) properties of water by changing the structure of the hydrogen bond network. In this Account, we provide an overview of the behavior of water molecules confined inside CNTs and slit pores between graphene and graphene oxide (GO) sheets. Water molecules confined in narrow CNTs are arranged in a single file and exhibit solidlike ordering at room temperature due to strong hydrogen bonding between nearest-neighbor molecules. Although molecules constrained to move along a line are expected to exhibit single-file diffusion in contrast to normal Fickian diffusion, we show, from a combination of molecular dynamics simulations and analytic calculations, that water molecules confined in short and narrow CNTs with open ends exhibit Fickian diffusion because of their collective motion as a single unit due to strong hydrogen bonding. Confinement leads to strong anisotropy in the orientational relaxation of water molecules. The time scale of relaxation of the dipolar correlations of water molecules arranged in a single file becomes ultraslow, of the order of several nanoseconds, compared with the value of 2.5 ps for bulk water. In contrast, the relaxation of the vector that joins the two hydrogens in a water molecule is much faster, with a time scale of about 150 fs, which is about 10 times shorter than the corresponding time scale for bulk water. This is a rare example of confinement leading to a speedup of orientational dynamics. The orientational relaxation of confined water molecules proceeds by angular jumps between two locally stable states, making the relaxation qualitatively different from that expected in the diffusive limit. The spontaneous entry of water inside the hydrophobic cavity of CNTs is primarily driven by an increase in the rotational entropy of water molecules inside the cavity, arising from a reduction in the average number of hydrogen bonds attached to a water molecule. From simulations using a variety of water models, we demonstrate that the relatively simple SPC/E water model yields results in close agreement with those obtained from polarizable water models. Finally, we provide an account of the structure and thermodynamics of water confined in the slit pore between two GO sheets with both oxidized and reduced parts. We show that the potential of mean force for the oxidized part of GO sheets in the presence of water exhibits two local minima, one corresponding to a dry cavity and the other corresponding to a fully hydrated cavity. The coexistence of these two regimes provides permeation pathways for water in GO membranes.

Dynamics and structural behavior of water in large confinement with planar amorphous walls

We study the structure and dynamics of liquid water confined between planar amorphous walls using molecular dynamics (MD) simulations. We report MD results for systems of more than 23 000 SPC/E water molecules confined between two hydrophilic or hydrophobic walls, separated by distances of about 15 nm. We find that the walls induce ordering of the liquid and slow down the dynamics, affecting the properties of the confined water up to distances of about 8 nm at 275 K. We quantify this influence by computing dynamic and static penetration lengths and studying their temperature dependence. Our results indicate that in the temperature range considered, hydrophobic walls perturb static properties over larger lengths compared to hydrophilic walls. We also find opposite temperature trends in the dynamic penetration lengths, with hydrophobic walls increasing their range of influence on increasing the temperature.

Surface dissimilarity affects critical distance of influence for confined water

In this study, the properties of nano-confined water, such as density, orientation etc., are monitored across varying confinement spacing to determine the critical distance of influence between dissimilar surfaces.

Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Using a molecular model for octamethylcyclotetrasiloxane (OMCTS), molecular dynamics simulations are carried out to probe the phase state of OMCTS confined between two mica surfaces in equilibrium with a reservoir. Molecular dynamics simulations are carried out for elevations ranging from 5 to 35 K above the melting point for the OMCTS model used in this study. The Helmholtz free energy is computed for a specific confinement using the two-phase thermodynamic (2PT) method. Analysis of the in-plane pair correlation functions did not reveal signatures of freezing even under an extreme confinement of two layers. OMCTS is found to orient with a wide distribution of orientations with respect to the mica surface, with a distinct preference for the surface parallel configuration in the contact layers. The self-intermediate scattering function is found to decay with increasing relaxation times as the surface separation is decreased, and the two-step relaxation in the scattering function, a signature of glassy dynamics, distinctly evolves as the temperature is lowered. However, even at 5 K above the melting point, we did not observe a freezing transition and the self-intermediate scattering functions relax within 200 ps for the seven-layered confined system. The self-diffusivity and relaxation times obtained from the Kohlrausch-Williams-Watts stretched exponential fits to the late α-relaxation exhibit power law scalings with the packing fraction as predicted by mode coupling theory. A distinct discontinuity in the Helmholtz free energy, potential energy, and a sharp change in the local bond order parameter, Q4, was observed at 230 K for a five-layered system upon cooling, indicative of a first-order transition. A freezing point depression of about 30 K was observed for this five-layered confined system, and at the lower temperatures, contact layers were found to be disordered with long-range order present only in the inner layers. These dynamical signatures indicate that confined OMCTS undergoes a slowdown akin to a fluid approaching a glass transition upon increasing confinement, and freezing under confinement would require substantial subcooling below the bulk melting point of OMCTS.

Molecular mechanism of water permeation in a helium impermeable graphene and graphene oxide membrane

Layers of graphene oxide (GO) are found to be good for the permeation of water but not for helium (Science, 2012, 335(6067), 442-444) suggesting that the GO layers are dynamic in the formation of a permeation route depending on the environment they are in (i.e., water or helium). To probe the microscopic origin of this observation we calculate the potential of mean force (PMF) of GO sheets (with oxidized and reduced parts), with the inter-planar distance as a reaction coordinate in helium and water. Our PMF calculation shows that the equilibrium interlayer distance between the oxidized part of the GO sheets in helium is at 4.8 Å leaving no space for helium permeation. In contrast, the PMF of the oxidized part of the GO in water shows two minima, one at 4.8 Å and another at 6.8 Å, corresponding to no water and a water filled region, thus giving rise to a permeation path. The increased electrostatic interaction between water with the oxidized part of the sheet helps the sheet open up and pushes water inside. Based on the entropy calculations for water trapped between graphene sheets and oxidized graphene sheets at different inter-sheet spacings, we also show the thermodynamics of filling.