Wavefunctions electronic

Differentiation of the basis set form of SCF, the Hartree-Fock-Roothaan equation, is complicated by the fact that the Fock operator is itself dependent on the orbital set. But this complication is not difficult to deal with. We will use C to represent the matrix of orbital expansion coefficients, S to be the matrix (operator) of the overlap of basis functions, F to be the Fock operator matrix, and E to be the orbital eigenvalue (orbital energy) matrix. The equation to be differentiated is 

To solve this equation, and all higher ofder equations, it is helpful (1) to write the derivative of C in terms of the zeroeth-order orbital coefficients, e.g., Co CoUo , (2) to multiply on the left by Cj, and (3) to substitute according to the following definitions  

All higher derivative equations have the same form, and so the general equation to be solved is 

Notice in this expression how the t factors discussed above are used and that the summations over p, v, and t are restricted according to the comparison rules for the integer lists. (This is what makes the computational implementation fairly straightforward.) 

As already mentioned, correlated wavefunctions present a major challenge for differentiation. The expansion coefficients of the configurations (or substitutions) have derivatives with respea to the parameters, and these derivatives require the orbital derivatives, meaning the Uo matrices. Although the computational approaches involve extensive coding, the example of SCF differentiation illustrates how the problem may be worked out. In the case of Cl, it is just the basic Cl eigenequation that must be differentiated  

---

These limitations lead to electron spin multiplicity restrictions and to differing nuclear spin statistical weights for the rotational levels. Writing the electronic wavefunction as the product of an orbital fiinction and a spin fiinction there are restrictions on how these functions can be combined. The restrictions are imposed by the fact that the complete function has to be of synnnetry... 

The magnitude of the perturbations can be calculated fairly quantitatively from high-quality electronic wavefunctions including configuration interaction [24]. 

A quantum mechanical treatment of molecular systems usually starts with the Bom-Oppenlieimer approximation, i.e., the separation of the electronic and nuclear degrees of freedom. This is a very good approximation for well separated electronic states. The expectation value of the total energy in this case is a fiinction of the nuclear coordinates and the parameters in the electronic wavefunction, e.g., orbital coefficients. The wavefiinction parameters are most often detennined by tire variation theorem the electronic energy is made stationary (in the most important ground-state case it is minimized) with respect to them. The... 

Z-matriccs arc commonly used as input to quantum mechanical ab initio and serai-empirical) calculations as they properly describe the spatial arrangement of the atoms of a molecule. Note that there is no explicit information on the connectivity present in the Z-matrix, as there is, c.g., in a connection table, but quantum mechanics derives the bonding and non-bonding intramolecular interactions from the molecular electronic wavefunction, starting from atomic wavefiinctions and a crude 3D structure. In contrast to that, most of the molecular mechanics packages require the initial molecular geometry as 3D Cartesian coordinates plus the connection table, as they have to assign appropriate force constants and potentials to each atom and each bond in order to relax and optimi-/e the molecular structure. Furthermore, Cartesian coordinates are preferable to internal coordinates if the spatial situations of ensembles of different molecules have to be compared. Of course, both representations are interconvertible. 

