To be clear, on my second attempt the bottom of my file looks like :<br><br>/<br> 5<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> 4S  2  0   2.00 0.00   2.0   2.0   0.00<br> 4P  3  1   1.00 0.00   2.0   2.0   0.00<br>
 3S  1  0   2.00 0.00   0.50  0.50  0.00<br><br><br>Thanks,<br>Brad<br><br><div class="gmail_quote">On Sat, Nov 14, 2009 at 6:19 PM, Brad Malone <span dir="ltr"><<a href="mailto:brad.malone@gmail.com">brad.malone@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">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>
</blockquote></div><br>