  We construct a Hamiltonian for a single domain protein where the contact
enthalpy and the chain entropy decrease linearly with the number of native
contacts. The hydration effect upon protein unfolding is included by modeling
water as ideal dipoles that are ordered around the unfolded surfaces, where the
influence of these surfaces, covered with an ``ice-like'' shell of water, is
represented by an effective field that directs the water dipoles. An
intermolecular pair interaction between water molecules is also introduced. The
heat capacity of the model exhibits the common feature of small globular
proteins, two peaks corresponding to cold and warm unfolding, respectively. By
introducing vibrational modes, we obtain quantitatively good accordance with
experiments.
