Repulsive solvent-induced interaction between C60 fullerenes in water

The Effects of External Interfaces on Hydrophobic Interactions I: Smooth Surface

External interfaces, such as the air-water and solid-liquid interfaces, are ubiquitous in nature. Hydrophobic interactions are considered the fundamental driving force in many physical and chemical processes occurring in aqueous solutions. It is important to understand the effects of external interfaces on hydrophobic interactions. According to the structural studies on liquid water and the air-water interface, the external interface primarily affects the structure of the topmost water layer (interfacial water). Therefore, an external interface may affect hydrophobic interactions. The effects of interfaces on hydrophobicity are related not only to surface molecular polarity but also to the geometric characteristics of the external interface, such as shape and surface roughness. This study is devoted to understanding the effects of a smooth interface on hydrophobicity. Due to hydrophobic interactions, the solutes tend to accumulate at external interfaces to maximize the hydrogen bonding of water. Additionally, these can be demonstrated by the calculated potential mean forces (PMFs) using molecular dynamic (MD) simulations.

Nuclear quantum effects in fullerene-fullerene aggregation in water

We studied the effects of the quantum delocalization in space of the hydrogen atoms of water in the aggregation process of two fullerene molecules. We considered a case using a purely repulsive water-fullerene interaction, as such a situation has shown that water-mediated effects play a key role in the aggregation process. This study becomes feasible, at a reduced computational price, by combining the path integral (PI) molecular dynamics (MD) method with a recently developed open-system MD technique. Specifically, only the mandatory solvation shell of the two fullerene molecules was considered at full quantum resolution, while the rest of the system was represented as a mean-field macroscopic reservoir of particles and energy. Our results showed that the quantum nature of the hydrogen atoms leads to a sizable difference in the curve of the free energy of aggregation; that is, that nuclear quantum effects play a relevant role.

Ethanol-Induced Aggregation of Nonpolar Nanoparticles in Water/Ethanol Mixed Solvents

The dispersity of nonpolar nanoparticles (NPs) in water/ethanol mixed solvents was studied using molecular dynamics simulations. Based on the rule of "like dissolves like," nonpolar NPs should be dispersed better in a solvent with a lower polarity. As the mole fraction of ethanol in a mixed solvent (R) increases from 0% (pure water) to 100% (pure ethanol), the polarity of the mixed solvent is indicated to decrease monotonically. However, the dispersity of nonpolar NP does not increase monotonically: it first decreases after the addition of a small fraction of ethanol (R < 8.0%) and then markedly increases as R further grows. When there is a small amount of ethanol, the ethanol molecules around aggregated NPs tend to simultaneously make contact with multiple NPs, which can increase the tendency of NP aggregation. Furthermore, with a considerable ethanol ratio, the interaction of the solvent with NPs becomes notably strong, which facilitates the dissolution of NPs. Our findings may help to better understand the mechanism of dispersion of NPs in mixed solvents and may provide a useful precipitation technology for NP production.

Water-mediated interactions destabilize proteins

Proteins are folded to avoid exposure of the nonpolar groups to water because water-mediated interactions between nonpolar groups are a promising factor in the thermodynamic stabilities of proteins-which is a well-accepted view as one of the unique effects of hydrophobic interactions. This article poses a critical question for this classical view by conducting an accurate solvation free-energy calculation for a thermodynamic cycle of a protein folding using a liquid-state density functional theory. Here, the solvation-free energy for a leucine zipper formation was examined in the coiled-coil protein GCN4-p1, a typical model for hydrophobic interactions, which demonstrated that water-mediated interactions were unfavorable for the association of nonpolar groups in the native state, while the dispersion forces between them were, instead, responsible for the association. Furthermore, the present analysis well predicted the isolated helical state stabilized by pressure, which was previously observed in an experiment. We reviewed the problems in the classical concept and semiempirical presumption that the energetic cost of the hydration of nonpolar groups is a driving force of folding.

