Density functional models

LS. In the LS phase the molecules are oriented normal to the surface in a hexagonal unit cell. It is identified with the hexatic smectic BH phase. Chains can rotate and have axial symmetry due to their lack of tilt. Cai and Rice developed a density functional model for the tilting transition between the L2 and LS phases [202]. Calculations with this model show that amphiphile-surface interactions play an important role in determining the tilt their conclusions support the lack of tilt found in fluorinated amphiphiles [203]. 

It is important to realize that whenever qualitative or frontier molecular orbital theory is invoked, the description is within the orbital (Hartree-Fock or Density Functional) model for the electronic wave function. In other words, rationalizing a trend in computational results by qualitative MO theory is only valid if the effect is present at the HF or DFT level. If the majority of the variation is due to electron correlation, an explanation in terms of interacting orbitals is not appropriate. 

Zhao, Y. Lynch, B. J. Truhlar, D. G. Development and assessment of a new hybrid density functional model for thermochemical kinetics. J. Phys. Chem. A 2004, 14, 2715-2719. 

Halls, M. D., Schlegel, H. B., 1998, Comparison of the Performance of Local, Gradient-Corrected, and Hybrid Density Functional Models in Predicting Infrared Intensities , J. Chem. Phys., 109, 10587. 

Clementi, E., and S. J. Chakravorty. 1990. Comparative study of density functional models to estimate molecular atomization energies. J. Chem. Phys. 93, 2591. 

Erratum Application of the velocity-dissipation probability density function model to inhomogeneous turbulent flows [Phys. Fluids A 3, 1947 (1991)]. Physics of Fluids A Fluid Dynamics 4, 1088. 

Pozorski, J., and J. P. Minier. 1999. Probability density function modeling of dispersed two-phase turbulent flow. Phys. Rev. E 59 855-63. 

M. R. Nyden, An orthogonality constrained generalization of the Weizacker density functional model. J. Chem. Phys. 78, 4048-4051 (1983). 

The Gordon-Kim interaction functions may be compared with empirical potential functions derived by energy- or net-force minimization methods using known crystal structures. The O—O Gordon-Kim potentials are more repulsive, as illustrated in Fig. 9.2. Spackman points out that the empirical potentials likely contain a significant attractive component because of the inadequate allowance for electrostatic interactions in their derivation. This attractive component is included in the electrostatic interaction in the density functional model. 

An alternative approach to improve upon Hartree-Fock models involves including an explicit term to account for the way in which electron motions affect each other. In practice, this account is based on an exacf solution for an idealized system, and is introduced using empirical parameters. As a class, the resulting models are referred to as density functional models. Density functional models have proven to be successful for determination of equilibrium geometries and conformations, and are (nearly) as successful as MP2 models for establishing the thermochemistry of reactions where bonds are broken or formed. Discussion is provided in Section II. 

This chapter reviews models based on quantum mechanics starting from the Schrodinger equation. Hartree-Fock models are addressed first, followed by models which account for electron correlation, with focus on density functional models, configuration interaction models and Moller-Plesset models. All-electron basis sets and pseudopotentials for use with Hartree-Fock and correlated models are described. Semi-empirical models are introduced next, followed by a discussion of models for solvation. 

Kohn-Sham Equations and Density Functional Models... 

One approach to the treatment of electron correlation is referred to as density functional theory. Density functional models have at their heart the electron density, p(r), as opposed to the many-electron wavefimction, F(ri, r2,...). There are both distinct similarities and distinct differences between traditional wavefunction-based approaches (see following two sections) and electron-density-based methodologies. First, the essential building blocks of a many-electron wavefunction are single-electron (molecular) orbitals, which are directly analogous to the orbitals used in density functional methodologies. Second, both the electron density and the many-electron wavefunction are constructed from an SCF approach which requires nearly identical matrix elements. 

For his discovery, leading up to the development of practical density functional models, Walter Kohn was awarded the Nobel Prize in Chemistry in 1998. 

