Derivatives wave functions

Every one of the four products represents an AO configuration. The first describes both electrons on nucleus a] the second, both on nucleus b and the two others, one electron on each of a and b. The first two terms describe ionic situations and might be expected to occur in a more accurately derived wave function with only very little weight. Ileitler and London simply used the function in which they were missing altogether, namely... 

The linear combination of atomic orbitals (LCAO) used as an interpolation method to derive the symmetry properties of energy bands in crystals [9] has as coefficients the values of plane waves, exp(ik.R), at the positions R of atoms in the unit cell. Rather than a LCAO the derived wave function therefore is a linear combination of Bloch sums with the full space-group symmetry of the crystal. Atomic orbitals, by contrast, are modes in a spherically symmetrical central Coulomb field. 

Figure A5. A schematic showing the comparison of the full electron wavefunction and the pseudopotential-derived wave function. The strongly bound core electrons are replaced by a smoother analytical function. This schematic was adapted from Payne et al.l3].

As mentioned in the introduction, the simplest way of approximately accounting for the geomehic or topological effects of a conical intersection incorporates a phase factor in the nuclear wave function. In this section, we shall consider some specific situations where this approach is used and furthermore give the vector potential that can be derived from the phase factor. 

Section VI shows the power of the modulus-phase formalism and is included in this chapter partly for methodological purposes. In this formalism, the equations of continuity and the Hamilton-Jacobi equations can be naturally derived in both the nonrelativistic and the relativistic (Dirac) theories of the electron. It is shown that in the four-component (spinor) theory of electrons, the two exha components in the spinor wave function will have only a minor effect on the topological phase, provided certain conditions are met (nearly nonrelativistic velocities and external fields that are not excessively large). 

The question of determination of the phase of a field (classical or quantal, as of a wave function) from the modulus (absolute value) of the field along a real parameter (for which alone experimental determination is possible) is known as the phase problem [28]. (True also in crystallography.) The reciprocal relations derived in Section III represent a formal scheme for the determination of phase given the modulus, and vice versa. The physical basis of these singular integral relations was described in [147] and in several companion articles in that volume a more recent account can be found in [148]. Thus, the reciprocal relations in the time domain provide, under certain conditions of analyticity, solutions to the phase problem. For electromagnetic fields, these were derived in [120,149,150] and reviewed in [28,148]. Matter or Schrodinger waves were... 

Projecting the nuclear solutions Xt( ) oti the Hilbert space of the electronic states (r, R) and working in the projected Hilbert space of the nuclear coordinates R. The equation of motion (the nuclear Schrddinger equation) is shown in Eq. (91) and the Lagrangean in Eq. (96). In either expression, the terms with represent couplings between the nuclear wave functions X (K) and X (R). that is, (virtual) transitions (or admixtures) between the nuclear states. (These may represent transitions also for the electronic states, which would get expressed in finite electionic lifetimes.) The expression for the transition matrix is not elementaiy, since the coupling terms are of a derivative type. 

This part of our chapter has shown that the use of the two variables, moduli and phases, leads in a direct way to the derivation of the continuity and Hamilton-Jacobi equations for both scalar and spinor wave functions. For the latter case, we show that the differential equations for each spinor component are (in the nearly nomelativistic limit) approximately decoupled. Because of this decoupling (mutual independence) it appears that the reciprocal relations between phases and moduli derived in Section III hold to a good approximation for each spinor component separately, too. For velocities and electromagnetic field strengths that ate nomrally below the relativistic scale, the Berry phase obtained from the Schrddinger equation (for scalar fields) will not be altered by consideration of the Dirac equation. 

