Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published July 10, 2020 | Published
Journal Article Open

THEMIS: A Parameter Estimation Framework for the Event Horizon Telescope

Broderick, Avery E. ORCID icon
Gold, Roman ORCID icon
Karami, Mansour ORCID icon
Preciado-López, Jorge A. ORCID icon
Tiede, Paul ORCID icon
Pu, Hung-Yi ORCID icon
Akiyama, Kazunori ORCID icon
Alberdi, Antxon ORCID icon
Alef, Walter
Asada, Keiichi
Azulay, Rebecca ORCID icon
Baczko, Anne-Kathrin ORCID icon
Baloković, Mislav ORCID icon
Barrett, John ORCID icon
Bintley, Dan
Blackburn, Lindy ORCID icon
Boland, Wilfred
Bouman, Katherine L. ORCID icon
Bower, Geoffrey C. ORCID icon
Bremer, Michael
Brinkerink, Christiaan D. ORCID icon
Brissenden, Roger ORCID icon
Britzen, Silke ORCID icon
Broguiere, Dominique
Bronzwaer, Thomas ORCID icon
Byun, Do-Young ORCID icon
Carlstrom, John E. ORCID icon
Chael, Andrew ORCID icon
Chatterjee, Shami ORCID icon
Chatterjee, Koushik ORCID icon
Chen, Ming-Tang ORCID icon
Chen, Yongjun ORCID icon
Cho, Ilje ORCID icon
Conway, John E. ORCID icon
Cordes, James M. ORCID icon
Crew, Geoffrey B. ORCID icon
Cui, Yuzhu ORCID icon
Davelaar, Jordy ORCID icon
De Laurentis, Mariafelicia ORCID icon
Deane, Roger ORCID icon
Dempsey, Jessica ORCID icon
Desvignes, Gregory ORCID icon
Doeleman, Sheperd S. ORCID icon
Eatough, Ralph P. ORCID icon
Falcke, Heino ORCID icon
Fish, Vincent L. ORCID icon
Fomalont, Ed ORCID icon
Fraga-Encinas, Raquel ORCID icon
Friberg, Per ORCID icon
Fromm, Christian M. ORCID icon
Galison, Peter ORCID icon
Gammie, Charles F. ORCID icon
García, Roberto ORCID icon
Gentaz, Olivier
Georgiev, Boris ORCID icon
Goddi, Ciriaco ORCID icon
Gómez, José L. ORCID icon
Gu, Minfeng ORCID icon
Gurwell, Mark ORCID icon
Hada, Kazuhiro ORCID icon
Hecht, Michael H.
Hesper, Ronald ORCID icon
Ho, Luis C. ORCID icon
Ho, Paul ORCID icon
Honma, Mareki ORCID icon
Huang, Chih-Wei L.
Huang, Lei ORCID icon
Hughes, David H.
Inoue, Makoto ORCID icon
Issaoun, Sara ORCID icon
James, David J. ORCID icon
Janssen, Michael ORCID icon
Jeter, Britton ORCID icon
Jiang, Wu ORCID icon
Jiménez-Rosales, Alejandra ORCID icon
Johnson, Michael D. ORCID icon
Jorstad, Svetlana ORCID icon
Jung, Taehyun ORCID icon
Karuppusamy, Ramesh ORCID icon
Kawashima, Tomohisa ORCID icon
Keating, Garrett K. ORCID icon
Kettenis, Mark ORCID icon
Kim, Jae-Young ORCID icon
Kim, Jongsoo ORCID icon
Kino, Motoki ORCID icon
Koay, Jun Yi ORCID icon
Koch, Patrick M. ORCID icon
Koyama, Shoko ORCID icon
Kramer, Michael ORCID icon
Kramer, Carsten ORCID icon
Krichbaum, Thomas P. ORCID icon
Kuo, Cheng-Yu ORCID icon
Lee, Sang-Sung ORCID icon
Li, Yan-Rong ORCID icon
Li, Zhiyuan ORCID icon
Lindqvist, Michael ORCID icon
Lico, Rocco ORCID icon
Liu, Kuo ORCID icon
Liuzzo, Elisabetta ORCID icon
Lo, Wen-Ping ORCID icon
Lobanov, Andrei P. ORCID icon
Loinard, Laurent ORCID icon
Lonsdale, Colin ORCID icon
Lu, Ru-Sen ORCID icon
MacDonald, Nicholas R. ORCID icon
Mao, Jirong ORCID icon
Marscher, Alan P. ORCID icon
Martí-Vidal, Iván ORCID icon
Matsushita, Satoki
Matthews, Lynn D. ORCID icon
Menten, Karl M. ORCID icon
Mizuno, Yosuke ORCID icon
Mizuno, Izumi ORCID icon
Moran, James M. ORCID icon
Moriyama, Kotaro ORCID icon
Moscibrodzka, Monika ORCID icon
Müller, Cornelia ORCID icon
Nagai, Hiroshi ORCID icon
Nagar, Neil M. ORCID icon
Nakamura, Masanori ORCID icon
Narayan, Ramesh ORCID icon
Narayanan, Gopal ORCID icon
Natarajan, Iniyan ORCID icon
Neri, Roberto ORCID icon
Ni, Chunchong ORCID icon
Noutsos, Aristeidis
Okino, Hiroki ORCID icon
Olivares, Héctor ORCID icon
Ortiz-León, Gisela N. ORCID icon
Oyama, Tomoaki ORCID icon
Palumbo, Daniel C. M. ORCID icon
Park, Jongho ORCID icon
Pen, Ue-Li ORCID icon
Pesce, Dominic W. ORCID icon
Piétu, Vincent
Plambeck, Richard ORCID icon
PopStefanija, Aleksandar
Porth, Oliver ORCID icon
Prather, Ben ORCID icon
Ramakrishnan, Venkatessh ORCID icon
Rao, Ramprasad ORCID icon
Rawlings, Mark G. ORCID icon
Raymond, Alexander W. ORCID icon
Rezzolla, Luciano ORCID icon
Ripperda, Bart ORCID icon
Roelofs, Freek ORCID icon
Rogers, Alan ORCID icon
Ros, Eduardo ORCID icon
Rose, Mel ORCID icon
Rottmann, Helge ORCID icon
Ruszczyk, Chet ORCID icon
Ryan, Benjamin R. ORCID icon
Rygl, Kazi L. J. ORCID icon
Sánchez, Salvador ORCID icon
Sánchez-Arguelles, David ORCID icon
Sasada, Mahito ORCID icon
Savolainen, Tuomas ORCID icon
Schloerb, F. Peter
Schuster, Karl-Friedrich ORCID icon
Shao, Lijing ORCID icon
Shen, Zhiqiang ORCID icon
Small, Des ORCID icon
Sohn, Bong Won ORCID icon
SooHoo, Jason ORCID icon
Tazaki, Fumie ORCID icon
Tilanus, Remo P. J. ORCID icon
Titus, Michael ORCID icon
Toma, Kenji ORCID icon
Torne, Pablo ORCID icon
Traianou, Efthalia ORCID icon
Trippe, Sascha ORCID icon
Tsuda, Shuichiro
van Bemmel, Ilse ORCID icon
van Langevelde, Huib Jan ORCID icon
van Rossum, Daniel R. ORCID icon
Wagner, Jan ORCID icon
Wardle, John ORCID icon
Weintroub, Jonathan ORCID icon
Wex, Norbert ORCID icon
Wharton, Robert ORCID icon
Wielgus, Maciek ORCID icon
Wong, George N. ORCID icon
Wu, Qingwen ORCID icon
Yoon, Doosoo ORCID icon
Young, André ORCID icon
Young, Ken ORCID icon
Younsi, Ziri ORCID icon
Yuan, Feng ORCID icon
Yuan, Ye-Fei ORCID icon
Zensus, J. Anton ORCID icon
Zhao, Guangyao ORCID icon
Zhao, Shan-Shan ORCID icon
Zhu, Ziyan ORCID icon
Event Horizon Telescope Collaboration

Abstract

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.

Additional Information

© 2020 The Author(s). Published by the American Astronomical Society. Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Received 2019 September 12; revised 2020 April 29; accepted 2020 May 7; published 2020 July 10. 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 Parallèle 2 from University of Sherbrooke, managed by Calcul Québec and Compute Canada. The operation of this supercomputer is funded by the Canada Foundation for Innovation (CFI), the ministère de l'Économie, de la science et de l'innovation du Québec (MESI), and the Fonds de recherche du Québec—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; Comisión Nacional de Investigación Científica y Tecnológica (CONICYT, Chile, via PIA ACT172033, Fondecyt 1171506, BASAL AFB-170002, ALMA-conicyt 31140007); Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico, projects 104497, 275201, 279006, 281692); the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Dirección General de Asuntos del Personal Académico–Universidad Nacional Autónoma de México (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 Economía 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 Astrofísica de Andalucía (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 Università 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.

Attached Files

Published - Broderick_2020_ApJ_897_139.pdf

Files

Broderick_2020_ApJ_897_139.pdf
Files (4.6 MB)
Name Size Download all
md5:8166ea7ce3a53585180a488fa60c2bfe
4.6 MB Preview Download

Additional details

Created:
August 22, 2023
Modified:
October 20, 2023