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.
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?
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_{3} is the center of mass of the three particles, so the "free particle" is actually the system taken as a whole.
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]
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.
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)}{\sqrt{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
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.
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
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]
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_{1} and y_{2} (and y_{3}) 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.
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.
You can swap x_{2} and x_{3} in the original coordinate transformation to get [tex]\begin{align*} y'_1 &= x_1-x_3 \\ y'_2 &= \frac{x_2+x_3}{2}-x_1 \\ y'_3 &= \frac{x_1+x_2+x_3}{3} \end{align*}[/tex] and find a solution. Other transformations derived from permuting the x_{i}'s will work as well, so the ground state solution would be [tex]\begin{align*} \psi(x_1,x_2,x_3) = \frac{1}{\sqrt{3}} [ & \psi_1(x_1-x_2)\psi_2\left(\frac{x_1+x_2}{2} - x_3\right)\psi_3\left(\frac{x_1+x_2+x_3}{3}\right) + \\ & \psi_1(x_2-x_3)\psi_2\left(\frac{x_2+x_3}{2} - x_1\right)\psi_3\left(\frac{x_1+x_2+x_3}{3}\right) + \\ & \psi_1(x_3-x_1)\psi_2\left(\frac{x_3+x_1}{2} - x_2\right)\psi_3\left(\frac{x_1+x_2+x_3}{3}\right) ] \end{align*}[/tex] where [itex]\psi_1[/itex] and [itex]\psi_2[/itex] are the harmonic oscillator ground state wavefunctions (which are even so sign changes in the arguments don't matter) and [itex]\psi_3[/itex] is the free particle wavefunction. The ground state is symmetric so the spin part of the state would have to take care of the overall symmetry or asymmetry of the state.
Ok, so I have to return back to the original coordinates x's and make the wave function symmetric with respect to the x's. In that case the spin part would have to be symmetric as well. Is this right?