Search

We present the results of optical photometric observations of three extreme TeV blazars, 1ES 0229$+$200, 1ES 0414$+$009, and 1ES 2344$+$514, taken with two telescopes (1.3 m Devasthal Fast Optical Telescope, and 1.04 m Sampuranand Telescope) in India and two (1.4 m Milankovi\'{c} telescope and 60 cm Nedeljkovi\'{c} telescope) in Serbia during 2013--2019. We investigated their flux and spectral variability on diverse timescales. We examined a total of 36 intraday $R-$band light curves of these blazars for flux variations using the power-enhanced {\it F}-test and the nested ANOVA test. No significant intraday variation was detected on 35 nights, and during the one positive detection the amplitude of variability was only 2.26 per cent. On yearly timescales, all three blazars showed clear flux variations in all optical wavebands. The weighted mean optical spectral index ($\alpha_{BR}$), calculated using $B - R$ color indices, for 1ES 0229$+$200 was 2.09 $\pm$ 0.01. We also estimated the weighted mean optical spectral indices of 0.67 $\pm$ 0.01 and 1.37 $\pm$ 0.01 for 1ES 0414$+$009, and 1ES 2344$+$514, respectively, by fitting a single power-law ($F_{\nu} \propto \nu^{-\alpha}$) in their optical ({\it VRI}) spectral energy distributions. A bluer-when-brighter trend was only detected in the blazar 1ES 0414$+$009. We briefly discuss different possible physical mechanisms responsible for the observed flux and spectral changes in these blazars on diverse timescales., Comment: 17 pages, 12 figures, 8 tables, accepted for publication in MNRAS journal

Full list of authors: Acciari, V. A.; Ansoldi, S.; Antonelli, L. A.; Arbet Engels, A.; Artero, M.; Asano, K.; Babić, A.; Baquero, A.; Barres de Almeida, U.; Barrio, J. A.; Batković, I.; Becerra González, J.; Bednarek, W.; Bellizzi, L.; Bernardini, E.; Bernardos, M.; Berti, A.; Besenrieder, J.; Bhattacharyya, W.; Bigongiari, C.; Blanch, O.; Bošnjak, Ž.; Busetto, G.; Carosi, R.; Ceribella, G.; Cerruti, M.; Chai, Y.; Chilingarian, A.; Cikota, S.; Colak, S. M.; Colombo, E.; Contreras, J. L.; Cortina, J.; Covino, S.; D'Amico, G.; D'Elia, V.; da Vela, P.; Dazzi, F.; de Angelis, A.; de Lotto, B.; Delfino, M.; Delgado, J.; Delgado Mendez, C.; Depaoli, D.; di Pierro, F.; di Venere, L.; Do Souto Espiñeira, E.; Dominis Prester, D.; Donini, A.; Doro, M.; Fallah Ramazani, V.; Fattorini, A.; Ferrara, G.; Fonseca, M. V.; Font, L.; Fruck, C.; Fukami, S.; García López, R. J.; Garczarczyk, M.; Gasparyan, S.; Gaug, M.; Giglietto, N.; Giordano, F.; Gliwny, P.; Godinović, N.; Green, J. G.; Green, D.; Hadasch, D.; Hahn, A.; Heckmann, L.; Herrera, J.; Hoang, J.; Hrupec, D.; Hütten, M.; Inada, T.; Inoue, S.; Ishio, K.; Iwamura, Y.; Jiménez, I.; Jormanainen, J.; Jouvin, L.; Kajiwara, Y.; Karjalainen, M.; Kerszberg, D.; Kobayashi, Y.; Kubo, H.; Kushida, J.; Lamastra, A.; Lelas, D.; Leone, F.; Lindfors, E.; Lombardi, S.; Longo, F.; López-Coto, R.; López-Moya, M.; López-Oramas, A.; Loporchio, S.; Machado de Oliveira Fraga, B.; Maggio, C.; Majumdar, P.; Makariev, M.; Mallamaci, M.; Maneva, G.; Manganaro, M.; Maraschi, L.; Mariotti, M.; Martínez, M.; Mazin, D.; Menchiari, S.; Mender, S.; Mićanović, S.; Miceli, D.; Miener, T.; Minev, M.; Miranda, J. M.; Mirzoyan, R.; Molina, E.; Moralejo, A.; Morcuende, D.; Moreno, V.; Moretti, E.; Neustroev, V.; Nigro, C.; Nilsson, K.; Nishijima, K.; Noda, K.; Nozaki, S.; Ohtani, Y.; Oka, T.; Otero-Santos, J.; Paiano, S.; Palatiello, M.; Paneque, D.; Paoletti, R.; Paredes, J. M.; Pavletić, L.; Peñil, P.; Perennes, C.; Persic, M.; Prada Moroni, P. G.; Prandini, E.; Priyadarshi, C.; Puljak, I.; Ribó, M.; Rico, J.; Righi, C.; Rugliancich, A.; Saha, L.; Sahakyan, N.; Saito, T.; Sakurai, S.; Satalecka, K.; Saturni, F. G.; Schmidt, K.; Schweizer, T.; Sitarek, J.; Šnidarić, I.; Sobczynska, D.; Spolon, A.; Stamerra, A.; Strom, D.; Strzys, M.; Suda, Y.; Surić, T.; Takahashi, M.; Tavecchio, F.; Temnikov, P.; Terzić, T.; Teshima, M.; Tosti, L.; Truzzi, S.; Tutone, A.; Ubach, S.; van Scherpenberg, J.; Vanzo, G.; Vazquez Acosta, M.; Ventura, S.; Verguilov, V.; Vigorito, C. F.; Vitale, V.; Vovk, I.; Will, M.; Wunderlich, C.; Zarić, D.; Baack, D.; Balbo, M.; Biederbeck, N.; Biland, A.; Bretz, T.; Buss, J.; Dorner, D.; Eisenberger, L.; Elsaesser, D.; Hildebrand, D.; Iotov, R.; Mannheim, K.; Neise, D.; Noethe, M.; Paravac, A.; Rhode, W.; Schleicher, B.; Sliusar, V.; Walter, R.; D'Ammando, F.; Horan, D.; Lien, A. Y.; Baloković, M.; Madejski, G. M.; Perri, M.; Verrecchia, F.; Leto, C.; Lähteenmäki, A.; Tornikoski, M.; Ramakrishnan, V.; Järvelä, E.; Vera, R. J. C.; Chamani, W.; Villata, M.; Raiteri, C. M.; Gupta, A. C.; Pandey, A.; Fuentes, A.; Agudo, I.; Casadio, C.; Semkov, E.; Ibryamov, S.; Marchini, A.; Bachev, R.; Strigachev, A.; Ovcharov, E.; Bozhilov, V.; Valcheva, A.; Zaharieva, E.; Damljanovic, G.; Vince, O.; Larionov, V. M.; Borman, G. A.; Grishina, T. S.; Hagen-Thorn, V. A.; Kopatskaya, E. N.; Larionova, E. G.; Larionova, L. V.; Morozova, D. A.; Nikiforova, A. A.; Savchenko, S. S.; Troitskiy, I. S.; Troitskaya, Y. V.; Vasilyev, A. A.; Merkulova, O. A.; Chen, W. P; Samal, M.; Lin, H. C.; Moody, J. W.; Sadun, A. C.; Jorstad, S. G.; Marscher, A. P.; Weaver, Z. R.; Feige, M.; Kania, J.; Kopp, M.; Kunkel, L.; Reinhart, D.; Scherbantin, A.; Schneider, L.; Lorey, C.; Acosta-Pulido, J. A.; Carnerero, M. I.; Carosati, D.; Kurtanidze, S. O.; Kurtanidze, O. M.; Nikolashvili, M. G.; Chigladze, R. A.; Ivanidze, R. Z.; Kimeridze, G. N.; Sigua, L. A.; Joner, M. D.; Spencer, M.; Giroletti, M.; Marchili, N.; Righini, S.; Rizzi, N.; Bonnoli, G.; MAGIC Collaboration; Fact Collaboration.-- This is an Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Aims. We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with a special focus on the multi-band flux correlations. Methods. The dataset has been collected through an extensive multi-wavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina, and Metsahovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. Results. The very-high-energy (VHE; E 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow dierent patterns, partly due to substantial changes in the Compton dominance for a few days without a simultaneous increase in the X-ray flux (i.e., orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. The peak frequency of the synchrotron bump varies by two orders of magnitude throughout the campaign. Our multi-band correlation study also hints at an anti-correlation between UV-optical and X-ray at a significance higher than 3. A VHE flare observed on MJD 57788 (2017 February 4) shows gamma-ray variability on multi-hour timescales, with a factor ten increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario.We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors, from 0 min = 2104 to 0 max = 6105. © 2021 Georg Thieme Verlag. All rights reserved., The MAGIC Collaboration would like to thank the Instituto de Astrofísica de Canarias for the excellent working conditions at the Observatorio del Roque de los Muchachos in La Palma. The financial support of the German BMBF, MPG and HGF; the Italian INFN and INAF; the Swiss National Fund SNF; the ERDF under the Spanish Ministerio de Ciencia e Innovación (MICINN) (FPA2017-87859-P, FPA2017-85668-P, FPA2017-82729-C6-5-R, FPA2017-90566-REDC, PID2019-104114RB-C31, PID2019-104114RB-C32, PID2019-105510GB-C31,PID2019-107847RB-C41, PID2019-107847RB-C42, PID2019-107847RB-C44, PID2019-107988GB-C22); the Indian Department of Atomic Energy; the Japanese ICRR, the University of Tokyo, JSPS, and MEXT; the Bulgarian Ministry of Education and Science, National RI Roadmap Project DO1-268/16.12.2019 and the Academy of Finland grant nr. 320045 is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia “Severo Ochoa” SEV-2016-0588, SEV-2017-0709 and CEX2019-000920-S, and “María de Maeztu” CEX2019-000918-M, the Unidad de Excelencia “María de Maeztu” MDM-2015-0509-18-2 and the “la Caixa” Foundation (fellowship LCF/BQ/PI18/11630012) and by the CERCA program of the Generalitat de Catalunya; by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02; by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3; the Polish National Research Centre grant UMO-2016/22/M/ST9/00382; and by the Brazilian MCTIC, CNPq and FAPERJ. The important contributions from ETH Zurich grants ETH-10.08-2 and ETH-27.12-1 as well as the funding by the Swiss SNF and the German BMBF (Verbundforschung Astro- und Astroteilchenphysik) and HAP (Helmoltz Alliance for Astroparticle Physics) are gratefully acknowledged. Part of this work is supported by Deutsche Forschungsgemeinschaft (DFG) within the Collaborative Research Center SFB 876 “Providing Information by Resource-Constrained Analysis”, project C3. We are thankful for the very valuable contributions from E. Lorenz, D. Renker and G. Viertel during the early phase of the project. We thank the Instituto de Astrofísica de Canarias for allowing us to operate the telescope at the Observatorio del Roque de los Muchachos in La Palma, the Max-Planck-Institut für Physik for providing us with the mount of the former HEGRA CT3 telescope, and the MAGIC collaboration for their support. The Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l’Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d’Études Spatiales in France. This work performed in part under DOE Contract DE-AC02-76SF00515. This work made use of data from the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. We thank the NuSTAR Operations, Software, and Calibration teams for support with the execution and analysis of these observations. This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Science Data Center (ASDC; Italy) and the California Institute of Technology (USA). This research has also made use of the XRT Data Analysis Software (XRTDAS) developed under the responsibility of the ASI Science Data Center (ASDC), Italy. A.A.E and D.P acknowledge support from the Deutsche Forschungs gemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 390783311. M. B. acknowledges support from the YCAA Prize Postdoctoral Fellowship and from the Black Hole Initiative at Harvard University, which is funded in part by the Gordon and Betty Moore Foundation (grant GBMF8273) and in part by the John Templeton Foundation. This publication makes use of data obtained at the Metsähovi Radio Observatory, operated by Aalto University in Finland. This research has made use of data from the OVRO 40-m monitoring program (Richards et al. 2011) which is supported in part by NASA grants NNX08AW31G, NNX11A043G, and NNX14AQ89G and NSF grants AST-0808050 and AST-1109911. I.A. acknowledges financial support from the Spanish “Ministerio de Ciencia e Innovación” (MCINN) through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía-CSIC (SEV-2017-0709). Acquisition and reduction of the MAPCAT data was supported in part by MICINN through grants AYA2016-80889-P and PID2019-107847RB-C44. The MAPCAT observations were carried out at the German-Spanish Calar Alto Observatory, which is jointly operated by Junta de Andalucía and Consejo Superior de Investigaciones Científicas. C.C. acknowledges support from the European Research Council (ERC) under the European Union Horizon 2020 research and innovation program under the grant agreement No 771282. This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants KP-06-H28/3 (2018), KP-06-H38/4 (2019) and KP-06-KITAJ/2 (2020). We acknowledge support by Bulgarian National Science Fund under grant DN18-10/2017 and National RI Roadmap Projects DO1-277/16.12.2019 and DO1-268/16.12.2019 of the Ministry of Education and Science of the Republic of Bulgaria. This research was supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia (contract No 451-03-68/2020-14/200002). G.D. acknowledges observing grant support from the Institute of Astronomy and Rozhen NAO BAS through the bilateral joint research project “Gaia Celestial Reference Frame (CRF) and fast variable astronomical objects” (2020–2022, head – G. Damljanovic). The BU group was supported in part by NASA Fermi guest investigator program grants 80NSSC19K1505 and 80NSSC20K1566. This study was based in part on observations conducted using the 1.8 m Perkins Telescope Observatory (PTO) in Arizona, which is owned and operated by Boston University. This article is partly based on observations made with the LCOGT Telescopes, one of whose nodes is located at the Observatorios de Canarias del IAC on the island of Tenerife in the Observatorio del Teide. This article is also based partly on data obtained with the STELLA robotic telescopes in Tenerife, an AIP facility jointly operated by AIP and IAC. The Abastumani team acknowledges financial support by the Shota Rustaveli National Science Foundation under contract FR-19-6174. Based on observations with the Medicina telescope operated by INAF – Istituto di Radioastronomia.

We present the variability study of the lowest-luminosity Seyfert 1 galaxy NGC 4395 based on the photometric monitoring campaigns in 2017 and 2018. Using 22 ground-based and space telescopes, we monitored NGC 4395 with a $\sim$5 minute cadence during a period of 10 days and obtained light curves in the UV, V, J, H, and K/Ks bands as well as the H$\alpha$ narrow-band. The RMS variability is $\sim$0.13 mag on \emph{Swift}-UVM2 and V filter light curves, decreasing down to $\sim$0.01 mag on K filter. After correcting for continuum contribution to the H$\alpha$ narrow-band, we measured the time lag of the H$\alpha$ emission line with respect to the V-band continuum as ${55}^{+27}_{-31}$ to ${122}^{+33}_{-67}$ min. in 2017 and ${49}^{+15}_{-14}$ to ${83}^{+13}_{-14}$ min. in 2018, depending on the assumption on the continuum variability amplitude in the H$\alpha$ narrow-band. We obtained no reliable measurements for the continuum-to-continuum lag between UV and V bands and among near-IR bands, due to the large flux uncertainty of UV observations and the limited time baseline. We determined the AGN monochromatic luminosity at 5100\AA\ $\lambda L_\lambda = \left(5.75\pm0.40\right)\times 10^{39}\,\mathrm{erg\,s^{-1}}$, after subtracting the contribution of the nuclear star cluster. While the optical luminosity of NGC 4395 is two orders of magnitude lower than that of other reverberation-mapped AGNs, NGC 4395 follows the size-luminosity relation, albeit with an offset of 0.48 dex ($\geq$2.5$\sigma$) from the previous best-fit relation of Bentz et al. (2013)., Comment: Accepted for publication in ApJ (Feb. 23rd, 2020). 18 pages, 10 figures

We report the results of our optical (VRI) photometric observations of the TeV blazar 1ES 0806$+$524 on 153 nights during 2011-2019 using seven optical telescopes in Europe and Asia. We investigated the variability of the blazar on intraday as well as on long-term timescales. We examined eighteen intraday light curves for flux and color variations using the most reliable power-enhanced F-test and the nested ANOVA test. Only on one night was a small, but significant, variation found, in both $V$ band and $R$ band light curves. The $V-R$ color index was constant on every one of those nights. Flux density changes of around 80 % were seen over the course of these eight years in multiple bands. We found a weighted mean optical spectral index of 0.639$\pm$0.002 during our monitoring period by fitting a power law ($F_{\nu} \propto \nu^{-\alpha}$) in 23 optical ($VRI$) spectral energy distributions of 1ES 0806$+$524. We discuss different possible mechanisms responsible for blazar variability on diverse timescales., Comment: 14 pages, 6 figures, 7 tables, accepted for publication in ApJ

In an effort to locate the sites of emission at different frequencies and physical processes causing variability in blazar jets, we have obtained high time-resolution observations of BL Lacertae over a wide wavelength range: with the Transiting Exoplanet Survey Satellite (TESS) at 6000–10000 Å with 2 minute cadence; with the Neil Gehrels Swift satellite at optical, UV, and X-ray bands; with the Nuclear Spectroscopic Telescope Array at hard X-ray bands; with the Fermi Large Area Telescope at γ-ray energies; and with the Whole Earth Blazar Telescope for measurement of the optical flux density and polarization. All light curves are correlated, with similar structure on timescales from hours to days. The shortest timescale of variability at optical frequencies observed with TESS is ∼0.5 hr. The most common timescale is 13 ± 1 hr, comparable with the minimum timescale of X-ray variability, 14.5 hr. The multiwavelength variability properties cannot be explained by a change solely in the Doppler factor of the emitting plasma. The polarization behavior implies that there are both ordered and turbulent components to the magnetic field in the jet. Correlation analysis indicates that the X-ray variations lag behind the γ-ray and optical light curves by up to ∼0.4 day. The timescales of variability, cross-frequency lags, and polarization properties can be explained by turbulent plasma that is energized by a shock in the jet and subsequently loses energy to synchrotron and inverse Compton radiation in a magnetic field of strength ∼3 G.

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability., 11 pages, 6 figures, 4 tables; Accepted for publication in MNRAS

We present the results of photometric observations of three TeV blazars, 3C 66A, S5 0954+658 and BL Lacertae, during the period 2013--2017. Our extensive observations were performed in a total of 360 nights which produced $\sim$6820 image frames in BVRI bands. We study flux and spectral variability of these blazars on these lengthy timescales. We also examine the optical Spectral Energy Distributions of these blazars, which are crucial in understanding the emission mechanism of long-term variability in blazars. All three TeV blazars exhibited strong flux variability during our observations. The colour variations are mildly chromatic on long timescales for two of them. The nature of the long-term variability of 3C 66A and S5 0954+658 is consistent with a model of a non-thermal variable component that has a continuous injection of relativistic electrons with power law distributions around 4.3 and 4.6, respectively. However, the long-term flux and colour variability of BL Lac suggests that these can arise from modest changes in velocities or viewing angle toward the emission region, leading to variations in the Doppler boosting of the radiation by a factor ~1.2 over the period of these observations., 13 pages, 15 figures, Accepted for publication in MNRAS

Context. The flat-spectrum radio quasar 4C +71.07 is a high-redshift (z = 2.172), γ-loud blazar whose optical emission is dominated by thermal radiation from the accretion disc. Aims. 4C +71.07 has been detected in outburst twice by the AGILE γ-ray satellite during the period from the end of October to mid-November 2015, when it reached a γ-ray flux of the order of F(E > 100 MeV)=(1.2 ± 0.3)×10 photons cm s and F(E > 100 MeV)=(3.1 ± 0.6)×10 photons cm s, respectively, allowing us to investigate the properties of the jet and the emission region. Methods. We investigated its spectral energy distribution by means of almost-simultaneous observations covering the cm, mm, near-infrared, optical, ultraviolet, X-ray, and γ-ray energy bands obtained by the GASP-WEBT Consortium and the Swift, AGILE, and Fermi satellites. Results. The spectral energy distribution of the second γ-ray flare (whose energy coverage is more dense) can be modelled by means of a one-zone leptonic model, yielding a total jet power of about 4 × 10 erg s. Conclusions. During the most prominent γ-ray flaring period our model is consistent with a dissipation region within the broad-line region. Moreover, this class of high-redshift, flat-spectrum radio quasars with high-mass black holes might be good targets for future γ-ray satellites such as e-ASTROGAM. © ESO 2019., AGILE is an ASI space mission developed with programmatic support by INAF and INFN. We acknowledge partial support through the ASI grant no. I/028/12/0. SV and PR acknowledge contract ASI-INAF I/004/11/0 and INAF/IASF Palermo where most of the work was carried out. SV acknowledges financial contribution from the agreement ASI-INAF no. 2017-14-H.0. Part of this work is based on archival data, software, or online services provided by the ASI SPACE SCIENCE DATA CENTER (ASI-SSDC). SV and PR thank Leonardo Barzaghi and Sara Baitieri for useful discussions. The Osservatorio di Torino team acknowledges the financial contribution from the agreement ASI-INAF No. 2017-14-H.0 and from the contract PRIN-SKA-CTA-INAF 2016. OMK acknowledges financial support from the Shota Rustaveli National Science Foundation under contract FR/217950/16 and grants NSFC11733001, NSFCU1531245. IA acknowledges support from a Ramón y Cajal grant of the Ministerio de Economía y Compet-itividad (MINECO) of Spain. The research at the IAA–CSIC was supported in part by the MINECO through grants AYA2016–80889–P, AYA2013–40825–P, and AYA2010–14844, and by the regional government of Andalucía through grant P09–FQM–4784. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). Calar Alto Observatory is jointly operated by the MPIA and the IAA-CSIC. This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grant DN 08-1/2016. The St. Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. AZT-24 observations are made within an agreement among the Pulkovo, Rome, and Teramo observatories. GD and OV gratefully acknowledge the observing grant support from the Institute of Astronomy and Rozhen National Astronomical Observatory, Bulgaria Academy of Sciences, via bilateral joint research project “Observations of ICRF radio-sources visible in optical domain” (PI G. Damljanovic). This work is a part of Project No. 176011 (“Dynamics and kinematics of celestial bodies and systems”), No. 176004 (“Stellar physics”) and No. 176021 (“Visible and invisible matter in nearby galaxies: theory and observations”) supported by the Ministry of Education, Science, and Technological Development of the Republic of Serbia. The Maidanak Observatory team acknowledges support from Uzbekistan Academy of Sciences grants No. F2-FA-F027 and F.4-16.

We report on a multiband variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 19. The study uses, among others, data from GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), Swift, Nuclear Spectroscopic Telescope Array (NuSTAR), Fermi Large Area Telescope, Very Energetic Radiation Imaging Telescope Array System (VERITAS), and Major Atmospheric Gamma Imaging Cherenkov (MAGIC). The large blazar activity and the 43 hr of simultaneous NuSTAR and MAGIC/VERITAS observations permitted variability studies on 15 minute time bins over three X-ray bands (3-7 keV, 7-30 keV, and 30-80 keV) and three very-high-energy (VHE; >0.1 TeV) 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 of 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-rays 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 characteristics of the correlation change substantially over time and across energy bands. Our findings favor 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. © 2020. The American Astronomical Society. All rights reserved., The financial support of the German BMBF and MPG, the Italian INFN and INAF, the Swiss National Fund SNF, the ERDF under the Spanish MINECO (FPA2017-87859-P, FPA2017-85668-P, FPA2017-82729-C6-2-R, FPA2017-827 29-C6-6-R, FPA2017-82729-C6-5-R, AYA2015-71042-P, AYA2016-76012-C3-1-P, ESP2017-87055-C2-2-P, and FPA 201790566REDC), the Indian Department of Atomic Energy, the Japanese JSPS and MEXT, the Bulgarian Ministry of Education and Science, National RI Roadmap Project DO-1153/28.08.2018 and the Academy of Finland grant No. 320045 is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia "Severo Ochoa" SEV-20160588 and SEV-2015-0548, and Unidad de Excelencia "Maria de Maeztu" MDM-2014-0369, by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, the Polish National Research Centre grant UMO-2016/22/M/ST9/00382, and by the Brazilian MCTIC, CNPq and FAPERJ. The Fermi-LAT.Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat a l'Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucleaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK), and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K..A..Wallenberg Foundation, the Swedish Research Council, and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Etudes Spatiales in France. This work made use of data from the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. We thank the NuSTAR Operations, Software, and Calibration teams for support with the execution and analysis of these observations. This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Science Data Center (ASDC; Italy) and the California Institute of Technology (USA). D.P. and A.B. are grateful to Amy Furniss and Wystan Benbow for providing the VERITAS VHE gamma-ray fluxes and for useful discussions about them. M.P. acknowledges support from the Lyman Jr. Spitzer Postdoctoral Fellowship and NASA Fermi grant No. 80NSSC18K1745., and J.F. is partially supported by NASA under contract S-15633Y. M.B. acknowledges support from NASA Headquarters under the NASA Earth and Space Science Fellowship Program (grant NNX14AQ07H) and the black hole Initiative at Harvard University, which is funded in part by the Gordon and Betty Moore Foundation (grant GBMF8273) and in part by the John Templeton Foundation. This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants DN 18-13/2017, DN 18-10/2017, KP06-H28/3 (2018), and KP-06-PN38/1 (2019), as well as for the Spanish MIMECO (AYA2016-80889-P, RYC-2013-14511) and the IAA-CSIC "Severo Ochoa" program SEV2017-0709. The St. Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. The Abastumani team acknowledges financial support of the project FR/638/6-320/12 by the Shota Rustaveli National Science Foundation under contract 31/77. T.G. acknowledges support from Istanbul University (Project numbers 49429 and 48285), Bilim Akademisi (BAGEP program), and TUBITAK (project numbers 13AT100-431, 13AT100-466, and 13AT60-430). The Boston University effort was supported in part by NASA grants NNX12AO90G and NNX14AQ58G. Data from the Steward Observatory spectropolarimetric monitoring project were used in this paper. This program is supported by Fermi Guest Investigator grants NNX08AW56G, NNX09 AU10G, NNX12AO93G, and NNX15AU81G. The OVRO 40 m monitoring program is supported in part by NASA grants NNX08AW31G and NNX11A043G and NSF grants AST-0808050 and AST-1109911. The Metsahovi team acknowledges the support from the Academy of Finland to our observing projects (numbers 212656, 210338, 121148, and others). W.M. acknowledges support from CONICYT project Basal AFB-170002.

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $��$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence., 13 pages, 4 figures, 2 tables, MNRAS accepted version

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability., Comment: 16 pages, 8 figures, 7 tables; Accepted for Publication to AJ

