Hi, I'm trying to generate a pseudopotential for Gallium with the entire third shell taken as valence (i.e., I want 3s2 3p6 3d10 4s2 4s1 all in valence). I didn't see anything in the /atomic_doc/pseudo_library/ that had an example where there were multiple occupied states of the same angular momentum.  My first attempt at doing so was with the following input file:<br>
<br>  &input<br>   atom = 'Ga'<br>   iswitch = 3<br>   rel = 1<br>   dft = 'PZ'<br>/<br> 8<br> 1S  1  0   2.00  1<br> 2S  2  0   2.00  1<br> 2P  2  1   6.00  1<br> 3S  3  0   2.00  1<br> 3P  3  1   6.00  1<br>
 3D  3  2  10.00  1<br> 4S  4  0   2.00  1<br> 4P  4  1   1.00  1<br>&inputp<br>   pseudotype = 2<br>   file_pseudopw = 'Ga.UPF'<br>   lloc = 0<br>   nlcc = .false.<br>   tm = .true.<br>/<br> 3<br> 3P  2  1   6.00  0.00  0.50  0.50  0.00<br>
 3D  3  2   10.00 0.00  0.50  0.50  0.00 <br> 3S  1  0   2.00 0.00   0.50  0.50  0.00<br><br>After generating the pseudopotential with this file, it became clear that I had the wrong valence (no 4s or 4p electrons presumably). So I tried to edit the file and add these in to the &inputp field. For the 4S and 4P I used a cutoff of 2 au, but the pseudopotential generation failed and said that I had the wrong number of nodes. Am I doing something wrong? I appreciate any help or suggestions. Thanks!<br>
<br>Best,<br>Brad <br>UC Berkeley<br><br>