Difference between revisions of "CamCASP/Programming/5"

From CUC3
Jump to navigation Jump to search
import>Am592
 
import>Am592
Line 18: Line 18:
 
# More consistent behaviour of aug-cc-pV''n''Z series.
 
# More consistent behaviour of aug-cc-pV''n''Z series.
 
# Don't have rotation matrices coded, but these have been derived.
 
# Don't have rotation matrices coded, but these have been derived.
  +
  +
===Simple example: 1===
  +
Here's an example that demonstrates the inadequacy of the spherical GTOs. Consider a ''1p1s'' (yes, in this order) main basis and a single occupied orbital. This would be the ''px'' orbital. The space the auxiliary basis needs to span is ''1p1s x 1p1s'' = ''1s1p1d''. So let us construct this auxiliary basis and perform the DF. Here's the question: '''Can the density be exactly fitted using this auxiliary basis?'''
  +
  +
The answer depends on which basis type is used for the auxiliary basis.
  +
  +
* Cartesian: The auxiliary basis consists of the following combinations:
  +
<math>
  +
s \ \ x \ \ y \ \ z \ \ x^2 \ \ xy \ \ xz \ \ y^2 \ \ yz \ \ z^2
  +
</math>
  +
  +
So the density, <math>x^2</math>, can be exactly described.
  +
  +
* Spherical: Now the auxiliary basis consists of the functions:
  +
<math>
  +
s x y z
  +
</math>

Revision as of 15:47, 12 May 2009

CamCASP => Programming => Auxiliary basis sets

Spherical versus Cartesian GTOs

  • What are the advantages of the two basis types: Spherical and Cartesian?
  • Which one results in more accurate energies and why?

Spherical GTOs:

  1. Smaller in size
  2. More well-defined set of linear equations during DF
  3. Lower accuracy.
  4. Seemingly erratic behaviour in aug-cc-pVnZ series.
  5. Rotations made easy as we have the Wigner rotation matrices.

Cartesian GTOs:

  1. Larger in size
  2. DF equations seem OK.
  3. Higher accuracy of energies
  4. More consistent behaviour of aug-cc-pVnZ series.
  5. Don't have rotation matrices coded, but these have been derived.

Simple example: 1

Here's an example that demonstrates the inadequacy of the spherical GTOs. Consider a 1p1s (yes, in this order) main basis and a single occupied orbital. This would be the px orbital. The space the auxiliary basis needs to span is 1p1s x 1p1s = 1s1p1d. So let us construct this auxiliary basis and perform the DF. Here's the question: Can the density be exactly fitted using this auxiliary basis?

The answer depends on which basis type is used for the auxiliary basis.

  • Cartesian: The auxiliary basis consists of the following combinations:

<math>

 s \ \  x \ \  y \ \  z \ \  x^2  \ \  xy  \ \ xz  \ \  y^2  \ \  yz  \ \  z^2

</math>

So the density, <math>x^2</math>, can be exactly described.

  • Spherical: Now the auxiliary basis consists of the functions:

<math>

 s   x   y   z   

</math>