Suchergebnisse

Context. While studies of large samples of jets of active galactic nuclei (AGN) are important in order to establish a global picture, dedicated single-source studies are an invaluable tool for probing crucial processes within jets on parsec scales. These processes involve in particular the formation and geometry of the jet magnetic field as well as the flow itself. Aims. We aim to better understand the dynamics within relativistic magneto-hydrodynamical flows in the extreme environment and close vicinity of supermassive black holes. Methods. We analyze the peculiar radio galaxy 3C 111, for which long-Term polarimetric observations are available. We make use of the high spatial resolution of the VLBA network and the MOJAVE monitoring program, which provides high data quality also for single sources and allows us to study jet dynamics on parsec scales in full polarization with an evenly sampled time-domain. While electric vectors can probe the underlying magnetic field, other properties of the jet such as the variable (polarized) flux density, feature size, and brightness temperature, can give valuable insights into the flow itself. We complement the VLBA data with data from the IRAM 30-m Telescope as well as the SMA. Results. We observe a complex evolution of the polarized jet. The electric vector position angles (EVPAs) of features traveling down the jet perform a large rotation of &180across a distance of about 20 pc. As opposed to this smooth swing, the EVPAs are strongly variable within the first parsecs of the jet.We find an overall tendency towards transverse EVPAs across the jet with a local anomaly of aligned vectors in between. The polarized flux density increases rapidly at that distance and eventually saturates towards the outermost observable regions. The transverse extent of the flow suddenly decreases simultaneously to a jump in brightness temperature around where we observe the EVPAs to turn into alignment with the jet flow. Also the gradient of the feature size and particle density with distance steepens significantly at that region. Conclusions. We interpret the propagating polarized features as shocks and the observed local anomalies as the interaction of these shocks with a localized recollimation shock of the underlying flow. Together with a sheared magnetic field, this shock-shock interaction can explain the large rotation of the EVPA. The superimposed variability of the EVPAs close to the core is likely related to a clumpy Faraday screen, which also contributes significantly to the observed EVPA rotation in that region.© ESO 2018. ; T.B. is grateful for fruitful discussions with M. Lyutikov and E. Kravchenko, and thanks D. Gabuzda for her input that helped to improve the paper. We thank the National Radio Astronomy Observatory, a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.; this research has made use of data from the MOJAVE database that is maintained by the MOJAVE team (Lister et al. 2009) and IRAM, which is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). This work has benefited from research funding from the European Community's Seventh Framework Programme. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. We also thank J. E. Davis for the development of the slxfig module that has been used to prepare the figures in this work. We made use of ISIS functions provided by ECAP/Remeis observatory and MIT (http://www.sternwarte.uni-erlangen.de/isis/) as well as the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. M.P. acknowledges financial support from the Spanish Ministerio de Economia, Industria y Competitividad (grants AYA2013-40979-P, and AYA2013-48226- C3-2-P) and from the local Valencian Government (Generalitat Valenciana, grant Prometeo-11/2014/069). I.A. acknowledges support by a Ramon y Cajal grant of the Ministerio de Economia y Competitividad (MINECO) of Spain. A.B.P. and Y.Y.K. acknowledge partial support by the Russian Foundation for Basic Research (grant 17-02-00197), by the government of the Russian Federation (agreement 05.Y09.21.0018), and by the Alexander von Humboldt Foundation. E.R. acknowledges partial support by the Spanish MINECO project AYA2012-38491-C02-01 and by the Generalitat Valenciana project PROMETEOII/2014/057. T.S. was funded by the Academy of Finland projects 274477 and 284495. The research at the IAA-CSIC was supported in part by the MINECO through grants AYA2016-80889-P, AYA2013-40825-P, and AYA2010-14844, and by the regional government of Andalucia through grant P09-FQM-4784.

