Flow of water through carbon nanotubes predicted by different atomistic water models

CONAN─Novel Tool to Create and Analyze Liquids in Confined Space

Modeling of complex liquids at solid surfaces and in confinement is gaining attention due to an increase in computer power and advancement of simulation techniques. Therefore, tools to set up structures and for analysis are needed. In this paper, we present CONAN─a Python code designed to facilitate the study of liquids interacting with solid structures, such as walls or pores. Among other things, the program provides the option to generate a variety of different structures, including carbon walls and nanotubes and their boron nitride analogs, as well as the ability to analyze various structural properties of confined and interfacial liquids. In the case of the ionic liquid 1-butyl-3-methylimidazolium acetate in carbon nanotubes of different sizes, we demonstrate the abilities of our tool. The average density within the confinement highly depends on the carbon nanotube size, and it is generally lower than the density of the bulk liquid. The arrangement of the individual species within the tube also depends on size, with radial layers forming within the tubular confinement. The density is largely increased in the respective layers, while it is drastically reduced between the layers.

Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models

In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.

Functionalized carbon nanocones performance in water harvesting

In this work, we investigate the water capture process for functionalized carbon nanocones (CNCs) through molecular dynamic simulations in the following three scenarios: a single CNC in contact with a reservoir containing liquid water, a single CNC in contact with a water vapor reservoir, and a combination of more than one CNC in contact with vapor. We found that water flows through the nanocones when in contact with the liquid reservoir if the nanocone tip presents hydrophilic functionalization. In contact with steam, we observed the formation of droplets at the base of the nanocone only when hydrophilic functionalization is present. Then, water flows through in a linear manner, a process that is more efficient than that in the liquid reservoir regime. The scalability of the process is tested by analyzing the water flow through more than one nanocone. The results suggest that the distance between the nanocones is a fundamental ingredient for the efficiency of water harvesting.

Horizontally Asymmetric Nanochannels of Graphene Oxide Membranes for Efficient Osmotic Energy Harvesting

Reverse electrodialysis (RED) directly harvests renewable energy from salinity gradients, and the achievable potential power heavily relies on the ion exchange membranes. Graphene oxides (GOs) are considered a solid candidate for the RED membrane because the laminated GO nanochannels with charged functional groups provide an excellent ionic selectivity and conductivity. Yet, a high internal resistance and poor stability in aqueous solutions limit the RED performance. Here, we develop a RED membrane that concurrently achieves high ion permeability and stable operation based on epoxy-confined GO nanochannels with asymmetric structures. The membrane is fabricated by reacting epoxy-wrapped GO membranes with ethylene diamine via vapor diffusion, overcoming the swelling properties in aqueous solutions. More importantly, the resultant membrane exhibits asymmetric GO nanochannels in terms of both channel geometry and electrostatic surface charges, leading to the rectified ion transport behavior. The demonstrated GO membrane exhibits the RED performance up to 5.32 W·m^-2 with >40% energy conversion efficiency across a 50-fold salinity gradient and 20.3 W·m^-2 across a 500-fold salinity gradient. Planck-Nernst continuum models coupled to molecular dynamics simulations rationalize the improved RED performance in terms of the asymmetric ionic concentration gradient within the GO nanochannel and the ionic resistance. The multiscale model also provides the design guidelines for ionic diode-type membranes configuring the optimum surface charge density and ionic diffusivity for efficient osmotic energy harvesting. The synthesized asymmetric nanochannels and their RED performance demonstrate the nanoscale tailoring of the membrane properties, establishing the potentials for 2D material-based asymmetric membranes.

Atmospheric water harvesting using functionalized carbon nanocones

In this work, we propose a method to harvest liquid water from water vapor using carbon nanocones. The condensation occurs due to the presence of hydrophilic sites at the nanocone entrance. The functionalization, together with the high mobility of water inside nanostructures, leads to a fast water flow through the nanostructure. We show using molecular dynamics simulations that this device is able to collect water if the surface functionalization is properly selected.

Coupling solute interactions with functionalized graphene membranes: towards facile membrane-level engineering

Optimizing ion transport through nanoporous graphene membranes with intricate engineering at nanoscale levels finds applications ranging from ion segregation to desalination. Such membrane-level engineering often requires futuristic and state-of-the-art micro- and nanofabrication infrastructure making it less accessible to widespread applications. In this study, the effective membrane pore size is modulated using macroscopic membrane functionalization, which, when combined with the solute concentration, can prove to be facile nanoscale engineering towards achieving selectivity. By performing robust molecular dynamics (MD) simulations of aqueous NaCl solution through a nanoporous graphene membrane, we demonstrate that varying membrane wettability influences the structural organization of ions and water molecules both in the vicinity and inside the nanopore, which is manifested in the form of altered permeation characteristics. Moreover, the disparate solvation characteristics of the ionic species in conjunction with the variable van der Waals interactive forces affect the ion-selective nature (Cl^- over Na⁺) of the membrane. The relative hydrophilization, resulting from the effective functionalization of the nanoporous graphene membrane, not only allows greater control over the permeation characteristics of ions and water molecules mediated by an altered depletion ratio but also gives rise to the ion-selective nature of the membrane, thus providing a sound understanding of the transport properties of ion-water solutions through nanoporous materials.

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.

