VHE gamma-ray emission from binary systems observed with the MAGIC telescopes Alicia López Oramas (IAC) O. Blanch, E. de Oña Wilhelmi, A. Fernández-Barral, D. Hadasch, J. Herrera, E. Moretti, P. Munar-Adrover, J.M. Paredes, M. Ribó (the MAGIC collaboration), P. Bordas, F. Brun and R. Zanin Sergi López/ Espai Astronomic 1
Transients with MAGIC A new gamma-ray binary: PSR J2032+4107 Searching VHE emission from massive microquasars SS 433 Cygnus X-1 Cygnus X-3 Searching VHE emission from LMXBs V404 Cygni NASA/Goddard/Feimer NASA/CXC/M.Weiss 2
The MAGIC telescopes La Palma Two telescopes, 17 m diameter Energy threshold (trigger) ~50 GeV Integral sensitivity E >290 GeV: (0.67 ± 0.04 )% of Crab Nebula flux in 50 hours (AlekiĆ et al. 2016) Energy resolution: 15-23 % Angular resolution: 0.1º AlekiĆ et al. 2016 VHE gamma-ray emission from binary systems observed with the MAGIC telescopes 3
Gamma-ray binaries: what are they? Definition: Bulk of the non-thermal emission lies in the γ-ray domain (E>1MeV) Only 7 (6+a new one!) display VHE emission Massive star + compact object (2 pulsars, 5 unidentified) Recently some microquasars detected at HE Pulsar wind Microquasar Mirabel 2006 Skilton et al, 2009 4
PSR J2032+4127/MT91 213 a new gamma-ray binary NASA/Goddard 5
VHE emission from PSR J2032 +4107 TeV 2032+4130 : unidentified source discovered by HEGRA (HEGRA 2002, Aharonian et al. 2005) and confirmed by Whipple (Konopelko et al. 2008), MAGIC, (Albert et al. 2008) Pulsar PSR J2032+4127 discovered in blind search by Fermi-LAT (Abdo et al. 2009). Young, high spin-down power TeV 2032+4130 possibly a wind nebula driven by the pulsar PSR J2032+4127 (Bednarek 2003, Aliu et al. 2014) Binary nature: PSR 2032+4127 associated to Be star MT91 213 (Lyne et al. 2005) Konopelko et al. 2008 NASA/Goddard 6
VHE emission from PSR J2032 +4107 Extremely eccentric binary: Orbital period: 50 years (Ho et al. 2017) Periastron: November 2017 (MJD 58070) Monitoring of the source: MAGIC +VERITAS joined campaign 53.7 h (2016 May - September) 34.2 h (2017 June - December) September 2017 (ATel #10810): TeV gamma-ray flux increased a factor 2 wrt June-August 2017 average November 2017 (ATel #10971) Periastron passage (MJD 58069.8): Flux increased almost a factor 10 wrt the average flux in June- August in only 1.9 h Dec (J2000) +42 00 0 30 0 +41 00 0 t: 53.7 h t: 1.9 h MAGIC J2031+4134 MAGIC J2032+4127 32 m 20 h 34 m RA (J2000) PSR J2032+4127 Cygnus X-3 Preliminary 30 m 20 15 10 5 0 5 Significance [ ] 7
VHE emission from PSR J2032 +4107 2016-Aug 2016-Dec 2017-Apr 2017-Aug 2017-Dec 2018-Apr 2017-Sep 2017-Oct 2017-Nov 2017-Dec 2018-Jan 6 X-ray Swift 6 X-ray Swift VERITAS and MAGIC. in prep E FX (0.3-10 kev) [ 10 12 erg cm 2 s 1 ] 4 2 Preliminary E FX (0.3-10 kev) [ 10 12 erg cm 2 s 1 ] 4 2 Preliminary FVHE (E > 200 GeV) [ 10 11 cm 2 s 1 ] 0 2.0 1.5 1.0 0.5 VHE -ray Date [mjd] VERITAS MAGIC Preliminary FVHE (E > 200 GeV) [ 10 11 cm 2 s 1 ] 0 2.5 2.0 2.0 1.5 1.5 1.0 1.0 0.5 0.5 EFVHE (E > 100 GeV) [ 10 11 erg cm 2 s 1 ] VHE -ray VERITAS Date [mjd] MAGIC Preliminary 2.5 2.0 1.5 1.0 0.5 EFVHE (E > 100 GeV) [ 10 11 erg cm 2 s 1 ] X-rays Swift LC: 0.0 57600 57750 57900 58050 58200 Date [MJD] Peaked about 30 days before periastron Gradually decreasing, minimum at periastron Recovery over the next 30 days (punctual flare 15 days after periastron) 0.0 0.0 VHE LC: 58000 58020 58040 58060 58080 58100 0.0 58120 Date [MJD] Flux peaked at periastron 7 days after periastron: sharp decrease of the flux compatible with the baseline emission. Flux recovered to periastron level few days later Sharp dip after periastron likely caused by γ-γ absorption 8
VHE emission from PSR J2032 +4107 2016: only steady emission component was present (associated to TeV 2032+4130) 2017: include contributions from both the steady and variable sources A joint fit was conducted to determine the spectral properties of the emission from the binary above the (baseline) background from the PWN (modeled as a PWL). E 2 df/de [TeV cm 2 s 1 ] SED compatible with predictions by Bednarek et al. (2018) 10 11 10 12 10 13 Preliminary Baseline 2016 2017 Average E 2 df/de [TeV cm 2 s 1 ] 10 11 10 12 10 13 VERITAS and MAGIC. in prep Baseline Low State High State Preliminary 10 14 10 1 10 0 10 1 Energy [TeV] Fall 2017 data: PWL with exponential cut-off (EPWL) at E =1.40 ± 0.97 TeV 10 14 10 1 10 0 10 1 Energy [TeV] Low state: EPWL with cut-off at E =0.58 ± 0.33 TeV High state: PWL 9
SS433, Cygnus X-1 & Cygnus X-3 Searching for VHE emission from massive microquasars NASA/CXC/M.Weiss 10
SS433/W50 complex NASA/CXC/U.Amsterdam/S.Migliari et al Only galactic super-critical accretor Persistent hadronic relativistic jets Porb ~13 days, Ppre ~162 days Embeded in W50 nebula Interaction regions: X-ray & radio non-thermal emission Detected at HE (Bordas et al. 2015) Prediction of gamma-ray fluxes at E>800GeV due to pp interactions (Reynoso et al. 2008) Absorption of putative emission in ~80% of the orbit Best opportunity for observations: Φpre=0.91-0.09 Goodall et al, 2011 11
SS 433/W50 at VHE Dec (J2000) dn/de [cm 2 s 1 TeV 1 ] +6 +5 +4 +3 10 7 10 8 10 9 10 10 10 11 10 12 10 13 10 14 MAGIC and H.E.S.S.., 2017, A&A, 612, A14 e3 e1 e2 SS 433 18 m 15 m 12 m RA (J2000) MAGIC H.E.S.S. MAGIC+H.E.S.S. 10 2 10 3 E [GeV] MGRO J1908+06 w1 w2 09 m + 19 h 06 m Ψ prec =0.9 0.1 Crab Nebula 6 5 4 3 2 1 0 1 2 Significance MAGIC + H.E.S.S. campaign No significant excess detected Upper limits (95% C.L.) compatible with predictions by Reynoso et al. 2008 Hadronic scenario Efficiency in transferring jet kinetic power to relativistic protons: We can constrain: qp 2.5 10-5 Leptonic scenario Interaction regions: X-rays (synchrotron origin ) -> presence of electrons up to 50 TeV Expected fluxes (Bordas et al. (2009): roughly at level of the reported ULs HE electrons may loose most of their E preventing an efficient production of gamma-rays ULs: constrains on the magnetic field: lower limit of 20-25μG 12
The microquasar Cygnus X-1 O9.7 Iab star + BH Highly collimated one-side jet (Stirling et al. 2001) Surrounded by radio/optical nebula (Gallo et al. 2005, Russell et al. 2007) Three transient episodes with AGILE during HS and IS(Bulgarelli et al. 2008, Sabatini et al. 2010, 2013) Detected at HE during HS: 7.5yr Fermi-LAT data (Zanin et al. 2016) Hint of emission with MAGIC: 4σ in 80 min (Albert et al. 2006) Simultaneously with hard X-ray flare During HS and SUPC Gallo et al. 2005 HS SS Zanin et al. 2016 13
Cygnus X-1: results 100 h (2007-2014) of MAGIC observations mainly at HS (83h) MAGIC coll., 2017, MNRAS, 472, 3474 No significant excess at either X- ray state for steady, orbital or daily basis emission No emission above 200 GeV due to interaction between jet and ISM Jet-medium interaction discarded as possible region for VHE emission above MAGIC sensitivity level: not affected by γ γ absorption Transient emission (Albert et al. 2007) still possible at binary scale 14
Galactic jets: Cygnus X-3 The system WR+compact object Short orbit: 4.8 h Major radio flares during soft state (Szostek et al. 2008) MAGIC coll, ICRC Proceedings 2017 (arxiv 1709.01725) HE emission detected: AGILE (Tavani et al. 2009) and Fermi-LAT (Abdo et al. 2009) Latest results August-September 2016 (radio) flare: 70 h observations during SS No steady/transient VHE emission (E>100 and E>300 GeV) Not orbital/daily modulation Explanation Extremely high absorption due to the WR MAGIC, ICRC 2017 VHE gamma rays, if produced, are originated inside the binary scale and not at the radio-emitting regions of the jets far from the compact object 15
V404 Cygni Searching for VHE emission from LMXBs ESA / ATG medialab 16
The LMXB V404 Cygni The system ~1 M star + 9-15 M BH (Khargharia et al. 2010) Orbital period: 6.5 days (Casares et al. 1992) What happened? Major outburst in June 2015 after 26 years in quiescent state MAGIC was triggered by INTEGRAL: ~11h between June 18-27 Selection of the flaring times for MAGIC: Bayesian block analysis on the INTEGRAL LC MAGIC coll. 2017, MNRAS, 471, 1688. Latest results 7h distributed in different nights No signal detected (0.08σ), UL (E= 200-1250 GeV): 4.8x10-12 ph cm -2 s -1 17
V404 Cygni: discussion MAGIC Coll. 2017, MNRAS, 471, 1688. Evidences of jet emission in optical observations on 26th June: + hint of detection (~4σ) in Fermi-LAT data (Loh et al. 2016) + giant radio flare + increase of hardness ratio in X-ray band + optical fast variability Jet environment dramatically changed No detection in MAGIC data (~1 hr taken simultaneously) Luminosity MAGIC upper limits 2 10 33 erg s 1, in contrast with extreme luminosity emitted in X-ray band 2 10 38 erg s 1 Estimated γ-ray opacity during flaring period & non-detection: inefficient γ-ray emission (0.003%) in V404 Cyg jets if VHE emitter is located r >1 10 10 cm from the compact object. 18
Summary PSR J2032+4127: a new gamma-ray binary! TeV emission of PSR J2032+4107 during periastron passage. Flux increased factor 10 wrt than baseline measured in June-August. SED described by a EPWL in low state (before periastron) and PWL during high state (during periastron passage). Baseline emission as a PWL Search for VHE emission from microquasars: No significant excess from SS433/W50 complex. Constrains on particle acceleration and magnetic field No detection of Cygnus X-1 at TeV energies: ULs for steady, daily and separated X-ray states (including phase-folded analysis) Cygnus X-3 shows extremely high absorption due to the WR. VHE emission possible at binary scale ULs on the LMXB V404 Cygni during the June 2015 major outburst: inefficient γ-ray emission 19
Gamma-ray binaries: state-of-the-art System Star spectral type Compact object Porb [days] HE emission VHE emission PSR B1259-53 Be 48ms pulsar 1236.72 yes yes LS 5039 O - 3.91 yes yes LS I +61 303 Be - 26.49 yes yes HESS J0632+057 Be - 315.50 yes yes extragalactic new binary FGL J1018.6-5856 O - 16.58 yes yes LMC P-3 O - 10.2 yes yes PSR J2032+4127 Be 143 ms pulsar 50 years yes yes HE emitters Cygnus X-1 O BH yes no (4σ hint once) Cygnus X-3 WR BH? yes no SS 433 A BH yes no 20
VHE gamma-ray emission from binary systems observed with the MAGIC telescopes Thanks for your attention Sergi López/ Espai Astronomic Alicia López Oramas (IAC) O. Blanch, E. de Oña Wilhelmi, A. Fernández-Barral, D. Hadasch, J. Herrera, E. Moretti, P. Munar-Adrover, J.M. Paredes, M. Ribó (the MAGIC collaboration), P. Bordas, F. Brun and R. Zanin 21