F-matrix

The solution of any such eigenvalue problem requires a number of computer operations that scales as the dimension of the F matrix to the third power. Since the indices on the F matrix label AOs, this means... 

In the Huckel theory of simple hydrocarbons, one assumes that the election density on a carbon atom and the order of bonds connected to it (which is an election density between atoms) are uninfluenced by election densities and bond orders elsewhere in the molecule. In PPP-SCF theory, exchange and electrostatic repulsion among electrons are specifically built into the method by including exchange and electrostatic terms in the elements of the F matrix. A simple example is the 1,3 element of the matrix for the allyl anion, which is zero in the Huckel method but is 1.44 eV due to election repulsion between the 1 and 3 carbon atoms in one implementation of the PPP-SCF method. 

The elements of the F matr ix depend on either the charge densities q or the bond orders p, which in turn depend on the elements of the F matrix. This circular dependence means that we must start with some initial F matrix, calculate eigenvectors, use the eigenvectors to calculate q and p, which lead to new elements in the F matr ix, calculate new eigenvectors leading to a new F matrix, and so on, until repeated iteration brings about no change in the results. The job now is to fill in the elements of the F matr ix. 

Beyond these approximations as to what the aetual integrals in the F matrix are, in some more reeent semiempirieal methods further adjustment of any or all of the elements of the F matrix is used to bring the ealeulated results into the best possible agreement with standard thermoehemieal results, largely The number of... 

Hence, the off-diagonal F matrix elements vanish, Fj y =0 for X i and Xv both on the same atom (a). 

The off-diagonal F matrix elements coupling orbitals from different atoms (a and b) are expressed as... 

In the PPP method, only the n- orbitals and the eorresponding 7i-eleetrons are eonsidered. The parameters ineluded in the F matrix... 

The S matrix contains the overlap elements between basis functions, and the F matrix contains the Fock matrix elements. Each element contains two parts from the Fock operator (eq. (3.36)), integrals involving the one-electron operators, and a sum over... 

The latter equation is now in a standard form for determining the eigenvalues of the F matrix. The eigenvectors contained in C can then be backtransformed to the original coordinate system (C = S C )-... 

Fig. 3. (a) Thermal expansion coefficients a for the inclusion (f), matrix (m), mesophase (i) and composite (c) of a typical iron-epoxy particulate composite, with 5 percent volume fraction for the inclusions, versus temperature, (b) the reduced longitudinal expansion of the same elements, normalized to the unit-length versus temperature (diameter of inclusions df = 150 pm)... 

Before returning to the non-BO rate expression, it is important to note that, in this spectroscopy case, the perturbation (i.e., the photon s vector potential) appears explicitly only in the p.i f matrix element because this external field is purely an electronic operator. In contrast, in the non-BO case, the perturbation involves a product of momentum operators, one acting on the electronic wavefimction and the second acting on the vibration/rotation wavefunction because the non-BO perturbation involves an explicit exchange of momentum between the electrons and the nuclei. As a result, one has matrix elements of the form (P/ t)Xf > in the non-BO case where one finds lXf > in the spectroscopy case. A primary difference is that derivatives of the vibration/rotation functions appear in the former case (in (P/(J.)x ) where only X appears in the latter. 

The F matrix elements in eqs. (15) and (16) are formally the same as for closed-shell systems, the only difference being the definition of the density matrix in eq. (17), where the singly occupied orbital (m) has also to be taken into account. The total electronic energy (not including core-core repulsions) is given by... 

In the unrestricted Hartree-Fock method, a single-determinant wave function is used with different molecular orbitals for a and jS spins, and the eigenvalue problem is solved with separate F and F matrices. With the zero differential overlap approximation, the F matrix elements (25) become... 

It should be stated that eqs. (3) through (21) are general in that they do not require a rr-electron approximation and do not represent directly any specific computational method. In the rr-electron approach with Pople s approximations, the F matrix elements become those below. 

From the point of view of the present state of computational possibilities, an extension of open-shell methods to a and a + rr electronic systems is rather tempting. This extension is easy for the method of Longuet-Higgins and Pople e.g., the CNDO/2 method is amenable to radicals having one nondegenerate open shell if in the original terms for F matrix elements (66),... 

In the Roothaan method, within the CNDO approximation the elements of the F matrix and of various matrices contributing to the F matrix are (68) ... 

From the above example it is apparent that there may be (and usually are) more unknown force constants than observed frequencies. If no additional sources of data are available, it is necessary to make some assumptions to simplify the force field. Often all or some of the off-diagonal elements in the F matrix [Eq. (63)] are neglected, leading to the so-called valence force field (VFF) or modified valence force field (MVFF), respectively. 

Karas, M. Bachmann, D. Bahr, U. Hillenkamp, F. Matrix-assisted ultraviolet laser desorption of non-volatile compounds. Int. J. Mass Spectrom. Ion Process. 1987,78, 53. 

Nordhoff, E. Ingendoh, A. Cramer, R. Overberg, A. Stahl, B. Karas, M. Hillenkamp, F. Crain, P. F. Matrix-assisted laser desorption/ionization mass spectrometry of nucleic acids with wavelengths in the UV and IR. Rapid Comm. Mass Spectrom. 1992, 6,771-776. 

The subscript e means that the second derivatives are taken for internal parameters of the known or assumed molecular geometry.) It then follows that, in principle, to every molecule a F-matrix is assigned with molecule-specific force constants Fd.Transferability can only be expected for molecules which contain qualitatively and quantitatively very similar internal coordinates For instance, a force field which reproduces very well frequencies of n-al-kanes (from which it was derived) gives significantly worse results for branched alkanes... 

In case the Newton-Raphson process is used for energy minimisation, the subsequent evaluation of the molecular vibrational frequencies is a fairly simple matter since the F-matrix is available (17). 

The actual calculation consists of minimizing the intramolecular potential energy, or steric energy, as a function of the nuclear coordinates. The potential-energy expressions derive from the force-field concept that features in vibrational spectroscopic analysis according to the G-F-matrix formalism [111]. The G-matrix contains as elements atomic masses suitably reduced to match the internal displacement coordinates (matrix D) in defining the vibrational kinetic energy T of a molecule ... 

Thus, the solution to the classical problem with the Hamiltonian (energy) function given by Equation 3.19 is a set of 3N harmonic oscillators with 3N frequencies Aj = 4tt2v2. The A s result from the diagonalization of the F matrix of Equation 3.15. 

---