Two independent particles in an infinite square well

9 replies · 3K views · last reply
gfd43tg
Gold Member
Messages
949
Reaction score
48

Homework Statement


Screenshot of a homework problem 2 based on Griffiths Problem 5.6, with parts a) and b) asking for the expectation value of (x1 - x2)^2 for two particles in an infinite square well, and for spatial wavefunction, spin and expectation value for electrons.


Homework Equations


Printed text introducing integral identities for integers m and n: integral from 0 to pi of x squared sin squared nx, integral of x sin squared nx, and integral of x sin mx sin nx, each with its closed-form right-hand side.

Boxed equation: psi_n(x) = sqrt(2/a) sin(n pi x / a), labelled [2.28].


The Attempt at a Solution


a) For this part, I know for distinguishable particles, the expectation value of the square distance

$$\langle (x_{1}^{2} - x_{2}^{2}) \rangle = \langle x^{2} \rangle_{2} + \langle x^{2} \rangle_{3} - 2 \langle x \rangle_{2} \langle x \rangle_{3}$$

This just gets really ugly...
$$\langle x^{2} \rangle_{2} = \int \psi_{2}^{*}x^{2}\psi_{2} dx = \frac {2}{a} \int x^{2} sin^{2}\Big (\frac {2 \pi}{a} x \Big ) dx $$
Using the equation given in the problem (is that even right? It seems suspicious that it is independent of ##n##)
$$\langle x^{2} \rangle_{2} = \frac {2}{a} \Big [ \frac {\pi^{3}}{6} - \frac {\pi}{4(\frac {2 \pi}{a})^{2}} \Big ] $$
$$\langle x^{2} \rangle_{2} = \frac {2}{a} \Big [ \frac {\pi^{3}}{6} - \frac {a^{2}}{16 \pi} \Big ] $$
$$ = \frac {\pi^{3}}{3a} - \frac {a}{8 \pi} $$

Similarly for ##\langle x^{2} \rangle_{3}##
$$\langle x^{2} \rangle_{3} = \frac {2}{a} \Big [ \frac {\pi^{3}}{6} - \frac {\pi}{4(\frac {3 \pi}{a})^{2}} \Big ]$$
$$\langle x^{2} \rangle_{3} = \frac {2}{a} \Big [ \frac {\pi^{3}}{6} - \frac {a^{2}}{36 \pi} \Big ] $$
$$=\frac {\pi^{3}}{3a} - \frac {a}{18 \pi}$$

For ##\langle x \rangle_{2}##
$$\langle x \rangle_{2} = \int \psi_{2}^{*}x \psi_{2} dx$$
$$\langle x \rangle_{2} = \frac {2}{a} \int x sin^{2}(\frac {2 \pi}{a} x) dx $$
From the equation given in the problem statement
$$\langle x \rangle_{2} = \frac {2}{a} \frac {\pi^{2}}{4} = \frac {\pi^{2}}{2a} $$
The same goes for ##\langle x \rangle_{3}##
$$\langle x \rangle_{3} = \frac {2}{a} \int x sin^{2}(\frac {3 \pi}{a} x) dx$$
$$\langle x \rangle_{3} = \frac {2}{a} \frac {\pi^{2}}{4} = \frac {\pi^{2}}{2a}$$

$$\langle (x_{1}^{2} - x_{2}^{2}) \rangle = \langle x^{2} \rangle_{2} + \langle x^{2} \rangle_{3} - 2 \langle x \rangle_{2} \langle x \rangle_{3}$$
$$\langle (x_{1}^{2} - x_{2}^{2}) \rangle = \Big (\frac {\pi^{3}}{3a} - \frac {a}{8 \pi} \Big ) + \Big (\frac {\pi^{3}}{3a} - \frac {a}{18 \pi} \Big ) - 2 \Big (\frac {\pi^{2}}{2a} \Big ) \Big ( \frac {\pi^{2}}{2a} \Big ) $$
$$\langle (x_{1}^{2} - x_{2}^{2}) \rangle = \frac {2 \pi^{3}}{3a} - \frac {13a}{72 \pi} - \frac {\pi^{4}}{2a^{2}} $$
----------------------------------------------------------------------------------------------------------------------------------------
(b) For this problem, I am basically looking at this webpage
http://galileo.phys.virginia.edu/classes/252/symmetry/Symmetry.html

symmetric:
$$\psi_{23} (x_{1},x_{2}) = A sin \Big (\frac {2 \pi}{a} x_{1} \Big ) cos \Big (\frac {3 \pi}{a} x_{2} \Big ) $$
Not sure how to get A, but anyways
For symmetric
$$ \psi_{23}^{S} = \frac {1}{\sqrt{2}} \Big [ sin \Big (\frac {2 \pi}{a} x_{1} \Big ) cos \Big (\frac {3 \pi}{a} x_{2} \Big ) +sin \Big (\frac {3 \pi}{a} x_{1} \Big ) cos \Big (\frac {2 \pi}{a} x_{2} \Big ) \Big ]$$

For antisymmetric
$$ \psi_{23}^{A} = \frac {1}{\sqrt{2}} \Big [ sin \Big (\frac {2 \pi}{a} x_{1} \Big ) cos \Big (\frac {3 \pi}{a} x_{2} \Big ) - sin \Big (\frac {3 \pi}{a} x_{1} \Big ) cos \Big (\frac {2 \pi}{a} x_{2} \Big ) \Big ] $$

This is slightly different from my notes, which had two sine terms instead of a cosine term, and I am not sure if that is because this time it has electrons. I am just so confused...
Handwritten-style typed notes showing the infinite square well single-particle wavefunction psi_n(x) = sqrt(2/a) sin(n pi x / a) and energy E_n = n^2 pi^2 hbar^2 / (2ma^2), then two-particle product wavefunction and total energy formulas for the ground state.

Handwritten-style derivation of symmetric and antisymmetric two-particle wavefunctions for indistinguishable particles in an infinite square well, with normalization A = 1/√2 and lowest-energy cases n1 = n2 = 1 and n1 = 1, n2 = 2.
 
Last edited:
For part (a), note that your potential well extends from x = 0 to x = a. Thus, you will integrate from x = 0 to x = a. Note, though, that Griffiths' integrals extend from 0 to ##\pi##. So, a change of integration variable will be required if you want to use Griffiths' integrals.

For part (b), the website that you link to assumes that the well extends from x = -L/2 to x = L/2. That causes some of the wavefunctions to be cosine functions. But, in your case where the well extends from x = 0 to x = a, all of the wavefunctions will be sine functions.
 
The answer you got in part (a) can't be correct because the units don't work out.

You should review how to construct the symmetric and antisymmetric spatial wave functions. It's pretty straightforward how to do so. You shouldn't need to copy a result off of a webpage.
 
I just posted my new work since writing out all the latex and editing it until it is right can take over an hour. But for part (a), I got

$$\langle (x_{1} - x_{2})^{2} \rangle = \frac {a^{2}}{27} - \frac {5a^{2}}{72 \pi^{2}}$$

At least the units are correct for this answer
Handwritten derivation on lined paper of wavefunctions and expectation values <(x1 - x2)^2> for two particles in an infinite square well, with integrals and a boxed final result.


(b) for part i. and ii.
Handwritten work on lined paper: part b defines the symmetric spatial wavefunction as one over root two times the sum of two product wavefunctions, with the spin-1 and spin-0 total spin values and their S-squared eigenvalues.


part iii. looks mathematically brutal =(
 
Attachments
  • Handwritten work on lined paper, rotated sideways: part b with symmetric and antisymmetric spatial wavefunctions psi_23 built from sin(pi x/a) terms, followed by spin results S=1 and S=0 with S^2 values and symmetric or antisymmetric labels.
    ImageUploadedByPhysics Forums1430714973.395573.jpg
    29.3 KB · Views: 545
Last edited:
Is this the correct integral to evaluate for part iii? This looks brutal

Handwritten notes on lined paper: spin results S=1 with S^2 = 2h-bar^2 labelled antisymmetric spatial, and S=0 with S^2 = 0 labelled symmetric spatial; below, an expectation value of (x1 - x2)^2 and a symmetric double integral over x1 and x2 of x^2 times a bracketed sum of sine-squared and sine product terms.
 
Attachments
  • Handwritten notes on lined paper: spin values S=1 and S=0 with S^2 formulas, labelled symmetric and antisymmetric spatial; then the expression <(x1-x2)^2> and a double integral over x1 and x2 from 0 to a of sin-squared and sine product terms.
    ImageUploadedByPhysics Forums1430717022.168788.jpg
    28.7 KB · Views: 514
Maylis said:
I just posted my new work since writing out all the latex and editing it until it is right can take over an hour. But for part (a), I got

$$\langle (x_{1} - x_{2})^{2} \rangle = \frac {a^{2}}{27} - \frac {5a^{2}}{72 \pi^{2}}$$

At least the units are correct for this answer

I don't believe your answer is correct. Note, for example, that in calculating ##\left <x^2 \right >_2## your change of variable led to an integral from 0 to ##2 \pi##. But Griffiths' integrals are from 0 to ##\pi##.
 
True, but that's why I just multiplied by 2, because the function is even, so I should be able to just double my answer, and mulitplied by 3 for the integral from 0 to ##3 \pi##
 
According to Mathematica, the expectation value in part (i) works out to ##\frac{a^2}{6}\left(1-\frac{13}{12\pi^2}\right)##.
 
Then it has something to do with my limits of integration. I went to the prof's office hours and he said to just double what you would get integrating from 0 to pi since it is an even function. I assumed the same for tripling.
 
Your prof must not have looked at your integral carefully. Changing Griffiths' limits of [0, ##\pi##] to [0, ## 2\pi##] will not double the result. A sketch of the integrand would help see that. I agree with vela's result (also using Mathematica).