Molecular modeling of freezing of simple fluids confined within carbon nanotubes

We report Monte Carlo simulation results for freezing of Lennard-Jones carbon tetrachloride confined within model multiwalled carbon nanotubes of different diameters. The structure and thermodynamic stability of the confined phases, as well as the transition temperatures, were determined from parallel tempering grand canonical Monte Carlo simulations and free-energy calculations. The simulations show that the adsorbate forms concentric molecular layers that solidify into defective quasi-two-dimensional hexagonal crystals. Freezing in such concentric layers occurs via intermediate phases that show remnants of hexatic behavior, similar to the freezing mechanism observed for slit pores in previous works. The adsorbate molecules in the inner regions of the pore also exhibit changes in their properties upon reduction of temperature. The structural changes in the different regions of adsorbate occur at temperatures above or below the bulk freezing point, depending on pore diameter and distance of the adsorbate molecules from the pore wall. The simulations show evidence of a rich phase behavior in confinement; a number of phases, some of them inhomogeneous, were observed for the pore sizes considered. The multiple transition temperatures obtained from the simulations were found to be in good agreement with recent dielectric relaxation spectroscopy experiments for CCl(4) confined within multiwalled carbon nanotubes.     

Melting behavior of water in cylindrical pores: carbon nanotubes and silica glasses

We report a study of the effects of confinement in multi-walled carbon nanotubes and mesoporous silica glasses (SBA-15) on the solid structure and melting of both H(2)O and D(2)O ice, using differential scanning calorimetry, dielectric relaxation spectroscopy, and neutron diffraction. Multi-walled nanotubes of 2.4, 3.9 and 10 nm are studied, and the SBA-15 studied has pores of mean diameter 3.9 nm; temperatures ranging from approximately 110 to 290 K were studied. We find that the melting point is depressed relative to the bulk water for all systems studied, with the depression being greater in the case of the silica mesopores. These results are shown to be consistent with molecular simulation studies of freezing in silica and carbon materials. The neutron diffraction data show that the cubic phase of ice is stabilized by the confinement in carbon nanotubes, as well as in silica mesopores, and persists up to temperatures of about 240 K, above which there is a transition to the hexagonal ice structure.     

Effects of nonionic surfactant C12E5 on the cooperative dynamics of water

A dielectric relaxation study of binary mixtures of nonionic surfactant C12E5 + water has been made as a function of temperature in the isotropic micellar, lamellar, and hexagonal regions of the phase diagram. Two dielectric dispersion steps were found and could be assigned to the intermolecular cooperative dynamics of water at the micellar interface and in the bulk water domains. A quantitative analysis is given. The relaxation amplitudes were used to determine effective hydration numbers. The activation energies of water relaxation were calculated from the relaxation times. The data indicate weaker surfactant-water and water-water interactions near the micellar interface compared to those of bulk liquid water. Further analysis revealed the presence of water clusters large enough to show a cooperative relaxation mode even at high surfactant concentrations. However, the relaxation time of this mode is larger compared to that of pure water. This points out the importance of confinement effects on water dynamics.     

An experimental study of melting of CCl4 in carbon nanotubes

We report dielectric relaxation spectroscopy measurements of the melting point of carbon tetrachloride confined within open-tip multi-walled carbon nanotubes with two different pore diameters, 4.0 and 2.8 nm. In both cases, a single transition temperature well above the bulk melting point was obtained for confined CCl4. These results contrast with what was obtained in our previous measurements using carbon nanotubes with a pore diameter of 5.0 nm, where multiple transition temperatures both above and below the bulk melting point of CCl4 were observed. Our experimental measurements are consistent with our recent molecular simulation results (F. R. Hung, B. Coasne, E. E. Santiso, K. E. Gubbins, F. R. Siperstein and M. Sliwinska-Bartkowiak, J. Chem. Phys., 2005, 122, 144706). Although the simulations overestimate the temperatures in which melting upon confinement occurs, both simulations and experiments suggest that all regions of adsorbate freeze at the same temperature, and that freezing occurs at higher temperatures upon reduction of the pore diameter.     

Thermodynamics of water entry in hydrophobic channels of carbon nanotubes

