Collective Proton Transfers in Cyclic Water-Ammonia Tetramers: A Path Integral Machine-Learning Study.

We present results from machine-learning-based path integral molecular dynamics simulations that describe isomerization paths articulated via collective proton transfers along mixed, cyclic tetramers combining water and ammonia at cryogenic conditions. The net result of such isomerizations is a reverse of the chirality of the global hydrogen-bonding architecture along the different cyclic moieties. In monocomponent tetramers, the classical free energy profiles associated with these isomerizations present the usual symmetric double-well characteristics whereas the reactive paths exhibit full concertedness among the different intermolecular transfer processes. Contrastingly, in mixed water/ammonia tetramers, the incorporation of a second component introduces imbalances in the strengths of the different hydrogen bonds leading to a partial loss of concertedness, most notably at the vicinity of the transition state. As such, the highest and lowest degrees of progression are registered along OH···N and O···HN coordinations, respectively. These characteristics lead to polarized transition state scenarios akin to solvent-separated ion-pair configurations. The explicit incorporation of nuclear quantum effects promotes drastic depletions in the activation free energies and modifications in the overall shape of the profiles which include central plateau-like stages, indicating the prevalence of deep tunneling regimes. On the other hand, the quantum treatment of the nuclei partially restores the degree of concertedness among the evolutions of the individual transfers.

Equilibrium and Dynamical Characteristics of Hydrogen Bond Bifurcations in Water-Water and Water-Ammonia Dimers: A Path Integral Molecular Dynamics Study.

We present path integral molecular dynamics results that describe the effects of nuclear quantum fluctuations on equilibrium and dynamical characteristics pertaining to bifurcation pathways in hydrogen bonded dimers combining water and ammonia, at cryogenic temperatures of the order of 20 K. Along these isomerizations, the hydrogen atoms in the molecules acting as hydrogen-bond donors interchange their original dangling/connective characters. Our results reveal that the resulting quantum transition paths comprise three stages: the initial and final ones involve overall rotations during which the two protons retain their classical-like characteristics. Effects from quantum fluctuation are clearly manifested in the changes operated at the intermediate passages over transition states, as the spatial extents of the protons stretch over typical lengths comparable to the distances between connective and dangling basins of attractions. Consequently, the classical over-the-hill path is replaced by a tunneling controlled mechanism which, within the path integral perspective, can be cast in terms of concerted inter-basin migrations of polymer beads from dangling-to-connective and from connective-to-dangling, at practically no energy costs. We also estimated the characteristic timescales describing such interconversions within the approximate ring polymer rate theory. Effects derived from full and partial deuteration are also discussed.

Structure and dynamics of nanoconfined water and aqueous solutions.

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Nuclear quantum effects on the hydrogen bond donor-acceptor exchange in water-water and water-methanol dimers.

We present results from path integral molecular dynamics simulations that describe effects from the explicit incorporation of nuclear quantum fluctuations on the topology of the free energy associated with the geared exchange of hydrogen bonds in the water-water dimer. Compared to the classical treatment, our results reveal important reductions in the free energy barriers and changes at a qualitative level in the overall profile. Most notable are those manifested by a plateau behavior, ascribed to nuclear tunneling, which bridges reactant and product states, contrasting with the usual symmetric double-well profile. The characteristics of the proton localizations along the pathway are examined. An imaginary time analysis of the rotational degrees of freedom of the partners in the dimer at the vicinities of transition states shows a clear "anticorrelation" between intermolecular interactions coupling beads localized in connective and dangling basins of attractions. As such, the transfer is operated by gradual concerted inter-basin migrations in opposite directions, at practically no energy costs. Modifications operated by partial deuteration and by the asymmetries in the hydrogen bonding characteristics prevailing in water-methanol heterodimers are also examined.

Combined Effects from Solvation and Nuclear Quantum Fluctuations on Autoionization Mechanisms in Aqueous Clusters.

