The Physics and Cosmology of TeV Blazars

Similar documents
The Physics and Cosmology of TeV Blazars

The Cosmological Impact of Blazars: from Plasma Instabilities to Structure Formation

The Plasma Physics and Cosmological Impact of TeV Blazars

Are TeV blazars heating the IGM?

Galaxies 626. Lecture 3: From the CMBR to the first star

Really, really, what universe do we live in?

Cosmic ray feedback in hydrodynamical simulations. simulations of galaxy and structure formation

Signal Model vs. Observed γ-ray Sky

Constraints on Extragalactic Background Light from Cherenkov telescopes: status and perspectives for the next 5 years

Cosmology and fundamental physics with extragalactic TeV γ-rays?

Cosmic Rays in Galaxy Clusters: Simulations and Perspectives

Planck meets the Lyman-α forest

Galaxy Formation and Evolution

High Energy Astrophysics

Hunting for Dark Matter in Anisotropies of Gamma-ray Sky: Theory and First Observational Results from Fermi-LAT

The Secondary Universe

EBL Studies with the Fermi Gamma-ray Space Telescope

Cosmic ray feedback in hydrodynamical simulations. simulations of galaxy and structure formation

Outline. Walls, Filaments, Voids. Cosmic epochs. Jeans length I. Jeans length II. Cosmology AS7009, 2008 Lecture 10. λ =

GRB history. Discovered 1967 Vela satellites. classified! Published 1973! Ruderman 1974 Texas: More theories than bursts!

ACTIVE GALACTIC NUCLEI: FROM THE CENTRAL BLACK HOLE TO THE GALACTIC ENVIRONMENT

Formation and growth of galaxies in the young Universe: progress & challenges

Fermi: Highlights of GeV Gamma-ray Astronomy

Cosmic Ray Astronomy. Qingling Ni

Constraining Dark Matter annihilation with the Fermi-LAT isotropic gamma-ray background

Probing the End of Dark Ages with High-redshift Quasars. Xiaohui Fan University of Arizona Dec 14, 2004

A New View of the High-Energy γ-ray Sky with the Fermi Telescope

Searching for Dark Matter in the Galactic Center with Fermi LAT: Challenges

Lecture 20 High-Energy Astronomy. HEA intro X-ray astrophysics a very brief run through. Swift & GRBs 6.4 kev Fe line and the Kerr metric

The Intergalactic Medium: Overview and Selected Aspects

Active Galactic Nuclei

Broadband X-ray emission from radio-quiet Active Galactic Nuclei

Accretion Disks. 1. Accretion Efficiency. 2. Eddington Luminosity. 3. Bondi-Hoyle Accretion. 4. Temperature profile and spectrum of accretion disk

AGN in hierarchical galaxy formation models

Observations of Active Galactic Nuclei at very high energies with H.E.S.S.

TEMA 6. Continuum Emission

Observational Cosmology

Can blazar flares be triggered by the VHE gamma-rays from the surrounding of a supermassive black hole?

Galaxies and Cosmology

The Black Hole in the Galactic Center. Eliot Quataert (UC Berkeley)

AST Cosmology and extragalactic astronomy. Lecture 20. Black Holes Part II

Dark Matter ASTR 2120 Sarazin. Bullet Cluster of Galaxies - Dark Matter Lab

Introduction to AGN. General Characteristics History Components of AGN The AGN Zoo

Cherenkov Telescope Array Status Report. Salvatore Mangano (CIEMAT) On behalf of the CTA consortium

Where are the missing baryons? Craig Hogan SLAC Summer Institute 2007

Cosmic Rays, Photons and Neutrinos

Rupert Croft. QuickTime and a decompressor are needed to see this picture.

The Inner Region of the Milky Way Galaxy in High Energy Gamma Rays

Active Galactic Nuclei - Zoology

Cosmological Evolution of Blazars

Radiation processes and mechanisms in astrophysics I. R Subrahmanyan Notes on ATA lectures at UWA, Perth 18 May 2009

Cherenkov Telescope Array ELINA LINDFORS, TUORLA OBSERVATORY ON BEHALF OF CTA CONSORTIUM, TAUP

Orianne ROOS CEA-Saclay Collaborators : F. Bournaud, J. Gabor, S. Juneau

Lower bound on the extragalactic magnetic field

Hubble Space Telescope ultraviolet spectroscopy of blazars: emission lines properties and black hole masses. E. Pian, R. Falomo, A.

Cross-Correlation of Cosmic Shear and Extragalactic Gamma-ray Background

Reionization. High-Redshift Galaxy UV Luminosity Function. Axion Dark Matter. Rosemary Wyse

Astr 2320 Thurs. April 27, 2017 Today s Topics. Chapter 21: Active Galaxies and Quasars

