What to do for processing the electronic/spin density?

For this calculation you can use the files of your ground state run. Therefore, follow this sequence:

1. Run SIESTA first; the total/spin densities are produced (the file <job>.RHO) provided that the SaveRho flag in the <job>.fdf file is set to TRUE
2. After that you should run get_param_siesta to get param.inc file
3. Compile lev00 using the script lev00.comp (provided)
4. Run lev00 as explained in Section 3.6.3

Thus, you will need the following SIESTA output files:

• <job>.XV file with the final geometry
• <job>.RHO file for the electron/spin density