Promoting Noncovalent Intermolecular Interactions Using a C60 Core Particle in Aqueous PC60s-Covered Colloids for Ultraefficient Photoinduced Particle Activity

Noncovalent intermolecular interactions in nanomaterials, such as van der Waals effects, allow adjustment of the nanoscopic size of compounds and their conformation in molecular crystal regimes. These strong interactions permit small particle sizes to be maintained as the crystals grow. In particular, these effects can be leveraged in the confined/reinforcing phase of molecules. With this in mind, we used C₆₀ molecules as a core particle in single-PC60 surfactant-covered colloid in a water-processable system. Compared with our previous results based on a PC₆₁BM core-PC60 shell particle, the PC60-C₆₀ colloid had a considerably smaller spherical structure due to the increased intermolecular interactions between C₆₀ (fullerene) molecules. Interestingly, the conformation of C₆₀ aggregates was altered depending on the mixed solvents and their volume fraction in the organic phase, which strongly affected the structural properties of the PC60-C₆₀ colloids. The particle facilitated strong interactions with a p-type core sphere when it was introduced as the shell part of a p-n heterojunction particle. This direct interaction provided effective electronic communication between p- and n-type particles, resulting in ultraefficient photonic properties, particularly in charge separation in aqueous heterostructured colloids. This enabled the development of an extremely efficient photovoltaic device with a 6.74% efficiency, which could provide the basis for creating high-performance water-processable solar cells based on p-n heterostructured NPs.

Theoretical analysis on thermodynamic stability of chignolin

Understanding the dominant factor in thermodynamic stability of proteins remains an open challenge. Kauzmann's hydrophobic interaction hypothesis, which considers hydrophobic interactions between nonpolar groups as the dominant factor, has been widely accepted for about sixty years and attracted many scientists. The hypothesis, however, has not been verified or disproved because it is difficult, both theoretically and experimentally, to quantify the solvent effects on the free energy change in protein folding. Here, we developed a computational method for extracting the dominant factor behind thermodynamic stability of proteins and applied it to a small, designed protein, chignolin. The resulting free energy profile quantitatively agreed with the molecular dynamics simulations. Decomposition of the free energy profile indicated that intramolecular interactions predominantly stabilized collapsed conformations, whereas solvent-induced interactions, including hydrophobic ones, destabilized them. These results obtained for chignolin were consistent with the site-directed mutagenesis and calorimetry experiments for globular proteins with hydrophobic interior cores.

Solvent-mediated interactions between nanostructures: From water to Lennard-Jones liquid

Solvent-mediated interactions emerge from complex mechanisms that depend on the solute structure, its wetting properties, and the nature of the liquid. While numerous studies have focused on the first two influences, here, we compare the results from water and Lennard-Jones liquid in order to reveal to what extent solvent-mediated interactions are universal with respect to the nature of the liquid. Besides the influence of the liquid, the results were obtained with classical density functional theory and brute-force molecular dynamics simulations which allow us to contrast these two numerical techniques.

Thermodynamics of hydration of fullerols [C60(OH)n] and hydrogen bond dynamics in their hydration shells

Molecular dynamics simulations of fullerene and fullerols [C₆₀(OH)_n, where n = 2-30] in aqueous solutions have been performed for the purpose of obtaining a detailed understanding of the structural and dynamic properties of these nanoparticles in water. The structures, dynamics and hydration free energies of the solute molecules in water have been analysed. Radial distribution functions, spatial density distribution functions and hydrogen bond analyses are employed to characterize the solvation shells of water around the central solute molecules. We have found that water molecules form two solvation shells around the central solute molecule. Hydrogen bonding in the bulk solvent is unaffected by increasing n. The large decrease in solvation enthalpies of these solute molecules for n > 14 enhances solubilisation. The diffusion constants of solute molecules decrease with increasing n. The solvation free energy of C₆₀ in water is positive (52.8 kJ/mol), whereas its value for C₆₀(OH)₃₀ is highly negative (-427.1 kJ/mol). The effects of surface hydroxylation become more dominant once the fullerols become soluble.