Using path-integral molecular dynamics simulations, we examine isomerization paths involving collective proton transfers in [H2O]5 and [H2O]8 clusters under cryogenic conditions. We focused attention on combined effects derived from solvation and nuclear quantum fluctuations on the characteristics of free energy barriers and relative stabilities of reactants and products. In particular, we analyzed two different processes: the first one involves the exchange of donor-acceptor hydrogen bond roles along cyclic moieties, whereas the second one corresponds to charge separation leading to stable [H3O]+[OH]- ion pairs. In the first case, the explicit incorporation of quantum tunneling introduces important modifications in the classical free energy profile. The resulting quantum profile presents two main contributions, one corresponding to compressions of O-O distances and a second one ascribed to nuclear tunneling of the light protons. Solvation effects promote a moderate polarization of the cyclic structures and a partial loss of concertedness in the collective modes, most notably, at the onset of tunneling. Still, the latter effects are also sufficiently strong to promote the stabilization of ion pairs along the classical trajectories. In contrast, the explicit incorporation of nuclear quantum fluctuations leads to charge separated configurations that are marginally stable. As such, the latter states could also be regarded as short-lived intermediate states along the reactive exchange path.

Isotopic equilibria in aqueous clusters at low temperatures: Insights from the MB-pol many-body potential.

By combining path-integrals molecular dynamics simulations with the accurate MB-pol potential energy surface, we investigate the role of alternative potential models on isotopic fractionation ratios between H and D atoms at dangling positions in water clusters at low temperatures. Our results show clear stabilizations of the lighter isotope at dangling sites, characterized by free energy differences ΔG that become comparable to or larger than kBT for temperatures below ∼75 K. The comparison between these results to those previously reported using the empirical q-TIP4P/F water model [P. E. Videla et al., J. Phys. Chem. Lett. 5, 2375 (2014)] reveals that the latter Hamiltonian overestimates the H stabilization by ∼25%. Moreover, predictions from the MB-pol model are in much better agreement with measured results reported for similar isotope equilibria at ice surfaces. The dissection of the quantum kinetic energies into orthogonal directions shows that the dominant differences between the two models are to be found in the anharmonic characteristics of the potential energy surfaces along OH bond directions involved in hydrogen bonds.

Ionic Transport and Speciation of Lithium Salts in Glymes: Experimental and Theoretical Results for Electrolytes of Interest for Lithium-Air Batteries.

Glycol ethers, or glymes, have been recognized as good candidates as solvents for lithium-air batteries because they exhibit relatively good stability in the presence of superoxide radicals. Diglyme (bis(2-methoxy-ethyl)ether), in spite of its low donor number, has been found to promote the solution mechanism for the formation of Li2O2 during the discharge reaction, leading to large deposits, that is, high capacities. It has been suggested that lithium salt association in these types of solvents could be responsible for this behavior. Thus, the knowledge of the speciation and transport behavior of lithium salts in these types of solvents is relevant for the optimization of the lithium-air battery performance. In this work, a comprehensive study of lithium trifluoromethanesulfonate (LiTf) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in 1,2-di-methoxyethane (DME) and diglyme, over a wide range of concentrations, have been performed. Consistent ion pairs and triplet ions formation constants have been obtained by resorting to well-known equations that describe the concentration dependence of the molar conductivities in highly associated electrolytes, and we found that the system LiTf/DME would be the best to promote bulky Li2O2 deposits. Unexpected differences are observed for the association constants of LiTf and, to a lesser extent, for LiTFSI, in DME and diglyme, whose dielectric constants are similar. Molecular dynamics (MD) simulations allowed us to rationalize these differences in terms of the competing interactions of the O-sites of the ethers and the SO x groups of the corresponding anions with Li+ ion. The limiting Li+ diffusivity derived from the fractional Walden rule agrees quite well with those obtained from MD simulations, when solvent viscosity is conveniently rescaled.

Positional Isotope Exchange in HX·(H2O)n (X = F, I) Clusters at Low Temperatures.

