Linker IT Software
menubar-top-links menubar-top-rechts
Home Help Search Login
Welcome, Guest. Please Login.
SQL*XL: Database to Excel bridge litLIB: Excel power functions pack ExcelLock: Locking and securing your valuable Excel spreadsheets encOffice: Protect your Excel file easy and safe encOffice: Protect your Excel file easy and safe
Pages: 1
CASSCF wavefunction for open shell singlet (Read 3934 times)
Gerrit-Jan Linker
YaBB Administrator

Posts: 75
CASSCF wavefunction for open shell singlet
17.12.10 at 11:50:32
CASSCF wavefunction for open shell singlet
An open shell singlet is a state in which two electrons with opposite spin occupy a different orbital. It is an excited state compared to the singlet ground state in which both electrons occupy the same orbital.
A ground state singlet wavefunction can be written as: Φ=|a↑a↓|, and the open shell singlet Φ=|a↑b↓|.
In this notation Φ is the wavefunction a and b are the orbitals, ↑ and ↓ denote the spin function (either up or down). The orbitals are actually spin orbitals. A spin orbital is the product of a and b as their spatial component and the up and down as the spin component.  
Wavefunctions are commonly built as Slater determinants so they fulfill the Pauli principle, it can be expanded as follows:
Φ=|a↑b↓|= a↑1b↓2 - b↓1a↑2,
where 1 and 2 denote the coordinates of electron 1 and 2.  
In CASSCF an open shell singlet can be constructed using 2 electrons in 2 orbitals: a CAS (2,2). In such a calculation the wavefunction can be built from 3 configuration state functions (CSF's) namely: Φ20, Φud, Φ02. In this notation the 0,u,d,2 denote the occupation of the basis function with respectively 0, 1, 1, or 2. u and d denote the spin coupling, either up or down.
To bring the two notations together, Φ20 can be written as Φ20=|a↑a↓|, Φud as Φud=|a↑b↓| or Φud=|b↑a↓| and Φ02 as Φ02=|b↑b↓|
A open shell singlet wavefunction can be built from these 3 CSF's as:
Ψ=c1Φud or Ψ=c2Φ20-c3Φ02 or written differently:
Ψ=c1|a↑b↓| or Ψ=c2|l↑l↑|-c3|r↓r↓|
Note that different sets of orbitals are used for the two wavefunctions. We can write the open shell singlet as the following linear combination Ψ=|a↑b↓|=|a↑b↓|-|a↓b↑| and then perform a transformation of the orbitals to be able to do this. So, in the [ab] basis the open shell singlet can be written as Ψ=|a↑b↓|-|a↓b↑| and a transformation exist of the basis [ab] to a basis [lr] so that we can write the open shell singlet as Ψ=c2|l↑l↑|-c3|r↓r↓|
To show that the open shell singlet wavefunction can be built in this way we use the following transformation of the orbitals a and b to the orbitals l and r:
l=(a+b)√2, r=(a-b)√2
This leads to:
[1] a=(l+r)/√2, b=(l-r)/√2
We use the following wavefunction for the open shell singlet:
[2] Ψ=|a↑b↓|-|a↓b↑|
Expanding the determinants in [2] gives:
[3] Ψ= ( (a↑1b↓2 - b↓1a↑2) - (a↓1b↑2 - b↑1a↓2) ) = ( a↑1b↓2 - b↓1a↑2 - a↓1b↑2 + b↑1a↓2 )
Using [1] we rewrite [3] as:
Ψ = ( (l↑1+r↑1) (l↓2-r↓2)/2 - (l↓1-r↓1) (l↑2+r↑2)/2 - (l↓1+r↓1)(l↑2-r↑2)/2 + (l↑1-r↑1)(l↓2+r↓2)/2 )
= ( (l↑1+r↑1) (l↓2-r↓2) - (l↓1-r↓1) (l↑2+r↑2) - (l↓1+r↓1)(l↑2-r↑2) + (l↑1-r↑1)(l↓2+r↓2) )
expanding gives:
[4] Ψ = ( ( (l↑1l↓2+r↑1l↓2) - (l↑1r↓2+r↑1r↓2) ) - ( (l↓1l↑2-r↓1l↑2) + (l↓1r↑2-r↓1r↑2) ) - ( (l↓1l↑2+r↓1l↑2) - (l↓1r↑2+r↓1r↑2) ) + ( (l↑1l↓2-r↑1l↓2) + (l↑1r↓2-r↑1r↓2) ) )
= ( l↑1l↓2 + r↑1l↓2 - l↑1r↓2 - r↑1r↓2 - l↓1l↑2 + r↓1l↑2 - l↓1r↑2 + r↓1r↑2 - l↓1l↑2 - r↓1l↑2 + l↓1r↑2 + r↓1r↑2 + l↑1l↓2 - r↑1l↓2 + l↑1r↓2 - r↑1r↓2 )
= ( l↑1l↓2 + l↑1l↓2 - l↓1l↑2 - l↓1l↑2 - l↓1r↑2 + l↓1r↑2 + l↑1r↓2 - l↑1r↓2 + r↓1l↑2 - r↓1l↑2 - r↑1l↓2 + r↑1l↓2 - r↑1r↓2 - r↑1r↓2 + r↓1r↑2 + r↓1r↑2 )
= ( 2l↑1l↓2 - 2l↓1l↑2 - 2r↑1r↓2 + 2r↓1r↑2 )
= ( l↑1l↓2 - l↓1l↑2 - r↑1r↓2 + r↓1r↑2 )
= ( (l↑1l↓2 - l↓1l↑2 )- (r↑1r↓2 - r↓1r↑2 ))
= |l↑l↓| - |r↑r↓|
Back to top
« Last Edit: 17.12.10 at 14:28:04 by Gerrit-Jan Linker »  

Gerrit-Jan Linker
Linker IT Software
Email WWW Gerrit-Jan Linker   IP Logged
Pages: 1