How to interpret molden AO/MO coefficients?

Question

A simple calculation with Molpro 2012 on H2 molecule in cc-pvtz basis produces a molden file with 28 molecular orbitals. Each MO is defined with 30 AO/MO coefficient, but the number of symmetric AOs is 32, and the number of contractions is 28.
memory, 200, m

basis = cc-pvtz;
symmetry nosym;
geomtyp = xyz;
geometry = {
2

H       0.00000000     0.00000000    -0.46904000
H       0.00000000     0.00000000     0.46904000
}

{rhf;
 WF, 2, 1, 0;
}

put, molden, h2.molden;
---

Part of Molpro output
 NUMBER OF PRIMITIVE AOS:          34
 NUMBER OF SYMMETRY AOS:           32
 NUMBER OF CONTRACTIONS:           28   (  28A   )

Why the number of AO/MO coefficients is not equal to the number of symmetric AOs or the number of contractions?

Susi Lehtola · Answer

Okay, so there are many layers to this question. cc-pVTZ for H is [5s2p1d/3s], which comes out to 3 + 2*3 + 5 = 3+6+5 = 14 basis functions per atom, which are composed of 16 primitives (the contracted s function).
Now, while there are 1 and 3 cartesians for the s and p shells, for the d shell you have 6 cartesian functions but only 5 spherical functions. This is important, since Gaussian-basis integrals are evaluated in the cartesian basis, and then contracted to get the spherical-basis integral. (There are closed-form relations for this translation.)
For H2, you then have

2*14 = 28 contracted basis functions; this is the "number of contractions".
2*15 = 30 contracted cartesian basis functions, which is thus probably the format which Molden is using.
2*16 = 32 primitive Gaussian basis functions with spherical D functions, this appears to be the "symmetry AOs".
2*18 = 34 primitive Gaussian cartesian basis functions; this is the "primitive AOs".

Hope this helps.

How to interpret molden AO/MO coefficients?

One Answer

Add your own answers!

Ask a Question