The structure and dynamic properties of interfacial water at the graphite and silica solid surfaces were investigated using molecular dynamics simulations. The effect of surface properties on the characteristics of interfacial water was quantified by computing density profiles, radial distribution functions, surface density distributions, orientation order parameters, and residence and reorientation correlation functions. In brief, our results show that the surface roughness, chemical heterogeneity, and surface heterogeneous charge distribution affect the structural and dynamic properties of the interfacial water molecules, as well as their rate of exchange with bulk water. Most importantly, our results indicate the formation of two distinct water layers at the SiO2 surface covered by a large density of hydroxyl groups. Further analysis of the data suggests a highly confined first layer where the water molecules assume preferential hydrogen-down orientation and a second layer whose behavior and characteristics are highly dependent on those of the first layer through a well-organized hydrogen bond network. The results suggest that water−water interactions, in particular hydrogen bonds, may be largely responsible for macroscopic interfacial properties such as adsorption and contact angle.

Although carbon nanotubes have attracted enormous research interest, their practical application is still hindered, primarily, by the difficulty of separating them into samples monodispersed in diameter, chirality, and length. Recent advances show that ultracentrifugating carbon nanotube dispersions stabilized by surfactants is a promising route for achieving the desired separation. For further perfectioning this procedure it is necessary to know how surfactants adsorb on nanotubes of different diameters, which determines the nanotube-surfactant aggregate effective density and the nanotube-nanotube potential of mean force. Because only limited experimental data are available to elucidate these phenomena, we report here an extensive all-atom molecular dynamics study on the morphology of sodium dodecyl sulfate (SDS) surfactant aggregates adsorbed on (6,6), (12,12), and (20,20) single walled carbon nanotubes at room conditions. Our calculations reveal that the nanotube diameter is the primary factor that determines the morphology of the aggregates because of a competition between the entropic and energetic advantage encountered by the surfactants when they wrap one nanotube, and the enthalpic penalty faced during this process due to bending of the surfactant molecule. The data are in qualitative agreement with the neutron scattering results reported by Yurekli et al. [J. Am. Chem. Soc. 2004, 126, 9902], and for the first time provide an atomic-level description helpful in designing better separation, as well as stabilization techniques for aqueous carbon nanotube dispersions.

The aggregate structure of sodium dodecyl sulfate (SDS) adsorbed at the graphite-water interface has been studied with the aid of molecular dynamics (MD) simulations. As expected, our results show that adsorbed SDS yields hemi-cylindrical micelles. The hemi-cylindrical aggregates in our simulations closely resemble all structural and morphological details provided by previous solution atomic force microscopy (AFM) experiments. More interestingly, our data indicate that SDS head groups do not provide a complete shield to the hydrophobic tails. Instead, we found regions in which the hydrophobic tails are exposed to the aqueous solution. By conducting a parametric study for SDS-like nonionic surfactants we show that electrostatic interactions between SDS head groups and counterions are responsible for the unexpected result. Our interpretation is corroborated by density profiles, analysis of the coordination states, and mean square displacement data for both the adsorbed SDS surfactants and the counterions in solution. Counterion condensation appears to be a physical phenomenon that could be exploited to direct the assembly of advanced nanostructured materials.

We use a combination of molecular dynamics simulations and density functional theory calculations to investigate the energetic disorder in fullerene systems. We show that the energetic disorder evaluated from an ensemble average contains contributions of both static origin (time-independent, due to loose packing) and dynamic origin (time-dependent, due to electron-vibration interactions). In order to differentiate between these two contributions, we compare the results obtained from an ensemble average approach with those derived from a time average approach. It is found that in both amorphous C60 and C70 bulk systems, the degrees of static and dynamic disorder are comparable, while in the amorphous PC61BM and PC71BM systems, static disorder is about twice as large as dynamic disorder.

We investigate the impact of electronic polarization, charge delocalization, and energetic disorder on the charge-transfer (CT) states formed at a planar C 60 /pentacene interface. The ability to examine large complexes containing up to seven pentacene molecules and three C 60 molecules allows us to take explicitly into account the electronic polarization effects. These complexes are extracted from a bilayer architecture modeled by molecular dynamics simulations and evaluated by means of electronic-structure calculations based on long range-separated functionals (ωB97XD and BNL) with optimized range-separation parameters. The energies of the lowest charge-transfer states derived for the large complexes are in very good agreement with the experimentally reported values. The average singlet-triplet energy splittings of the lowest CT states are calculated not to exceed 10 meV. The rates of geminate recombination as well as of dissociation of the triplet excitons are also evaluated. In line with experiment, our results indicate that the pentacene triplet excitons generated through singlet fission can dissociate into separated charges on a picosecond time scale, despite the fact that their energy in C 60 /pentacene heterojunctions is slightly lower than the energies of the lowest CT triplet states.3