skip to main content

Title: Multiband optical variability of the TeV blazar PG 1553 + 113 in 2019
ABSTRACT

We report the flux and spectral variability of PG 1553 + 113 on intra-night (IDV) to short-term time-scales using BVRI data collected over 91 nights from 28 February to 8 November 2019 employing 10 optical telescopes: three in Bulgaria, two each in India and Serbia, and one each in Greece, Georgia, and Latvia. We monitored the blazar quasi-simultaneously for 16 nights in the V and R bands and 8 nights in the V, R, I bands and examined the light curves (LCs) for intra-day flux and colour variations using two powerful tests: the power-enhanced F-test and the nested ANOVA test. The source was found to be significantly (>99 per cent) variable in 4 nights out of 27 in R-band, 1 out of 16 in V-band, and 1 out of 6 nights in I-band. No temporal variations in the colours were observed on IDV time-scale. During the course of these observations the total variation in R-band was 0.89 mag observed. We also investigated the spectral energy distribution (SED) using B-, V-, R-, and I-band data. We found optical spectral indices in the range of 0.878 ± 0.029 to 1.106 ± 0.065 by fitting a power law (Fν∝ν−α) to these SEDs of PG 1553 + 113. We found that the more » source follows a bluer-when-brighter trend on IDV time-scales. We discuss possible physical causes of the observed spectral variability.

