Function with adapted symmetry

In the symmetry-adapted formulation, the 43- term no longer occurs because the d-orbital density contains a vertical mirror plane even if such a plane is absent in the point group. This is illustrated as follows. Point groups without vertical mirror planes differ from those with vertical mirror planes by the occurrence of both dlm+ and d(m functions, with m being restricted to n, the order of the rotation axis. But the coordinate system can be rotated around the main symmetry axis such that P4 becomes zero. As proof, we write the (p dependence as... 

For systems with high symmetry, in particular for atoms, symmetry properties can be used to reduce the matrix of the //-electron Hamiltonian to separate noninteracting blocks characterized by global symmetry quantum numbers. A particular method will be outlined here [263], to complete the discussion of basis-set expansions. A symmetry-adapted function is defined by 0 = 04>, where O is an Hermitian projection operator (O2 = O) that characterizes a particular irreducible representation of the symmetry group of the electronic Hamiltonian. Thus H commutes with O. This implies the turnover rule (0 > II 0 >) = (), which removes the projection operator from one side of the matrix element. Since the expansion of OT may run to many individual terms, this can greatly simplify formulas and computing algorithms. Matrix elements (0/x H ) simplify to (4 H v) or... 

To introduce the perturbation expansions corresponding to perturbation theories with different symmetry adaptation, it is useful to introduce a general concept of the interpolation function ve(t) defined such that86 ... 

In this equation, N is equal to the number of unit cells in the crystal. Note how the function in Eq. 5.27 is the same as that of Eq. 5.19 for cyclic tt molecules, if a new index is defined ask = liij/Na. Bloch sums are simply symmetry-adapted linear combinations of atomic orbitals. However, whereas the exponential term in Eq. 5.19 is the character of the yth irreducible representation of the cychc group to which the molecule belongs, in Eq. 5.27 the exponential term is related to the character of the Mi irreducible representation of the cychc group of infinite order (Albright, 1985). This, in turn, may be replaced with the infinite linear translation group because of the periodic boundary conditions. It turns out that SALCs for any system with translational symmetry are con-stmcted in this same manner. Thus, as with cychc tt systems, there should never be a need to use the projection operators referred to earher to generate a Bloch sum. 

The symmetry properties of spontaneous strains are most conveniently understood by referring to the irreducible representations and basis functions for the point group of the high symmetry phase of a crystal of interest. These are given in Table 2 for the point group Almmm as an example. Basis functions x + y ) and are associated with the identity representation and are equivalent to (ei + ei) and C3 respectively. This is the same as saying that both strains are consistent with Almmm symmetry e = ei). The shear strain e - ei) is equivalent to the basis function (x - y ) which is associated with the Big representation, the shear strain e is equivalent to xy (B2g) and shear strains e and e to xz, yz respectively (Eg). The combinations (ci + 62) and (ci - ei) are referred to as symmetry-adapted strains because they have the form of specific basis functions of the... 

This defines a set of equations for the mean field Hamiltonians HPF. These equations have to be solved self-consistently since the thermodynamic values within the angle brackets in (109) involve the mean field Hamiltonians // F. In principle, all // F can be different in practice, we impose symmetry relations. Therefore, we choose a unit cell, compatible with the symmetry of the lattice introduced in Section II,D, and we put Hpf equal to // F whenever P and P belong to the same sublattice. Moreover, we apply unit cell symmetry that relates the mean field Hamiltonians on different sublattices. By using the symmetry-adapted functions introduced in Section II,B, the latter symmetry can be imposed as follows. We select a set of molecules constituting the asymmetric part of the unit cell. Then we assign to all other molecules P Euler angles tip-through which the mean field. Hamiltonian of some molecule P in the asymmetric part has to be rotated in order to obtain HrF. As a result, we... 

At this point let us make a remark concerning the size of the basis. In order to obtain convergence, one must sometimes include (Briels et al., 1984) basis functions with high values of / and n. High values of l are needed in particular when the orientations of the molecules are fairly well localized. This leads to a rapidly increasing size of the basis. Two measures can be taken to simplify the problem. First, one can adapt the basis of molecule P to the site symmetry at P, which block-diagonalizes the secular problem. If this does not sufficiently reduce the problem, the mean field model Hamiltonian (96) can be further separated by writing... 