We report the results of decade-long (2008-2018) γ -ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, Fermi and Swift data, as well as polarimetric and spectroscopic data. The X-ray and γ -ray light curves correlate well, with no delay ≳ 3 h, implying general cospatiality of the emission regions. The γ -ray-optical flux-flux relation changes with activity state, ranging from a linear to amore complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz Very Long Baseline Array images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain γ -ray variability on very short time-scales. The MgII emission line flux in the 'blue' and 'red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands, we find progressive delays of the most prominent light-curve maxima with decreasing frequency, as expected from the frequency dependence of the τ= 1 surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at ∼5 GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet. © 2020 The Author(s). ; We thank the referee for attentive reading and comments that helped to improve presentation of the manuscript. The data collected by the WEBT collaboration are stored in the WEBT archive at the Osservatorio Astrofisico di Torino -INAF (ht tp://www.oato.inaf.it/blazars/webt/); for questions regarding their availability, please contact the WEBT President Massimo Villata(massimo.villata@inaf.it).TheSt.Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. The research at BU was supported in part by National Science Foundation grant AST-1615796 and NASA Fermi Guest Investigator grants 80NSSC17K0649, 80NSSC19K1504, and 80NSSC19K1505. The PRISM camera at Lowell Observatory was developed by K. Janes et al. at BU and Lowell Observatory, with funding from the NSF, BU, and Lowell Observatory. The emission-line observations made use of the DCT at Lowell Observatory, supported by Discovery Communications, Inc., BU, the University of Maryland, the University of Toledo, and Northern Arizona University. The VLBA is an instrument of the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the US NSF, operated under cooperative agreement by Associated Universities, Inc. This research has used data from the UMRAO which was supported by the University of Michigan; research at this facility was supported by NASA under awards NNX09AU16G, NNX10AP16G, NNX11AO13G, and NNX13AP18G, and by the NSF under award AST-0607523. The Steward Observatory spectropolarimetric monitoring project was supported by NASA Fermi Guest Investigator grants NNX08AW56G, NNX09AU10G, NNX12AO93G, and NNX15AU81G. The Torino group acknowledges financial contribution from agreement ASI-INAF n.2017-14-H.0 and from contract PRIN-SKA-CTA-INAF 2016. I.A. acknowledges support by a Ramon y Cajal grant (RYC-2013-14511) of the 'Ministerio de Ciencia, Innovacion, y Universidades (MICIU)' of Spain and from MCIU through the 'Center of Excellence Severo Ochoa' award for the Instituto de Astrofisica de Andalucia-CSIC (SEV-20170709). Acquisition and reduction of the POLAMI and MAPCAT data were supported by MICIU through grant AYA2016-80889-P. The POLAMI observations were carried out at the IRAM 30-m Telescope, supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). The MAPCAT observations were carried out at theGerman-Spanish Calar Alto Observatory, jointly operated by the Max-Plank-Institut fur Astronomie and the Instituto de Astrofisica de Andalucia-CSIC. The study is based partly on data obtained with the STELLA robotic telescopes in Tenerife, an AIP facility jointly operated by AIP and IAC. The OVRO 40-m monitoring program is supported in part by NASA grants NNX08AW31G, NNX11A043G, and NNX14AQ89G, and NSF grants AST-0808050 and AST-1109911. TH was supported by the Academy of Finland projects 317383 and 320085. AZT-24 observations were made within an agreement between Pulkovo, Rome and Teramo observatories. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. The Abastumani team acknowledges financial support by the Shota Rustaveli National Science Foundation under contract FR/217950/16. r This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants DN 081/2016, DN 18-13/2017, KP-06-H28/3 (2018), and KP-06-PN38/1 (2019), Bulgarian National Science Programme 'Young Scientists and Postdoctoral Students 2019', Bulgarian National Science Fund under grant DN18-10/2017 and National RI Roadmap Projects DO1-157/28.08.2018 and DO1-153/28.08.2018 of the Ministry of Education and Science of the Republic of Bulgaria. GD and OV gratefully acknowledge observing grant support from the Institute of Astronomy and Rozhen National Astronomical Observatory via bilateral joint research project `Study of ICRF radio-sources and fast variable astronomical objects' (head -G. Damljanovic). This work was partly supported by the National Science Fund of the Ministry of Education and Science of Bulgaria under grant DN 08-20/2016, and by project RD-08-37/2019 of the University of Shumen. This work is a part of projects nos 176011, 176004, and 176021, supported by theMinistry of Education, Science and Technological Development of the Republic of Serbia. MGM acknowledges support through the Russian Government Program of Competitive Growth of Kazan Federal University. The Astronomical Observatory of the Autonomous Region of the Aosta Valley (OAVdA) is managed by the Fondazione Clement Fillietroz-ONLUS, which is supported by the Regional Government of the Aosta Valley, the Town Municipality of Nus and the 'Unite des Communes vald 'otainesMont-Emilius'. The research at the OAVdA was partially funded by several `Research and Education' annual grants from Fondazione CRT. This article is partly based on observations made with the IAC80 and TCS telescopes operated by the Instituto de Astrofisica de Canarias in the Spanish Observatorio del Teide on the island of Tenerife. A part of the observations were carried out using theRATAN-600 scientific equipment (SAO of the Russian Academy of Sciences).

