<div>Thanks for the quick response Paolo. I missed it in the documentation, thinking incorrectly that the &test section wouldn't have anything to do with the construction of the pseudopotential. <br><br>Best,<br>Brad<br>
UC Berkeley<br><br> </div><blockquote style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;" class="gmail_quote">you need to specify a test configuration (containing all valence  <br>
states) in the &test section<br><br>P.<br>---<br>Paolo Giannozzi, Dept of Physics, University of Udine<br>via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></blockquote><pre>
</pre><br><br><br><blockquote style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;" class="gmail_quote"><br><br></blockquote><div class="gmail_quote">On Sat, Nov 14, 2009 at 6:35 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;"><div class="im">To be clear, on my second attempt the bottom of my file looks like :<br><br>
/<br> 5<br></div><div class="im"> 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></div> 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<div class="im">
<br>
 3S  1  0   2.00 0.00   0.50  0.50  0.00<br><br><br></div>Thanks,<br><font color="#888888">Brad<br><br></font><div class="gmail_quote"><div class="im">On Sat, Nov 14, 2009 at 6:19 PM, Brad Malone <span dir="ltr"><<a href="mailto:brad.malone@gmail.com" target="_blank">brad.malone@gmail.com</a>></span> wrote:<br>

</div><div><div></div><div class="h5"><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></div></div><br>
</blockquote></div><br>