BEC Paper Q: Eigenvectors & Eigenvalues of Operator

  • Thread starter Coffee_
  • Start date
  • Tags
    Bec Paper
In summary, the author is explaining how to calculate the eigenstates of an operator for a single particle system. The author gives an example of how to do this, but says that it is not a unique procedure.
  • #1
Coffee_
259
2
http://bec.science.unitn.it/infm-bec/papers/preprints/becenc.pdf

It's in regards to expression (4)

I don't understand why the eigenvectors and eigenvalues of this operator happen to be the single particle eigenstates and the occupation numbers as they mention just below this expression.

I do understand that in general a hermitian operator can be written like this in terms of its eigenfunctions and eigenvalues - I just don't understand how to identify these eigenfunctions as the single particle eigenstates.

By single particle eigenstates I understand them being the eigenstates of the Hamiltonian for a single particle.
 
Physics news on Phys.org
  • #2
Think about how you construct the many-body state vector for a system of N identical particles. If the particles are all identical then you might be tempted to simply add their single-particle state vectors. But that won’t do because if the particles are identical then they must be genuinely indistinguishable. A simple addition of state vectors leaves the particles distinguishable by their generalised co-ordinates.

There are two ways you can construct the many-body state vector: by taking a linear combination of the single-particle state vectors which is either symmetric or anti-symmetric with respect to particle exchange.

Google “exchange interaction” to see many worked examples.

That particular exchange-symmetry, which you are forced to impose on the many-body wavefunction when you build it up from the individual particle’s wave functions means that when you construct the density matrix in step (4) of their calculation, the interference terms vanish and the single-particle eigenstates pop back out.
 
  • Like
Likes Greg Bernhardt and Nugatory
  • #3
So basically we are talking about the wave function of the system is:

##\Psi(r1,r2,...) = \Psi_1 (r1) \Psi_2 (r2) ... + \Psi_1(r2) \Psi_2 (r1) ... + ...## all permutations

Where ##\Psi_i## are the single particle eigenstates.

By plugging this into expression (3) of the paper I should get (4)?
 
  • #4
Yes... but... if you’re thinking of trying it yourself, there’s a subtlety here that I’m sure you’re aware of, although I don’t suppose it would do any harm to mention it.

You need to remember that the density matrix is fundamentally a statistical beast. It’s conventionally defined to be the ensemble average of the product of the ket and bra state vectors. You won’t be able to make the off-diagonal elements vanish without averaging over random phase in the integrands.

The assumption behind this is that the system is in thermal equilibrium. If the system is thermal then the relative phase of the integrands is randomised by contact with the thermal reservoir. (Absolute phase is, of course, irrelevant and un-observable.) When you take the ensemble average you are effectively averaging an ensemble of random numbers which yields zero.

(As an aside, if the system is coherently excited, it is indeed possible to obtain nonvanishing off-diagonal elements in the density matrix. The relaxation of this non-equilibrium state is called dephasing, or decoherence.)
 
  • #5
Daz said:
You need to remember that the density matrix is fundamentally a statistical beast. It’s conventionally defined to be the ensemble average of the product of the ket and bra state vectors.
The one density matrix has here little to do with statistical mechanics. Specifically, it does not involve an ensemble average.
The one density matrix is a hermitian one particle operator and as such there exists a diagonal representation in terms of its eigenfunctions, which are usually called natural orbitals.
Specifically, the eigenvalues and vectors fulfill
## n_i \phi_i(r)=\int dr' n^{(1)}(r',r) \phi_i(r') ##.
 
  • #6
Hi DrDu - thanks for your input. Clearly you know more about this than I do. I can see how your expression there works if the system is in a pure state but I’m struggling to see how it would work in a real system at finite temperature. I’m sure you’re right but any light you can shed on this would be greatly appreciated. The way I see it is that if you have a system consisting of a large number of interacting particles then surely the only way to describe that is statistically?
 
  • #7
I suppose that for a statistical ensemble, the one density matrix is just a partial trace over the full density matrix, i.e. tracing over all but one degree of freedom: ##n^{(1)}=\int dx_2\ldots dx_N \langle x_1 x_2\ldots x_N| \rho |x_N\ldots x_2 x'_1\rangle ##.
 
  • Like
Likes Daz
  • #8
OK Thanks for that. I’m still struggling, to be honest, but I think I need to sleep on this and maybe I’ll have more questions later.

To the OP (Coffee_): I think you might be well advised to follow DrDu’s guidance, rather than mine.
 
  • #9
