You can try a multiwfn software. This is really cool tool for processing data after QM-calculation. Multiwfn works with .wfn files (if my memory is correct), so you have to add keyword "AIMPAC=.true." in $CONTRL section to generate this file. 

Mr. Gromov,

Please pay attention to the links below. There are shown examples of Gamess-US test computations. (Example 34)

http://myweb.liu.edu/~nmatsuna/gamess/tests.html

Please read carefully AIP Advances, 4, 127104 (2014), that is available free of charges. A Google Search for the article or of "Understanding Density Functional Theory and Completing It in Practice" will lead to the article. URL of the Article: http://dx.doi.org/10.1063/1.4903408

To get correct band gaps, low laying excited state energies and related properties (i.e., density of excited states), one must first ascertain that electronic structure calculations search for and attain the absolute minima of the occupied energies.  This search requires beginning the calculations with a relatively small basis set (no smaller than the minimum basis set), and carrying out successive, self-consistent calculations. Except for the first calculation, every other calculation is done with the basis set of the one preceding it plus one orbital (representing an excited states in the atomic or ionic species present in the system under study). While two consecutive calculations leading to the same occupied energies indicate a minimum (i.e., minima of occupied energies), only three consecutive calculations leading to the same occupied energies indicate the attainment of the absolute minima of the occupied energies. The minimum with just two calculation is a local one and not the absolute one (i.e., the ground state).

Now, of the three calculations giving the same occupied energies, only the first one (with the smallest basis set), provides the DFT description of the material under study.  The third of these calculations and other calculations done with basis sets obtained by augmenting that of the third calculation lead to some unoccupied energies lower than their counterpart resulting from the calculation that gives the DFT description.

Once the above steps are followed, then the density of excited states, generally up to 10 eV, are found to agree with experiment. Oh, the band gap also agrees with experiment for the calculation giving the DFT description of the material.  The physical content of the low energy excited states is destroyed when basis sets much larger than that of the first one to reach the absolute minima of the occupied energies.

When only one basis set is employed to carry out self-consistent calculations of electronic structures, the selected basis set is often found to be over-complete for the description of the ground state [resulting in a mathematical artifact (stemming from the Rayleigh theorem for eigenvalues) consisting of lowering some unoccupied energies below their values from the first of the above three calculations that lead to the same occupied energies].

It seems, I formulated question a bit ambiguously. I need to plot electron (and if possible spin) density in excited state (calculated in TDDFT formalism). 

I'm getting a couple of states of interest from Spin Flip TDDFT calculation, so they are formally excited with respect to reference state. Generally, can get these densities manually, but that will be quite long job.

Dear Oleg,

Using $td tdprp=true iroot=1 $end

Gamess writes both MO for ground state and natural orbitals of that special excited state (iroot=n) in log file,

For visualuize density of excited state you need wxmacmolplt and choosing 3d plot of density and natural orbitals of excited state

Bests,

Yavar, 

Thank you, seems promising, will try it.

______________________________________________________________

Maksym,

It wouldn't work as aimpac option writes only occupied reference orbitals

______________________________________________________________

Bojidarka, Diola

Thanks for reply, unfortunately the topic was reformulated by editor uncorrectly, here some misunderstanding arose. 

Yavar

tdprp option did not allow to plot densities (at least in my hands, may be I do not understand something), I have not found any natural orbitals. Though it prints charges, and for now it is enough, so many thanks!

Best regards

Oleg