[molpro-user] LCC2 population analysis

[Borys Szefczyk]

Dear Molpro users,

I'm trying to do an MPA based on the density from LCC2 calculations.
I am able to save the density and compute some other needed properties
with the command

lt-df-lcc2;eom,-3.1,singlet=1,tranes=-2.1,propes=-2.1,startle=1,trans=1,densave=3333.2;expec,dm

Then, I do

{pop;density,3333.2,type=CC2,state=1.1}
{pop;density,3333.2,type=FULL,state=2.1}

to get the charges in the ground- and excited state. The problem is that they
sum-up to 0, whereas my molecule has +1 net charge. Not to mention, that the 'Total'
column looks like the nucleus charge was already subtracted...
I've tried to do MPA "by hand", using MATROP program and I've got exactly the same
result. The charges look fine at the Hartree-Fock level. What am I doing wrong?

Best regards,
Borys Szefczyk

-- 
  Molecular Modelling & Quantum Chemistry Group,
  Institute of Physical & Theoretical Chemistry,
  Wroclaw University of Technology
  http://ichfit.ch.pwr.wroc.pl/people/szefczyk