Flow-Induced Autonomic Ordering of Hydrogen Molecules under a Non-Equilibrium Flow

A non-equilibrium molecular flow through a carbon nanotube (CNT) serves as a key system for revealing molecular transport and establishing nanofluidics. It has been challenging to simulate a non-equilibrium flow of hydrogen molecules exhibiting strong nuclear quantumness. Taking advantage of the quantum molecular dynamics method that can calculate real-time trajectories of hydrogen molecules even under a non-equilibrium flow, we found that the non-equilibrium flow makes hydrogen molecules more condensed and accelerates their adsorption near a CNT surface, letting the molecules flow more smoothly by propagating velocity momenta more efficiently along the CNT axis and by suppressing transverse molecular dynamics on the CNT cross section. Such flow-induced autonomic ordering indicates the importance of monitoring and investigating dynamics and adsorption of hydrogen molecules under a non-equilibrium circumstance as well as in a quiet equilibrium state, opening a new strategy for efficient hydrogen liquefaction and storage.

From a bulk to nanoconfined water chain: bridge water at the pore of the (6,6) carbon nanotube

Hydrophobic porous materials with nano-pores are critical in many processes such as water desalination and biological membrane transportation. Herein, we performed molecular dynamics (MD) simulations on a prototypical hydrophobic nanochannel consisting of a (6,6) carbon nanotube (CNT) of 4.12 Å in radius and 13.72 Å in length immersed in water. The simulation shows that there are two major filling numbers of water N = 5 and N = 6, with the former being the most stable one. The confined waters form a single-file water chain with two hydrogen bonds per water. An extending water chain is formed for N = 5, with a bridge water near the pore of the CNT linking the water confined inside the CNT and hydration layer around the pore of the CNT. The bridge water can be considered as intermediate water characterized by three hydrogen bonds that distinguish from the confined water and bulk water. On the other hand, the hydration layer is depleted from the pore when N = 6. The analyses of the correlation of the bond order for the adjacent hydrogen bond pair of the hydration layer around the pore of the CNT does not show apparent difference from that of bulk water, though the former is slightly ordered. van Hove analysis of the bridge water shows that it tends to move inside the CNT when N < 5, in order to maintain the chemical equilibrium between the confined water and bulk water. This study highlights the unique structure of water around the hydrophobic pore of a sub-nanometer nanochannel.

Formation of a Stable Bridge between Two Disjoint Nanotubes with Single-File Chains of Water

It was recently demonstrated that stable water bridges can form between two relatively large disjoint nanochannels, such as carbon nanotubes (CNTs), under an applied pressure drop. Such bridges are relevant to fabrication of nanostructured materials, drug delivery, water desalination devices, hydrogen fuel cells, dip-pen nanolithography, and several other applications. If the nanotubes are small enough, however, then one has only single-file hydrogen-bonded chains of water molecules. The distribution of water in such nanotubes manifests unusual physical properties that are attributed to the low number of hydrogen bonds (HBs) formed in the channel since, on average, each water molecule in a single-file chain forms only 1.7 HBs, almost half of the value for bulk water. Using extensive molecular dynamics simulations, we demonstrate that stable bridges can form even between two small disjoint CNTs that contain single-file chains of water. The structure, stability, and properties of such bridges and their dependence on the applied pressure drop and the length of the gap between the two CNTs are studied in detail, as is the distribution of the HBs. We demonstrate, in particular, that the efficiency of flow through the bridge is at maximum at a specific pressure difference.

Metallicity-Dependent Ultrafast Water Transport in Carbon Nanotubes

Carbon nanotubes (CNTs) with hydrophobic and atomically smooth inner channels are promising for building ultrahigh-flux nanofluidic platforms for energy harvesting, health monitoring, and water purification. Conventional wisdom is that nanoconfinement effects determine water transport in CNTs. Here, using full-atomistic molecular dynamics simulations, it is shown that water transport behavior in CNTs strongly correlates with the electronic properties of single-walled CNTs (metallic (met) vs semiconducting (s/c)), which is as dominant as the effect of nanoconfinement. Three pairs of CNTs (i.e., (8,8)_met , 10.85 Å vs (9,7)_s/c , 10.88 Å; (9,8)_s/c , 11.53 Å vs (10,7)_met , 11.59 Å; and (9,9)_met , 12.20 Å vs (10,8)_s/c , 12.23 Å) are used to investigate the roles of diameter and metallicity. Specifically, the (9,8)_s/c can restrict the hydrogen-bonding-mediated structuring of water and give the highest reduction in carbon-water interaction energy, providing an extraordinarily high water flux, around 250 times that of the commercial reverse osmosis membranes and approximately fourfold higher than the flux of the state-of-the-art boron nitrate nanotubes. Further, the high performance of (9,8)_s/c is also reproducible when embedded in lipid bilayers as synthetic high-water flux porins. Given the increasing availability of high-purity CNTs, these findings provide valuable guides for realizing novel CNT-enhanced nanofluidic systems.

Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS2 nanopores

Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based on low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based on nanoporous materials have been suggested in recent years. One of the materials suggested is MoS₂, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct force fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were built based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, and the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through nanopores of two sizes were used - one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter - and two rigid water models from the TIP4P family - TIP4P/2005 and TIP4P/ε. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only can the water model influence the outcomes, but the ion model plays a crucial role when the pore is small enough.