74 CHINESE OPTICS LETTERS / Vol. 8, No. 1 / January 10, 2010
First-principles study on the optical properties
for CsI crystal with a pair of V
1−
Cs
-V
1+
I
Tingyu Liu (444ÉÉÉ)
1,2∗
, Feinan Yan (îîîIII)
1,2
, Jun Chen ( ddd)
1,2
,
and Liping Liang (ùùùwww±±±)
1,2
1
College of Science, University of Shanghai for Science and Technology, Shanghai 200093, China
2
Shanghai Key Laboratory of Contemporary Optics System, Shanghai 200093, China
∗
E-mail: liutyyxj@163.com
Received March 16, 2009
The electronic structures and optical properties of both the perfect CsI crystal and the crystal containing
a pair of V
1−
Cs
-V
1+
I
are calculated using CASTEP code with the lattice structure optimized. The calculated
results indicate that the optical symmetry of the CsI crystal coincides with the lattice structure geometry
of the CsI crystal. The absorption spectrum of the CsI crystal containing a pair of V
1−
Cs
-V
1+
I
also does not
o ccur in the visible and near-ultraviolet range. It reveals that the existence of the pair of V
1−
Cs
-V
1+
I
in CsI
crystal has no visible effects on the optical properties of the CsI crystal.
OCIS co des: 160.4670, 160.2220, 160.4760.
doi: 10.3788/COL20100801.0074.
Cesium iodide crystal is widely used as a scintillator
material in particle detectors due to its high emission
efficiency and fast decay time, especially the crystal
doped with thallium
[1−4]
. In general, this kind of sim-
ple cubic structured crystal would be easily damaged by
high-energy irradiation, such as X- or γ-ray, resulting in
the creation of color centers and making the crystal occur
absorption bands in the visible range.
The Kr ion irradiated CsI
[5]
exhibits several absorption
bands peaking at 240, 335, 410, 560, and 751 nm with
a shoulder at 716 nm. The absorption bands of color
centers are harmful to the scintillate properties. Many
researchers have discussed the origins of these absorp-
tion bands
[5−8]
. The CsI crystal after the irradiation
exhibits six absorption bands peaking at 375, 560, 750,
1020, 1180, and 285 nm
[9]
. The former five bands are
assigned to I
1−
2
color center, H band, F band, I
2−
3
color
center, and F
2
color center, respectively. The origin of
the absorption band at 285 nm is still unclear
[9]
. The
growth of optical absorption spectra in a high-Tl crystal
over a wide range of γ-ray dose show a group of six peaks,
which are at 355, 390, 430, 460, 520, and 560 nm. In
the infrared region, 840-nm absorption band originating
from the F-band and 980-nm absorption band related to
the F
3
-band are also observed
[10]
. The short wavelength
absorption bands are generally assigned to V band and
the long wavelength absorption bands are assigned to F
bands
[6]
.
On the other hand, there are a large number of suc-
cessful simulations on the optical prop erties of crystals
using the density functional theory (DFT)
[11,12]
with the
code called CASTEP. In order to further study the phys-
ical properties of the pair of V
1−
Cs
-V
1+
I
related defects, the
electronic structures, dielectric function, complex refrac-
tive index, and absorption spectra of both the perfect CsI
cluster and the CsI cluster containing a pair of V
1−
Cs
-V
1+
I
under the irradiation of a polarized light are calculated
using CASTEP code with the lattice structure optimized.
The lattice structure optimization is required before we
study the optical properties of the crystal due to the lat-
tice distortion caused by the existence of V
1−
Cs
-V
1+
I
.
The CsI single crystal is simple cubic structured. The
CsI supercell used here consists of 27 Cs and 27 I centered
at I, having the chemical formula of [Cs
27
I
27
]. The CsI
supercell containing a pair of V
1−
Cs
-V
1+
I
has a Cs vacancy
V
1−
Cs
nearest centered I ion and an iodine V
1+
I
replacing
the centered I
1−
ions.
The CASTEP code is used for the lattice optimiza-
tion for the CsI containing one or two vacancies and the
calculation of the ground-state electronic structure us-
ing a plane-wave pseudo-potential formulation within the
framework of DFT with generalized gradient correction
in the form of PW91
[13−15]
. We use ultrasoft pseudo-
potentials for the iodine and cesium atoms and a plane-
wave cutoff energy of 340 eV. The basic parameters are
chosen as follows: kinetic energy cutoff of 340 eV, fast
Fourier transform (FFT) grid dimensions of 60×60×64,
space representation being reciprocal, and self-consistent
field (SCF) tolerance 1.0×10
−6
of eV/atom. Optimal
atomic positions are determined until satisfying the con-
ditions: 1) the maximal force on them is smaller than
0.5 eV/nm; 2) the maximal change of energy per atom
is smaller than 0.00001 eV; and 3) the maximal displace-
ment is smaller than 0.0001 nm. All other calculations
are performed on basis of the lattice structure being op-
timized.
The dielectric function of an anisotropic material is a
complex symmetric second-order tensor which describes
the linear response of an electronic system to an applied
external electric field. The imaginary part of the dielec-
tric tensor is directly related to the electronic band struc-
ture of a solid, so it can be computed from the knowledge
of single-particle orbitals and energies approximated by
the solutions of the Kohn-Sham equations. However, it is
a well known fact that DFT calculations underestimate
the band gap. To take this into account, a ‘scissors oper-
ator’ is used, allowing a shift of the bands situated above
the valence band and a rescaling of the matrix elements.
1671-7694/2010/010074-04
c
° 2010 Chinese Optics Letters