способ определения скорости распространения гравитационных волн

Формула изобретения

Способ определения скорости распространения гравитационных волн, основанный на измерении разности фаз гравитационного сигнала от источника излучения, принятого в двух пространственно разнесенных точках, отличающийся тем, что при годовом орбитальном обращении Земли с гравитационно-волновым детектором последовательно в произвольно фиксированные моменты времени t₁ и t₂ измеряют фазы гравитационного сигнала ₁ и ₂ от двойных астрофизических объектов, вычисляют непосредственно измеренную разность фаз _O, расчетную разность фаз _C, определяемую в предположении неподвижности Земли относительно центра масс Солнечной системы, и вычисляют скорость распространения гравитационной волны V по формуле



где _g - частота принимаемого гравитационно-волнового сигнала, D - изменение расстояния до источника гравитационно-волнового излучения, обусловленное орбитальным движением Земли, _O= ₂- ₁ - наблюденная разность фаз, _C - расчетная разность фаз.

Описание изобретения к патенту

Изобретение относится к области гравитационно-волновой физики и может быть использовано для определения скорости распространения гравитационных волн.

Известен способ [1] определения скорости электромагнитных волн, основанный на измерении разности моментов регистрации сигнала от одного и того же источника излучения, принятого в двух пространственно разнесенных точках, который наиболее близок к заявляемому способу и может быть принят в качестве ПРОТОТИПА. Способ по прототипу предполагает регистрацию сигнала от определенного источника излучения в двух разнесенных в пространстве точках и по разности между измеренными моментами времени и известной разности расстояний между источником и точками регистрации вычисление скорости распространения электромагнитной волны. Точность результата зависит от расстояния между точками регистрации и становится приемлемой при астрономических величинах порядка расстояний между планетами Солнечной системы. Это предполагает наличие станций регистрации не только на Земле, но и на других планетах или космических аппаратах, что и является существенным недостатком способа.

Задача, на решение которой направлено заявляемое изобретение, заключается в разработке способа определения скорости распространения гравитационных волн при использовании одного гравитационно-волнового детектора вместо двух, что возможно при регистрации сигналов от источника стабильного периодического излучения. Это позволит существенно сократить расходы на создание и обеспечение функционирования внеземных гравитационно-волновых детекторов. Кроме того, работа с периодическими источниками дает возможность повышать точность результатов путем накопления данных в течение произвольных интервалов измерений.

Сущность изобретения заключается в том, что скорость гравитационной волны определяется из измерений фазы сигнала, регистрируемого.от двойного релятивистского астрофизического объекта низкочастотного периодического гравитационного излучения детектором, размещенным на Земле и перемещающимся вместе с ней в процессе ее орбитального движения вокруг центра масс Солнечной системы, фазы гравитационного сигнала ₁ и ₂ от двойных астрофизических объектов измеряют в произвольно фиксированные моменты времени t₁ и t ₂ вычисляют разность фаз _o= ₂- ₁, и расчетную разность фаз _c, определяемую в предположении неподвижности Земли относительно центра масс Солнечной системы, и вычисляют скорость распространения гравитационной волны по формуле

где _g - частота принимаемого гравитационно-волнового сигнала, D - изменение расстояния до источника гравитационно-волнового излучения за время t₂-t₁, обусловленное орбитальным движением Земли.

В отличие от известного способа, где разность моментов регистрации принятого от источника (одного) импульсного сигнала измеряют в двух разнесенных точках, в заявляемом изобретении разность фаз периодического сигнала от источника, принятого одним детектором, измеряют в различных точках пространства при его орбитальном движении.

При разработке заявляемого способа по определению скорости распространения гравитационных волн предполагалось, что безразмерная амплитуда на несколько порядков ниже уровня шумов. Отсюда вытекает необходимость их достаточно длительного накопления до требуемого уровня соотношения сигнал/шум. Это приводит к погрешностям определения измеренных значений фазы ₁ и ₂, из которых вытекает ошибка определения наблюдаемой разности фаз _O. В предположении погрешности ошибок _i, заданной техническими характеристиками ГВ-детектора и принятой методикой выделения слабых сигналов из шумов [2, 3], влияние ошибки может быть сведено к минимуму оптимизацией выбора источника ГВ-излучения (астрофизический объект должен лежать вблизи плоскости эклиптики) и времени наблюдений (измерения следует выполнять. с интервалом в полгода, когда объект лежит на одной прямой линии с Солнцем).

Представим вывод основных формул, на которых основывается заявляемый способ.

Запишем фазу сигнала, локализованного в начале координатной системы, которой является барицентр (т.е. центр масс) Солнечной системы:

где _g - случайная начальная фаза ГВ-сигнала, _g - его частота:

где T - период обращения двойной системы.

Предположим, что выполнены два измерения фазы гравитационной волны в моменты времени t₁ и t₂ в точках 1 и 2, разнесенных в пространстве. Тогда фазы сигнала, принимаемого в моменты времени t₁ и t₂ в точках 1 и 2 соответственно, запишутся так:

