Minimum Basis Iterative Stockholder (MBIS) scheme
[1]:
from setup import prepare_argument_dict, prepare_grid_and_dens, print_results
from horton_part import MBISWPart
mol, grid, rho = prepare_grid_and_dens("data/h2o.fchk")
part = MBISWPart(**prepare_argument_dict(mol, grid, rho))
part.do_all()
print_results(part)
The number of electrons: 10.000003764139395
Coordinates of the atoms
[[ 0. 0. 0.224]
[-0. 1.457 -0.896]
[-0. -1.457 -0.896]]
Atomic numbers of the atom
[8 1 1]
================================================================================
Information of integral grids.
--------------------------------------------------------------------------------
Grid size of molecular grid: 18460
************************************ Atom 0 ************************************
|-- Radial grid size: 120
|-- Angular grid sizes:
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18
18 18 18 18 18 18 26 26 26 26 38 38 38 38 38 50 50 50 50 50
86 86 110 110 110 110 170 194 194 434 590 590 434 434 434 302 302 302 194 194
170 110 110 110 110 110 110 110 110 110 110 110 86 50 50 18 18 18 18 18
************************************ Atom 1 ************************************
|-- Radial grid size: 120
|-- Angular grid sizes:
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 18 18 18 18 18
18 18 18 18 18 18 18 18 18 26 26 26 38 38 38 38 38 38 38 50
50 50 50 50 50 86 86 86 110 110 110 170 170 170 194 194 194 194 170 170
170 110 110 110 110 110 110 110 110 110 110 86 86 50 50 50 18 18 18 18
************************************ Atom 2 ************************************
|-- Radial grid size: 120
|-- Angular grid sizes:
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 18 18 18 18 18
18 18 18 18 18 18 18 18 18 26 26 26 38 38 38 38 38 38 38 50
50 50 50 50 50 86 86 86 110 110 110 170 170 170 194 194 194 194 170 170
170 110 110 110 110 110 110 110 110 110 110 86 86 50 50 50 18 18 18 18
--------------------------------------------------------------------------------
================================================================================
Performing a density-based AIM analysis with a wavefunction as input.
Molecular grid : MolGrid
Using local grids : True
Iteration Change Entropy
1 7.03232e-01 6.49265e-01
2 7.33878e-02 1.10215e-01
3 1.63712e-02 8.93342e-02
4 1.53463e-02 8.03457e-02
5 1.17546e-02 7.51961e-02
6 8.91499e-03 7.21761e-02
7 6.81170e-03 7.03768e-02
8 5.24340e-03 6.92914e-02
9 4.06100e-03 6.86304e-02
10 3.16086e-03 6.82248e-02
11 2.47010e-03 6.79746e-02
12 1.93651e-03 6.78194e-02
13 1.52209e-03 6.77229e-02
14 1.19882e-03 6.76627e-02
15 9.45749e-04 6.76251e-02
16 7.47079e-04 6.76016e-02
17 5.90751e-04 6.75868e-02
18 4.67524e-04 6.75775e-02
19 3.70239e-04 6.75717e-02
20 2.93355e-04 6.75681e-02
21 2.32527e-04 6.75658e-02
22 1.84378e-04 6.75643e-02
23 1.46232e-04 6.75634e-02
24 1.16007e-04 6.75628e-02
25 9.20394e-05 6.75625e-02
26 7.30378e-05 6.75623e-02
27 5.79605e-05 6.75621e-02
28 4.60040e-05 6.75620e-02
29 3.65118e-05 6.75620e-02
30 2.89842e-05 6.75619e-02
31 2.30049e-05 6.75619e-02
32 1.82643e-05 6.75619e-02
33 1.44964e-05 6.75619e-02
34 1.15106e-05 6.75619e-02
35 9.13538e-06 6.75619e-02
36 7.25496e-06 6.75619e-02
37 5.75706e-06 6.75619e-02
38 4.57307e-06 6.75619e-02
39 3.62798e-06 6.75619e-02
40 2.88284e-06 6.75619e-02
41 2.28613e-06 6.75619e-02
42 1.81756e-06 6.75619e-02
43 1.44040e-06 6.75619e-02
44 1.14613e-06 6.75619e-02
45 9.07345e-07 6.75619e-02
Computing atomic populations.
Computing atomic charges.
Computing density decomposition for atom 0
Computing density decomposition for atom 1
Computing density decomposition for atom 2
Computing cartesian and pure AIM multipoles and radial AIM moments.
Storing proatom density spline for atom 0.
Storing proatom density spline for atom 1.
Storing proatom density spline for atom 2.
charges:
[-0.907 0.454 0.454]
cartesian multipoles:
[[-0.907 0. -0. 0.157 -5.377 -0. -0. -4.824 -0. -5.102 0. 0. 0.013 0. -0. 0. 0. 0.392 0. 0.302]
[ 0.454 0. 0.019 0.002 -0.277 -0. -0. -0.264 0.008 -0.266 0. -0.001 0.013 0. 0. 0. 0.008 0.011 0.004 0.042]
[ 0.454 0. -0.019 0.002 -0.277 0. -0. -0.264 -0.008 -0.266 0. 0.001 0.013 0. -0. 0. -0.008 0.011 -0.004 0.042]]
[ ]: