Context. We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Aims. A summary of the contents of Gaia EDR3 is presented, accompanied by a discussion on the differences with respect to Gaia DR2 and an overview of the main limitations which are present in the survey. Recommendations are made on the responsible use of Gaia EDR3 results. Methods. The raw data collected with the Gaia instruments during the first 34 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium and turned into this early third data release, which represents a major advance with respect to Gaia DR2 in terms of astrometric and photometric precision, accuracy, and homogeneity. Results. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the (GBP − GRP) colour are also available. The passbands for G, GBP, and GRP are provided as part of the release. For ease of use, the 7 million radial velocities from Gaia DR2 are included in this release, after the removal of a small number of spurious values. New radial velocities will appear as part of Gaia DR3. Finally, Gaia EDR3 represents an updated materialisation of the celestial reference frame (CRF) in the optical, the Gaia-CRF3, which is based solely on extragalactic sources. The creation of the source list for Gaia EDR3 includes enhancements that make it more robust with respect to high proper motion stars, and the disturbing effects of spurious and partially resolved sources. The source list is largely the same as that for Gaia DR2, but it does feature new sources and there are some notable changes. The source list will not change for Gaia DR3. Conclusions. Gaia EDR3 represents a significant advance over Gaia DR2, with parallax precisions increased by 30 per cent, proper motion precisions increased by a factor of 2, and the systematic errors in the astrometry suppressed by 30–40% for the parallaxes and by a factor ~2.5 for the proper motions. The photometry also features increased precision, but above all much better homogeneity across colour, magnitude, and celestial position. A single passband for G, GBP, and GRP is valid over the entire magnitude and colour range, with no systematics above the 1% level

Since the launch of the Fermi satellite, BL Lacertae has been moderately active at γrays and optical frequencies until May 2011, when the source started a series of strong flares. The exceptional optical sampling achieved by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) in collaboration with the Steward Observatory allows us to perform a detailed comparison with the daily γ-ray observations by Fermi. Discrete correlation analysis between the optical and γ-ray emission reveals correlation with a time lag of 0 ± 1 d, which suggests cospatiality of the corresponding jet emitting regions. A better definition of the time lag is hindered by the daily gaps in the sampling of the extremely fast flux variations. In general, optical flares present more structure and develop on longer time scales than corresponding γ-ray flares. Observations at X-rays and at millimetre wavelengths reveal a common trend, which suggests that the region producing the mm and X-ray radiation is located downstream from the optical and γ-ray-emitting zone in the jet. The mean optical degree of polarisation slightly decreases over the considered period and in general it is higher when the flux is lower. The optical electric vector polarisation angle (EVPA) shows a preferred orientation of about 15 • , nearly aligned with the radio core EVPA and mean jet direction. Oscillations around it increase during the 2011-2012 outburst. We investigate the effects of a geometrical interpretation of the long-term flux variability on the polarisation. A helical magnetic field model predicts an evolution of the mean polarisation that is in reasonable agreement with the observations. These can be fully explained by introducing slight variations in the compression factor in a transverse shock waves model.

Blazars are active galactic nuclei, which are powerful sources of radiation whose central engine is located in the core of the host galaxy. Blazar emission is dominated by non-thermal radiation from a jet that moves relativistically towards us, and therefore undergoes Doppler beaming. This beaming causes flux enhancement and contraction of the variability timescales, so that most blazars appear as luminous sources characterized by noticeable and fast changes in brightness at all frequencies. The mechanism that produces this unpredictable variability is under debate, but proposed mechanisms include injection, acceleration and cooling of particles, with possible intervention of shock waves or turbulence. Changes in the viewing angle of the observed emitting knots or jet regions have also been suggested as an explanation of flaring events and can also explain specific properties of blazar emission, such as intra-day variability, quasi-periodicity and the delay of radio flux variations relative to optical changes. Such a geometric interpretation, however, is not universally accepted because alternative explanations based on changes in physical conditions-such as the size and speed of the emitting zone, the magnetic field, the number of emitting particles and their energy distribution-can explain snapshots of the spectral behaviour of blazars in many cases. Here we report the results of optical-to-radio-wavelength monitoring of the blazar CTA 102 and show that the observed long-term trends of the flux and spectral variability are best explained by an inhomogeneous, curved jet that undergoes changes in orientation over time. We propose that magnetohydrodynamic instabilities or rotation of the twisted jet cause different jet regions to change their orientation and hence their relative Doppler factors. In particular, the extreme optical outburst of 2016-2017 (brightness increase of six magnitudes) occurred when the corresponding emitting region had a small viewing angle. The agreement between observations and theoretical predictions can be seen as further validation of the relativistic beaming theory.

We present coordinated multiwavelength observations of the bright, nearby BL Lacertae object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very high energy(VHE) γ-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep 3 G » power law, with no evidence for an exponential cutoff or additional hard components up to ∼80keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure that relates to the two bumps of the broadband SED. In each bump, the variability increases with energy, which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in situ electron acceleration and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

