Getting radial Kohn-Sham orbital from the output of an electronic structure code using gaussian basis sets

As we know that in solving the Schrodinger equation (Kohn-Sham equation), we expand the wavefunction in terms of some basis sets and them optimize the energy (the expectation value of Hamiltonian) by varying the coefficients of the basis in the expansion. Most of the electronic structure codes use the basis sets. What I need is that if I use gaussian basis sets in NWCHEM for example, how would I get radial Kohn-Sham orbitals or the optimized coefficients of the basis? I will be grateful if someone can tell any other option to get these orbitals from different electronic structure code but the one which uses Gaussian basis sets.

• Every qc code should have some documented option to print this information. Did you check the manual? May 21 '19 at 6:46
• I am aware of getting molecular orbitals in NWCHEM ("nwchem-sw.org/index.php/…) in Gaussian cube format but do not know how to get the radial orbitals directly or from the corresponding cube file. I know how to get the radial density. This I learned from a similar question "chemistry.stackexchange.com/questions/70021/…" but for radial density. May 21 '19 at 11:18
• I have trouble understanding the question, particularly I do not know what is meant by radial Kohn-Sham orbital. With the optimised coefficients don't you just mean one of the density matrices? May 21 '19 at 13:17
• Okay, I understand your point. The systems I am interested in have spherical symmetry (like closed-shell atoms, so the wavefunction can be separated in radial and angular part). Even if it's not spherically symmetric, it should be the wavefunction which is expanded in terms of the basis set chosen (I am interested in knowing the optimized coefficients of this expansion) and the density matrices then be calculated from the wavefunction. May 21 '19 at 14:14
• For example "physics.stackexchange.com/questions/420899/…". May 21 '19 at 14:23