« less
Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; « less
Publication Date:
NSF-PAR ID:
10390203
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
519
Issue:
2
Page Range or eLocation-ID:
p. 2796-2811
ISSN:
0035-8711
Publisher:
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Resonant tunneling diodes (RTDs) have come full-circle in the past 10 years after their demonstration in the early 1990s as the fastest room-temperature semiconductor oscillator, displaying experimental results up to 712 GHz and fmax values exceeding 1.0 THz [1]. Now the RTD is once again the preeminent electronic oscillator above 1.0 THz and is being implemented as a coherent source [2] and a self-oscillating mixer [3], amongst other applications. This paper concerns RTD electroluminescence – an effect that has been studied very little in the past 30+ years of RTD development, and not at room temperature. We present experiments and modeling of an n-type In0.53Ga0.47As/AlAs double-barrier RTD operating as a cross-gap light emitter at ~300K. The MBE-growth stack is shown in Fig. 1(a). A 15-μm-diam-mesa device was defined by standard planar processing including a top annular ohmic contact with a 5-μm-diam pinhole in the center to couple out enough of the internal emission for accurate free-space power measurements [4]. The emission spectra have the behavior displayed in Fig. 1(b), parameterized by bias voltage (VB). The long wavelength emission edge is at  = 1684 nm - close to the In0.53Ga0.47As bandgap energy of Ug ≈ 0.75 eV at 300 K.more »The spectral peaks for VB = 2.8 and 3.0 V both occur around  = 1550 nm (h = 0.75 eV), so blue-shifted relative to the peak of the “ideal”, bulk InGaAs emission spectrum shown in Fig. 1(b) [5]. These results are consistent with the model displayed in Fig. 1(c), whereby the broad emission peak is attributed to the radiative recombination between electrons accumulated on the emitter side, and holes generated on the emitter side by interband tunneling with current density Jinter. The blue-shifted main peak is attributed to the quantum-size effect on the emitter side, which creates a radiative recombination rate RN,2 comparable to the band-edge cross-gap rate RN,1. Further support for this model is provided by the shorter wavelength and weaker emission peak shown in Fig. 1(b) around = 1148 nm. Our quantum mechanical calculations attribute this to radiative recombination RR,3 in the RTD quantum well between the electron ground-state level E1,e, and the hole level E1,h. To further test the model and estimate quantum efficiencies, we conducted optical power measurements using a large-area Ge photodiode located ≈3 mm away from the RTD pinhole, and having spectral response between 800 and 1800 nm with a peak responsivity of ≈0.85 A/W at  =1550 nm. Simultaneous I-V and L-V plots were obtained and are plotted in Fig. 2(a) with positive bias on the top contact (emitter on the bottom). The I-V curve displays a pronounced NDR region having a current peak-to-valley current ratio of 10.7 (typical for In0.53Ga0.47As RTDs). The external quantum efficiency (EQE) was calculated from EQE = e∙IP/(∙IE∙h) where IP is the photodiode dc current and IE the RTD current. The plot of EQE is shown in Fig. 2(b) where we see a very rapid rise with VB, but a maximum value (at VB= 3.0 V) of only ≈2×10-5. To extract the internal quantum efficiency (IQE), we use the expression EQE= c ∙i ∙r ≡ c∙IQE where ci, and r are the optical-coupling, electrical-injection, and radiative recombination efficiencies, respectively [6]. Our separate optical calculations yield c≈3.4×10-4 (limited primarily by the small pinhole) from which we obtain the curve of IQE plotted in Fig. 2(b) (right-hand scale). The maximum value of IQE (again at VB = 3.0 V) is 6.0%. From the implicit definition of IQE in terms of i and r given above, and the fact that the recombination efficiency in In0.53Ga0.47As is likely limited by Auger scattering, this result for IQE suggests that i might be significantly high. To estimate i, we have used the experimental total current of Fig. 2(a), the Kane two-band model of interband tunneling [7] computed in conjunction with a solution to Poisson’s equation across the entire structure, and a rate-equation model of Auger recombination on the emitter side [6] assuming a free-electron density of 2×1018 cm3. We focus on the high-bias regime above VB = 2.5 V of Fig. 2(a) where most of the interband tunneling should occur in the depletion region on the collector side [Jinter,2 in Fig. 1(c)]. And because of the high-quality of the InGaAs/AlAs heterostructure (very few traps or deep levels), most of the holes should reach the emitter side by some combination of drift, diffusion, and tunneling through the valence-band double barriers (Type-I offset) between InGaAs and AlAs. The computed interband current density Jinter is shown in Fig. 3(a) along with the total current density Jtot. At the maximum Jinter (at VB=3.0 V) of 7.4×102 A/cm2, we get i = Jinter/Jtot = 0.18, which is surprisingly high considering there is no p-type doping in the device. When combined with the Auger-limited r of 0.41 and c ≈ 3.4×10-4, we find a model value of IQE = 7.4% in good agreement with experiment. This leads to the model values for EQE plotted in Fig. 2(b) - also in good agreement with experiment. Finally, we address the high Jinter and consider a possible universal nature of the light-emission mechanism. Fig. 3(b) shows the tunneling probability T according to the Kane two-band model in the three materials, In0.53Ga0.47As, GaAs, and GaN, following our observation of a similar electroluminescence mechanism in GaN/AlN RTDs (due to strong polarization field of wurtzite structures) [8]. The expression is Tinter = (2/9)∙exp[(-2 ∙Ug 2 ∙me)/(2h∙P∙E)], where Ug is the bandgap energy, P is the valence-to-conduction-band momentum matrix element, and E is the electric field. Values for the highest calculated internal E fields for the InGaAs and GaN are also shown, indicating that Tinter in those structures approaches values of ~10-5. As shown, a GaAs RTD would require an internal field of ~6×105 V/cm, which is rarely realized in standard GaAs RTDs, perhaps explaining why there have been few if any reports of room-temperature electroluminescence in the GaAs devices. [1] E.R. Brown,et al., Appl. Phys. Lett., vol. 58, 2291, 1991. [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [2] M. Feiginov et al., Appl. Phys. Lett., 99, 233506, 2011. [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [3] Y. Nishida et al., Nature Sci. Reports, 9, 18125, 2019. [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [4] P. Fakhimi, et al., 2019 DRC Conference Digest. [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018).« less
  2. ABSTRACT

    In 2021 BL Lacertae underwent an extraordinary activity phase, which was intensively followed by the Whole Earth Blazar Telescope (WEBT) Collaboration. We present the WEBT optical data in the BVRI bands acquired at 36 observatories around the world. In mid-2021 the source showed its historical maximum, with R = 11.14. The light curves display many episodes of intraday variability, whose amplitude increases with source brightness, in agreement with a geometrical interpretation of the long-term flux behaviour. This is also supported by the long-term spectral variability, with an almost achromatic trend with brightness. In contrast, short-term variations are found to be strongly chromatic and are ascribed to energetic processes in the jet. We also analyse the optical polarimetric behaviour, finding evidence of a strong correlation between the intrinsic fast variations in flux density and those in polarization degree, with a time delay of about 13 h. This suggests a common physical origin. The overall behaviour of the source can be interpreted as the result of two mechanisms: variability on time-scales greater than several days is likely produced by orientation effects, while either shock waves propagating in the jet, or magnetic reconnection, possibly induced by kink instabilities in the jet, can explain variability onmore »shorter time-scales. The latter scenario could also account for the appearance of quasi-periodic oscillations, with periods from a few days to a few hours, during outbursts, when the jet is more closely aligned with our line of sight and the time-scales are shortened by relativistic effects.

    « less
  3. Abstract In order to constrain the size of the optical continuum emission region in the dwarf Seyfert 1 galaxy NGC 4395 through reverberation mapping, we carried out high-cadence photometric monitoring in the griz filter bands on two consecutive nights in 2022 April using the four-channel MuSCAT3 camera on the Faulkes Telescope North at Haleakalā Observatory. Correlated variability across the griz bands is clearly detected, and the r -, i -, and z -band light curves show lags of 7.72 − 1.09 + 1.01 , 14.16 − 1.25 + 1.22 , and 20.78 − 2.09 + 1.99 minutes with respect to the g band when measured using the full-duration light curves. When lags are measured for each night separately, the Night 2 data exhibit lower cross-correlation amplitudes and shorter lags than the Night 1 light curves. Using the full-duration lags, we find that the lag–wavelength relationship is consistent with the τ ∝ λ 4/3 dependence found for more luminous active galactic nuclei. Combining our results with continuum lags measured for other objects, the lag between g and z band scales with optical continuum luminosity as τ gz ∝ L 0.56±0.05 , similar to the scaling of broad-line region size with luminosity,more »reinforcing recent evidence that diffuse continuum emission from the broad-line region may contribute substantially to optical continuum variability and reverberation lags.« less
  4. MeerKAT’s large number (64) of 13.5 m diameter antennas, spanning 8 km with a densely packed 1 km core, create a powerful instrument for wide-area surveys, with high sensitivity over a wide range of angular scales. The MeerKAT Galaxy Cluster Legacy Survey (MGCLS) is a programme of long-track MeerKAT L -band (900−1670 MHz) observations of 115 galaxy clusters, observed for ∼6−10 h each in full polarisation. The first legacy product data release (DR1), made available with this paper, includes the MeerKAT visibilities, basic image cubes at ∼8″ resolution, and enhanced spectral and polarisation image cubes at ∼8″ and 15″ resolutions. Typical sensitivities for the full-resolution MGCLS image products range from ∼3−5 μJy beam −1 . The basic cubes are full-field and span 2° × 2°. The enhanced products consist of the inner 1.2° × 1.2° field of view, corrected for the primary beam. The survey is fully sensitive to structures up to ∼10′ scales, and the wide bandwidth allows spectral and Faraday rotation mapping. Relatively narrow frequency channels (209 kHz) are also used to provide H  I mapping in windows of 0 <  z  < 0.09 and 0.19 <  z  < 0.48. In this paper, we provide an overview of the survey and the DR1 products, including caveatsmore »for usage. We present some initial results from the survey, both for their intrinsic scientific value and to highlight the capabilities for further exploration with these data. These include a primary-beam-corrected compact source catalogue of ∼626 000 sources for the full survey and an optical and infrared cross-matched catalogue for compact sources in the primary-beam-corrected areas of Abell 209 and Abell S295. We examine dust unbiased star-formation rates as a function of cluster-centric radius in Abell 209, extending out to 3.5 R 200 . We find no dependence of the star-formation rate on distance from the cluster centre, and we observe a small excess of the radio-to-100 μm flux ratio towards the centre of Abell 209 that may reflect a ram pressure enhancement in the denser environment. We detect diffuse cluster radio emission in 62 of the surveyed systems and present a catalogue of the 99 diffuse cluster emission structures, of which 56 are new. These include mini-halos, halos, relics, and other diffuse structures for which no suitable characterisation currently exists. We highlight some of the radio galaxies that challenge current paradigms, such as trident-shaped structures, jets that remain well collimated far beyond their bending radius, and filamentary features linked to radio galaxies that likely illuminate magnetic flux tubes in the intracluster medium. We also present early results from the H  I analysis of four clusters, which show a wide variety of H  I mass distributions that reflect both sensitivity and intrinsic cluster effects, and the serendipitous discovery of a group in the foreground of Abell 3365.« less
  5. Context. Stellar evolution models are highly dependent on accurate mass estimates, especially for highly massive stars in the early stages of stellar evolution. The most direct method for obtaining model-independent stellar masses is derivation from the orbit of close binaries. Aims. Our aim was to derive the first astrometric plus radial velocity orbit solution for the single-lined spectroscopic binary star MWC 166 A, based on near-infrared interferometry over multiple epochs and ∼100 archival radial velocity measurements, and to derive fundamental stellar parameters from this orbit. A supplementary aim was to model the circumstellar activity in the system from K band spectral lines. Methods. The data used include interferometric observations from the VLTI instruments GRAVITY and PIONIER, as well as the MIRC-X instrument at the CHARA Array. We geometrically modelled the dust continuum to derive relative astrometry at 13 epochs, determine the orbital elements, and constrain individual stellar parameters at five different age estimates. We used the continuum models as a base to examine differential phases, visibilities, and closure phases over the Br γ and He  I emission lines in order to characterise the nature of the circumstellar emission. Results. Our orbit solution suggests a period of P  = 367.7 ± 0.1 d, approximatelymore »twice as long as found with previous radial velocity orbit fits. We derive a semi-major axis of 2.61 ± 0.04 au at d  = 990 ± 50 pc, an eccentricity of 0.498 ± 0.001, and an orbital inclination of 53.6 ± 0.3°. This allowed the component masses to be constrained to M 1  = 12.2 ± 2.2  M ⊙ and M 2  = 4.9 ± 0.5  M ⊙ . The line-emitting gas was found to be localised around the primary and is spatially resolved on scales of ∼11 stellar radii, where the spatial displacement between the line wings is consistent with a rotating disc. Conclusions. The large spatial extent and stable rotation axis orientation measured for the Br γ and He  I line emission are inconsistent with an origin in magnetospheric accretion or boundary-layer accretion, but indicate an ionised inner gas disc around this Herbig Be star. We observe line variability that could be explained either with generic line variability in a Herbig star disc or V/R variations in a decretion disc scenario. We have also constrained the age of the system, with relative flux ratios suggesting an age of ∼(7 ± 2)×10 5 yr, consistent with the system being composed of a main-sequence primary and a secondary still contracting towards the main-sequence stage.« less