Structure du fichier définissant le système

Un fichier est absolumment nécessaire : la définition du système à étudier. Le nom de fichier n'a pas d'importance, puisque la suite est produite par la commande

geninput < fichier_du_systeme

Voilà alors sa structure (en format libre)

1) NATOM: nombre d'atomes dans la molécule

Pour chaque atome on donne des cartes avec

2) NAT, NSHL, POS: numéro atomique, nombre de couches atomiques, position

Pour chaque couche :

3) ISTYPE, NPRIM: type de couche (0,1,2 pour $s,p,d$) et nombre de primitives dans cette couche.

Pour chaque primitive :

4) EXX, COEFF: Exposant et coefficient

Comme exemple peut servir NH$_3$ dans une petite base

4                                  4 atoms in this molecule
 7 6   -2.481795969   -0.023514305   0.000006953    the N
  0 5
 4177.5687     0.001829476
 627.23610     0.013995972
 142.77926     0.06836863
 40.271825     0.23023378
 12.893158     0.46369183
  0  1
 4.4175095     0.35840086
  0  1
 .76130776    -0.54135141
  0  1
 .22402545    -0.58812336
  1  3
 14.089373     0.037761819
 3.0322858     0.20922098
 .82469052     0.50670505
  1  1
 .22435351     0.47458390
 1 2  -2.865352733    1.859673616   -0.000029210    the first  H
  0  3
 13.361500     0.01906000
 2.0133000     0.13424000
 .45375700     0.47449000
  0  1
 .12331700      1.0
 1 2  -3.312917613   -0.792592254    1.552599025    the second H
  0  3
 13.361500     0.01906000
 2.0133000     0.13424000
 .45375700     0.47449000
  0  1
 .12331700     1.0
 1 2  -3.312654127   -0.792598221   -1.552723180    the third H
  0  3
 13.361500     0.01906000
 2.0133000     0.13424000
 .45375700     0.47449000
  0  1
 .12331700      1.0

Il est nécessaire d'inclure une base séparée pour chaque atome, même si on utilise la même base pour plusieurs atomes. Des pseudopotentiels ne sont pas utilisables.