Blazar emission is dominated by non-thermal radiation from a relativistic jet pointing toward us, therefore undergoing Doppler beaming. This is responsible for flux enhancement and contraction of the variability time scales, so that most blazars appear as luminous sources characterized by noticeable and fast flux changes at all frequencies. The mechanisms producing their unpredictable variability are debated and include injection, acceleration and cooling of particles, with possible intervention of shock waves or turbulence. Changes in the viewing angle of the emitting knots or jet regions have also been suggested to explain flaring events or specific properties such as intraday variability, quasi-periodicities, or the delay of radio flux variations relative to optical changes. However, such a geometric interpretation has not been universally accepted because alternative explanations based on changes of physical conditions can also work in many cases. Here we report the results of optical-to-radio monitoring of the blazar CTA 102 by the Whole Earth Blazar Telescope Collaboration and show that the observed long-term flux and spectral variability is best explained by an inhomogeneous, curved jet that undergoes orientation changes. We propose that magnetohydrodynamic instabilities or rotation of a twisted jet cause different jet regions to change their orientation and hence their relative Doppler factors. In particular, the recent extreme optical outburst (six magnitudes) occurred when the corresponding jet emitting region acquired a minimum viewing angle., Comment: 31 pages, 11 figures. This is the accepted version (before Nature editing) of the paper published by Nature, online version at http://dx.doi.org/10.1038/nature24623, Nature 2017/12/04/online

