Simulations of AGN jets on 100 kpc scale Salvatore Cielo, MPIA - - PowerPoint PPT Presentation

▶
simulations of agn jets on 100 kpc scale
SMART_READER_LITE
LIVE PREVIEW

Simulations of AGN jets on 100 kpc scale Salvatore Cielo, MPIA - - PowerPoint PPT Presentation

Apr 02, 2024 15 likes •279 views

Simulations of AGN jets on 100 kpc scale Salvatore Cielo, MPIA Heidelberg, DE Andrea V. Macci, MPIA Heidelberg, DE Vincenzo Antonuccio-Delogu, INAF, Catania, IT Salvatore Cielo - Q&Q - 2014-07-15 Why are Jets so mysterious? Jet powers

slide-1
SLIDE 1

Salvatore Cielo - Q&Q - 2014-07-15

Simulations of AGN jets

  • n 100 kpc scale

Salvatore Cielo, MPIA Heidelberg, DE Andrea

  • V. Macciò, MPIA Heidelberg, DE

Vincenzo Antonuccio-Delogu, INAF, Catania, IT

slide-2
SLIDE 2

Salvatore Cielo - Q&Q - 2014-07-15

Why are Jets so mysterious?

2

Jet powers are well constrained, but ...

slide-3
SLIDE 3

Salvatore Cielo - Q&Q - 2014-07-15

Why are Jets so mysterious?

2

?

Several unknowns

Jet powers are well constrained, but ...

slide-4
SLIDE 4

Salvatore Cielo - Q&Q - 2014-07-15

Why are Jets so mysterious?

2

?

Several unknowns

Duty cycles link to BH accretion Jet powers are well constrained, but ...

slide-5
SLIDE 5

Salvatore Cielo - Q&Q - 2014-07-15

Why are Jets so mysterious?

2

?

Several unknowns

Duty cycles link to BH accretion FRI / FR II intrinsic or environment? Jet powers are well constrained, but ...

slide-6
SLIDE 6

Salvatore Cielo - Q&Q - 2014-07-15

Why are Jets so mysterious?

2

?

Several unknowns

Duty cycles link to BH accretion FRI / FR II intrinsic or environment? Jet powers are well constrained, but ... Composition / EoS light/hard: Blandford-Znajek VS Blandford-Payne

Blandford, Znajek 1977

slide-7
SLIDE 7

Salvatore Cielo - Q&Q - 2014-07-15

Why are Jets so mysterious?

2

?

Several unknowns

Duty cycles link to BH accretion FRI / FR II intrinsic or environment? Jet powers are well constrained, but ... Observations: obscured, sparse, small scale gas. Simulations: need to link several physical scales. Composition / EoS light/hard: Blandford-Znajek VS Blandford-Payne

Blandford, Znajek 1977

slide-8
SLIDE 8

Salvatore Cielo - Q&Q - 2014-07-15

So, can jets from AGNs...

3

IN GALAXIES, few Myr ...directly reduce SF? (positive VS negative FB)

slide-9
SLIDE 9

Salvatore Cielo - Q&Q - 2014-07-15

So, can jets from AGNs...

3

IN GALAXIES, few Myr ...directly reduce SF? (positive VS negative FB)

  • V. Gaibler (+), 2014 (12), ...

Silk + 2013, Dubois + 2013, ...

slide-10
SLIDE 10

Salvatore Cielo - Q&Q - 2014-07-15

So, can jets from AGNs...

3

IN GALAXIES, few Myr ...directly reduce SF? IN HALOS, >100 Myr ...prevent cool flows by heat? (positive VS negative FB)

  • V. Gaibler (+), 2014 (12), ...

Silk + 2013, Dubois + 2013, ...

slide-11
SLIDE 11

Salvatore Cielo - Q&Q - 2014-07-15

So, can jets from AGNs...

3

IN GALAXIES, few Myr ...directly reduce SF? IN HALOS, >100 Myr ...prevent cool flows by heat? (positive VS negative FB)

  • V. Gaibler (+), 2014 (12), ...

Silk + 2013, Dubois + 2013, ... Perseus cluster ~150 kpc CHANDRA nasa.gov

Complex heating mechanisms (bubbles, etc.) invoked with little testing

slide-12
SLIDE 12

Salvatore Cielo - Q&Q - 2014-07-15

  • Dark Matter

NFW potential (static, spherical, external)

  • Hot halo

Isothermal, self-grav + DM hydrostatic equilibrium, rad. cooling (e.g. Mhalo = 2e12 Msun => T ~ 5e6 K)

My setup, in the FLASH code

4

slide-13
SLIDE 13

Salvatore Cielo - Q&Q - 2014-07-15

  • Dark Matter

NFW potential (static, spherical, external)

  • Hot halo

Isothermal, self-grav + DM hydrostatic equilibrium, rad. cooling (e.g. Mhalo = 2e12 Msun => T ~ 5e6 K)

My setup, in the FLASH code

4

  • Jet

density contrast 1:100 Pjet = ρjet A vjet3 = 1e42-1e46 erg/s FRI/FRII: Mach 5 - Mach 100

  • n/off cycles: 1-100 Myr

Box: 60-640kpc Resolution:10-150 pc

FRII: recollimation shocks

