Multiwfn is able to plot RDF of density (norm of orbital wavefunction) for any kind of orbital via subfunction 5 of main function 200. However, if you plot RDF for orbital wavefunction of p, d, f... type, then is expected that the result should be zero, because in the angular integral, the positive and negative parts exactly cancel with each other.
If you want to plot RDF of density of an orbital, the correct way is: Enter subfunction 26 of main function 6, set occupation number of all orbitals to zero first, and then set occupation number of the orbital of interest to 1.0. Then go to the RDF module, select electron density as the function to be studied, you should be able to get expected result.
Best regards,
Tian
]]>Does anyone have any advise as to how to address this?
Regards
Tavar