For the four thermodynamic states: temperature T ) 283. 15, 298.15, 313.15, and 328.15 K and the corresponding bulk water density F ) 0.9997, 0.9970, 0.9922, and 0.9875 g cm -3 , for which experimental data are available, we have studied hydration structure, hydration thermodynamics, and molecular conformational equilibria for n-butane in water at infinite dilution, by means of the hypernetted chain closure reference interaction site model (HNC-RISM) theory with an all-atom solute model. The hydration structures of the trans and the gauche conformers of n-butane are presented and analyzed at the atomic level in terms of the atomic solute-solvent radial distribution functions. With these radial distribution functions as input, the n-butane conformational average hydration free energies, energies, enthalpies, and entropies are calculated. At room temperature, the normalized equilibrium distribution of n-butane conformers, the water solvent-induced rotational free energy surface and the trans-gauche and trans-cis cavity thermodynamic properties are calculated. With the optimized nonbonded potential parameters based on the CHARMM96 all-atom model for alkanes (Yin, D.; Mackerell, A. D., Jr. J. Comput. Chem. 1998, 19, 334), n-butane hydration thermodynamics and its conformational equilibria in water are well described by the HNC-RISM theory in comparison with the available experimental and computer simulation results. We also calculated the solute density derivatives of the water-water radial distribution functions δh vv , with the optimized CHARMM96 all-atom model, the united-atom OPLS (optimized potentials for liquid simulations), and the all-atom OPLS models for n-butane, respectively. The δh vv (r) reflect the effect of increased pressure disrupting the hydrogen bonding between water molecules. The all-atom model seems to enhance such an effect due to the well-documented shortcoming of the RISM theory in the treatment of the excluded volume of so-called auxiliary sites.
I. IntroductionTwo important problems in the study of the aqueous solution behavior of biomolecules are the restructuring of water around the molecules and the solvent influence on the conformational structure of the molecules. Small flexible hydrocarbon molecules have been used as models for studying those problems. 1 Among them, n-butane, which has only one conformational degree of freedom, is perhaps the simplest molecule available for studying the two problems at the same time. n-butane has been studied extensively in the gas phase, 1-8 in the neat liquid, 9-17 and in apolar 15,18-23 and aqueous solutions. 15,16,23-33 There are two kinds of models for representing the n-butane molecule. One is a united-atom model, namely that the methyl and methylene groups are approximately considered as one site of pair interactions, respectively. In this model, the hydrogens of the CH n groups are implicitly treated, the interaction sites for the CH n groups are centered on the carbons, and the CH n groups are all taken as neutral. The other one is an all-...