Effect of Water Hydrogen Bonding on the Solvent-Mediated "Oscillatory" Repulsion of C60 Fullerenes in Water

The solvent-mediated interaction of C60 fullerenes in liquid water is examined by using the combination of the probabilistic hydrogen bond model with the density functional theory. This combination allows one to take into account the effect of hydrogen bonding between water molecules on their interaction with fullerenes and to construct an approximation for the distribution of water molecules in the system, which provides an efficient foundation for studying hydrophobic phenomena. Our numerical evaluations predict the solvent-induced interaction of two C60 fullerenes in water at 293 K to have an oscillatory-repulsive character (previously observed in molecular dynamics simulations) only when the vicinal water-water hydrogen bonds are slightly weaker than bulk ones. Besides indicating the direction of the energetic alteration of water-water hydrogen bonds near C60 fullerenes, our model also suggests that the hydrogen bonding ability of water plays a defining role in the solvent-mediated C60-C60 repulsion.

Interaction of C₆₀ with water: first-principles modeling and environmental implications

The nature of fullerene-water interactions has been the subject of much research and debate. Specifically, the presence of a stabilizing, negative surface potential on colloidal aggregates of C60 in water is unexpected, given the neutral nature of pure carbon, and is not well understood. Previous simulation efforts have focused on the C60-water interaction using molecular dynamics simulations that lacked the ability to account for charge transfer and distribution interactions. In this study, first-principles density functional theory was used to analyze the fundamental electronic interactions to elucidate the polarization and charge transfer between water and C60. Simulations show that charge is inductively transferred to the C60 from water molecules, with subsequent polarization of the C60 molecule. In a case with two neighboring C60 molecules, the charge polarization induces a charge onto the second C60. Simulation suggests that this charge transfer and polarization may contribute at least partly to the observed negative surface potential of fullerene aggregates and, combined with hydrogen bonding network formation around C60, provides a fundamental driving force for aggregate formation in water.

Dynamic solvation shell and solubility of C60 in organic solvents

The notion of (static) solvation shells has recently proved fruitful in revealing key molecular factors that dictate the solubility and aggregation properties of fullerene species in polar or ionic solvent media. Using molecular dynamics schemes with carefully evaluated force fields, we have scrutinized both the static and the dynamic features of the solvation shells of single C60 particle for three nonpolar organic solvents (i.e., chloroform, toluene, and chlorobenzene) and a range of system temperatures (i.e., T = 250-330 K). The central findings have been that, while the static structures of the solvation shell remain, in general, insensitive to the effects of changing solvent type or system temperature, the dynamic behavior of solvent molecules within the shell exhibits prominent dependence on both factors. Detailed analyses led us to propose the notion of dynamically stable solvation shell, effectiveness of which can be characterized by a new physical parameter defined as the ratio of two fundamental time constants representing, respectively, the solvent relaxation (or residence) time within the first solvation shell and the characteristic time required for the fullerene particle to diffuse a distance comparable to the shell thickness. We show that, for the five (two from the literature) different solvent media and the range of system temperatures examined herein, this parameter bears a value around unity and, in particular, correlates intimately with known trends of solubility for C60 solutions. We also provide evidence revealing that, in addition to fullerene-solvent interactions, solvent-solvent interactions play an important role, too, in shaping the dynamic solvation shell, as implied by recent experimental trends.

Fullerene C60: surface energy and interfacial interactions in aqueous systems

