05.03.2014 Views

Calculation of reduced matrix elements - IPNL

Calculation of reduced matrix elements - IPNL

Calculation of reduced matrix elements - IPNL

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!