Lecture 27 The Intergalactic Medium

Set 4: Active Galaxies

The visible constituents of the Universe: Non-relativistic particles ( baryons ): Relativistic particles: 1. radiation 2.

CTA SKA Synergies. Stefan Wagner Landessternwarte (CTA Project Office) Heidelberg

Feedback and Galaxy Formation

Exploring the Ends of the Rainbow: Cosmic Rays in Star-Forming Galaxies

Galaxy formation and evolution II. The physics of galaxy formation

TeV gamma-rays from UHECR sources 22 radio log10(e /ev ) 16 photon horizon γγ e + e CMB 14 IR kpc 10kpc 100kpc M pc Virgo 10M pc 100M pc G

2. Active Galaxies. 2.1 Taxonomy 2.2 The mass of the central engine 2.3 Models of AGNs 2.4 Quasars as cosmological probes.

Non-Blazar Gamma-ray Active Galactic Nuclei seen by Fermi-LAT. C.C. Teddy Cheung Naval Research Lab/NRC on behalf of the Fermi-LAT Collaboration

Gamma-ray echo (pair echo) from TeV sources and extragalactic magnetic fields

High Energy Emission. Brenda Dingus, LANL HAWC

Lecture 9. Quasars, Active Galaxies and AGN

PERSPECTIVES of HIGH ENERGY NEUTRINO ASTRONOMY. Paolo Lipari Vulcano 27 may 2006

SIMILARITY AND DIVERSITY OF BLACK HOLE SYSTEMS View from the Very High Energies

Galaxies 626. Lecture 5

Cosmic Ray Electrons and GC Observations with H.E.S.S.

The Inner Region of the Milky Way Galaxy in High Energy Gamma Rays

Radio Observations of TeV and GeV emitting Supernova Remnants

Quasars ASTR 2120 Sarazin. Quintuple Gravitational Lens Quasar

Jun. Prof. Dr. Alessandro MIRIZZI (Universität Hamburg)

80 2 Observational Cosmology L and the mean energy

Low-Energy Cosmic Rays

v Characteristics v Possible Interpretations L X = erg s -1

Probing dark matter and the physical state of the IGM with the Lyα forest

AGN and Radio Galaxy Studies with LOFAR and SKA

Diagnostics of Leptonic vs. Hadronic Emission Models for Blazars Prospects for H.E.S.S.-II and CTA Markus Böttcher North-West University

Princeton December 2009 The fine-scale structure of dark matter halos

Emmanuel Moulin! on behalf of the CTA Consortium!!! Rencontres de Moriond 2013! Very High Energy Phenomena in the Universe! March 9-16, La Thuile,

Testing a DM explanation of the positron excess with the Inverse Compton scattering

Gamma-ray Astrophysics

GAMMA-RAYS FROM MASSIVE BINARIES

Discovery of TeV Gamma-ray Emission Towards Supernova Remnant SNR G Last Updated Tuesday, 30 July :01

Active Galactic Alexander David M Nuclei

Evidence for BH: Active Galaxies

Age-redshift relation. The time since the big bang depends on the cosmological parameters.

High Redshift Universe

WFIRST and JWST synergies in the study of First Light. M. Stiavelli STScI, Baltimore

Modelling the Variability of Blazar Jets

1. GAMMA-RAY BURSTS & 2. FAST RADIO BURSTS

Explosive reconnection of the double tearing mode in relativistic plasmas

The High-Energy Interstellar Medium

Transcription:

1 in collaboration with Avery E. Broderick 2, Phil Chang 3, Ewald Puchwein 1, Volker Springel 1 1 Heidelberg Institute for Theoretical Studies, Germany 2 Perimeter Institute/University of Waterloo, Canada 3 University of Wisconsin-Milwaukee, USA Jul 11, 2012 / Meeting, DFG Research Unit FOR 1254

The Hitchhiker s Guide to... Blazar Heating the extragalactic TeV Universe plasma physics for cosmologists

The Hitchhiker s Guide to... Blazar Heating the extragalactic TeV Universe plasma physics for cosmologists consequences for intergalactic magnetic fields extragalactic gamma-ray background thermal history of the Universe Lyman-α forest missing dwarf galaxies H I mass function galaxy cluster bimodality

Outline Physics of blazar heating TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 1 Physics of blazar heating TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 2 Properties of blazar heating Thermal history of the IGM The Lyman-α forest 3 Formation of dwarf galaxies Puzzles in galaxy formation Conclusions

TeV gamma-ray astronomy TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background H.E.S.S. MAGIC I VERITAS MAGIC II

The TeV gamma-ray sky TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background There are several classes of TeV sources: Galactic - pulsars, BH binaries, supernova remnants Extragalactic - mostly blazars, two starburst galaxies

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Unified model of active galactic nuclei narrow line region relativistic jet broad line region dusty torus central SMBH

