agn_obs_met - INAF-Osservatorio Astronomico di Roma

Metabolics: feeding the
monster
•Energetics
•Accretion
disks
•Emitted
spectrum
•Relativistic
discs
Fueling AGNs
Conversion of mass to energy with some efficiency 
The major problem for fueling an AGN is not the energy
requirement but angular momentum considerations.
Consider the angular momentum of a particle in the solar
circle around the Galactic centre:
Accretion: Basic concepts
Tidal disruption of a star
Accretion disk structure
• Rotating mass of gas in a cylindrical potential well
• Time-scale for processes to redistribute angular
momentum >> dynamical time-scales or radiative timescales  each gas element looses energy via collision and
radiative cooling but retains its angular momentum.
• Circular orbits
Viscosity present among gas anuli rotating with different
 tends to reduce the difference in velocity, dumping
out the shearing motion, and therefore tends to
dissipate energy as heat and radiation. Thus viscosity
converts gravitational potential energy into radiation in
an efficient manner.
Equation of motion of the gas.
Disk has surface density (R,t) and radial velocity vR(R,t).
Consider motion of an annulus of gas with inner radius R and extent R
Mass is: M = 2R R 
Angular momentum is L= MR2  = 2R R  R2
Equation of continuity
The variation of the mass of the annulus with the time must be equal
to the difference between the mass entering the radius R and the
mass going out.
Conservation of angular momentum
If viscosity is introduced among contiguous gas annuli the rate of change
of L will depend also in the viscous force:
If the rate of change of L would not depend on G then the disk
surface density time derivative would be zero and = ( R).
Conversely,  is a compicate function of R and t.
In general = (R,t, ). If = const. the previous eq. can be
integrated analitically.
If we assume an initial
configuration with all mass
confined at R0 at t0
(t0)=(R- R0), then while time
passes  broadens until all the
mass is distributed toward the
centre and all the angular
momentum is flown toward the
outer disk.
Dissipative processes act to
smooth the differences and
bring the mass toward the
inner disk.
Energy dissipation in the disk
Energy dissipation in the disk
During radial motion half of the liberated potential
energy goes into increasing the kinetic energy and half
goes into heat.
Viscous stresses transfer mechanical energy besides
momentum.
Hydrostatic equilibrium
along z-direction
If the motion in the disk in the z-direction is subsonic, then the
disk is in hydrostatic equilibrium. Gas and radiation pressure
gradient is balanced by the z component of the grav. attraction.
Assuming that the disk is homogeneus (=const along z).
Viscosity
Steady disks
Energy does not depend on viscosity! It is only the
gravitational potential that pays.
1/2
¥æ
æ æ R ö1/2 ö
·
æ R0 ö ö -2
2
0
L = 4p ò Q RdR = 3M ò çç1- ç ÷ ÷÷W RdR = 3M GM ò çç1- ç ÷ ÷÷R dR
èRø ø
èRø ø
R0
R0 è
R0 è
¥
·
+
¥
·
M GM
=
= accretion rate times the binding energy of the last stable orbit!
R0
independent on dissipative forces
The heat loss per unit area is given by the transfer equation:
Q - = acTc4 / t
a= Stefan's constant, t = r Hk( r, Tc )
if the disk is optically thick, each element radiates as a BB with TS Þ
Q- = s T 4
where s is the Stefan-Bolzman constant.
putting Q+ = Q -
Þ
·
é
æ æ R ö1/2 öù
3GM
M
ê
çç1- ç 0 ÷ ÷÷ú
Ts =
3
ê 8p R s è è R ø øú
ë
û
1/4
At each radius the spectrum is a BB with TS so the total spectrum is:
Sn µ
ò
Rout
R0
Bn (TS (R))2p RdR =
ò
Rout
R0
n3
2p RdR
exp(hn / kTS ) -1
· ö
æ
49
3GM M ÷
occurres at R= R0 Þ Tmax = 0.488 ç
ç 8p R 3s ÷
36
è
ø
1/4
Tmax
X-ray spectrum
•
•
Dissipate energy in optically thick disk – cool, no hard X-rays
MUST dissipate in optically thin material so that E >> kT (Compton)
Optically thin accretion flow –
low L/LEdd only!
Magnetic reconnection above
disk – no known alternatives at
high L/LEdd!
Inverse Compton scattering of lower energy photons by
energetic electrons in a corona surrounding the disk
Thomson diffusion
A monochromatic wave interacts with an
e-. The e- will be accelerated and will emit
radiation. The direction of the emitted
radiation will be in general different from
the direction of the incident radiation. If the
particle is not relativistic the frequency of
the emitted radiation will be the same of
the incident radiation.
Scattering cross section
Compton effect
E1 = hn1
p1 =
quantum mechanical particle approach
hn1 h
=
c
l1
E 2 = hn 2
p2 =
pe =
hn 2 h
=
c
l2
1
2
E 2 - me c 4 mec 2 = initial e - energy
c
E = final e - energy
Momentum conservation : p1 = p2 + pe Þ pe2 = p12 + p22 - 2 p1 · p2 = p12 + p22 - 2 p1 p2 cos q
Energy conservation : E 0 = mec 2 , E = E 02 + pe2c 2 Þ p1c + me c 2 = p2c + E 02 + pe2c 2
c( p1 - p2 ) + me c 2 = E 02 + pe2c 2 squaring : c 2 ( p1 - p2 ) 2 + 2cE 0 ( p1 - p2 ) = pe2c 2
pe2 = p12 + p22 - 2 p1 p2 + 2E 0 ( p1 - p2 ) /c
- p1 p2 + E 0 ( p1 - p2 ) /c = - p1 p2 cos q
p1 p2 (1- cosq ) = E 0 ( p1 - p2 ) /c
multiply each term by hc / p1 p2 E 0 :
hc
( p - p2 )h
(1- cos q ) = 1
now use l = h / p
E0
p1 p2
Compton effect
Inverse Compton
e1 = initial photon energy in the lab frame
e1 ' = initial photon energy in the e - frame
e2 = photon energy in the lab frame after the scattering
e2 ' = photon energy in the e - frame after the scattering
e1 '= e1g (1- b cos q ) from the lab frame to the e- frame
if e1'<< me c 2 one can apply the Thompson scattering, i.e. e2 '= e1 '
let's go back to the lab frame : e2 = e2 ' g (1+ b cosq ) = e1g 2 (1- b cosq )(1+ b cos q )
The photon energy has been incremented by a factor g 2 . This is for 1 photon and 1 e - .
Now we want to find the energy emitted per unit time by an isotropic distribution of
photons scattered by a isotropic distribution of e - .
N(e) = number of photons with energy e
U rad = N(e)e = N(hn )hn
Inverse Compton
Let's begin considering non - relativistic Thomson scattering. If the Poynting flux
2
c
(power per unit area) of a plane wave incident on the e is S =
E the E will
4p
accelerate the e -, and the accelerated e - will in turn emit radiation according to the
Larmour formula. The net result is to scatter a portion of the incoming radiation with
no net tranfer of energy. The scattered radiation had power :
P = S s T = cs T U rad
2
U rad = S /c = E /4 p
Let's consider now the radiation scattered by a relativistic e- . The Thomson formula
is valid in the primed frame if e'<< mec 2 . In this case :
P'= cs T U'rad
We need to go back in the lab frame. We know that P is a Lorentz invariant, therefore :
P = P'= cs T U'rad
and we only need to transform U'rad in U rad .
Inverse Compton
We know that e'= eg (1+ b cosq ). The rate at which successive photons arrive is multiplied by the
same factor, so N'(e') = N(e)g (1+ b cos q ). In the e - frame :
U' rad = N' hn '= Ng (1+ b cosq )hng (1+ b cos q ) = U rad g 2 (1+ b cos q ) 2
Thus the transformation between U' rad and U rad depends on the angle q between the direction of
the photons and the e - motion. The total energy density in the e- frame of a radiation field that is
isotropic in the observer frame is obtained by integrating over all directions:
U' rad = U rad g
2p
2
p
ò dj ò sinqdq (1+ b cosq )
0
2
0
where is the azimuthal angle around the x axis.
é 4g 2 1 2
ù
4(g 2 -1/4)
2
U' rad = U rad ê
- g (1- b )ú = U rad
3
3
ë 3
û
[ g 2 (1- b 2 ) = 1]
4
cs T U rad (g 2 -1/4) = Total power in the radiation field after IC upscattering.
3
The initial power of the photons was cs T U rad , so the net power added to the radiation field is :
4
4
4
PIC = cs T U rad (g 2 -1/4) - cs T U rad = cs T U rad (g 2 -1) = cs T U rad g 2b 2 [(g 2 -1) = g 2b 2 ]
3
3
3
P'= P = cs T U'rad =
Thermal Comptonization
éaverage fractional energyù émean number of ù
ê
ú´ê
ú º Compton y parameter
ëchange per scattering
û ëscatterings
û
ef 4 2
æ kT ö2
A = » g = 16ç 2 ÷ mean amplification per scattering
è mc ø
ei 3
assuming a thermal electron distribution : N(E) » E 2 exp(-E /kT)
after k scatterings : ek = ei A k
the probability of a photon undergoing k scatterings before escaping is:
pk (t es ) ~ eit esk
The intensity of the emergent radiation at ek is proportional to pk (t es ) :
-a
æ
ö
e
I(ek ) ~ I(ei )t esk ~ I(ei )ç k ÷
è ei ø
Þ a=
-log(t es )
log(A)
Thermal Comptonization
•  depends on 2 parameters: the optical
depth of the medium and the
temperature. Any spectral shape can be
produced with ad hoc choices of these
parameters.
• problem: in AGN  ~ 1 [0.5-1.5]
Two phases disc
(Haardt & Maraschi 90’)
Optically thick emission from the cool layer provides soft
photons input for Comptonization and the hard Comptonized
photons contribute to the heating of the thick phase. The
feedback between the two phases determines the fraction of
power emitted in three main components: a BB from the thick
phase, a power law from Comptonization in the hot layer and
a reflection component. The resulting spectrum is
~independent of the coronal parameters.
Two phase disc
• A fraction f of the gravitational power PG is
dissipated in the hot layer of optical dept <1,
while (1-f)PG is dissipated in the optically thick
phase.
• The total luminosity of the hot phase is
LT=ALS, where LS is the luminosity of the thick
phase. The luminosity added by the hot
phase is LC=(A-1)LS
• LC=LCU+LCD with LCD=LC ~0.5
• Photons directed downward are partly
absorbed and partly reflected:
– Lrfl=aLCD a~0.1-0.2
– Labs=(1-a)LCD will contribute to LS
Energy balance
LS = (1- f )PG + (1- a)LCD for phase 1
ALS = fPG + LS
for phase 2
solving for A and L S
LS = {1- f [1- (1- a)h]} PG
f
1- f [1- (1- a)h]
The outgoing luminosity is given by :
A = 1+
Lout = LS + LCU + Lrfl
LCU =
where :
(1- h) f
ah f
L S and Lrfl =
LS
1- f + fh - fah
1- f + fh - fah
For small f LCU and Lrfl are proportional to LS.
For f~1 LCU and Lrfl are determined by a and 
Energy balance
f
A = 1+
1- f [1- (1- a)h]
-log(t es )
a=
log(A)
æ kT ö 2
A = 16ç 2 ÷
è mc ø
The energy balance in the first
ew. Implies a relationship
between optical depth and
temperature, and therefore the
spectral shape of the
Comptonized component!
Emitted spectra
Spettro dei raggi X
Il disco di accrescimento che circonda il buco nero è una sorgente di radiazione UV
e X di bassa energia (soft X-rays). Lo spettro di emissione è di tipo termico
(radiazione di corpo nero). La “Comptonizzazione”(*) della componente soft X-rays
in una corona che circonda il disco di accrescimento è una possibile causa dello
spettro a legge di potenza per la componente di raggi X di energia 1-100 KeV (hard
X-rays)
La componente continua di hard Xrays incide sul disco di accrescimento
e produce uno spettro riflesso,
caratterizzato
dai
fenomeni
di
scattering
Compton
e
di
assorbimento fotoelettrico, seguito da
emissione di righe di flourescenza o
diseccitazione di tipo Auger.
Lo spettro di emissione dei metalli è
sovrapposto allo spettro continuo
(*) Si tratta del fenomeno di ICS (Inverse Compton Scattering) per cui un fotone aumenta la propria energia a
seguito di un processo di diffusione su un elettrone.
La riga di emissione Fe Kα
Quali sono le evidenze osservative del BH-Paradigma?
Quali informazioni possono essere estratte dalle misure astronomiche
sulle proprietà fisiche e geometriche del buco nero e del disco di accrescimento?
A causa dell’elevato valore dell’abbondanza cosmica del ferro, la riga Fe Kα è la
componente principale dello spettro di emissione.
Per assorbimento fotoelettrico uno dei due elettroni della
shell K (la shell più interna, con n=1) viene etratto
dall’atomo, lasciando una lacuna. Un elettrone della shell L
(n=2) occupa il “posto vacante”, rilasciando 6.4 KeV di
energia.
Lo studio del profilo di riga è un importante
strumento di diagnostica delle proprietà fisiche e
geometriche del buco nero e del disco di
accrescimento
[A. C. Fabian, G. Miniutti, astro-ph/0507409 v1 18 Jul 2005]
Profilo di riga
Numerosi effetti modificano il profilo della riga del ferro:
• effetti RELATIVISTICI
• effetti di ORIENTAZIONE
• posizione dell’ultima orbita stabile (ISCO)
• effetti di IONIZZAZIONE
• profilo di EMISSIVITÀ del disco
La simulazione dei diversi effetti permette di calcolare la forma di riga in
funzione dei parametri fisici e geometrici del sistema costituito dal buco nero
centrale e dal disco di accrescimento. Il confronto dei dati sperimentali con le
previsioni teoriche fornisce importanti informazioni sulla natura degli oggetti
astrofisici.
Effetti relativistici
EFFETTI RELATIVISTICI producono l’allargamento del profilo di riga e lo
spostamento verso il rosso (red-shift) del picco di emissione.
NEWTONIANO
RELATIVITA’
SPECIALE
EFFETTO DOPPLER
EFFETTO DOPPLER
TRASVERSO
BEAMING RELATIVISTICO
RELATIVITA’
GENERALE
REDSHIFT
GRAVITAZIONALE
PROFILO DI RIGA
CONVOLUZIONE DEI
DIFFERENTI EFFETTI
Broad lines from relativistic discs
Axisymmetric disc orbiting a
Schwarzschild BH. Disk
extend from r0 to r1 and it is
observed at inclination i.
The ratio of the emitted
energies of the photons from
a point in the disk to the
observed energies is given
by:
a
u
k
(a )
E em
em
1+ z =
=
E obs ( ua k a )
obs
Second order contri butions to Doppler effect
The expression of the Doppler effect in General Relativity is the
ratio between the observer frame and emitted frame product between
the quadri - velocity
and the quadri - wave vector.
(1+ z) =
' a
u
( ak )
E em n1
obs
=
=
E obs n 0 ( ua k a )
em
'
dx
dx
where ua' = a = (-1;0,0,0) and ua = a are the emitted and observed frame
ds
ds
w
quadri - velocities and k a = (- ; k ).
c
k is the direction of propagation of the wave with angular velocity
2p w
w = 2pn ; k =
=
where v p is the phase - velocity
l vp
ds is the metric. Let us consider the expression for ds provided by Roberson :
æ 2MG
ö 2 æ
M 2G 2
MG ö 2 2
2
2
2
2
2
ds = -c ç1+
2
b
g
+
....
dt
+
1+
2
g
dr
+
r
sin
q
d
j
+
d
q
(
)
ç
÷
÷
(
)
2
4 2
2
è
ø
rc
c
r
rc
è
ø
ds2 = -c 2Bdt 2 + Adr 2 + r 2 (sin 2 qdj 2 + dq 2 )
2MG
and A = 1+ 2g
2
rc
A
v
ds = cdt B 1- b 2 where b =
B
c
cdt
ua = (u 0 ;u i ) where u 0 =
=
ds
B
where B = 1-
MG
rc 2
and v is the three dimensional velocity.
1
1- b 2
A
B
w
and u i =
dx i
=
ds
recall that ua = (-1;0,0,0) and that ka = (- ; k ) then :
c
w
w
1
v
1
ua ka = ; ua ka = + i ·k
c
c
c
2 A
2 A
B 1- b
B 1- b
B
B
vi
c B 1- b 2
A
B
n1 (u a' ka' ) em
= a
=
n 0 (u ka ) obs - w
c
=
-
w
c
=
1
vi
1
+ ·k
A c
A
B 1- b 2
B 1- b 2
B
B
A
1
B =
=
1
v
c
(1- b cosq )
(1- i · k )
c
w
A
B 1- b 2
B
B 1- b 2
é
MG ù
1/ 2
1+ 2g 2 ú
æ 2GM ö ê
2
rc
1b
ç1÷
ê
2GM ú
è
rc 2 ø ê
ú
12
ë
rc û
=
(1- b cos q )
1/ 2
=
æ 2GM ö1/ 2 é
MG öæ 2GM öù
2æ
= ç11b
1+
2
g
1(1+ b cosq )
ç
2 ÷ ê
2 ֍
2 ÷ú
è
è
rc ø ë
rc øè
rc øû
1/ 2
1/ 2
ù
æ 2GM ö1/ 2 é
æ
ö
æ
ö
MG
2GM
2
1b
(1+ b cosq ) =
ç1÷
ç1+ 2g 2 ÷ç1ê
2
2 ÷ú
è
è
rc ø ë
rc øè
rc øû
v
1 v 2 GM
1+ cosq + .....
c
2 c 2 rc 2
leaving only the contributions up to second order in
1 v2
- 2 = transverse Doppler effect
2c
GM
= gravitational redshift
2
rc
r
2GM
if R =
; Rg = 2 and cosq = + /-1 than :
3Rg
c
GM
v
r = 1 ; and GM = 1
=
c
c
rc 2 6R
6R
n1
1
1
1
= 1+
n0
6R 12R 6R
v
c
Effetto di orientazione
La forma della riga dipende dall’angolo di inclinazione dell’asse del disco rispetto
alla linea di vista. Al crescere dell’angolo di inclinazione, l’effetto principale è
l’allargamento della riga che si estende verso le più alte energie.
Il contenuto di riga nel blu è una misura dell’angolo di inclinazione del disco.
i [deg]
(Metrica di Schwarzschild)
Posizione della ISCO
Schwarzschild vs Kerr
Il contenuto di riga nel rosso è una stima del raggio della ISCO (Innermost Stable
Circular Orbit) e permette di distinguire un buco nero di tipo Schwarzschild
(statico) da un buco nero di tipo Kerr (rotante).
Rg =
GM
c2
Schwarzschild  RISCO = 6Rg
Kerr  RISCO = 1.24Rg
parametro di spin
a/M
Effetti di ionizzazione
x (r ) =
4p × Fx (r )
n( r )
parametro di ionizazione
Fx (r )
n(r )
Flusso di raggi X incidente per unità di area alla distanza r
Densità di elettroni
Se la materia del disco è
fortemente ionizzata, aumenta
l ’ energia di soglia del
processo
di
assorbimento
fotoelettrico, e di conseguenza
si riduce l ’ efficienza per
l ’ emissione della riga di
fluorescenza.
Profilo di emissività del disco
Il profilo di emissività del disco definisce l’efficienza con cui la luce è emessa in
funzione della coordinata radiale del disco. Si assume una legge di potenza:
q = indice di emissività
e (r ) = r - q
emissività uniforme
emissività standard
emissività “steep”
Il caso “ steep ” implica
un ’ illuminazione
più
efficiente a piccoli raggi, cioè
nelle regioni più interne del
disco di accrescimento.
In questo caso, il profilo di
linea si allarga e si estende
verso il rosso: maggiore peso
è dato, infatti, alle zone
centrali del disco, dove
dominano gli effetti di redshift
gravitazionale.
MGC – 6 - 30 -15
Galassia di tipo Seyfert I
z = 0.00775
Osservazione della riga Fe Kα
CONCLUSIONI
• i = 33° ± 1°
XMM – Newton
Chandra
• rin = 1.8 ± 0.1 Rg (ISCO)
• a/M = 0.93 ± 0.01  KERR
• ξ < 30 erg cm s-1
• profilo di emissività
rout
rbr
rin
qout = 3.0 ± 0.1
standard
qin = 6.9 ± 0.6
steep
Dischi relativistici