DFT method

The projector augmented-wave (PAW) DFT method was invented by Blochl to generalize both the pseudopotential and the LAPW DFT teclmiques [M]- PAW, however, provides all-electron one-particle wavefiinctions not accessible with the pseudopotential approach. The central idea of the PAW is to express the all-electron quantities in tenns of a pseudo-wavefiinction (easily expanded in plane waves) tenn that describes mterstitial contributions well, and one-centre corrections expanded in tenns of atom-centred fiinctions, that allow for the recovery of the all-electron quantities. The LAPW method is a special case of the PAW method and the pseudopotential fonnalism is obtained by an approximation. Comparisons of the PAW method to other all-electron methods show an accuracy similar to the FLAPW results and an efficiency comparable to plane wave pseudopotential calculations [, ]. PAW is also fonnulated to carry out DFT dynamics, where the forces on nuclei and wavefiinctions are calculated from the PAW wavefiinctions. (Another all-electron DFT molecular dynamics teclmique using a mixed-basis approach is applied in [84].)... 

Computational solid-state physics and chemistry are vibrant areas of research. The all-electron methods for high-accuracy electronic stnicture calculations mentioned in section B3.2.3.2 are in active development, and with PAW, an efficient new all-electron method has recently been introduced. Ever more powerfiil computers enable more detailed predictions on systems of increasing size. At the same time, new, more complex materials require methods that are able to describe their large unit cells and diverse atomic make-up. Here, the new orbital-free DFT method may lead the way. More powerful teclmiques are also necessary for the accurate treatment of surfaces and their interaction with atoms and, possibly complex, molecules. Combined with recent progress in embedding theory, these developments make possible increasingly sophisticated predictions of the quantum structural properties of solids and solid surfaces. 

To use direct dynamics for the study of non-adiabatic systems it is necessary to be able to efficiently and accurately calculate electronic wave functions for excited states. In recent years, density functional theory (DFT) has been gaining ground over traditional Hartree-Fock based SCF calculations for the treatment of the ground state of large molecules. Recent advances mean that so-called time-dependent DFT methods are now also being applied to excited states. Even so, at present, the best general methods for the treatment of the photochemistry of polyatomic organic molecules are MCSCF methods, of which the CASSCF method is particularly powerful. 

The first point to remark is that methods that are to be incorporated in MD, and thus require frequent updates, must be both accurate and efficient. It is likely that only semi-empirical and density functional (DFT) methods are suitable for embedding. Semi-empirical methods include MO (molecular orbital) [90] and valence-bond methods [89], both being dependent on suitable parametrizations that can be validated by high-level ab initio QM. The quality of DFT has improved recently by refinements of the exchange density functional to such an extent that its accuracy rivals that of the best ab initio calculations [91]. DFT is quite suitable for embedding into a classical environment [92]. Therefore DFT is expected to have the best potential for future incorporation in embedded QM/MD. 

Molecular dipole moments are often used as descriptors in QPSR models. They are calculated reliably by most quantum mechanical techniques, not least because they are part of the parameterization data for semi-empirical MO techniques. Higher multipole moments are especially easily available from semi-empirical calculations using the natural atomic orbital-point charge (NAO-PC) technique [40], but can also be calculated rehably using ab-initio or DFT methods. They have been used for some QSPR models. 

The molecular electronic polarizability is one of the most important descriptors used in QSPR models. Paradoxically, although it is an electronic property, it is often easier to calculate the polarizability by an additive method (see Section 7.1) than quantum mechanically. Ah-initio and DFT methods need very large basis sets before they give accurate polarizabilities. Accurate molecular polarizabilities are available from semi-empirical MO calculations very easily using a modified version of a simple variational technique proposed by Rivail and co-workers [41]. The molecular electronic polarizability correlates quite strongly with the molecular volume, although there are many cases where both descriptors are useful in QSPR models. 

In formulating a mathematical representation of molecules, it is necessary to define a reference system that is defined as having zero energy. This zero of energy is different from one approximation to the next. For ah initio or density functional theory (DFT) methods, which model all the electrons in a system, zero energy corresponds to having all nuclei and electrons at an infinite distance from one another. Most semiempirical methods use a valence energy that cor-... 

The conductor-like screening model (COSMO) is a continuum method designed to be fast and robust. This method uses a simpler, more approximate equation for the electrostatic interaction between the solvent and solute. Line the SMx methods, it is based on a solvent accessible surface. Because of this, COSMO calculations require less CPU time than PCM calculations and are less likely to fail to converge. COSMO can be used with a variety of semiempirical, ah initio, and DFT methods. There is also some loss of accuracy as a result of this approximation. 

If spin contamination is small, continue to use unrestricted methods, preferably with spin-annihilated wave functions and spin projected energies. Do not use spin projection with DFT methods. When the amount of spin contamination is more significant, use restricted open-shell methods. If all else fails, use highly correlated methods. 