We present molecular dynamics simulation results describing proton/deuteron exchange equilibria along hydrogen bonds at the vicinity of HX acids (X = F, I) in aqueous clusters at low temperatures. To allow for an adequate description of proton transfer processes, our simulation scheme resorted on the implementation of a multistate empirical valence bond hamiltonian coupled to a path integral scheme to account for effects derived from nuclear quantum fluctuations. We focused attention on clusters comprising a number of water molecules close to the threshold values necessary to stabilize contact-ion-pairs. For X = F, our results reveal a clear propensity of the heavy isotope to lie at the bond bridging the halide to the nearest water molecule. Contrasting, for X = I, the thermodynamic stability is reversed and the former connectivity is preferentially articulated via the light isotope. These trends remain valid for undissociated and ionic descriptions of the stable valence bond states. The preferences are rationalized in terms of differences in the quantum kinetic energies of the isotopes which, in turn, reflect the extent of the local spatial confinements prevailing along the different hydrogen bonds in the clusters. In most cases, these features are also clearly reflected in the characteristics of the corresponding stretching bands of the simulated infrared spectra. This opens interesting possibilities to gauge the extent of the isotopic thermodynamic stabilizations and the strengths of the different hydrogen bonds by following the magnitudes and shifts of the spectral signals in temperature-controlled experiments, performed on mixed clusters combining H2O and HOD.

Isotopic Preferential Solvation of I(-) in Low-Temperature Water Nanoclusters.

We present results from ring polymer molecular dynamics experiments that provide microscopic insights into the characteristics of the isotopic stabilizations of H and D aqueous species in the first solvation shell of a halide I(-) ion in water nanoclusters at low temperatures. The analysis of the simplest I(-)·(HOD) dimer shows a clear propensity for the light isotope to lie at the non-hydrogen-bonded dangling position. Our results predict that, at T ∼ 50 K, I(-)·(DOH) isomers are three times more abundant than I(-)·(HOD) ones. The reasons for such stabilization can be traced back to differences in the nuclear kinetic energy projected along directions perpendicular to the plane of the water molecule. Dynamical implications of these imbalances are shown to be reflected in the characteristics of the corresponding bands of the infrared spectroscopic signals. A similar analysis performed in larger aggregates containing ∼20 water molecules reveals, in contrast, a stabilization of the light isotope along I(-)···HO hydrogen bonds. Effects derived from the consideration of smaller halide anions with larger electric fields at the surface are also examined.

Lithium solvation in dimethyl sulfoxide-acetonitrile mixtures.

We present molecular dynamics simulation results pertaining to the solvation of Li(+) in dimethyl sulfoxide-acetonitrile binary mixtures. The results are potentially relevant in the design of Li-air batteries that rely on aprotic mixtures as solvent media. To analyze effects derived from differences in ionic size and charge sign, the solvation of Li(+) is compared to the ones observed for infinitely diluted K(+) and Cl(-) species, in similar solutions. At all compositions, the cations are preferentially solvated by dimethyl sulfoxide. Contrasting, the first solvation shell of Cl(-) shows a gradual modification in its composition, which varies linearly with the global concentrations of the two solvents in the mixtures. Moreover, the energetics of the solvation, described in terms of the corresponding solute-solvent coupling, presents a clear non-ideal concentration dependence. Similar nonlinear trends were found for the stabilization of different ionic species in solution, compared to the ones exhibited by their electrically neutral counterparts. These tendencies account for the characteristics of the free energy associated to the stabilization of Li(+)Cl(-), contact-ion-pairs in these solutions. Ionic transport is also analyzed. Dynamical results show concentration trends similar to those recently obtained from direct experimental measurements.

Surface Isotope Segregation as a Probe of Temperature in Water Nanoclusters.

Using ring polymer molecular dynamics simulations, we examine equilibrium and dynamical characteristics of solid-like, aqueous clusters that combine isotopic mixtures of HDO dilute in H2O, at temperatures intermediate between 50 and 175 K. In particular, we focus attention on the relative thermodynamic stabilities of the two isotopes at dangling hydrogen bond sites. The water octamer is analyzed as a reference system. For this aggregate, decreasing temperature yields a gradual stabilization of the light isotope at dangling sites in molecules acting as single-donor-double-acceptors of hydrogen bonds. At T ∼ 50 K, the imbalance between the corresponding quantum kinetic energies leads to a free energy difference between dangling and hydrogen bonded sites of the order of ∼2kBT. Similar free energy differences were found at dangling sites in Nw = 50 water clusters. The extent of the H/D segregation can be adequately monitored by modifications in the peak intensity of the high frequency shoulder of the stretching band of the infrared spectrum. These signals, in turn, represent a potential experimental signature of the elusive temperature of clusters in molecular beams.

