Natural Bond Orbital (NBO) Analysis

Gaussian contains version 3.1 of the NBO program by F. Weinhold and coworkers. NBO analysis is based on a method for optimally transforming a given wave function into localized form, corresponding to the one-center ("lone pairs") and two-center ("bonds") elements of the chemist's Lewis structure picture. In NBO analysis, the input atomic orbital basis set is transformed via natural atomic orbitals (NAOs) and natural hybrid orbitals (NHOs) into natural bond orbitals (NBOs). The NBOs obtained in this fashion correspond to the widely used Lewis picture, in which two-center bonds and lone pairs are localized.

The Weinhold group offers updated versions of the NBO program. The corresponding web page offers a wealth of "natural" information and should be considered the most authoritative source of NBO-related information. The conceptual world of natural orbitals has been laid out recently in a very detailed manner in "Valency and Bonding - A Natural Bond Orbital Donor-Acceptor Perspective" by F. Weinhold, C. R. Landis, Cambridge University Press, 2005.

An online-manual for the 3.0 version of the NBO program as combined with MOPAC 6.0 explaining some of the underlying concepts and some of the keywords can be found here.

A full NBO analysis is obtained in Gaussian when using the POP=NBO keyword, while POP=NPA requests just the Natural Population Analysis (NPA) phase of NBO, in which atomic partial charges are obtained through summation over NAOs.

Details of the NBO analysis procedure can best be illustrated using a specific example and formaldehyde will in the following be used for this purpose. Results of the NBO analysis are given with respect to the standard orientation of the system as given at the beginning of the Gaussian output file. For formaldehyde in its RHF/STO-3G structure the following input file performs a standard NBO analysis at the RHF/STO-3G level:

#P HF/STO-3G scf=tight pop=nbo

HF/STO-3G//HF/STO-3G sp formaldehyde

0 1
C1
O2  1  r2
H3  1  r3  2  a3
H4  1  r3  2  a3  3  180.0

r2=1.21672286
r3=1.10137241
a3=122.73666566

 


The following output is obtaind for this example:

 ******************************Gaussian NBO Version 3.1******************************
             N A T U R A L   A T O M I C   O R B I T A L   A N D
          N A T U R A L   B O N D   O R B I T A L   A N A L Y S I S
 ******************************Gaussian NBO Version 3.1******************************
       /RESON  / : Allow strongly delocalized NBO set

 Analyzing the SCF density

 Job title: HF/STO-3G//HF/STO-3G sp formaldehyde

 Storage needed:       510 in NPA,       633 in NBO (   8388541 available)


Note that all NBO output is formatted to a maximum 80-character width for convenient display on a computer terminal. The NBO heading echoes any requested keywords (none for the present default case) and shows an estimate of the memory requirements (in double precision words) for the separate steps of the NBO process, compared to the total allocated memory. Increase the memory allocated to your current job with the "%mem" directive if the estimated NBO requests exceed the available storage.

The next four NBO output segments summarize the results of natural population analysis (NPA). The first segment is the main NAO table, as shown below:

 NATURAL POPULATIONS:  Natural atomic orbital occupancies
 
  NAO Atom No lang   Type(AO)    Occupancy      Energy
 ---------------------------------------------------------
   1    C  1  s      Cor( 1s)     2.00000     -11.05947
   2    C  1  s      Val( 2s)     1.08461      -0.26145
   3    C  1  px     Val( 2p)     0.90165      -0.04494
   4    C  1  py     Val( 2p)     0.99304       0.07032
   5    C  1  pz     Val( 2p)     0.85379       0.14739

   6    O  2  s      Cor( 1s)     1.99999     -20.15602
   7    O  2  s      Val( 2s)     1.79982      -1.09031
   8    O  2  px     Val( 2p)     1.09835      -0.11626
   9    O  2  py     Val( 2p)     1.91857      -0.38126
  10    O  2  pz     Val( 2p)     1.37055      -0.15404

  11    H  3  s      Val( 1s)     0.98982       0.02250

  12    H  4  s      Val( 1s)     0.98982       0.02250