Experiments and computer simulations demonstrate that water spontaneously fills the hydrophobic cavity of a carbon nanotube. To gain a quantitative thermodynamic understanding of this phenomenon, we use the recently developed two phase thermodynamics method to compute translational and rotational entropies of confined water molecules inside single-walled carbon nanotubes and show that the increase in energy of a water molecule inside the nanotube is compensated by the gain in its rotational entropy. The confined water is in equilibrium with the bulk water and the Helmholtz free energy per water molecule of confined water is the same as that in the bulk within the accuracy of the simulation results. A comparison of translational and rotational spectra of water molecules confined in carbon nanotubes with that of bulk water shows significant shifts in the positions of the spectral peaks that are directly related to the tube radius.     

Relaxation dynamics of a polymer in a 2D confinement

The molecular dynamics of oligomeric poly(propylene glycol) (PPG) liquids (MW=1000, 2000, and 4000 g/mol) confined in a two-dimensional layer-structured Na-vermiculite clay has been studied by broadband dielectric spectroscopy. The alpha-relaxation and the normal mode relaxation processes were studied for all samples in bulk and confinement. The most prominent experimental observation was that for the normal mode process: the relaxation rate in the clay is drastically shifted to lower frequencies compared to that of the bulk material. This slowing down is probably caused by the strongly reduced number of accessible chain conformations in two dimensions. Also the temperature dependence of the relaxation time for the normal mode process is strongly affected by the confinement. In contrast, for the alpha-relaxation of the confined polymers we observed only a slight increase of the relaxation rate at high temperatures compared to the corresponding bulk samples, and a decrease of its relaxation strength relative to the beta relaxation. Thus, the glass transition is unaffected by the 2D confinement, suggesting that the underlying phenomena responsible for the glass transition is the same as in bulk. Moreover, in the clay the intensity of the normal mode is stronger than that of the alpha-process, in contrast to the bulk samples where the opposite behavior is observed.     

Fast Water Relaxation through One‐Dimensional Channels by Rapid Energy Transfer

Water in carbon nanotubes is surrounded by hydrophobic carbon surfaces and shows anomalous structural and fast transport properties. However, the dynamics of water in hydrophobic nanospaces is only phenomenologically understood. In this study, water dynamics in hydrophobic carbon nanotubes is evaluated based on water relaxation using nuclear magnetic resonance spectroscopy and molecular dynamics simulations. Extremely fast relaxation (0.001 s) of water confined in carbon nanotubes of 1 nm in diameter on average is observed; the relaxation times of water confined in carbon nanotubes with an average diameter of 2 nm (0.40 s) is similar to that of bulk water (0.44 s). The extremely fast relaxation time of water confined in carbon nanotubes with an average diameter of 1 nm is a result of frequent energy transfer between water and carbon surfaces. Water relaxation in carbon nanotubes of average diameter 2 nm is slow because of the limited number of collisions between water molecules. The dynamics of interfacial water can therefore be controlled by varying the size of the hydrophobic nanospace.     

Dielectric relaxation in bulk and cylindrically confined octylcyanobiphenyl (8CB)

The dielectric relaxation spectrum over the frequency range 10² to 1.8×10⁹ Hz of 4-octyl-4'-cyanobiphenyl, 8CB, in bulk and confined to 200 nm diameter cylindrical pores is reported. We used matrices with parallel cylindrical pores, obtaining different alignments of the molecular director depending on the treatment. Results show that there are two relaxations in the isotropic phase and in the mesophases for parallel alignment and three for perpendicular alignment. The molecular origin of theses modes and the effect of the confinement on their dynamics are discussed. To compare properly the results for bulk and confined 8CB, a re-scaling of the experimental data is proposed.     

Optical absorption of semiconducting and metallic nanospheres with the confined electron-phonon coupling

We study the optical absorption, especially the (far-) infrared absorption by phonons, of semiconducting and metallic nanospheres. In the nanoscopic sphere, phonons as well as states of electronic excitations are quantized by confinement. It is also known that in the nanoscopic geometry, the confined electron-phonon interaction has a different form from the usual one in the bulk. First, we analyze the phonon and electron contributions to the dielectric response of nanospheres like epsilon(q,omega)=epsilon(ph)(q,omega)+epsilon(el)(q,omega) or 1epsilon(q,omega)=1epsilon(sc-ph)(q,omega)+1epsilon(el)(q,omega) from the confined electron-phonon interaction for three cases: the intrinsic semiconductor, the doped semiconductor, and the metal. From the dielectric response, the optical absorption spectra are calculated within the semiclassical framework concentrating on the (far-) infrared region and compared to the spectra without imposing confinement. Nontrivial differences of the spectra with confined phonons stem from two features: the electron-phonon coupling matrix has a different form and the phase space q of the confined phonon is reduced because of its quantization to q(n). Finally, size distribution effects in an ensemble of isolated nanospheres are briefly discussed. Those effects are found to be important in metallic spheres with rapid sweepings of resonances by a small change of the sphere size.     