Polarizabilities and hyperpolarizabilities have been calculated with semi-empirical, ah initio, and DFT methods. The general conclusion from these studies is that a high level of theory is necessary to correctly predict nonlinear optical properties. 

Ah initio methods can yield reliable, quantitatively correct results. It is important to use basis sets with diffrise functions and high-angular-momentum polarization functions. Hyperpolarizabilities seem to be relatively insensitive to the core electron description. Good agreement has been obtained between ECP basis sets and all electron basis sets. DFT methods have not yet been used widely enough to make generalizations about their accuracy. 

The HE, GVB, local MP2, and DFT methods are available, as well as local, gradient-corrected, and hybrid density functionals. The GVB-RCI (restricted configuration interaction) method is available to give correlation and correct bond dissociation with a minimum amount of CPU time. There is also a GVB-DFT calculation available, which is a GVB-SCF calculation with a post-SCF DFT calculation. In addition, GVB-MP2 calculations are possible. Geometry optimizations can be performed with constraints. Both quasi-Newton and QST transition structure finding algorithms are available, as well as the SCRF solvation method. 

LCAO (linear combination of atomic orbitals) refers to construction of a wave function from atomic basis functions LDA (local density approximation) approximation used in some of the more approximate DFT methods... 

LORG (localized orbital-local origin) technique for removing dependence on the coordinate system when computing NMR chemical shifts LSDA (local spin-density approximation) approximation used in more approximate DFT methods for open-shell systems LSER (linear solvent energy relationships) method for computing solvation energy... 

OPW (orthogonalized plane wave) a band-structure computation method P89 (Perdew 1986) a gradient corrected DFT method parallel computer a computer with more than one CPU Pariser-Parr-Pople (PPP) a simple semiempirical method PCM (polarized continuum method) method for including solvation effects in ah initio calculations... 

PW91 (Perdew, Wang 1991) a gradient corrected DFT method QCI (quadratic conhguration interaction) a correlated ah initio method QMC (quantum Monte Carlo) an explicitly correlated ah initio method QM/MM a technique in which orbital-based calculations and molecular mechanics calculations are combined into one calculation QSAR (quantitative structure-activity relationship) a technique for computing chemical properties, particularly as applied to biological activity QSPR (quantitative structure-property relationship) a technique for computing chemical properties... 

Recently, a third class of electronic structure methods have come into wide use density functional methods. These DFT methods are similar to ab initio methods in many ways. DFT calculations require about the same amount of computation resources as Hartree-Fock theory, the least expensive ab initio method. 

DFT methods are attractive because they include the effects of electron correlation—the fact that electrons in a molecular system react to one another s motion and attempt to keep out of one another s way—in their model. Hartree-Fock calculations consider this effect only in an average sense—each electron sees and... 

In the last few years, methods based on Density Functional Theory have gained steadily in popularity. The best DFT methods achieve significantly greater accuracy than Harttee-Fock theory at only a modest increase in cost (far less than MP2 for medium-size and larger molecular systems). They do so by including some of the effects of electron correlation much less expensively than traditional correlated methods. 

DFT methods compute electron correlation via general functionals of the electron density (see Appendix A for details). DFT functionals partition the electronic energy into several components which are computed separately the kinetic energy, the electron-nuclear interaction, the Coulomb repulsion, and an exchange-correlation term accounting for the remainder of the electron-electron interaction (which is itself... 

These SVWN5 results are somewhat fortuitous. Be careful not to overgeneralize from their agreement to experiment. We will see a different result in Exercise 6.7. Several other excerises will also include comparisons of DFT methods to Hartree-Fock theory, MP2 and other electron correlation methods. 

Optimizing water dimer can be challenging in general, and DFT methods are known to have difficulty with weakly-bound complexes. When your optimization succeeds, make sure that you have found a minimum and not a transition structure by verifying that there are no imaginary frequencies. In the course of developing this exercise, we needed to restart our initial optimization from an improved intermediate step and to use Opt=CalcAII to reach a minimum. 

Following on the work of Kohn and Sham, the approximate functionals employed by current DFT methods partition the electronic energy into several terms ... 

Pure DFT methods are defined by pairing an exchange functional with a correlation functional. For example, the well-known BLYP functional pairs Becke s gradient-corrected exchange functional with the gradient-corrected correlation functional of Lee, Yang and Parr. 

In general, DFT calculations proceed in the same way as Hartree-Fock calculations, with the addition of the evaluation of the extra term, This term cannot be evaluated analytically for DFT methods, so it is computed via numerical integration. 

The foundation for the use of DFT methods in computational chemistry was the introduction of orbitals by Kohn and Sham. 5 The main problem in Thomas-Fermi models is that the kinetic energy is represented poorly. The basic idea in the Kohn and Sham (KS) formalism is splitting the kinetic energy functional into two parts, one of which can be calculated exactly, and a small correction term. 

---