where H 0 ({r, ṙ} i ) is the Hamiltonian for the i th molecule, comprising the bonded potential energy of a molecule U 0 ({r} i ), and the kinetic energy ∑ N i j=1 m j v 2 j /2, where, v j is the velocity of the j th particle in the i th molecule. The functional of interaction energy W [{ φ }] is the usual Flory-Huggins mixing term plus the incompressibility term that is usually employed in HhPF. In terms of a local energy density functional w( φ (r)), this functional can be written as: