# Spontaneous Emission of an Excited Atom in a Dusty Unmagnetized Plasma Medium

Authors: Naser Alinejad; Noushin Pishbin
Journal: Advances in Materials Science and Engineering
 (2014)

---

## Abstract

Investigation of spontaneous decay of an excited atom in dusty unmagnetized plasma is presented in this paper. The transverse contribution to the decay rate is normally associated with spontaneous emission. The rate of spontaneous emission can be obtained by Fermi’s golden rule. In this calculation, the transverse contribution to dielectric permittivity and Green function technique are used. Calculation of the decay rate of atoms is applicable to understand the particular structure of the vacuum state of the electromagnetic field.


---

## Body

## 1. Introduction

Spontaneous decay of an excited state arises from the interaction between an excited atom or molecule and the ground state of the quantized electromagnetic field [1]. The rate of spontaneous emission depends partly on the environment of a light source. This means that, by placing the light source in a special environment, the rate of spontaneous emission can be modified [2]. Spontaneous emission is a quantum effect, which in a semiclassical picture can be described as an emission which is stimulated by vacuum noise, that is, by the zero point fluctuations of the electromagnetic field [3]. In a well-known paper published in 1948, Welton [4] wrote that spontaneous emission “can be thought of as forced emission taking place under the action of the fluctuating field.”Spontaneous emission of light or luminescence is a fundamental process that plays an essential role in many phenomena in nature and forms the basis of many applications, such as fluorescent tubes, older television screens (cathode ray tubes), plasma display panels, lasers (for startup normal continues operation works by stimulated emission instead), and light emitting diodes.Spontaneous emission is ultimately responsible for most of the light around us. We would not be here without it. The first person to derive the rate of spontaneous emission directly from the first principles was Dirac [5], who used the newly formulated quantum theory of radiation.The rate of spontaneous emission (i.e., the radiative rate) can be described by Fermi’s golden rule [6]; the rate of emission depends on two factors: an “atomic part,” which describes the internal structure of light source, and a “field part,” which describes the density of electromagnetic modes of the environment. The atomic part describes the strength of a transition between two states in terms of transition moments. In a homogeneous medium, such as free space, the rate of spontaneous emission in the dipole approximation is given by [7]
(1)
Γ
0
=
ω
3
d
2
3
π
ℏ
ε
0
c
3
.Clearly, the rate of spontaneous emission in free space increases withω
3. In contrast with atom, which has a discrete emission spectrum, quantum dots can be tuned continuously by changing their size. This property has been used to check the ω
3-frequency dependence of the spontaneous emission rate as described by Fermi’s golden rule.The purpose of present work is to investigate the decay rate of an excited atom in dusty unmagnetized plasma, using Green function technique. Dusty plasma is loosely defined as normal electron-ion plasma with an additional charged component which increases the complexity of the system even further. Dusty plasma is low temperature fully or partially ionized electrically conducting gases whose constituents are electron, ions, charged dust grains, and neutral atoms. Dust grains are massive and their sizes range from nanometers to millimeters. Dust grains may be metallic, conducting, or made of ice particulates [8].The plasma frequencyω
p in uniform, cold, dusty unmagnetized plasma can be described as follows [8]:
(2)
ω
P
2
=
∑
s
4
π
n
s
0
q
s
2
m
s
=
∑
s
ω
p
s
2
,
where n
s
0, q
s, and m
s are, respectively, the unperturbed number density, the charge, and the mass of the plasma species s (s equals e for electrons, i for ions, and d for dust grains). And ω
p
s
=
(
4
π
n
s
0
q
s
2
/
m
s
)
1
/
2 represents the plasma frequency associated with the plasma species s.Considering the usual expression for the decay rate of an excited atom in dipole approximation [9] and since we are interested in the decay rate of an excited atom in dusty unmagnetized plasma, at the first stage, we calculate Green’s function in this medium.


## 2. Green Function