Dielectric relaxation in bulk and cylindrically confined octylcyanobiphenyl (8CB)

The dielectric relaxation spectrum over the frequency range 10² to 1.8×10⁹ Hz of 4‐octyl‐4′‐cyanobiphenyl, 8CB, in bulk and confined to 200 nm diameter cylindrical pores is reported. We used matrices with parallel cylindrical pores, obtaining different alignments of the molecular director depending on the treatment. Results show that there are two relaxations in the isotropic phase and in the mesophases for parallel alignment and three for perpendicular alignment. The molecular origin of theses modes and the effect of the confinement on their dynamics are discussed. To compare properly the results for bulk and confined 8CB, a re‐scaling of the experimental data is proposed.     

Influence of cylindrical submicrometer confinement on the static and dynamic properties in nonyloxycyanobiphenyl (9OCB)

Broadband dielectric spectroscopy (10(2)-1.9 x 10(9) Hz) and specific heat measurements have been performed on nonyloxycyanobiphenyl (9OCB) in the isotropic (I), nematic (N), and smectic A (SmA) phases confined to 200 nm diameter parallel cylindrical pores of Anopore membranes. Untreated and HTBA-treated membranes have been found to obtain axial and radial confinements, respectively. However, structural or configurational transitions in untreated membranes have been reported to exist in the SmA-mesophase of 9OCB. Both confinements clearly affect the N-I and SmA-N phase transitions. In the axial confinement, the analysis of the specific heat and static dielectric permittivity data leads to a second order SmA-N phase transition, which is known to be weakly first order for bulk 9OCB. Dynamic dielectric measurements have accounted for the different molecular motions in both confinements. On both mesophases, either N or SmA, the relaxation processes in axial configuration are faster than in the bulk. However, in radial confinement, they are either equal or slower than in the bulk. Additionally, there are no differences in the energy barrier hindering the molecular motions between the axial and radial confinements and even in relation to bulk. Likewise, dielectric results suggest that the extension inside the pores of the surface pinned molecular layer (proved to be temperature-dependent) persists at high enough temperature as a residual-thin layer adjacent to the pore wall.     

Dynamics of small-molecule glass formers confined in nanopores

We report a comparative neutron scattering study of the molecular mobility and nonexponential relaxation of three structurally similar glass-forming liquids, isopropanol, propylene glycol, and glycerol, both in bulk and confined in porous Vycor glass. Confinement reduces molecular mobility in all three liquids, and suppresses crystallization in isopropanol. High-resolution quasielastic neutron scattering spectra were fit to Fourier transformed Kohlrausch functions exp[-(t∕τ)(β)], describing the α-relaxation processes in these liquids. The stretching parameter β is roughly constant with wavevector Q and over the temperature range explored in bulk glycerol and propylene glycol, but varies both with Q and temperature in confinement. Average relaxation times <τ(Q)> are longer at lower temperatures and in confinement. They obey a power law <τ(Q)> ∝ Q(-γ), where the exponent γ is modified by confinement. Comparison of the bulk and confined liquids lends support to the idea that structural and∕or dynamical heterogeneity underlies the nonexponential relaxation of glass formers, as widely hypothesized in the literature.     

水中受拉伸负载下单壁纳米碳管机械性质的分子动力学研究

利用分子动力学方法研究了(5,5)扶手椅型和(10,10)锯齿型纳米碳管在水中受拉伸负载下的机械性质.通过计算纳米碳管中氧和氢原子的局部密度分布研究了限制效应.结果表明,碳管在水中的杨式系数与在真空下相同,而碳管在水中的拉伸应力小于在真空中的.     

Local dielectric spectroscopy of near-surface glassy polymer dynamics