In 2015 July 29 - September 1 the satellite XMM-Newton pointed at the BL Lac object PG 1553+133 six times, collecting data for 218 hours. During one of these epochs, simultaneous observations by the Swift satellite were requested to compare the results of the X-ray and optical-UV instruments. Optical, near-infrared and radio monitoring was carried out by the Whole Earth Blazar Telescope (WEBT) collaboration for the whole observing season. We here present the results of the analysis of all these data, together with an investigation of the source photometric and polarimetric behaviour over the last three years. The 2015 EPIC spectra show slight curvature and the corresponding light curves display fast X-ray variability with a time scale of the order of 1 hour. In contrast to previous results, during the brightest X-ray states detected in 2015 the simple log-parabolic model that best-fits the XMM-Newton data also reproduces reasonably well the whole synchrotron bump, suggesting a peak in the near-UV band. We found evidence of a wide rotation of the polarization angle in 2014, when the polarization degree was variable, but the flux remained almost constant. This is difficult to interpret with deterministic jet emission models, while it can be easily reproduced by assuming some turbulence of the magnetic field., Comment: 15 pages, 12 figures, accepted by Monthly Notices of the Royal Astronomical Society

The astrometric European Space Agency (ESA) Gaia mission was launched in December 19, 2013. One of the tasks of the Gaia mission is production of an astrometric catalog of over one billion stars and more than 500000 extragalactic sources. The quasars (QSOs), as extragalactic sources and radio emitters, are active galactic nuclei objects (AGNs) whose coordinates are well determined via Very Long Baseline Interferometry (VLBI) technique and may reach sub-milliarcsecond accuracy. The QSOs are the defining sources of the quasi-inertial International Celestial Reference Frame (ICRF) because of their core radio morphology, negligible proper motions (until sub-milliarcsecond per year), and apparent point-like nature. Compact AGNs, visible in optical domain, are useful for a direct link of the future Gaia optical reference frame with the most accurate radio one. Apart from the above mentioned activities, Gaia has other goals such as follow-up of transient objects. One of the most important Gaia's requirements for photometric alerts is a fast observation and reduction response, that is, submition of observations within 24 hours. For this reason we have developed a pipeline. In line with possibilities of our new telescope (D(cm)/F(cm)=60/600) at the Astronomical Station Vidojevica (ASV, of the Astronomical Observatory in Belgrade), we joined the Gaia-Follow-Up Network for Transients Objects (Gaia-FUN-TO) for the photometric alerts. Moreover, in view of the cooperation with Bulgarian colleagues (in the frst place, SV), one of us (GD) initiated a local mini-network of Serbian { Bulgarian telescopes useful for the Gaia-FUN-TO and other astronomical purposes. During the next year we expect a new 1.4 m telescope at ASV site. The speed of data processing (from observation to calibration server) could be one day. Here, we present an overview of our activities in the Gaia-FUN-TO which includes establishing Serbian { Bulgarian mini-network (of five telescopes at three sites, ASV in Serbia, Belogradchik and Rozhen in Bulgaria), the Gaia-FUN-TO test observations, and some results. [Projekat Ministarstva nauke Republike Srbije, br. 176011: Dynamics and kinematics of celestial bodies and systems, br. 176004: Stellar physics, i br. 176021: Visible and invisible matter in nearby galaxies: theory and observations]

