Magnetic and Antimagnetic rotation in Covariant Density Functional Theory Jie Meng 孟 杰 北京航空航天大学物理与核能工程学院 School of Physics and Nuclear Energy Engineering Beihang University (BUAA) 北京大学物理学院 School of Physics Peking University (PKU)
Magnetic and Antimagnetic rotation in Covariant Density Functional Theory
Jie Meng 孟 杰
北京航空航天大学物理与核能工程学院 School of Physics and Nuclear Energy Engineering
Beihang University (BUAA)
北京大学物理学院 School of Physics Peking University (PKU)
Outline
2
2011-9-21
Introduction
Theoretical framework
Magnetic rotation
Antimagnetic rotation
Summary & Perspectives
Angular Momentum World of Nucleus
3
2011-9-21
Electric Rotation
4
2011-9-21
Substantial quadrupole deformation
Strong electric quadrupole (E2) transitions
Coherent collective rotation of many nucleons Bohr PR1951
Twin PRL1986
Electric & Magnetic Rotation 2011-9-21
∆I = 2
E2 Transitions
∆I = 1
M1 Transitions
Hübel PPNP2005
Twin PRL1986
Magnetic rotation
6
2011-9-21
near spherical or weakly deformed nuclei
strong M1 and very weak E2 transitions
rotational bands with ∆I = 1
shears mechanism
∆I=1 regular bands
First attempt in verify MR
ΔI = 1 Enhanced magnetic dipole transition
How does B(M1) change with spin I ?
Good agreement between TAC and PRM
Good agreement with prediction for BM1 versus I
Experiment: MR
11
2011-9-21
Experiment: MR
Magnetic rotation: 78 nuclei
A~60
Antimagnetic Rotation (AMR) Magnetic rotation Ferromagnet rotational bands with ΔI = 1
near spherical nuclei; weak E2 transitions
strong M1 transitions
B(M1) decrease with spin
shears mechanism
Antimagnetic rotation Antiferromagnet rotational bands with ΔI = 2
near spherical nuclei; weak E2 transitions
no M1 transitions
B(E2) decrease with spin
two “shears-like” mechanism
Experiment: AMR Antimagnetic rotation: 3 nuclei
Other mass regions
Simons PRL2003; Simons PRC2005
Small B(E2) Decrease tendency
Large J(2)/B(E2) Increase tendency
Theory
Semiclassical particle plus rotor model Clark ARNPS2000
simple geometry for the energies and transition probabilities
2
2 2( )( ) ( )
2I j jE I V Pπ ν θ− −
= +ℑ
jπ
jν
I
R
Core Rotor Particle shears mechanism
Theory
Semiclassical particle plus rotor model Clark ARNPS2000
simple geometry for the energies and transition probabilities
Pairing-plus-quadrupole tilted axis cranking (TAC) model Frauendorf NPA1993; Frauendorf NPA2000
semi-phenomenological Hamiltonian
H H Jω′ = − ⋅
2
22sphH H Q Q GP P Nµ µµ
χ λ+ +
=−
= − − −∑
Theory
Semiclassical particle plus rotor model Clark ARNPS2000
simple geometry for the energies and transition probabilities
Pairing-plus-quadrupole tilted axis cranking (TAC) model Frauendorf NPA1993; Frauendorf NPA2000
semi-phenomenological Hamiltonian
Phenomenological investigations
polarization effects are neglected or only partially considered
nuclear currents are treated without self-consistency
adjusted to show MR/AMR in some way or another
A fully self-consistent microscopic investigation?
DFT: Cranking version
TAC based on Covariant Density Functional Theory Meson exchange version:
3-D Cranking: Madokoro, Meng, Matsuzaki, Yamaji, PRC 62, 061301 (2000)
2-D Cranking: Peng, Meng, Ring, Zhang, PRC 78, 024313 (2008)
Point coupling version: Simple and more suitable for systematic investigations
2-D Cranking: Zhao, Zhang, Peng, Liang, Ring, Meng, PLB 699, 181 (2011)
TAC based on Skyrme Density Functional Theory
Fully self-consistent microscopic investigations
fully taken into account polarization effects
self-consistently treated the nuclear currents
without any adjustable parameters for rotational excitations
3-D Cranking: Olbratowski, Dobaczewski, Dudek, Płóciennik, PRL 93, 052501(2004)
2-D Cranking: Olbratowski, Dobaczewski, Dudek, Rzaca-Urban, Marcinkowska, Lieder, APPB 33, 389(2002)
Outline
19
2011-9-21
Introduction
Theoretical framework
Magnetic rotation
Antimagnetic rotation
Summary & Perspectives
ˆ
δδEh =
ρ iiih ϕεϕ =ˆMean field: Eigenfunctions:
ˆ
2
δδδ EV =ρ ρ
Interaction:
Skyrme Gogny
Extensions: Pairing correlations, Covariance Relativistic Hartree Bogoliubov (RHB) theory
Density functional theory in nuclear physics
Starting point of CDFT
21
2011-9-21
Nucleons are coupled by exchange of mesons via an effective Lagrangian with all relativistic symmetries, used in a mean field concept and no-sea approximation
σ ω ρ
π
meson Jπ T π 0− 1 σ 0+ 0 ω 1− 0 ρ 1− 1
22
Brief introduction of CDFT 2011-9-21
3152
2 2 2
2 2
[ ( ) ]
1 1 1 1 1 2 2 4 2 41 1 1 1 2 2 2 4
fmL i M g g g e A
m m R R
m m F F
π
π
τµ µ µµ σ ω µ ρ µ µ µ
µ µν µ µνµ σ µν ω µ µν
µ µ µνρ µ µ π µν
ψ γ σ γ ω τ ρ γ γ π τ ψ
σ σ σ ω ω
ρ ρ π π π π
−= ∂ − − − + • + − ∂ •
+ ∂ ∂ − − Ω Ω + − •
+ + ∂ • ∂ − • −
RF A A
µν µ ν ν µ
µν µ ν ν µ
µν µ ν ν µ
ω ω
ρ ρ
Ω = ∂ − ∂
= ∂ − ∂
= ∂ − ∂
4
, , , ,
1( ) ( ) ( ) ( , ) ( ) ( )2
i ii A
H i M d y x y D x y y x
T Vσ ω ρ π
ψ ψ ψ ψ ψ ψ=
= − • ∇ + + Γ
= +
∑∫γ
1 2
1 2 5 1 5 2
2
em 3 1 3 2
(1, 2) (1) (2), (1, 2) ( ) ( ) ,
(1, 2) ( ) ( ) , (1, 2) ( ) ( )
(1,2) ( (1 )) ( (1 ))4
f fm m
g g g g
g g
e
π π
π π
µσ σ σ ρ ρ µ ρ
µ νω ω µ ω µ π µ ν
µµ
γ τ γ τ
γ γ τγ γ τγ γ
γ τ γ τ
Γ ≡ − Γ ≡ +
Γ ≡ + Γ ≡ − ∂ ∂
Γ ≡ + − −
Hamiltonian:
Lagrangian:
Brief introduction of CDFT
23
2011-9-21
'
'† ††
†
†
( )( ) [ ( ) ]
( ) ([ ( ) ])
i
i
i
i
i
i
i ti i
i
i t
i ti i
i tii
i
x f e c
x
g e d
ef e dc g
εε
ε ε
ψ
ψ −
−= +
= +
∑
∑
xx
xx
†0 0cα
α
Φ = ∏
0 0 0 0 0 0, , , ,
emem
ii A
D D D Dk
E E E E
E H T V
E E E E EE E E E Eσ ω ρ π
σ ω ρσ ω ρ π
=
= Φ Φ = Φ Φ + Φ Φ
= + + + + + + + + +
∑
Energy density functional:
, , , ,i
i AH T V
σ ω ρ π=
= + ∑
†
††1 2 ' ' ' '
; ' '
( ) ,
1 (1) (2) (1,2) (1,2) (2) (1)2i i i
T d f i M f c c
V d d c c c c f f D f f
α β α βαβ
α β β α α β β ααβ α β
γ= − ⋅∇ +
= Γ
∑∫
∑∫
x
x xHartree
Fock
Equations of motion in RMF theory
24
2011-9-21
Same footing for
Deformation Rotation Pairing (RHB,BCS,SLAP)
…
For system with time invariance:
( )( ) ( ) i i iV M Sα β ψ ε ψ⋅ + + + = p r r
33
1( ) ( ) ( ) ( )2
( ) ( )
V g g e A
S r g
ω ρ
σ
τω τ ρ
σ
− = + + =
r r r r
r
( ) ( )
2 2 32 3
2 33
2 33
s
b
n pb b
m g g g
m g c
m g d
σ σ
ω ω
ρ ρ
σ ρ σ σ
ω ρ ω
ρ ρ ρ ρ
−∆ + = − − − −∆ + = −
−∆ + = − −
1
1
3 31
31
( ) ( ) ( )
( ) ( ) ( )
( ) ( ) ( )
1( ) ( ) ( )2
i
i
i
As i ii
Av ii
Aii
Ac ii
ρ ψ ψ
ρ ψ ψ
ρ ψ τ ψ
τρ ψ ψ
=
+=
+=
+=
= = = − =
∑∑∑
∑
r r r
r r r
r r r
r r r
RMF theory with Point-Coupling interaction
25
2011-9-21
( )
2 2 2 3 42 3
0 30
14
1 1 1((
12
) )2 3 4
12
ii i i
v
s
H i M F F
m g g g
gg
νν
σ σ
ρω
ψ γ ψ
σ σ σρ σ σ
ω ρ ρρ
= − • ∇ + +
+ ∇ + + +
++
+
ω
( )
2 3 4
22
14
1 1 1 12 2 3 4
1 12 2
1 12 2
ii i i
S S S S S S
V
S S S
TV V TV TV TV V V VV
H i M F Fνν
α ρ
ψ γ ψ
α ρ δ ρ ρ β ρ γ ρ
α ρ δ ρ ρδ ρ ρ
= − ∇ + +
+ + ∆ + +
+ +∆ ∆++
ω
2 2
2
2
2 2 4
11 / v v v vv v v
ggm m
g gm m
ω
ωω ω
ω
ω
ωωω ρ ρ ρ ρα ρδ= = + ∆ + ≈ + ∆
− ∆
Equations of motion in RMF-PC theory
26
2011-9-21
For system with time invariance:
( )( ) ( ) i i iV M Sα β ψ ε ψ⋅ + + + = p r r
3 33 3
2 3
1( ) ( ) ( ) ( ) ( ) ( ) ( )2
( )
V V V V V V TV TV TV TV
S S S S S S S S
V e A
S r
τα ρ γ ρ δ ρ τ α ρ τ δ ρ
α ρ β ρ γ ρ δ ρ
− = + + ∆ + + ∆ + = + + + ∆
r r r r r r r
1
1
3 31
31
( ) ( ) ( )
( ) ( ) ( )
( ) ( ) ( )
1( ) ( ) ( )2
i
i
i
As i ii
Av ii
Aii
Ac ii
ρ ψ ψ
ρ ψ ψ
ρ ψ τ ψ
τρ ψ ψ
=
+=
+=
+=
= = = − =
∑∑∑
∑
r r r
r r r
r r r
r r r
Without Klein-Gordon equation
Tilted axis cranking CDFT
27
2011-9-21
3 4 2
3
( )1 1 1( )( ) ( )( ) ( )( )2 2 21 1 1( ) ( ) [( )( )]3 4 41 1 1( ) ( ) ( ) ( ) ( ) ( )2 2 21 1
2 4
S V TV
S S V
S V TV
L i m
e A F F
µµ
µ µµ µ
µµ
ν ν µ νν ν µ ν µ µ
µ µνµ µν
ψ γ ψ
α ψψ ψψ α ψγ ψ ψγ ψ α ψτγ ψ ψτγ ψ
β ψψ γ ψψ γ ψγ ψ ψγ ψ
δ ψψ ψψ δ ψγ ψ ψγ ψ δ ψτγ ψ ψτγ ψ
τ ψγ ψ
= ∂ −
− − −
− − −
− ∂ ∂ − ∂ ∂ − ∂ ∂
−− −
General Lagrangian density
Transformed to the frame rotating with the uniform velocity
Koepf NPA1989; Kaneko PLB1993; Madokoro PRC1997
( ,0, ) ( cos ,0, sin )x z θ θΩ ΩΩ = Ω Ω = Ω Ω
1 00
t ttx xα α
= → = = x R xx
TAC RMF:equations of motion
28
2011-9-21
Dirac Equation
Potential
Spatial components of vector field are involved due to the time-reversal invariance broken
( )( ) ( )( ) ( ) i i ii M S VV Jα β ψ ε ψΩ⋅ ⋅ − ∇ − + + + − =
rr r
2 3
3 33 3
( )1( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
2
S S S S S S S S
V V V V V V TV TV TV TV
S r
V j j j j j e Aµ µ µ µ µ µ µ
α ρ β ρ γ ρ δ ρτα γ δ τ α τ δ
= + + + ∆ −
= + + ∆ + + ∆ +r r r r r r r
Observables 2011-9-21
3 2tot
1
3 4 2
00
1
1 1 1( ) ( )2 2 2
2 3 3 1 1( ( ) ) ( )3 4 4 2 21 1( ) | |2 2
A
k S S V V V TV TV TVk
S S S S V V V S S S V V V
A
TV TV TV pk
E d r j j j j
j j j j
j j ej A k J k
µ µµ µ
µ µµ µ
µµ
α ρ α α
β ρ γ ρ γ δ ρ ρ δ
δ
=
=
= − + +
+ + + + ∆ + ∆
+ ∆ + + ⟨ Ω ⟩
∑ ∫
∑
Binding energy
Angular momentum 2ˆ ( 1)J I I⟨ ⟩ = +
Quadrupole moments and magnetic moments
2 2 2 220 22
5 153 ,16 32
Q z r Q x yπ π
= ⟨ − ⟩ = ⟨ − ⟩
23 ?†[ ( ) ( ) ( ) ( )]
A
i i i i i ii
mcd r Q r r r r rc
µ ψ αψ κψ β ψ= × + Σ∑∫
Where 1, 0, 1.793, 1.913p n p nQ Q κ κ= = = = −
Observables 2011-9-21
3 2tot
1
3 4 2
00
1
1 1 1( ) ( )2 2 2
2 3 3 1 1( ( ) ) ( )3 4 4 2 21 1( ) | |2 2
A
k S S V V V TV TV TVk
S S S S V V V S S S V V V
A
TV TV TV pk
E d r j j j j
j j j j
j j ej A k J k
µ µµ µ
µ µµ µ
µµ
α ρ α α
β ρ γ ρ γ δ ρ ρ δ
δ
=
=
= − + +
+ + + + ∆ + ∆
+ ∆ + + ⟨ Ω ⟩
∑ ∫
∑
Binding energy
Angular momentum 2ˆ ( 1)J I I⟨ ⟩ = +
Quadrupole moments and magnetic moments
2 23 3( 1) ( sin cos )8 8 x J z JB M µ µ θ µ θπ π⊥= = −
B(M1) and B(E2) transition probabilites
2
2 220 22
3 2( 2) cos (1 sin )8 3J JB E Q Qθ θ
= + +
RMF parameterizations
31
2011-9-21
Meson Exchange
2
NL3, NLSH, TM1, TM2, PK1, …
Density dependent parameterizations:
2 3 3 3, , , , , , , , , ,M m m m g g g g g c dσ ω ρ σ ω ρ
Point Coupling
, , , , , , , , ,S V TV S V TV S S VM α α α δ δ δ β γ γ
PC-LA, PC-F1, PC-PK1 …
Density dependent parameterizations:
Nonlinear parameterizations: Nonlinear parameterizations:
TW99, DD-ME1, DD-ME2, PKDD, … DD-PC1, …
, , , , ( ), ( ), ( )M m m m g g gσ ω ρ σ ω ρρ ρ ρ , , ( ), ( ), ( )S S V TVM δ α ρ α ρ α ρ
Parameterizations: PC-PK1
32
2011-9-21
Zhao, Li, Yao, Meng, PRC 82, 054319 (2010)
Spherical nuclei
33
2011-9-21
Zhao, Li, Yao, Meng, PRC 82, 054319 (2010)
Deformed nuclei
34
2011-9-21
Zhao, Li, Yao, Meng, PRC 82, 054319 (2010)
Nuclear Mass
35
2011-9-21
Exp value for 2149 nuclei from Audi et al. NPA2003
Outline
36
2011-9-21
Introduction
Theoretical framework
Magnetic rotation
Antimagnetic rotation
Summary & Perspectives
MR: 60Ni lightest nucleus with magnetic rotation
Torres PRC 2008
MR: 60Ni
Harmonic oscillator shells: Nf = 10
Parameter set: PC-PK1
Configurations:
Zhao, Zhang, Peng, Liang, Ring, Meng, PLB 699, 181 (2011)
Numerical Details
(Py=PRy(π))
P PxT Px PzT Pz PyT Py
Jx √ × √ √ × √ ×
Jz √ √ × × √ √ ×
Identification of the energy levels
Symmetry
|nx ny nz ns⟩ → |nljmz⟩ → |nljmx⟩
Spherical basis Cartesian basis
Constrained intrinsic framework
( )22 1 2 1
12
H H C Q a− −′ = + −
Enforce principal axes to be identical with the x, y, and z axis.
2 1 0a − =
Numerical Details
Parallel transport principle
(δ) ( ) 1j iφ φΩ + Ω Ω ≈
Ω : Φi(Ω) Ω+δΩ : Φj(Ω+δΩ)
MR: 60Ni
Zhao, Zhang, Peng, Liang, Ring, Meng, PLB 699, 181 (2011)
MR: 60Ni
Zhao, Zhang, Peng, Liang, Ring, Meng, PLB 699, 181 (2011)
MR: 60Ni
Shears mechanism
Magnetic Rotation
Electric Rotation
Zhao, Zhang, Peng, Liang, Ring, Meng, PLB 699, 181 (2011)
MR: 60Ni
Electromagnetic transition properties Zhao, Zhang, Peng, Liang, Ring, Meng, PLB 699, 181 (2011)
Outline
45
2011-9-21
Introduction
Theoretical framework
Magnetic rotation
Antimagnetic rotation
Summary & Perspectives
AMR: 105Cd first odd-A nucleus with antimagnetic rotation Choudhury PRC2010
AMR: 105Cd Harmonic oscillator shells: Nf = 10
Parameter set: PC-PK1
Configurations:
Polarizations:
107Sn 107Sn + π[(g9/2)-2] 105Cd Zhao, Peng, Liang, Ring, Meng PRL 177, 122501(2011)
Choudhury PRC2010
2011-9-21
AMR: Single particle routhians
Neutron Proton
• Time reversal symmetry broken energy splitting
• For proton, two holes in the top of g9/2 shell
• For neutron, one particle in the bottom of h11/2 shell, the other six are
distributed over the (gd) shell with strong mixing
• This configuration is similar to , but not exactly
AMR: 105Cd
Zhao, Peng, Liang, Ring, Meng PRL 177, 122501(2011)
AMR: 105Cd
Polarization effects play important roles in the description of AMR, especially
for E2 transitions.
Zhao, Peng, Liang, Ring, Meng PRL 177, 122501(2011)
AMR: 105Cd
Two “shears-like” mechanism
Zhao, Peng, Liang, Ring, Meng PRL 177, 122501(2011)
AMR: 105Cd
In the microscopic point of view, increasing angular momentum results from
the alignment of the proton holes and the mixing within the neutron orbitals.
Zhao, Peng, Liang, Ring, Meng PRL 177, 122501(2011)
Outline
53
2011-9-21
Introduction
Theoretical framework
Magnetic rotation
Antimagnetic rotation
Summary & Perspectives
Summary and Perspectives 2011-9-21
Covariant density functional theory has been extended to describe rotational excitations including MR and AMR.
PC-PK1: fitted with BE & Rc for 60 nuclei provide better description for the isospin dependence of BE 60Ni: MR reproduce well the E, I, B(M1), and B(E2) values, inclduing the transition from electric rotation to magnetic rotation
105Cd: AMR reproduce well the AMR pictures, E, I, and B(E2) values in a fully self-consistent microscopic way for the first time and find that polarization effects play important roles
Summary and Perspectives 2011-9-21
? Pairing effects
? Transition from the electric to magnetic rotation
? Anti-magnetic rotation in other mass regions
? High-K isomer
Thank You!