Context. The sudden optical brightening of two young stellar objects, HBC 722 and VSX J205126.1+440523, located in the North America/Pelican Nebula Complex, was announced in August 2010. Early photometric and spectroscopic observations of these objects indicated that they may belong to the FUor or EXor class of young eruptive stars. The eruptions of FUors and EXors are often explained by enhanced accretion of material from the circumstellar disk to the protostar. Aims. In order to determine the true nature of these two objects, we started an optical and near-infrared monitoring program, and complemented our data with archival observations and data from the literature.Methods. We plot and analyze pre-outburst and outburst spectral energy distributions (SEDs), multi-filter light curves, and colorcolor diagrams. Results. The quiescent SED of HBC 722 is consistent with that of a slightly reddened normal T Tauri-type star. The source brightened monotonically in about two months, and the SED obtained during maximum brightness indicates the appearance of a hot, singletemperature blackbody. The current fading rate implies that the star will return to quiescence in about a year, which questions its classification as a bone fide FUor. The quiescent SED of VSX J205126.1+440523 looks like that of a highly embedded Class I source. The outburst of this source happened more gradually, but reached an unprecedentedly high amplitude. Its light curves showed a deep minimum two and a half months after the peak, when the object was close to its pre-outburst optical brightness. Further monitoring indicates that it is still far from being quiescent. Conclusions. The shape of the light curves as well as the bolometric luminosities and accretion rates suggest that these objects do not fit into the classic FUor group. Although HBC 722 exhibits all spectral characteristics of a bona fide FUor, its luminosity and accretion rate is too low and its timescale is too fast compared with classical FUors. VSX J205126.1+440523 seems to be an example where quick extinction changes modulate the light curve.

Aims. Amongst more than fifty blazars detected in very high energy (VHE, E > 100 GeV) γ rays, only three belong to the subclass of flat spectrum radio quasars (FSRQs). The detection of FSRQs in the VHE range is challenging, mainly because of their soft spectra in the GeV-TeV regime. MAGIC observed PKS 1510−089 (z = 0.36) starting 2012 February 3 until April 3 during a high activity state in the high energy (HE, E > 100 MeV) γ-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 standard deviations (σ). We study the multi-frequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO, and VLBA telescopes), X-ray (Swift satellite), and HE γ-ray frequencies. Methods. We study the VHE γ-ray emission, together with the multi-frequency light curves, 43 GHz radio maps, and spectral energy distribution (SED) of the source. The quasi-simultaneous multi-frequency SED from the millimetre radio band to VHE γ rays is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infrared torus and a slow sheath surrounding the jet around the Very Long Baseline Array (VLBA) core. Results. We find that the VHE γ-ray emission detected from PKS 1510−089 in 2012 February-April agrees with the previous VHE observations of the source from 2009 March-April. We find no statistically significant variability during the MAGIC observations on daily, weekly, or monthly time scales, while the other two known VHE FSRQs (3C 279 and PKS 1222+216) have shown daily scale to sub-hour variability. The γ-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multi-frequency light curves suggest a common origin for the millimetre radio and HE γ-ray emission, and the HE γ-ray flaring starts when the new component is ejected from the 43 GHz VLBA core and the studied SED models fit the data well. However, the fast HE γ-ray variability requires that within the modelled large emitting region, more compact regions must exist. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

Context. This work is part of the Gaia Data Processing and Analysis Consortium papers published with the Gaia Early Data Release 3 (EDR3). It is one of the demonstration papers aiming to highlight the improvements and quality of the newly published data by applying them to a scientific case. Aims. We use the Gaia EDR3 data to study the structure and kinematics of the Magellanic Clouds. The large distance to the Clouds is a challenge for the Gaia astrometry. The Clouds lie at the very limits of the usability of the Gaia data, which makes the Clouds an excellent case study for evaluating the quality and properties of the Gaia data. Methods. The basis of our work are two samples selected to provide a representation as clean as possible of the stars of the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC). The selection used criteria based on position, parallax, and proper motions to remove foreground contamination from the Milky Way, and allowed the separation of the stars of both Clouds. From these two samples we defined a series of subsamples based on cuts in the colour-magnitude diagram; these subsamples were used to select stars in a common evolutionary phase and can also be used as approximate proxies of a selection by age. Results. We compared the Gaia Data Release 2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasible with the use of additional external data. We derive radial and tangential velocity maps and global profiles for the LMC for the several subsamples we defined. To our knowledge, this is the first time that the two planar components of the ordered and random motions are derived for multiple stellar evolutionary phases in a galactic disc outside the Milky Way, showing the differences between younger and older phases. We also analyse the spatial structure and motions in the central region, the bar, and the disc, providing new insightsinto features and kinematics. Finally, we show that the Gaia EDR3 data allows clearly resolving the Magellanic Bridge, and we trace the density and velocity flow of the stars from the SMC towards the LMC not only globally, but also separately for young and evolved populations. This allows us to confirm an evolved population in the Bridge that is slightly shift from the younger population. Additionally, we were able to study the outskirts of both Magellanic Clouds, in which we detected some well-known features and indications of new ones.

Aims. We aim to characterize the multiwavelength emission from Markarian 501 (Mrk 501), quantify the energy-dependent variability, study the potential multiband correlations, and describe the temporal evolution of the broadband emission within leptonic theoretical scenarios. Methods. We organized a multiwavelength campaign 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. Results. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of ∼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 ∼3 CU, and the peak of the high-energy spectral component was found to be at ∼2 TeV. Both the X-ray and VHE gamma-ray spectral slopes were measured to be extremely hard, with spectral indices < 2 during most of the observing campaign, regardless of the X-ray and VHE flux. This study reports the hardest Mrk 501 VHE spectra measured to date. The fractional variability was found to increase with energy, with the highest variability occurring at VHE. Using the complete data set, we found correlation between the X-ray and VHE bands; however, if the June 9 flare is excluded, the correlation disappears (significance < 3σ) despite the existence of substantial variability in the X-ray and VHE bands throughout the campaign. Conclusions. 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 data set acquired shows that the broadband spectral energy distribution (SED) of Mrk 501, and its transient evolution, is very complex, requiring, within the framework of synchrotron self-Compton (SSC) models, various emission regions for a satisfactory description. Nevertheless the one-zone 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 behavior 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.