The file .hesswk looks as

  30 2 2 2 2 2 2 2 2 2 2 3 3 3 3 4 4 4 5 5 6 6 6 6 6 6 6 6 6 6-1
  20 1 1 1 1 1 2 2 3 3 4 5 5 5 5 5 5 5 5 5 6
  20 1 1 1 1 1 2 2 3 3 4 5 5 5 5 5 5 5 5 5 6

    0    0    0         0.000000000        -0.368369197

         0.000000000         0.000000000        -0.000000292
        -0.000000158         0.000000000         0.000000146
         0.000000158         0.000000000         0.00000014
         0.000000000         0.000000000         0.000000000
    1    1    1         0.010000000        -0.368338591
        0.006116027         0.000000000         0.000084854
       -0.003027212         0.000000000         0.002424326

       -0.003088815         0.000000000        -0.002509179

       -0.004992946         0.000000000         0.765405092

etc etc

the last line represents dipole moment, but when I run a crystal, DMol3 will put zeros on this line, how can I to compute by hand using Mulliken?