Researcher profile

M. Tornikoski

M. Tornikoski contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

A multiband study and exploration of the radio wave - $γ$-ray connection in 3C 84

Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $γ$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process regression. We find that the jet particle and magnetic field energy densities are in equipartition ($k_\textrm{r} = 1.08\pm0.18$). The jet apex is located $z_\textrm{91.5 GHz}= 22 - 645$ $R_\textrm{s}$ ($2 - 20 \times 10^{-3}$ pc) upstream of the 3 mm radio core; at that position, the magnetic field amplitude is $B_\textrm{core}^\textrm{91.5 GHz}= 3 - 10$ G. Our results are in good agreement with earlier studies, which utilised very-long-baseline interferometry. Furthermore, we investigate the temporal relation between the ejection of radio and $γ$-ray flares. Our results are in favour of the $γ$-ray emission being associated with the radio emission. We are able to tentatively connect the ejection of features identified at 43 and 86 GHz to prominent $γ$-ray flares. Finally, we compute the multiplicity parameter $λ$ and the Michel magnetisation $σ_\textrm{M}$ and find that they are consistent with a jet launched by the Blandford & Znajek 1977 mechanism.

preprint2022arXiv

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.

preprint2022arXiv

Jet kinematics in the transversely stratified jet of 3C 84 A two-decade overview

3C84 (NGC1275) is one of the brightest radio sources in the mm radio-bands, which led to a plethora of VLBI observations at numerous frequencies over the years. They reveal a two-sided jet structure, with an expanding but not well-collimated parsec-scale jet, pointing southward. High resolution mm-VLBI observations allow the study and imaging of the jet base on the sub-parsec scale. This could facilitate the investigation of the nature of the jet origin, also in view of the previously detected two-railed jet structure and east-west oriented core region seen with RadioAstron at 22 GHz. We produce VLBI images of this core and inner jet region, observed during the past twenty years at 15, 43, and 86 GHz. We determine the kinematics of the inner jet and ejected features at 43 and 86 GHz and compare their ejection times with radio and $γ$-ray variability. For the moving jet features, we find an average velocity of $β^\textrm{avg}_\textrm{app} = 0.055-0.22$c ($μ^\textrm{avg} = 0.04-0.18\,$mas/yr). From the time-averaged VLBI images at the three frequencies, we measure the transverse jet width along the bulk flow. On the $\leq 1.5$ parsec scale, we find a clear trend of the jet width being frequency dependent, with the jet being narrower at higher frequencies. This stratification is discussed in the context of a spine-sheath scenario, and is compared to other possible interpretations. From quasi-simultaneous observations at 43 and 86\,GHz we obtain spectral index maps, revealing a time-variable orientation of the spectral index gradient, due to structural variability of the inner jet.

preprint2022arXiv

New Tests of Millilensing in the Blazar PKS 1413+135

Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational millilensing, PKS 1413+135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the millilens itself. We discuss two possible candidates for the putative millilens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989-day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in August 2022, May 2025, and February 2028.

preprint2022arXiv

The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate

Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise&#39;&#39; of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6σ$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.

preprint2021arXiv

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet.

preprint2020arXiv

A decade of multi-wavelength observations of the TeV blazar 1ES 1215+303: Extreme shift of the synchrotron peak frequency and long-term optical-gamma-ray flux increase

Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed.

preprint2020arXiv

Absorbed relativistic jets in radio-quiet narrow-line Seyfert 1 galaxies

Narrow-line Seyfert 1 (NLS1) galaxies are peculiar active galactic nuclei (AGN). Most of them do not show strong radio emission, but recently seven radio-quiet (or -silent) NLS1s have been detected flaring multiple times at 37 GHz by the Metsähovi Radio Telescope, indicating the presence of relativistic jets in these peculiar sources. We observed them with the Karl G. Jansky Very Large Array (JVLA) in A configuration at 1.6, 5.2, and 9.0 GHz. Our results show that these sources are either extremely faint or not detected in the JVLA bands. At those frequencies, the radio emission from their relativistic jet must be absorbed, either via synchrotron self-absorption as it occurs in gigahertz-peaked sources or, more likely, via free-free absorption by a screen of ionized gas associated with starburst activity or shocks. Our findings cast new shadows on the radio-loudness criterion, which seems to be more and more frequently a misleading parameter. New high-frequency and high-resolution radio observations are essential to test our hypotheses.

