simulations of agn jets on 100 kpc scale
play

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

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


  1. 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

  2. Why are Jets so mysterious? Jet powers are well constrained, but ... Salvatore Cielo - Q&Q - 2014-07-15 2

  3. Why are Jets so mysterious? Jet powers are well constrained, but ... ? Several unknowns Salvatore Cielo - Q&Q - 2014-07-15 2

  4. Why are Jets so mysterious? Jet powers are well constrained, but ... ? Duty cycles link to BH accretion Several unknowns Salvatore Cielo - Q&Q - 2014-07-15 2

  5. Why are Jets so mysterious? Jet powers are well constrained, but ... ? FRI / FR II Duty cycles intrinsic or environment? link to BH accretion Several unknowns Salvatore Cielo - Q&Q - 2014-07-15 2

  6. Why are Jets so mysterious? Jet powers are well constrained, but ... ? FRI / FR II Duty cycles intrinsic or environment? link to BH accretion Several unknowns Composition / EoS light/hard: Blandford, Znajek 1977 Blandford-Znajek VS Blandford-Payne Salvatore Cielo - Q&Q - 2014-07-15 2

  7. Why are Jets so mysterious? Jet powers are well constrained, but ... ? FRI / FR II Duty cycles intrinsic or environment? link to BH accretion Several unknowns Composition / EoS Observations: obscured, light/hard: Blandford, Znajek 1977 sparse, small scale gas. Blandford-Znajek Simulations: need to link VS several physical scales. Blandford-Payne Salvatore Cielo - Q&Q - 2014-07-15 2

  8. So, can jets from AGNs... IN GALAXIES, few Myr ...directly reduce SF? (positive VS negative FB) Salvatore Cielo - Q&Q - 2014-07-15 3

  9. So, can jets from AGNs... IN GALAXIES, few Myr ...directly reduce SF? (positive VS negative FB) V. Gaibler (+), 2014 (12), ... Silk + 2013, Dubois + 2013, ... Salvatore Cielo - Q&Q - 2014-07-15 3

  10. So, can jets from AGNs... IN GALAXIES, few Myr IN HALOS, >100 Myr ...directly reduce SF? ...prevent cool flows by heat? (positive VS negative FB) V. Gaibler (+), 2014 (12), ... Silk + 2013, Dubois + 2013, ... Salvatore Cielo - Q&Q - 2014-07-15 3

  11. So, can jets from AGNs... IN GALAXIES, few Myr IN HALOS, >100 Myr ...directly reduce SF? ...prevent cool flows by heat? Perseus cluster (positive VS negative FB) ~150 kpc CHANDRA nasa.gov Complex heating V. Gaibler (+), 2014 (12), ... mechanisms (bubbles, etc.) Silk + 2013, invoked with little testing Dubois + 2013, ... Salvatore Cielo - Q&Q - 2014-07-15 3

  12. My setup, in the FLASH code • Dark Matter NFW potential (static, spherical, external) • Hot halo Isothermal, self-grav + DM hydrostatic equilibrium, rad. cooling (e.g. M halo = 2e12 M sun => T ~ 5e6 K) Salvatore Cielo - Q&Q - 2014-07-15 4

  13. My setup, in the FLASH code • Dark Matter NFW potential (static, spherical, external) FRII: recollimation • Hot halo shocks Isothermal, self-grav + DM hydrostatic equilibrium, rad. cooling (e.g. M halo = 2e12 M sun => T ~ 5e6 K) • Jet density contrast 1:100 P jet = ρ jet A v jet3 = 1e42-1e46 erg/s FRI/FRII: Mach 5 - Mach 100 on/off cycles: 1-100 Myr Box: 60-640kpc Resolution:10-150 pc Salvatore Cielo - Q&Q - 2014-07-15 4

  14. My setup, in the FLASH code • Dark Matter NFW potential (static, spherical, external) FRII: recollimation • Hot halo shocks Isothermal, self-grav + DM hydrostatic equilibrium, rad. cooling (e.g. M halo = 2e12 M sun => T ~ 5e6 K) • Jet density contrast 1:100 P jet = ρ jet A v jet3 = 1e42-1e46 erg/s FRI/FRII: Mach 5 - Mach 100 on/off cycles: 1-100 Myr Box: 60-640kpc • Warm gas clouds (1e4 K) Resolution:10-150 pc in a central galaxy + SF with sink particles. One day. Salvatore Cielo - Q&Q - 2014-07-15 4

  15. 10 42 erg/s 120 kpc Results: Evolutionary stages Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 5

  16. 1. Cocoon 10 42 erg/s and HS 120 kpc formation Results: Evolutionary stages Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 5

  17. 1. Cocoon 10 42 erg/s and HS 120 kpc formation 2. Forward propagation Results: Evolutionary stages Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 5

  18. 1. Cocoon 10 42 erg/s and HS 120 kpc formation 2. Forward propagation 3. Lobes inflation Results: Evolutionary stages Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 5

  19. 1. Cocoon 10 42 erg/s and HS 120 kpc formation 2. Forward propagation 3. Lobes inflation 4. Bubbles (jet off) Results: Evolutionary stages Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 5

  20. Bubbles still there after .6 Gyr Salvatore Cielo - Q&Q - 2014-07-15 6

  21. Results: Energy and Volume Coverage “p dV” Mechanical Work as fraction of P jet 1 1E45 erg/s All gas 1E44 erg/s Coupling: (W pdV / t age ) / P jet 0.8 1E43 erg/s 1E42 erg/s 0.6 0.4 Lobes/Bubbles only 0.2 0 0 2 10 10 t age / Myr Cielo +, in prep. Salvatore Cielo - Q&Q - 2014-07-15 7

  22. Results: Energy and Volume Coverage “p dV” Mechanical Work Volume coverage within as fraction of P jet r Cool ~30kpc (2e12 M Sun ) 1 1 1E45 erg/s 1E45 erg/s Low � res All gas 1E44 erg/s 1E44 erg/s Coupling: (W pdV / t age ) / P jet 0.8 1E43 erg/s 0.8 1E43 erg/s 1E42 erg/s 30 kpc / V sphere 30 kpc 1E42 erg/s 0.6 0.6 0.4 0.4 Lobes/Bubbles only V C 0.2 0.2 0 0 0 2 1 2 10 10 10 10 t age / Myr t age / Myr Cielo +, in prep. Salvatore Cielo - Q&Q - 2014-07-15 7

  23. Results: Energy and Volume Coverage “p dV” Mechanical Work Volume coverage within as fraction of P jet r Cool ~30kpc (2e12 M Sun ) 1 1 1E45 erg/s 1E45 erg/s Low � res All gas 1E44 erg/s 1E44 erg/s Coupling: (W pdV / t age ) / P jet 0.8 1E43 erg/s 0.8 1E43 erg/s 1E42 erg/s 30 kpc / V sphere 30 kpc 1E42 erg/s 0.6 0.6 FRII can make it... FRI not large enough! 0.4 0.4 Lobes/Bubbles only V C 0.2 0.2 0 0 0 2 1 2 10 10 10 10 t age / Myr t age / Myr Cielo +, in prep. Salvatore Cielo - Q&Q - 2014-07-15 7

  24. Results: angle flip 10 45 erg/s ~150 kpc Can the volume fraction increase with multiple injections? Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 8

  25. Results: angle flip 10 45 erg/s ~150 kpc Can the volume fraction Not really... :( increase with multiple injections? Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 8

  26. Summary • Can we simulate AGNs radio mode with AMR? • self-consistently? • Do jets have enough energy? • Do they affect enough volume in the ISM? Salvatore Cielo - Q&Q - 2014-07-15 9

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend


More recommend