arXiv:2005.09122v1 [astro-ph.IM] 18 May 2020 ino ua itne rdcn mgn eouinat resolution imaging producing distance, with lunar GHz frac- substantial a 22 a of and or km, tion 4.8, 240,000 to 1.7, up at baselines radio projected fringes space VLBI arrays meter VLBI useful terrestrial 10 found with conjunction a in used with telescope mission, RadioAstron The referencing al- phase [4]. another, done same-beam to been of has source amount one limited from a i slew telescopes though rapidly large to of space inability wide-angle the free produce to to due of intended data, None astrometric was fringes. missions and produced three [2] three acquired (VSOP) these all have HALCA and data [1], [3], to TDRS RadioAstron missions space, in three data these, VLBI Of antennas VLBI space. of in deployment led the produc- has for which baselines proposals fringes), longest numerous (i.e., the correlations of significant length frequency ing the given by a at limited observations are VLBI of resolution The Introduction 1 the in ICRF the of future. stabilty near and relatively accuracy the in im- two-orders-of-magnitude SMBH provement a equipped of to observation lead properly could the A horizons to event devoted future. mission VLBI forseeable space the for CRF any anchored be to CRF rectly resulting “emis- astrometri the as enable SMBH rings) will these photon references associated of the Use (and rings” sion smaller shadow. is black-hole region emitting the the than that showing mm, 1.3 at M87 RIGS 00 oe tl,2 uut-5Spebr2020 September 5 - August 29 Italy, Rome, 2020, GASS URSI a bevdte40 the observed a has array, (EHT), (VLBI) Telescope Interferometry Baseline Horizon Long Very mm-wave Event Frames the Reference Recently, Celestial resulting (CRFs). the sy and in differ- instabilities errors at to tematic different leading appear frequencies, will observing and ent inher change are to jets These subject quasars. ently power- and (SMBH) pro- nuclei Holes galactic jets Black active Massive ing of Frame Super observations the Reference by from duced Celestial derived Frame Gaia are Reference optical (Gaia-CRF2) the Celestial and International (ICRF) radio the Both nhrdi hdw:TigteClsilRfrneFaeDi Frame Reference Celestial the Tying Shadows: in Anchored nSB hdw;teutmt eeec onsfor points reference ultimate the shadows; SMBH in µ sdaee hdwo h MHin SMBH the of shadow diameter as Abstract [email protected] amBy lrd,USA Florida, Bay, Palm pc ntaie Inc Initiatives Space .M Eubanks M. T. di- to s- n c - , aeegho . m(3 H)adhdarslto of resolution a had and GHz) M87 (230 ∼ of mm image 1.3 EHT of first wavelength The a used. the frequency of the diameter and the a Earth by or limited EHT, be will existing EHT, the augmented of LEO existing resolution the the However, over [7]. sensitivity EHT and range substantial dynamic a in provide would increase the (LEO) Orbit of limiting Earth extension Low an sites, to Even EHT high-altitude coverage. imaging and plane dry (u,v) EHT the very the limits mostly which using and to used frequencies at- which with the bands, caused both possible sub-mm and limits limitations is mm the the what in to absorption of mospheric due limits VLBI, the terrestrial pushing purely is EHT The [6]. (EHT) Telescope of Horizon shadow Event the by re- of M87 imaging been successful has the mission frequencies. by VLBI provided space higher cently new even a to at impetus telescopes Further with extended radio and space continued be larger to needs clearly work This ore a xii rgtestmeaue 10 > radio temperatures natural brightness that exhibit showing can and level, sources second arc micro the aeiea htwvlnt.Fgr hw htms of of most masses that have shows which 1 known, Figure AGN the wavelength. that at baseline evto foesuc n nosraino nte.Any another. of observation an and source ob- an one for of delay servation baseline between difference d the signal on astrometric pending astrometr the doubly-differential, all is that effectively data of and arrival VLBI telescope) of time another the at minus arriva telescope signal of one time at the signal is a delay of (the differential is data VLBI All me and delays. position propagation telescope dia the in of errors fraction in- errors, thus timing a clude delays within Non-astrometric to wavelength. ideally observing system, the of apparent VLBI time or the the delays in at non-geometric delays all geometry also network and the observation both of knowledge modeling requires or astrometry Accurate and therein). [8] e.g., quan- erences (see, these sources in different change from the obtained tities) and observ- delay VLBI phase and the (group comparing ables on based is astrometry VLBI planet. of the use from removed the well require telescopes will VLBI resolution substan- EHT A in Earth. unless improvement the EHT to tial current close very the cosmologically by are resolved they be not will masses) lar 20 µ s ls otefig pcn fa1Erhdiameter Earth 1 a of spacing fringe the to close as, etyt lc oeEetHorizons Event Hole Black to rectly . 10 10 13 M J [5]. K (so- ref- e- ic - l 100000 non-geometric delays that do not change between different 11 Angular Size for 10 MSun 10 observations, or whose change can be adequately modeled, Angular Size for 10 MSun 10000 9 can be removed from the astrometric solution using doubly Angular Size for 10 MSun M87 (as measured) differential data. 1000

In most global (or wide-field) VLBI astrometry, such as the 100 data used to produce the various radio CRFs, observations 10 EHT are not simultaneous, but are as close together in time as possible using fast-slewing telescopes. Also, the geometri- Angular Size (microas) 1 cal changes in the network, such as Earth rotation changes Earth−Moon VLBI and continental drift, are typically slow in an Earth-fixed 0.1 ELS−2 VLBI frame and can be estimated from external models or from 0.01 the VLBI data itself. However, space VLBI telescopes typ- 0.0001 0.001 0.01 0.1 1 10 ically have very low slew rates and are of course in orbit, Red shift (z) moving rapidly with respect to the Earth, and their orbits Figure 1. The angular size as a function of red shift and frequently have unmodeled or hard to model perturbations. SMBH mass for a reasonable range of AGN masses, af- The fast-slewing doubly differential techniques of terres- ter [11], together the resolution of the EHT and Earth- trial astrometry cannot be applied. Moon and Earth-Sun Lagrange Point 2 (ESL-2) baselines. TON618, the heaviest SMBH known, has a mass of 6.6 10 11 ∼ × 2 Multi-Beam Multi-Antenna VLBI 10 solar masses, and so the upper curve, at 10 MJ, rep- resents a optimistic upper limit for the sizes of any SMBH Although there is a long history of highly accurate same- shadows in the universe [12]. Note that space baselines will beam VLBI [9], with standard radio telescopes this tech- be required to resolve most AGN SMBH shadows out to nique can only be applied in the unusual cases where two or cosmological distances. more sources are within the same telescope beam (i.e., have separations of a few arc-minutes at most). The Japanese VERA (VLBI Exploration of Radio Astrometry) project is can cover a total beamwidth of ε radians. For a given pri- mary source, this secondary antenna can cover an array of multi-beam VLBI antennas used mostly at 22 and 48 GHz for differential for phase referencing extra- A = 2 π ε sinθ steradians (1) galactic sources and astronomical water and silicon oxide masers. This system has four receiver bands, at 2, 8, 22, of the sky. For θ 90 degrees and ε 1 degree, for any and 43 GHz and can observe different objects with maxi- given primary source∼ the secondary antenna∼ will be able to 3 mum 2.2 degree separation [10], which allows for a reason- reach secondary sources in about 4 10− of the sky, and able probability of finding a phase reference source near the thus, even given only a few thousand× suitable sources, one desired target. Because the two beams are close together or more simultaneous source pairs would be generally ac- in the sky, and the differential VLBI is done with pairs of cessible for any given primary source. Each observation sources with the same clock, with very close paths through would thus determine the relative positions of 2 widely- the atmosphere, and similar instrumental and antenna de- separated sources; by repeating this with many pairs of lays, VERA, with an accuracy goal of 10 µ as, is capable sources a wide-field CRF could be built up over time. A ∼ of highly accurate differential astrometry and shows what multibeam space telescope or telescopes could be used to be done with narrow-. This dual-beam technique could also construct a space-Earth or space-space CRF), providing a be used to extend the sensitivity limit of space VLBI; phase substantial increase in both global astrometric accuracy and changes observed with a strong source in one beam could differential phase VLBI. be used to correct the phase errors on a much weaker source in the adjacent beam. 3 Anchoring on Rings and Shadows