preprint2020arXiv

Investigating the mini and giant radio flare episodes of Cygnus X-3

The microquasar Cygnus X-3 underwent a giant radio flare in April 2017, reaching a maximum flux of $\sim 16.5$ Jy at 8.5 GHz. We present results from a long monitoring campaign carried out with Medicina at 8.5, 18.6 and 24.1 GHz, in parallel to the Metsähovi radio telescope at 37 GHz, from 4 to 11 April 2017. We observe a spectral steepening from $α= 0.2$ to 0.5 (with $S_ν \propto ν^{-α}$) within 6 h around the epoch of the peak maximum of the flare, and rapid changes in the spectral slope in the following days during brief enhanced emission episodes while the general trend of the radio flux density indicated the decay of the giant flare. We further study the radio orbital modulation of Cyg X-3 emission associated with the 2017 giant flare and with six mini-flares observed in 1983, 1985, 1994, 1995, 2002 and 2016. The enhanced emission episodes observed during the decline of the giant flare at 8.5 GHz coincide with the orbital phase $ϕ\sim 0.5$ (orbital inferior conjunction). On the other hand the light curves of the mini-flares observed at $15-22$ GHz peak at $ϕ\sim 0$, except for the 2016 light curve which is shifted of 0.5 w.r.t. the other ones. We attribute the apparent phase shift to the variable location of the emitting region along the bent jet. This might be explained by the different accretion states of the flaring episodes (the 2016 mini-flare occurred in the hypersoft X-ray state).

preprint2020arXiv

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $γ$-ray (VERITAS, MAGIC), high-energy (HE) $γ$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Metsähovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline&#39; epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($δ\gtrsim 33$) and the size of the emission region ($ δ^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

preprint2020arXiv

The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135

PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a &#34;multizone&#34; model is needed to account for the high-energy emission.

preprint2020arXiv

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

preprint2019arXiv

A Fast Very High Energy $γ$-ray Flare from BL Lacertae during a Period of Multiwavelength activity in June 2015

The mechanisms producing fast variability of the $γ$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $γ$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $γ$-ray bands. In this paper we present the MAGIC VHE $γ$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $γ$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $γ$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $γ$-rays is very similar to two other observed VHE $γ$-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43\,GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution, based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the $γ$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $γ$-ray flare, but all have some tension with the multiwavelength observations.

preprint2019arXiv

Long term variability and correlation study of the blazar 3C 454.3 in radio, NIR and optical wavebands

We performed a long-term optical (B, V, R bands), infra-red (J and K bands) and radio band (15, 22, 37 GHz band) study on the flat spectrum radio quasar, 3C 454.3, using the data collected over a period of more than 8 years (MJD 54500--57500). The temporal variability, spectral properties and inter-waveband correlations were studied by dividing the available data into smaller segments with more regular sampling. This helped us constrain the size and the relative locations of the emission regions for different wavebands. Spectral analysis of the source revealed the interplay between the accretion disk and jet emission. The source predominantly showed a redder-when-brighter trend, though we observed a bluer-when-brighter trend at high flux levels which could be signatures of particle acceleration and radiative cooling. Significant correlations with near-zero lag were seen between various optical/infra-red bands, indicating that these emission regions are co-spatial. Correlations with a time lag of about 10--100 days are seen between optical/infra-red and radio bands indicating these emissions arise from different regions. We also observe the DCF peak lag change from year to year. We try to explain these differences using a curved jet model where the different emission regions have different viewing angles resulting in a frequency dependent Doppler factor. This variable Doppler factor model explains the variability timescales and the variation in DCF peak lag between the radio and optical emissions in different segments. Lags of 6-180 days are seen between emissions in various radio bands, indicating a core-shift effect

preprint2018arXiv

Multi-wavelength characterization of the blazar S5~0716+714 during an unprecedented outburst phase

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $γ$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

preprint2018arXiv

The extreme HBL behaviour of Markarian 501 during 2012

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

preprint2017arXiv

Multiwavelength observations of a VHE gamma-ray flare from PKS 1510-089 in 2015

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.