Charge Models#

Minimal Basis Iterative Stockholder (MBIS)#

The Minimal Basis Iterative Stockholder (MBIS) [Ayers2016] method is a variant of the Hirshfeld method. The idea behind this approach is to expand the pro-atomic density in a minimal set of atom-centered s-type Slater functions. MBIS is used in various applications like Force-Field parameterization or rescaling atomic polarizabilities.

MBIS charges can easily be requested by the simple input keyword MBIS:

!wB97X-V DEF2-TZVP MBIS

Within ORCA its also possible to compute atomic dipole-, quadrupole-, and octupolemoments within the MBIS scheme. These can be printed via the extended output option:

%method
MBIS_LARGEPRINT TRUE  # default = FALSE 
end

There are various further options to use with MBIS that can be found in the ORCA manual.

Example: 5,6-diaminobenzene-2,3-dicarbonitrile#

In this example, we will compute MBIS properties for 5,6-diaminobenzene-2,3-dicarbonitrile. This molecules has, for a neutral compound, a high experimental dipole moment of 10.5±0.2 Debye [Müllen2016].

../_images/mbis-example-2.png

To obtain the MBIS charges and atomic multipoles we can employ the following input:

!wB97X-V DEF2-TZVP MBIS

%method
MBIS_LARGEPRINT TRUE  # default = FALSE
end

We can now find the MBIS charges in the output file basename.out:

------------------
MBIS ANALYSIS
------------------

Convergence threshold (charges)      ...        1.000e-06
Number of iterations                 ...        74

Total integrated alpha density       ...     52.999978264
Total integrated beta density        ...     52.999978264

  ATOM     CHARGE    POPULATION     SPIN
   0 N   -0.843178    7.843178    0.000000
   1 C    0.263789    5.736211    0.000000
   2 C    0.263801    5.736199    0.000000
   3 N   -0.843087    7.843087    0.000000
   4 C   -0.257154    6.257154    0.000000
   5 H    0.386423    0.613577    0.000000
   6 C   -0.027760    6.027760    0.000000
   7 C    0.356352    5.643648    0.000000
   8 C   -0.027763    6.027763    0.000000
   9 C    0.356292    5.643708    0.000000
  10 C   -0.257237    6.257237    0.000000
  11 H    0.386406    0.613594    0.000000
  12 H    0.384292    0.615708    0.000000
  13 N   -0.433365    7.433365    0.000000
  14 N   -0.433376    7.433376    0.000000
  15 H    0.384314    0.615686    0.000000
  16 H    0.170620    0.829380    0.000000
  17 H    0.170678    0.829322    0.000000

  TOTAL   0.000048   81.999952    0.000000

The charges mapped to to respective atoms are:

../_images/mbis-example-2-charges.png

As we requested the large printout we also find the atomic dipole-

 MBIS ATOMIC DIPOLE MOMENT (A.U.):
  ATOM           X             Y             Z
   0 N      -0.038066     -0.015118     -0.103029
[...]

quadrupole-, and

 MBIS ATOMIC QUADRUPOLE MOMENT (A.U.):
  ATOM          XX            YY            ZZ            XY            XZ            YZ
   0 N      -5.193059     -5.361218     -5.246343     -0.051584      0.008768      0.095351
[...]

octupole moments

 MBIS ATOMIC OCTUPOLE MOMENT (A.U.):
  ATOM          XXX           YYY           ZZZ           XXY           XXZ           XYY           XYZ           XZZ           YYZ           YZZ
   0 N      -0.327912      0.056948     -0.069312     -0.055253     -0.039152     -0.006909      0.012998      0.314205     -0.394240     -0.035321
[...]

Note

Note that there are various other charge models available within ORCA including Mulliken, Hirshfeld, and more.

Starting structure#

5,6-diaminobenzene-2,3-dicarbonitrile
18

  N       1.34721241437119     -0.19875204578868      2.39728170838216
  C       0.70867203692614     -0.08409568583235      1.17165794089898
  C      -0.69557523765992     -0.03503410755084      1.15547373801566
  N      -1.36695245837416     -0.05357620226644      2.36862799274125
  C      -1.35463448913564      0.08033827321651     -0.05146262039998
  H       2.34810629260485     -0.10903909384438      2.36022439042418
  C      -0.65492341506459      0.11650354440257     -1.24923904175946
  C      -1.38024608777480      0.22912849642218     -2.47568708881038
  C       0.73267914111713      0.02684352514116     -1.23535246715296
  C       1.49072510922012      0.04700848828691     -2.44697511349353
  C       1.40013842906872     -0.06847776855218     -0.02246437719730
  H      -2.36630779122682     -0.13781135092574      2.29511415417341
  H       0.95194597461754      0.38337577789381      3.11975943676764
  N      -1.98904585311221      0.31503849289177     -3.44573900381660
  N       2.12487320894712      0.06661771348260     -3.40431289691542
  H      -0.98996937874613     -0.71100834890760      3.03401439822709
  H      -2.43381979990049      0.12883159310508     -0.07188432898727
  H       2.47966190412196     -0.11602130117439     -0.01964682109748