The blazar sequence TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background continuous sequence from LBL IBL HBL TeV blazars are dim (very sub-eddington) TeV blazars have rising spectra in the Fermi band (α < 2) define TeV blazar = hard IBL + HBL Ghisellini (2011), arxiv:1104.0006

Propagation of TeV photons TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 1 TeV photons can pair produce with 1 ev EBL photons: γ TeV + γ ev e + + e

Propagation of TeV photons TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 1 TeV photons can pair produce with 1 ev EBL photons: γ TeV + γ ev e + + e mean free path for this depends on the density of 1 ev photons: λ γγ (35... 700) Mpc for z = 1... 0 pairs produced with energy of 0.5 TeV (γ = 10 6 )

Propagation of TeV photons TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 1 TeV photons can pair produce with 1 ev EBL photons: γ TeV + γ ev e + + e mean free path for this depends on the density of 1 ev photons: λ γγ (35... 700) Mpc for z = 1... 0 pairs produced with energy of 0.5 TeV (γ = 10 6 ) these pairs inverse Compton scatter off the CMB photons: mean free path is λ IC λ γγ /1000 producing gamma-rays of 1 GeV E γ 2 E CMB 1 GeV

Propagation of TeV photons TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 1 TeV photons can pair produce with 1 ev EBL photons: γ TeV + γ ev e + + e mean free path for this depends on the density of 1 ev photons: λ γγ (35... 700) Mpc for z = 1... 0 pairs produced with energy of 0.5 TeV (γ = 10 6 ) these pairs inverse Compton scatter off the CMB photons: mean free path is λ IC λ γγ /1000 producing gamma-rays of 1 GeV E γ 2 E CMB 1 GeV each TeV point source should also be a GeV point source

What about the cascade emission? TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Every TeV source should be associated with a 1-100 GeV gamma-ray halo not seen! TeV spectra expected cascade emission TeV detections Fermi constraints Neronov & Vovk (2010)

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Measuring IGM B-fields from TeV/GeV observations TeV beam of e + /e are deflected out of the line of sight reducing the GeV IC flux lower limit on B Larmor radius r L = E ( ) ( ) 1 E eb 30 B 3 TeV 10 16 Mpc G

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Measuring IGM B-fields from TeV/GeV observations TeV beam of e + /e are deflected out of the line of sight reducing the GeV IC flux lower limit on B Larmor radius r L = E ( ) ( ) 1 E eb 30 B 3 TeV 10 16 Mpc G IC mean free path x IC 0.1 ( ) 1 E Mpc 3 TeV for the associated 10 GeV IC photons the Fermi angular resolution is 0.2 or θ 3 10 3 rad x IC r L > θ B 10 16 G

Missing plasma physics? TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background How do beams of e + /e propagate through the IGM? plasma processes are important interpenetrating beams of charged particles are unstable consider the two-stream instability: one frequency (timescale) and one length in the problem: ω p γ = 4πe 2 n e γ 2, λ p = γc ( γ ) 10 14 cm m e ω p 10 ρ(z=0) 6

Two-stream instability: mechanism TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background wave-like perturbation with k v beam, longitudinal charge oscillations in background plasma (Langmuir wave): initially homogeneous beam-e : attractive (repulsive) force by potential maxima (minima) e attain lowest velocity in potential minima bunching up e + attain lowest velocity in potential maxima bunching up Φ p p e e

Two-stream instability: mechanism TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background wave-like perturbation with k v beam, longitudinal charge oscillations in background plasma (Langmuir wave): beam-e + /e couple in phase with the background perturbation: enhances background potential stronger forces on beam-e + /e positive feedback exponential wave-growth instability Φ e + e + e p e p e e

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Two-stream instability: energy transfer particles with v v phase : pair energy plasma waves growing modes particles with v v phase : plasma wave energy pairs damped modes

Oblique instability Physics of blazar heating TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background k oblique to v beam : real word perturbations don t choose easy alignment = all orientations Bret (2009), Bret+ (2010)

Beam physics growth rates TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background excluded for collective plasma phenomena oblique consider a light beam penetrating into relatively dense plasma maximum growth rate 0.4 γ n beam n IGM ω p IC oblique instability beats IC by two orders of magnitude Broderick, Chang, C.P. (2012)

Beam physics complications... TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background non-linear saturation: non-linear evolution of these instabilities at these density contrasts is not known expectation from PIC simulations suggest substantial isotropization of the beam assume that they grow at linear rate up to saturation