At the beginning of 2013, we have a team of scientists at the Astronomical Observatory of Belgrade, focusing on observations and research of blazars. At the same time, we established a cooperation with the Astronomical Observatory of Torino (Osservatorio Astrofisico di Torino), center for the Whole Earth Blazar Telescope (WEBT) international collaboration. In this paper, we describe the WEBT project, present our equipment as well as some preliminary results. We also point to some issues that we have encountered during the one year experience and found important from the observational point of view. [Projekat Ministarstva nauke Republike Srbije, br. 176011: Dynamics and Kinematics of Celestial Bodies and Systems, br. 176021: Visible and Invisible Matter in Nearby Galaxies: Theory and Observations i br. 176004: Stellar Physics]

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models., 16 pages, 17 Figures; ApJ accepted

We present our optical photometric observations of three TeV blazars, PKS 1510-089, PG 1553+113 and Mrk 501 taken using two telescopes in India, one in Bulgaria, one in Greece and one in Serbia during 2012 - 2014. These observations covered a total of 95 nights with a total of 202 B filter frames, 247 images in V band, 817 in R band while 229 images were taken in the I filter. This work is focused on multi-band flux and colour variability studies of these blazars on diverse timescales which are useful in understanding the emission mechanisms. We studied the variability characteristics of above three blazars and found all to be active over our entire observational campaigns. We also searched for any correlation between the brightness of the sources and their colour indices. During the times of variability, no significant evidence for the sources to display spectral changes correlated with magnitude was found on timescales of a few months. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux variability., Comment: 11 pages, 7 figures, 9 tables. Accepted for publication in MNRAS

