Unphysical results with custom Mg basis set

[Richard Fogarty]

Dear all,

I'm trying to use a custom basis-set for Mg in CP2K but am getting 
unphysical results; I suspect the reason is that i'm making a mistake with 
the basis set formatting. The simplest unphysical result is for a 
calculation on a single atom, whereby the ground state has p-orbitals 
(rather than s) fully occupied (see attached *.inp and *.cpout files).

The basis set is generated in an external program using the following steps:
1) A calculation is carried out on an Mg atom with the s-shell occupied 
(neutral atom, no net spin). This calculation yields a list of grid points 
(distances from the nucleus, due to the spherical symmetry), and the value 
of each shell at these grid points (which I'll call *\phi_{r}*). I'm using 
the same pseudopotential ("Mg GTH-PBE-q2") in this calculation as in the 
subsequent CP2K calculations.
2) The value of *\phi_{r}* is multiplied by sqrt((2l+1))*(r**l)
3) Gaussians are fit to the data from step 2. 

Im loading the exponents and coefficients from step 3) into CP2K. One 
example is attached (file MG_TEST_BASIS). Note that the basis set attached 
leads to reasonable results for a single atom in the external program 
(where i know im entering it in the correct format).

Any help would be greatly appreciated.

Thanks in advance,
Richard


Misc info:
This is using cp2k-6.1 in serial mode (cp2k.sopt executable)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180911/b2d6f2e0/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Mg_atom.inp
Type: chemical/x-gamess-input
Size: 1357 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180911/b2d6f2e0/attachment.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mg_singl_atom.cpout
Type: application/octet-stream
Size: 120438 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180911/b2d6f2e0/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MG_TEST_BASIS
Type: application/octet-stream
Size: 991 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180911/b2d6f2e0/attachment-0001.obj>