Tight binding model (covalent bonding)

[amjad-sh]

Homework Statement

(linear combination of atomic orbitals):
Lets consider two atoms which are bond together with a covalent bond. Let's consider any sets of wavefunctions |n\rangle for n=1,2,...,N. Let's call orbital |1\rangle around nucleus 1 and orbital |2\rangle around nucleus 2 and so on.
for simplicity we assume this basis is orthonormal \langle n|m\rangle=\delta_{n,m}
Let us write a trial wavefunction for our ground state as |\psi\rangle=\sum_n\phi_n|n\rangle.
This is known as linear combination of atomic orbitals, LCAO,or tight binding.

We would like to find the lowest -energy wavefunction we can construct in this form.i.e, the best approximation to the actual ground state wavefunction. ( The more states we use in our basis, generally, the more accurate our results will be).

We claim that the ground state is given by the solution of the effective Schroedinger equation H\phi=E\phi where \phi is the vector of N coefficients \phi_{n} and H is the N by N matrix H_{n,m}=\langle n|H|m \rangle with H the full system we are considering.

we want to prove this by constructing the energy E=\frac{\langle\psi|H|\psi\rangle}<br /> {\langle\psi|\psi\rangle}

The question is : Show that minimizing this energy E w.r.t to each \phi_{n} gives the same eigenvalue equation above.(Caution:\phi_n is generally complex! )

Homework Equations

//////////////

The Attempt at a Solution

[/B]
after expanding
E=\frac{\langle\psi|H|\phi\rangle}<br /> {\langle\psi|\psi\rangle}
I tried to minimize the energy E w.r.t to \phi_{n} by letting the derivative of it w.r.t to\phi_{n} equals to zero, but this didn't lead me to the eigenvalue equation, if anybody can give me more hints? thanks!

 

[amjad-sh]

The Attempt at a Solution

E=\frac{\langle\psi|H|\psi\rangle}{\langle\psi|\psi\rangle}=

\frac{\phi^*_1\phi_1\langle1|H|1\rangle+...+\phi^*_1\phi_n\langle1|H|1\rangle+\phi^*_2\phi_1\langle 2|H|1\rangle+...+\phi^*_2\phi_n\langle 2|H| n\rangle+...+\phi^*_n\phi_1\langle n|H|1\rangle+...+\phi^*_n\phi_n\langle n|H|n\rangle}{ \left| \phi_1 \right|^2+\left| \phi_2 \right|^2+...+\left| \phi_n \right|^2 }

Let A=\left| \phi_1 \right|^2+\left| \phi_2 \right|^2+...+\left| \phi_n \right|^2

Let's minimize E by differentiating it w.r.t to \phi_1

dE=\frac{\partial E}{\partial \phi_1}d\phi_1+...+\frac{\partial E}{\partial \phi_n}d\phi_n
then\frac{dE}{d\phi_1}=\frac{\partial E}{\partial \phi_1}\frac{d\phi_1}{d\phi_1}+...+\frac{\partial E}{\partial \phi_n}\frac{d\phi_n}{d\phi_1}as all the \phi_i &#039;s are independent then \frac{dE}{d\phi_1}=\frac{\partial E}{\partial \phi_1}\frac{\partial E}{\partial \phi_1}=\frac{\phi_1^*A\langle1|H|1\rangle-\phi^*_1\phi_1\phi^*_1\langle1|H|1\rangle}{A^2}+\frac{\phi_2^*A\langle2|H|1\rangle-\phi^*_2\phi_1\phi^*_1\langle2|H|1\rangle}{A^2}+...+\frac{\phi_n^*A\langle n|H|1\rangle-\phi^*_n\phi_1\phi^*_1\langle n|H|1\rangle}{A^2}=0

thenA(\phi^*_1\langle 1|H|1\rangle+\phi^*_2\langle 2|H|1\rangle+...+\phi^*_n\langle n|H|1\rangle)-\left| \phi_1\right|^2(\phi^*_1\langle 1|H|1\rangle+\phi^*_2\langle2|H|1\rangle+...+\phi^*_n\langle n|H|1\rangle)=0

this means that \phi^*_1\langle1|H|1\rangle+...+\phi^*_n\langle n|H|1\rangle=0
lets complex conjugate we will get \phi_1\langle1|H|1\rangle+...+\phi_n\langle 1|H|n\rangle=0

moving in this with every variable \phi_i we will get H\phi=0 where H is the N by N Hamiltonian matrix in the |1> ,|2>,...|n> basis and \phi is the column vector containg \phi_1,\phi_2,...\phi_n

I reached this part and I don't know how to proceed.

 

[IanBerkman]

We have $|\psi\rangle = \sum_n \phi_n |n\rangle$ and define $H|n\rangle = E_n|n\rangle$ which gives
$$\langle \psi| H | \psi \rangle = \langle m | \sum_m \phi_m^* \sum_n E_n \phi_n| n \rangle\\
= \sum_{m,n}\phi_m^*\phi_n E_n \langle m| n \rangle\\
= \sum_{m,n}\phi_m^*\phi_n E_n\delta_{m,n}\\
= \sum_n E_n |\phi_n|^2$$
The denominator gives
$$ \langle \psi | \psi \rangle = \sum_n |\phi_n|^2$$
I think from here on it is a bit easier.

 

[amjad-sh]

We have $|\psi\rangle = \sum_n \phi_n |n\rangle$ and define $H|n\rangle = E_n|n\rangle$ which gives
$$\langle \psi| H | \psi \rangle = \langle m | \sum_m \phi_m^* \sum_n E_n \phi_n| n \rangle\\
= \sum_{m,n}\phi_m^*\phi_n E_n \langle m| n \rangle\\
= \sum_{m,n}\phi_m^*\phi_n E_n\delta_{m,n}\\
= \sum_n E_n |\phi_n|^2$$
The denominator gives
$$ \langle \psi | \psi \rangle = \sum_n |\phi_n|^2$$
I think from here on it is a bit easier.

But I think this approach will not work since the orbitals|1>,|2>,...|n> are not the eigenstates of the Hamiltonian H which is the Hamiltonian of the electron when the two atoms are combined.

 

[amjad-sh]

I found the solution. :)

 

[IanBerkman]

Great :). Do you mind to give the solution? I find this problem also quite interesting.

 

[amjad-sh]

Great :). Do you mind to give the solution? I find this problem also quite interesting.

It is problem 6.2. My approach was right but I differentiated the Energy E w.r.t to phi1 wrongly.