A noncontact scanning probe microscopy method was used to probe local near-surface dielectric susceptibility and dielectric relaxation in polyvinyl acetate near the glass transition. Dielectric spectra were measured from 10(-4) to 10(2) Hz as a function of temperature. The measurements probed a 20 nm thick layer below the free surface of a bulk film. A significant change in the fragility index and moderate narrowing of the distribution of relaxation times were found in the near-surface layer. In contrast to results for ultrathin films confined on or between metallic electrodes, no reduction in the dielectric strength was found, inconsistent with the immobilization of slower modes.     

Ultrafast energy transfer in water-AOT reverse micelles

A spectroscopic investigation of the vibrational dynamics of water in a geometrically confined environment is presented. Reverse micelles of the ternary microemulsion H2O/AOT/n-octane (AOT = bis-2-ethylhexyl sulfosuccinate or aerosol-OT) with diameters ranging from 1 to 10 nm are used as a model system for nanoscopic water droplets surrounded by a soft-matter boundary. Femtosecond nonlinear infrared spectroscopy in the OH-stretching region of H2O fully confirms the core/shell model, in which the entrapped water molecules partition onto two molecular subensembles: a bulk-like water core and a hydration layer near the ionic surfactant headgroups. These two distinct water species display different relaxation kinetics, as they do not exchange vibrational energy. The observed spectrotemporal ultrafast response exhibits a local character, indicating that the spatial confinement influences approximately one molecular layer located near the water-amphiphile boundary. The core of the encapsulated water droplet is similar in its spectroscopic properties to the bulk phase of liquid water, i.e., it does not display any true confinement effects such as droplet-size-dependent vibrational lifetimes or rotational correlation times. Unlike in bulk water, no intermolecular transfer of OH-stretching quanta occurs among the interfacial water molecules or from the hydration shell to the bulk-like core, indicating that the hydrogen bond network near the H2O/AOT interface is strongly disrupted.     

Effect of confinement on the liquid-liquid phase transition of supercooled water

We report on an observation of the phase transition between two liquid phases of supercooled confined water in simulations. The temperature of the liquid-liquid transition of water at zero pressure slightly decreases due to confinement in the hydrophobic pore. The hydrophilic confinement affects this temperature in the opposite direction and shifts the critical point of the liquid-liquid transition to a higher pressure. As a result, in a strongly hydrophilic pore the liquid-liquid phase transition becomes continuous at zero pressure, indicating the shift of its critical point from negative to a positive pressure. These findings indicate that experimental studies of water confined in the pores of various hydrophobicity/hydrophilicity may clarify the location of the liquid-liquid critical point of bulk water.     

Testing the core/shell model of nanoconfined water in reverse micelles using linear and nonlinear IR spectroscopy

A core/shell model has often been used to describe water confined to the interior of reverse micelles. The validity of this model for water encapsulated in AOT/isooctane reverse micelles ranging in diameter from 1.7 to 28 nm (w0 = 2-60) and bulk water is investigated using four experimental observables: the hydroxyl stretch absorption spectra, vibrational population relaxation times, orientational relaxation rates, and spectral diffusion dynamics. The time dependent observables are measured with ultrafast infrared spectrally resolved pump-probe and vibrational echo spectroscopies. Major progressive changes appear in all observables as the system moves from bulk water to the smallest water nanopool, w0 = 2. The dynamics are readily distinguishable for reverse micelle sizes smaller than 7 nm in diameter (w0 = 20) compared to the response of bulk water. The results also demonstrate that the size dependent absorption spectra and population relaxation times can be quantitatively predicted using a core-shell model in which the properties of the core (interior of the nanopool) are taken to be those of bulk water and the properties of the shell (water associated with the headgroups) are taken to be those of w0 = 2. A weighted sum of the core and shell components reproduces the size dependent spectra and the nonexponential population relaxation dynamics. However, the same model does not reproduce the spectral diffusion and the orientational relaxation experiments. It is proposed that, when hydrogen bond structural rearrangement is involved (orientational relaxation and spectral diffusion), dynamical coupling between the shell and the core cause the water nanopool to display more homogeneous dynamics. Therefore, the absorption spectra and vibrational lifetime decays can discern different hydrogen bonding environments whereas orientational and spectral diffusion correlation functions predict that the dynamics are size dependent but not as strongly spatially dependent within a reverse micelle.     

Confinement and surface effects on the molecular dynamics of a nematic mixture investigated by dielectric relaxation spectroscopy

