plotting density of open shell

[cavallo]

Hi,

working with an open-shell C(S) symmetric system, I would like to plot the total (alpha+beta) density of the a' and a'' irriducible representations in the xz-plane separately.  From the manual it seems I've to uste TM-5.10, and put some lines into the control file.

It is unclear what lines... what I've understood so far is something like :

$pointval geo=plane
  grid1 vector 1 0 0 range -10,10 points 501
  grid2 vector 0 0 1 range -10,10 points 501
  origin 0  0  0   

however, I think this way I get the total (a'+a'') density into the td.plt file.  Should I modify the first line something like:

$pointval geo=plane mo 1a' 2a' ... up to the last occupied a' MO

to get the total a' (alpha-beta) density ? Smells wrong to me, but this is the best guess I can imagine.
TIA
Luigi

[Arnim]

Hi,

I think, there is no automatic way to do that.
If you use the option
$pointval mo
then each MO is written on seperated files, ordered by symmetry and alpha-beta.
If you want to visualize the density of all the MOs belonging to
the same symmetry, you would have to sum them up somehow.

With
$pointval geo=plane mo 1-10 fmt=xyz
the informations are written as plane 3d coordinates.
Then you have to sum up the values in 4th column, probably using a script.

I hope this helps a bit
Arnim

[cavallo]

Oh, it helps a lot.  It seems this is the way to go.
Summing up density from the various MO is easy.
Thanks,
Luigi

[mpjohans]

Hello!

Another way that should work just fine is to edit the control file and comment out first the a' and then the a" orbital occupations, and run "dscf/ridft -proper". This way, only the density of the a'/a" orbitals will be considered in the analysis part.

Have a nice day,
    Mikael J.
    http://www.iki.fi/~mpjohans