Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.asc.rssi.ru/RadioAstron/publications/articles/aa,583,a100_2015.pdf Äàòà èçìåíåíèÿ: Mon Nov 2 16:03:06 2015 Äàòà èíäåêñèðîâàíèÿ: Sun Apr 10 00:07:11 2016 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: m 87 jet

A&A 583, A100 (2015) DOI: 10.1051/0004-6361/201526335
c ESO 2015

Astronomy & Astrophysics

RadioAstron space VLBI imaging of polarized radio emission in the high-redshift quasar 0642+449 at 1.6 GHz
A. P. Lobanov1,2 , J. L. GÑmez3 , G. Bruni1 , Y. Y. Kovalev4,1 , J. Anderson1,5 , U. Bach1 , A. Kraus1 , J. A. Zensus1 , M. M. Lisakov4 , K. V. Sokolovsky4,6 , and P. A. Voytsik4
1 2 3 4 5 6

Max-Planck-Institut fÝr Radioastronomie, Auf dem HÝgel 69, 53121 Bonn, Germany e-mail: alobanov@mpifr-bonn.mpg.de Institut fÝr Experimentalphysik, UniversitÄt Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany Instituto de AstrofÌsica de AndalucÌa ­ CSIC, Glorieta de la AstronomÌa s/n, 18008 Granada, Spain Astro Space Center of Lebedev Physical Institute, Profsoyuznaya 84/32, 117997 Moscow, Russia Helmholtz-Zentrum Potsdam, Deutsches GeoForschungsZentrum GFZ, Telegrafenberg A6, 14473 Potsdam, Germany Sternberg Astronomical Institute, Moscow State University, Universitetskii pr. 13, 119992 Moscow, Russia

Received 16 April 2015 / Accepted 7 July 2015
ABSTRACT Context. Polarization of radio emission in extragalactic jets at a sub-milliarcsecond angular resolution holds important clues for understanding the structure of the magnetic field in the inner regions of the jets and in close vicinity of the supermassive black holes in the centers of active galaxies. Aims. Space VLBI observations provide a unique tool for polarimetric imaging at a sub-milliarcsecond angular resolution and studying the properties of magnetic field in active galactic nuclei on scales of less than 104 gravitational radii. Methods. A space VLBI observation of high-redshift quasar TXS 0642+449 (OH 471), made at a wavelength of 18 cm (frequency of 1.6 GHz) as part of the early science programme (ESP) of the RadioAstron mission, is used here to test the polarimetric performance of the orbiting Space Radio Telescope (SRT) employed by the mission, to establish a methodology for making full Stokes polarimetry with space VLBI at 1.6 GHz, and to study the polarized emission in the target object on sub-milliarcsecond scales. Results. Polarization leakage of the SRT at 18 cm is found to be within 9% in amplitude, demonstrating the feasibility of high fidelity polarization imaging with RadioAstron at this wavelength. A polarimetric image of 0642+449 with a resolution of 0.8 mas (signifying an 4 times improvement over ground VLBI observations at the same wavelength) is obtained. The image shows a compact corejet structure with low (2%) polarization and predominantly transverse magnetic field in the nuclear region. The VLBI data also uncover a complex structure of the nuclear region, with two prominent features possibly corresponding to the jet base and a strong recollimation shock. The maximum brightness temperature at the jet base can be as high as 4 â 1013 K. Key words. galaxies: jets ­ galaxies: nuclei ­ quasars: individual: 0642+449

1. Introduction
Very long baseline interferometry (VLBI) observations in which one of the antennas is placed on board of an Earth satellite (space VLBI) have provided a capability of reaching unprecedentedly high angular resolution of astronomical observations (Arsentev et al. 1982; Levy et al. 1986; Hirabayashi et al. 2000). The space VLBI mission RadioAstron combines ground-based radio antennas operating at frequencies 0.32, 1.6, 5, and 22 Gigahertz (GHz) with a 10 m antenna (SRT Space Radio Telescope) on board of the satellite Spektr-R launched on 18 July 2011 (Kardashev et al. 2013). At 0.32, 1.6, and 22 GHz, the SRT delivers data in both left (LCP) and right (RCP) circular polarization, which enables fully reconstructing the linearly polarized emission from target objects at these frequencies. At 5 GHz, the SRT provides only LCP data because of a failure of the on-board hardware. After its launch, the SRT was first tested in a singleantenna mode during the in-orbit-checkout (IOC) period (Kovalev et al. 2014) and later in the interferometry mode, in combination with ground-based antenna (fringe search, FS), successfully delivering fringes at each of the four observing bands (Kardashev et al. 2013, 2014a).

The ensuing early science program (ESP) of RadioAstron started in February 2012 and continued through June 2013 (Kardashev et al. 2014a). One of the prime objectives of the program was to provide a bridge between the initial experimental mode of operations, observing, and data processing, and routine operations that started after completion of the ESP. One of the main goals of the ESP was to establish the feasibility of imaging experiments with RadioAstron. The first imaging experiment with RadioAstron, targeting the bright and compact radio source 0716+714, was made during the ESP period at 5 GHz (Kardashev et al. 2014a). Polarization imaging was tested with the observation presented in this paper. Polarization imaging at high angular resolution, both along and transverse to the jet direction, is crucial for measuring the magnetic field structures of jets to distinguish between different possible magnetic field configurations. The improvement of angular resolution provided by space VLBI observations makes it particularly interesting to compare the brightness temperature and polarization measurements made with RadioAstron to the ground VLBI with a similar resolution obtained at higher observing frequencies. It has been shown that a minimum of several, independent, resolution elements across a
A100, page 1 of 10

Article published by EDP Sciences


A&A 583, A100 (2015)

