Analytical expressions for the effective local force constants, correlated Einstein frequency and temperature, first cumulant or net thermal expansion, second cum ulant or Debye Waller f
Trang 1C A L C U L A T IO N O F X A F S CƯ M ULANTS F O R F C C CRYSTALS
C O N T A IN IN G I M P U R I T Y A TOM
N g u y e n V a n H u n g , N g u y e n T h i T h u H o a i , L e H a i H u n g
D e p a r tm e n t o f P h y sic s, College o f Science, V N U
A b s tr a c t: A new procedure for calculation and evaluation of XAFS cumulants of fee crystals containing impurity atom has been developed based on the quantum statistical theory with correlated Einstein model Analytical expressions for the effective local force constants, correlated Einstein frequency and temperature, first cumulant or net thermal expansion, second cum ulant or Debye Waller factor and third cumulant of fee crystals containing impurity atom have been derived Morse potential param eters of pure crystals and those with impurity included in the derived expressions have been calculated Numerical results for Cu Ni Ni-Cu are found to be
in good agreement with experiment
I n t r o d u c t i o n
T h e c r y s t a l s w ith fee s t r u c t u r e occupies a b o u t 25 % of e le m e n ts in the
)erio d ical M e n d e le e v s y s te m T h a t is w hy th e c a lc u la tio n of p h y s ic a l p a r a m e t e r s of
h e s e c r y s ta ls is v ery i m p o r t a n t To s tu d y t h e r m o d y n a m i c p r o p e r t i e s of a s u b s ta n c e
t is n e c e s s a r y to i n v e s t i g a t e its effective local force c o n s t a n t s , c o r r e la te d E in s te in
r e q u e n c y a n d t e m p e r a t u r e , n e t t h e r m a l e x p a n s io n , m e a n s q u a r e re la tiv e lis p la c e m e n t (MSRD) or Debye W a lle r fac to r a n d t h i r d c u m u l a n t [1-13] w hich arc*
o n t a i n e d in th e X -ra y a b s o r p tio n fine s t r u c t u r e (XAFS) [9J M o reo ver, th e im p u rity
>r d o p a n t a to m c a n in fu e n c e on t h e p h y s ic a l p a r a m e t e r s t a k e n from th e XAFS
p e c t r a 110] a n d on t h e efficiency of u s in g t h e s e s u b s t a n c e s T h e rm o d y n a m ic
>roperties of a lk a li m e t a l s u n d e r in flu e n c e of i m p u r i t y h a s b e e n s tu d i e d [1 1]
T he p u r p o s e of t h i s w ork is to develop a m e th o d for c a lc u la tio n a n d e v a lu a tio n
>f t h e effective local force c o n s t a n t s , c o r r e l a t e d E i n s t e i n freq u e n cy an d
e m p e r a t u r e , f ir s t c u m u l a n t or n e t t h e r m a l e x p a n s io n , seco n d c u m u l a n t w hich IS
•qual to M SRD or Debve W a lle r fa c to r a n d t h i r d c u m u l a n t of fee c r y s t a l s c o n ta in in g
I d o p a n t or i m p u r i t y (I) a to m as a b s o rb e r in th e X A FS p ro cess Its n e a r e s t
le ig h b o rs a r e th e h o s t (H) a to m s T h e d e r iv a tio n is b a s e d on th e q u a n t u m statistic a l th e o ry w ith t h e c o r r e l a t e d E i n s t e i n m odel [7] w h ic h is c o n sid e re d at
j r e s e n t as “th e b e s t t h e o r e t ic a l f ra m e w o r k w ith w hich th e e x p e r i m e n t a l i s t can
e la te force c o n s t a n t s to t e m p e r a t u r e d e p e n d e n t XAFS" [1 0] Kor c o m p le tin g th e ab nitio c a lc u la tio n p r o c e d u re th e p a r a m e t e r s of M orse p o t e n t i a l of p u r e c r y s ta ls a n d hose w ith i m p u r i t y h a v e b ee n also c a lc u la te d N u m e r ic a l c a l c u l a t i o n s for Cu Ni
m d Cu doped by Ni a to m h av e b e e n c a r r ie d o u t in c o m p a r is o n to th o se of Lh(' p u re
n a t e r i a l s to show t h e r m o d y n a m i c a l effects of fee c r y s t a l u n d e r in flu en c e of th e
m p u r ity ato m T h e c a l c u la te d r e s u l t s a r e found to be in good a g r e e m e n t w ith 'X perim ent for M orse p o t e n t i a l a n d for th e o t h e r t h e r m o d y n a m i c p a r a m e t e r s [13]
Trang 2C al cu la ti on o f XAFS c u m u l a n t s for fee. 9
2 F o r m a l i s m
T h e e x p r e s s io n for th e M SRD in XAFS th e o r y is d e r iv e d b a s e d on th e
a n h a r m o n i c c o r r e la te d E i n s t e i n model [7] acco rd in g to w h ic h th e effective
i n te r a c tio n E i n s t e i n p o t e n t i a l of th e s y s te m c o n s is tin g of a n i m p u r i t y (I) a to m as
a b s o rb e r a n d th e o t h e r h o s t (H) a to m s is given by
i' c/f(x ) - - l<cirx 2 + /vvY? +■■■ - ~ a T x R ' 2 'R, i
M ị +M„
H ere X is d e v ia tio n b e tw e e n th e i n s t a n t a n e o u s bond le n g t h r a n d its e q u ilib riu m
v alu e r „ , kcfJ is effective local force c o n s ta n t, a n d ky t h e cubic p a r a m e t e r giving an
a s y m m e t r y in th e p a i r d i s t r i b u t i o n fun ctio n , R is bo nd u n i t vector T h e c o rre la te d
E i n s t e i n m odel m a y be d efin ed as a o s c illatio n of a p a i r of a t o m s w ith m a s s e s M ị
a n d MH (e.g., of i m p u r i t y a to m a s a b s o r b e r a n d of h o s t a to m as b a c k s c a t t e r e r ) in a
given sy stem T h e ir o s c illa tio n is in flu e n c e d by t h e i r n e ig h b o r s given by th e last
te r m in th e l e f t - h a n d side of Eq (1), w h e re t h e s u m / is o v er a b s o r b e r ( / = 1) and
b a r k s c a t t e r e r 0 = 2), a n d t h e s u m j is over all t h e i r n e a r e s t n e ig h b o rs , excluding
th e a b s o rb e r a n d b a c k s c t t e r e r th e m s e lv e s T h e l a t t e r c o n t r i b u t i o n s a r e d e s c rib e d Ly
th e te r m vm (v).
For w eak a n h a r m o n i c i t y in t h e XAFS p ro cess t h e M o rse p o t e n t i a l is given by
I ho e x p a n s io n
V(x) = d {Lr2 u ' - 2 c ) = d ( - 1 + a 2.x2 - e r V + ■ • ■) (2)
for th e p u r e m a t e r i a l a n d
V/If ( v) “ D jii (-1 + ơ J ịị -V —aj Hx H ) 0 ) for th e case w ith i m p u r ity , w h e r e M o rse p o t e n t i a l p a r a m e t e r s h a v e b e e n obtained
by a v e r a g in g th o s e of t h e p u r e m a t e r i a l s a n d a r e giv en by
r, D l + D ll „ 2 D l a i + D H a H 3 _ D l a ) + D H a hl (1\
U s in g th e d e f in itio n [2, 7] V = A' - ứ as th e d e v ia tio n from t h e e q uilib rium
v a lu e of X th e Eq (1) is r e w r i t t e n in th e su m of th e h a r m o n ic c o n t r ib u t i o n a n d t i e (in h a rm o n ic c o n t r ib u t i o n ỔV a s a p u r t u r b a t i o n
T a k in g in to a c c o u n t th e a to m ic d i s t r i b u t i o n of fee c r y s t a l a n d u s in g th e abo/e
e q u a t i o n s we o b ta in th e effective local force c o n s ta n t
Kir = 2 D m a ///[1 + 3(//,2 +/jị)] + ị o Ha ị = /.ưoị, 6 ,
th e cubic a n h a r m o n i c p a r a m e t e r
Trang 3ky — Dm ct'm (1 + //|' + //2),
th e a n h a r m o n i c c o n t r ib u t i o n to th e effective p o t e n t i a l of t h e s y s te m
SV(y) = 2 D IHaj H(\ + 3( rì + v l ) ) + ^ D flaj f ay - DIHcc]H{\ + r ì + f i ị ) Ý
ơ )
(8)
th e c o r r e la te d E i n s t e i n fre q u e n c y
[f*
a n d th e c o r r e la te d E i n s t e i n t e m p e r a t u r e
Dm a m [1 + 1(MỈ + MÌ )] + ị D Ha l
I 2
(9)
k B [// Dịh&Ĩh D + 3(//f + M i)]+ 2 ^ H a ~H (10)
w here
M , +M// , / / 2 M j + M h / / 2 =
/7
T h e c u m u l a n t s h a v e b een d e riv e d by a v e r a g i n g p ro c e d u re , u s in g th e
s ta t i s t i c a l d e n s ity m a t r i x /9 a n d th e c a n o n ic a l p a r t i t i o n fu n c tio n z in th e form
< y m > = - T r ( p y m), m = 1 ,2 ,3 ,- ,
'Zd
Z = Trp, p = p 0 +ổp, Z * Z 0 =Trpa ,
p 0 =e - fiHo, H 0 = ậ - + ị - k eíry 2, fi = \ / k BT
(13)
(14)
2// 2
where k B is Boltzmann constant and ỏp is neglected due to small anharmonicity in XAFS [2]
U sing th e above r e s u lts we calcu late the second c u m u l a n t or D ebye-W aller factor
(15)
Ơ
Í above r e s u lts we calcu late the second c u m u l a n t or
Debye-w h e re Debye-we e x p r e s s y in t e r m s of a n i h i l a t i o n a n d c r e a tio n o p e r a t o r s , à a n d ( 1+, i e
{ « + « •); ’ *
\JLIO) e
an d use h a r m o n ic o s c illa to r s t a t e I n) w ith e ig e n v a lu e En = nh( 0 E (ig n o rin g th e zero
p o in t e n e rg y for c o n v e n ie n ce )
T h e refo re , th e e x p r e s s io n for seco n d c u m u l a n t (M SRD) or D e b y e -W a lle r fac to r
is r e s u lte d as
ơ ~ = Ơ (1 + z)
(1-z)
PiCOr
Trang 4C al cu la ti on of XAFS c u m u l a n t s for fee. 11
Now we c a lc u la te th e odd c u m u l a n t s
! ^ e- 0 E _ e-pE„.
Ơ
0
(18)
U s in g th e c a lc u la te d m a tr ix e le m e n ts a n d m a t h e m a t i c a l f o r m u la s for d iffe re n t
t r a n s f o r m a t i o n s we o b ta in t h e e x p r e s s io n s for th e firs t c u m u l a n t (m — 1)
i1 + r) J (l)- 3Dw g // / (l + //|3+/<23) ^ - (19)
uyi ) — u () _ x , V () Ị- !
( I )
( l - z )
DiHa ĨH (1 + 3/'|2 + Ml ) + 2
a n d for th e t h i r d c u m u l a n t (m=3)
n O) = (3) (l + ioz + z2) (3) = D ih GIH M\ /^2 ) (20)
( 1 - 2 )
DiHa ìn (1 + 3/'|2 + ) + 9 D Ha~H
In th e above e x p r e s s io n s ơ i '1, <702, ơ i’ 1 a r e z e r o -p o in t c o n t r i b u t i o n s to t h e first, second a n d t h i r d c u m u l a n t s , re sp ectiv e ly T h e y c h a r a c t e r i z e q u a n t u m effects
o ccu rred by u s in g q u a n t u m t h e o r y in o u r c a lc u la tio n a n d in flu e n c e on t h e o b ta in e d
r e s u lts T h e above d e r iv e d c u m u l a n t s a r e c o n ta in e d in t h e XAFS in c lu d in g
a n h a r m o n i c c o n t r ib u t i o n s [9]
3 N u m e r i c a l r e s u l t s a n d c o m p a r i s o n t o e x p e r i m e n t
Now we a p p ly th e d e riv e d e x p r e s s io n s to n u m e r ic a l c a l c u l a t i o n s for Cu, Ni
a n d Cu doped by Ni a to m a s a b s o r b e r in th e XAFS p rocess T h e i r M o rse p o t e n t i a l
p a r a m e t e r s h a v e b e e n c a l c u l a t e d u s in g t h e p r o c e d u re p r e s e n t e d in [16] T h e
c a lc u la te d v a l u e s of M o rse p o t e n t i a l p a r a m e t e r s ; c o r r e la te d E i n s t e i n fre q u e n c y a n d
t e m p e r a t u r e ; effective local force c o n s ta n t for t h e p u r e Cu, Ni a n d th o s e doped by
N i a to m a r e w r i t t e n in T a b le I T h e y a r e found to be in good a g r e e m e n t w ith
e x p e r i m e n t [13]
T a b l e I C a lc u la te d v a l u e s of M o rse p o t e n t i a l p a r a m e t e r s D, a ; c o r r e la te d E i n s t e i n fre q u en cy z , a n d t e m p e r a t u r e eE\effective local force c o n s t a n t keff for N i-N i, C
u-Cu, N i-C u in c o m p a r is o n to e x p e r i m e n t [13] _
Bond D(eV) a ( A ’) r„ ( Ả ) k ef f { N / m ) o > e ( x 1013 Hz )
Ni-Ni, present 0.4263 1.3819 2.8033 65.2158 3' 64 73 278.6038
Ni-Ni, exp.[ 13] 0.4100 1.3900 2.9035 63.4596 3.5979 274.8271
Cu-Cu, present 0.3367 1.3549 2.8701 49.5156 3.0544 233.3151
Cu-Cu, exp.[13] 0.3300 1.3800 2.9802 50.3450 3.0799 235.2661
Ni-Cu, exp.[ 13] 0.3700 1.3855 2.9337 57.8621 3.3693 257.3708
Trang 5F i g u r e 1 i l l u s t r a t e s o u r c a l c u l a t e d M o rse p o t e n t i a l s of C u, Ni a n d of Cu
d o p p e d by Ni a t o m w h ic h a g r e e s w ell w ith t h e e x p e r i m e n t a l r e s u l t [13]
r(A°)
Figure 1: Morse potentials of Cu, Ni and Cu doped by Ni atom com pared to experim ent [13].
Figure 2: T e m p e r a t u r e d ep e n d e n c e of o u r ca lc u la te d first c u m u l a n t c r ^ ( r ) of Cu (dash-
dot), Ni (dash) a n d of Cu doped by Ni ato m (solid) co m p a re d to e x p e r im e n t (dot) [13]
Trang 6Ca lc u la ti o n o f XA FS c u m u l a n t s f o r fee. 13
F ig u re 3: T e m p e ratu re dependence of our calculated second c u m u la n t ơ { t ) of Cu (dash-
dot), Ni (dash) and of Cu doped by Ni atom (solid) compared to ex perim ent (dot) [13]
X 1Q'4
T(K)
Figure 4: T e m p e ra tu re dependence of our calculated th ird c u m u la n t cr '^ (r) of Cu (dash”
dot), Ni (dash) and of Cu doped by Ni atom (solid) compared to ex perim ent (dot) [13]
Trang 7T he t e m p e r a t u r e d e p e n d e n c e of o u r c a l c u l a t e d f i r s t c u m u l a n t or n e t t h e r m a l
e x p a n s io n ơ-(l)( ĩ ) ( F ig u re 2), second c u m u l a n t or D e b y e - W a l l e r fa c to r ơ 2( t ) (F igure
3) a n d t h i r d c u m u l a n t c r ^ ( r ) (F ig u re 4) s h o w s s i g n i f i c a n t c h a n g e s of t h e s e v a lu es
w h en Cu is d o p p ed by Ni a to m a n d a r e a s o n a b l e a g r e e m e n t b e t w e e n t h e c a lc u la te d
by th e p r e s e n t th e o r y a n d e x p e r i m e n t a l v a lu e s F i g u r e s 2, 3, 4 show t h a t th e
c u m u l a n t s of Cu becom e w e a k e r due to t h e d o p a n t by Ni a to m T h e s e i m p u 'r i t y effects a r e very i m p o r t a n t a n d th e y h a v e to be t a k e n in to a c c o u n t in th e e v a lu a tio n
of t h e r m o d y n a m ic p r o p e r t i e s of th e s u b s t a n c e s T h e c a l c u l a t e d f ir s t, seco n d a n d
t h i r d c u m u l a n t c o n t a i n i n g i m p u r i t y a to m a lso s a t i s f y a l l i m p o r t a n t p r o p e rtie s disco v ered in t h e o r y [7, 17] a n d e x p e r i m e n t [1] T h e y c o n t a i n z e ro -p o in t
c o n tr ib u tio n a t low t e m p e r a t u r e a n d a p p r o a c h t h e c l a s s i c a l t h e o r y r e s u l t s a t h ig h
t e m p e r a t u r e , i., e., t h e p r o p o r t i o n a l i t y to th e h ig h t e m p e r a t u r e is l i n e a r for th e firs t
a n d second c u m u l a n t a n d q u a d r a t i c a l for th e t h i r d c u m u l a n t
4 C o n c l u s i o n s
A new a n a l y t i c a l m e th o d for c a lc u l a t i o n a n d e v a l u a t i o n of th e r m o d y n a m ic
p r o p e r tie s of th e fee c r y s t a l s c o n t a i n i n g i m p u r i t y a t o m h a s b e e n d e v e lo p e d b ased on
th e q u a n t u m s t a t i s t i c a l th e o r y w ith c o r r e l a t e d E i n s t e i n m odel
O u r d e v e lo p m e n t is th e d e r iv a tio n of t h e a n a l y t i c a l e x p r e s s i o n s for t h e local effective force c o n s ta n ts , c o r r e la te d E i n s t e i n f r e q u e n c y a n d t e m p e r a t u r e , t h e first, second a n d t h i r d X AFS c u m u l a n t of fee c r y s t a l s c o n t a i n i n g i m p u r i t y ato m T h ey a re sig n ific a n tly d iff e r e n t from th o se of t h e p u r e m a t e r i a l s , b u t s a t i s f y all s t a n d a r d
p r o p e r tie s of t h e s e q u a n t i t i e s T h e s e d iff e r e n c e s d e n o t e t h e i m p u r i t y effects w hich discov ered in e x p e r i m e n t a n d th e y h a v e to be t a k e n in to a c c o u n t in th e e v a lu a tio n
of th e r m o d y n a m ic p r o p e r t i e s of th e s u b s t a n c e s
M orse p o t e n t i a l p a r a m e t e r s h a v e b e e n a ls o a n a l y t i c a l l y c a l c u la te d th u s
c o m p le tin g th e ab in itio c a l c u la tio n p r o c e d u r e of t h e c o n s id e r e d v a lu e s
T h e good a g r e e m e n t b e tw e e n t h e c a l c u l a t e d a n d t h e e x p e r i m e n t a l r e s u lts
d e m o n s t r a t e s th e efficiency a n d p o s s ib ility of u s i n g t h e p r e s e n t d e v e lo p e d p ro c e d u re
in XAFS d a t a a n a ly s is
A c k n o w l e d g m e n t s One of th e a u t h o r s (N V H u n g ) t h a n k s Dr I V Pirog for s e n d in g Ref 13 a n d p ro v id in g th e e x p e r i m e n t a l d a t a T h i s w o rk is s u p p o rte d in
p a r t by th e b a sic science r e s e a r c h n a t i o n a l p r o g r a m p r o v id e d by t h e M in is try of
S cience a n d T ech n o lo g y No 41.10.04
R e f e r e n c e s
1 E A S te rn , p L iv ins, z Z h a n g , P h y s R ev B, 4 3 ( 1 9 9 1 ) 8550
2 A.I F r a n k e l , J J R e h r, P hys R e v B, 4 8 (1 9 9 3 ) 585.
3 T M iy a n a g a , T F u jik a w a , J Phys Soc J p n , 6 3 ( 1 9 9 4 ) 1036 a n d 3683.
4 N V H u n g , R F r a h m , P h y s ic a B, 2 0 8 - 2 0 9 ( 1 9 9 5 ) 91.
Trang 8Ca lc u la ti o n o f XAFS c u r n u la n t s for f e e 15
5 N V H u n g , R F r a h m , H K a m its u b o , J Phys Soc J p n , 65(1996) 3571.
6 N V H u n g , J de P h y s iq u e , IV (1 9 9 7 ) C2 : 279.
7 N V H u n g , J J R e h r , P h y s Rev B, 56 (1 9 9 7) 43.
8 N V H u n g , C o m m u n P h y s., 8, No 1(1998) 46-54
9 N V H u n g , N B Due, R F r a h m , J P h y s Soc J p n , 72 (2 0 03 ) 1254.
10 M D a n ie l, D M P e a s e , e t al, a c c e p te d for p u b lic a tio n in P h y s Rev B.
1 1 N V H u n g , V N U - J o u r o f S cien c e, Vol 19, No 4(2003).
12 See X - r a y a b s o r p t i o n , e d i t e d by D c K o n in g s b e rg e r a n d R P r i n s (Wiley,
N ew Y o rk ,1 9 8 8 )
13 I V P iro g , T I N e d o s e i k i n a , I A Z a r u b i n a n d A T S h u v a e v , J Phys.:
C ondens M a t t e r , 1 4 (2 0 0 2 ) 1825.
14 L A G irifa lc o a n d V G W e iz e r, P h y s Rev 114(1959) 687.
15 N V H u n g , V N U - J o u r S c i e n c e , Vol 18, No 3, (2002), 17-23.
16 N V H u n g , D X V iet, V N U - J o u r S cien c e 19, N o 2(2003) 19.
17 J.M Z im a n , P r i n c i p l e s o f th e T h e o r y of S o lid s , 2 nd ed by C a m b rid g e
U n iv e r s ity P r e s s , 1972
18 R F F e y n m a n , S t a t i s t i c a l M e c h a n ic s (b e n ja m in , R e a d in g , 1972).