Population

DESCRIPTION

This properties keyword controls printing of molecular orbitals and several types of population analysis and atomic charge assignments. The default is to print just the total atomic charges and orbital energies, except for Guess=Only jobs, for which the default is Pop=Full (see below). Populations are done once for single-point calculations and at the first and last points of geometry optimizations. Note that the Population keyword requires an option.

The density that is used for the population analysis is controlled by the Density keyword. Note that only one density and method of charge fitting can be used in a job step. If several combinations are of interest, additional jobs steps can be added by specifying Guess=Only Density=Check, to avoid repeating any costly calculations.

Population analysis results are given in the standard orientation.

Output controlled by the Pop keyword includes:

APT charges are also computed by default during vibrational frequency calculations [Cioslowski89].

OPTIONS CONTROLLING OUTPUT FILE CONTENTS

None
No orbitals are printed, and no population analysis is done. This is the default for all calculations using the ZIndo method.

Minimal
Total atomic charges and orbital energies are printed. This is the default for all job types and methods except Guess=Only and/or ZIndo.

Regular
The five highest occupied and five lowest virtual orbitals are printed, along with the density matrices and a full (orbital by orbital and atom by atom) Mulliken population analysis. Since the size of the output depends on the square of the size of the molecule, it can become quite substantial for larger molecules.

Full
Same as the Regular population analysis, except that all orbitals are printed. This is the default for Guess=Only jobs.

Always
Perform a population analysis at every optimization step rather than just the initial and final ones.

Orbitals[=N]
Perform a population analysis of the highest N occupied and lowest N virtual orbitals (see the examples below). N defaults to 10. AllOrbitals may also be specified instead of Orbitals=N to request analysis of all orbitals. For open shell calculations, both alpha and beta orbitals are included.

ThreshOrbitals=N
Set the minimum contribution percentage to include in individual orbital population analysis. The default is 10.

POPULATION ANALYSIS OPTIONS

Bonding
Do a bonding population analysis in addition to the standard analysis. This is a Mulliken population analysis in which only density terms involving pairs of basis functions on different centers are retained.

MBS
Perform minimum basis set Mulliken population analysis [Montgomery99, Montgomery00]. NoMBS tells the program not to perform MBS Mulliken population analysis.

Hirshfeld
Perform Hirshfeld population analysis [Hirshfeld77, Ritchie85, Ritchie87]. The HirshfeldEE option requests interatomic electrostatic interactions as well.

Biorthogonalize
Biorthogonalize unrestricted molecular orbitals in order to maximally align electron pairs.

NaturalTransitionOrbitals
Requests Natural Transition Orbital analysis [Martin03] of a CI-Singles or TD-DFT excited state. Must be accompanied by Density=Transition=N or Density=(Read,Transition=N) in order to specify which transition density is to be used to generate the orbitals. To print the orbitals from several states of interest, run successive Pop=NTO Density=(Read,Transition=N) Guess=(Read,Only) jobs after the initial excited state calculation. NTO is a synonym for this option.

SaveNaturalTransitionOrbitals
Save the generated orbitals in the checkpoint file, replacing the canonical ones if the density was read-in from there. SaveNTO is a synonym for this option. If you want to visualize the orbitals, you need to write them back to the checkpoint file. It is a good idea to do so with a copy of the checkpoint file for each state. After the initial excited state calculation (using %Chk ex.chk), use a technique like the following to generate visualization data for each state:

$ cp ex.chk staten.chk                     Copy original checkpoint file.
$ g09 <<END                                Run a Guess=Only job just to generate orbitals.
%Chk=staten                                Use the checkpoint file copy.
# Geom=AllCheck ChkBas Guess=(Read,Only)
  Density=Transition=n Pop=SaveNTO         Repeat process for each state of interest.
END                                        End of Gaussian 09 input file.

NATURAL ORBITAL-RELATED OPTIONS

NaturalOrbitals
Do a natural orbital analysis of the total density. NO is a synonym for NaturalOrbitals.

NOAB
Do separate natural orbital analyses for the α and β densities. NaturalSpinOrbitals is a synonym for NOAB.

AlphaNatural
Do separate natural orbital analyses for the α and β densities, but store only the α densities for use in a .wfn file (see Output=WFN). NOA is a synonym for AlphaNatural.

BetaNatural
Do separate natural orbital analyses for the α and β densities, but store only the β densities for use in a .wfn file (see Output=WFN). NOB is a synonym for BetaNatural.