It was stated above that the Schrodinger equation cannot be solved exactly for any molecular systems. However, it is possible to solve the equation exactly for the simplest molecular species, Hj (and isotopically equivalent species such as ITD" ), when the motion of the electrons is decoupled from the motion of the nuclei in accordance with the Bom-Oppenheimer approximation. The masses of the nuclei are much greater than the masses of the electrons (the resting mass of the lightest nucleus, the proton, is 1836 times heavier than the resting mass of the electron). This means that the electrons can adjust almost instantaneously to any changes in the positions of the nuclei. The electronic wavefunction thus depends only on the positions of the nuclei and not on their momenta. Under the Bom-Oppenheimer approximation the total wavefunction for the molecule can be written in the following form ... 

A determinant is the most convenient way to write down the permitted functional forms of a polv electronic wavefunction that satisfies the antisymmetry principle. In general, if we have electrons in spin orbitals Xi,X2, , Xn (where each spin orbital is the product of a spatial function and a spin function) then an acceptable form of the wavefunction is ... 

The most commonly employed tool for introducing such spatial correlations into electronic wavefunctions is called configuration interaction (Cl) this approach is described briefly later in this Section and in considerable detail in Section 6. 

Electronic Wavefunctions Must Also Possess Proper Symmetry. These Include Angular Momentum and Point Group Symmetries... 

We have extended the linear combination of Gaussian-type orbitals local-density functional approach to calculate the total energies and electronic structures of helical chain polymers[35]. This method was originally developed for molecular systems[36-40], and extended to two-dimensionally periodic sys-tems[41,42] and chain polymers[34j. The one-electron wavefunctions here are constructed from a linear combination of Bloch functions c>>, which are in turn constructed from a linear combination of nuclear-centered Gaussian-type orbitals Xylr) (in ihis case, products of Gaussians and the real solid spherical harmonics). The one-electron density matrix is given by... 

A basis set is the mathematical description of the orbitals within a system (which in turn combine to approximate the total electronic wavefunction) used to perform the theoretical calculation. Larger basis sets more accurately approximate the orbitals by imposing fewer restrictions on the locations of the electrons in space. In the true quantum mechanical picture, electrons have a finite probability of existing anywhere in space this limit corresponds to the infinite basis set expansion in the chart we looked at previously. 

Solving this equation for the electronic wavefunction will produce the effective nuclear potential function It depends on the nuclear coordinates and describes the potential energy surface for the system. 

For an electronic wavefunction, antisymmetry is a physical requirement following from the fact that electrons are fermions. It is essentially a requirement that y agree with the results of experimental physics. More specifically, this requirement means that any valid wavefunction must satisfy the following condition ... 

Multiplying a molecular orbital function by a or P will include electron spin as part of the overall electronic wavefunction i /. The product of the molecular orbital and a spin function is defined as a spin orbital, a function of both the electron s location and its spin. Note that these spin orbitals are also orthonormal when the component molecular orbitals are. 

We might reasonably expect that the electronic wavefunction would depend on the particular values of Ra and 7 b at which the nuclei were fixed, and I have indicated this in the expression above. 

The electronic wavefunction is thus given as solution of = ggiAe and the total energy is given by... 

As before, the nuclei are to be thought of as being clamped in position for the purpose of evaluating the electronic energy and electronic wavefunction. The electronic wavefunction depends implicitly on the nuclear coordinates, which is why I have shown the functional dependence. 

Imagine a model hydrogen molecule with non-interacting electrons, such that their Coulomb repulsion is zero. Each electron in our model still has kinetic energy and is still attracted to both nuclei, but the electron motions are completely independent of each other because the electron-electron interaction term is zero. We would, therefore, expect that the electronic wavefunction for the pair of electrons would be a product of the wavefunctions for two independent electrons in H2+ (Figure 4.1), which I will write X(rO and F(r2). Thus X(ri) and T(r2) are molecular orbitals which describe independently the two electrons in our non-interacting electron model. 

I ain going to leave you to prove for yourself that the wavefunction corresponding to this infinite-distance H2 problem is a product of two hydrogen atom wavefunc-tions. Physically, you might have expected this the two atoms are independent so the electronic wavefunctions multiply to give the molecular electronic... 

If I write possible atomic orbitals for hydrogen atom A as Xa possible atomic orbitals for hydrogen atom B as the molecular electronic wavefunction will be... 

Electrons are indistinguishable, they simply cannot be labelled. This means that an acceptable electronic wavefunction has to treat all electrons on an equal footing. Thus, although 1 have so far implied that electron 1 is to be associated with nucleus Ha, and electron 2 with nucleus Hb, 1 must also cater for the alternative description where electron 1 is associated with nucleus Hb and electron 2 with nucleus Ha. 1 therefore have to modify Table 4.1 to Table 4.2. 

We can construct a total electronic wavefunction as the product of a spatial part and a spin part. For the electronic ground state of H2 we can consider combinations of the two spatial terms... 

As computational facilities improve, electronic wavefunctions tend to become more and more complicated. A configuration interaction (Cl) calculation on a medium-sized molecule might be a linear combination of a million Slater determinants, and it is very easy to lose sight of the chemistry and the chemical intuition , to say nothing of the visualization of the results. Such wavefunctions seem to give no simple physical picture of the electron distribution, and so we must seek to find ways of extracting information that is chemically useful. 

The total electronic wavefunction is the product of a spatial part and a spin part it is it(r) times a(s) or /3(s) for this one-electron molecule. There are thus two different quantum states having the same spatial part i/r(r). In the absence of a magnetic field, these are degenerate. 

What happens when we have a many-electron wavefunction, such as the one below which relates to the simple valence-bond treatment of dihydrogen ... 

Integration of P(r) with respect to the coordinates of this electron (now written r) gives the number of electrons, 2 in this case. In the case of a many-electron wavefunction that depends on the spatial coordinates of electrons 1,2,..., m, we define the electron density as... 

For every electronic wavefunction that is an eigenfunction of the electron spin operator S, the one-electron density function always comprises an spin part... 

Once electron repulsion is taken into account, this separation of a many-electron wavefunction into a product of one-electron wavefunctions (orbitals) is no longer possible. This is not a failing of quanmm mechanics scientists and engineers reach similar conclusions whenever they have to deal with problems involving three or more mutually interacting particles. We speak of the three-body problem. 

There are m doubly occupied molecular orbitals, and the number of electrons is 2m because we have allocated an a and a spin electron to each. In the original Hartree model, the many-electron wavefunction was written as a straightforward product of one-electron orbitals i/p, i/ and so on... 

Don t confuse the state wavefunction with a molecular orbital we might well want to build the state wavefunction, which describes all the 16 electrons, from molecular orbitals each of which describe a single electron. But the two are not the same. We would have to find some suitable one-electron wavefunctions and then combine them into a slater determinant in order to take account of the Pauli principle. 

---