Charges

Default Atomic Charge Analysis

By default BAND computes the following atomic charge analyses:

  • Hirshfeld Charges [1] [2]
  • Voronoi Deformation Charges (VDD, Voronoi Deformation Density)
  • Mulliken Charges (note: not calculated for Spin-Orbit calculations)
  • CM5 (Charge Model 5) [3] [4]

These atomic charges are printed to the output file and can be visualized using the AMSview GUI module.

A more detailed output of the atomic charges can be printed by specifying following print option (note: in Band 2017 and previous versions this detailed output was printed by default):

Print AtomicChargesDetails

Bader Analysis (AIM)

The QTAIM (Quantum Theory of Atoms in Molecules), also known as Bader Analysis can be enabled in the GridBasedAIM input block:

GridBasedAIM
   Enabled Yes/No
   Iterations integer
   SmallDensity float
   UseStartDensity Yes/No
End
GridBasedAIM
Type:Block
Description:Invoke the ultra fast grid based Bader analysis.
Enabled
Type:Bool
Default value:No
GUI name:Bader (AIM): Atomic properties
Description:Invoke the ultra fast grid based Bader analysis.
Iterations
Type:Integer
Default value:40
Description:The maximum number of steps that may be taken to find the nuclear attractor for a grid point.
SmallDensity
Type:Float
Default value:1e-06
Description:Value below which the density is ignored. This should not be chosen too small because it may lead to unassignable grid points.
UseStartDensity
Type:Bool
Default value:No
Description:Whether the analysis is performed on the startup density (True) or on the final density (False).
AIMCriticalPoints
   Enabled Yes/No
   EqvPointsTol float
   GridPadding float
   GridSpacing float
End
AIMCriticalPoints
Type:Block
Description:Compute the critical points of the density (Atoms In Molecules). The algorithm starts from a regular mesh of points, and from each of these it walks towards its corresponding critical point.
Enabled
Type:Bool
Default value:No
GUI name:: Critical points and bond paths
Description:Compute the critical points of the density (Atoms In Molecules). The algorithm starts from a regular mesh of points, and from each of these it walks towards its corresponding critical point.
EqvPointsTol
Type:Float
Default value:0.27
Unit:Bohr
Description:If the distance between two critical points is smaller than this value, the two critical points are considered to be the same point.
GridPadding
Type:Float
Default value:0.7
Unit:Bohr
Description:How much extra space is added to the starting guess domain in the search for the critical points
GridSpacing
Type:Float
Default value:0.5
Unit:Bohr
Description:The distance between the initial trial points.

Note

The Bader (AIM) analysis is performed on the fitted density (see ZlmFit). We advise to use a Good (or better) ZlmFit quality.

References

[1]F.L. Hirshfeld, Bonded-atom fragments for describing molecular charge densities, Theoretica Chimica Acta 44, 129 (1977)
[2]K.B. Wiberg and P.R. Rablen, Comparison of atomic charges derived via different procedures, Journal of Computational Chemistry 14, 1504 (1993)
[3]A.V. Marenich, S.V. Jerome, C.J. Cramer, D.G. Truhlar, Charge Model 5: An Extension of Hirshfeld Population Analysis for the Accurate Description of Molecular Interactions in Gaseous and Condensed Phases, Journal of Chemical Theory and Computation 8, 527 (2012)
[4]C.A. Peeples and G. Schreckenbach, Implementation of the SM12 Solvation Model into ADF and Comparison with COSMO, Journal of Chemical Theory and Computation 12, 4033 (2016)