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 '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$

then$A(\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.