Perhaps you should consider consulting the literature about various basis set kinds and their peculiarities: Gaussian functions, all sorts of plane waves and wavelets, localized muffin-tin orbitals, numerical atomic orbitals, finite grids and splines... -- just to understand why you don't want Gaussian bases in SIESTA or why you can't do an all-electron calculation with it.
The answers to your questions are: (1) SIESTA does support any basis functions that are a product of a radial part and a spherical harmonic AND have a finite support, and hence, (2) it IS possible to set them manually in input (almost exactly), but it's like the last thing anybody would want to do with bases in SIESTA. If you need both linear scaling and Gaussian basis sets, your choice is pretty much limited to Gaussian 03 (which is, to my opinion, not quite good) and CP2k (they use a combined Gaussian--plane wave basis set for KS orbitals and density, respectively). 2007/2/5, Oleksandr Loboda <[EMAIL PROTECTED]>:
If I understood you correctly SIESTA does not support any of Dunning-Hay, Huzinaga basis sets? And it is impossible to set them manually in input.