Я зашел на https://eventhorizontelescope.org/array и прочитал о десяти сайтах, перечисленных как часть EHT, у меня есть их скриншоты ниже.
Я сделал небольшой скрипт, также показанный ниже, с приблизительными координатами, и увидел, что Sgr A* находится под склонением около -29,5 градусов (южнее). Я вращал его вокруг Земли, вычислял скалярное произведение между векторами положения каждой точки и всеми нормалями, указывающими на Sgr A*, когда он вращается вокруг Земли.
Затем я нанес угол возвышения над горизонтом для всех десяти участков. Я был очень удивлен!
Я сделал ошибку? Если нет, то как эти десять мест объединяют свои данные для изображения Sgr A*?
class Site(object):
def __init__(self, name, lat, lon, alt):
self.name = str(name)
self.lat = rads * float(lat)
self.lon = rads * float(lon)
self.alt = km * float(alt)
(clat, slat), (clon, slon) = [[f(x) for f in (np.cos, np.sin)] for x in (self.lat, self.lon)]
self.norm = np.array([clon*clat, slon*clat, slat])
data = (('Northern Extended Millimeter Array', (44.634, 5.908, 2550.)),
('IRAM 30 meter telescope', (37.066, -3.393, 2850.)),
('The Greenland Telescope now near Thule Air Base', (76.531, -68.703, 10.)),
('Combined Array for Research in Millimeter-wave Astronomy CARMA', (37.280, -118.142, 2196.)),
('Kitt Peak National Observatory 12 meter Submillimeter Telescope (SMT)', (1.9583, -111.5967, 2096.)),
('Mt. Graham International Observatory 12 meter ALMA prototye', (32.701, -109.892, 3191.)),
('The Large Millimeter Telescope Alfonso Serrano', (18.985, -97.315, 4600.)),
('ALMA', (-22.971, -67.703, 4800.)),
('Caltech Submillimeter Observatory', (19.823, -155.476, 4140.)),
('South Pole Telescope', (-90.0, 0.0, 2800.)))
# https://eventhorizontelescope.org/array
# https://astronomy.stackexchange.com/questions/26413/math-behind-a-uv-plot-in-interferometry
datadict = dict(data)
import numpy as np
import matplotlib.pyplot as plt
from skyfield.api import Topos, Loader, EarthSatellite
from mpl_toolkits.mplot3d import Axes3D
halfpi, pi, twopi = [f*np.pi for f in (0.5, 1, 2)]
degs, rads = 180/pi, pi/180
km = 0.001
sites = [Site(a, *b) for a, b in datadict.items()]
Dec = rads * -29.5 # SgrA_star
cDec, sDec = [f(Dec) for f in (np.cos, np.sin)]
th = twopi * np.linspace(0, 1, 100001)[:-1]
cth, sth = [f(th) for f in (np.cos, np.sin)]
zth, oth = np.zeros_like(th), np.ones_like(th)
SgrA_star = np.stack([cDec*cth, cDec*sth, sDec*oth], axis=1)
for site in sites:
site.elev = np.arcsin(np.dot(SgrA_star, site.norm))
if True:
for site in sites:
plt.plot(degs*site.elev)
plt.plot(zth, '-k', linewidth=2)
plt.ylim(-90, 90)
plt.ylabel('elevation (deg)', fontsize=16)
plt.show()
Я не думаю, что необходимо, чтобы все телескопы одновременно наблюдали цель, если вы можете сделать предположение, что наблюдаемый вами источник изменяется только в масштабах времени, превышающих время, необходимое для того, чтобы источник был виден всеми телескопами. Конечно, максимальная базовая линия, которую вы можете получить, будет определяться самой удаленной парой одновременно наблюдающих телескопов, но это не мешает вам складывать вместе изображения, сделанные в другое время.
Вероятно, именно поэтому M87 является первым изображением, представленным EHT. Временная шкала значительной изменчивости вокруг этой (большой) черной дыры будет порядка нескольких раз. , что для черной дыры M87 составляет несколько дней. Получение изображения Sgr A* будет более сложным (или, по крайней мере, более размытым), поскольку временная шкала изменчивости составляет минуты.
Похоже (из документа по сокращению данных ), что первые шаги включают обработку каждой базовой линии попарно (как я предложил выше), а затем объединение их по сети с любыми базовыми линиями, которые наблюдались в то время, чтобы получить данные. с удивительно (для меня) высоким временным разрешением 10 с.
Теперь, когда данные/изображения Sgr A* были опубликованы ( Акияма и др., 2022 г. ), появилась следующая информация:
Наблюдения EHT проводились с помощью восьми обсерваторий в шести местах: ALMA и Atacama Pathfinder Experiment (APEX) на Льяно-де-Чахнантор в Чили, Большой миллиметровый телескоп Альфонсо Серрано (LMT) на вулкане Сьерра-Негра в Мексике, телескоп Джеймса Клерка Максвелла. (JCMT) и субмиллиметровая решетка (SMA) на Маунакеа на Гавайях, 30-метровый телескоп IRAM на Пико Велета (PV) в Испании, субмиллиметровый телескоп (SMT) на горе Грэм в Аризоне и телескоп Южного полюса (SPT) ) в Антарктиде.
В недавно объявленных наблюдениях M87 использовались только 8 из 10 сайтов. Этот документ включает эту диаграмму. Кажется, даже прибор Южного полюса использовался только для калибровки (пунктирные линии).
- Sgr A* никогда не находится над горизонтом для Гренландского телескопа!
Что ж, Sgr A* — не единственная цель EHT; на веб-сайте упоминается, что он также изучает другие объекты, такие как галактика M87 в Деве , которая при склонении +12 ° видна из Гренландии.
- Sgr A * никогда не виден одновременно даже во все остальные девять телескопов.
Интерферометрия со сверхдлинной базой просто объединяет результаты телескопов, которые могут (и были) направлены на конкретную цель в определенное время (плюс-минус несколько наносекунд, потому что некоторые телескопы могут быть дальше от цели, чем другие, и свет требуется немного больше или меньше времени, чтобы добраться до них). Таким образом, он никогда не будет использовать весь потенциал всех телескопов массива одновременно.
В месте расположения коррелятора данные воспроизводятся. Время воспроизведения регулируется в соответствии с сигналами атомных часов на (лентах / дисках / оптоволоконном сигнале) и расчетным временем прихода радиосигнала на каждый из телескопов. Обычно проверяется диапазон времени воспроизведения в диапазоне наносекунд, пока не будет найдено правильное время.
ПрофРоб
ооо