Vibrational normal coordinates
The potential energy is approximated by a second-order Taylor expansion around the stationary geometry x0. V x V xo x - xo 2 x - xo t IdV - Xo 1637 The energy for the expansion point, V x0 , may be chosen as zero, and the first derivative is zero since x0 is a stationary point. Here F is a 3Natom x 3Natom force constant matrix containing the second derivatives of the energy with respect to the coordinates. The nuclear Schr dinger equation for an Natom system is given by eq. 16.39 . Eq. 16.39 is...
Distance geometry methods
The idea in Distance Geometry methods is that trial geometries can be generated from a set of lower and upper bounds on distances between all pairs of atoms.55 The method was originally developed for generating possible geometries based on experimental information such as NMR NOE effects, which place certain constraints on the distance between protons that may be far from each other in terms of bonding. The bonding information itself, however, also places restrictions on distances between all...
Isogyric and Isodesmic Reactions
The most difficult part in calculating absolute stabilities heat of formation is the correlation energy. For calculating energies relative to isolated atoms, the goal of the Gn CBS models, essentially all the correlation energy of the bond being broken must be recovered. This in turn necessitates large basis sets and sophisticated correlation methods. This is also the reason why ab initio energies are not converted into heat of formation, as is normally done for semi-empirical methods eq. 3.97...
H T V V V
tot R, r n R e R, r 1.23 The above approximation corresponds to neglecting the coupling between the nuclear and electronic velocities, i.e. the nuclei are stationary from the electronic point of view. The electronic wave function thus depends parametrically on the nuclear coordinates, since it only depends on the position of the nuclei, not on their momentum. To a good approximation, the electronic wave function thus provides a potential energy surface upon which the nuclei move, and this...
Contracted Basis Sets
One disadvantage of all energy-optimized basis sets is the fact that they primarily depend on the wave function in the region of the inner-shell electrons. The 1s-electrons account for a large part of the total energy, and minimizing the energy will tend to make the basis set optimum for the core electrons, and less so for the valence electrons. However, chemistry is mainly dependent on the valence electrons. Furthermore, many properties for example polarizability depend mainly on the wave...
Classification of Basis Sets
Having decided on the type of function STO GTO and the location nuclei , the most important factor is the number of functions to be used. The smallest number of functions possible is a minimum basis set. Only enough functions are employed to contain all the electrons of the neutral atom s . For hydrogen and helium this means a single s-function. For the first row in the periodic system it means two s-functions 1s and 2s and one set of p-functions 2px, 2py and 2pz . Lithium and beryllium...
Onestructure interpolation methods coordinate driving linear and quadratic
The intuitively simple approach for locating a TS is to select one or a few internal reaction coordinates, i.e. those that describes the main difference between the reac-tant and product structures. A typical example is a torsional angle for describing a con-formational TS, or two bond distances for a bond breaking forming reaction. The selected coordinate s is are fixed at certain values, while the remaining variables are optimized, thereby adiabatically mapping the energy as a function of the...
B The Variational Principle
The Variational Principle states that an approximate wave function has an energy that is above or equal to the exact energy. The equality holds only if the wave function is exact. The proof is as follows. Assume that we know the exact solutions to the Schr dinger equation. H E, v i i 0,1,2, , - B.1 There are infinitely many solutions and we assume that they are labelled according to their energies, E0 being the lowest. Since the H operator is Hermitian, the solutions form a complete basis. We...
Slater and Gaussian Type Orbitals
There are two types of basis functions also called Atomic Orbitals AO , although they in general are not solutions to an atomic Schr dinger equation commonly used Introduction to Computational Chemistry, Second Edition. Frank Jensen. 2007 John Wiley amp Sons, Ltd in electronic structure calculations Slater Type Orbitals STO and Gaussian Type Orbitals GTO . Slater type orbitals2 have the functional form shown in eq. 5.1 . Xc,n,im r, e, j NYm 0, j rn -1e-Zr 5.1 Here N is a normalization constant...
Force fields for metal coordination compounds
Coordination chemistry is an area that is especially plagued with the problems of assigning suitable functions for describing the individual energy terms and deriving good parameters.42 The bonding around metals is much more varied than for organic molecules, where there are just two, three or four bonds. Furthermore, for a given number of ligands, more than one geometrical arrangement is usually possible. A four-coordinated metal, for example, may either be tetrahedral or square planar, and a...
Concepts from Density Functional Theory
The success of FMO theory is not because the neglected terms in the second-order perturbation expansion eq. 15.1 are especially small an actual calculation will reveal that they completely swamp the HOMO-LUMO contribution. The deeper reason is that the shapes of the HOMO and LUMO resemble features in the total electron density, which determines the reactivity. There are also other quantities derived from density functional theory that directly relate to the properties and reactivity of...
Localized Orbitals
A Hartree-Fock wave function can be written as a single Slater determinant, composed of a set of orthonormal MOs eqs 9.19 and 3.20 . fl l f 2 1 L fN 1 For computational purposes, it is convenient to work with canonical MOs, i.e. those that make the matrix of Lagrange multipliers diagonal, and that are eigenfunctions of the Fock operator at convergence eq. 3.42 . This corresponds to a specific choice of a unitary orthogonal transformation of the occupied MOs. Once the SCF procedure has...
Examples
Tables 9.1 and 9.2 give some examples of atomic charges and bond orders calculated by various methods as a function of the basis set at the HF level of theory. It is evident that the Mulliken and Lowdin methods do not converge as the basis set is increased, and the values in general behave unpredictably. In particular, the presence of diffuse functions leads to absurd behaviours, as the aug-cc-pVXZ basis sets illustrate for CH4. Note also that for sufficiently large basis sets, the charge on...
Final Considerations
Should DFT methods be considered ab initio or semi-empirical If ab initio is taken to mean the absence of fitting parameters, LSDA methods are ab initio but gradient-corrected methods may or may not be. The LSDA exchange energy contains no parameters and the correlation functional is known accurately as a tabulated function of the density. The use of a parameterized interpolation formula in practical calculations does not represent fitting in order to improve the performance for atomic and...
Ts 1
Except for reactions with low barriers i.e. lt 40kJ mol at T 300K or at high temperatures, the quantity AE kT is large and the last series can be neglected. The tunnelling correction is then given completely in terms of the magnitude of the imaginary frequency. For small values of u the first term may be Taylor expanded to give eq. 14.26 . k Bell 1 u 2 14.26 The first-order term is known as the Wigner correction.31 It is possible to derive tunnelling corrections for functional forms of the...
Ccsdtq 1
We have so far been careful to use the wording formal scaling. As already discussed, HF is formally an M4 method but in practice the scaling may be reduced all the way down to M1. Similarly, MP2 is formally an M5 method. However, an MP2 calculation consists of three main parts the HF calculation, the AO to MO integral transformation, and the MP2 energy calculation. Only the second part has a formal scaling of M5, the others are formal M4 steps. In the large system limit, the transformation...
Intrinsic Reaction Coordinate Methods
The optimization methods described in Sections 12.2-12.4 concentrate on locating stationary points on an energy surface. The important points for discussing chemical reactions are minima, corresponding to reactant s and product s , and saddle points, corresponding to transition structures. Once a TS has been located, it should be verified that it indeed connects the desired minima. At the TS the vibrational normal coordinate associated with the imaginary frequency is the reaction coordinate...
Conrotatory jta
Figure 15.23 State correlation diagram for the dis- and conrotatory ring-closure of butadiene Figure 15.23 State correlation diagram for the dis- and conrotatory ring-closure of butadiene The same conclusion may again be reached by considering only the HOMO orbital. For the conrotatory path the orbital interaction leads directly to a bonding orbital, while the orbital phases for the disrotatory motion lead to an antibonding orbital. Figure 15.24 HOMO orbital for the ring-closure of butadiene...
Excited States
The development of HF and correlated methods in the previous chapters has focused on the electronic ground state. In some cases it is also of interest to consider electronically excited states. It is useful to distinguish between two cases, depending on whether the excited state has the same or a different symmetry than the lower state s . The different symmetry case is easy to handle, as the lowest energy state of a given symmetry may be handled completely analogously to the ground state. An...
Restricted and Unrestricted HartreeFock
So far there has not been any restriction on the MOs used to build the determinantal trial wave function. The Slater determinant has been written in terms of spin-orbitals, eq. 3.20 , being products of a spatial orbital and a spin function a or b . If there are no restrictions on the form of the spatial orbitals, the trial function is an Unrestricted Hartree-Fock UHF wave function.11 The term different orbitals for different spins DODS is also sometimes used. If the interest is in systems with...
Basis set effect at the HartreeFock level
Figure 11.1 shows the bond dissociation curves at the HF level with the STO-3G, 3-21G, 6-31G d,p , cc-pVDZ and cc-pVQZ basis sets. The total energy drops considerably upon going from the STO-3G to the 3-21G and again to the 6-31G d,p basis. This is primarily due to the improved description of the oxygen 1s-orbital. The two different types of DZP basis sets, 6-31G d,p and cc-pVDZ, give very similar results, and the improvement upon going to the cc-pVQZ basis is relatively minor. More important...
Qualitative Molecular Orbital Theory
Frontier molecular orbital theory is closely related to various schemes of qualitative orbital theory where interactions between fragment MOs are considered.14 Ligand field theory, as commonly used in systems involving coordination to metal atoms, can be considered as a special case where only the d-orbitals on the metal and selected orbitals of the ligands are considered. Two interacting orbitals will in general produce two new orbitals, having lower and higher energies than the...
Ahlrichs type basis sets
The group centred around R. Ahlrichs has designed basis sets of DZ, TZ and QZ quality for the elements up to Kr. The Split Valence Polarized SVP basis set is a 3s2p contraction of a 7s4p set of primitive functions contraction 5,1,1 and 3,1 , while the Triple Zeta Valence TZV basis set is a 5s3p contraction of an 11s6p set of primitive functions contraction 6,2,1,1,1 and 4,1,1 .23 More recently, the series has been extended by a Quadruple Zeta Valence QZV basis set, being a 7s4p contraction of a...
Translational degrees of freedom
The translational degrees of freedom can be exactly separated from the other 3N - 3 coordinates. The allowed quantum states for the translational energy are determined by placing the molecule in a box, i.e. the potential is zero inside the box but infinite outside. The only purpose of the box is to allow normalization of the translational wave function, i.e. the exact size is not important. The solutions to the Schrodinger equation for such a particle in a box are standing waves, cosine and...
Correlation consistent basis sets
The primary disadvantage of ANO basis sets is that a very large number of primitive GTOs are necessary for converging towards the basis set limit. Dunning and cowork-ers have proposed a somewhat smaller set of primitives that yields comparable results to the ANO basis sets.26 The correlation consistent cc the convention is to use lower case letters as the acronym, to distinguish it from coupled cluster CC basis sets are geared towards recovering the correlation energy of the valence electrons....