Functionals, wave function calculations

Luce T A and Bennemann K H 1998 Nonlinear optical response of noble metals determined from first-principles electronic structures and wave functions calculation of transition matrix elements P/rys. Rev. B 58 15 821-6... 

Floizwarth N A W, Matthews G E, Dunning R B, Tackett A R and Zeng Y 1997 Comparison of the pro]ector augmented-wave, pseudopotentlal and linearized augmented-plane-wave formalisms for density-functional calculations of solids Phys. Rev. B 55 2005... 

View the contour map m several planes to see the general Torm of the distiibiiiioii. As long as you don t alter the molecular coordinates, you don t need to repeat th e wave function calculation. Use the left mouse button and the IlyperChem Rotation or Translation tools (or Tool icons ) to change the view of amolecnle without changing its atomic coordinates. 

If a covalent bond is broken, as in the simple case of dissociation of the hydrogen molecule into atoms, then theRHFwave function without the Configuration Interaction option (see Extending the Wave Function Calculation on page 37) is inappropriate. This is because the doubly occupied RHFmolecular orbital includes spurious terms that place both electrons on the same hydrogen atom, even when they are separated by an infinite distance. 

See also Extending the Wave Function Calculation on page 37. 

The quantity, V(R), the sum of the electronic energy Egjg. computed in a wave function calculation and the nuclear-nuclear coulomb interaction V(R,R), constitutes a potential energy surface having 3N independent variables (the coordinates R). The independent variables are the coordinates of the nuclei but having made the Born-Oppenheimer approximation, we can think of them as the coordinates of the atoms in a molecule. 

Table 1.10. Coefficients of Wave Functions Calculated for Methyl Cation by the CNDO/2...

Several models are just different ways of representing the same physical phenomenon. For example, the lattiee-gas model, eapillary waves, and the density functional calculations presented above basically have the same view of the interface. However,... 

We should also mention that basis sets which do not actually comply with the LCAO scheme are employed under certain circumstances in density functional calculations, i. e., plane waves. These are the solutions of the Schrodinger equation of a free particle and are simple exponential functions of the general form... 

As noted above, in the reference model the dependence of the electron wave functions A and B on the nuclear coordinates was entirely neglected and the wave functions A and B for the isolated ions or the wave functions calculated for corresponding equilibrium nuclear configurations Qk0i and Qkof according to Eqs. (8) were usually used in the calculations. 

The second approach typically involves expanding the wave functions in terms of atomic or atomic-like orbitals. Frequently s- and p-symmetry functions suffice for silicon and impurities up through the third period. A minimal basis set for silicon would consist of four basis functions, one 5-function and three p-functions on each atom. Some approaches supplement the minimal basis either with more atomic-like functions or with additional types of functions, such as plane waves. Some calculations use only plane waves for the basis. 

An unambiguous identification of anomalous muonium with the bond-center site became possible based on pseudopotential-spin-density-functional calculations (Van de Walle, 1990). For an axially symmetric defect such as anomalous muonium the hyperfine tensor can be written in terms of an isotropic and an anisotropic hyperfine interaction. The isotropic part (labeled a) is related to the spin density at the nucleus, ip(0) [2 it is often compared to the corresponding value in vacuum, leading to the ratio i7s = a/Afee = j i (O) Hi/) / (O) vac- The anisotropic part (labeled b) describes the p-like contribution to the defect wave function. 

It is not possible to use normal AO basis sets in relativistic calculations The relativistic contraction of the inner shells makes it necessary to design new basis sets to account for this effect. Specially designed basis sets have therefore been constructed using the DKH Flamiltonian. These basis sets are of the atomic natural orbital (ANO) type and are constructed such that semi-core electrons can also be correlated. They have been given the name ANO-RCC (relativistic with core correlation) and cover all atoms of the Periodic Table.36-38 They have been used in most applications presented in this review. ANO-RCC are all-electron basis sets. Deep core orbitals are described by a minimal basis set and are kept frozen in the wave function calculations. The extra cost compared with using effective core potentials (ECPs) is therefore limited. ECPs, however, have been used in some studies, and more details will be given in connection with the specific application. The ANO-RCC basis sets can be downloaded from the home page of the MOLCAS quantum chemistry software (http //www.teokem.lu.se/molcas). 

Non-Bom-Oppenheimer wave functions calculated in this way look more like their Born-Oppenheimer counterparts in the smaller basis set limits, and thus a good starting guess for these may be taken from Bom-Oppenheimer calculations in the same basis. Thus we calculate the electronic part first (this requires much fewer basis functions than does a full non-Bom-Oppenhimer calculation) and then form the total basis function by multiplying each electronic portion by a guess for the nuclear portion ... 

By calculating A.U (R) and Al/ (i ) separately, we can straightforwardly calculate the total adiabatic correction V (R) for any isotopes of A and B. The adiabatic corrections are calculated by numerical differentiation of the multi-configurational self-consistent field (MCSCF) wave functions calculated with Dalton [23]. The nurnerical differentiation was performed with the Westa program developed 1986 by Agren, Flores-Riveros and Jensen [22],... 

Recall from Chapter 3 that the basic idea of a plane-wave DFT calculation is to express the electron density in functions of the form exp (/k r). Electrons associated with plane waves of this form have energy E = (hk)2/2m. As a result, once a DFT calculation has been performed, the electronic DOS can... 

In the remainder of this section, we give a brief overview of some of the functionals that are most widely used in plane-wave DFT calculations by examining each of the different approaches identified in Fig. 10.2 in turn. The simplest approximation to the true Kohn-Sham functional is the local density approximation (LDA). In the LDA, the local exchange-correlation potential in the Kohn-Sham equations [Eq. (1.5)] is defined as the exchange potential for the spatially uniform electron gas with the same density as the local electron density ... 

Here, is the distance between atoms i andj, C(/ is a dispersion coefficient for atoms i andj, which can be calculated directly from tabulated properties of the individual atoms, and /dampF y) is a damping function to avoid unphysical behavior of the dispersion term for small distances. The only empirical parameter in this expression is S, a scaling factor that is applied uniformly to all pairs of atoms. In applications of DFT-D, this scaling factor has been estimated separately for each functional of interest by optimizing its value with respect to collections of molecular complexes in which dispersion interactions are important. There are no fundamental barriers to applying the ideas of DFT-D within plane-wave DFT calculations. In the work by Neumann and Perrin mentioned above, they showed that adding dispersion corrections to forces... 

---