Excess protons in mesoscopic water-acetone nanoclusters.

We carried out molecular dynamics simulation experiments to examine equilibrium and dynamical characteristics of the solvation of excess protons in mesoscopic, [m:n] binary polar clusters comprising m = 50 water molecules and n = 6, 25, and 100 acetone molecules. Contrasting from what is found in conventional macroscopic phases, the characteristics of the proton solvation are dictated, to a large extent, by the nature of the concentration fluctuations prevailing within the clusters. At low acetone contents, the overall cluster morphology corresponds to a segregated aqueous nucleus coated by an external aprotic phase. Under these circumstances, the proton remains localized at the surface of the water core, in a region locally deprived from acetone molecules. At higher acetone concentrations, we found clear evidence of the onset of the mixing process. The cluster structures present aqueous domains with irregular shape, fully embedded within the acetone phase. Still, the proton remains coordinated to the aqueous phase, with its closest solvation shell composed exclusively by three water molecules. As the relative concentration of acetone increases, the time scales characterizing proton transfer events between neighboring water molecules show considerable retardations, stretching into the nanosecond time domain already for n ~ 25. In water-rich aggregates, and similarly to what is found in the bulk, proton transfers are controlled by acetone/water exchange processes taking place at the second solvation shell of the proton. As a distinctive feature of the transfer mechanism, translocation pathways also include diffusive motions of the proton from the surface down into inner regions of the underlying water domain.

Excess protons in water-acetone mixtures.

Using molecular dynamics experiments, we analyze equilibrium and dynamical characteristics related to the solvation of excess protons in water-acetone mixtures. Our approach is based on the implementation of an extended valence-bond Hamiltonian, which incorporates translocation of the excess charge between neighboring water molecules. Different mixtures have been analyzed, starting from the pure water case down to solutions with a water molar fraction x(w) = 0.25. In all cases, we have verified that the structure of the first solvation shell of the H(3)O(+) moiety remains practically unchanged, compared to the one observed in pure water. This shell is composed by three water molecules acting as hydrogen bond acceptors, with no evidence of hydrogen bond donor-like connectivity. Moreover, the increment in the acetone concentration leads to a gradual stabilization of Eigen-like [H(3)O[middle dot](H(2)O)(3)](+) configurations, in detriment of Zundel-like [H[middle dot](H(2)O)(2)](+) ones. Rates of proton transfer and proton diffusion coefficients have been recorded at various water-acetone relative concentrations. In both cases, we have found a transition region, in the vicinity of x(w) ∼ 0.8, where the concentration dependences of the two magnitudes change at a quantitative level. A crude estimate shows that, at this tagged concentration, the volumes "occupied" by the two solvents become comparable. The origins of this transition separating water-rich from acetone-rich realms is rationalized in terms of modifications operated in the nearby, second solvation shell, which in the latter solutions, normally includes at least, one acetone molecule. Our results would suggest that one possible mechanism controlling the proton transfer in acetone-rich solutions is the exchange of one of these tagged acetone molecules, by nearby water ones. This exchange would give rise to Zundel-like structures, exhibiting a symmetric, first solvation shell composed exclusively by water molecules, and would facilitate the transfer between neighboring water molecules along the resonant complex.

Liquid methanol confined within functionalized silica nanopores. 2. Solvation dynamics of coumarin 153.

