Calculation of reduced matrix elements - IPNL
Calculation of reduced matrix elements - IPNL
Calculation of reduced matrix elements - IPNL
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
•〈l a ‖Y (L) ‖l b 〉= ?<br />
The spherical harmonics Y (L)<br />
M<br />
(ˆr) are the coordinate representation <strong>of</strong> the 2L+1 spherical<br />
components <strong>of</strong> a tensor <strong>of</strong> rank L, the Wigner-Eckart theorem can be written (for example<br />
for M= 0)<br />
( )<br />
〈l a m a |Y (L)<br />
0 |l bm b 〉=(−1) l a−m a la L l b<br />
〈l<br />
-m a 0 m a ‖Y (L) ‖l b 〉. (6)<br />
b<br />
Furthermore using〈l b m b |ˆr〉=Y (l a)∗<br />
m a<br />
(ˆr) and〈ˆr|l a m a 〉=Y (l b)<br />
m b<br />
(ˆr) one has<br />
∫<br />
〈l a m a |Y (L)<br />
0 |l bm b 〉= dˆr Y (l a)∗<br />
m a<br />
(ˆr)Y (L)<br />
0 (ˆr)Y(l b)<br />
m b<br />
(ˆr)<br />
where L andl b can be coupled tolusing<br />
√<br />
∑<br />
Y (L)<br />
(2L+1)(2lb + 1)(2l+1)<br />
0 (ˆr)Y(l b)<br />
m b<br />
(ˆr)=<br />
4π<br />
so one has<br />
∑<br />
〈l a m a |Y (L)<br />
0 |l bm b 〉=<br />
lm<br />
lm<br />
√<br />
(2L+1)(2lb + 1)(2l+1)<br />
4π<br />
The integral gives<br />
∫ ∫<br />
dˆr Y (l a)∗<br />
m a<br />
(ˆr)Y m<br />
(l)∗ (ˆr)=(−1)m<br />
so<br />
〈l a m a |Y (L)<br />
0 |l bm b 〉=(−1) −m a<br />
Using the properties <strong>of</strong> the 3− j symbols<br />
( )( )<br />
L la l b L la l b<br />
=<br />
0 -m a m b 0 0 0<br />
so<br />
〈l a m a |Y (L)<br />
0 |l bm b 〉=(−1) −m a<br />
( )( )<br />
L l lb L l lb<br />
0 m m b 0 0 0<br />
( )( )∫<br />
L l lb L l lb<br />
0 m m b 0 0 0<br />
dˆr Y (l a)∗<br />
m a<br />
(ˆr)Y -m (l)∗ (ˆr)=(−1)m δ la lδ ma -m<br />
√<br />
(2L+1)(2la + 1)(2l b + 1)<br />
4π<br />
√<br />
(2L+1)(2la + 1)(2l b + 1)<br />
4π<br />
Y m<br />
(l)∗ (ˆr)<br />
dˆr Y (l a)∗<br />
m a<br />
( )( )<br />
L la l b L la l b<br />
.<br />
0 -m a m b 0 0 0<br />
( )( )<br />
la L l b la L l b<br />
δ<br />
-m a 0 m b 0 0 0<br />
ma m b<br />
Using eq. (6) and eq. (7) we obtain<br />
√<br />
〈l a ‖Y (L) ‖l b 〉=(−1) l (2L+1)(2la + 1)(2l<br />
a<br />
b + 1)<br />
4π<br />
(ˆr)Y (l)∗ (ˆr).<br />
( )( )<br />
la L l b la L l b<br />
δ<br />
-m a 0 m b 0 0 0<br />
ma m b<br />
.<br />
(7)<br />
( )<br />
la L l b<br />
.<br />
0 0 0<br />
m<br />
3