add model II of the PK1985 paper
suppress k-space output by default
compute the dielectric constant