<span class="gI">Dear CP2k users,<br><br> I am learning how to generate GTH pseudopotential using atom.x. I ran the atom.x program with the input file shown below for Ne:<br><br>File: atom.dat:<br><br> Ne<br> PADE<br> r<br>
60 10 30.0<br> 1.3d0 50.0<br> 1 3<br> 2 0 2.00 0.00<br> 2 1 6.00 0.00<br> 3 0 0.00 0.00<br><br><br><br> I got the following screen output. It seems that the program did run, and created three files, namely <a href="http://atom.ae">atom.ae</a>, psp.par and weights.par. Unfortunately, they are all EMPTY with nothing in it. <br>
<br> Has anyone encountered such error, or my input file is somewhat incorrect ?<br><br> Thanks.<br><br>Hanning Chen<br>Department of Chemistry<br>Northwestern University<br>Evanston, IL 60208<br><br>Screen output:<br><br>
Ne all electron calculation<br> ---------------------------<br><br> r e l a t i v i s t i c ! !<br><br> correlation = PADE non spin-polarized<br><br> nuclear charge = 10.000000<br> number of core orbitals = 1<br>
number of valence orbitals = 4<br> electronic charge = 10.000000<br> ionic charge = 0.000000<br><br><br> input data for orbitals<br><br> i n l s j occ<br><br> 1 1 0 0.5 0.5 2.0000<br>
2 2 0 0.5 0.5 2.0000<br> 3 2 1 -0.5 0.5 2.0000<br> 4 2 1 0.5 1.5 4.0000<br> 5 3 0 0.5 0.5 0.0000<br><br><br> radial grid parameters<br><br> r(1) = .0 , r(2) =0.153620E-06 , ... , r( 492) = 58.11173790<br>
a = 10.00171382 b = 30.00000000<br><br> iter = 1 dvmax = 0.368E+00 xmixo =0.353E+00<br> iter = 2 dvmax = 0.202E+00 xmixo =0.353E+00<br> iter = 3 dvmax = 0.972E-01 xmixo =0.353E+00<br> iter = 4 dvmax = 0.414E-01 xmixo =0.353E+00<br>
iter = 5 dvmax = 0.189E-01 xmixo =0.353E+00<br> iter = 6 dvmax = 0.104E-01 xmixo =0.353E+00<br> ..................<br> iter = 51 dvmax = 0.129E-10 xmixo =0.353E+00<br> iter = 52 dvmax = 0.849E-11 xmixo =0.353E+00<br>
spliq<br><br><br><br> <br><br></span>