The underlying mechanisms of fullerene-fullerene, fullerene-water, and fullerene-soil surface interactions in aqueous systems are not well understood. To advance our understanding of these interfacial interactions, the surface properties of Buckminsterfullerene (C60) and quartz surfaces were investigated. From application of the van Oss-Chaudhury-Good model and the Young-Dupre equation, the Lifshitz-van der Waals, acid-base, and the total surface energies of C60 powder and quartz surfaces were calculated from contact angle measurements using the sessile drop technique. C60 powder measurements indicate low to medium energy surfaces of 41.7 mJ/m2 with a dominant Lifshitz-van der Waals component. In aqueous systems, hydrophobic attraction due to the high cohesion of water is the driving force for C60 aggregation. The high free energy of hydration (DeltaG(pw)(total) = -90.5 mJ/m2) indicates the high affinity of C60 particles for water. Hamaker constants of 4.02 x 10(-21) J (A(pwp)) and 2.59 x 10(-21) J (A(pws)) were derived for C60-C60 and C60-quartz interactions in aqueous systems. The results of this study indicate that surface energy is an important physical parameter that should be considered as a basic characterization property of fullerene nanomaterials.

Potential of mean force of association of large hydrophobic particles: toward the nanoscale limit

The potentials of mean force (PMFs) were determined, in both water with the TIP3P water model and in vacuo, for systems involving formation of nonpolar dimers composed of bicyclooctane, adamantane (both an all-atom model and a sphere with the radius of 3.4 A representing adamantane), and fullerene, respectively. A series of umbrella-sampling molecular dynamics simulations with the AMBER force field were carried out for each pair under both environmental conditions. The PMFs were calculated by using the weighted histogram analysis method. The results were compared with our previously determined PMF for neopentane. The shape of the PMFs for dimers of all four nonpolar molecules is characteristic of hydrophobic interactions with contact and solvent-separated minima and desolvation maxima. The positions of all these minima and maxima change with the size of the nonpolar molecule; for larger molecules they shift toward larger distances. Comparison of the PMFs of the bicyclooctane, adamantane, and fullerene dimers in water and in vacuo shows that hydrophobic interactions in each dimer are different from that for the dimer of neopentane. Interactions in the bicyclooctane, adamantane, and fullerene dimers are stronger in vacuo than in water. These dimers cannot be treated as classical, spherical, hydrophobic objects. The solvent contribution to the PMF was also computed by subtracting the PMF determined in vacuo from that in explicit solvent. The solvent contribution to the PMFs of bicyclooctane, adamantane, and fullerene is positive, as opposed to that of neopentane. The water molecules in the first solvation sphere of both adamantane and neopentane dimers are more ordered as compared to bulk water, with their dipole moments pointing away from the surface of the dimers. The average number of hydrogen bonds per water molecule in the first hydration shell of adamantane is smaller compared to that in bulk water, but this shell is thicker for all-atom adamantane than for neopentane or a spherical model of adamantane. In the second hydration shell, the average number of hydrogen bonds is greater compared to that in bulk water only for neopentane and a spherical model of adamantane but not for the all-atom model. The strength of the hydrophobic interactions shows a linear dependence on the number of carbon atoms both in water and in vacuo. Smaller nonpolar particles interact more strongly in water than in vacuo. For larger molecules, such as bicyclooctane, adamantane and fullerene, the reversed tendency is observed.

Computer simulation study of fullerene translocation through lipid membranes

Recent toxicology studies suggest that nanosized aggregates of fullerene molecules can enter cells and alter their functions, and also cross the blood-brain barrier. However, the mechanisms by which fullerenes penetrate and disrupt cell membranes are still poorly understood. Here we use computer simulations to explore the translocation of fullerene clusters through a model lipid membrane and the effect of high fullerene concentrations on membrane properties. The fullerene molecules rapidly aggregate in water but disaggregate after entering the membrane interior. The permeation of a solid-like fullerene aggregate into the lipid bilayer is thermodynamically favoured and occurs on the microsecond timescale. High concentrations of fullerene induce changes in the structural and elastic properties of the lipid bilayer, but these are not large enough to mechanically damage the membrane. Our results suggest that mechanical damage is an unlikely mechanism for membrane disruption and fullerene toxicity.

