Electron densities

For the above mentioned programs setting of keyword

$pointval dens

or simply

$pointval

yields calculation of densities

$$\rho(\vec{R}_P)=\sum_{\nu\mu} D_{\nu\mu}\phi_{\nu}(\vec{R}_P)\phi_{\mu}(\vec{R}_P)$$ (12.1)

dens

on an orthogonal grid ($R_P$ mathend000#), the size of which is automatically adjusted to the size of the molecule and the resolution is adjusted to yield acceptable gOpenMol plots (for specification of non-default grid types (planes, lines) and non-default output formats see Section 14.2.17).

Names of output files are:

td.plt

total density (UHF: $\alpha$ mathend000# density plus $\beta$ mathend000# density )

sd.plt

spin density ($\alpha$ mathend000# density minus $\beta$ mathend000# density )

mp2d.plt

MP2 density

mp2sd.plt

MP2 spin density

ed.plt

differential density for excited state

esd.plt

differential spin density for excited state

<myname>.plt

general density passed e.g. by the RICC2 program.

The .plt files may directly be visualized by gOpenMol; the file coord.xyz, which is also necessary for gOpenMol, is generated by the above programs, if $pointval is set in the control-file.

Two-component wave functions (only module RIDFT and only if $soghf is set): Total density is on file td.plt like for one-component wave functions; this is also true for all other quantities depending only on the density matrix (electrostatic potential etc.). sd.plt contains the absolute value of the spin vector density, which is the absolute value of the following vector:

$$s_i(\mathbf{r})=\left(\begin{array}{cc}\Psi_\alpha^*&\Psi_\beta^*... ..._i\left(\begin{array}{c}\Psi_\alpha\\ \Psi_\beta\end{array}\right) ~~~~~i=x,y,z$$

$pointval fmt=txt

leads to a file containing the spin vector density vectors, which can be used by gOpenMol. It is advisable to choose ca. one Bohr as the distance between two gridpoints.