We present the theory of the analytical first and second derivatives of the Hartree-Fock energy for a molecular solute with respect to the nuclear coordinates, within the framework of the polarizable continuum model. The formulation refers to a cavity with an accurately modelled molecular shape.