Quantum mechanics exercise

eoghan

1. The problem statement, all variables and given/known data
Let there be 3 particles with mass m moving in the 1D potential:
[tex]\frac{k}{2}[(x_1-x_2)^2 + (x_2-x_3)^2 + (x_1-x_3)^2][/tex]
where [tex]x_i[/tex] is the coordinate of the particle i.

1)Show that with the following coordinat change the Schroedinger equation is easy to solve:
[tex]y_1=x_1-x_2[/tex]
[tex]y_2=\frac{1}{2}(x_1+x_2)-x_3[/tex]
[tex]y_3=\frac{1}{3}(x_1+x_2+x_3)[/tex]


2) Find the eigenstates and the energies of the equation you got in point 1)

2. Relevant equations



3. The attempt at a solution
[tex]x_1-x_2=y_1[/tex]
[tex]x_2-x_3=y_2-\frac{1}{2}y_1[/tex]
[tex]x_1-x_3=y_2+\frac{1}{2}y_1[/tex]

[tex]V=\frac{k}{2}\left[\frac{3}{2}y_1^2+2y_2^2\right][/tex]
[tex]H=\frac{P_1^2}{2m}+\frac{P_2^2}{2m}+\frac{P_3^2}{2m}+V[/tex]

So I have 2 independent harmonic oscillators with angular frequencies [tex]\sqrt{\frac{3k}{2m}}[/tex] and [tex]\sqrt{\frac{2k}{m}}[/tex]
and a free particle whose eigenfunction is [tex]exp\left[\frac{i}{\hbar}\vec P\vec r\right][/tex]
So the eigenstates are the tensor product of the eigenstates of two harmonic oscillators and an exponential.
The energies are [tex](a+\frac{1}{2})\hbar w_1+(b+\frac{1}{2})\hbar w_2 + E[/tex]
where w1 and w2 are the two frequencies of the two harmonic oscillators and E is the energy of the free particle.

vela

You need to convert the momenta to the new coordinates as well.

eoghan

ok, if I convert the momenta I get:
[tex]H=\frac{2P_1^2}{2m}+\frac{3}{4m}P_2^2+\frac{5}{18m}P_3^2+V[/tex]

Where [tex]P_i[/tex] now refers to the momenta in the new basis.

Now the new Hamiltonian can be written as [tex]H=H_1+H_2+H_3[/tex]
So I have the hamiltonian of 3 independent particles:
[tex]H_1=2\left[\frac{P_1^2}{2m}+\frac{k}{2}\left(\frac{3}{4}y_1^2\right)\right][/tex]
[tex]H_2=\frac{3}{2}\left[\frac{P_2^2}{2m}+\frac{k}{2}\left(\frac{9}{4}y_2^2\right)\right][/tex]
[tex]H_3=\frac{5}{9}\left(\frac{P_3^2}{2m}\right)[/tex]

H1 and H2 are two independent harmonic oscillators with their own frequency, and H3 is an independent particle whose energy is multiplied by 5/9
Is this right?

vela

I got the same potential in terms of the y's, but for the kinetic term, I found

[tex]\frac{1}{2m}\left(\frac{1}{2}p_1^2+\frac{2}{3}p_2^2+3 p_3^2\right)[/tex]

Besides the constant factors, though, you're correct in that you get what looks like two oscillators and a free particle. The coordinate y₃ is the center of mass of the three particles, so the "free particle" is actually the system taken as a whole.

eoghan

Uhm... I calculate the momenta in the new coordinates again and I get:
[tex]\frac{1}{2m}\left(2p_1^2+\frac{3}{2}p_2^2+\frac{1}{3}p_3^2\right)[/tex]
Anyway... in the end, the energy is the sum of 3 energies:
E1: energy of a harmonic oscillator with mass m/2 (or 2m according to your result) and frequency [tex]\left(\omega=\sqrt{\frac{m}{k}}\right)[/tex] w=3 (or 3/4 according to your results);
E2: the energy of the second harmonic oscillator
E3: the energy of a free particle of mass 3m, i.e. E3=[tex]\frac{P_3^2}{3m}[/tex]

vela

Your answer is probably right. A mass of 3m for the whole system makes more sense than m/3.

EDIT: I redid the calculation and now get the same result you did.

eoghan

And if the 3 particles are bosons what is the ground state energy?
I guess it is:
3*(E1+E2+E3)
and the wavefunction is:
[tex]\frac{\psi_1(x_1)\psi_2(x_2)\psi_3(x_3)+\psi_1(x_1)\psi_2(x_3)\psi_3(x_ 2)+\psi_1(x_2)\psi_2(x_1)\psi_3(x_3)+\psi_1(x_2)\psi_2(x_3)\psi_3(x_1)+ \psi_1(x_3)\psi_2(x_2)\psi_3(x_1)+\psi_1(x_3)\psi_2(x_1)\psi_3(x_2)}{\s qrt{6}}[/tex]

But then, wouldn't the energy be 6*(E1+E2+E3)?

Where [tex]\psi_1[/tex] is the wavefunction of the first harmonic oscillator,
[tex]\psi_1[/tex] is the wavefunction of the second HO
and [tex]\psi_1[/tex] is the wavefunction of the free particle

eoghan

Ah.. no...the energy is still 3*(E1+E2+E3) because I have to divide by 6

vela

What are E1, E2, and E3? The argument of the wavefunctions should be the y's, not the x's. The Hamiltonian isn't symmetric with respect to exchanging the y's, so if [itex]\psi_1(y_1)\psi_2(y_2)\psi_3(y3)[/itex] is a solution, [itex]\psi_1(y_2)\psi_2(y_1)\psi_3(y3)[/itex] won't generally be one.

eoghan

E1: energy of a harmonic oscillator with mass m/2 (or 2m according to your result) and frequency [tex]\left(\omega=\sqrt{\frac{m}{k}}\right)[/tex] w=3;
E2: the energy of the second harmonic oscillator
E3: the energy of a free particle of mass 3m.

So the wave function of the ground state is simply [itex]\psi_1(y_1)\psi_2(y_2)\psi_3(y_3)[/itex]?
Where psi1 is the wavefunction of a HO with frequency w=3, psi2 is the wavefunction of the second HO and psi3 is the wavefunction of a free particle

vela

So why would the energy be 3*(E1+E2+E3)? Wouldn't it be just E1+E2+E3?

eoghan

Ah, yes, the energy is E1+E2+E3.
So, I don't care if the particles are bosons or fermions? In both cases the wave-function is
[itex]
\psi_1(y_1)\psi_2(y_2)\psi_3(y_3)
[/itex]

vela

I'm not sure actually. I'll see if someone else can offer some help on this point.

thegtoman

plug and chug

diazona

Yeah, I think it doesn't matter. As far as I remember, the Pauli exclusion principle really comes from the antisymmetrization of the wavefunctions of identical fermionic particles, but in this case, the masses corresponding to y₁ and y₂ (and y₃) are different, so the particles are distinguishable and we don't have that issue. And in any case, these are bosons you're talking about, so there's definitely no exclusion principle to worry about. Just drop everything in its lowest energy state and add them up.

eoghan

Thanks to everybody for your help!

vela

I'm wondering about the ground state wavefunctions for the fermion and boson cases. Because of the symmetry of the original Hamiltonian, I'd expect the individual particle wavefunctions to be the same and the ground state would be symmetric and antisymmetric combinations depending on the type of particle. Since the ground state wavefunctions in terms of the the x's are different, I'd expect the ground state wavefunctions in terms of the y's to be different.