Beam physics complications... TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background non-linear saturation: non-linear evolution of these instabilities at these density contrasts is not known expectation from PIC simulations suggest substantial isotropization of the beam assume that they grow at linear rate up to saturation plasma instabilities dissipate the beam s energy, no (little) energy left over for inverse Compton scattering off the CMB

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background TeV emission from blazars a new paradigm γ TeV + γ ev e + + e IC off CMB γ GeV plasma instabilities heating IGM

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background TeV emission from blazars a new paradigm γ TeV + γ ev e + + e IC off CMB γ GeV plasma instabilities heating IGM absence of γ GeV s has significant implications for... intergalactic B-field estimates γ-ray emission from blazars: spectra, background

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background TeV emission from blazars a new paradigm γ TeV + γ ev e + + e IC off CMB γ GeV plasma instabilities heating IGM absence of γ GeV s has significant implications for... intergalactic B-field estimates γ-ray emission from blazars: spectra, background additional IGM heating has significant implications for... thermal history of the IGM: Lyman-α forest late time structure formation: dwarfs, galaxy clusters

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Implications for B-field measurements Fraction of the pair energy lost to inverse-compton on the CMB: f IC = Γ IC /(Γ IC + Γ oblique ) Broderick, Chang, C.P. (2012)

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Conclusions on B-field constraints from blazar spectra it is thought that TeV blazar spectra might constrain IGM B-fields this assumes that cooling mechanism is IC off the CMB + deflection from magnetic fields beam instabilities may allow high-energy e + /e pairs to self scatter and/or lose energy isotropizes the beam no need for B-field 1 10% of beam energy to IC CMB photons

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Conclusions on B-field constraints from blazar spectra it is thought that TeV blazar spectra might constrain IGM B-fields this assumes that cooling mechanism is IC off the CMB + deflection from magnetic fields beam instabilities may allow high-energy e + /e pairs to self scatter and/or lose energy isotropizes the beam no need for B-field 1 10% of beam energy to IC CMB photons TeV blazar spectra are not suitable to measure IGM B-fields!

TeV blazar luminosity density: today TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background collect luminosity of all 23 TeV blazars with good spectral measurements account for the selection effects (sky coverage, duty cycle, galactic occultation, TeV flux limit) TeV blazar luminosity density is a scaled version (η B 0.2%) of that of quasars! Broderick, Chang, C.P. (2012)

Unified TeV blazar-quasar model TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Quasars and TeV blazars are: regulated by the same mechanism contemporaneous elements of a single AGN population: TeV-blazar activity does not lag quasar activity Broderick, Chang, C.P. (2012)

Unified TeV blazar-quasar model TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Quasars and TeV blazars are: regulated by the same mechanism contemporaneous elements of a single AGN population: TeV-blazar activity does not lag quasar activity assume that they trace each other for all redshifts! Broderick, Chang, C.P. (2012)

How many TeV blazars are there? TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Hopkins+ (2007)

How many TeV blazars are there? TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Hopkins+ (2007)

How many TeV blazars are there? TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Hopkins+ (2007)

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Fermi number count of TeV blazars η = 0.8, 1.6, 3.1 number evolution of TeV blazars that are expected to have been observed by Fermi vs. observed evolution colors: different flux (luminosity) limits connecting the Fermi and the TeV band: L TeV,min (z) = η L Fermi,min (z) Broderick, Chang, C.P. (2012)

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Fermi number count of TeV blazars η = 0.8, 1.6, 3.1 number evolution of TeV blazars that are expected to have been observed by Fermi vs. observed evolution colors: different flux (luminosity) limits connecting the Fermi and the TeV band: L TeV,min (z) = η L Fermi,min (z) Broderick, Chang, C.P. (2012) evolving (increasing) blazar population consistent with observed declining evolution (Fermi flux limit)!

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background How many TeV blazars are there at high-z? Hopkins+ (2007)

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Fermi probes dragons of the gamma-ray sky Fermi LAT Extragalactic Gamma-ray Background Intensity (GeV photons per cm 2 per sec per steradian) 10-6 10-7 Unknown contributors Background accounted for by unresolved AGN 10-8 0.1 1 10 Energy (GeV) 100

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Extragalactic gamma-ray background assume all TeV blazars have identical intrinsic spectra: E b is break energy, F E = LˆF E 1 (E/E b ) α L 1 + (E/E b ) α 1, α L < α are low and high-energy spectral indexes extragalactic gamma-ray background (EGRB): E 2 dn 1 (E, z) = de 4π z dv (z ) η B Λ Q (z )ˆF E 4πD 2 L E = E(1 + z ) is gamma-ray energy at emission, Λ Q is physical quasar luminosity density, η B 0.2% is blazar fraction, τ is optical depth e τ E (E,z ),

TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background Extragalactic gamma-ray background dotted: unabsorbed EGRB due to TeV blazars dashed: absorbed EGRB due to TeV blazars solid: absorbed EGRB, after subtracting the resolved TeV blazars (z < 0.25) Broderick, Chang, C.P. (2012)