In addition to a multiple-beam approach, a multiple- The black hole event horizon, of course, basically emits no telescope approach seems promising, where 2 antennas are photons; SMBH shadows are observed through the photons mounted together at a fixed or nearly fixed relative angle. emitted from near, but outside the even horizon. While The space multiple-telescope observing mode would be to SMBH shadows may or may not be suitable “standard point one antenna at a source, while the entire assemblage is rulers” for cosmology [11, 12], the photon emission region rotated about the vector pointing at the first source to point will tied to the mass and spin rate of the host SMBH and the second antenna at a second source. This would pro- will make it possible (in many cases) to parametrically fit vide for instantaneous mutual observations, with the same for the shadow centroid. Figure 1 shows that for many clock and spacecraft position, between two widely sepa- SMBH at cosmological distances the shadow size will be rated sources without require any rapid slewing at any time. .1 µas.

Suppose that the second antenna rotatesat an angle θ and General Relativity predicts that inside the bright emission 10000 [15] This effort has in fact already begun; the ICRF-3 in- 300 km / s (typical Peculiar Velocity) 6000 km / s (extreme merger kick) cludes a model for the acceleration caused by the Galaxy, 1000 i.e., the aberration rate caused by the rotation of the So- )

−1 100 lar System Barycenter (SSB) about the center of the Milky Way galaxy, with the VLBI estimates of the aberration am- 10 1 plitude being in the range of 5.1 to 6.4 µas yr− [16]; this effect is near the accuracy limit of the current VLBI CRF. 1 Near−Earth VLBI