This factorization amounts to the statement thatEq. (2.25) breaks down into two separate linear systems, one for the determination of a orbitals, and the other for n orbitals. In the Hartree-Fock scheme, self-consistent field equations (SCF equations 2.25) have as solutions symmetry-adapted functions (i.e. in the case of planar unsaturated molecules symmetric or antisymmetric functions with respect ot the molecular plane), at least for closed-shell ground states iM8,20,2i) ... 

In fact, symmetry requirements on molecular orbitals introduce in a variational calculation certain constraints, which raise the total energy 18>. This problem, called the symmetry dilemma , has been studied for some n electron systems 19). It is not important for the present discussion because for a system of closed-shell type the NO s associated with a total wave function of correct symmetry are automatically symmetry-adapted. 

This form of the Hamiltonian shows explicitly the couplings between wave functions with different Mk values and makes possible to factorize occupation vectors in alpha- and beta-strings like done in non-relativistic Cl theory. The difference with non-relativistic theory is that calculations are not restricted to one value of Mk. Applied without further approximation the formalism gives therefore no dramatic reduction in operation count over the symmetry-adapted unrestricted scheme described in the previous section. An advantage of the formalism is, however, that it facilitates incorporation of the relativistic scheme in non-relativistic Cl or MCSCF implementations [35] and that the scheme gives a natural subdivision of the full Cl matrix. 

Combining these rotation functions with the usual vibrational basis functions provides a symmetry-adapted rotation-vibration basis. 

We can now determine the symmetry-adapted coordinates by applying the projection operators, but the results can be written down almost immediately by again using the criterion of overlap with central symmetry functions. The A g, T u, and Eg -I- l2g SALCs reflect the nodal patterns of central p, and d functions, respectively. The T g mode corresponds to the rotation and evidently consists of tangential displacements of ligands in the equator perpendicular to the rotation axis. Finally, the T2u is a buckling mode, which has the symmetry of central f orbitals, viz. 

Several of the above-described publications extracted rotational spectra from inverse Laplace transforms of imaginary-time autocorrelation functions, quantities readily calculated with RQMC. The utility of defining a larger set of correlation functions, so-called symmetry-adapted imaginary-time autocorrelation functions was explored in a recent paper [50]. Computational efficiency in the calculation of weak spectral features was demonstrated by a study of He-CO binary complex. Some preliminary results of an analysis of a recently observed satellite band in the IR spectrum of CO2 doped He clusters were presented. 

It is usually expedient to set up the symmetry-adapted s at the outset of any calculation, using the methods to be discussed presently, and the number of terms in the Cl expansion (3.1.4) will then be drastically reduced. It must be stressed that, in principle, there is no need to work with symmetry-adapted functions, and that the solutions of the full secular problem would automatically come out with correct symmetry properties in practice, however, the secular equations are so big that there are compelling reasons for seeking every possible reduction. 

Table 5.1 The density-matrix elements and molecular integrals for the hydrogen molecule in a symmetry-adapted basis of hydrogenic Is functions with exponents I (atomic units). Rows containing only zero elements and rows with elements that are related to those of other rows by permutational symmetry are not listed...

To identify the states which arise from a given atomic configuration and to construct properly symmetry-adapted determinental wave functions corresponding to these symmetries, one must employ L and Ml and S and Ms angular momentum tools. One first identifies those determinants with maximum Ms (this then defines the maximum S value that occurs) within that set of determinants, one then identifies the determinant(s) with maximum Ml (this identifies the highest L value). This determinant has S and L equal to its Ms and Ml values (this can be verified, for example for L, by acting on this determinant with f2 in the form... 

Molecular orbitals are not unique. The same exact wave function could be expressed an infinite number of ways with different, but equivalent orbitals. Two commonly used sets of orbitals are localized orbitals and symmetry-adapted orbitals (also called canonical orbitals). Localized orbitals are sometimes used because they look very much like a chemist s qualitative models of molecular bonds, lone-pair electrons, core electrons, and the like. Symmetry-adapted orbitals are more commonly used because they allow the calculation to be executed much more quickly for high-symmetry molecules. Localized orbitals can give the fastest calculations for very large molecules without symmetry due to many long-distance interactions becoming negligible. 

ADF uses a STO basis set along with STO fit functions to improve the efficiency of calculating multicenter integrals. It uses a fragment orbital approach. This is, in essence, a set of localized orbitals that have been symmetry-adapted. This approach is designed to make it possible to analyze molecular properties in terms of functional groups. Frozen core calculations can also be performed. 

---