We present a new observational campaign, DWARF, aimed at detection of circumbinary extrasolar planets using the timing of the minima of low-mass eclipsing binaries. The observations will be performed within an extensive network of relatively small to medium-size telescopes with apertures of ~20-200 cm. The starting sample of the objects to be monitored contains (i) low-mass eclipsing binaries with M and K components, (ii) short-period binaries with sdB or sdO component, and (iii) post-common-envelope systems containing a WD, which enable to determine minima with high precision. Since the amplitude of the timing signal increases with the orbital period of an invisible third component, the timescale of project is long, at least 5-10 years. The paper gives simple formulas to estimate suitability of individual eclipsing binaries for the circumbinary planet detection. Intrinsic variability of the binaries (photospheric spots, flares, pulsation etc.) limiting the accuracy of the minima timing is also discussed. The manuscript also describes the best observing strategy and methods to detect cyclic timing variability in the minima times indicating presence of circumbinary planets. First test observation of the selected targets are presented.

Monthly Notices of the Royal Astronomical Society, 454 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

The spectral region in the vicinity of 5394 A contains three prominent photospheric spectral lines, which can be used as a solar plasma diagnostic tool. The occurrence of telluric lines in this region is a potential source of systematic and random errors in these solar spectral lines. The goal of our investigation was to determine the telluric line contamination of this interesting spectral region. Several series of high-resolution solar spectra within an interval of about 4 A around the 5394 A wavelength were observed at different zenith distances of the Sun. Comparison of these spectra has permitted identification of telluric lines in this spectral interval. The observations were carried out with the horizontal solar spectrograph of the Heliophysical Observatory in Debrecen. Telluric feature blending was identified in the blue and red wings of the Fe I 5393.2 A line, and in the local continuum of the Mn I 5394.7 A line. The blue wing of the Fe I 5395.2 A line is contaminated by a weak telluric feature too. The red continuum of this line has a more prominent telluric contamination. A dozen of water vapor telluric lines that determined the observed telluric features were identified in this spectral interval. The profiles of three telluric lines that have a significant influence on both the profiles of solar spectral lines and the level of local continuum were derived, and the variation of their parameters (equivalent width and central depth) with air mass were analyzed.

