Cube format and ECPs

[antti_karttunen]

Hello,

I noticed a small issue when plotting molecular orbitals using the new cube format (a very useful addition, by the way!). If a system contains atoms that use an ECP, dscf and ridft do not take the ECP electrons into account when writing the cube header. For example, gold (79 electrons, 60-electron ECP) is written as "19" in the header. Due to this, visualization programs interpret the gold atoms as potassium atoms (19 electrons). This is not a big issue, as it can be easily remedied with a simple sed command wrapped in a shell script. Maybe it can be corrected in a future release...

[christof.haettig]

Can send this: turbomole@cosmologic.de ?

Do you have such a script which can provide in this forum for other users?

Regards,
Christof

[antti_karttunen]

Hello,

I'll also report the issue to turbomole@cosmologic.de. The script I was talking about is available at
http://cc.joensuu.fi/~ankarttu/tm/tmcube.tar.gz
I have not tested it extensively, but it worked for all my cube files. It also includes another feature that is useful for those who visualize the cubes with GaussView.