The associated electromagnetic fields are determined by Maxwell’s equations(3)
∇
×
E
=
-
μ
0
∂
H
∂
t
,
∇
×
H
=
j
+
∂
D
∂
t
,
where E, H, and the electric displacement D have positive and negative frequency decomposition as follows:
(4)
E
^
(
r
,
t
)
=
∫
0
∞
d
ω
{
E
^
+
(
r
,
ω
)
e
-
i
ω
t
+
E
^
-
(
r
,
ω
)
e
+
i
ω
t
}
.
The positive and negative frequency parts in the integrand involve only annihilation and creation operators, respectively.Taking the Fourier transform of these equations with respect tot, we obtain for positive frequency part of the fields
(5)
∇
×
E
+
(
r
,
ω
)
=
i
ω
B
+
(
r
,
ω
)
,
∇
×
H
+
(
r
,
ω
)
=
-
i
ω
D
+
(
r
,
ω
)
+
J
+
(
r
,
ω
)
.
Using the relationship between the displacement and the electric field
(6)
D
i
+
(
r
,
t
)
=
∫
-
∞
t
d
t
′
∫
d
r
′
ε
i
j
(
r
-
r
′
,
t
-
t
′
)
E
j
(
r
′
,
t
′
)
,
D
i
(
r
,
ω
)
=
ε
i
j
(
r
,
ω
)
E
j
(
r
,
ω
)
.
And utilize the gauge in which the scalar potential vanishes:
(7)
E
+
(
r
,
ω
)
=
i
ω
A
+
(
r
,
ω
)
.We find the inhomogeneous wave equation forA
+
(
r
,
ω
) in the frequency domain
(8)
∇
×
∇
×
A
+
(
r
,
ω
)
-
q
2
n
2
A
+
(
r
,
ω
)
=
μ
(
ω
)
ε
0
c
2
J
+
(
r
,
ω
)
,
where n
=
ε
(
r
,
ω
)
μ
(
ω
), c
2
=
1
/
ε
0
μ
0, and q
=
ω
/
c. Since the medium is unmagnetized, then μ
(
ω
)
=
1.This equation can be solved by using Fourier time transform Green function as(9)
A
+
(
r
,
ω
)
=
∫
d
3
r
′
G
(
r
-
r
′
,
ω
)
·
J
+
(
r
′
,
ω
)
J
+
(
r
,
ω
)
=
∫
d
3
r
′
δ
(
r
-
r
′
)
J
+
(
r
′
,
ω
)
.Substitution of (9) into (8) shows that the classical Green function G
(
r
,
r
′
,
ω
) satisfies the differential equation:
(10)
∇
×
∇
×
G
(
r
-
r
′
,
ω
)
-
q
2
ε
(
k
,
ω
)
G
(
r
-
r
′
,
ω
)
=
I
ε
0
c
2
δ
(
r
-
r
′
)
,
where I is the unit tensor.For solving this equation, we used the formal approach that relies on Fourier transform, where the differential equations are transformed into reciprocal space. At the end, after algebraic calculation and separating Green function to transverse and longitudinal component in coordinate of wave vector [7], we obtain
(11)
G
i
j
T
(
k
,
ω
)
=
(
δ
i
j
-
k
i
k
j
/
k
2
)
ε
0
ω
2
[
(
k
c
/
ω
)
2
-
ε
t
]
,
(12)
G
i
j
L
(
k
,
ω
)
=
-
k
i
k
j
k
2
ε
0
ε
l
ω
2
.The dielectric tensor for dusty unmagnetized plasma is obtained by Poisson-Maxwell equation [10]. Consider
(13)
ε
i
j
(
k
,
ω
)
=
k
i
k
j
k
2
ε
l
(
k
,
ω
)
+
(
δ
i
j
-
k
i
k
j
k
2
)
ε
t
(
k
,
ω
)
+
(
δ
i
j
-
Ω
i
Ω
j
Ω
2
)
ε
d
(
k
,
ω
)
,
where Ω
i is the angular velocity of the ith grain. The influence of the dust grain rotation is described by ε
d
(
k
,
ω
). Since in this paper the unmagnetized plasma (B
0
=
0) is investigated, the influence of the dust grain rotation is neglected. Therefore ε
d
(
ω
,
k
)
=
0.Alsoε
l
(
k
,
ω
) and ε
t
(
k
,
ω
) are the longitudinal and the transverse dielectric permittivity, respectively. They are given by
(14)
ε
l
(
k
,
ω
)
=
1
+
∑
s
ω
p
s
2
k
2
V
T
s
2
[
1
-
I
+
(
ω
k
V
T
s
)
]
,
ε
t
(
k
,
ω
)
=
1
-
∑
s
ω
p
s
2
ω
2
I
+
(
ω
k
V
T
s
)
,
where V
T
s
=
(
k
B
T
s
/
m
s
)
1
/
2 is the thermal speed of the species s and s includes ion, electron, and dust grain.The functionI
+
(
x
) is
(15)
I
+
(
x
)
=
x
2
π
∫
-
∞
∞
d
z
exp
⁡
(
-
z
2
/
2
)
x
-
z
.The asymptotic forms of (15) are as follows.For|
x
|
≫
1, |
Re
x
|
≫
|
Im
⁡
x
| and Im
⁡
x
<
0,
(16)
I
+
(
x
)
≈
1
+
1
x
2
+
⋯
-
i
π
2
x
exp
⁡
⁡
(
-
x
2
2
)
.And for|
x
|
≪
1,
(17)
I
+
(
x
)
≈
-
i
π
2
x
.Substituting the transverse dielectric permittivity in (11), we find that the transverse Green function has the following expansion:
(18)
G
i
j
T
(
k
,
ω
)
=
(
δ
i
j
-
k
i
k
j
/
k
2
)
ε
0
ω
2
[
(
k
c
/
ω
)
2
-
(
1
-
∑
s
(
ω
p
s
2
/
ω
2
)
I
+
(
ω
/
k
V
T
s
)
)
]
.Here, we study the range of high frequencies; then(19)
|
ω
k
V
T
s
|
≫
1
,
and
then
I
+
(
ω
k
V
T
s
)
≈
1
.By using this approximation and taking inverse Fourier transformation, the transverse Green function in coordinate space is obtained:(20)
G
i
j
T
(
R
,
ω
)
=
1
4
π
ε
0
ω
2
(
1
-
ω
p
2
/
ω
2
)
×
{
[
δ
i
j
-
3
R
i
R
j
R
2
]
+
q
′
3
[
(
1
q
′
R
+
i
(
q
′
R
)
2
-
1
(
q
′
R
)
3
)
δ
i
j
-
(
1
q
′
R

+

3
i
(
q
′
R
)
2

-

3
(
q
′
R
)
3
)
R
i
R
j
R
2
]
×
e
i
q
′
R
[
δ
i
j
-
3
R
i
R
j
R
2
]
}
,
in which R
=
r
-
r
′ and q
′
=
(
ω
/
c
)
1
-
ω
p
2
/
ω
2.Expanding the exponential term in (20), it is straightforward to show that, in the limit R
→
0, the transverse Green function takes the form
(21)
lim
⁡
R
→
0
⁡
G
i
j
T
(
R
,
ω
)
=
1
4
π
ε
0
c
2
{
R
i
R
j
2
R
3
+
δ
i
j
2
R
+
2
i
ω
3
c
(
1
-
ω
p
2
ω
2
)
δ
i
j
+
o
(
R
)
}
.The singular behavior of the real part arises from the transverse delta function associated with the canonical field commutation relation [1]. But the imaginary part of this Green function is well behaved in the limit as r tends to zero:
(22)
Im
⁡
G
i
j
T
(
0
,
ω
)
=
Real
(
(
1
-
ω
p
2
/
ω
2
)
)
ω
6
π
ε
0
c
3
δ
i
j
.Substituting the longitudinal dielectric permittivity in Green function, the longitudinal Green function is obtained. But because of considering only the transverse radiative modes of the electromagnetic field, the calculation of longitudinal component is omitted here.