SpinNatural
Generate natural orbitals for the spin density (with α considered positive).

By default, natural orbitals are not included in the checkpoint file. Use a second job step of this form to place the natural orbitals into the checkpoint file:

--Link1--
%Chk=name
# Guess=(Save,Only,NaturalOrbitals) Geom=AllCheck ChkBasis

Run the formchk utility on the resulting checkpoint file to prepare the orbitals for visualization.

OPTIONS FOR GENERATING ELECTROSTATIC POTENTIAL-DERIVED CHARGES

MK
Produce charges fit to the electrostatic potential at points selected according to the Merz-Singh-Kollman scheme [Singh84, Besler90]. ESP and MerzKollman are synonyms for MK.

CHelp
Produce charges fit to the electrostatic potential at points selected according to the CHelp scheme [Chirlian87].

CHelpG
Produce charges fit to the electrostatic potential at points selected according to the CHelpG scheme [Breneman90].

Dipole
When fitting charges to the potential, constrain them to reproduce the dipole moment. ESPDipole is a synonym for Dipole.

AtomDipole
When fitting charges to the potential, also fit a point dipole at each atomic center.

ReadRadii
Read in alternative radii (in Angstroms) for each element for use in fitting potentials. These are read as pairs of atomic symbol and radius, terminated by a blank line.

ReadAtRadii
Read in alternative radii (in Angstroms) for each atom for use in fitting potentials. These are read as pairs of atom number and radius, terminated by a blank line.

NBO-RELATED OPTIONS

NBO
Requests a full Natural Bond Orbital analysis, using NBO version 3 [Foster80, Reed83a, Reed85, Reed85a, Carpenter87, Carpenter88, Reed88, Weinhold88].

NCS
Requests a partitioning of the NMR shielding tensors (computed using GIAOs) into magnetic contributions from bonds and lone pairs using the Natural Chemical Shielding analysis of Bohmann et al. [Bohmann97], which is based upon the NBO analysis method. By default, an analysis of the isotropic shielding is performed. NoNCS skips this analysis.

NCSDiag
Requests an NCS analysis of the diagonal tensor elements.

NCSAll
Requests an NCS analysis of all tensor components.

NPA
Requests just the Natural Population Analysis phase of NBO.

NBORead
Requests a full NBO analysis, with input controlling the analysis read from the input stream. Use this option to specify keywords for NBO. Refer to the NBO documentation for details on this input.

NBODel
Requests NBO analysis of the effects of deletion of some interactions. Only possible with SCF methods. Implies that NBO input will be read; refer to the NBO documentation for details. Note that NBO input starts in column 2 so that the UNIX shell does not interpret the initial $.

SaveNBOs
Save natural bond orbitals in the checkpoint file (for later visualization).

SaveNLMOs
Save natural localized molecular orbitals in the checkpoint file (for later visualization).

SaveMixed
Save the NBOs for the occupied orbitals and the NLMOs for the unoccupied orbitals in the checkpoint file (for later visualization).

RELATED KEYWORDS

Density, Output=WFN

EXAMPLES

Orbital-by-Orbital Population Analysis. The following route section requests population analysis of the lowest 3 virtual orbitals and the highest 3 occupied orbitals:

# UHF/6-311+G(d) Pop=Orbitals=3

Here is the resulting output from a calculation on FeO+ quartet:

 Atomic contributions to Alpha molecular orbitals:
 Alpha occ 16 OE=-0.923 is Fe1-d=1.00
 Alpha occ 17 OE=-0.699 is O2-p=0.88
 Alpha occ 18 OE=-0.690 is O2-p=0.68 Fe1-s=0.21
 Alpha vir 19 OE=-0.253 is Fe1-s=0.70 Fe1-p=0.27
 Alpha vir 20 OE=-0.188 is Fe1-p=0.71 O2-p=0.29
 Alpha vir 21 OE=-0.133 is Fe1-p=1.04
 
 Atomic contributions to Beta  molecular orbitals:
 Beta  occ 13 OE=-0.801 is O2-p=0.79
 Beta  occ 14 OE=-0.783 is Fe1-d=1.00
 Beta  occ 15 OE=-0.758 is O2-p=0.89
 Beta  vir 16 OE=-0.241 is Fe1-s=0.81 Fe1-p=0.17
 Beta  vir 17 OE=-0.139 is Fe1-p=0.91 Fe1-d=0.14

