An electric dipole transition is the dominant effect of an interaction of an electron in an atom with the electromagnetic field.
Following reference,[1] consider an electron in an atom with quantum Hamiltonian
H0
{E}({r},t)=E0{\hat{z}}\cos(ky-\omegat), {B}({r},t)=B0{\hat{x}}\cos(ky-\omegat).
Write the Hamiltonian of the electron in this electromagnetic field as
H(t) = H0+W(t).
Treating this system by means of time-dependent perturbation theory, one finds that the most likely transitions of the electron from one state to the other occur due to the summand of
W(t)
WDE(t)=
qE0 | |
me\omega |
pz\sin\omegat,
where
q
me
H0+WDE(t)
Between certain electron states the electric dipole transition rate may be zero due to one or more selection rules, particularly the angular momentum selection rule. In such a case, the transition is termed electric dipole forbidden, and the transitions between such levels must be approximated by higher-order transitions.
The next order summand in
W(t)
WDM(t)=
q | |
2me |
(Lx+2Sx)B0\cos\omegat
and describes magnetic dipole transitions.
Even smaller contributions to transition rates are given by higher electric and magnetic multipole transitions.
One way of modelling and understanding the effect of light (mainly electric field) on an atom is to look at a simpler model consisting of three energy levels. In this model, we have simplified our atom to a transition between a state of 0 angular momentum (
Jg=0)
Je=1
n=1,\ell=1
In order to understand the effect of the electric field on this simplified atom we are going to choose the electric field linearly polarized with the polarization axis to be parallel with the axis of the
|g\rangle
|e\rangle
\hat{z}
With this in mind, we can limit ourselves to just the transition from
|g\rangle
|e\rangle
\vec{E}(\vec{r},t)=\hat{z}l{E}(\vec{r})\cos(wt-\phi(\vec{r}))
\hat{z}
w
\phi
l{E}
Now, the main question we want to solve is what is the average force felt by the atom under this kind of light? We are interested by
f=\overline{\langle\hat{F}\rangle}=\overline{\langle-\nabla\hat{V}e.d.(\hat{r},t)\rangle}
\hat{V}e.d.
This potential can be further be written as
\hat{V}e.d=-\hat{D} ⋅ \vec{E}
\hat{D}
The reason we use a two state model is that it allow us to write explicitly the dipole transition operator as
\hat{D}=d0(|e\rangle\langleg|+|g\rangle\langlee|)
\hat{V}e.d=-d0l{E}(\hat{\vec{r}})(|e\rangle\langleg|+|g\rangle\langlee|)\cos(\omegat-\phi(\hat{\vec{r}}))=-d0l{E}(\hat{\vec{r}})(\hat{\sigma}++\hat{\sigma}-)\cos(\omegat-\phi(\hat{\vec{r}}))
f(\vec{r})=\overline{\langled0(\hat{\sigma}++\hat{\sigma}-)\nabla[l{E}(\vec{r})\cos(\omegat-\phi(\hat{\vec{r}})]\rangle}=\overline{d(\vec{r},t)\nablal{E}(\vec{r},t)}=\overline{f(\vec{r},t)}
Now, the semi-classical approach means that we write the dipole moment as the polarizability of the atom times the electric field:
d(\vec{r},t)=\alpha(\omega)l{E}(\vec{r},t)
And as such
f(\vec{r},t)=\alpha(\omega)l{E}(\vec{r},t)\nablal{E}(\vec{r},t)=
1 | |
2 |
\alpha(\omega)\nablal{E}2(\vec{r},t)
f(\vec{r})=\nabla[
1 | |
2 |
\alpha(\omega)\overline{l{E}2(\vec{r},t)}]=\nabla[
1 | |
4 |
\alpha(\omega)l{E}2(\vec{r})]
Ve.d=-
1 | |
4 |
\alpha(\omega)l{E}2(\vec{r})
Before progressing in the math, and trying to find a more explicit expression for the proportionality constant
\alpha(\omega)
I\propto\overline{l{E}2(\vec{r},t)}
Now, let's look at how to get the expression of the polarizability
\alpha(\omega)
We will use the density matrix formalism, and the optical Bloch equations for this.
The main idea here is that the non-diagonal density matrix elements can be written as
\rhoeg=Tr(\hat{\sigma}-\hat{\rho})
\rhoge=Tr(\hat{\sigma}+\hat{\rho})
d(\vec{r},t)=Tr(\hat{D}\hat{\rho}(t))=Tr(d0(\hat{\sigma}++\hat{\sigma}-)\rho(t))=d0[Tr(\hat{\sigma}+\rho(t))+Tr(\hat{\sigma}-\rho(t))]=d0(\rhoge+\rhoeg)
Here is where the optical Bloch equations will come in handy, they give us an equation to understand the dynamics of the density matrix.
Indeed, we have:
d\hat{\rho | |
(t)}{dt} |
=
1 | |
i\hbar |
[\hat{H},\rho]
and another term that describes the spontaneous emissions of the atom:
d\hat{\rho | |
eg |
(t)}{dt}=-
\Gamma | |
2 |
\hat{\rho}eg,
d\hat{\rho | |
ge |
(t)}{dt}=-
\Gamma | |
2 |
\hat{\rho}ge
Where
\hat{H}
\hat{H}=\hbar\omega0|e\rangle\langlee|+Ve.d.
\Gamma=
| |||||||
3\pi\epsilon0\hbarc3 |
\Gamma
\Gamma-1
Let us introduce the Rabi frequency
\Omega
\hbar\Omega=-d0l{E}(\vec{r})
Then we can write the optical Bloch equations for
\rhoeg
\rhoge
For this part we take the equation of the evolution of the
\rho
i\hbar
d\hat{\rho | |
eg |
(t)}{dt}=\hbar(\omega0-i
\Gamma | |
2 |
)\rhoeg-d0l{E}(\vec{r},t)(\rhogg-\rhoee)
We can get the equation for
\rhoge
We can then repeat the process for all 4 matrix elements, but in our study we will apply a small field approximation, so that the electric field is small enough that we can uncouple the 4 equations. This approximation is written mathematically using the Rabi frequency as:
\Omega\ll|\Delta|,\Gamma
\Delta=\omega-\omega0
Then we can neglect
\rhoee=0
\rhogg=1
l{E}
\rhoee=0
\rhogg=1
We can then rewrite the evolution equation to:
d\hat{\rho | |
eg |
(t)}{dt}=
-(i\omega | ||||
|
)\rhoeg-i\Omega\cos(wt-\phi)
This is an ordinary first order differential equation with an inhomogeneous term in cosines. This can easily be solved by using the Euler's formula for the cosine.
We get the following solution:
\rhoeg(t)=\Omega/2(
e-i(\omega | |
\omega-\omega0+i\Gamma/2 |
-
ei(\omega | |
\omega+\omega0-i\Gamma/2 |
)
Furthermore, if we say that the detuning
\Delta
\Gamma
\omega,\omega0
\Gamma
\rhoeg(t)=\Omega/2(
e-i(\omega | |
\omega-\omega0 |
-
ei(\omega | |
\omega+\omega0 |
)
and
\rhoge(t)=\Omega/2(
ei(\omega | |
\omega-\omega0 |
-
e-i(\omega | |
\omega+\omega0 |
)
And coming back to our average dipole moment:
d(t)=d0\Omega/\overline{\Delta}\cos(\omegat-\phi)=-
| ||||||||||
(\vec{r})}{\hbar |
\overline{\Delta}}\cos(\omegat-\phi)
1 | |
\overline{\Delta |
Then it is clear that
d(t)=-
| |||||||
\hbar\overline{\Delta |
\alpha(\omega)=-
| |||||||
\hbar\overline{\Delta |
Finally, we can write the potential felt by the atom in due to the electric dipole interaction as:
Ve.d.(\vec{r})=
| |||||||
4\hbar\overline{\Delta |
The essential points worth discussing here are as said previously that the light intensity
I\propto\overline{l{E}2(\vec{r},t)}
\overline{\Delta}
\omega<\omega0
\omega>\omega0