That would be great, really. I did that in this way (for what I called DZP-DZP):
%block PAO.Basis Fe 2 n=4 0 2 P 1 0.000 0.000 n=3 2 2 P 1 0.000 0.000 %endblock PAO.Basis Then I varied EnergyShift and SplitNorm iteratively within a shell script to allow Siesta find the combination of those two parameters that results in the minimum energy for the single-atom Fe unit cell. Finally, I saved the generated PAO.Basis to be used in other calculations. I do not know neither if that is the correct way to generate this DZP-DZP basis set nor if this is the best protocol to minimize basis sets. Roberto ________________________________ From: Oleksandr Voznyy <[email protected]> To: [email protected] Sent: Wednesday, January 7, 2009 7:38:02 PM Subject: Re: [SIESTA-L] Still energy differences and BSSE corrections Dear Victor, I knew that one can make more than one polarized orbital but didn't think that one might ever need to generate three polarized orbitals per each angular momentum. How do you generate them? Could you show your basis block to make things more clear? Sincerely, Alexander