Note that both alpha and beta orbital information is included. For each orbital, the output reports the orbital energy (labeled OE and given in atomic units), followed by the fractional contribution of all basis functions of a given angular momentum for each relevant atom.

This is an example of a system where it is hard to tell the spin state of the system, because the canonical α and β orbitals are quite different. If you subsequently run a Guess=(Read,Only,BiOrthogonalize) Pop=Orbital calculation to analyze the results, then the program transforms the α and β orbitals to match up as much as possible (occupied and virtual separately). In this case, orbital energies for the transformed orbitals are nor given. Rather, the Pop=Orbital analysis reports the overlaps between the pairs of corresponding orbitals (i.e., α 19 with β 19), with a value of 1 indicating 100% correspondence. Instead of labeling the orbitals as occupied or virtual, they are labeled:

Docc

  

Doubly occupied: alpha and beta are occupied, with a match of 90% or better.

Asing, Bsing

  

Singly occupied: alpha without matching beta occupied, or lower occupieds which do not match up with any orbital of the opposite spin.

Dvir

  

Virtual orbital with nearly the same alpha and beta orbitals.

AVir, BVir

  

Virtual orbital which does not match up with any of the other opposite-spin virtuals.

The program lists the orbitals which match another orbital only once within the output. Also, for each unpaired excess alpha spin occupied orbital, there is always some beta virtual orbital which matches it, and the latter are also omitted.

Here is the output for FeO+ quartet, which has 3 more alpha electrons than beta but which turns out to really have 4 unpaired alpha occupieds and one unpaired beta occupied (the lowest doubly occupied and highest unoccupied orbitals have been cut from the output):

 Atomic contributions to molecular orbitals:
 Docc. orb 13 abOv=0.999 is Fe1-d=1.00
 Docc. orb 14 abOv=0.990 is O2-p=0.72 Fe1-s=0.14 Fe1-d=0.14
 Asing orb 15 abOv=0.316 is Fe1-d=0.99
 Asing orb 16 abOv=1.000 is Fe1-d=0.94
 Asing orb 17 abOv=1.000 is Fe1-d=0.91
 Asing orb 18 abOv=1.000 is Fe1-d=0.91
 Dvirt orb 19 abOv=1.000 is O2-s=1.92 Fe1-s=-0.67 O2-p=-0.43 Fe1-p=0.14
 Dvirt orb 20 abOv=1.000 is Fe1-s=0.52 Fe1-p=0.37
 Bsing orb 15 abOv=0.316 is O2-p=0.92

Orbitals 1-14 are doubly occupied. Alpha orbital 15 and beta orbital 15 are different singly occupied orbitals; all 4 unpaired alpha spins are on the Fe, and the one unpaired beta spin is on the oxygen. There are no unmatched alpha and beta virtuals within the default range of orbitals analyzed.

Fragment Level Decomposition. If fragment information is present, the output also reports populations over fragments. For Pop=Orbital jobs, the decomposition of each orbital by fragment is reported.

 Default output including fragment populations:
 Mulliken charges with hydrogens summed into heavy atoms:
              1
     1  Pd  -0.265855
     2  P    0.346314
     3  P    0.346314
     4  Cl  -0.168156
     5  Cl  -0.168156
     6  C    0.060982
     7  C    0.060982
     8  C   -0.106213
     9  C   -0.106213
 Sum of Mulliken charges with hydrogens summed into heavy atoms =   0.00000
          Condensed to fragments (all electrons):
     1      -0.265855
     2      -0.168156
     3      -0.168156
     4       0.060982
     5       0.060982
     6       0.480203
 
 Pop=Orbital output:
 Alpha occ 60 OE=-0.247 is Cl4-p=0.22 Cl5-p=0.22 P3-p=0.12 P2-p=0.12
                           Fr6=0.36 Fr2=0.22 Fr3=0.22

Sample NBO 3 Input. The following input file requests a bond order analysis using NBO 3:

# B3LYP/6-31G(d,p) Pop=NBORead 

Example of NBO bond orders

 0  1
C     0.000000    0.665676    0.000000
H     0.919278    1.237739    0.000000
H    -0.919239    1.237787    0.000000
C     0.000000   -0.665676    0.000000
H    -0.919278   -1.237739    0.000000
H     0.919239   -1.237787    0.000000

$nbo bndidx $end

Last updated on: 9 Jun 2009