We use optical color indices (colors) from the SDSS database to study the effect of dust in starburst galaxies by mea-suring the dependence of colors on galaxy inclination. Starburst galaxies with ongoing star formation, are rich with metals/dust and are, therefore, an excellent objects for studying the effect of dust in galaxies. They are selected using the [O III]λ 5007/Hα vs. [N II]λ 6584/Hβ diagram, that is, the BPT-diagram. We use Kauffmann's empirical demarcation line in the BPT-diagram to exclude galaxies with active galactic nuclei (AGN) from the sample because they have different physical and dust properties from normal galaxies. The sample is divided into bins according to galaxy stellar mass and 4000 A break (which is a coarse measure of a galaxy star formation history; SFH) and the reddening with inclination is studied as a function of these two physical parameters. Assuming that the dust effect is negligible in the SDSS z -band, we derive the attenuation curves for these galaxies. We fit the attenuation curves with a simple power law and use power law index to interpret the relative distribution of dust and stars in the starburst galaxies (© 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Attenuation curves (laws) in the optical part of the spectrum (from about 3500 A to 9000 A) are determined for galaxies with different star formation histories (SFHs) using images in five photometrical bands from the Sloan Digital Sky Survey (SDSS). Owing to large surveys like SDSS, it is possible to define several subgroups of galaxies with similar SFHs using two spectral indices that trace SFH. Attenuation curves are analyzed in terms of SFH and compared to the curves that are often used for dust correction of the observed data. .

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events., Comment: 9

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings., Comment: 9 pages, 4 Figures, ApJ Letter accepted

