Gradients wave functions

This makes it desirable to define other representations in addition to the electronically adiabatic one [Eqs. (9)-(12)], in which the adiabatic electronic wave function basis set used in the Bom-Huang expansion (12) is replaced by another basis set of functions of the electronic coordinates. Such a different electronic basis set can be chosen so as to minimize the above mentioned gradient term. This term can initially be neglected in the solution of the / -electionic-state nuclear motion Schrodinger equation and reintroduced later using perturbative or other methods, if desired. This new basis set of electronic wave functions can also be made to depend parametrically, like their adiabatic counterparts, on the internal nuclear coordinates q that were defined after Eq. (8). This new electronic basis set is henceforth refened to as diabatic and, as is obvious, leads to an electronically diabatic representation that is not unique unlike the adiabatic one, which is unique by definition. 

The gradient of the PES (force) can in principle be calculated by finite difference methods. This is, however, extremely inefficient, requiring many evaluations of the wave function. Gradient methods in quantum chemistiy are fortunately now very advanced, and analytic gradients are available for a wide variety of ab initio methods [123-127]. Note that if the wave function depends on a set of parameters X], for example, the expansion coefficients of the basis functions used to build the orbitals in molecular orbital (MO) theory. 

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. 

For a quantum mechanical calculation, the single point calculation leads to a wave function for the molecular system and considerably more information than just the energy and gradient are available. In principle, any expectation value might be computed. You can get plots of the individual orbitals, the total (or spin) electron density and the electrostatic field around the molecule. You can see the orbital energies in the status line when you plot an orbital. Finally, the log file contains additional information including the dipole moment of the molecule. The level of detail may be controlled by the PrintLevel entry in the chem.ini file. 

Pseudo-Newton-Raphson methods have traditionally been the preferred algorithms with ab initio wave function. The interpolation methods tend to have a somewhat poor convergence characteristic, requiring many function and gradient evaluations, and have consequently primarily been used in connection with semi-empirical and force field methods. 

If there is more than one constraint, one additional multiplier term is added for each constraint. The optimization is then performed on the Lagrange function by requiring that the gradient with respect to the x- and A-variable(s) is equal to zero. In many cases the multipliers A can be given a physical interpretation at the end. In the variational treatment of an HF wave function (Section 3.3), the MO orthogonality constraints turn out to be MO energies, and the multiplier associated with normalization of the total Cl wave function (Section 4.2) becomes the total energy. 

The pseudopotential density-functional technique is used to calculate total energies, forces on atoms and stress tensors as described in Ref. 13 and implemented in the computer code CASTEP. CASTEP uses a plane-wave basis set to expand wave-functions and a preconditioned conjugate gradient scheme to solve the density-functional theory (DFT) equations iteratively. Brillouin zone integration is carried out via the special points scheme by Monkhorst and Pack. The nonlocal pseudopotentials in Kleynman-Bylander form were optimized in order to achieve the best convergence with respect to the basis set size. 5... 

The interaction between the third magnetic field gradient and each atom collapses the wave function a to either a or p with equal probabilities. 

Unlike the wave function, the electron density is an observable and can be measured experimentally, e. g. by X-ray diffraction. One of its important features is that at any position of an atom, p(r) exhibits a maximum with a finite value, due to the attractive force exerted by the positive charge of the nuclei. However, at these positions the gradient of the density has a discontinuity and a cusp results. This cusp is a consequence of the singularity ZA... 

Shepard R (1995) The analytic gradient method for configuration interaction wave functions. Yarkony DR (ed) In Modern electronic structure theory part I, World Scientific, Singapore, p 345... 

In this diabatic Schrodinger equation, the only terms that couple the nuclear wave functions Xd(R-/v) are the elements of the W RjJ and zd q%) matrices. The —(fi2/2p)W i(Rx) matrix does not have poles at conical intersection geometries [as opposed to W(2 ad(R>.) and furthermore it only appears as an additive term to the diabatic energy matrix cd(q>.) and does not increase the computational effort for the solution of Eq. (55). Since the neglected gradient term is expected to be small, it can be reintroduced as a first-order perturbation afterward, if desired. 

As discussed in Section II. A, the adiabatic electronic wave functions, a and / 1,ad depend on the nuclear coordinates R> only through the subset (which in the triatomic case consists of a nuclear coordinate hyperradius p and a set of two internal hyperangles this permits one to relate the 6D vector W(1)ad(Rx) to another one w(1 ad(q J that is 3D. For a triatomic system, let aIX = (a1 -. blk, crx) be the Euler angles that rotate the space-fixed Cartesian frame into the body-fixed principal axis of inertia frame IX, and let be the 6D gradient vector in this rotated frame. The relation between the space-fixed VRi and is given by... 

---