Water in carbon nanotubes: adsorption isotherms and thermodynamic properties from molecular simulation

Grand canonical Monte Carlo simulations are performed to study the adsorption of water in single-walled (6:6), (8:8), (10:10), (12:12), and (20:20) carbon nanotubes in the 248-548 K temperature range. At room temperature the resulting adsorption isotherms in (10:10) and wider single-walled carbon nanotubes (SWCNs) are characterized by negligible water uptake at low pressures, sudden and complete pore filling once a threshold pressure is reached, and wide adsorption/desorption hysteresis loops. The width of the hysteresis loops decreases as pore diameter narrows and it becomes negligible for water adsorption in (8:8) and (6:6) SWCNs. Results for the isosteric heat of adsorption, density profiles along the pore axis and across the pore radii, order parameter across the pore radii, and x-ray diffraction patterns are presented. Layered structures are observed when the internal diameter of the nanotubes is commensurate to the establishment of a hydrogen-bonded network. The structure of water in (8:8) and (10:10) SWCNs is ordered when the temperature is 298 and 248 K, respectively. By simulating adsorption isotherms at various temperatures, the hysteresis critical temperature, e.g., the lowest temperature at which no hysteresis can be detected, is determined for water adsorbed in (20:20), (12:12), and (10:10) SWCNs. The hysteresis critical temperature is lower than the vapor-liquid critical temperature for bulk Simple Point Charge-Extended (SPCE) water model.     

Influence of surface roughness on cetyltrimethylammonium bromide adsorption from aqueous solution

The influence of surface roughness on surfactant adsorption was studied using a quartz crystal microbalance with dissipation (QCM-D). The sensors employed had root-mean-square (R) roughness values of 2.3, 3.1, and 5.8 nm, corresponding to fractal-calculated surface area ratios (actual/nominal) of 1.13, 1.73, and 2.53, respectively. Adsorption isotherms measured at 25 °C showed that adsorbed mass of cetyltrimethylammonium bromide per unit of actual surface area below 0.8 cmc, or above 1.2 cmc, decreases as the surface roughness increases. At the cmc, both the measured adsorbed amount and the measured dissipation increased dramatically on the rougher surfaces. These results are consistent with the presence of impurities, suggesting that roughness exacerbates well-known phenomena reported in the literature of peak impurity-related adsorption at the cmc. The magnitude of the increase, especially in dissipation, suggests that changes in adsorbed amount may not be the only reason for the observed results, as aggregates at the cmc on rougher surfaces are more flexible and likely contain larger amounts of solvent. Differences in adsorption kinetics were also found as a function of surface roughness, with data showing a second, slower adsorption rate after rapid initial adsorption. A two-rate Langmuir model was used to further examine this effect. Although adsorption completes faster on the smoother surfaces, initial adsorption at zero surface coverage is faster on the rougher surfaces, suggesting the presence of more high-energy sites on the rougher surfaces.     

Aggregation of POSS monomers in liquid hexane: a molecular-simulation study

Polyhedral oligomeric silsesquioxanes (POSS) are multifunctional molecules that can be employed as building blocks to develop nanocomposite materials whose mechanical properties often improve upon those of traditional polymeric materials. We report here molecular simulation results for the effective potential of mean force between octamethyl POSS monomers and between POSS monomers in which one methyl group has been substituted by a linear alkane chain of nine carbon atoms in liquid normal hexane at 300 and 400 K. The results are discussed and compared to available data for the effective interactions between octamethyl POSS monomers in normal hexadecane. Our results show that the effective short-ranged POSS-POSS attraction is significantly weaker in hexane than it is in hexadecane, perhaps explaining why normal hexane is often the solvent of choice for the preparation of POSS-containing materials. Additionally, we provide results for the radial distribution functions between selected sites in the POSS monomers that can be used both to understand the association between POSS monomers in solution and to parametrize coarse-grained simulation models. Such models will be used to study the formation of POSS-containing supramolecular structures such as lamellae or micelles that are currently not accessible by atomistic simulation and can be compared to experimental observations.     