Equilibrium and dynamical characteristics pertaining to the solvation of the fluorescent probe coumarin 153 in liquid methanol confined within cylindrical silica pores are investigated using molecular dynamics techniques. Three kinds of pores are examined: (i) Soft hydrophobic cavities, in which wall-solvent interactions were exclusively of the Lennard-Jones type; (ii) Hydrophilic cavities, in which unsaturated oxygen sites at the wall were transformed into hydroxyl groups; (iii) Rugged pores, in which 60% of the polar groups were transformed into bulkier and mobile trimethylsilyl moieties. Equilibrium solvation structures in the three pores differ considerably: In hydrophobic environments, the solute remains adsorbed to the pore wall, with its molecular plane mostly parallel to the interface. Upon hydroxylation, the solid interface becomes preferentially coated by methanol, leading to a bistable solvation state of the probe, with alternation of "wall-like" and "bulk-like" events. An increment in the interface roughness promotes a solvation structure characterized by the embedding of the probe within a wall domain surrounded by trimethylsilyl groups. In hydrophobic environments, the relevant dynamical modes of the probe can be cast in terms of in-the-wall rotations, whereas in hydrophilic pores, out-of-the-wall evolutions are also present. The embedding of the probe at wall domains in more rugged pores, leads to restrained angular motions, with maximum amplitudes of the order of 20°. Results of early stages of the solvation response of the environment following a vertical excitation of the probe are also presented. During the initial 30 ps, we found no evidence of modifications in the spatial localizations of the probe. The overall responses are found to be between 2 and 4.5 times slower than the one observed in the bulk, being the fastest relaxation the one associated to rugged pores whereas the slowest one corresponds to hydrophilic cavities. These features are rationalized in terms of the composition of the first solvation shells and the local dynamical inhomogeneities prevailing within the different regions of the pores.

Aqueous electrolytes confined within functionalized silica nanopores.

Molecular dynamics simulations have been carried out to investigate structural and dynamical characteristics of NaCl aqueous solutions confined within silica nanopores in contact with a "bulk-like" reservoir. Two types of pores, with diameters intermediate between 20 Å and 37.5 Å, were investigated: The first one corresponded to hydrophobic cavities, in which the prevailing wall-solution interactions were of the Lennard-Jones type. In addition, we also examined the behavior of solutions trapped within hydrophilic cavities, in which a set of unsaturated O-sites at the wall were transformed in polar silanol Si-OH groups. In all cases, the overall concentrations of the trapped electrolytes exhibited important reductions that, in the case of the narrowest pores, attained 50% of the bulk value. Local concentrations within the pores also showed important fluctuations. In hydrophobic cavities, the close vicinity of the pore wall was coated exclusively by the solvent, whereas in hydrophilic pores, selective adsorption of Na(+) ions was also observed. Mass and charge transport were also investigated. Individual diffusion coefficients did not present large modifications from what is perceived in the bulk; contrasting, the electrical conductivity exhibited important reductions. The qualitative differences are rationalized in terms of simple geometrical considerations.

Structure and dynamics of liquid methanol confined within functionalized silica nanopores.

Molecular dynamics simulations have been carried out to investigate the structure and dynamics of liquid methanol confined in 3.3 nm diameter cylindrical silica pores. Three cavities differing in the characteristics of the functional groups at their walls have been examined: (i) smooth hydrophobic pores in which dispersive forces prevail, (ii) hydrophilic cavities with surfaces covered by polar silanol groups, and (iii) a much more rugged pore in which 60% of the previous interfacial hydroxyl groups were replaced by the bulkier trimethylsilyl ones. Confinement promotes a considerable structure at the vicinity of the pore walls which is enhanced in the case of hydroxylated surfaces. Moreover, in the presence of the trimethylsilyl groups, the propagation of this interface-induced spatial ordering extends down to the central region of the pore. Concerning the dynamical modes, we observed an overall slowdown in both the translational and rotational motions. An analysis of these mobilities from a local perspective shows that the largest retardations operate at the vicinity of the interfaces. The gross features of the rotational dynamics were analyzed in terms of contributions arising from bulk and surface states. Compared to the bulk dynamical behavior, the characteristic timescales associated with the rotational motions show the most dramatic increments. A dynamical analysis of hydrogen bond formation and breaking processes is also included.

Confined polar mixtures within cylindrical nanocavities.