For each of the 12 NAO functions, this table lists the atom to which the NAO is attached, the angular momentum type 'lang' (s, px, etc.), the orbital type (whether core, valence, or Rydberg, and a conventional hydrogen-type label), the orbital occupancy (number of electrons, or 'natural population' of the orbital), and the orbital energy (in Hartree). Taking NAO 5 as an example we see that this NAO describes the energetically highest lying atomic orbital located at carbon and has the character of a 2pz atomic orbital, occupied by 0.85379 electrons. The principal quantum numbers for the NAO labels (1s, 2s, 3s, etc.) are assigned on the basis of the energy order if a Fock matrix is available, or on the basis of occupancy otherwise. A message is printed warning of a 'population inversion' if the occupancy and energy ordering do not coincide.

The next segment is an atomic summary showing the natural atomic charges (nuclear charge minus summed natural populations of NAOs on the atom) and total core, valence, and Rydberg populations on each atom:

 Summary of Natural Population Analysis:

                                       Natural Population
               Natural   -----------------------------------------------
    Atom No    Charge        Core      Valence    Rydberg      Total
 -----------------------------------------------------------------------
      C  1    0.16692      2.00000     3.83308    0.00000     5.83308
      O  2   -0.18728      1.99999     6.18729    0.00000     8.18728
      H  3    0.01018      0.00000     0.98982    0.00000     0.98982
      H  4    0.01018      0.00000     0.98982    0.00000     0.98982
 =======================================================================
   * Total *  0.00000      3.99999    12.00001    0.00000    16.00000

                                 Natural Population
 --------------------------------------------------------
   Core                       3.99999 ( 99.9997% of   4)
   Valence                   12.00001 (100.0001% of  12)
   Natural Minimal Basis     16.00000 (100.0000% of  16)
   Natural Rydberg Basis      0.00000 (  0.0000% of  16)
 --------------------------------------------------------

    Atom No         Natural Electron Configuration
 ----------------------------------------------------------------------------
      C  1      [core]2s( 1.08)2p( 2.75)
      O  2      [core]2s( 1.80)2p( 4.39)
      H  3            1s( 0.99)
      H  4            1s( 0.99)


This table succinctly describes the molecular charge distribution in terms of NPA charges. The carbon atom of formaldyde, for example, is assigned a net NPA charge of +0.16692 e at this level. Note also the slightly less positive charge on H3 and H4. Next follows a summary of the populations in the Natural Minimal Basis (NMB) and the Natural Rydberg Basis (NRB) for the composite system, summed over atoms.

The natural populations are summarized as an effective valence electron configuration ("Natural Electron Configuration") for each atom. Although the occupancies of the atomic orbitals are non-integer in the molecular environment, the effective atomic configurations can be related to idealized atomic states in 'promoted' configurations. For example, the carbon atom in the above table is most nearly described by an 1s22s12p2.75 (or idealized sp2) electron configuration.

The next segments of the output summarize the results of the NBO analysis. The first segment reports on details of the search for an NBO natural Lewis structure:

NATURAL BOND ORBITAL ANALYSIS:

               Occupancies        Lewis Structure  Low  High
    Occ.       -----------------  ---------------  occ  occ
Cycle Thresh.  Lewis   Non-Lewis  CR  BD  3C  LP   (L)  (NL) Dev

===================================================================
 1(1)  1.90    15.90533 0.09467    2   4   0   2    0    0   0.04
 
Structure accepted: No low occupancy Lewis orbitals


Normally there is a single cycle in the NBO search. The table summarizes a variety of information for each cycle: the occupancy threshold for a 'good' pair in the NBO search; the total populations of Lewis and non-Lewis NBOs; the number of core (CR), 2-center bond (BD), 3-center bond (3C), and lone pair (LP) NBOs in the natural Lewis structure; the number of low-occupancy Lewis (L) and 'high-occupancy' (> 0.1e) non-Lewis (NL) orbitals; and the maximum deviation ('Dev') of any formal bond order from a nominal estimate (NAO Wiberg bond index) for the structure. If the latter exceeds 0.1, additional NBO searches are initiated (indicated by the parenthesized number under 'Cycle') for alternative Lewis structures. The Lewis structure is accepted if all orbitals of the formal Lewis structure exceed the occupancy threshold (default: 1.90 electrons).

