The rotationally invariant bond order  between all pairs of atoms is printed. In this context a bond is defined as the sum of the squares of the density matrix elements connecting any two atoms. For ethane, ethylene, and acetylene, the carbon-carbon bond orders are roughly 1.00, 2.00 and 3.00 respectively. The diagonal terms are the valencies calculated from the atomic terms only and are defined as the sum of the bonds the atom makes with other atoms. In RHF calculations, the total density matrix (alpha plus beta density matrices) is perfectly duodempotent, that is, the square of the density matrix, P, is exactly two times the density matrix, or P*P = 2*P (see Bond Orders), and valencies will be correct. In UHF and non-variationally optimized wavefunctions the calculated valency will be incorrect, the degree of error being proportional to the non-duodempotency of the total density matrix. (In UHF work, the individual alpha and beta density matrices are idempotent, that is, P(alpha)*P(alpha) = 1.0*P(alpha) and P(beta)*P(beta) = 1.0*P(beta), but, in general, the sum of these two matrices is not duodempotent, i.e., (P(alpha)+P(beta))*(P(alpha)+P(beta)) ≠. 2.0*(P(alpha)+P(beta)).)
The bonding contributions of all M.O.s in the system are printed immediately before the bonds matrix. The idea of molecular orbital valency was developed by Gopinathan, Siddarth, and Ravimohan . Just as an atomic orbital has a 'valency', so does a molecular orbital. This leads to the following relations: The sum of the bonding contributions of all occupied M.O.s is the same as the sum of all valencies which, in turn is equal to two times the sum of all bonds. The sum of the bonding contributions of all M.O.s is zero.
In MOZYME calculations, only bonds of order greater than 0.01 are printed. If ALLBONDS is present, only bonds of order greater than 0.001 are printed.
If LARGE is present, then the Medrano-Bochicchio-Reale population analysis is printed. For each atom, the following quantities are generated:
Note that the last two terms are expressed in units of the electron, not the proton charge.
For more information, see [24,25,26,27], and also see J. A. Medrano, R. C. Bochicchio, S. G. Das, "The ROHF Extension of the Statistical Population Analysis of Electron and Spin Densities", J. Phys. B.