Simulated water adsorption in chemically heterogeneous carbon nanotubes

Grand canonical Monte Carlo simulations are used to study the adsorption of water in single-walled (10:10), (12:12), and (20:20) carbon nanotubes at 298 K. Water is represented by the extended simple point charge model and the carbon atoms as Lennard-Jones spheres. The nanotubes are decorated with different amounts of oxygenated sites, represented as carbonyl groups. In the absence of carbonyl groups the simulated isotherms are characterized by negligible amounts of water uptake at low pressures, sudden and complete pore filling once a threshold pressure is reached, and wide adsorption-desorption hysteresis loops. In the presence of a few carbonyl groups the simulated adsorption isotherms are characterized by pore filling at lower pressures and by narrower adsorption-desorption hysteresis loops compared to the results obtained in the absence of carbonyl groups. Our results show that the distribution of the carbonyl groups has a strong effect on the adsorption isotherms. For carbonyl groups localized in a narrow section the adsorption of water may be gradual because a cluster of adsorbed water forms at low pressures and grows as the pressure increases. For carbonyl groups distributed along the nanotube the adsorption isotherm is of type V.     

Aqueous NaCl and CsCl solutions confined in crystalline slit-shaped silica nanopores of varying degree of protonation

All-atom molecular dynamics simulations were conducted to study the dynamics of aqueous electrolyte solutions confined in slit-shaped silica nanopores of various degrees of protonation. Five degrees of protonation were prepared by randomly removing surface hydrogen atoms from fully protonated crystalline silica surfaces. Aqueous electrolyte solutions containing NaCl or CsCl salt were simulated at ambient conditions. In all cases, the ionic concentration was 1 M. The results were quantified in terms of atomic density distributions within the pores, and the self-diffusion coefficient along the direction parallel to the pore surface. We found evidence for ion-specific properties that depend on ion-surface, water-ion, and only in some cases ion-ion correlations. The degree of protonation strongly affects the structure, distribution, and the dynamic behavior of confined water and electrolytes. Cl(-) ions adsorb on the surface at large degrees of protonation, and their behavior does not depend significantly on the cation type (either Na(+) or Cs(+) ions are present in the systems considered). The cations show significant ion-specific behavior. Na(+) ions occupy different positions within the pore as the degree of protonation changes, while Cs(+) ions mainly remain near the pore center at all conditions considered. For a given degree of protonation, the planar self-diffusion coefficient of Cs(+) is always greater than that of Na(+) ions. The results are useful for better understanding transport under confinement, including brine behavior in the subsurface, with important applications such as environmental remediation.     

Adsorption of comb copolymers on weakly attractive solid surfaces

In this work continuum and lattice Monte Carlo simulation methods are used to study the adsorption of linear and comb polymers on flat surfaces. Selected polymer segments, located at the tips of the side chains in comb polymers or equally spaced along the linear polymers, are attracted to each other and to the surface via square-well potentials. The rest of the polymer segments are modeled as tangent hard spheres in the continuum model and as self-avoiding random walks in the lattice model. Results are presented in terms of segment-density profiles, distribution functions, and radii of gyration of the adsorbed polymers. At infinite dilution the presence of short side chains promotes the adsorption of polymers favoring both a decrease in the depletion-layer thickness and a spreading of the polymer molecule on the surface. The presence of long side chains favors the adsorption of polymers on the surface, but does not permit the spreading of the polymers. At finite concentration linear polymers and comb polymers with long side chains readily adsorb on the solid surface, while comb polymers with short side chains are unlikely to adsorb. The simple models of comb copolymers with short side chains used here show properties similar to those of associating polymers and of globular proteins in aqueous solutions, and can be used as a first approximation to investigate the mechanism of adsorption of proteins onto hydrophobic surfaces.     

Effect of temperature on the adsorption of water in porous carbons