Outline Physics of blazar heating Properties of blazar heating Thermal history of the IGM The Lyman-α forest 1 Physics of blazar heating TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 2 Properties of blazar heating Thermal history of the IGM The Lyman-α forest 3 Formation of dwarf galaxies Puzzles in galaxy formation Conclusions

Evolution of the heating rates HI,HeI /HeII reionization Properties of blazar heating Thermal history of the IGM The Lyman-α forest Heating Rates [ev Gyr 1 ] 10 3 10 2 10 1 0.1 2 10 photoheating photoheating standard blazar standard blazar fit optimistic blazar optimistic blazar fit blazar heating 10x larger heating 10 5 2 1 1 +z Chang, Broderick, C.P. (2012)

Blazar heating vs. photoheating Properties of blazar heating Thermal history of the IGM The Lyman-α forest total power from AGN/stars vastly exceeds the TeV power of blazars

Blazar heating vs. photoheating Properties of blazar heating Thermal history of the IGM The Lyman-α forest total power from AGN/stars vastly exceeds the TeV power of blazars T IGM 10 4 K (1 ev) at mean density (z 2) ε th = kt m pc 2 10 9

Blazar heating vs. photoheating Properties of blazar heating Thermal history of the IGM The Lyman-α forest total power from AGN/stars vastly exceeds the TeV power of blazars T IGM 10 4 K (1 ev) at mean density (z 2) ε th = kt m pc 2 10 9 radiative energy ratio emitted by BHs in the Universe (Fukugita & Peebles 2004) ε rad = η Ω bh 0.1 10 4 10 5

Blazar heating vs. photoheating Properties of blazar heating Thermal history of the IGM The Lyman-α forest total power from AGN/stars vastly exceeds the TeV power of blazars T IGM 10 4 K (1 ev) at mean density (z 2) ε th = kt m pc 2 10 9 radiative energy ratio emitted by BHs in the Universe (Fukugita & Peebles 2004) ε rad = η Ω bh 0.1 10 4 10 5 fraction of the energy energetic enough to ionize H I is 0.1: ε UV 0.1ε rad 10 6 kt kev

Blazar heating vs. photoheating Properties of blazar heating Thermal history of the IGM The Lyman-α forest total power from AGN/stars vastly exceeds the TeV power of blazars T IGM 10 4 K (1 ev) at mean density (z 2) ε th = kt m pc 2 10 9 radiative energy ratio emitted by BHs in the Universe (Fukugita & Peebles 2004) ε rad = η Ω bh 0.1 10 4 10 5 fraction of the energy energetic enough to ionize H I is 0.1: ε UV 0.1ε rad 10 6 kt kev photoheating efficiency η ph 10 3 kt η ph ε UV m pc 2 ev (limited by the abundance of H I/He II due to the small recombination rate)

Blazar heating vs. photoheating Properties of blazar heating Thermal history of the IGM The Lyman-α forest total power from AGN/stars vastly exceeds the TeV power of blazars T IGM 10 4 K (1 ev) at mean density (z 2) ε th = kt m pc 2 10 9 radiative energy ratio emitted by BHs in the Universe (Fukugita & Peebles 2004) ε rad = η Ω bh 0.1 10 4 10 5 fraction of the energy energetic enough to ionize H I is 0.1: ε UV 0.1ε rad 10 6 kt kev photoheating efficiency η ph 10 3 kt η ph ε UV m pc 2 ev (limited by the abundance of H I/He II due to the small recombination rate) blazar heating efficiency η bh 10 3 kt η bh ε rad m pc 2 10 ev (limited by the total power of TeV sources)

Thermal history of the IGM Properties of blazar heating Thermal history of the IGM The Lyman-α forest 10 5 only photoheating standard BLF optimistic BLF T [K] 10 4 10 3 20 10 5 2 1 1 +z Chang, Broderick, C.P. (2012)

Properties of blazar heating Thermal history of the IGM The Lyman-α forest Evolution of the temperature-density relation no blazar heating

Properties of blazar heating Thermal history of the IGM The Lyman-α forest Evolution of the temperature-density relation no blazar heating blazars and extragalactic background light are uniform: blazar heating rate independent of density

Properties of blazar heating Thermal history of the IGM The Lyman-α forest Evolution of the temperature-density relation no blazar heating blazars and extragalactic background light are uniform: blazar heating rate independent of density makes low density regions hot causes inverted temperature-density relation, T 1/δ

Properties of blazar heating Thermal history of the IGM The Lyman-α forest Evolution of the temperature-density relation no blazar heating with blazar heating Chang, Broderick, C.P. (2012) blazars and extragalactic background light are uniform: blazar heating rate independent of density makes low density regions hot causes inverted temperature-density relation, T 1/δ

