THEMIS: A Parameter Estimation Framework for the Event Horizon Telescope
This is an open accBroderick, Avery E.; Gold, Roman; Karami, Mansour; Preciado-López, Jorge A.; Tiede, Paul; Pu, Hung-Yi; Akiyama, Kazunori; Alberdi, Antxon; Alef, Walter; Asada, Keiichi; Azulay, Rebecca; Baczko, Anne-Kathrin; Baloković, Mislav; Barrett, John; Bintley, Dan; Blackburn, Lindy; Boland, Wilfred; Bouman, Katherine L.; Bower, Geoffrey C.; Bremer, Michael; Brinkerink, Christiaan D.; Brissenden, Roger; Britzen, Silke; Broguiere, Dominique; Bronzwaer, Thomas; Byun, Do-Young; Carlstrom, John E.; Chael, Andrew; Chatterjee, Shami; Chatterjee, Koushik; Chen, Ming-Tang; Chen, Yongjun; Cho, Ilje; Conway, John E.; Cordes, James M.; Crew, Geoffrey B.; Cui, Yuzhu; Davelaar, Jordy; De Laurentis, Mariafelicia; Deane, Roger; Dempsey, Jessica; Desvignes, Gregory; Doeleman, Sheperd S.; Eatough, Ralph P.; Falcke, Heino; Fish, Vincent L.; Fomalont, Ed; Fraga-Encinas, Raquel; Friberg, Per; Fromm, Christian M.; Galison, Peter; Gammie, Charles F.; García, Roberto; Gentaz, Olivier; Georgiev, Boris; Goddi, Ciriaco; Gómez, José L.; Gu, Minfeng; Gurwell, Mark; Hada, Kazuhiro; Hecht, Michael H.; Hesper, Ronald; Ho, Luis C.; Ho, Paul; Honma, Mareki; Huang, Chih-Wei L.; Huang, Lei; Hughes, David H.; Inoue, Makoto; Issaoun, Sara; James, David J.; Janssen, Michael; Jeter, Britton; Jiang, Wu; Jiménez-Rosales, Alejandra; Johnson, Michael D.; Jorstad, Svetlana; Jung, Taehyun; Karuppusamy, Ramesh; Kawashima, Tomohisa; Keating, Garrett K.; Kettenis, Mark; Kim, Jae-Young; Kim, Jongsoo; Kino, Motoki; Koay, Jun Yi; Koch, Patrick M.; Koyama, Shoko; Kramer, Michael; Kramer, Carsten; Krichbaum, Thomas P.; Kuo, Cheng-Yu; Lee, Sang-Sung; Li, Yan-Rong; Li, Zhiyuan; Lindqvist, Michael; Lico, Rocco; Liu, Kuo; Liuzzo, Elisabetta; Lo, Wen-Ping; Lobanov, Andrei P.; Loinard, Laurent; Lonsdale, Colin; Lu, Ru-Sen; MacDonald, Nicholas R.; Mao, Jirong; Marscher, Alan P.; Martí-Vidal, Iván; Matsushita, Satoki; Matthews, Lynn D.; Menten, Karl M.; Mizuno, Yosuke; Mizuno, Izumi; Moran, James M.; Moriyama, Kotaro; Moscibrodzka, Monika; Müller, Cornelia; Nagai, Hiroshi; Nagar, Neil M.; Nakamura, Masanori; Narayan, Ramesh; Narayanan, Gopal; Natarajan, Iniyan; Neri, Roberto; Ni, Chunchong; Noutsos, Aristeidis; Okino, Hiroki; Olivares, Héctor; Ortiz-León, Gisela N.; Oyama, Tomoaki; Palumbo, Daniel C. M.; Park, Jongho; Pen, Ue-Li; Pesce, Dominic W.; Piétu, Vincent; Plambeck, Richard; PopStefanija, Aleksandar; Porth, Oliver; Prather, Ben; Ramakrishnan, Venkatessh; Rao, Ramprasad; Rawlings, Mark G.; Raymond, Alexander W.; Rezzolla, Luciano; Ripperda, Bart; Roelofs, Freek; Rogers, Alan; Ros, Eduardo; Rose, Mel; Rottmann, Helge; Ruszczyk, Chet; Ryan, Benjamin R.; Rygl, Kazi L. J.; Sánchez, Salvador; Sánchez-Arguelles, David; Sasada, Mahito; Savolainen, Tuomas; Schloerb, F. Peter; Schuster, Karl-Friedrich; Shao, Lijing; Shen, Zhiqiang; Small, Des; Sohn, Bong Won; SooHoo, Jason; Tazaki, Fumie; Tilanus, Remo P. J.; Titus, Michael; Toma, Kenji; Torne, Pablo; Traianou, Efthalia; Trippe, Sascha; Tsuda, Shuichiro; van Bemmel, Ilse; van Langevelde, Huib Jan; van Rossum, Daniel R.; Wagner, Jan; Wardle, John; Weintroub, Jonathan; Wex, Norbert; Wharton, Robert; Wielgus, Maciek; Wong, George N.; Wu, Qingwen; Yoon, Doosoo; Young, André; Young, Ken; Younsi, Ziri; Yuan, Feng; Yuan, Ye-Fei; Zensus, J. Anton; Zhao, Guangyao; Zhao, Shan-Shan; Zhu, Ziyan; Event Horizon Telescope Collaboration ess article.-- Full list of authors: ; The Event Horizon Telescope (EHT) provides the unprecedented ability to directly resolve the structure and dynamics of black hole emission regions on scales smaller than their horizons. This has the potential to critically probe the mechanisms by which black holes accrete and launch outflows, and the structure of supermassive black hole spacetimes. However, accessing this information is a formidable analysis challenge for two reasons. First, the EHT natively produces a variety of data types that encode information about the image structure in nontrivial ways; these are subject to a variety of systematic effects associated with very long baseline interferometry and are supplemented by a wide variety of auxiliary data on the primary EHT targets from decades of other observations. Second, models of the emission regions and their interaction with the black hole are complex, highly uncertain, and computationally expensive to construct. As a result, the scientific utilization of EHT observations requires a flexible, extensible, and powerful analysis framework. We present such a framework, Themis, which defines a set of interfaces between models, data, and sampling algorithms that facilitates future development. We describe the design and currently existing components of Themis, how Themis has been validated thus far, and present additional analyses made possible by Themis that illustrate its capabilities. Importantly, we demonstrate that Themis is able to reproduce prior EHT analyses, extend these, and do so in a computationally efficient manner that can efficiently exploit modern high-performance computing facilities. Themis has already been used extensively in the scientific analysis and interpretation of the first EHT observations of M87. © 2020. The Author(s). Published by the American Astronomical Society. ; This work was made possible by the facilities of the Shared Hierarchical Academic Research Computing Network (SHARCNET: www.sharcnet.ca) and Compute/Calcul Canada (www.computecanada.ca). Computations were made on the supercomputer Mammouth Parallele 2 from University of Sherbrooke, managed by Calcul Quebec and Compute Canada. The operation of this supercomputer is funded by the Canada Foundation for Innovation (CFI), the ministere de l'Economie, de la science et de l'innovation du Quebec (MESI), and the Fonds de recherche du Quebec-Nature et technologies (FRQ-NT). This work was supported in part by 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 Canada and by the Province of Ontario through the Ministry of Economic Development, Job Creation and Trade. A.E.B. thanks the Delaney Family for their generous financial support via the Delaney Family John A. Wheeler Chair at Perimeter Institute. A.E.B., P.T., and M.K. receive additional financial support from the Natural Sciences and Engineering Research Council of Canada through a Discovery Grant. R.G. receives additional support from the ERC synergy grant "BlackHoleCam: Imaging the Event Horizon of Black Holes" (grant No. 610058). We thank Jason Dexter for helpful comments and suggestions. We 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; Black Hole Initiative at Harvard University, through grants from the John Templeton Foundation and the Gordon and Betty Moore 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-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 PGC2018-098915-B-C21, 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 program 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 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. ; Peer reviewed