## 3. Calculation of Decay Rate

According to the expression of decay rate [7],
(23)
Γ
=
2
ω
2
ℏ
Im
⁡
{
d
i
G
i
j
(
0
,
ω
)
d
j
}
.
The transverse contribution to the total transition rate can be written as follows.For|
ω
/
k
V
T
s
|
≫
1,
(24)
Γ
T
=
Real
(
1
-
ω
p
2
ω
2
)
Γ
0
.In this relation, if(25)
ω
p
2
>
ω
2
,
then (
1
-
ω
p
2
/
ω
2
) is imaginary and Γ
T
=
0. This is linked to the nonexistence of propagating modes in the plasma at this range of frequency into which a photon can be emitted.And if(26)
ω
p
2
<
ω
2
,
then (
1
-
ω
p
2
/
ω
2
)is real and
(27)
Γ
T
=
(
(
1
-
ω
p
2
ω
2
)
)
Γ
0
.The result of calculation of transverse decay rate is applied in high energy astrophysical plasma, collisional radioactive model in plasma, fluorescence quenching, and other options.In Figure1 we plot the transverse decay rate as given in (24).Figure 1
The transverse spontaneous emission rate (24). Transverse decay rate is zero for the range of frequency ω
<
ω
P and for ω
>
ω
P; the diagram corresponds to (27). The plasma frequency in this model is considered 5.5
×
10
10 Hz.


## 4. Conclusions

In this paper, we have calculated the transverse contribution of decay by using Green function technique and fluctuation dissipation theorem [11] in dusty unmagnetized plasma. Spontaneous emission is the process by which a light source such as an atom, molecule, or nucleus in an excited state undergoes a transition to a state with a lower energy and emits a photon.The transverse decay rate in dusty plasma is obtained by modifying the dielectric permittivity.This behavior arises from the existence of the grain and radiating interaction of them with electrical field.



---

*Source: 746742-2014-03-31.xml*