Supramolecular self-organization in PEO-modified C60 fullerene/water solutions: influence of polymer molecular weight and nanoparticle concentration

Utilizing a first-principles-based coarse-grained implicit solvent model, we have investigated the self-association of C(60) fullerenes that have been symmetrically modified with six grafted poly(ethylene oxide) (PEO) chains in aqueous solution. Despite the highly symmetric nature of the pair interactions between PEO-grafted fullerenes, their supramolecular assemblies are highly anisotropic and resemble the linear clusters formed in Stockmayer fluids. The dipole-like interaction between these symmetrically modified fullerenes results from the shielding of the C(60) fullerenes by PEO, favoring the addition of more PEO-grafted fullerenes to the linear clusters at the relatively unprotected ends. At low nanoparticle concentrations, self-association is dominated by the formation of stable dimers and trimers resulting from fullerene-fullerene contact and favorable PEO-fullerene interactions. With increasing nanoparticle concentration, larger clusters become increasingly probable. The molecular weight of the PEO tethers can be treated as a temperature-like analogue, with a reduction in average cluster size with increasing chain length due to increased steric repulsion, which is qualitatively similar to effects observed in Stockmayer fluids with increasing temperature. The role of PEO in supramolecular self-organization in PEO-modified C(60) fullerene/water solutions is complex, contributing not only to steric stabilization but also to favorable energetic interactions, nanoparticle shielding, and depletion-driven aggregation.

Structure and UV−Vis Spectrum of C60 Fullerene in Ethanol: A Sequential Molecular Dynamics/Quantum Mechanics Study

A molecular dynamics simulation combined with semiempirical quantum mechanics calculations has been performed to investigate the structure, dynamical, and electronic properties of pure C60 in liquid ethanol. The behavior of the fullerene alcoholic solution was obtained by using the NPT ensemble under ambient conditions, including one C60 fullerene immersed in 1000 ethanol molecules. Our analyzed center-of-mass pairwise radial distribution function indicated that, on average, there are 32, 72, 132, and 187 ethanol molecules around, respectively, the first, second, third, and fourth solvation shells of the C60 molecule. To investigate the UV-vis transition energies of C60 in the presence of ethanol, we have considered constituents of the time uncorrelated supramolecular structures of the first solvation shell, i.e., clusters of C60@{EtOH}32 types. The semiempirical calculations were performed at the intermediate neglect of differential overlap level with configuration interaction singles (INDO/CIS). Our results have pointed out that the characteristic C60 UV-vis absorbance peaks are slightly shifted to longer wavelengths, as compared to the isolated molecule. These findings are in connection with the weak donor-acceptor character of the interactions involving electron lone pairs of oxygen atoms on the solvent and the fullerene surface.

Application of the level-set method to the implicit solvation of nonpolar molecules

A level-set method is developed for numerically capturing the equilibrium solute-solvent interface that is defined by the recently proposed variational implicit solvent model [Dzubiella, Swanson, and McCammon, Phys. Rev. Lett. 104, 527 (2006); J. Chem. Phys. 124, 084905 (2006)]. In the level-set method, a possible solute-solvent interface is represented by the zero level set (i.e., the zero level surface) of a level-set function and is eventually evolved into the equilibrium solute-solvent interface. The evolution law is determined by minimization of a solvation free energy functional that couples both the interfacial energy and the van der Waals type solute-solvent interaction energy. The surface evolution is thus an energy minimizing process, and the equilibrium solute-solvent interface is an output of this process. The method is implemented and applied to the solvation of nonpolar molecules such as two xenon atoms, two parallel paraffin plates, helical alkane chains, and a single fullerence C(60). The level-set solutions show good agreement for the solvation energies when compared to available molecular dynamics simulations. In particular, the method captures solvent dewetting (nanobubble formation) and quantitatively describes the interaction in the strongly hydrophobic plate system.

Effects of lengthscales and attractions on the collapse of hydrophobic polymers in water