It is a well-established fact that the Mn I 539.47 nm line exhibits significant cycle dependence similar to lines arising in the chromosphere which are affected by non-thermal heating in the chromospheric plages. Among these lines the case of the Mn I 539.47 nm line seems unique. It is of photospheric origin on the basis of theoretical calculations but its cyclic dependence hints at a chromospheric nature. The present work provides further evidence to both connections. The line exhibits asymmetry features and a center-to-limb variation as though influenced by the photospheric granulation. On the other hand an enhancement of the central intensity has been detected in a well identifiable Ca plage area. The line seems to be a promising candidate as an irradiance variation indicator.

We present an analysis of the multiwavelength behaviour of the blazar OJ 248 at z = 0.939 in the period 2006-2013. We use low-energy data (optical, near-infrared, and radio) obtained by 21 observatories participating in the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), as well as data from the Swift (optical-UV and X-rays) and Fermi (gamma-rays) satellites, to study flux and spectral variability and correlations among emissions in different bands. We take into account the effect of absorption by the Damped Lyman Alpha intervening system at z = 0.525. Two major outbursts were observed in 2006-2007 and in 2012-2013 at optical and near-IR wavelengths, while in the high-frequency radio light curves prominent radio outbursts are visible peaking at the end of 2010 and beginning of 2013, revealing a complex radio-optical correlation. Cross-correlation analysis suggests a delay of the optical variations after the gamma-ray ones of about a month, which is a peculiar behaviour in blazars. We also analyse optical polarimetric and spectroscopic data. The average polarization percentage P is less than 3 per cent, but it reaches about 19 per cent during the early stage of the 2012-2013 outburst. A vague correlation of P with brightness is observed. There is no preferred electric vector polarisation angle and during the outburst the linear polarization vector shows wide rotations in both directions, suggesting a complex behaviour or structure of the jet and possible turbulence. The analysis of 140 optical spectra acquired at the Steward Observatory reveals a strong Mg II broad emission line with an essentially stable flux of 6.2 e-15 erg cm-2 s-1 and a full width at half-maximum of 2053 km s-1., Comment: 16 pages, 20 figures, 2 tables, in press for MNRAS

Using the 2 m telescope of the Bulgarian National Astronomical Observatory at Rozhen, observations of 222 double or multiple stars were carried out during three nights in 2011. This is the sixth series of measurements of CCD frames of double and multiple stars obtained at Rozhen. Also in 2011, using the 0.6 m telescope of the Serbian Astronomical Station on the mountain of Vidojevica, observations of 208 double or multiple stars were carried out during six nights. This is the first series of measurements of CCD frames of double and multiple stars obtained at this station. In this paper, we present the results for the position angle and angular separation for 337 pairs and residuals for 72 pairs with published orbital elements or linear solutions. These observations have angular separations in the range from 1.''37 to 172.''81, with a median angular separation of 7.''66. We also present the recalculated linear solutions for four pairs and one linear solution that has been calculated for the first time.