Blazars cause hot voids Properties of blazar heating Thermal history of the IGM The Lyman-α forest no blazar heating with blazar heating Chang, Broderick, C.P. (2012)

Blazars cause hot voids Properties of blazar heating Thermal history of the IGM The Lyman-α forest no blazar heating with blazar heating Chang, Broderick, C.P. (2012) blazars completely change the thermal history of the diffuse IGM and late-time structure formation

Simulations with blazar heating Properties of blazar heating Thermal history of the IGM The Lyman-α forest Puchwein, C.P., Springel, Broderick, Chang (2012): L = 15h 1 Mpc boxes with 2 384 3 particles one reference run without blazar heating three with blazar heating at different levels of efficiency (address uncertainty) used an up-to-date model of the UV background (Faucher-Giguère+ 2009)

-2-2 Physics of blazar heating Temperature-density relation Properties of blazar heating Thermal history of the IGM The Lyman-α forest 7.0 no blazar heating intermediate blazar heating 6.5 6.0 log 10 ( N HI N H ) = 8-7 -6-8 -7-6 log 10 (T /K) 5.5 5.0 4.5-5 -4-5 -4 4.0-3 -3 3.5 Viel et al. 2009, F=0.1-0.8 log 10 (Mpix/(h 1 M )) Viel et al. 2009, F=0-0.9 5 6 7 8 9 10 3.0 2 1 0 1 2 3 2 1 0 1 2 3 log 10 (ρ/ ρ ) Puchwein, C.P., Springel, Broderick, Chang (2012)

Ly-α spectra Physics of blazar heating Properties of blazar heating Thermal history of the IGM The Lyman-α forest 1.0 0.8 transmitted flux fraction e τ 0.6 0.4 e τ 0.2 0.0 0.25 0.20 0.15 0.10 0.05 no blazar heating intermediate b. h. 0.00 0 1000 2000 3000 4000 velocity [km s 1 ] Puchwein+ (2012)

Optical depths and temperatures Properties of blazar heating Thermal history of the IGM The Lyman-α forest effective optical depth τeff 0.5 0.4 0.3 0.2 no blazar heating weak blazar heating intermediate blazar heating strong blazar heating Viel et al. 2004 Tytler et al. 2004 FG 08 0.1 1.8 2.0 2.2 2.4 2.6 2.8 3.0 redshift z

Optical depths and temperatures Properties of blazar heating Thermal history of the IGM The Lyman-α forest effective optical depth τeff 0.5 0.4 0.3 0.2 no blazar heating weak blazar heating intermediate blazar heating strong blazar heating Viel et al. 2004 Tytler et al. 2004 FG 08 temperature T ( )/(10 4 K) 3.5 3.0 2.5 no blazar heating weak blazar heating intermediate blazar heating strong blazar heating Becker et al. 2011 2.0 0.1 1.8 2.0 2.2 2.4 2.6 2.8 3.0 redshift z 2.0 2.2 2.4 2.6 2.8 3.0 redshift z Puchwein+ (2012) Redshift evolutions of effective optical depth and IGM temperature match data only with additional heating, e.g., provided by blazars!

Properties of blazar heating Thermal history of the IGM The Lyman-α forest Ly-α flux PDFs and power spectra PDF of transmitted flux fraction 10 1 10 0 10 1 10 1 tuned UV background no blazar heating weak blazar heating intermediate blazar heating strong blazar heating Kim et al. 2007 z = 2.52 10 0 z = 2.94 10 1 0.0 0.2 0.4 0.6 0.8 1.0 transmitted flux fraction Puchwein+ (2012)

Voigt profile decomposition Properties of blazar heating Thermal history of the IGM The Lyman-α forest decomposing Lyman-α forest into individual Voigt profiles allows studying the thermal broadening of absorption lines

Properties of blazar heating Thermal history of the IGM The Lyman-α forest Voigt profile decomposition line width distribution 0.08 0.07 0.06 N HI > 10 13 cm 2 2.75 < z < 3.05 no blazar heating weak blazar heating intermediate blazar h. strong blazar heating Kirkman & Tytler 97 PDF of b [skm 1 ] 0.05 0.04 0.03 0.02 0.01 0.00 0 10 20 30 40 50 60 b[kms 1 ] Puchwein+ (2012)

Properties of blazar heating Thermal history of the IGM The Lyman-α forest Lyman-α forest in a blazar heated Universe improvement in modelling the Lyman-α forest is a direct consequence of the peculiar properties of blazar heating: heating rate independent of IGM density naturally produces the inverted T ρ relation that Lyman-α forest data demand recent and continuous nature of the heating needed to match the redshift evolutions of all Lyman-α forest statistics magnitude of the heating rate required by Lyman-α forest data the total energy output of TeV blazars (or equivalently 0.2% of that of quasars)