0.1 Figure 2 shows that with a network extending to the ESL- Proper Motion (microas yr Earth−Moon VLBI 2 Lagrange point the galactic aberration would be two or- 0.01 ELS−2 VLBI ders of magnitude above the noise and the expected galac- tic peculiar velocities (typically order 300 km s 1) could 0.001 − 0.0001 0.001 0.01 0.1 1 10 be detected over much of the recent cosmological past. The Red shift (z) “kicks” resulting from equal mass SMBH binary mergers Figure 2. Kinematic proper motions for the velocity range can be significantly larger and, although rare, should be eas- expected for non-binary SMBH, together with expected ily detectable [17]. limits from one decade of observing with near-Earth, Earth- Moon and ESL-2 baselines, respectively. For VLBI with 5 Conclusions ESL-2 baselines, kinematic proper motions should be ob- servable out to cosmological distances. A space-Earth or space-space multibeam multiantenna VLBI network would tie the Celestial Reference System to SMBH shadows, providing a substantial increase in the ring will be “photon rings,” narrow arcs representing pho- global astrometric accuracy in the resulting CRF, and a tons escaping after one or more orbits around the eventhori- number of new cosmological observables. In a uniform zon. The photon ring is described as photons escaping from Hubble expansion there are no transverse motions and thus near the critical impact parameter, b = bc, where (for a non- “all transverse motion is peculiar,” [18]. Any transverse rotating black hole) motions detected, whether kinematic [19], or due to grav- itational lensing [14], or due to gravitational waves [20] bc = 3 √3 RS 5.196 RS, (2) ∼ will, provide information about the large scale distribution of matter in the cosmos. The resulting CRF will undoubt- and R is the Schwarzschild radius. Photons that are tan- S edly also find uses in terrestrial and spacecraft navigation gent at the critical distance are brought closer to the black and in studies of solar system and galactic dynamics. hole and are trapped, orbiting an infinite number of times at the bound orbit radius (b = 3 RS). Just outside the critical distance unstable orbits exist, where photons can orbit the References black hole multiple times and then escape, causing narrow rings in the image plane of the observer. Although these [1] R. P. Linfield, G. S. Levy, J. S. Ulvestad, and et al., rings can be very bright caustics, Gralla et al. [13] conclude “Very Long Baseline Interferometry Observations Us- that the lensed photon ring and lensing ring (representing ing the Tracking and Data Relay Satellite as an Orbit- different numbers of black hole orbits) are not likely to ing Radio Telescope,” in The Impact of VLBI on Astro- dominate the total emission, which will mostly come from physics and Geophysics (M. J. Reid and J. M. Moran, the emission ring. Although the emission ring can be of eds.), vol. 129 of IAU Symposium, p. 457, Jan 1988. variable distance from the center of the SMBH shadow, it will still surround and demark it, providing a suitable refer- [2] H. Hirabayashi, H. Hirosawa, H. Kobayashi, and et ence for accuracy astrometry. al., “The VLBI Space Observatory Programme and the Radio-Astronomical Satellite HALCA,” Pub. As- 4 Astrometric Cosmology tron. Soc. Japan, vol. 52, pp. 955–L965, Dec 2000. [3] N. S. Kardashev, V. V. Khartov, V. V. Abramov, and A multibeam multiantenna VLBI network would enable po- et al., ““RadioAstron”-A telescope with a size of 300 sitioning at the µas level and below and would, if contin- 000 km: Main parameters and first observational re- ued over a suitable duration, also provide information about sults,” Astronomy Reports, vol. 57, pp. 153–194, Mar the angular motions of the SMBH on the sky. While grav- 2013. itational lensing events and microlensing will cause tran- sient source motions [14], most temporal changes at cos- [4] R. W. Porcas, M. J. Rioja, J. Machalski, and mological distances will be secular, appearing as proper H. Hirabayashi, “Phase-Reference Observations with motions. SMBH shadow proper motions from networks VSOP,” in Astrophysical Phenomena Revealed by extending out the Moon and beyond will provide a sound Space VLBI (H. Hirabayashi, P. G. Edwards, and foundation for the development of astrometric cosmology D. W. Murphy, eds.), pp. 245–252, Apr 2000. [5] Y. Y. Kovalev, N. S. Kardashev, K. I. Kellermann, [15] T. M. Eubanks, “Astrometric Cosmology Using Ra- A. P. Lobanov, M. D. Johnson, L. I. Gurvits, P. A. dio Interferometry,” in Bulletin of the American As- Voitsik, J. A. Zensus, J. M. Anderson, U. Bach, D. L. tronomical Society, vol. 23, p. 1465, Sep 1991. Jauncey, F. Ghigo, T. Ghosh, A. Kraus, Y. A. Ko- valev, M. M. Lisakov, L. Y. Petrov, J. D. Romney, [16] D. MacMillan, A. Fey, J. Gipson, D. Gordon, C. Ja- C. J. Salter, and K. V. Sokolovsky, “RadioAstron Ob- cobs, H. Krásná, S. Lambert, C. Ma, Z. Malkin, servations of the Quasar 3C273: A Challenge to the O. Titov, G. Wang, M. Xu, and N. Zacharias, “Galac- Brightness Temperature Limit,” Ap. J. Lett., vol. 820, tic Aberration in VLBI Analysis:,” in International p. L9, Mar 2016. VLBI Service for Geodesy and Astrometry 2018 Gen- eral Meeting Proceedings: ”Global Geodesy and the [6] Event Horizon Telescope Collaboration, “First M87 Role of VGOS - Fundamentalto Sustainable Develop- Event Horizon Telescope Results. I. The Shadow of ment, pp. 163–168, Sep 2019. the Supermassive Black Hole,” Ap. J. Lett., vol. 875, [17] U. Sperhake, R. Rosca-Mead, D. Gerosa, and E. Berti, p. L1, Apr 2019. “Amplification of superkicks in black-hole binaries [7] V. L. Fish, M. Shea, and K. Akiyama, “Imaging black through orbital eccentricity,” Phys. Rev. D, vol. 101, holes and jets with a VLBI array including multiple p. 024044, Jan 2020. space-based telescopes,” Advances in Space Research, [18] J. Darling and A. E. Truebenbach, “All Transverse vol. 65, pp. 821–830, Jan 2020. Motion Is Peculiar: Connecting the Proper Motions of Galaxies to the Matter Power Spectrum,” Ap. J., [8] A. L. Fey, D. Gordon, C. S. Jacobs, C. Ma, R. A. vol. 864, p. 37, Sep 2018. Gaume, E. F. Arias, G. Bianco, D. A. Boboltz, S. Böckmann, S. Bolotin, P. Charlot, A. Collioud, [19] J. Paine, J. Darling, R. Graziani, and H. M. Cour- G. Engelhardt, J. Gipson, A.-M. Gontier, R. Heinkel- tois, “Secular Extragalactic Parallax: Measurement mann, S. Kurdubov, S. Lambert, S. Lytvyn, D. S. Methods and Predictions for Gaia,” arXiv e-prints, MacMillan, Z. Malkin, A. Nothnagel, R. Ojha, p. arXiv:1912.11935, Dec 2019. E. Skurikhina, J. Sokolova, J. Souchay, O. J. Sovers, V. Tesmer, O. Titov, G. Wang, and V. Zharov, “The [20] C. R. Gwinn, T. M. Eubanks, T. Pyne, M. Birkin- Second Realization of the International Celestial Ref- shaw, and D. N. Matsakis, “Quasar Proper Motions erence Frame by Very Long Baseline Interferometry,” and Low-Frequency Gravitational Waves,” Ap. J., Astron. J., vol. 150, p. 58, Aug. 2015. vol. 485, pp. 87–91, Aug 1997.

[9] J. M. Marcaide, I. I. Shapiro, J. Ball, N. Bartel, T. Clark, B. E. Corey, M. V. Gorenstein, R. A. Pre- ston, M. I. Ratner, A. E. E. Rogers, and A. R. Whit- ney, “Nearly Microarcsecond Precision Differential Astrometry,” in VLBI and Compact Radio Sources (R. Fanti, K. I. Kellermann, and G. Setti, eds.), vol. 110 of IAU Symposium, p. 361, Jan 1984.

[10] H. Kobayashi, “VERA and East Asian VLBI net- work,” in Resolving The Sky - Radio Interferometry: Past, Present and Future, p. 29, Jan 2012.

[11] O. Y. Tsupko, Z. Fan, and G. S. Bisnovatyi-Kogan, “Black hole shadow as a standard ruler in cosmology,” arXiv e-prints, p. arXiv:1905.10509, May 2019.

[12] S. Vagnozzi, C. Bambi, and L. Visinelli, “Concerns regarding the use of black hole shadows as stan- dard rulers,” arXiv e-prints, p. arXiv:2001.02986, Jan 2020.

[13] S. E. Gralla, D. E. Holz, and R. M. Wald, “Black hole shadows, photon rings, and lensing rings,” Phys. Rev. D, vol. 100, p. 024018, Jul 2019.

[14] T. Yano, “Cosmic Error Caused by the Gravitational Microlensing Effect in High-precision Astrometry,” Ap. J., vol. 757, p. 189, Oct 2012.