Next follows a more detailed breakdown of the Lewis and non-Lewis occupancies into core, valence, and Rydberg shell contributions:

Core           3.99999 (100.000% of  4)
Valence Lewis 11.90534 ( 99.211% of 12)
================== ============================
Total Lewis 15.90533 ( 99.408% of 16)

Valence non-Lewis 0.09467 ( 0.592% of 16)
Rydberg non-Lewis 0.00000 ( 0.000% of 16)
================== ============================
Total non-Lewis 0.09467 ( 0.592% of 16)


This shows the general quality of the natural Lewis structure description in terms of the percentage of the total electron density (in the above case about 99.4 %). The table also exhibits the relatively important role of the valence non-Lewis orbitals relative to the extra-valence orbitals (e.g. Rydberg NBOs which occur when using a non-minimal basis set) in the slight departures from a localized Lewis structure model.

Next follows the main listing of NBOs, displaying the form and occupancy of the complete set of NBOs that span the input AO space:

(Occupancy) Bond orbital/ Coefficients/ Hybrids
---------------------------------------------------------------------------------
1. (1.99777) BD ( 1) C 1 - O 2
      ( 41.41%) 0.6435* C 1 s( 34.36%)p 1.91( 65.64%)
      0.0000 0.5862 0.0000 0.0000 0.8102
      1s     2s     2px    2py    2pz

      ( 58.59%) 0.7654* O 2 s( 24.04%)p 3.16( 75.96%)
      0.0000 0.4903 0.0000 0.0000 -0.8716

This first NBO corresponds to a sigma(C-O) bond with approximate composition of 0.6435 C(sp1.91) + 0.7654 O(sp3.16). The weights are obtained from the squares of the coefficients as (0.6435)2 = 0.4141, corresponding to 41.41 % localization on carbon C1. In a similar way the 58.59 % localization on oxygen O2 is obtained. Overall, this describes a polar sigma(C-O) bond.

2. (2.00000) BD ( 2) C 1 - O 2
      ( 45.08%) 0.6714* C 1 s( 0.00%)p 1.00(100.00%)
      0.0000 0.0000 1.0000 0.0000 0.0000

      ( 54.92%) 0.7411* O 2 s( 0.00%)p 1.00(100.00%)
      0.0000 0.0000 1.0000 0.0000 0.0000

3. (1.99522) BD ( 1) C 1 - H 3
      ( 51.57%) 0.7181* C 1 s( 32.82%)p 2.05( 67.18%)
      0.0000 0.5729 0.0000 0.7071 -0.4145

      ( 48.43%) 0.6959* H 3 s(100.00%)
      1.0000

4. (1.99522) BD ( 1) C 1 - H 4
      ( 51.57%) 0.7181* C 1 s( 32.82%)p 2.05( 67.18%)
      0.0000 0.5729 0.0000 -0.7071 -0.4145

      ( 48.43%) 0.6959* H 4 s(100.00%)
      1.0000

5. (2.00000) CR ( 1) C 1 s(100.00%)
      1.0000 0.0000 0.0000 0.0000 0.0000

6. (1.99999) CR ( 1) O 2 s(100.00%)
      1.0000 0.0000 0.0000 0.0000 0.0000

7. (1.99857) LP ( 1) O 2 s( 75.96%)p 0.32( 24.04%)
      0.0000 0.8716 0.0000 0.0000 0.4903

8. (1.91857) LP ( 2) O 2 s( 0.00%)p 1.00(100.00%)
      0.0000 0.0000 0.0000 1.0000 0.0000