где V - скорость распространения гравитационных волн, а величины s₁ и s₂ определяют, насколько детектор ближе к источнику, чем барицентр; они вычисляются как скалярные произведения векторов положения детектора в точках 1 и 2 и вектора направления на источник.

Назовем измеренной разность фаз _o, которая получится как разница измерений:

где D=s₂-s₁ введена для обозначения разности расстояний точек 1 и 2 от источника.

Назовем расчетной разность фаз _c, определяемую формулой:

Расчетная разность фаз может быть интерпретирована двояко. Во-первых, это разность фаз, которая наблюдалась бы неподвижным детектором, находящимся в барицентре Солнечной системы. Во-вторых, это разность фаз, которая наблюдалась бы в случае мгновенного распространения гравитационной волны, т.е. если бы величина V была бы равна бесконечности или была бы столь велика, что наши средства измерения не позволили бы отличить ее от бесконечной.

Измеренная и расчетная разности фаз отличаются на величину, определяемую скоростью гравитационной волны:

Из этого выражения выведем формулу для вычисления скорости распространения гравитационной волны:

Приведем также формулы для оценки точности получаемого результата:

где m_V - погрешность определения скорости гравитационной волны, m - погрешность определения разности фаз.

Определение скорости распространения гравитационной волны рассмотрим на примере низкочастотного периодического ГВ-источника PSR 1756-2251, представляющего собой двойную систему с периодом обращения Т=0,3196339 суток, расположенного вблизи плоскости эклиптики (эклиптическая широта =0,5707153°). Расчетная амплитуда ГВ-сигнала h ₀=9,3·10^-23. Рассчитаем частоту гравитационного сигнала:

Выберем систему координат таким образом, чтобы ось Х лежала в плоскости эклиптики и была направлена в сторону источника. При этом направляющий вектор источника будет (cos , 0, sin ). Выполним первое измерение фазы гравитационной волны ₁ в точке земной орбиты, наиболее удаленной от источника, а второе измерение ₂ - в точке, наиболее близкой к нему. В первом случае вектор координат ГВ-детектора будет близок к (-a, 0, 0), во втором (а, 0, 0). Здесь а=1,4959787·10¹¹ м - радиус земной орбиты.

Пусть измерены следующие величины (в долях циклов):

₁=0,924; ₂=0,716

и сосчитаны целые циклы (когда =0) в количестве 1142.

Шаг первый: вычислим разность фаз _o:

_o= ₂- ₁=1143+0,716-0,924=1142,792 цикла 7180,374 рад.

Шаг второй: вычислим расчетную разность фаз _c. В нашем модельном случае моменты времени t₁ и t₂ разделяет полугодовой интервал времени:

где 3,1558196·10⁷ - период обращения Земли в секундах. Подставим эту величину в формулу:

_c= _g(t₂-t₁)=4,550334·10 ^-4·1,5779098·10⁷ 7180,017 рад.

Шаг третий: вычислим расхождение между фазами:

_o- _c=7180,374-7180,017 0,357 рад.

Далее вычислим разность расстояний точек от источника:

s₁=-acos =-1,4959787·10¹¹ м·0.9999504 -1,4959045·10¹¹ м;

s ₂=+acos =+1,4959787·10¹¹ м·0.9999504 +1,4959045·10¹¹ м;

D=s₂-s₁ 2,9918090·10¹¹ м.

В реальных условиях координаты Земли и ГВ-детектора на моменты измерений будут вычисляться в соответствии с точными эфемеридами теории DE405/LE405 [4].

Шаг четвертый: вычислим скорость гравитационной волны:

Оценим точность полученного результата. Инструментальная точность измерения фазы может быть достигнута порядка 0,01 цикла, или 0,06 рад. Тогда относительная погрешность измеренной скорости составит:

Абсолютная величина погрешности равна

Результат записывается как V=3,81±0,64·10 ⁸ м/с.

Точность может быть повышена путем накопления повторных измерений.

Кроме того, для уменьшения погрешности результата можно:

- увеличить базу разноса точек, увеличив размер орбиты ГВ-детектора;

- подбирать для работы источники ГВ-излучения с наименьшими периодами.

Таким образом, заявляемый способ выгодно отличается от прототипа возможностью определения скорости распространения гравитационных волн путем сравнения разности фаз _o, между сигналами, принимаемыми от ГВ-источников последовательно при перемещении ГВ-детектора по орбите в течение произвольно долгих интервалов измерений, и расчетной разности фаз _c, определяемой в предположении неподвижности Земли относительно центра масс Солнечной системы.

Источники информации

1. Пестряков В.Б. Фазовые радиотехнические системы. - М.: Сов. Радио, 1968, 466с. (прототип).

2. Даишев Р.А., Левин С.Ф., Мурзаханов З.Г., Скочилов А.Ф. // Измерительная техника, 2001, № 2, с.6.

3. Константинов В.Б., Мурзаханов З.Г. // Измерительная техника, 2001, № 3, с.6.

4. Standish, E.M.: 1998, "JPL Planetary and Lunar Ephemerides, DE405/LE405", JPL IOM 312.F-98-048.