jet are necessary to prevent spurious results from being obtained (Hovatta et al. 2012). Space VLBI observations are essential for resolving transversely a significant number of jets and studying the structure and properties of their magnetic field. Feasibility of polarization measurements with space VLBI has been demonstrated with the VSOP (Hirabayashi et al. 2000) observations made with a satellite antenna recording data only in single (left circular) polarization (Kemball et al. 2000; Gabuzda & GÑmez 2001). Results from the first polarization observation with RadioAstron at a frequency of 1.6 GHz are reported here. This observation was made in March 2013 as part of the ESP and it targeted a bright, compact radio source TXS 0642+449 (OH 471). The radio source 0642+449 is a distant low-polarization quasar (LPQ) located at a redshift of 3.396 (Osmer et al. 1994) which corresponds to the luminosity distance of 30.2 Gigaparsecs (Gpc) and a linear scale of 7.58 kpc/arcsecond, assuming the standard CDM cosmology (Planck Collaboration XIII 2015). The optical spectrum of 0642+449 does not show strong broad lines (Torrealba et al. 2012). In radio, its radio structure has a compact "core­jet" morphology (Gurvits et al. 1992; Xu et al. 1995), with a flat spectral index -0.1 (Hovatta et al. 2014) of the milliarcsecond (mas) scale core (bright narrow end of the jet). The maximum brightness temperature, T b , measured in 0642+449 decreases with observing frequency. Previous VSOP space VLBI observations of 0642+449 at 5 GHz put a lower limit of T b 3.3 â 1012 K (Dodson et al. 2008), while ground VLBI observations at 15 and 86 GHz resolved the compact core and yielded somewhat lower brightness temperatures of 1.2­4.3 â 1012 K (Kovalev et al. 2005) and 1.1­1.6 â 1011 K (Lobanov et al. 2000; Lee et al. 2008). This implies that the compact core becomes gradually resolved out at higher frequencies, calling for comparison with the new space VLBI observations by RadioAstron. Ground VLBI monitoring of the source has revealed a range of proper motions of 0.012­0.078 mas/yr in the jet (Lister et al. 2013) and provided estimates of the jet Lorentz factor, j = 5.4, viewing angle j = 0.8 and intrinsic opening angle j = 0.3 (Pushkarev et al. 2009). With these parameters, an equipartition magnetic field of 0.3 G can be estimated at the location of the jet core (Lobanov 1998) observed at 1.6 GHz, based on measurements of frequency dependent shift of the core position which have been reported for 0642+449 between frequencies of 8.1, 8.4, 12.1, and 15.4 GHz (Pushkarev et al. 2012). Ground VLBI polarization observations at 5.0 and 8.4 GHz (O'Sullivan et al. 2011) have shown a small (1.1­2.6%) degree of polarization in the core, with no conclusive evidence for Faraday rotation ( = 15.6 and 15.3 at 5.0 and 8.4 GHz, respectively). At higher frequencies, between 8.4 and 15.1 GHz, a rotation measure (RM) of -280 rad/m2 has been reported in the core (Hovatta et al. 2012). The Galactic RM in the direction of 0642+449 is estimated to be -5.4 ± 18.3 rad/m2 (Taylor et al. 2009; Oppermann et al. 2012). The polarization observations of 0642+449 with RadioAstron signal the first successful space VLBI polarization experiment using dual circular polarization measurements made at the space-borne antenna. The methodology for RadioAstron polarimetric observations and data correlation is presented in Sect. 2. The RadioAstron observation of 0642+449 and reduction of the polarimetric data are described in Sect. 3. The polarization image of 0642+449 is discussed in Sect. 4 and compared with previous ground VLBI observations of this object.
A100, page 2 of 10

2. Imaging with RadioAstron
RadioAstron observations are performed by combining together ground-based radio telescopes with the orbiting 10 m antenna of the SRT. The data are recorded at a rate of 128 Megabits per second, providing a total bandwidth of up to 32 MHz per circular polarization channel (Andreyanov et al. 2014). The satellite telemetry and science data are received and recorded at tracking stations in Puschino (Andreyanov et al. 2014) and Green Bank (Ford et al. 2014) and later transferred for correlation at one of the RadioAstron correlation facilities. Accurate reconstruction of the orbit and of the momentary state vector of the Spektr-R spacecraft is achieved by combining together radiometric range and radial velocity measurements made at the satellite control stations in Bear Lakes and Ussurijsk (Russia), Doppler measurements obtained at each of the tracking stations, and laser ranging and optical astrometric measurements of the spacecraft position in the sky (Khartov et al. 2014; Zakhvatkin et al. 2014). The accuracy can be further enhanced by simultaneous VLBI measurements performed on the narrow-band downlink signal from the spacecraft (Duev et al. 2012). This was demonstrated (Duev et al. 2015) using the data from the ESP observations of 0642+449 discussed here. The Spektr-R satellite has a highly elliptical orbit, with the apogee reaching up to 350 000 km and the perigee varying in the 7000­80 000 km range. The orbital period of the satellite varies between 8.3 and 9.0 days (Kardashev et al. 2013, 2014b). These factors make imaging observations with RadioAstron particularly challenging, with optimal conditions realized typically only for a limited fraction of sky and for limited periods of time. The imaging observations are planned by simulating the observing conditions in the FakeRaT software package1 (Zhuravlev 2015), which is based on the FakeSat software (Murphy 1991; Murphy et al. 1994). The Fourier spacings provided by the ground-space baselines to the SRT are typically limited to a very narrow range of position angle and often does not reach down to the spacings of the baselines between ground-based telescopes. In order to reduce the adverse effect of these two factors, imaging with RadioAstron is performed during the perigee passages of the SRT (perigee imaging), typically with 12­24 h-long observing segments. Whenever feasible, these perigee imaging segments are augmented by several 1­2 h-long segments (longbaseline tracking) allocated during the following days and including a smaller number (3­4) of ground antennas. The long-baseline tracking provides information on largest Fourier spacings, enabling detecting the emission on the finest angular scales.
2.1. Correlation of RadioAstron experiments

Correlation of RadioAstron experiments is performed at three different correlator facilities: the ASC RadioAstron correlator (Kardashev et al. 2013; Andrianov et al. 2014), the JIVE SFXC correlator (Kettenis 2010), and the DiFX software correlator (Deller et al. 2007, 2011). In order to enable the correlation of RadioAstron data to be made in the DiFX correlator, the DiFX code has been upgraded by the Max-Planck-Institute for Radio Astronomy (MPIfR) in Bonn, addressing the specific aspects of data and telemetry formats of the SRT telescope on-board the Spektr-R spacecraft (Bruni 2014; Bruni et al. 2015).
1

http://www.asc.rssi.ru/radioastron/software/soft. html


A. P. Lobanov et al.: RadioAstron polarization image of 0642+449 Table 1. Radio telescopes participating in RadioAstron observation of 0642+449. Telescope Spektr-R SRT (RU) Effelsberg (DE) Jodrell Bank (UK) WSRT (NL) Torun (PL) Urumqi (CH) Shanghai (CH) Noto (I) Hartebeesthoek (SA) Green Bank (USA) Zelenchukskaya (RU) Code RA EF JB WB TR UR SH NT HH GB ZC D [m] 10 100 76 66 32 25 65 32 26 100 32 SEFD [Jy] 2840 19 65 40 300 300 670 784 430 10 300 Observing GK047A 10:00­01:00 10:00­01:00 10:00­01:00 13:00­01:00 10:00­01:00 10:00­23:00 10:00­19:00 10:30­01:00 14:30­20:00 17:00­01:00 19:30­23:30 time [UT] GK047C 14:00­15:00 14:00­15:00 14:50­15:00 14:00­15:00 Bmax [DE ] ... 5.87 5.93 5.90 2.76 3.02 3.26 2.69 1.91 1.63 1.67

Notes. D ­ antenna diameter (() ­ equivalent diameter); SEFD ­ system equivalent flux density (system noise), describing the antenna sensitivity; Bmax ­ largest projected baseline in the data, in units of the Earth diameter.

The upgrade has comprised a new routine to read the native RadioAstron data format (RDF) into the DiFX correlator and a number of modifications to the correlator code enabling it to process data recorded at the SRT. The delay model server CALC from the CALC/SOLVE package2 used by the DiFX correlator is modified to be able to calculate delay information for telescopes with arbitrary coordinates and velocities, accounting for the general relativistic effects in the gravitational field of the Earth. The DiFX metadata system is upgraded to deal with the changing position and velocity of the spacecraft as a function of time. Satellite orientation parameters provided in the telemetry are used in the upgraded version of the code for calculating the equivalent of parallactic angle correction for the spaceborne antenna. The estimated accuracy of the orientation parameters is 1 (Lisakov et al. 2014). This enables accurate calibration of the instrumental polarization of the SRT. The timestamps for the SRT data (set by the tracking station clock at the moment it begins recording individual VLBI scans) are modified to refer to the relative time of arrival of the astronomical signal at the spacecraft. This is done by calculating the delay for the transmission of the signal from the spacecraft to the tracking station. These corrections are inserted to the metadata handler and to the delay model server, which provides the full compatibility of the RadioAstron data stream with data streams from ground antennas. The upgraded DiFX correlator is now used at the MPIfR to correlate RadioAstron data and being prepared for merging with the trunk version of the DiFX correlator code.

telemetry were recorded by the RadioAstron satellite tracking station in Puschino. Basic parameters of the ground telescopes participating in the segments GK047A and GK047C are given in Table 1. The data were recorded in two polarization channels (left and right circularly polarized, LCP and RCP), with a bandwidth of 32 MHz per channel, centered at a frequency of 1.6601875 GHz. In each of the channels, the total bandwidth was split into two independent intermediate frequency (IF) bands of 16 MHz. The space VLBI segments on 0642+449 were observed in a 40/120 min on/off cycle, with the gaps between the segments required for cooling the motor drive of the on-board high-gain antenna of the Spektr-R. During these gaps, the strong radio sources 1328+207 (3C 286) and 0851+202 (OJ 287) were also observed with the ground antennas as polarization calibrators. Figure 1 shows the resulting visibility coverage of the Fourier domain (uv coverage, expressed in units of spatial frequency defined as a ratio of instantaneous projected baseline length to the observing wavelength, ) of the data recorded for 0642+449.

3.1. Data correlation

3. RadioAstron observation of 0642+449
The 1.6 GHz RadioAstron observation of 0642+449 (RadioAstron project code raes12; global VLBI project code GK047) was made on 9­10 March 2013, with three observing segments allocated for the experiment. The main perigee imaging segment, GK047A (raes12a; 10:00 UT/09/03­ 01:00 UT/10/03) was complemented by VLBI observations of the satellite downlink signal at 8.4 GHz (segment GK047B) made for the purpose of improving the Spektr-R orbit determination (Duev et al. 2015). The third segment, GK047C (raes12b; 14:00 UT/10/03­15:00 UT/10/03) was introduced as a long-baseline tracking segment. The SRT data and the satellite
2

http://gemini.gsfc.nasa.gov/solve/

The data correlation was performed at the DiFX software correlator of the MPIfR. Fringe search for RadioAstron was performed for every scan, and ad-hoc values for delay offsets and rates were determined and applied, in order to compensate for the acceleration terms of the spacecraft. The initial fringe-search was performed with 1024 spectral channels per IF and an integration time of 0.1 s (providing respective search ranges of 10 microseconds and 10-9 s/s for the delay and the delay rate), in order to accommodate for potential large residual delays at the SRT. Fringes between the SRT and the reference antenna of Effelsberg were found with a residual delay offset of only about 1 microsecond or less, and a delay rate of about 2 â 10-11 s/s (60 cm/s). The fringe solutions for the SRT were typically found at high signal-to-noise ratio (SNR): all but three scans showed fringes with SNR > 10. For the three scans with low SNR, values of delay offset and rate were interpolated from adjacent scans. Final correlation was performed with 32 spectral channels per IF (with the corresponding channel width of 500 kHz) and 0.5 s of integration time.
A100, page 3 of 10


A&A 583, A100 (2015)

Fig. 2. Total delay rates of the SRT obtained from baseline-based fringe fitting of the correlated data on the baseline to Effelsberg and accounting for the clock rates at each of the two telescopes. For each SRT observing segment, the rates are obtained from fringe fitting solutions with a time interval of 10 min. The rates are plotted separately for each IF and polarization channels. The plot reflects the accuracy of the velocity determination in the reconstructed orbit of the SRT. All of the velocity and acceleration values correspond to the baseline projections of the true residual offset, velocity, and acceleration of the spacecraft.

The segments GK047A and GK047C were calibrated independently and combined together for the final imaging. Polarization calibration was performed using the segment GK047A (the parallactic angle coverage of GK047C was insufficient for using this segment for the polarization calibration). The a priori amplitude calibration was applied using the default antenna gains and the system temperature measurements made at each antenna during the observation. For the SRT, the sensitivity parameters measured during the IOC (Kovalev et al. 2014) were used. The data were edited using the flagging information from the station log data. Parallactic angle correction was applied, to correct for feed rotation with respect to the target sources.
3.2.1. Fringe fitting

Fig. 1. Coverage of the Fourier domain (uv coverage) of the RadioAstron observation of 0642+449 at 1.6 GHz ( = 18 cm), plotted in units of M. The central part corresponds to the perigee imaging segment (GK047A) and the short segment at v 400 M represents the long-baseline tracking segment (GK047C). Color marks the upper limit of brightness temperature obtained from visibility amplitudes (for details, see discussion in Sect. 4.1).

The data were fringe fitted in two steps, first by applying manual phase-cal corrections and then using the global fringe search for antenna rates and single and multi-band delays. The group delay difference between the two polarization channels was corrected. The resulting post-fringe residual delay rates obtained for the SRT are plotted in Fig. 2. These slowly evolving residuals agree well with the expected accuracy of the orbital velocity determination for the SRT (Kardashev et al. 2013), but they still should be viewed only as an indicator of the fidelity of the fringe fitted data, while a more detailed discussion of the accuracy of the orbit determination for the SRT is presented elsewhere (Duev et al. 2015).
3.2.2. Bandpass calibration

3.2. Post-correlation data reduction

The correlated data were reduced in several steps using AIPS3 and Difmap (Shepherd 1997, 2011). The initial calibration was performed in AIPS, imaging was done using Difmap, and both packages were used for calibrating the instrumental polarization and providing absolute calibration of the polarization vectors.
3

Astronomical Image Processing Software of the National Radio Astronomy Observatory, USA. A100, page 4 of 10

After the fringe fitting, the receiver bandpasses were corrected. The bandpass-calibrated data on the baselines between the SRT and the ground antennas in Effelsberg and Urumqi is shown in Fig. 3. The interferometric visibility signal is clearly detected both in the total intensity (LL and RR) and the cross-polarization (LR and RL) channels, demonstrating the excellent quality of the data. Table 2 lists typical amplitude and phase errors of the


A. P. Lobanov et al.: RadioAstron polarization image of 0642+449 Table 2. Typical amplitude and phase errors on the baselines to the SRT. Ant. Code EF JB WB TR UR SH NT HH GB ZC Polarization LL,RR LR,RL amp ph amp ph [%] [ ] [%] [ ] 5.5 2.9 22.6 50.2 9.3 3.0 30.7 51.9 6.0 3.1 30.5 46.8 11.6 4.6 38.3 54.2 13.6 6.5 45.0 60.5 53.9 56.6 53.9 84.6 18.3 9.3 54.1 80.0 52.8 65.1 50.8 67.2 6.5 2.6 33.4 42.7 25.1 13.8 50.6 57.7

Notes. Columns represent average fractional amplitude, amp and phase, ph , errors of the total (LL, RR) and cross-polarized (LR, RL) correlated flux density detected on the SRT baselines to the ground radio telescopes.

Fig. 3. Bandpass amplitudes and phases on ground-space baselines between the SRT and Effelsberg (top two rows) and Urumqi (bottom two rows), after applying global fringe fitting solutions. The bandpasses for total intensity (parallel hands: LL, RR) and polarized emission (crosshands: LR, RL) are shown.

visibility data on the longest baselines to the SRT observed in the segment GK047A, demonstrating that the source was detected on baselines to all participating antennas (the errors were calculated by vector averaging the calibrated data over two IFs and over a scan length of 10 min, and hence these errors are conservative estimates containing also contributions from residual bandpasses and offsets between the IFs, and from amplitude and phase variations over the scan time).
3.2.3. Polarization calibration

After applying the bandpass calibration, the calibrated Stokes I data were exported to Difmap and imaged using hybrid imaging with CLEAN deconvolution and self-calibration (with full phase self-calibration and limited amplitude calibration using a single amplitude gain correction per antenna). The resulting selfcalibrated visibility amplitudes and phases (and their CLEAN model representations) are shown in Fig. 4 for all of the baselines as a function of Fourier spacing (uv distance). The visibilities indicate that the source structure is clearly detected up to the longest ground-space baselines of the observation. The polarization calibration proceeded then with the selfcalibrated model of the source structure obtained. The model was imported into AIPS and used as the calibration model for the full Stokes data, with amplitude and phase calibration applied on time intervals of 1 min.

The instrumental polarization (polarization leakage terms or D-terms, comprising the fractional amplitude and phase) was obtained using the LPCAL method (Leppanen et al. 1995), with Effelsberg used as the reference antenna. The resulting polarization leakages determined through this procedure are listed in Table 3. Errors of the D-terms were estimated from averaging the respective values determined for each of the calibrators and for the target itself. To estimate the errors for the SRT (which did not observe the calibrators), the D-terms obtained from each of the calibrators with the ground antennas were applied to the data on 0642+449 and respective solutions for the D-terms of the SRT were obtained and compared with the solution obtained directly from the target data (hence the D-terms errors for the SRT might be underestimated). The amplitude of instrumental polarization of the SRT is found to be within 9% and with excellent phase consistency between the IFs. This figure cannot be directly compared with the results from pre-launch measurements, as those measurements were made only for the antenna feeds of the SRT, yielding an 2% leakage at 1.6 GHz (Turygin 2014). The polarization leakage of the SRT is comparable to the instrumental polarization obtained for the ground antennas, demonstrating robust polarization performance of the SRT and ensuring reliable polarization measurements with data on all ground-space baselines. Similar estimates of the L-band polarization leakage were obtained for the SRT from the analysis of the RadioAstron AGN survey observations (Pashchenko et al. 2015). After applying the D-terms determined for the SRT and the ground antennas, the calibrated full Stokes data were exported into Difmap, imaged in Stokes I , Q, and U , and combined together, providing images of the total intensity, linearly polarized intensity, and polarization vectors (electric vector position angle, EVPA).
3.2.4. EVPA calibration

The absolute EVPA calibration was performed using single-dish measurements of the total flux density polarization of the target and calibrator sources. These measurements (listed in Table 4) were made concurrently with the RadioAstron observation at the Effelsberg 100-m telescope. The averaged EVPA of the polarized emission detected in the VLBI data was related to the polarization vectors that were measured with the Effelsberg antenna,
A100, page 5 of 10


A&A 583, A100 (2015)

GK047A GK047C

SRT baselines

Fig. 4. Visibility amplitude (top) and phase (bottom) distributions as a function of uv radius, overplotted with the CLEAN model (red color) obtained during the hybrid imaging of the source structure. Table 3. Instrumental polarization (D-terms). Antenna RA EF JB WB TR UR SH NT HH GB ZC 6. 7. 2. 3. 1. 0. 4. 4. 8. 8. 11. 12. 6. 6. 5. 5. 21. 18. 5. 4. 8. 9. 8 0 6 1 3 8 6 6 0 7 6 5 3 3 8 3 7 1 4 8 0 5 RCP m [%] [ ] ± 0.3 82 ± 6 ± 0.2 81 ± 8 ± 0.3 10 ± 2 ± 0.3 8±3 ± 0.4 140 ± 7 ± 0.4 306 ± 7 ± 0.5 176 ± 3 ± 0.9 182 ± 3 ± 0.5 160 ± 3 ± 0.4 164 ± 3 ± 1.0 271 ± 12 ± 1.1 294 ± 12 ± 0.5 173 ± 31 ± 1.0 170 ± 32 ± 0.3 247 ± 5 ± 0.2 184 ± 5 ± 1.8 151 ± 23 ± 1.8 198 ± 23 ± 0.3 248 ± 6 ± 0.3 236 ± 6 ± 0.7 160 ± 7 ± 0.7 146 ± 7 LCP m [%] [ ] ± 0.8 98 ± 1 ± 0.9 98 ± 2 ± 0.1 167 ± 9 ± 0.1 185 ± 9 ± 0.6 336 ± 3 ± 0.6 330 ± 3 ± 0.2 166 ± 13 ± 0.4 199 ± 12 ± 0.6 14 ± 8 ± 0.3 31 ± 8 ± 1.1 56 ± 18 ± 1.2 93 ± 18 ± 0.4 13 ± 30 ± 0.8 349 ± 30 ± 0.2 239 ± 6 ± 0.2 178 ± 8 ± 5.2 80 ± 35 ± 5.2 149 ± 35 ± 0.8 295 ± 10 ± 0.8 315 ± 10 ± 1.0 146 ± 3 ± 1.0 141 ± 3 Table 4. Effelsberg measurements of total and polarized flux densities of the targets. Source 0642+449 0851+202 1328+307 S tot [Jy] 1.32 ± 0.02 2.77 ± 0.03 13.68 ± 0.21 S pol [Jy] 0.022 ± 0.008 0.069 ± 0.016 1.37 ± 0.10 [deg] 77.2 ± 12.0 32.1 ± 1.7 33.3 ± 1.4

8. 8. 2. 2. 3. 4. 2. 2. 7. 7. 10. 12. 5. 4. 5. 5. 19. 9. 5. 3. 5. 7.

2 7 3 4 5 3 5 4 2 6 2 2 5 5 9 4 9 5 1 5 2 2

Notes. Measurements were made on 9­10 March 2013. Column designation: S tot ­ total flux density; S pol ­ polarized flux density; m ­ percentage of polarization; ­ position angle of the electric vector position angle (EVPA).

and the resulting EVPA offsets were applied to the VLBI images. The resulting ground array images of the calibrator sources are shown in Fig. 5, and the image parameters are listed in Table 6.

4. Polarization image of 0642+449
The RadioAstron polarization image of 0642+449 is shown in Fig. 6, and the basic parameters of the image are listed in Table 6. The source shows a compact core-dominated structure, with a resolved nuclear region which contains about 98% of the total flux density in the image. The weak and extended jet feature, J1, located at about 4 mas separation from the nucleus has a peak brightness of 15 mJy/beam, which is not surprising considering energy losses due to adiabatic expansion of the flow up to the estimated deprojected linear separation of 2.3 kpc of this feature from the core of the source (assuming the jet viewing angle of 0.8 ; Pushkarev et al. 2009). The linear separation can

Notes. For each antenna, listed are the fractional amplitude, m, and phase, , of the instrumental polarization (polarization leakage) in the right (RCP) and left (LCP) circular polarization channel and in the first (IF1, top row) and second (IF2, bottom row) intermediate frequency channel. The D-term phases have been corrected for the phase offsets obtained from the EVPA calibration described in Sect. 3.2.4. A100, page 6 of 10


A. P. Lobanov et al.: RadioAstron polarization image of 0642+449

Fig. 5. Ground and 0851+202 parameters are colors indicate orientation.

array images of the calibrator sources 1328+307 (left) (right), made using the uniform data weighting. Image listed in Table 6. Contours show the total intensity; linearly polarized intensity; vectors show the EVPA

be a factor of 3.5 smaller, if the observed misalignment of 2 between the position angles of C2 and J1 reflects the physical change of the jet direction. The observed concentration of the polarized flux density along the southern edge of this feature may reflect interaction with the ambient medium occurring in this region, although it may also result from the errors of the D-terms determination or from other instrumental effects (cf. Gabuzda & GÑmez 2001; Hovatta et al. 2012). The core region is clearly resolved and extended along a PA of 81 , It can be decomposed into two circular Gaussian components. The resulting three-component decomposition of the source structure is presented in Table 5. The two nuclear components, C1 and C2, are separated by 0.76 mas (410 pc, deprojected). Similar structural composition of the nuclear region is also detected in the ground VLBI from the MOJAVE survey at 15 GHz (Lister et al. 2013). The same structure is also present in the images of 0642+449 made at 24 and 43 GHz (Fey et al. 2000). Overall, the structure observed in 0642+449 may be similar to the main features of the jet in M 87 (cf. Sparks et al. 1996), in which a prominent standing shock (component HST-1) is located at 200­280 pc (considering the uncertainty of the jet viewing angle) and a prominent feature A located at 2.8­3.7 kpc distance from the core is believed either to be an oblique shock (Bicknell & Begelman 1996) or to result from the helical mode of KelvinHelmholtz instability (Lobanov et al. 2003). This may also explain the apparent displacement between the total intensity and linearly polarized emission in the region J1, which is similar to the polarization asymmetries observed in M 87 at these linear scales (Owen et al. 1989). The polarized emission detected in 0642+449 is dominated by the two nuclear components, with fractional polarization reaching 15.6 mJy (2%), and the magnetic field showing predominantly transverse orientation with respect to the general jet direction (assuming that the innermost component C1 corresponds to the canonical "VLBI core" region in which the

emission becomes optically thin at the frequency of observation; Konigl 1981). The polarization vectors vary by about 10 across the nucleus, with the polarization angle observed in the most compact region C1 (80 ) coming close to the polarization angle measured with Effelsberg. A substantial rotation of the EVPA vectors is seen in some of the MOJAVE images of 0642+449, which is likely to result from plasma propagation in a jet observed at an extremely small viewing angle. The weak (2 mJy/beam) linearly polarized emission in the jet component J1 corresponds to a fractional polarization of 15%, with a large error margin on this estimate. The magnetic field is roughly longitudinal in this region. The MOJAVE database4 lists, for the frequency of 15 GHz, an EVPA of 122 in March 2011. Analysis of the multifrequency MOJAVE observations (Hovatta et al. 2012) yields, for 0642+449, a median rotation measure of -268 rad/m2 over the whole source structure, with a slightly higher rotation measure of -280 rad/m2 over the compact core region. Assuming a linear trend in the EVPA in time as derived from the two latest MOJAVE observations, the polarization angle should be 105 in March 2013 at 15 GHz. We measure = 80 ­90 over the core region. Formally, this corresponds to a rotation measure of -15 rad/m2 , under assumption that the rotation measure is not affected by the core shift. It would then require about three turns of phase to get this value to be reconciled with the RM reported from the MOJAVE measurements. Alternatively, the smaller RM estimated from the RadioAstron measurements may reflect the reported trend of decreasing RM with decreasing frequency (Kravchenko et al. 2015). It would still require about one turn of phase to bring the estimated RM value to a better agreement with the empirical dependence |RM|() 1.6 derived in Kravchenko et al. (2015).
4.1. Brightness temperature

The Gaussian modelfit described in Table 5 can be used for obtaining estimates of brightness temperature (also presented in the same Table), with the highest brightness temperature of 1.3 â 1012 K associated with the most compact core feature C1. This value is somewhat lower than the one that should be expected at 1.6 GHz, based on the trend log10 T b [K] = 13.3 - 1.1 log10 [GHz] that can be obtained from the existing measurements at 5, 15, and 86 GHz (Lobanov et al. 2000; Kovalev et al. 2005; Dodson et al. 2008; Lee et al. 2008). However, since the morphology of the emitting region can be different from a simple Gaussian shape, the actual brightness temperature of the emission may differ from the modelfit-based estimates, for instance, if the brightness distribution is smooth and the respective brightness temperature is largely determined mostly by the transverse dimension of the flow (Lobanov 2015). To investigate this possibility, the modelfit-based estimates of T b can be compared with estimates obtained directly from visibility amplitudes and their errors, which yields the absolute minimum brightness temperature, T b,min , and an estimate of the limiting brightness temperature, T b,lim that can obtained from the data under the requirements that the structural detail sampled by the given visibility is resolved (Lobanov 2015). These two estimates are compared in Fig. 7 with the brightness temperatures calculated from the Gaussian modelfit. This comparison shows that the visibility amplitudes on the space baselines longer that 140 M require the brightness temperature to be larger than 1012 K, which agrees with T b obtained
4

http://www.physics.purdue.edu/MOJAVE/ A100, page 7 of 10


A&A 583, A100 (2015) Table 5. Decomposition of compact structure in 0642+449. Comp. C1 C2 J1 S tot [mJy] 462 ± 16 721 ± 19 20 ± 9 r [mas] 0.38 ± 0.01 0.32 ± 0.01 3.80 ± 0.20 [ ] -88.3 ± 0.7 70.1 ± 1.1 92.0 ± 3.2 [mas] 0.40 ± 0.01 0.60 ± 0.01 0.89 ± 0.43 min [mas] 0.16 0.13 0.75 Tb [1011 K] 12.8 ± 0.6 8.8 ± 0.3 0.11 ± 0.04 S pol [mJy] 7.7 ± 0.8 9.9 ± 0.7 1.7 ± 0.7 mpol [%] 1.7 ± 0.2 1.4 ± 0.1 8.5 ± 5.2 pol [ ] 79 ± 3 89 ± 2 39 ± 11

Notes. Gaussian model description: S tot ­ total flux density; component position in polar coordinates (r, ) with respect to the map center; component size, ; minimum resolvable size, min , of the respective component; brightness temperature, T b , in the observer's frame, derived for the parameters of the Gaussian fit. Polarization properties (measured from the Stokes U and Q images: S pol ­ flux density of linearly polarized emission; mpol ­ fractional polarization; pol ­ polarization position angle in the sky plane (measured north-through-east).

C2 C1 J1

10 pc

Fig. 6. RadioAstron space VLBI image of 0642+449 of total (contours) and linearly polarized (color shades) emission, made using the uniform data weighting. The orientation of linear polarization is indicated by vectors. Basic parameters of the total intensity image are listed in Table 6. Labels indicate the modelfit decomposition of the source structure described in Table 5. Table 6. Parameters of total intensity images. Source 0851+202 1328+307 0642+449 S tot 2.06 7.87 1.21 S peak 1.73 0.89 0.81 S neg ­19.0 ­47.1 3.6 rms 6.5 8.6 0.7 Beam 12.2, 2.93, 12.6 7.48, 2.73, 14.9 1.38, 0.81, 66.5

Notes. Column designation: S tot [Jy] ­ total flux density; S peak [Jy/beam] ­ peak flux density; S neg [mJy/beam] ­ maximum negative flux density in the image; rms [mJy/beam] ­ rms noise in the image; Beam: major axis, minor axis, position angle of major axis [mas, mas, ].

for the most compact feature C1 of the Gaussian modelfit. The limiting brightness temperature T b,lim = 3.6 â 1012 K can be estimated at the longest space baselines. These visibilities concentrate however only within a narrow range of position angles (-10 2 ± 0 1), with corresponding structural sensitivity lim. . ited to PA 80 , which essentially coincides with direction of the inner jet (hence resulting in lower visibility amplitudes and, consequently, lower estimates of T b,lim ).
A100, page 8 of 10

Indeed, some of the T b,lim estimated from the visibility amplitudes at shorter space baselines and different position angles reach above 1014 K, notably with a large spread of the estimates obtained at similar spatial frequencies. The observed spread reflects the difference of signal-to-noise ratios (SNR) in the measurements on different baselines, and it also may indicate potential biases due to a priori gain calibration and amplitude selfcalibration. The three particularly outstanding clusters of points in Fig. 7, which have T b,lim > 1014 K, come from visibilities on the baselines to Westerbork, Torun, Noto, and Hartebeesthoek. As the visibility SNR on these baselines is typically lower than on baselines to Effelsberg or Jodrell Bank (which do not yield such high T b,lim at similar radial distances and position angles), it is most likely that the exceptionally high brightness temperature limits are caused by problems with the amplitude calibration of the antenna gains or specific observing scans with these antennas. If the outliers are ignored and only the bulk of the T b,lim estimates are considered in the range of 100­200 M, the average T b,lim increases to 9.1 â 1012 K, and this can be viewed as a viable upper limit on the brightness temperature of the most


A. P. Lobanov et al.: RadioAstron polarization image of 0642+449

15 Tb 14 Tb 13 Tb 12 Tb 11 10 9 8
,min ,lim ,lim

2.4 0.4 1014 K 3.6 1.8 1012 K

log10 Tb K

,mod

1.3 1012 K 1.2 0.1 1012 K

Space baselines

0

100

200 q M

300

400

Fig. 7. Visibility-based estimates of the brightness temperature in 0642+449, calculated in the observer's frame. Data points represent the minimum (dark gray), T b,min , and maximum limiting (light gray), T b,lim , brightness temperatures derived from individual visibilities measured at different uv distances, q (Lobanov 2015). The two-dimensional distribution of T b,lim is also presented in Fig. 1 with a color wedge. Braces above and below the brightness temperature data indicate the range of uv distances covered by the space baselines of the RadioAstron observation. The histograms show the respective brightness temperatures averaged over bins of 5 M in size. The thick dashed line indicates the brightness temperature, T b,mod = 1.3 â 1012 K, derived for the most compact feature, C1, of the Gaussian fit of the source structure described in Table 5. The dotted lines indicate the respective brightness temperatures of the other two Gaussian components. The minimum brightness temperature, T b,min = 1.2 â 1012 K, is constrained by the visibilities with q > 380 M. The maximum limiting brightness temperature T b,lim = 2.4 â 1014 K can be estimated from the data on the SRT baselines to Noto and Hartebeesthoek, although potential problems with amplitude calibration cannot be ruled out. A more conservative estimate of T b,lim = 3.6 â 1012 K can be obtained from the data on the longest space-ground baselines.

compact structure in the jet of 0642+449. This corresponds to brightness temperature of 4.0 â 1013 K in the rest frame of the source. If the brightness temperature is generally determined by the transverse dimension of the flow, as suggested by the analysis of the MOJAVE visibility data (Lobanov 2015), the inner jet of 0642+449 should have a width of 0.15 mas to satisfy the limiting brightness temperature estimated above.

resolution of ground VLBI images the physical conditions in the jet as 0.2 mas (linear scale of 1.5 SNR-driven resolution limits of the

at 18 cm. The image probes on angular scales as small pc) taking into account the observation.

5. Summary
The RadioAstron Early Science observations of the high-redshift quasar TXS 0642+449 have tested the polarization performance of the Space Radio Telescope (SRT) at 1.6 GHz ( = 18 cm) and provided a framework for establishing a set of procedures for correlation, post-processing and calibration of space VLBI polarization data. The instrumental polarization of the SRT is found to be less than 9%, which enables robust reconstruction of the polarization signal on the ground-space baselines of RadioAstron observations and provides a basis for high-fidelity imaging of polarized emission. The total and linearly polarized emission from the quasar 0642+449 is detected on the SRT baselines of up to 75 560 km (5.93 Earth diameters) in length, corresponding to angular scales of 0.5 milliarcseconds (mas). The resulting RadioAstron image of linearly polarized emission in 0642+449 has a resolution of 0.8 mas, which is 4 times better than the

The structure of total intensity and linearly polarized emission recovered from the RadioAstron image of 0642+449 suggests that this object is likely to be a "canonical" quasar with a powerful relativistic jet observed at an extremely small angle of sight and featuring a bright "core" of the jet (component C1) and a likely recollimation shock (component C2) located at a deprojected distance of about 400 pc away from the core. Each of these regions is only weakly polarized, with fractional polarization not exceeding 2% and polarization vector implying a predominantly transverse magnetic field (in agreement with expectations for strong shocks in jet). The weak and more distant feature J1 observed at a separation of 3.8 mas (2 kpc, deprojected) has a much stronger polarization, possibly exceeding a 10% level. The dominant magnetic field components is likely to be poloidal in this region, implying that it either represents a weak shock or may even result from plasma instability developing in the jet. The RadioAstron data have also been used for making brightness temperature measurements based on the Gaussian decomposition (modelfit) of the structure and on estimates made directly from the visibility data. The most robust estimate
A100, page 9 of 10


A&A 583, A100 (2015)

provided by these measurements yields a brightness temperature of 3.6 â 1012 K for the most compact region in the jet. Further analysis of the RadioAstron data indicate that the the brightness temperature of the radio emission from this region cannot be lower that 1012 K and is not likely to exceed 9 â 1012 K (corresponding to 4 â 1013 K in the reference frame of the source).
Acknowledgements. The RadioAstron project is led by the Astro Space Center of the Lebedev Physical Institute of the Russian Academy of Sciences and the Lavochkin Scientific and Production Association under a contract with the Russian Federal Space Agency, in collaboration with partner organizations in Russia and other countries. This research is based on observations correlated at the Bonn Correlator, jointly operated by the Max Planck Institute for Radio Astronomy (MPIfR), and the Federal Agency for Cartography and Geodesy (BKG). Y.Y.K., M.M.L., K.V.S., P.A.V. are supported by the Russian Foundation for Basic Research (RFBR) grant 13-02-12103. K.V.S. is also supported by the RFBR grant 14-02-31789. J.L.G. acknowledges support from the Spanish Ministry of Economy and Competitiveness grant AYA2013-40825P. The European VLBI Network is a joint facility of European, Chinese, South African and other radio astronomy institutes funded by their national research councils. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

References
Andreyanov, V. V., Kardashev, N. S., & Khartov, V. V. 2014, Cosmic Res., 52, 319 Andrianov, A. S., Girin, I. A., Zharov, V. E., et al. 2014, Proc. of S.A. Lavochkin Association, 3, 55 Arsentev, V. M., Berzhatyi, V. I., Blagov, V. D., et al. 1982, Akademiia Nauk SSSR Doklady, 264, 588 Bicknell, G. V., & Begelman, M. C. 1996, ApJ, 467, 597 Bruni, G. 2014, in COSPAR Meeting, Vol. 40, 40th COSPAR Scientific Assembly, Held 2­10 August 2014, in Moscow, Russia, Abstract E1.10-2214., 418 Bruni, G., Anderson, J. M., Alef, W., Lobanov, A. P., & Zensus, J. A. 2015, in Proc. 12th European VLBI Network Symp., PoS(EVN 2014)119 Deller, A. T., Tingay, S. J., Bailes, M., & West, C. 2007, PASP, 119, 318 Deller, A. T., Brisken, W. F., Phillips, C. J., et al. 2011, PASP, 123, 275 Dodson, R., Fomalont, E. B., Wiik, K., et al. 2008, ApJS, 175, 314 Duev, D. A., Molera CalvÈs, G., Pogrebenko, S. V., et al. 2012, A&A, 541, A43 Duev, D. A., Zakhvatkin, M. V., Stepanyants, V. A., et al. 2015, A&A, 573, A99 Fey, A. L., Boboltz, D. A., Gaume, R. A., & Johnston, K. J. 2000, in International VLBI Service for Geodesy and Astrometry, General Meeting Proc., eds. N. R. Vandenberg, & K. D. Baver, 285 Ford, H. A., Anderson, R., Belousov, K., et al. 2014, in Proc. SPIE, 9145, 91450 Gabuzda, D. C., & GÑmez, J. L. 2001, MNRAS, 320, L49 Gurvits, L. I., Kardashev, N. S., Popov, M. V., et al. 1992, A&A, 260, 82 Hirabayashi, H., Hirosawa, H., Kobayashi, H., et al. 2000, PASJ, 52, 955 Hovatta, T., Lister, M. L., Aller, M. F., et al. 2012, AJ, 144, 105 Hovatta, T., Aller, M. F., Aller, H. D., et al. 2014, AJ, 147, 143

Kardashev, N. S., Khartov, V. V., Abramov, V. V., et al. 2013, Astron. Rep., 57, 153 Kardashev, N. S., Alakoz, A. V., Kovalev, Y. Y., et al. 2014a, Proc. S.A. Lavochkin Association, 3, 4 Kardashev, N. S., Kreisman, B. B., Pogodin, A. V., et al. 2014b, Cosmic Res., 52, 332 Kemball, A., Flatters, C., Gabuzda, D., et al. 2000, PASJ, 52, 1055 Kettenis, M. 2010, in 10th European VLBI Network Symposium and EVN Users Meeting: VLBI and the New Generation of Radio Arrays, 86 Khartov, V. V., Shirshakov, A. E., Artyukhov, M. I., et al. 2014, Cosmic Res., 52, 326 Konigl, A. 1981, ApJ, 243, 700 Kovalev, Y. Y., Kellermann, K. I., Lister, M. L., et al. 2005, AJ, 130, 2473 Kovalev, Y. A., Vasil'kov, V. I., Popov, M. V., et al. 2014, Cosmic Res., 52, 393 Kravchenko, E. V., Cotton, W. D., & Kovalev, Y. Y. 2015, in IAU Symp. 313, eds. F. Massaro, C. C. Cheung, E. Lopez, & A. Siemiginowska, 128 Lee, S.-S., Lobanov, A. P., Krichbaum, T. P., et al. 2008, AJ, 136, 159 Leppanen, K. J., Zensus, J. A., & Diamond, P. J. 1995, AJ, 110, 2479 Levy, G. S., Linfield, R. P., Ulvestad, J. S., et al. 1986, Science, 234, 187 Lisakov, M. M., Voinakov, S. M., Syrov, A. S., et al. 2014, Cosmic Res., 52, 365 Lister, M. L., Aller, M. F., Aller, H. D., et al. 2013, AJ, 146, 120 Lobanov, A. P. 1998, A&A, 330, 79 Lobanov, A. P. 2015, A&A, 574, A84 Lobanov, A. P., Krichbaum, T. P., Graham, D. A., et al. 2000, A&A, 364, 391 Lobanov, A. P., Hardee, P., & Eilek, J. 2003, New Astron. Rev., 47, 629 Murphy, D. W. 1991, in Radio Interferometry. Theory, Techniques, and Applications, IAU Colloq. 131, eds. T. J. Cornwell, & R. A. Perley, ASP Conf. Ser., 19, 107 Murphy, D. W., Yakimov, V., Kobayashi, H., Taylor, A. R., & Fejes, I. 1994, in VLBI TecnologyY: Progress and Future Observational Possibilities, eds. T. Sasao, S. Manabe, O. Kameya, & M. Inoue, 34 Oppermann, N., Junklewitz, H., Robbers, G., et al. 2012, A&A, 542, A93 Osmer, P. S., Porter, A. C., & Green, R. F. 1994, ApJ, 436, 678 O'Sullivan, S. P., Gabuzda, D. C., & Gurvits, L. I. 2011, MNRAS, 415, 3049 Owen, F. N., Hardee, P. E., & Cornwell, T. J. 1989, ApJ, 340, 698 Pashchenko, I. N., Kovalev, Y. Y., & Voitsik, P. A. 2015, Cosmic Res., 53, 199 Planck Collaboration XIII. 2015, A&A, submitted [arXiv:1502.01589] Pushkarev, A. B., Kovalev, Y. Y., Lister, M. L., & Savolainen, T. 2009, A&A, 507, L33 Pushkarev, A. B., Hovatta, T., Kovalev, Y. Y., et al. 2012, A&A, 545, A113 Shepherd, M. C. 1997, in Astronomical Data Analysis Software and Systems VI, eds. G. Hunt & H. Payne, ASP Conf. Ser., 125, 77 Shepherd, M. 2011, Astrophysics Source Code Library [record ascl:1103.001] Sparks, W. B., Biretta, J. A., & Macchetto, F. 1996, ApJ, 473, 254 Taylor, A. R., Stil, J. M., & Sunstrum, C. 2009, ApJ, 702, 1230 Torrealba, J., Chavushyan, V., Cruz-GonzÀlez, I., et al. 2012, Rev. Mex. Astron. Astrofis., 48, 9 Turygin, M. S. 2014, Cosmic Res., 52, 403 Xu, W., Readhead, A. C. S., Pearson, T. J., Polatidis, A. G., & Wilkinson, P. N. 1995, ApJS, 99, 297 Zakhvatkin, M. V., Ponomarev, Y. N., Stepan'yants, V. A., Tuchin, A. G., & Zaslavskiy, G. S. 2014, Cosmic Res., 52, 342 Zhuravlev, V. I. 2015, Cosmic Res., 53, 216

A100, page 10 of 10