NBO 7 and NBO 8 correspond to lone pairs on oxygen O2 with spz0.32 and with pure py character, respectively.


9. (0.00338) BD*( 1) C 1 - O 2
      ( 58.59%) 0.7654* C 1 s( 34.36%)p 1.91( 65.64%)
      0.0000 0.5862 0.0000 0.0000 0.8102

      ( 41.41%) -0.6435* O 2 s( 24.04%)p 3.16( 75.96%)
      0.0000 0.4903 0.0000 0.0000 -0.8716

10. (0.00000) BD*( 2) C 1 - O 2
     ( 54.92%) 0.7411* C 1 s( 0.00%)p 1.00(100.00%)
     ( 45.08%) -0.6714* O 2 s( 0.00%)p 1.00(100.00%)

11. (0.04564) BD*( 1) C 1 - H 3
     ( 48.43%) 0.6959* C 1 s( 32.82%)p 2.05( 67.18%)
     0.0000 -0.5729 0.0000 -0.7071 0.4145
     
     ( 51.57%) -0.7181* H 3 s(100.00%)
     -1.0000

12. (0.04564) BD*( 1) C 1 - H 4
     ( 48.43%) 0.6959* C 1 s( 32.82%)p 2.05( 67.18%)
     0.0000 -0.5729 0.0000 0.7071 0.4145
     
     ( 51.57%) -0.7181* H 4 s(100.00%)
     -1.0000


For each NBO (1-12), the first line of the printout shows the occupancy (between 0 and 2.0000 electrons) and unique label of the NBO. This label gives the type ("BD" for 2-center bond, "CR" for 1-center core pair, "LP" for 1-center valence lone pair, "RY*" for 1-center Rydberg, and "BD*" for 2-center antibond, the unstarred and starred labels corresponding to Lewis and non-Lewis NBOs, respectively), a serial number (1, 2,... if there is a single, double,... bond between the pair of atoms), and the atom(s) to which the NBO is affixed. The next lines summarize the natural atomic hybrids hA of which the NBO is composed, giving the percentage (100xcA2) of the NBO on each hybrid (in parentheses), the polarization coefficient cA, the atom label, and a hybrid label showing the sp composition (percentage s-character, p-character, etc.) of each hA.

The next segment summarizes the second-order perturbative estimates of 'donor- acceptor' (bond-antibond) interactions in the NBO basis:

Second Order Perturbation Theory Analysis of Fock Matrix in NBO Basis

Threshold for printing: 0.50 kcal/mol

                                               E(2)   E(j) -E(i) F(i,j)
     Donor NBO (i)        Acceptor NBO (j)   kcal/mol    a.u.     a.u.
=======================================================================
within unit 1
1. BD ( 1) C 1 - O 2 / 11. BD*( 1) C 1 - H 3   0.92      1.76    0.036
1. BD ( 1) C 1 - O 2 / 12. BD*( 1) C 1 - H 4   0.92      1.76    0.036
3. BD ( 1) C 1 - H 3 /  9. BD*( 1) C 1 - O 2   1.85      1.59    0.048
3. BD ( 1) C 1 - H 3 / 12. BD*( 1) C 1 - H 4   1.71      1.36    0.043
4. BD ( 1) C 1 - H 4 /  9. BD*( 1) C 1 - O 2   1.85      1.59    0.048
4. BD ( 1) C 1 - H 4 / 11. BD*( 1) C 1 - H 3   1.71      1.36    0.043
8. LP ( 2) O 2       / 11. BD*( 1) C 1 - H 3  29.85      1.04    0.159
8. LP ( 2) O 2       / 12. BD*( 1) C 1 - H 4  29.85      1.04    0.159


