1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33

pub const SMALL: f64 = 1.0e-14;
pub static SG1RADII: [f64; 19] = [0.0, 1.0000, 0.5882,   //H, He
3.0769, 2.0513, 1.5385, 1.2308, 1.0256, 0.8791, 0.7692, 0.6838,  //2nd Period
4.0909, 3.1579, 2.5714, 2.1687, 1.8750, 1.6514, 1.4754, 1.3333];  //3rd Period

pub const BOHR: f64 = 0.52917721092;

pub static BRAGG0: [f64; 131] = [0.0,  // Ghost atom
0.35,                                     1.40,             // 1s
1.45, 1.05, 0.85, 0.70, 0.65, 0.60, 0.50, 1.50,             // 2s2p
1.80, 1.50, 1.25, 1.10, 1.00, 1.00, 1.00, 1.80,             // 3s3p
2.20, 1.80,                                                 // 4s
1.60, 1.40, 1.35, 1.40, 1.40, 1.40, 1.35, 1.35, 1.35, 1.35, // 3d
            1.30, 1.25, 1.15, 1.15, 1.15, 1.90,             // 4p
2.35, 2.00,                                                 // 5s
1.80, 1.55, 1.45, 1.45, 1.35, 1.30, 1.35, 1.40, 1.60, 1.55, // 4d
            1.55, 1.45, 1.45, 1.40, 1.40, 2.10,             // 5p
2.60, 2.15,                                                 // 6s
1.95, 1.85, 1.85, 1.85, 1.85, 1.85, 1.85,                   // La, Ce-Eu
1.80, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75,             // Gd, Tb-Lu
      1.55, 1.45, 1.35, 1.35, 1.30, 1.35, 1.35, 1.35, 1.50, // 5d
            1.90, 1.80, 1.60, 1.90, 1.45, 2.10,             // 6p
1.80, 2.15,                                                 // 7s
1.95, 1.80, 1.80, 1.75, 1.75, 1.75, 1.75,
1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75,
1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75,
            1.75, 1.75, 1.75, 1.75, 1.75, 1.75,
1.75, 1.75,
1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75, 1.75];

pub static LEBEDEV_NGRID: [usize; 33] = [1, 6, 14, 26, 38, 50, 74, 86, 110, 146, 170,
194,  230,  266,  302,  350,  434,  590,  770,  974, 1202, 1454,
1730, 2030, 2354, 2702, 3074, 3470, 3890, 4334, 4802, 5294, 5810];