The very high energy (VHE ¿ 100 GeV) -ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z = 0.368 or z ¿ 0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. Methods. We investigated the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and have put it in context with possible emission scenarios. We collected photometric data in radio, optical, X-ray, and ¿-ray. We studied both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Results. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy ¿-ray(HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100. Conclusions. The high emission state during the flare allows us to compile the simultaneous broadband SED and to characterize it in the scope of blazar jet emission models. The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars (FSRQs), that is, inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to FSRQs and low-frequency peaked BL Lac objects (LBL). © ESO 2018. ; The financial support of the German BMBF and MPG, the Italian INFN and INAF, the Swiss National Fund SNF, the ERDF under the Spanish MINECO (FPA2015-69818-P, FPA2012-36668, FPA2015-68378-P, FPA2015-69210-C6-2-R, FPA2015-69210-C6-4-R, FPA2015-69210-C6-6-R, AYA2015-71042-P, AYA2016-76012-C3-1-P, ESP2015-71662-C2-2-P, CSD2009-00064), and the Japanese JSPS and MEXT is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia "Severo Ochoa" SEV-2012-0234 and SEV-2015-0548, and Unidad de Excelencia "María de Maeztu" MDM-2014-0369, by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, the Polish National Research Centre grant UMO-2016/22/M/ST9/00382 and by the Brazilian MCTIC, CNPq, and FAPERJ. IA acknowledges support by a Ramón y Cajal grant of the Ministerio de Economía, Industria y Competitividad (MINECO) of Spain. The research at the IAA–CSIC was supported in part by the MINECO through grants AYA2016–80889–P, AYA2013–40825–P, and AYA2010–14844, and by the regional government of Andalucía through grant P09–FQM–4784. ; Peer Reviewed

Context. Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are the most important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a direct comparison with observational data. Aims. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. Methods. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a millimetre VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects, compared to the addition of thermal noise only. Using synthetic data based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M 87, we performed case studies to assess the image quality that can be obtained with the current and future EHT array for different weather conditions. Results. Our synthetic observations show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of our GRMHD source models can be recovered robustly with the EHT2017 array after performing calibration steps, which include fringe fitting, a priori amplitude and network calibration, and self-calibration. With the planned addition of new stations to the EHT array in the coming years, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. © 2020 ESO. ; This work is supported by the ERC Synergy Grant "BlackHoleCam: Imaging the Event Horizon of Black Holes" (Grant 610058). I. Natarajan and R. Deane are grateful for the support from the New Scientific Frontiers with Precision Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory (SARAO), which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Technology (DST) of South Africa. The authors of the present paper further thank the following organizations and programmes: the Academy of Finland (projects 274477, 284495, 312496); the Advanced European Network of E-infrastructures for Astronomy with the SKA (AENEAS) project, supported by the European Commission Framework Programme Horizon 2020 Research and Innovation action under grant agreement 731016; the Alexander von Humboldt Stiftung; the Black Hole Initiative at Harvard University, through a grant (60477) from the John Templeton Foundation; the China Scholarship Council; Comision Nacional de Investigacio Cientifica y Tecnologica (CONICYT, Chile, via PIA ACT172033, Fondecyt 1171506, BASAL AFB-170002, ALMAconicyt 31140007); Consejo Nacional de Ciencia y Tecnologia (CONACYT, Mexico, projects 104497, 275201, 279006, 281692); the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Direccion General de Asuntos del Personal Academico-Universidad Nacional Autonoma de Mexico (DGAPA-UNAM, project IN112417); the Generalitat Valenciana postdoctoral grant APOSTD/2018/177; the Gordon and Betty Moore Foundation (grants GBMF-3561, GBMF-5278); the Istituto Nazionale di Fisica Nucleare (INFN) sezione di Napoli, iniziative specifiche TEONGRAV; the GenT Program (Generalitat Valenciana) under project CIDEGENT/2018/021; the International Max Planck Research School for Astronomy and Astrophysics at the Universities of Bonn and Cologne; the Jansky Fellowship program of the National Radio Astronomy Observatory (NRAO); the Japanese Government (Monbukagakusho: MEXT) Scholarship; the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for JSPS Research Fellowship (JP17J08829); the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (CAS, grants QYZDJ-SSW-SLH057, QYZDJ-SSW-SYS008); the Leverhulme Trust Early Career Research Fellowship; the Max-Planck-Gesellschaft (MPG); the Max Planck Partner Group of the MPG and the CAS; the MEXT/JSPS KAKENHI (grants 18KK0090, JP18K13594, JP18K03656, JP18H03721, 18K03709, 18H01245, 25120007); the MIT International Science and Technology Initiatives (MISTI) Funds; the Ministry of Science and Technology (MOST) of Taiwan (105-2112-M-001-025-MY3, 106-2112-M-001-011, 106-2119-M-001027, 107-2119-M-001-017, 107-2119-M-001-020, and 107-2119-M-110-005); the National Aeronautics and Space Administration (NASA, Fermi Guest Investigator grant 80NSSC17K0649); NASA through the NASA Hubble Fellowship grant #HST-HF2-51431.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. , for NASA, under contract NAS5-26555; the National Institute of Natural Sciences (NINS) of Japan; the National Key Research and Development Program of China (grant 2016YFA0400704, 2016YFA0400702); the National Science Foundation (NSF, grants AST-0096454, AST-0352953, AST-0521233, AST-0705062, AST-0905844, AST-0922984, AST-1126433, AST-1140030, DGE-1144085, AST-1207704, AST-1207730, AST-1207752, MRI-1228509, OPP-1248097, AST-1310896, AST-1312651, AST-1337663, AST-1440254, AST-1555365, AST-1715061, AST-1615796, AST-1716327, OISE-1743747, AST-1816420); the Natural Science Foundation of China (grants 11573051, 11633006, 11650110427, 10625314, 11721303, 11725312, 11933007); the Natural Sciences and Engineering Research Council of Canada (NSERC, including a Discovery Grant and the NSERC Alexander Graham Bell Canada Graduate Scholarships-Doctoral Program); the National Youth Thousand Talents Program of China; the National Research Foundation of Korea (the Global PhD Fellowship Grant: grants NRF-2015H1A2A1033752, 2015-R1D1A1A01056807, the Korea Research Fellowship Program: NRF-2015H1D3A1066561); the Netherlands Organization for Scientific Research (NWO) VICI award (grant 639.043.513) and Spinoza Prize SPI 78-409; the New Scientific Frontiers with Precision Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory (SARAO), which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Technology (DST) of South Africa; the Onsala Space Observatory (OSO) national infrastructure, for the provisioning of its facilities/observational support (OSO receives funding through the Swedish Research Council under grant 2017-00648) the Perimeter Institute for Theoretical Physics (research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development and by the Province of Ontario through the Ministry of Research, Innovation and Science); the Princeton/Flatiron Postdoctoral Prize Fellowship; the Russian Science Foundation (grant 17-12-01029); the Spanish Ministerio de Economia y Competitividad (grants AYA2015-63939-C21-P, AYA2016-80889-P); the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709); the Toray Science Foundation; the US Department of Energy (USDOE) through the Los Alamos National Laboratory (operated by Triad National Security, LLC, for the National Nuclear Security Administration of the USDOE (Contract 89233218CNA000001)); the Italian Ministero dell'Istruzione Universita e Ricerca through the grant Progetti Premiali 2012-iALMA (CUP C52I13000140001); the European Union's Horizon 2020 research and innovation programme under grant agreement No 730562 RadioNet; ALMA North America Development Fund; the Academia Sinica; Chandra TM6-17006X. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), supported by NSF grant ACI-1548562, and CyVerse, supported by NSF grants DBI-0735191, DBI-1265383, and DBI1743442. XSEDE Stampede2 resource at TACC was allocated through TGAST170024 and TG-AST080026N. XSEDE JetStream resource at PTI and TACC was allocated through AST170028. The simulations were performed in part on the SuperMUC cluster at the LRZ in Garching, on the LOEWE cluster in CSC in Frankfurt, and on the HazelHen cluster at the HLRS in Stuttgart. This research was enabled in part by support provided by Compute Ontario (http://computeontario.ca), Calcul Quebec (http://www. calculquebec.ca) and Compute Canada (http://www.computecanada.ca).We thank the sta ff at the participating observatories, correlation centers, and institutions for their enthusiastic support. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2017.1.00841.V. ALMA is a partnership of the European Southern Observatory (ESO; Europe, representing its member states), NSF, and National Institutes of Natural Sciences of Japan, together with National Research Council (Canada), Ministry of Science and Technology (MOST; Taiwan), Academia Sinica Institute of Astronomy and Astrophysics (ASIAA; Taiwan), and Korea Astronomy and Space Science Institute (KASI; Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, Associated Universities, Inc. (AUI)/NRAO, and the National Astronomical Observatory of Japan (NAOJ). The NRAO is a facility of the NSF operated under cooperative agreement by AUI. APEX is a collaboration between the Max-Planck-Institut fur Radioastronomie (Germany), ESO, and the Onsala Space Observatory (Sweden). The SMA is a joint project between the SAO and ASIAA and is funded by the Smithsonian Institution and the Academia Sinica. The JCMT is operated by the East Asian Observatory on behalf of the NAOJ, ASIAA, and KASI, as well as the Ministry of Finance of China, Chinese Academy of Sciences, and the National Key R&D Program (No. 2017YFA0402700) of China. Additional funding support for the JCMT is provided by the Science and Technologies Facility Council (UK) and participating universities in the UK and Canada. The LMT is a project operated by the Instituto Nacional de Astrofisica, Optica, y Electronica (Mexico) and the University of Massachusetts at Amherst (USA). The IRAM 30m telescope on Pico Veleta, Spain is operated by IRAM and supported by CNRS (Centre National de la Recherche Scientifique, France), MPG (Max-Planck-Gesellschaft, Germany) and IGN (Instituto Geografico Nacional, Spain). The SMT is operated by the Arizona Radio Observatory, a part of the Steward Observatory of the University of Arizona, with financial support of operations from the State of Arizona and financial support for instrumentation development from the NSF. The SPT is supported by the National Science Foundation through grant PLR-1248097. Partial support is also provided by the NSF Physics Frontier Center grant PHY-1125897 to the Kavli Institute of Cosmological Physics at the University of Chicago, the Kavli Foundation and the Gordon and Betty Moore Foundation grant GBMF 947. The SPT hydrogen maser was provided on loan from the GLT, courtesy of ASIAA. The EHTC has received generous donations of FPGA chips from Xilinx Inc., under the Xilinx University Program. The EHTC has benefited from technology shared under open-source license by the Collaboration for Astronomy Signal Processing and Electronics Research (CASPER). The EHT project is grateful to T4Science and Microsemi for their assistance with Hydrogen Masers. This research has made use of NASA's Astrophysics Data System. We gratefully acknowledge the support provided by the extended staff of the ALMA, both from the inception of the ALMA Phasing Project through the observational campaigns of 2017 and 2018. We would like to thank A. Deller and W. Brisken for EHT-specific support with the use of DiFX. We acknowledge the significance that Maunakea, where the SMA and JCMT EHT stations are located, has for the indigenous Hawaiian people. The software presented in this work makes use of the Numpy (van derWalt et al. 2011), Scipy (Jones et al. 2001), Astropy (Astropy Collaboration 2013, 2018) libraries and the KERN software bundle (Molenaar & Smirnov 2018).

The Event Horizon Telescope Collaboration ; Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope, and high- precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD, and KORAL. Agreement among the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. © 2019. The American Astronomical Society. All rights reserved. ; R.N. thanks the National Science Foundation (NSF; grants OISE-1743747, AST-1816420) and acknowledges computational support from the NSF via XSEDE resources (grant TG-AST080026N). L.D.Z. acknowledges support from the PRIN-MIUR project Multi-scale Simulations of High-Energy Astrophysical Plasmas (Prot. 2015L5EE2Y) and from the INFN-TEONGRAV initiative. C.J.W. made use of thComision Nacional de Investigacion Cientifica y Tecnologica (CONICYT, ChileComision Nacional de Investigacion Cientifica y Tecnologica (CONICYT, Chilee Extreme Science and Engineering Discovery Environment (XSEDE) Comet at the San Diego Supercomputer Center through allocation AST170012. The H-AMR high-resolution simulation was made possible by NSF PRAC award Nos. 1615281 and OAC-1811605 at the Blue Waters sustained-petascale computing project and supported in part under grant No. NSF PHY-1125915 (PI A. Tchekhovskoy). K.C. and S.M. are supported by the Netherlands Organization for Scientific Research (NWO) VICI grant (No. 639.043.513); M.L. is supported by the NWO Spinoza Prize (PI M.B.M. van der Klis). The HARM-Noble simulations were made possible by NSF PRAC award No. NSF OAC-1515969, OAC-1811228 at the Blue Waters sustained-petascale computing project, and supported in part under grant No. NSF PHY-1125915. The BHAC CKS-GRMHD simulations were performed on the Dutch National Supercomputing cluster Cartesius and are funded by the NWO computing grant 16431. S.C.N. was supported by an appointment to the NASA Postdoctoral Program at the Goddard Space Flight Center administered by USRA through a contract with NASA. Y.M., H.O., O.P., and L.R. acknowledge support from the ERC synergy grant >BlackHoleCam: Imaging the Event Horizon of Black Holes> (grant No. 610058). M.B. acknowledges support from the European Research Council (grant No. 715368-MagBURST) and from the Gauss Centre for Supercomputing e.V. (www.Gauss-centre.eu) for funding this project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (www.lrz.de).P.C.F.was supported by NSF grant AST-1616185 and used resources from the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by NSF grant No. ACI-1548562. Work by P.A. was performed in part under the auspices of the US Department of Energy by Lawrence Livermore National Laboratory under contract DE-AC52-07NA27344. The authors of the present paper further thank the following organizations and programs: the Academy of Finland (projects 274477, 284495, 312496); the Advanced European Network of E-infrastructures for Astronomy with the SKA (AENEAS) project, supported by the European Commission Framework Programme Horizon 2020 Research and Innovation action under grant agreement 731016; the Alexander von Humboldt Stiftung; the Black Hole Initiative at Harvard University, through a grant (60477) from the John Templeton Foundation; the China Scholarship Council; Comision Nacional de Investigacion Cientifica y Tecnologica (CONICYT, Chile, via PIA ACT172033, Fondecyt 1171506, BASAL AFB-170002, ALMA-conicyt 31140007); Consejo Nacional de Ciencia y Tecnologia (CONACYT, Mexico, projects 104497, 275201, 279006, 281692); the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Direccion General de Asuntos del Personal Academico-Universidad Nacional Autonoma de Mexico (DGAPA-UNAM, project IN112417); the European Research Council Synergy Grant >BlackHoleCam: Imaging the Event Horizon of Black Holes> (grant 610058); the Generalitat Valenciana postdoctoral grant APOSTD/2018/177; the Gordon and Betty Moore Foundation (grants GBMF-3561, GBMF-5278); the Istituto Nazionale di Fisica Nucleare (INFN) sezione di Napoli, iniziative specifiche TEONGRAV; the International Max Planck Research School for Astronomy and Astrophysics at the Universities of Bonn and Cologne; the Jansky Fellowship program of the National Radio Astronomy Observatory (NRAO); the Japanese Government (Monbukagakusho: MEXT) Scholarship; the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for JSPS Research Fellowship (JP17J08829); the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (CAS, grants QYZDJ-SSW-SLH057, QYZDJ-SSW-SYS008); the Leverhulme Trust Early Career Research Fellowship; the Max-Planck-Gesellschaft (MPG); the Max Planck Partner Group of the MPG and the CAS; the MEXT/JSPS KAKENHI (grants 18KK0090, JP18K13594, JP18K03656, JP18H03721, 18K03709, 18H01245, 25120007); the MIT International Science and Technology Initiatives (MISTI) Funds; the Ministry of Science and Technology (MOST) of Taiwan (105-2112-M-001-025-MY3, 106-2112-M-001-011, 106-2119-M-001-027, 107-2119-M-001-017, 107-2119-M-001-020, and 107-2119-M-110-005); the National Aeronautics and Space Administration (NASA, Fermi Guest Investigator grant 80NSSC17K0649); the National Institute of Natural Sciences (NINS) of Japan; the National Key Research and Development Program of China (grant 2016YFA0400704, 2016YFA0400702); the National Science Foundation (NSF, grants AST-0096454, AST-0352953, AST-0521233, AST-0705062, AST-0905844, AST-0922984, AST-1126433, AST-1140030, DGE-1144085, AST-1207704, AST-1207730, AST-1207752, MRI-1228509, OPP-1248097, AST-1310896, AST-1312651, AST-1337663, AST-1440254, AST-1555365, AST-1715061, AST-1615796, AST-1716327, OISE-1743747, AST-1816420); the Natural Science Foundation of China (grants 11573051, 11633006, 11650110427, 10625314, 11721303, 11725312); the Natural Sciences and Engineering Research Council of Canada (NSERC, including a Discovery Grant and the NSERC Alexander Graham Bell Canada Graduate Scholarships Doctoral Program); the National Youth Thousand Talents Program of China; the National Research Foundation of Korea (the Global PhD Fellowship Grant: grants NRF-2015H1A2A1033752, 2015-R1D1A1A01056807, the Korea Research Fellowship Program: NRF-2015H1D3A1066561); the Netherlands Organization for Scientific Research (NWO) VICI award (grant 639.043. 513) and Spinoza Prize SPI 78-409; the New Scientific Frontiers with Precision Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory (SARAO), which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Technology (DST) of South Africa; the Onsala Space Observatory (OSO) national infrastructure, for the provisioning of its facilities/observational support (OSO receives funding through the Swedish Research Council under grant 2017-00648); the Perimeter Institute for Theoretical Physics (research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development, and by the Province of Ontario through the Ministry of Research, Innovation and Science); the Russian Science Foundation (grant 17-12-01029); the Spanish Ministerio de Economia y Competitividad (grants AYA2015-63939-C2-1-P, AYA2016-80889-P); the State Agency for Research of the Spanish MCIU through the >Center of Excellence Severo Ochoa> award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709); the Toray Science Foundation; the US Department of Energy (USDOE) through the Los Alamos National Laboratory (operated by Triad National Security, LLC, for the National Nuclear Security Administration of the USDOE (Contract 89233218CNA000001)); the Italian Ministero dell'Istruzione Universita e Ricerca through the grant Progetti Premiali 2012-iALMA (CUP C52I13000140001); the European Union's Horizon 2020 research and innovation programme under grant agreement No. 730562 RadioNet; ALMA North America Development Fund; the Academia Sinica; and Chandra TM6-17006X. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), supported by NSF grant ACI-1548562, and CyVerse, supported by NSF grants DBI-0735191, DBI-1265383, and DBI-1743442. XSEDE Stampede2 resource at TACC was allocated through TG-AST170024 and TG-AST080026N. XSEDE JetStream resource at PTI and TACC was allocated through AST170028. The simulations were performed in part on the SuperMUC cluster at the LRZ in Garching, on the LOEWE cluster in CSC in Frankfurt, and on the HazelHen cluster at the HLRS in Stuttgart. This research was enabled in part by support provided by Compute Ontario (http://computeontario.ca), Calcul Quebec (http://www.calculquebec.ca), and Compute Canada (http://www.computecanada.ca).