This analysis is carried out by examining all possible interactions between 'filled' (donor) Lewis-type NBOs and 'empty' (acceptor) non-Lewis NBOs, and estimating their energetic importance by 2nd-order perturbation theory. Since these interactions lead to loss of occupancy from the localized NBOs of the idealized Lewis structure into the empty non-Lewis orbitals (and thus, to departures from the idealized Lewis structure description), they are referred to as 'delocalization' corrections to the zeroth-order natural Lewis structure. For each donor NBO (i) and acceptor NBO (j), the stabilization energy E(2) associated with delocalization ("2e-stabilization") is estimated. As the heading indicates, entries are included in this table only when the interaction energy exceeds a default threshold of 0.5 kcal/mol. The strongest interaction is in this example identified for the interaction of lone pair orbital NBO 8 localized on O2 with the adjacent sigma*(C-H) bonds.

Next appears a condensed summary of the principal NBOs, showing the occupancy, orbital energy, and the qualitative pattern of delocalization interactions associated with each:

Natural Bond Orbitals (Summary):

                                                   Principal Delocalizations
       NBO                 Occupancy    Energy     (geminal,vicinal,remote)
==========================================================================
Molecular unit 1 (CH2O)
 1. BD ( 1) C 1 - O 2    1.99777      -1.09402      11(g),12(g)
 2. BD ( 2) C 1 - O 2    2.00000      -0.44319
 3. BD ( 1) C 1 - H 3    1.99522      -0.69465      9(g),12(g)
 4. BD ( 1) C 1 - H 4    1.99522      -0.69465      9(g),11(g)
 5. CR ( 1) C 1          2.00000     -11.05946
 6. CR ( 1) O 2          1.99999     -20.15611
 7. LP ( 1) O 2          1.99857      -0.98034
 8. LP ( 2) O 2          1.91857      -0.38126      11(v),12(v)
 9. BD*( 1) C 1 - O 2    0.00338       0.89756
10. BD*( 2) C 1 - O 2    0.00000       0.28199
11. BD*( 1) C 1 - H 3    0.04564       0.66156
12. BD*( 1) C 1 - H 4    0.04564       0.66156
----------------------------------------
      Total Lewis   15.90533 ( 99.4083%)
Valence non-Lewis    0.09467 (  0.5917%)
Rydberg non-Lewis    0.00000 (  0.0000%)
----------------------------------------
Total unit 1        16.00000 (100.0000%)
Charge unit 1        0.00000


Details of the NBO analysis can be specified with a number of additional keywords. These latter keywords are read from the input stream after the end of the Z-Matrix specification, provided that the POP=NBORead keyword is used. Some of the more helpful keywords are:

NBOSUM requests printing of the NBO summary table. This combines elements of the NBO table and 2nd-order perturbation theory analysis table in a convenient form for recognizing the principal delocalization patterns.

RESONANCE requests search for highly delocalized structures. The NBO search normally aborts when one or more Lewis NBOs has less than the default occupancy threshold of 1.90 electrons for a 'good' electron pair. When the RESONANCE keyword is activated, this threshold is successively lowered in 0.10 decrements to 1.50, and the NBO search repeated to find the best Lewis structure within each occupancy threshold. The program returns with the best overall Lewis structure (lowest total non-Lewis occupancy) found in these searches. This is useful for benzene and other highly delocalized molecules.

NOBOND requests that no bonds (2-center NBOs) are to be formed in the NBO procedure. The resulting NBOs will then simply be 1-center atomic hybrids. This is useful for highly ionic species.

3CBOND requests search for 3-center bonds. The normal default is to search for only 1- and 2-center NBOs. This is useful for diborane and other electron-deficient 'bridged' species.

NLMO computes and prints out the summary table of Natural Localized Molecular Orbitals (NLMOs). The latter are similar to Boys or Edmiston-Ruedenberg LMOs, but more efficiently calculated. This is useful for a 'semi-localized' description of an SCF or correlated wavefunction and is activated automatically by all keywords that pertain to NLMOs (e.g., AONLMO, SPNLMO, DIPOLE).

BNDIDX calculates and prints Wiberg bond indices.

A detailed overview of the performance of the NPA method in describing the charge distribution in water can be found in: F. Martin, H. Zipse, J. Comp. Chem. 2005, 26, 97 - 105. A copy of the pdf file is available here.