Coffee_ said:
http://bec.science.unitn.it/infm-bec/papers/preprints/becenc.pdf
I just don't understand how to identify these eigenfunctions as the single particle eigenstates.
By single particle eigenstates I understand them being the eigenstates of the Hamiltonian for a single particle.
That is in deed a subtle question.
The problem is that there is no unique one particle hamiltonian for a system of interacting electrons.
If the system is non-interacting, then you (yes, you) can show that the electronic wavefunction is a Slater determinant formed from the eigenstates ##\phi_i## of the one-particle hamiltonian h. The natural orbitals are identical to the ##\phi_i## and the ##n_i = 1## for the occupied orbitals and 0 for the unoccupied ones.

Edit: Ups, this holds for a system of Fermions. For a system of Bosons, as we are interested here, instead of a Slater determinant you have a permanent and the occupation numbers are not restricted to be 1. Rather, for the ground state, n_0=N.
 
  • #10
DrDu thanks for your input.

I have tried what you mention, and even for a 2 particle state I seem to fail. Let me show you my attempt:

http://imgur.com/LiUDLnA
 
  • #11
DrDu said:
...

Hey, sorry if I'm persisting but I'm quite eager to lift this confusion/misunderstanding. Could you at least tell me vaguely where in my attempt I go wrong? I have overlooked it a few times and don't see the mistake.
 
  • #12
Sorry, I couldn't read what you posted on my smartphone. If psi_1 and psi_2 are orthogonal, then the integrals vanish and it seems you are done.
 
  • #13
I feel so stupid... That's what you get from position representation, should have seen it's just <f1 | f2 >

Oh really how silly, but still thanks a lot.

It would seem that this occupation level ##n_l## they speak of is just ##\sqrt(2)## is I consider the initial function to have another ##\frac{1}{\sqrt{2}}##
 
  • #14
Coffee_ said:
I feel so stupid... That's what you get from position representation, should have seen it's just <f1 | f2 >

Oh really how silly, but still thanks a lot.

It would seem that this occupation level ##n_l## they speak of is just ##\sqrt(2)## is I consider the initial function to have another ##\frac{1}{\sqrt{2}}##
Sounds good. You should also look at the more trivial case where the ground state wavefunction is a simple product of identical functions ##\phi_1(r_1)\ldots \phi_1(r_i)\ldots\phi_1(r_N)##.
 
  • #15
Wait something must be wrong here, consider eq(4) from the paper. Now set r=r' and integrate over the volume. The left term is 1 and the right term is a sum over ##n_l##. This doesn't result in the normalization condition they mention.
 

1. What is an eigenvector and eigenvalue?

An eigenvector is a vector that does not change its direction when multiplied by a specific linear transformation, while an eigenvalue is the scalar value that represents the amount by which the eigenvector is scaled.

2. What is the significance of eigenvectors and eigenvalues in mathematics and science?

Eigenvectors and eigenvalues are essential in solving many mathematical and scientific problems, particularly in linear algebra, quantum mechanics, and data analysis. They provide valuable information about the behavior and properties of linear transformations and operators.

3. How are eigenvectors and eigenvalues calculated?

Eigenvectors and eigenvalues can be determined by solving the characteristic equation (det(A-λI)=0) of a square matrix, where A is the matrix and λ is the eigenvalue. The eigenvectors are then found by solving the system of linear equations (A-λI)x=0.

4. What is the relationship between eigenvectors and eigenvalues?

The eigenvectors and eigenvalues of a linear operator are closely related. The eigenvalues represent the scaling factor of the eigenvectors, and each eigenvalue corresponds to a unique eigenvector. The eigenvectors form a basis for the vector space in which the linear operator acts.

5. How are eigenvectors and eigenvalues used in quantum mechanics?

In quantum mechanics, eigenvectors and eigenvalues are used to describe the state of a quantum system and the possible outcomes of measurements. The eigenvectors of an operator represent the states of the system, while the eigenvalues represent the possible results of measuring that operator on the system.

Similar threads

  • Atomic and Condensed Matter
Replies
1
Views
3K
  • Atomic and Condensed Matter
Replies
0
Views
247
  • Quantum Physics
Replies
2
Views
907
  • Quantum Physics
Replies
10
Views
2K
  • Atomic and Condensed Matter
Replies
4
Views
1K
Replies
18
Views
2K
  • Atomic and Condensed Matter
Replies
5
Views
1K
  • Advanced Physics Homework Help
Replies
1
Views
1K
Replies
67
Views
4K
Replies
9
Views
1K
Back
Top