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 are well constrained, but ... Salvatore Cielo - Q&Q - 2014-07-15 2
Why are Jets so mysterious? Jet powers are well constrained, but ... ? Several unknowns Salvatore Cielo - Q&Q - 2014-07-15 2
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
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
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
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
So, can jets from AGNs... IN GALAXIES, few Myr ...directly reduce SF? (positive VS negative FB) Salvatore Cielo - Q&Q - 2014-07-15 3
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
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
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
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
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
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
10 42 erg/s 120 kpc Results: Evolutionary stages Salvatore Cielo - Q&Q - 2014-07-15 Salvatore Cielo - Q&Q - 2014-07-15 5
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
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
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
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
Bubbles still there after .6 Gyr Salvatore Cielo - Q&Q - 2014-07-15 6
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
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
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
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
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
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
Recommend
More recommend
