1 / 7

Jack Simons , Henry Eyring Scientist and Professor Chemistry Department University of Utah

Electronic Structure Theory Session 8. Jack Simons , Henry Eyring Scientist and Professor Chemistry Department University of Utah. Coupled-Cluster Theory (CC): one expresses the wave function as  = exp(T)  ,

kaylee
Download Presentation

Jack Simons , Henry Eyring Scientist and Professor Chemistry Department University of Utah

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Electronic Structure Theory Session 8 Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

  2. Coupled-Cluster Theory (CC): one expresses the wave function as  = exp(T) , where is a single CSF (usually a single determinant) used in the SCF process to generate a set of spin-orbitals. The operator T is given in terms of operators that generate spin-orbital excitations: T = i,m tim m+ i + i,j,m,n ti,jm,n m+ n+ j i + ..., Here m+ i denotes creation of an electron in spin-orbitalmand removal of an electron from spin-orbitalito generate a single excitation. The operation m+ n+ j i represents a double excitation fromij tom n.

  3. When including in T only double excitations { m+ n+ j i}, the CC wave function exp(T) contains contributions from double, quadruple, sextuple, etc. excited determinants: exp(T)  = {1 + m,n,Iij tm,n,i,j m+ n+ j i + 1/2 (m,n,Iij tm,n,i,j m+ n+ j i) (m,n,Iij tm,n,i,j m+ n+ j i) + 1/6 (m,n,Iij tm,n,i,j m+ n+ j i) (m,n,Iij tm,n,i,j m+ n+ j i) (m,n,Iij tm,n,i,j m+ n+ j i) + …}. But note that the amplitudes of the higher excitations are given as products of amplitudes of lower excitations (unlinked).

  4. You may recall from your studies in statistical mechanics the so-called Mayer-Mayer cluster expansion in which the potential energy’s contribution to the partition function Q = exp(-V/kT)dr1 dr2 ... drN is approximated. Q = exp(-V/kT)dr1 dr2 ... drN = exp(-J,KV(|rJ-rK|) /kT)dr1 dr2 ... drN = J,Kexp(-V(|rJ-rK|) /kT)dr1 dr2 ... drN = J,K(1+ {exp(-V(|rJ-rK|) /kT)-1})dr1 dr2 ... drN = (1+J,K(exp(-V(|rJ-rK|) /kT)}-1) dr1 dr2 ... drN + J,K I,L (exp(-V(|rJ-rK|) /kT)}-1) (exp(-V(|rI-rL|) /kT)}-1) dr1 dr2 ... drN +... The unlinked terms are more in number and are important at lower densities. VN VN-1N(N-1)/2∫ (exp(-V(r) /kT)}-1) dr VN-2 N(N-1)(N-2)(N-3)/4∫(exp(-V(r) /kT)}-1)dr ∫(exp(-V(r’) /kT)}-1)dr’ VN-3 N(N-1)(N-2)/2∫(exp(-V(r1,2) /kT)}-1) (exp(-V(r1,3) /kT)}-1)dr1dr2dr3

  5. To obtain the equations of CC theory, one writes: Hexp(T) exp(T) then exp(-T) Hexp(T) then uses the Baker-Campbell-Hausdorf expansion: exp(-T) H exp(T) = H -[T,H] + 1/2 [[T,T,H]] - 1/6 [[[T,T,T,T,H]]] +… The equations one must solve for the t amplitudes are quartic: < im | H + [H,T] + 1/2 [[H,T],T] + 1/6 [[[H,T],T],T] + 1/24 [[[[H,T],T],T],T] |  > = 0; < i,jm,n |H + [H,T] + 1/2 [[H,T],T] + 1/6 [[[H,T],T],T] + 1/24 [[[[H,T],T],T],T] |> =0; < i,j,km,n,p|H + [H,T] + 1/2[[H,T],T] + 1/6 [[[H,T],T],T] + 1/24 [[[[H,T],T],T],T] |> = 0, The amplitudes of the double excitations that arise in the lowest approximation are identical to those of MP2 ti,jm,n = - < i,j | e2/r1,2 | m,n >'/ [ m-i +n -j ].

  6. Summary: Basis sets should be used that are flexible in the valence region to allow for the different radial extents of the neutral and anion’s orbitals, include polarization functions to allow for good treatment of electron correlations, and, include extra diffuse functions if very weak electron binding is anticipated. For high precision, it is useful to carry out basis set extrapolations using results calculated with a range of basis sets (e.g., VDZ, VTZ, VQZ). Electron correlation should be included because correlation energies are significant (e.g., 0.5 eV per electron pair). Correlation allows the electrons to avoid one another by forming polarized orbital pairs. There are many ways to handle electron correlation (e.g., CI, MPn, CC, DFT, MCSCF). Single determinant zeroth order wave functions may not be adequate if the spin and space symmetry adapted wave function requires more than one determinant. Open-shell singlet wave functions are the most common examples for which a single determinant can not be employed. In such cases, methods that assume dominance of a single determinant should be avoided. The computational cost involved in various electronic structure calculations scales in a highly non-linear fashion with the size of the AO basis, so careful basis set choices must be made.

  7. <ab|g|cd> M4/8 to be calculated and stored- linear scaling is being pursued to reduce this. F, Ci, = i S, Ci, operations to find all M i and Ci, <ij|g|kl> M4 of them each needing M5 steps. J HI,J CJ = E CI CI requires NC2 operations to get one E MP2- scales as M5

More Related