Jtdylktuy

I'm currently working on the following problem for my Quantum Mechanics course, I could use some help with one of the integrals involved.

A particle of mass $m$ and charge $q$ is confined in a one-dimensional infinite square well of width $a$ in the $x$-direction. The particle starts at its ground state $(n=1)$. Under electromagnetic wave radiation with frequency $omega$ and an electric field strength $Ehat x$, the particle has a probability to transit to higher energy levels.

What is the transition probability to the $n=2$ state? Express your answer in terms of frequency and time. Is the transition to $n=3$ allowed? Can you give a general transition rule?

The following is the work I have done. I will point out my question when the time comes. If you see a mistake in my work prior to my question, I'd appreciate some help with that as well.

So,
$$ vec E=Ehat x Rightarrow H' = -qE_0 x $$
$$ Rightarrow H'_{ba}=qlanglepsi_n|x|psi_mrangle E_0 = -rho E_0$$

For $rho=qlanglepsi_n|x|psi_mrangle$. The transition probability between states $n$ and $m$ is:
$$ c^{(1)}_b approx frac{-i}{hbar} int_0^t H'_{ba}e^{iomega_0 t'}dt'= frac{i}{hbar} rho E_0 int_0^t e^{iomega_0 t'}dt'= frac{q}{hbar omega_0}rho E_0(e^{iomega_0 t} -1)$$

The following part is where I could use some help. I'm not sure how to evaluate the integral for $rho$.

$$langlepsi_n|x|psi_mrangle = int_0^a bigg ( sqrt{frac{2}{a}} sin(frac{npi x}{a})bigg) x bigg ( sqrt{frac{2}{a}} sin(frac{mpi x}{a})bigg) \ = frac{2}{a} int_0^a xsin(frac{npi x}{a})sin(frac{mpi x}{a})dx $$

I know that this integral is zero for $m$ is odd, which corresponds to even transition jumps in the physics. I'm just not sure how to calculate the general case for when $m$ is even (odd transition jumps). For the $n=2$ and $n=3$ states, it is easy enough to plug these values into WolframAlpha and get the correct probabilities, but I'd like to show the general equation before moving onto those parts.

Much appreciated.

It's easiest to do this kind of integral via tabular integration because one factor is a polynomial and the other is easy to integrate.
$$frac2aint_0^axsinleft(frac{npi x}aright)sinleft(frac{mpi x}aright)dx=frac1aint_0^axleft[cosleft(frac{(n-m)pi x}aright)-cosleft(frac{(n+m)pi x}aright)right]dx$$
Now for our table:
$$begin{array}{ccc}frac xa&&cosleft(frac{(n-m)pi x}aright)-cosleft(frac{(n+m)pi x}aright)\
&searrow+&\
frac1a&&frac a{(n-m)pi}sinleft(frac{(n-m)pi x}aright)-frac a{(n+m)pi}sinleft(frac{(n+m)pi x}aright)\
&searrow-&\
0&&-frac{a^2}{(n-m)^2pi^2}cosleft(frac{(n-m)pi x}aright)+frac{a^2}{(n+m)^2pi^2}cosleft(frac{(n+m)pi x}aright)end{array}$$
So you start with the polynomial in the left column and the other factor on the right. You differentiate the polynomial until it's zero and integrate the other factor as many times. Then you multiply pairs diagonally downwards with alternating signs as indicated by the arrows. After adding up results we get
$$begin{align}frac2aint_0^axsinleft(frac{npi x}aright)sinleft(frac{mpi x}aright)dx&=left[frac x{(n-m)pi}sinleft(frac{(n-m)pi x}aright)-frac x{(n+m)pi}sinleft(frac{(n+m)pi x}aright)right.\
&left.+frac{a}{(n-m)^2pi^2}cosleft(frac{(n-m)pi x}aright)-frac{a}{(n+m)^2pi^2}cosleft(frac{(n+m)pi x}aright)right]_0^a\
&=left[frac{a}{(n-m)^2pi^2}-frac{a}{(n+m)^2pi^2}right]left[(-1)^{n+m}-1right]\
&=frac{4nmaleft[(-1)^{n+m}-1right]}{(n^2-m^2)^2pi^2}end{align}$$
If you try doing that integral by parts without a way to efficiently lay out computations like tabular integration, it gets very messy.

The symmetry is broken when temperature is lowered. What I am trying to say that probability of transition should consider temperature effect as well.