U(qJ is referred to as an adiabatic-to-diabatic transformation (ADT) matrix. Its mathematical sbucture is discussed in detail in Section in.C. If the electronic wave functions in the adiabatic and diabatic representations are chosen to be real, as is normally the case, U(q ) is orthogonal and therefore has n n — l)/2 independent elements (or degrees of freedom). This transformation mabix U(qO can be chosen so as to yield a diabatic electronic basis set with desired properties, which can then be used to derive the diabatic nuclear motion Schrodinger equation. By using Eqs. (27) and (28) and the orthonormality of the diabatic and adiabatic electronic basis sets, we can relate the adiabatic and diabatic nuclear wave functions through the same n-dimensional unitary transformation matrix U(qx) according to... 

In the -electronic-state adiabatic representation involving real electronic wave functions, the skew-symmetiic first-derivative coupling vector mahix... 

Importantly for direct dynamics calculations, analytic gradients for MCSCF methods [124-126] are available in many standard quantum chemistiy packages. This is a big advantage as numerical gradients require many evaluations of the wave function. The evaluation of the non-Hellmann-Feynman forces is the major effort, and requires the solution of what are termed the coupled-perturbed MCSCF (CP-MCSCF) equations. The large memory requirements of these equations can be bypassed if a direct method is used [233]. Modem computer architectures and codes then make the evaluation of first and second derivatives relatively straightforward in this theoretical framework. 

In the derivation used here, it is clear that two approximations have been made—the configurations are incoherent, and the nuclear functions remain localized. Without these approximations, the wave function fonn Eq. (C.l) could be an exact solution of the Schrddinger equation, as it is in 2D MCTDH form (in fact is in what is termed a natural orbital form as only diagonal configurations are included [20]). 

If the reaction is elementary, there is only a single transition state between A and B. At this point the derivative of the total electronic wave function with respect to the reaction coordinate Qa b vanishes ... 

In modem quantum chemistry packages, one can obtain moleculai basis set at the optimized geometry, in which the wave functions of the molecular basis are expanded in terms of a set of orthogonal Gaussian basis set. Therefore, we need to derive efficient fomiulas for calculating the above-mentioned matrix elements, between Gaussian functions of the first and second derivatives of the Coulomb potential ternis, especially the second derivative term that is not available in quantum chemistry packages. Section TV is devoted to the evaluation of these matrix elements. 

Since the form of the electronic wave functions depends also on the coordinate p (in the usual, parametric way), the matrix elements (21) are functions of it too. Thus it looks at first sight as if a lot of cumbersome computations of derivatives of the electronic wave functions have to be carried out. In this case, however, nature was merciful the matrix elements in (21) enter the Hamiltonian matrix weighted with the rotational constant A, which tends to infinity when the molecule reaches linear geometry. This means that only the form of the wave functions, that is, of the matrix elements in (21), in the p 0 limit are really needed. In the above mentioned one-elecbon approximation... 

A convenience of electronic basis functions (53) is that they reduce at infinitesimal-amplitude bending to (28) with the same meaning of the angle 9 we may employ these asymptotic forms in the computation of the matrix elements of the kinetic energy operator and in this way avoid the necessity of carrying out calculations of the derivatives of the electronic wave functions with respect to the nuclear coordinates. The electronic part of the Hamiltonian is represented in the basis (53) by... 

By using Eq. (A. 17), the first derivative coupling for the complex electronic wave function assumes the fomi [4]... 

Drukker, K., Hammes-Schiffer, S. An analytical derivation of MC-SCF vibrational wave functions for the quantum dynamical simulation of multiple proton transfer reactions Initial application to protonated water chains. J. Chem. Phys. 107 (1997) 363-374. 

For large molecules, computation time becomes a consideration. Orbital-based techniques, such as Mulliken, Lowdin, and NBO, take a negligible amount of CPU time relative to the time required to obtain the wave function. Techniques based on the charge distribution, such as AIM and ESP, require a sig-nihcant amount of CPU time. The GAPT method, which was not mentioned above, requires a second derivative evaluation, which can be prohibitively expensive. 

When you perform a single point semi-empirical or ab initio calculation, you obtain the energy and the first derivatives of the energy with respect to Cartesian displacement of the atoms. Since the wave function for the molecule is computed in the process, there are a number of other molecular properties that could be available to you. Molecular properties are basically an average over the wave function of certain operators describing the property. For example, the electronic dipole operator is basically just the operator for the position of an electron and the electronic contribution to the dipole moment is... 

Just as the electrical behaviour of a real diatomic molecule is not accurately harmonic, neither is its mechanical behaviour. The potential function, vibrational energy levels and wave functions shown in Figure f.i3 were derived by assuming that vibrational motion obeys Hooke s law, as expressed by Equation (1.63), but this assumption is reasonable only... 

The spin part pl can be derived by labelling the electrons 1 and 2 and remembering that, in general, each can have an a or /i spin wave function giving four possible combinations a(l)P(2), P(l)a(2), a(l)a(2) and P(l)P(2). Because the first two are neither symmetric nor antisymmetric to the exchange of electrons, which is equivalent to the exchange of the labels 1 and 2, they must be replaced by linear combinations giving... 

---