slide-14
SLIDE 14

Salvatore Cielo - Q&Q - 2014-07-15

  • Dark Matter

NFW potential (static, spherical, external)

  • Hot halo

Isothermal, self-grav + DM hydrostatic equilibrium, rad. cooling (e.g. Mhalo = 2e12 Msun => T ~ 5e6 K)

My setup, in the FLASH code

4

  • Warm gas clouds (1e4 K)

in a central galaxy + SF with sink particles. One day.

  • Jet

density contrast 1:100 Pjet = ρjet A vjet3 = 1e42-1e46 erg/s FRI/FRII: Mach 5 - Mach 100

  • n/off cycles: 1-100 Myr

Box: 60-640kpc Resolution:10-150 pc

FRII: recollimation shocks

slide-15
SLIDE 15

Salvatore Cielo - Q&Q - 2014-07-15 5

Results: Evolutionary stages

Salvatore Cielo - Q&Q - 2014-07-15

1042 erg/s 120 kpc

slide-16
SLIDE 16

Salvatore Cielo - Q&Q - 2014-07-15 5

Results: Evolutionary stages

  • 1. Cocoon

and HS formation

Salvatore Cielo - Q&Q - 2014-07-15

1042 erg/s 120 kpc

slide-17
SLIDE 17

Salvatore Cielo - Q&Q - 2014-07-15 5

Results: Evolutionary stages

  • 1. Cocoon

and HS formation

  • 2. Forward

propagation

Salvatore Cielo - Q&Q - 2014-07-15

1042 erg/s 120 kpc

slide-18
SLIDE 18

Salvatore Cielo - Q&Q - 2014-07-15 5

Results: Evolutionary stages

  • 1. Cocoon

and HS formation

  • 2. Forward

propagation

  • 3. Lobes

inflation

Salvatore Cielo - Q&Q - 2014-07-15

1042 erg/s 120 kpc

slide-19
SLIDE 19

Salvatore Cielo - Q&Q - 2014-07-15 5

Results: Evolutionary stages

  • 1. Cocoon

and HS formation

  • 2. Forward

propagation

  • 3. Lobes

inflation

  • 4. Bubbles

(jet off)

Salvatore Cielo - Q&Q - 2014-07-15

1042 erg/s 120 kpc

slide-20
SLIDE 20

Salvatore Cielo - Q&Q - 2014-07-15

Bubbles still there after .6 Gyr

6

slide-21
SLIDE 21

Salvatore Cielo - Q&Q - 2014-07-15

Results: Energy and Volume Coverage

7

“p dV” Mechanical Work as fraction of Pjet

10 10

2

0.2 0.4 0.6 0.8 1 tage / Myr Coupling: (WpdV / tage) / Pjet

1E45 erg/s 1E44 erg/s 1E43 erg/s 1E42 erg/s

All gas Lobes/Bubbles only

Cielo +, in prep.

slide-22
SLIDE 22

Salvatore Cielo - Q&Q - 2014-07-15

Results: Energy and Volume Coverage

7

“p dV” Mechanical Work as fraction of Pjet

10 10

2

0.2 0.4 0.6 0.8 1 tage / Myr Coupling: (WpdV / tage) / Pjet

1E45 erg/s 1E44 erg/s 1E43 erg/s 1E42 erg/s

All gas Lobes/Bubbles only

10

1

10

2

0.2 0.4 0.6 0.8 1 tage / Myr VC

30 kpc / Vsphere 30 kpc

Lowres

1E45 erg/s 1E44 erg/s 1E43 erg/s 1E42 erg/s

Volume coverage within rCool ~30kpc (2e12 MSun) Cielo +, in prep.

slide-23
SLIDE 23

Salvatore Cielo - Q&Q - 2014-07-15

Results: Energy and Volume Coverage

7

“p dV” Mechanical Work as fraction of Pjet

10 10

2

0.2 0.4 0.6 0.8 1 tage / Myr Coupling: (WpdV / tage) / Pjet

1E45 erg/s 1E44 erg/s 1E43 erg/s 1E42 erg/s

All gas Lobes/Bubbles only

10

1

10

2

0.2 0.4 0.6 0.8 1 tage / Myr VC

30 kpc / Vsphere 30 kpc

Lowres

1E45 erg/s 1E44 erg/s 1E43 erg/s 1E42 erg/s

Volume coverage within rCool ~30kpc (2e12 MSun) FRII can make it... FRI not large enough! Cielo +, in prep.

slide-24
SLIDE 24

Salvatore Cielo - Q&Q - 2014-07-15

Results: angle flip

8

Can the volume fraction increase with multiple injections?

Salvatore Cielo - Q&Q - 2014-07-15

1045 erg/s ~150 kpc

slide-25
SLIDE 25

Salvatore Cielo - Q&Q - 2014-07-15

Results: angle flip

8

Can the volume fraction increase with multiple injections?

Not really... :(

Salvatore Cielo - Q&Q - 2014-07-15

1045 erg/s ~150 kpc

slide-26
SLIDE 26

Salvatore Cielo - Q&Q - 2014-07-15

  • Can we simulate AGNs radio mode with AMR?
  • self-consistently?
  • Do jets have enough energy?
  • Do they affect enough volume in the ISM?

Summary

9