"+1" is summed for epsilon for each atom group
The dielectric module sums (\varepsilon-1) + 1
for each atom group in epsilon.py (line 340 and line 375).
I.e. for a single atom group results are correct, but for multiple atomgroups the resulting dielectric constant is shifted by 1 for each further atomgroup.
This also affects 0.6!