The origins of the hydrophobic effect are widely thought to lie in structural changes of the water molecules surrounding a nonpolar solute. The spatial distribution functions of the water molecules surrounding benzene and cyclohexane computed previously from molecular dynamics simulations show a high density first hydration shell surrounding both solutes. In addition, benzene showed a strong preference for hydrogen bonding with two water molecules, one to each face of the benzene ring. The position data alone, however, do not describe the majority of orientational changes in the water molecules in the first hydration shells surrounding these solutes. In this paper, we measure the changes in orientation of the water molecules with respect to the solute through spatial orientation functions as well as radial͞angular distribution functions. These data show that the water molecules hydrogen bonded to benzene have a strong orientation preference, whereas those around cyclohexane show a weaker tendency. In addition, the water-water interactions within and between the first two hydration shells were measured as a function of distance and ''best'' hydrogen bonding angle. Water molecules within the first hydration shell have increased hydrogen bonding structure; water molecules interacting across shell 1 and shell 2 have reduced hydrogen bonding structure.hydrophobic effect ͉ orientation I t is widely understood that the structural changes in water induced by exposure to nonpolar surfaces are responsible for the phenomenon known as the hydrophobic effect (1). The small size, tetrahedral geometry, and hydrogen-bonding (H-bonding) ability of water all play a role in the interactions water takes part in with itself as well as with polar and nonpolar solutes. The structure of liquid water is well characterized experimentally in terms of radial distribution functions, g(r), derived from diffraction experiments (2, 3). Similar experiments on solutions of alcohols and tetraalkylammonium ions, however, show very little structural differences in g(r) compared with bulk water (2, 4). Yet the thermodynamics of insertion of nonpolar solutes into water show an unfavorable free energy, a positive entropy at room temperature, and a large positive heat capacity (5, 6). There are several possible reasons why the g(r) may lack sensitivity in detecting structural changes in the hydration shells of hydrophobic solutes: (i) solubility of such solutes is small, so the majority of the observed signal arises from bulk water; (ii) solutes may be incompletely mixed (7); and (iii) the g(r) is inherently insensitive to orientational structure.The tetrahedral geometry of the water molecule and the directional nature of the H-bond impart a highly oriented structure to water. The orientational correlation function of liquid water has been generated by applying maximum entropy methods to g(r) functions (8,9). The resulting density distribution shows clear density peaks beyond the H-bond donor locations (along the O-H bonds), with a broader distribution...