| 10:01:07 AM PDT - Mon, Aug 26th 2013 |
|
Hi Alexx,
The design of Gaussian basis is a rather complex science/art. The overall aim is to construct a, preferentially, number of functions that span the space of shapes of orbitals accurately for a target set of properties. Aside from what the target properties are of a particular basis set and how the required accuracy is selected, there is also some freedom in how a basis set is represented. In principle one could choose basis functions such that there is a one-to-one correspondence between a basis function and a particular atomic orbital (see e.g. ANO basis sets), but one does not have to do that. Basis sets that use atomic orbitals have as a disadvantage that each basis function contains many exponents (often they are generally contracted basis sets) which leads to high integral evaluation costs unless the integral code has been designed to specifically support these kinds of basis sets. In segmented contracted basis sets one exponent contributes to only one basis function. This makes the integral evaluation cheaper but eliminates the one-to-one correspondence between basis functions and orbitals. For the physics these differences in representation do not matter as long as both representations span the same space of functions.
NWChem currently uses integral codes that have been designed for segmented contracted basis sets. Therefore it is beneficial to use basis sets where the basis functions are not atomic orbitals. The basis set you are using is a segmented contracted basis set and hence you got the result you showed. If you were to use the ANO-RCC basis set for Hydrogen you should get one coefficient that is much closer to one.
Best wishes, Huub
|