We report experimental and simulation studies to investigate the effect of temperature on the adsorption isotherms for water in carbons. Adsorption isotherms are measured by a gravimetric technique in carbon-fiber monoliths at 378 and 423 K and studied by molecular simulation in ideal carbon pores in the temperature range 298-600 K. Experimental adsorption isotherms show a gradual water uptake, as the pressure increases, and narrow adsorption-desorption hysteresis loops. In contrast, simulated adsorption isotherms at room temperature are characterized by negligible uptake at low pressures, sudden and complete pore filling once a threshold pressure is reached, and wide adsorption-desorption hysteresis loops. As the temperature increases, the relative pressure at which pore filling occurs increases and the size of the hysteresis loop decreases. Experimental adsorption-desorption hysteresis loops are narrower than those from simulation. Discrepancies between simulation and experimental results are attributed to heterogeneities in chemical composition, pore connectivity, and nonuniform pore-size distribution, which are not accounted for in the simulation model. The hysteresis phase diagram for confined water is obtained by recording the pressure-density conditions that bound the simulated hysteresis loop at each temperature. We find that the hysteresis critical temperature, i.e., the lowest temperature at which no hysteresis is detected, can be hundreds of degrees lower than the vapor-liquid critical temperature for bulk model water. The properties of confined water are discussed with the aid of simulation snapshots and by analyzing the structure of the confined fluid.     

Interfacial water studies and their relevance for the energy sector

ABSTRACT

Interfacial water is ubiquitous, and its investigation has attracted the interest of many for both fundamental and applied purposes. This perspective provides a few highlights concerning how molecular and sometimes multi-scale investigations on the properties of interfacial water could be of practical relevance for the energy sector. The discussion focuses on the transport of electrolytes through narrow pores, and on the solubility and transport properties of confined hydrocarbons, with particular attention in attempting to describe how interfacial water can affect such properties, especially within narrow pores. Recent results on the self-assembly of surfactants on heterogeneous surfaces are also briefly discussed. Finally, a few possible directions for future research are suggested. It should be stressed that this is not a comprehensive review on the possible impact of interfacial water studies on the energy sector, but rather a brief commentary on some personal interests, with the intent of stimulating debate.     

Hydrogen-bond dynamics for water confined in carbon tetrachloride-acetone mixtures

In a variety of biological scenarios water is found trapped within hydrophobic environments (e.g., ion channels). Its behavior under such conditions is not well understood and therefore is attracting enormous scientific attention. It is of particular interest to understand how the confining environment affects both the structure and dynamics of water. Within this scenario, we report molecular dynamics simulation results for water trapped in a mixture of acetone and carbon tetrachloride whose composition mimics the one employed in recently reported experiments [Gilijamse, J. J.; et al. Proc. Natl. Acad. Sci. U.S.A. 2005, 102, 3202]. We show here that the water molecules dissolved in the carbon tetrachloride-acetone mixture assemble in clusters of varying sizes, that the longevity of hydrogen bonds between confined water molecules strongly depends on the cluster size, and that hydrogen bonds last longer for small water clusters in confined water than they do in bulk water. The simulated FT-IR spectra for the confined water are shifted at longer frequencies compared to those observed for bulk liquid water. We discuss the dependence of the FT-IR spectrum on the size of the water clusters dispersed in the carbon tetrachloride-acetone matrix. We also study in detail the rotational orientation of the dispersed water molecules, and we discuss how the composition of the organic matrix affects the results. By enhancing the interpretation of the experimental data, our results contribute to developing a molecular-based understanding of the relationship between environment and water properties.     

Steric stabilization of spherical colloidal particles: Implicit and explicit solvent

We present the results of Monte Carlo simulations and density functional theory treatment of interactions between spherical colloidal brushes both in implicit (good) solvent and in an explicit polymeric solution. Overall, theory is seen to be in good agreement with simulations. We find that interactions between hard-sphere particles grafted with hard-sphere chains are always repulsive in implicit solvent. The range and steepness of the repulsive interaction is sensitive to the grafting density and the length of the grafted chains. When the brushes are immersed in an explicit solvent of hard-sphere chains, a weak mid-range attraction arises, provided the length of the free chains exceeds that of the grafted chains.