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 Options). 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.
Output controlled by the Population keyword includes:
- Molecular orbitals and orbital energies. By default, all orbitals are included, but the output can be limited to a specific orbital range with the Orbitals option.
- Atomic charge distribution. The total charge per fragment is also reported if applicable.
- Multipole moments: dipole through hexadecapole.
APT charges are also computed by default during vibrational frequency calculations [Cioslowski89].
Options Controlling Output File Contents
No orbitals are printed, and no population analysis is done. This is the default for all calculations using the ZIndo method.
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.
Same as the Regular population analysis, except that all orbitals are printed. This is the default for Guess=Only jobs.
Perform a population analysis at every optimization step rather than just the initial and final ones.
Perform a population analysis of the highest N occupied and lowest N virtual orbitals, including the atomic contributions to each MO (see the examples). 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.
Set the minimum contribution percentage to include in individual orbital population analysis. The default is 10.
Population Analysis Options
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.
Perform Hirshfeld population analysis [Hirshfeld77, Ritchie85, Ritchie87]. CM5 atomic charges [Marenich12] are computed along with Hirshfeld charges. CM5 is a synonym for Hirshfeld. The HirshfeldEE option requests interatomic electrostatic interactions as well.
Biorthogonalize unrestricted molecular orbitals in order to maximally align electron pairs.
Save biorthogonalized orbitals in the checkpoint file over the canonical MOs.
Requests Natural Transition Orbital analysis [Martin03] of a CI-Singles or TD-DFT excited state. Must be accompanied by Density=(Check,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=(Check,Transition=N) Guess=Only jobs after the initial excited state calculation. NTO is a synonym for this option.
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:
|$ g16 <<END||Run a Guess=Only job just to generate orbitals.|
|%OldChk=excited.chk||Work on copy of original checkpoint file.|
|%Chk=staten.chk||Checkpoint file for this state.|
|# Geom=AllCheck ChkBas Guess=Only|
|Density=(Check,Transition=n) Pop=SaveNTO||Repeat for each state.|
|END||End of Gaussian input file.|
Charge Save Options
Suppress orthogonalization of saved Pop=(SaveMixed,NoOrthogonalize) orbitals (primarily for visualizing raw NBO results).
Computed atomic charges can be stored on the checkpoint file for use in a later MM calculation with Geom=Check. The options Pop=SaveMulliken, Pop=SaveESP, Pop=SaveNPA, Pop=SaveCM5, and so on can be used to save the corresponding charges. In the case of a multilayer ONIOM calculation, only explicitly computed charges are saved by default: i.e., charges for atoms in the QM layer(s). Any atomic charges present in the input file will not be used and they will be replaced by the newly fitted charges saved by these options.
The additional option Uncharged will keep the atomic charges present in the input file and only fit charges on uncharged atoms in the QM layers(s). The option combinations Pop=(Uncharged, SaveMulliken), Pop=(Uncharged, SaveCM5) and so on will save both the original atomic charges plus the newly fitted charges for the atoms that were originally uncharged.
Natural Orbital-Related Options
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.
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.
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
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. The data file for Antechamber (the AMBER program for generating RESP charges) can be generated using Pop=MK IOp(6/50=1) and specifying the file name on a separate line at the end of the Gaussian input file.
Uses the MK fitting but using UFF radii, which are defined for the full periodic table.
Produce charges fit to the electrostatic potential at points selected according to the CHelp scheme [Chirlian87].
Produce charges fit to the electrostatic potential at points selected according to the CHelpG scheme [Breneman90].
Specifies the Hu, Lu, and Yang charge fitting method [Hu07].
Specifies the Hu, Lu, and Yang charge fitting method, but using Gaussian’s standard atomic densities instead of those of HLY. The authors of HLY only parametrized the atomic densities required for the model for the first 18 elements. This is an alternative version that uses the HLY fitting scheme but with Gaussian’s standard atomic densities, which are available for the entire periodic table. For systems which can be done either way, the difference in atomic charges is usually between 1% and 5%.
Include the RESP (restrained electrostatic potential) constraint [Bayly97, Cornell93, Cieplak95] in computing potential-derived charges. Specify the weight factor as its argument: Pop=(MK,Resp=N) applies a weight of N x 10-6 Hartrees to the squared charges when computing Merz-Singh-Kollman charges. Other electrostatic potential-derived charge schemes also accept this option. N defaults to 2.
When fitting charges to the potential, also fit a point dipole at each atomic center.
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.
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 (Version 3)
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.
Requests an NCS analysis of the diagonal tensor elements.
Requests an NCS analysis of all tensor components.
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 $.
Save natural bond orbitals in the checkpoint file (for later visualization).
Save natural localized molecular orbitals in the checkpoint file (for later visualization).
Save the NBOs for the occupied orbitals and the NLMOs for the unoccupied orbitals in the checkpoint file (for later visualization).
NBO-Related Options (Versions 6 and 7)
An interface to NBO versions 6 and 7 is provided. Pop=NPA6, Pop=NBO6, Pop=NBO6Read and Pop=NBO6Delete request Natural Population Analysis, full Natural Bond Orbital Analysis, full NBO with NBO input read from the input stream and NBO analysis of the effects of deletion of some interactions (respectively), using the separate NBO6 program via the external interface. The options NPA7, NBO7, NBO7Read and NBO7Delete request the same using NBO7. For NBO7, the additional option Pop=NEDA is used to perform Natural Energy Decomposition Analysis. The analysis uses the same input information about fragments as counterpoise calculations.
The script required to run NBO6/NBO7 and the program itself must be obtained from Frank Weinhold (nbo6.chem.wisc.edu).
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
# 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: 23 July 2019. [G16 Rev. C.01]