Using molecular dynamics experiments, we have extended our previous analysis of equimolar mixtures of water and acetonitrile confined between silica walls [J. Phys. Chem. B 2009, 113, 12744] to examine similar solutions trapped within carbon nanotubes and cylindrical silica pores. Two different carbon tube sizes were investigated, (8,8) tubes, with radius R(cnt) = 0.55 nm, and (16,16) ones, with R(cnt) = 1.1 nm. In the narrowest tubes, we found that the cylindrical cavity is filled exclusively by acetonitrile; as the radius of the tube reaches approximately 1 nm, water begins to get incorporated within the inner cavities. In (16,16) tubes, the analysis of global and local concentration fluctuations shows a net increment of the global acetonitrile concentration; in addition, the aprotic solvent is also the prevailing species at the vicinity of the tube walls. Mixtures confined within silica nanopores of radius approximately 1.5 nm were also investigated. Three pores, differing in the effective wall/solvent interactions, were analyzed, (i) a first class, in which dispersive forces prevail (hydrophobic cavities), (ii) a second type, where oxygen sites at the pore walls are transformed into polar silanol groups (hydrophilic cavities), and (iii) finally, an intermediate scenario, in which 60% of the OH groups are replaced by mobile trimethylsilyl groups. Within the different pores, we found clear distinctions between the solvent layers that lie in close contact with the silica substrate and those with more central locations. Dynamical modes of the confined liquid phases were investigated in terms of diffusive and rotational time correlation functions. Compared to bulk results, the characteristic time scales describing different solvent motions exhibit significant increments. In carbon nanotubes, the most prominent modifications operate in the narrower tubes, where translations and rotations become severely hindered. In silica nanopores, the manifestations of the overall retardations are more dramatic for solvent species lying at the vicinity of trimethylsilyl groups.

Coaxial cross-diffusion through carbon nanotubes.

We present results from nonequilibrium molecular dynamics experiments describing the relaxation of local concentrations at two reservoirs, initially filled with water (W) and acetonitrile (ACN), as they become connected through a membrane composed of (16,16) carbon nanotubes. Within the hydrophobic nanotube cavities, the equilibrium concentrations contrast sharply to those observed at the reservoirs, with a clear enhancement of ACN, in detriment of W. From the dynamical side, the relaxation involves three well-differentiated stages; the first one corresponds to the equilibration of individual concentrations within the nanotubes. An intermediate interval with Fickian characteristics follows, during which the overall transport can be cast in terms of coaxial opposite fluxes, with a central water domain segregated from an external ACN shell, in close contact with the tube walls. We also found evidence of a third, much slower, mechanism to reach equilibration, which involves structural modifications of tightly bound solvation shells, in close contact with the nanotube rims.

Polar mixtures under nanoconfinement.

We present results from molecular dynamics simulations describing structural and dynamical characteristics of equimolar mixtures of water and acetonitrile, confined between two silica walls separated at interplate distances of d=0.6, 1, and 1.5 nm. Two different environments were investigated: a first one where wall-solvent dispersion forces prevail (hydrophobic confinement) and a second one in which the terminal O atoms at the silica surface are transformed into silanol groups (hydrophilic confinement). For the former case, we found that, at the shortest interplate distance examined, the confined region is devoid of water molecules. At an interplate distance of the order of 1 nm, water moves into the confined region, although, in all cases, there is a clear enhancement of the local concentration of acetonitrile in detriment of that of water. Within hydrophilic environments, we found clear distinctions between a layer of bound water lying in close contact with the silica substrates and a minority of confined water that occupies the inner liquid slab. The bound aqueous layer is fully coordinated to the silanol groups and exhibits minimal hydrogen bonding with the second solvation layer, which exclusively includes acetonitrile molecules. Dynamical characteristics of the solvent mixture are analyzed in terms of diffusive and rotational motions in both environments. Compared to bulk mixtures, we found significant retardations in all dynamical modes, with those ascribed to water molecules bound to the hydrophilic plates being the most dramatic.

Building up nanotubes: docking of "Janus" cyclodextrins in solution.

Using molecular dynamics experiments, we analyze the association of "Janus" 6-amino-6-deoxy-2O-carboxymethyl-beta-cyclodextrins (JCD) in aqueous solutions. In JCD dimers, the free energy associated with the primary-rim-secondary-rim docking shows a stable minimum of approximately -45 kcal mol(-1). Trimers in solution are also remarkably stable, exhibiting minimal distortions in their spatial and orientational distributions. The resulting geometrical docking shows the incipient characteristics of flexible nanotubes in solution, with eventual water interchange between the central channel and the bulk at the junctions between monomers. Structural and dynamical properties of the trapped water filling the nanotube are dictated to a large extent by the charge density at the rims.