Broadband dielectric spectroscopy (10(-2)-10(9) Hz) was employed to investigate the molecular dynamics of the liquid crystalline mixture E7 confined in both untreated and lecithin-treated 20 nm Anopore membranes. Because E7 does not crystallize, it was possible to cover a temperature range of more than 200 K, providing an exhaustive dielectric characterization of a liquid crystal confined to Anopore membranes for the first time. In the nematic state, the tumbling (alpha-) and the delta-relaxation are observed, also under confinement conditions. The analysis of their relative intensities give that the orientation of the E7 molecules is preferentially axial in untreated but opposite radial in lecithin-treated pores. The radial alignment of the liquid crystals in the modified membrane is understood as a tail-to-tail conformation of E7 molecules imposed by the adsorbed lecithin molecules. The relaxation rate of the alpha-process is enhanced for E7 confined in native Anopore compared with the bulk and E7 in treated pores. This is interpreted as resulting from a less dense molecular packing of E7 in the middle of the pore compared to the bulk. In both untreated and treated membranes, the relaxation rate of the delta-process is higher than in the bulk, and the values of the respective Vogel-Fulcher-Tammann temperatures depend on the actual surface treatment. Additionally, a surface process, due to molecular fluctuations of molecules within an adsorbed layer at the pore wall, was detected.     

Dielectric relaxation of aqueous solutions of hydrophilic versus amphiphilic peptides

We report on molecular dynamics simulations of the frequency-dependent dielectric relaxation spectra at room temperature for aqueous solutions of a hydrophilic peptide and an amphiphilic peptide at two concentrations. We find that only the high-concentration amphiphilic peptide solution exhibits an anomalous dielectric increment over that of pure water, while the hydrophilic peptide exhibits a significant dielectric decrement. The dielectric component analysis carried out by decomposing these peptide solutions into peptide, hydration layer, and outer layer(s) of water clearly shows the presence of a unique dipolar component with a relaxation time scale on the order of approximately 25 ps (compared to the bulk water time scale of approximately 11 ps) that originates from the interaction between the hydration layer water and the outer layer(s) of water. Results obtained from the dielectric component analysis further show the emergence of a distinct and much lower frequency relaxation process for the high-concentration amphiphilic peptide compared to the hydrophilic peptide due to strong peptide dipolar couplings to all constituents, accompanied by a slowing of the structural relaxation in all water layers, giving rise to time scales close to approximately 1 ns. We suggest that the molecular origin of the dielectric relaxation anomalies is due to frustration in the water network arising from the amphiphilic chemistry of the peptide that does not allow it to reorient on the picosecond time scale of bulk water motions. This explanation is consistent with the idea of the "slaving" of residue side chain motions to protein surface water, and furthermore offers the possibility that the anomalous dynamics observed from a number of spectroscopies arises at the interface of hydrophobic and hydrophilic domains on the protein surface.     

Dynamics of propylene glycol and its oligomers confined to a single molecular layer

The dynamics of propylene glycol (PG) and its oligomers 7-PG and poly-propylene glycol (PPG), with M(w) = 4000 (approximately 70 monomers), confined in a Na-vermiculite clay have been investigated by quasielastic neutron scattering. The liquids are confined to single molecular layers between clay platelets, giving a true two-dimensional liquid. Data from three different spectrometers of different resolutions were Fourier transformed to S(Q,t) and combined to give an extended dynamical time range of 0.3-2000 ps. An attempt was made to distinguish the diffusive motion from the methyl group rotation and a fast local motion of hydrogen in the polymer backbone. The results show that the average relaxation time tau(d) of this diffusive process is, as expected, larger than the relaxation time tau averaged over all dynamical processes observed in the experimental time window. More interesting, it is evident that the severe confinement has a relatively small effect on tau(d) at T = 300 K, this holds particularly for the longest oligomer, PPG. The most significant difference is that the chain-length dependence of tau(d) is weaker for the confined liquids, although the slowing down in bulk PG due to the formation of a three-dimensional network of OH-bonded end groups reduces this difference. The estimated average relaxation time tau at Q = 0.92 Angstroms(-1) for all the observed processes is in excellent agreement with the previously reported dielectric alpha relaxation time in the studied temperature range of 260-380 K. The average relaxation time tau (as well as the dielectric alpha relaxation time) is also almost unaffected by the confinement to a single molecular layer, suggesting that the interaction with the clay surfaces is weak and that the reduced dimensionality has only a weak influence on the time scale of all the dynamical processes observed in this study.