Despite the fact that numerical integration is involved, pseudoanalytical procedures have been developed for calculation of first and second energy derivatives. This means that density functional models, like Hartree-Fock models are routinely applicable to determination of equilibrium and transition-state geometries and of vibrational frequencies. 

Calculated equilibrium geometries for hydrogen and main-group hydrides containing one and two heavy (non-hydrogen) atoms are provided in Appendix A5 (Tables A5-1 and A5-10 for molecular mechanics models, A5-2 and A5-11 for Hartree-Fock models, A5-3 and A5-12 for local density models, A5-4 to A5-7 and A5-13 to A5-16 for BP, BLYP, EDFl and B3LYP density functional models, A5-8 and A5-17 for MP2 models and A5-9 and A5-18 for MNDO, AMI and PM3 semi-empirical models). Mean absolute errors in bond lengths are provided in Tables 5-1 and 5-2 for one and two-heavy-atom systems, respectively. 

B3LYP density functional models provide somewhat better bond length results than the other density functional models, generally very close to experimental distances and to those from MP2 calculations. As with the other density functional models, the errors are largest where one (or two) second-row elements are involved. This is apparent from Figure 5-4, which compares B3LYP/6-311+G and experimental bond distances. 

As was the case with hydrocarbons, 6-3IG and 6-311+G basis sets lead to similar bond lengths for all density functional models as well as for the MP2 model. This is reflected in the mean absolute errors. It is difficult to justify use of the larger basis set models for routine structure determinations. 

As with hydrocarbons, accurate descriptions of equilibrium structures for molecules with heteroatoms from density functional and MP2 models requires polarization basis sets. As shown in Table A5-20 (Appendix A5), bond distances in these compounds obtained from (EDF 1 and B3LYP) density functional models and from MP2 models... 

Calculated heavy-atom bond distances in molecules with three or more first and/or second-row atoms are tabulated in Appendix A5 molecular mechanics models (Table A5-21), Hartree-Fock models (Table A5-22), local density models (Table A5-23), BP, BLYP, EDFl and B3LYP density functional models (Tables A5-24 to A5-27), MP2 models (Table A5-28), and MNDO, AMI and PM3 semi-empirical models (Table A5-29). Results for STO-3G, 3-21G, 6-31G and 6-311+G basis sets are provided for Hartree-Fock models, but as in previous comparisons, only 6-3IG and 6-311+G basis sets are employed for local density, density functional and MP2 models. 

None of the semi-empirical models perform as well as Hartree-Fock models (except STO-3G), local density models, density functional models or MP2 models. PM3 provides the best overall description, although on the basis of mean absolute errors alone, all three models perform to an acceptable standard. Given the large difference in cost of application, semi-empirical models clearly have a role to play in structure determination. 

Results from local density models and BP, BLYP and EDF 1 density functional models are, broadly speaking, comparable to those from 6-3IG models, consistent with similarity in mean absolute errors. As with bond length comparisons, BLYP models stand out as inferior to the other non-local models. Both B3LYP/6-31G and MP2/6-31G models provide superior results, and either would appear to be a suitable choice where improved quality is required. 

Consistent with earlier remarks made for bond length comparisons, little if any improvement results in moving from the 6-3IG to the 6-311+G basis set for Hartree-Fock, local density and density functional models, but significant improvement results for MP2 models. 

Comparative data for heavy-atom bond lengths and skeletal bond angles for molecules incorporating one or more third or fourth-row, main-group elements are provided in Appendix A5 Table A5-39 for Hartree-Fock models with STO-3G, 3-2IG and 6-3IG basis sets. Table A5-40 for the local density model, BP, BLYP, EDFl andB3LYP density functional models and the MP2 model, all with the 6-3IG basis set, and in Table A5-41 for MNDO, AMI and PM3 semi-empirical models. 6-31G, local density, density functional and MP2 calculations have been restricted to molecules with third-row elements only. Also, molecular mechanics models have been excluded from the comparison. A summary of errors in bond distances is provided in Table 5-8. 

---