Matrix elements of one-body operator

O|n'\rangle a_n^\dagger a_{n'} a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle = \sum_{nn'}\delta_{n_1n'}\langle n|O|n'\rangle a_{n_3}^\dagger |0\rangle = \langle 0 | a_{n_3}^\dagger \langle n_1 | O | n_1 \rangle |0\rangle = \langle n_1 | O | n_1 \rangle .Note that the second
  • #1
Bill Foster
338
0

Homework Statement



Consider two fermion states:

[tex]|n_1 n_2\rangle = a_{n_1}^\dagger a_{n_2}^\dagger|0\rangle[/tex]

and

[tex]|n_1 n_3\rangle = a_{n_1}^\dagger a_{n_3}^\dagger|0\rangle[/tex]

where [tex]a_n^\dagger[/tex] denotes the fermion creation operator in the single-particle state [tex]u_n\left(\vec{r},s\right)[/tex]. Evaluate the matrix elements of the one-body operator [tex]\hat{O}=\sum_{nn'}\langle n|O|n'\rangle a_n^\dagger a_{n'}^\dagger[/tex] between these two states.

The Attempt at a Solution



Not sure if I should start out like this:

[tex]\langle 0| a_{n_1}^\dagger a_{n_2}^\dagger \sum_{nn'} \langle n|O|n'\rangle a_n^\dagger a_{n'}^\dagger a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]

of if I should start out like this:

[tex]\langle n_1 n_2 | \sum_{nn'} \langle n|O|n'\rangle a_n^\dagger a_{n'}^\dagger|n_1 n_3\rangle[/tex]
 
Physics news on Phys.org
  • #2
You'll need to use the anticommutation relations to compute the matrix element, so you might as well start with the form with all of the operators

[tex]\langle 0| a_{n_1}^\dagger a_{n_2}^\dagger \sum_{nn'} \langle n|O|n'\rangle a_n^\dagger a_{n'}^\dagger a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]

However, there is an error with this formula, since taking the conjugate of [tex]|n_1n_2\rangle[/tex] leads to

[tex]\langle n_1 n_2 | = \langle 0 | a_{n_2} a_{n_1}. [/tex]

I'd also question whether or not there's a typo in that one-body operator, since it doesn't conserve the fermion number and it's matrix element between even particle states vanishes.

[tex]\sum_{nn'} \langle n|O|n'\rangle a_n^\dagger a_{n'}[/tex]

would make more sense.
 
  • #3
fzero said:
You'll need to use the anticommutation relations to compute the matrix element, so you might as well start with the form with all of the operators
However, there is an error with this formula, since taking the conjugate of [tex]|n_1n_2\rangle[/tex] leads to

[tex]\langle n_1 n_2 | = \langle 0 | a_{n_2} a_{n_1}. [/tex]

I'd also question whether or not there's a typo in that one-body operator, since it doesn't conserve the fermion number and it's matrix element between even particle states vanishes.

[tex]\sum_{nn'} \langle n|O|n'\rangle a_n^\dagger a_{n'}[/tex]

would make more sense.

You are correct. The typo was on my part.

So I start with this:

[tex]
\langle 0| a_{n_1} a_{n_2} \sum_{nn'} \langle n|O|n'\rangle a_n^\dagger a_{n'} a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle
[/tex]

And then start applying the anti-commutation relations. I'll give that a shot.

Thanks.
 
  • #4
[tex]\{a_{n'},a_{n}^\dagger\}=a_{n'}a_{n}^\dagger + a_{n}^\dagger a_{n'}=\delta_{nn'}[/tex]

[tex]a_{n}^\dagger a_{n'}=\delta_{nn'}-a_{n'}a_{n}^\dagger[/tex]

Plugging that in...

[tex]\langle 0| a_{n_1} a_{n_2} \sum_{nn'} \langle n|O|n'\rangle \left(\delta_{nn'}-a_{n'}a_{n}^\dagger\right) a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]
[tex]=\langle 0| a_{n_1} a_{n_2} \sum_{nn'} \langle n|O|n'\rangle \delta_{nn'} a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle-\langle 0| a_{n_1} a_{n_2} \sum_{nn'} \langle n|O|n'\rangle a_{n'}a_{n}^\dagger a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]
[tex]=\langle 0| a_{n_1} a_{n_2} \langle n|O|n\rangle a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle-\langle 0| a_{n_1} a_{n_2} \sum_{nn'} \langle n|O|n'\rangle a_{n'}a_{n}^\dagger a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]

So my question now is this:

Is the following true?

[tex]\langle n|O|n\rangle = O[/tex]

If so, then I have this:

[tex]=\langle 0| a_{n_1} a_{n_2} O a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle-\langle 0| a_{n_1} a_{n_2} \sum_{nn'} \langle n|O|n'\rangle a_{n'}a_{n}^\dagger a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]
[tex]=\langle n_1 n_2| O |n_1 n_3\rangle-\langle 0| a_{n_1} a_{n_2} \sum_{nn'} \langle n|O|n'\rangle a_{n'}a_{n}^\dagger a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]

Now I'm not quite sure how to deal with that second term.
 
  • #5
Bill Foster said:
Is the following true?

[tex]\langle n|O|n\rangle = O[/tex]

No, but [tex]\langle n|O|n'\rangle [/tex] is a c-number and can be brought out of the expectation value that you're computing.

As for how to do these computations in general, note that

[tex]a_i |0\rangle =0.[/tex]

Therefore, we should use the anticommutation relations to move the [tex]a^\dagger[/tex]s to the left.
 
  • #6
fzero said:
No, but [tex]\langle n|O|n'\rangle [/tex] is a c-number and can be brought out of the expectation value that you're computing.

I don't know what you mean by "c-number".
 
  • #7
Bill Foster said:
I don't know what you mean by "c-number".

A c-number is an ordinary scalar quantity. It doesn't depend on the state vectors or other operators in an expectation value, so

[tex]\langle \alpha | c A | \beta \rangle = c \langle \alpha | A | \beta \rangle .[/tex]
 
  • #8
Which is true:

[tex]\langle 0|a_{n_3}^\dagger \langle n_2|O|n_1\rangle a_{n_1}|0\rangle = 0[/tex]

or


[tex]\langle 0|a_{n_3}^\dagger \langle n_2|O|n_1\rangle a_{n_1}|0\rangle = \langle n_2|O|n_1\rangle[/tex]
 
  • #9
Bill Foster said:
Which is true:

[tex]\langle 0|a_{n_3}^\dagger \langle n_2|O|n_1\rangle a_{n_1}|0\rangle = 0[/tex]

or


[tex]\langle 0|a_{n_3}^\dagger \langle n_2|O|n_1\rangle a_{n_1}|0\rangle = \langle n_2|O|n_1\rangle[/tex]

Those are both zero because the a annihilates the vacuum. On the other hand,


[tex]\langle 0|a_{n_3} \langle n_2|O|n_1\rangle a_{n_1}^\dagger|0\rangle = \delta_{n_1 n_3}\langle n_2|O|n_1\rangle .[/tex]
 
  • #10
So once I get one creation operator all the way to the left (or one annihilation operator all the way to the right), then that entire term goes to zero, correct?

For example, suppose I have:

[tex]\langle 0|a_{n_1} \sum_{nn'}\langle n|O|n' \rangle a_n^\dagger a_{n_2} a_{n'} a_{n_1}^\dagger a_{n_3}^\dagger |0\rangle[/tex]

Then all I need to do is swap [tex]a_{n_1}[/tex] with [tex]a_n^\dagger[/tex] and then that term will go to zero? (Except, of course), the additional term resulting from the delta function).
 
  • #11
Yes, you just have to get one annihilation operator to the right (or vice versa). The remaining delta terms contain all of the information about the original matrix element.
 
  • #12
Consider the product:

[tex]
\left\langle p_{1}, \ldots, p_{m} | q_{1}, \ldots q_{n} \right\rangle \equiv\left\langle 0 | a_{p_{1}} \, \ldots a_{p_{m}} \, a^{\dagger}_{q_{n}} \ldots a^{\dagger}_{q_{1}} | 0 \right\rangle
[/tex]

Let us start moving [itex]a_{p_{m}}[/itex] to the right. With each creation operator it gives a contraction:

[tex]
a_{p_{m}} \, a^{\dagger}_{q_{k}} = \delta_{p_{m}, q_{k}} + \zeta \, a^{\dagger}_{q_{k}} \, a_{p_{m}}
[/tex]

and one factor of [itex]\zeta[/itex] when the flipping occurs. When it reaches the furthest right, it annihilates the vacuum giving zero [itex]a_{p_{m}} \, | 0 \rangle = 0[/itex]. Finally, we are left with the following sum:

[tex]
\sum_{k_{m} = 1}^{n}{\zeta^{k_{m}} \, \delta_{p_{m}, q_{k}} \, \langle 0 | a_{p_{1}} \, \ldots a_{p_{m - 1}} \,\ldots a^{\dagger}_{q_{k_{m} + 1}} \, a^{\dagger}_{q_{k_{m} - 1}} \ldots | 0 \rangle}
[/tex]

The brackets in the sum differ from the original one by a [itex]a_{p_{m}}[/itex] and a [itex]a_{q_{k_{m}}}[/itex] and each one has a factor of [itex]\zeta^{k_{m}}[/itex], where [itex]\zeta = +1[/itex] for bosons and [itex]\zeta = -1[/itex] for fermions.

We may repeat this procedure for [itex]a_{p_{m - 1}}[/itex] for each of the terms in the sum. Then, we will get a double sum with one pair of creation-annihilation operators less and some power of [itex]\zeta[/itex]. Continuing this procedure, we will be left with either a bracket of the form:

[tex]
\langle 0 | a_{p_{1}} | 0 \rangle = 0, \; m > n
[/tex]

or:

[tex]
\langle 0 | a^{\dagger}_{q_{j_{n - m}}} \, a^{\dagger}_{q_{j_{n - m - 1}}} \, a^{\dagger}_{q_{j_{1}}} | 0 \rangle = 0, \; n > m
[/tex]

So, we can only get a non-zero matrix if [itex]m = n[/itex], i.e. we have the same number of annihilation operators as creation operators. Furthermore, due to the Kronecker deltas, we must have the single-particle states [itex]p_{1}, \ldots , p_{n}[/itex] be some permutation of the states [itex]q_{1}, \ldots q_{n}[/itex]. Finally, the power of [itex]\zeta[/itex] is equal to the number of transopsitions necessary to bring the sequence [itex]\{q_{1}, \ldots , q_{n}\}[/itex] to be identical to the sequence [itex]\{p_{1}, \ldots, p_{n}[/itex] (notice that [itex]\zeta^{2} = 1[/itex] for any kind of particle).

You could have used this rule if you wrote the one-body operator as:

[tex]
O = \sum_{m, m'}{(m | O | m') a^{\dagger}_{m} \, a_{m'}} = \zeta \, \sum_{m, m'}{(m' | O | m) a^{\dagger}_{m} \, a_{m'}} - \zeta \, \sum_{m}{(m | O |m)}
[/tex]

EDIT:
I noticed that you wrote the one-body operator with two creation operators. If this is so, then the matrix element is zero.
 
Last edited:
  • #13
Here's my final answer:

[tex]\langle 0|\langle n_2|O|n_1\rangle \delta_{n_1 n_3}|0\rangle - \langle 0|\langle n_1|O|n_1\rangle \delta_{n_2 n_3}|0\rangle - \langle 0|\langle n_1|O|n_3\rangle \delta_{n_1 n_1}|0\rangle + \langle 0|\langle n_1|O|n_3\rangle \delta_{n_2 n_1}|0\rangle[/tex]
 
  • #14
Dickfore said:
I noticed that you wrote the one-body operator with two creation operators. If this is so, then the matrix element is zero.

That was a typo.

I don't suppose you all could help me with this problem: https://www.physicsforums.com/showthread.php?t=431118

It's another matrix element problem.
 
  • #15
Bill Foster said:
Here's my final answer:

[tex]\langle 0|\langle n_2|O|n_1\rangle \delta_{n_1 n_3}|0\rangle - \langle 0|\langle n_1|O|n_1\rangle \delta_{n_2 n_3}|0\rangle - \langle 0|\langle n_1|O|n_3\rangle \delta_{n_1 n_1}|0\rangle + \langle 0|\langle n_1|O|n_3\rangle \delta_{n_2 n_1}|0\rangle[/tex]

The 3rd term should have [tex]\langle n_2|O|n_3\rangle[/tex]. I'd simplify this by pulling the [tex]O[/tex] matrix elements out and using [tex]\langle 0 | 0 \rangle=1.[/tex]
 
  • #16
Bill Foster said:
That was a typo.

I don't suppose you all could help me with this problem: https://www.physicsforums.com/showthread.php?t=431118

It's another matrix element problem.

I don't understand what that problem is trying to get at. The Coulomb potential commutes with the spin operators, so it doesn't make any sense for the matrix element to depend on spin. To see spin effects in a two electron system, you have to consider the magnetic fields of the electrons, which results in the spin-orbit coupling. The Coulomb potential is electrostatic.

There's a discussion of the expectation value of the Coulomb Hamiltonian in Baym's QM book, but I can't view the entire argument on Google Docs. He certainly doesn't have any spin dependence in his result.

Edit: I figured it out. There's an energy splitting due to Fermi statistics. I'll post in the other thread.
 
Last edited:

FAQ: Matrix elements of one-body operator

What are matrix elements of one-body operator?

Matrix elements of one-body operator are the elements of a matrix that represent the action of a one-body operator on a particular state or basis vector. They provide information about the expectation value of a one-body operator in a quantum mechanical system.

How are matrix elements of one-body operator calculated?

Matrix elements of one-body operator are calculated by taking the inner product of the operator with the basis vectors of the system. This involves taking the complex conjugate of the operator and then multiplying it by the basis vector and integrating over all space.

What is the significance of matrix elements of one-body operator in quantum mechanics?

Matrix elements of one-body operator are important in quantum mechanics because they allow us to calculate the expectation value of a one-body operator in a quantum mechanical system. This provides information about the physical properties of the system and helps us understand its behavior.

Can matrix elements of one-body operator be used to describe the time evolution of a quantum system?

Yes, matrix elements of one-body operator can be used to describe the time evolution of a quantum system. By taking the matrix elements of the Hamiltonian operator, we can calculate the time evolution of the system and determine how it changes over time.

How do matrix elements of one-body operator change when the basis set is changed?

Matrix elements of one-body operator can change when the basis set is changed. This is because the basis vectors are different and therefore the inner product will be different. However, the overall physical properties of the system should not change, as long as the new basis set is complete and orthonormal.

Similar threads

Back
Top