Optimization of PAO by yourself

Starting from the primitive basis functions, you can optimize the radial shape variationally so that the accuracy can be increased. See the details in the Section 'Orbital optimization'.