Outline Physics of blazar heating Formation of dwarf galaxies Puzzles in galaxy formation Conclusions 1 Physics of blazar heating TeV emission from blazars Plasma instabilities and magnetic fields Extragalactic gamma-ray background 2 Properties of blazar heating Thermal history of the IGM The Lyman-α forest 3 Formation of dwarf galaxies Puzzles in galaxy formation Conclusions

Entropy evolution Physics of blazar heating Formation of dwarf galaxies Puzzles in galaxy formation Conclusions temperature evolution 10 5 only photoheating standard BLF optimistic BLF T [K] 10 4 10 3 20 10 5 2 1 1 +z

Entropy evolution Physics of blazar heating Formation of dwarf galaxies Puzzles in galaxy formation Conclusions temperature evolution entropy evolution 10 5 only photoheating standard BLF 10 2 only photoheating standard BLF optimistic BLF optimistic BLF T [K] 10 4 K e [kev cm 2 ] 10 1 10 3 20 10 5 2 1 1 +z 0.1 20 10 5 2 1 1 +z C.P., Chang, Broderick (2012) evolution of entropy, K e = ktn 2/3 e, governs structure formation blazar heating: late-time, evolving, modest entropy floor

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Dwarf galaxy formation Jeans mass thermal pressure opposes gravitational collapse on small scales characteristic length/mass scale below which objects do not form

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Dwarf galaxy formation Jeans mass thermal pressure opposes gravitational collapse on small scales characteristic length/mass scale below which objects do not form hotter IGM higher IGM pressure higher Jeans mass: M J c3 s ρ 1/2 ( ) 1/2 T 3 IGM M J,blazar ρ M J,photo ( Tblazar T photo ) 3/2 30 depends on instantaneous value of c s

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Dwarf galaxy formation Jeans mass thermal pressure opposes gravitational collapse on small scales characteristic length/mass scale below which objects do not form hotter IGM higher IGM pressure higher Jeans mass: M J c3 s ρ 1/2 ( ) 1/2 T 3 IGM M J,blazar ρ M J,photo ( Tblazar T photo ) 3/2 30 depends on instantaneous value of c s filtering mass depends on full thermal history of the gas: accounts for delayed response of pressure in counteracting gravitational collapse in the expanding universe apply corrections for non-linear collapse

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Dwarf galaxy formation Filtering mass M F [ h -1 M O ] 10 12 10 10 10 8 10 6 1 + δ = 1, z reion = 10 blazar heating only photoheating linear theory non linear theory non-linear theory optimistic blazar standard blazar only photoheating M ~ 10 M F M ~ 10 M F 11 10 o o M F, blazar / M F 10 1 10 1 1 + z C.P., Chang, Broderick (2012)

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Peebles void phenomenon explained? mean density void, 1 + δ = 0.5 10 12 1 + δ = 1, z reion = 10 10 12 1 + δ = 0.5, z reion = 10 M F [ h -1 M O ] M F, blazar / M F 10 10 10 8 10 6 10 1 linear theory non-linear theory optimistic blazar standard blazar only photoheating 10 1 1 + z M F [ h -1 M O ] M F, blazar / M F 10 10 10 8 10 6 10 1 linear theory non-linear theory optimistic blazar standard blazar only photoheating 10 1 1 + z C.P., Chang, Broderick (2012) blazar heating efficiently suppresses the formation of void dwarfs within existing DM halos of masses < 3 10 11 M (z = 0) may reconcile the number of void dwarfs in simulations and the paucity of those in observations

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Missing satellite problem in the Milky Way Springel+ (2008) Dolphin+ (2005) Substructures in cold DM simulations much more numerous than observed number of Milky Way satellites!

When do dwarfs form? Formation of dwarf galaxies Puzzles in galaxy formation Conclusions 10 Myr M110 100 Myr 10 Gyr 1 Gyr Dolphin+ (2005) isochrone fitting for different metallicities star formation histories

When do dwarfs form? Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Dolphin+ (2005) red: τ form > 10 Gyr, z > 2

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Milky Way satellites: formation history and abundance satellite formation time Maccio & Fontanot (2010) late forming satellites (< 10 Gyr) not observed!

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Milky Way satellites: formation history and abundance satellite formation time satellite luminosity function Maccio & Fontanot (2010) late forming satellites (< 10 Gyr) not observed! no blazar heating: linear theory non linear theory Maccio+ (2010) blazar heating suppresses late satellite formation, may reconcile low observed dwarf abundances with CDM simulations

Galactic H I-mass function Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Mo+ (2005) H I-mass function is too flat (i.e., gas version of missing dwarf problem!) photoheating and SN feedback too inefficient IGM entropy floor of K 15 kev cm 2 at z 2 3 successful!

