Accurate description of water structure affects simulation of protein folding, substrate binding, macromolecular recognition, and complex formation. We study the hydration of buckminsterfullerene, the smallest hydrophobic nanosphere, by molecular dynamics simulations using a state-of-the-art quantum mechanical polarizable force field (QMPFF3), derived from quantum mechanical data at the MP2/aug-cc-pVTZ(-hp) level augmented by CCSD(T). QMPFF3 calculation of the hydrophobic effect is compared to that obtained with empirical force fields. Using a novel and highly sensitive method, we see polarization increases ordered water structure so that the imprint of the hydrophobic surface atoms on the surrounding waters is stronger and extends to long-range. We see less water order for empirical force fields. The greater order seen with QMPFF3 will affect biological processes through a stronger hydrophobic effect.B uckminsterfullerene (C 60 ), commonly known as the buckyball, has been the focus of many experimental and theoretical studies because of its promise in both biomedicine and nanomaterials (1-8). Carbon nanotubes, a nanomaterial similar to fullerenes, exhibit unique mechanical properties, electrical and heat conductivity; they are often hailed as a twenty-first century material that could revolutionize a number of industries (2). The behavior of fullerenes in water is of special interest for biosensors (9). Moreover, water-soluble buckyball derivatives could potentially be used for a wide variety of medical applications, including HIV protease inhibition, DNA photoclevage, antibacterial activity, and photocytotoxicity for cancer treatment (10).The interaction of water with hydrophobic surfaces is responsible for many biological processes, such as protein folding, substrate binding and biological self-assembly of micelles, lyotropic mesophases, and lipid membranes. The molecular surface of proteins can have extended hydrophobic patches, so that buckyballs, graphite surfaces, and carbon nanotubes provide model systems for study of water molecules near large hydrophobic surfaces. An accurate description of interactions between water and its surrounding environment is essential for understanding and simulating such processes. Molecular dynamics (MD) simulations are used to study the behavior of water molecules next to hydrophobic surfaces at atomic detail and subpicosecond time resolution. The hydration structure of nonpolar solutes of varying size has been studied extensively in the past (11-22) using empirical force fields. Work from our lab has used empirical potential energy functions with popular water models to model the effect of a single molecule of benzene, cyclohexane (18, 21), methane, and C 60 (22) on the structural details of water around these solutes. All these studies used force fields with nonpolarizable fixed point charge water models and analyzed water structure with a low-sensitivity method unable to detect the effect of solute surface "roughness" on water structure.Polarization, which allows the ef...