We present results from extensive molecular dynamics simulations of collapse transitions of hydrophobic polymers in explicit water focused on understanding effects of lengthscale of the hydrophobic surface and of attractive interactions on folding. Hydrophobic polymers display parabolic, protein-like, temperature-dependent free energy of unfolding. Folded states of small attractive polymers are marginally stable at 300 K and can be unfolded by heating or cooling. Increasing the lengthscale or decreasing the polymer-water attractions stabilizes folded states significantly, the former dominated by the hydration contribution. That hydration contribution can be described by the surface tension model, DeltaG = gamma(T)DeltaA, where the surface tension, gamma, is lengthscale-dependent and decreases monotonically with temperature. The resulting variation of the hydration entropy with polymer lengthscale is consistent with theoretical predictions of Huang and Chandler [Huang DM, Chandler D (2000) Proc Natl Acad Sci USA 97:] that explain the blurring of entropy convergence observed in protein folding thermodynamics. Analysis of water structure shows that the polymer-water hydrophobic interface is soft and weakly dewetted, and is characterized by enhanced interfacial density fluctuations. Formation of this interface, which induces polymer folding, is strongly opposed by enthalpy and favored by entropy, similar to the vapor-liquid interface.

A molecular dynamics simulation study of buckyballs in water: Atomistic versus coarse-grained models of C60

A simple coarse-grained water-C(60) model potential is proposed to study the solvation behavior of C(60) in water. Extensive molecular dynamics simulations are performed to compare the hydration behavior of the coarse-grained C(60) with that of its atomistic counterpart. In the coarse-grained description, we model the C(60) as a spherical solute and its interaction with water is represented by a simple two-body central potential as obtained from a coarse graining of the interactions of a water molecule with all the atoms of the C(60). On the other hand, the atomistic model takes into account a discrete nature of all the atoms of C(60) explicitly. Molecular dynamics simulations are carried out in an isothermal-isobaric ensemble at normal temperature and pressure. Various spatial and orientational correlations of the water around the two model C(60)s are compared. The coarse-grained model is shown to reproduce the hydration behavior of the C(60) quite well. The simplicity and the computational economy of the coarse-grained model will allow for simulations of self-assembly processes of a much larger system over a longer period of time.

Molecular dynamics simulation study of the structure of poly(ethylene oxide) brushes on nonpolar surfaces in aqueous solution

The structure of poly(ethylene oxide) (PEO, M(w) = 526) brushes of various grafting density (sigma) on nonpolar graphite and hydrophobic (oily) surfaces in aqueous solution has been studied using atomistic molecular dynamics simulations. Additionally, the influence of PEO-surface interactions on the brush structure was investigated by systematically reducing the strength of the (dispersion) attraction between PEO and the surfaces. PEO chains were found to adsorb strongly to the graphite surface due primarily to the relative strength of dispersion interactions between PEO and the atomically dense graphite compared to those between water and graphite. For the oily surface, PEO-surface and water-surface dispersion interactions are much weaker, greatly reducing the energetic driving force for PEO adsorption. This reduction is mediated to some extent by a hydrophobic driving force for PEO adsorption on the oily surface. Reduction in the strength of PEO-surface attraction results in reduced adsorption of PEO for both surfaces, with the effect being much greater for the graphite surface where the strong PEO-surface dispersion interactions dominate. At high grafting density (sigma approximately 1/R(g)(2)), the PEO density profiles exhibited classical brush behavior and were largely independent of the strength of the PEO-surface interaction. With decreasing grafting density (sigma < 1/R(g)(2)), coverage of the surface by PEO requires an increasingly large fraction of PEO segments resulting in a strong dependence of the PEO density profile on the nature of the PEO-surface interaction.

Formalism for calculation of polymer-solvent-mediated potential