Conclusions on blazar heating Formation of dwarf galaxies Puzzles in galaxy formation Conclusions explains puzzles in high-energy astrophysics: lack of GeV bumps in blazar spectra without IGM B-fields unified TeV blazar-quasar model explains Fermi source counts and extragalactic gamma-ray background

Conclusions on blazar heating Formation of dwarf galaxies Puzzles in galaxy formation Conclusions explains puzzles in high-energy astrophysics: lack of GeV bumps in blazar spectra without IGM B-fields unified TeV blazar-quasar model explains Fermi source counts and extragalactic gamma-ray background novel mechanism; dramatically alters thermal history of the IGM: uniform and z-dependent preheating rate independent of density inverted T ρ relation quantitative self-consistent picture of high-z Lyman-α forest

Conclusions on blazar heating Formation of dwarf galaxies Puzzles in galaxy formation Conclusions explains puzzles in high-energy astrophysics: lack of GeV bumps in blazar spectra without IGM B-fields unified TeV blazar-quasar model explains Fermi source counts and extragalactic gamma-ray background novel mechanism; dramatically alters thermal history of the IGM: uniform and z-dependent preheating rate independent of density inverted T ρ relation quantitative self-consistent picture of high-z Lyman-α forest significantly modifies late-time structure formation: suppresses late dwarf formation (in accordance with SFHs): missing satellites, void phenomenon, H I-mass function group/cluster bimodality of core entropy values

When do clusters form? Formation of dwarf galaxies Puzzles in galaxy formation Conclusions 10 15 mass accretion history 5 mass accretion rates M 200c [ M O ] 10 14 10 13 z 0.5 = 0.42 z 0.5 = 0.55 z 0.5 = 0.70 z 0.5 = 0.86 z 0.5 = 0.95 S = d log M 200c / d log a (z) 4 3 2 1 10 12 1 + z 1 0 1 + z 1 C.P., Chang, Broderick (2012) most cluster gas accretes after z = 1, when blazar heating can have a large effect (for late forming objects)!

Entropy floor in clusters Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Cluster entropy profiles Cavagnolo+ (2009)

Entropy floor in clusters Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Cluster entropy profiles Planck stacking of optical clusters Cavagnolo+ (2009) Planck Collaboration (2011) Do optical and X-ray/Sunyaev-Zel dovich cluster observations probe the same population? (Hicks+ 2008, Planck Collaboration 2011)

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Entropy profiles: effect of blazar heating varying formation time varying cluster mass 1000 1000 z = 0 z = 0.5 z = 1 z = 2 M 200 = 3 x 10 13 M O M 200 = 1 x 10 14 M O, z = 0.5 M 200 = 3 x 10 13 M O, z = 0.5 M 200 = 1 x 10 13 M O, z = 0.5 K e [kev cm 2 ] 100 K e [kev cm 2 ] 100 10 optimistic blazar 10 optimistic blazar 0.01 0.10 1.00 r / R 200 0.01 0.10 1.00 r / R 200 C.P., Chang, Broderick (2012) assume big fraction of intra-cluster medium collapses from IGM: redshift-dependent entropy excess in cores greatest effect for late forming groups/small clusters

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Gravitational reprocessing of entropy floors Borgani+ (2005) 2 K 0= 25 kev cm preheated no preheating + floor no preheating greater initial entropy K 0 more shock heating greater increase in K 0 over entropy floor net K 0 amplification of 3-5 expect: median K e,0 150 kev cm 2 max. K e,0 600 kev cm 2 Borgani+ (2005)

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Cool-core versus non-cool core clusters cool hot Cavagnolo+ (2009)

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Cool-core versus non-cool core clusters cool hot early forming, t > t merger cool late forming, t merger < t cool Cavagnolo+ (2009)

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions Cool-core versus non-cool core clusters cool hot early forming, t > t merger cool late forming, t merger < t cool Cavagnolo+ (2009) time-dependent preheating + gravitational reprocessing CC-NCC bifurcation (two attractor solutions) need hydrodynamic simulations to confirm this scenario

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ]

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ]

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ]

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ]

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ]

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ]

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ]

Formation of dwarf galaxies Puzzles in galaxy formation Conclusions How efficient is heating by AGN feedback? 10 4 C.P., Chang, Broderick (2011) E b, 2500 (kt X = 5.9 kev) E b, 2500 (kt X = 3.5 kev) E b, 2500 (kt X = 2.0 kev) Ecav= 4PVtot [10 58 erg] 10 2 10 0 E b, 2500 (kt X = 1.2 kev) E b, 2500 (kt X = 0.7 kev) max K 0 10-2 cool cores non-cool cores 1 10 100 K e,0 [kev cm 2 ] AGNs cannot transform CC to NCC clusters (on a buoyancy timescale)