Dielectric constant: diel.f

This utility computes static dielectric permittivity (dielectric constant) from fluctuations of the total molecular momenta using expression (S.W.deLeeuw, J.W.Perram, E.R.Smith, Ann.Rev.Phys.Chem., 37, 245 (1986)):

$$\epsilon = 1 + \frac{4\pi}{3kTV} \langle \vert \sum_i \vec{\mu}_i \vert^2 \rangle$$ (13)

where $\vec{\mu}_i$ is the dipole moment of molecule $i$ and the sum is taken over all the molecules. Note that this formula is applicable only for small and rigid (or almost rigid) molecules. It is important that the molecules were not broken by the periodic boundary conditions.

Compilation:

f77 -O3 -o diel diel.f tranal_base.f

This utility needs only a single parameter - the temperature, which should be written directly after the trajectory parameters.