A simple theoretical approach is proposed for calculation of a solvent-mediated potential (SMP) between two colloid particles immersed in a polymer solvent bath in which the polymer is modeled as a chain with intramolecular degrees of freedom. The present recipe is only concerned with the estimation of the density profile of a polymer site around a single solute colloid particle instead of two solute colloid particles separated by a varying distance as done in existing calculational methods for polymer-SMP. Therefore the present recipe is far simpler for numerical implementation than the existing methods. The resultant predictions for the polymer-SMP and polymer solvent-mediated mean force (polymer-SMMF) are in very good agreement with available simulation data. With the present recipe, change tendencies of the contact value and second virial coefficiency of the SMP as a function of size ratio between the colloid particle and polymer site, the number of sites per chain, and the polymer concentration are investigated in detail. The metastable critical polymer concentration as a function of size ratio and the number of sites per chain is also reported for the first time. To yield the numerical solution of the present recipe at less than on a personal computer, a rapid and accurate algorithm for the numerical solution of the classical density functional theory is proposed to supply rapid and accurate estimation of the density profile of the polymer site as an input into the present formalism.

Water-Induced Interactions between Carbon Nanoparticles

Molecular dynamics simulations were carried out in order to study the hydration of C60 fullerenes, carbon nanotubes, and graphene sheets in aqueous solution and the nature of water-induced interactions between these carbon nanoparticles. The hydration of these nonpolar carbon nanoparticles does not exhibit classical hydrophobic character due to the high density of surface atoms (carbon) resulting in strong water-surface dispersion interactions. Water was found to wet the nanoparticle surfaces independent of nanoparticle surface curvature, with the decrease in the extent of water-water hydrogen bonding with decreasing surface curvature being offset by stronger water-surface interactions. While all carbon nanoparticles investigated are anticipated to aggregate in water due to strong direct nanoparticle-nanoparticle interactions, the water-induced interactions between nanoparticles were found to be repulsive and, in contrast to the wetting behavior, were observed to exhibit strong dependence on surface curvature. The strength of the water-induced interaction between carbon nanoparticles was found to correlate well with the number of hydration water molecules displaced upon particle aggregation, which, relative to the amount of direct nanoparticle-nanoparticle contact engendered upon aggregation, decreases with decreasing surface curvature.

A molecular-dynamics simulation study of solvent-induced repulsion between C60 fullerenes in water

Molecular-dynamics simulations of a single C(60) fullerene and pairs of C(60) fullerenes in aqueous solution have been performed for the purpose of obtaining improved understanding of the nature of solvent-induced interactions between C(60) fullerenes in water. Our simulations reveal repulsive solvent-induced interactions between two C(60) fullerenes in aqueous solution in contrast to the associative effects observed for conventional nonpolar solutes. A decomposition of the solvent-induced potential of mean force between fullerenes into entropy and energy (enthalpy) contributions reveals that the water-induced repulsion between fullerenes is energetic in origin, contrasting strongly to entropy-driven association observed for conventional nonpolar solutes. The dominance of energy in the solvent-induced interactions between C(60) fullerenes arises primarily from the high atomic density of the C(60) molecule, resulting in strong C(60)-water van der Waals attraction that is reduced upon association of the fullerenes. The water-induced repulsion is found to decrease with increasing temperature due largely to an increasing contribution from a relatively weak entropy-driven association.

Molecular dynamics simulation study of the role of evenly spaced poly(ethylene oxide) tethers on the aggregation of C60 fullerenes in water

The aggregation behavior of C60 fullerenes and C60 fullerenes with six symmetrically tethered poly(ethylene oxide) oligomers [(PEO)-6-C60] in aqueous solutions has been studied using implicit solvent molecular dynamics simulations. Our simulations reveal that while the attraction between two (PEO)-6-C60 fullerenes in aqueous solution is stronger and longer range than that between two bare C60 fullerenes, the (PEO)-6-C60 fullerenes do not phase-separate in water but rather aggregate in chain-like clusters at concentrations where unmodified fullerenes completely phase-separate.