<div dir="ltr"><div><div><div><div><div><div><div><div>Dear All<br><br></div>I am trying to generate PPs with a core hole.<br></div>My approach so far was to use input files from pslib script (us_ps_high.job) and just change the electronic configuration.<br></div>I have got some reasonable results in comparison with the experimental values for XPS shift for molecules with C 1sstar and O 1sstar PPs generated in this way. However whenever I run a calculation with the Ca 2pstar PPs it fails with<br><br>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     from cdiaghg : error #        27<br>     diagonalization (ZHEGV*) failed<br>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br></div>error message. Even when I calculate the energy of Ca 2pstar atom in a box I get the same error. I have increased Ecut off up to 150 Ry but no change.<br></div><div>Can anybody tell me why this naive way of generating PPs for C 1sstar and O 1sstar works but not C 2pstar? <br></div><div><br></div>Input file I use is the following:<br><br>-------<br> &input<br>   title='Ca',<br>   zed=20.,<br>   rel=1,<br>   config='1s2 2s2 2p5 3s2 3p6 4s2 4p0 3d-1',<br>   iswitch=3,<br>   dft='pbe'<br> /<br> &inputp<br>   lpaw=.false.,<br>   pseudotype=3,<br>   file_pseudopw='Ca.ion.pbe-spn-rrkjus.UPF',<br>   author='ADC',<br>   lloc=-1,<br>   rcloc=1.5,<br>   which_augfun='PSQ',<br>   rmatch_augfun_nc=.true.,<br>   nlcc=.true.,<br>   new_core_ps=.true.,<br>   rcore=1.0,<br>   tm=.true.<br> /<br>6<br>3S  1  0  2.00  0.00  1.20  1.30  0.0<br>4S  2  0  2.00  0.00  1.20  1.30  0.0<br>3P  2  1  6.00  0.00  1.40  1.60  0.0<br>4P  3  1  0.00  0.00  1.40  1.60  0.0<br>3D  3  2  0.00  0.05  1.20  1.70  0.0<br>3D  3  2  0.00  1.00  1.20  1.70  0.0<br>-------<br><br></div>Regards<br></div>Evren Ataman<br></div>University of Copenhagen / Department of Chemistry<br></div>