Предложите экспериментальную стратегию для поиска непрерывных гравитационных волн от быстро вращающегося нейтронного звезды: оцените порядок величины требуемой амплитуды деформации (ellipticity), частоты и чувствительности детектора, и какие астрономические наблюдения (радио, рентген) помогут уточнить цель и повысить шанс детекции
Кратко — стратегия + порядок величин + какие наблюдения нужны. 1) Физика и ключевая формула - Для триаксиальной деформации (излучение на частоте fgw=2frotf_{\rm gw}=2f_{\rm rot}fgw=2frot) амплитуда волн (на полюсе) даётся стандартно h0=4π2Gc4I ϵ fgw2d,
h_0=\frac{4\pi^2 G}{c^4}\frac{I\,\epsilon\,f_{\rm gw}^2}{d}, h0=c44π2GdIϵfgw2,
где III — момент инерции (типично I∼1038 kg m2I\sim10^{38}\,\mathrm{kg\,m^2}I∼1038kgm2), ϵ\epsilonϵ — ellipticity, ddd — расстояние. - Удобный упрощённый масштабный вид (для I=1038 kg m2I=10^{38}\,\mathrm{kg\,m^2}I=1038kgm2, d=1 kpcd=1\,\mathrm{kpc}d=1kpc): h0≈1.05×10−26(ϵ10−6)(fgw100 Hz)2.
h_0\approx 1.05\times10^{-26}\left(\frac{\epsilon}{10^{-6}}\right)\left(\frac{f_{\rm gw}}{100\,\mathrm{Hz}}\right)^2. h0≈1.05×10−26(10−6ϵ)(100Hzfgw)2. 2) Примеры чисел (порядки величин) - Миллисекундные пульсары: frot∼200 − 700f_{\rm rot}\sim 200\!-\!700frot∼200−700 Hz ⇒fgw∼400 − 1400\Rightarrow f_{\rm gw}\sim 400\!-\!1400⇒fgw∼400−1400 Hz. Молодые пульсары: frot∼10 − 100f_{\rm rot}\sim 10\!-\!100frot∼10−100 Hz ⇒fgw∼20 − 200\Rightarrow f_{\rm gw}\sim 20\!-\!200⇒fgw∼20−200 Hz. Хорошая чувствительность детекторов — в диапазоне 20 − 150020\!-\!150020−1500 Hz. - Типичная амплитуда для d=1d=1d=1 kpc: - ϵ=10−6, fgw=500\epsilon=10^{-6},\,f_{\rm gw}=500ϵ=10−6,fgw=500 Hz ⇒h0∼3×10−25\Rightarrow h_0\sim 3\times10^{-25}⇒h0∼3×10−25. - ϵ=10−7, fgw=500\epsilon=10^{-7},\,f_{\rm gw}=500ϵ=10−7,fgw=500 Hz ⇒h0∼3×10−26\Rightarrow h_0\sim 3\times10^{-26}⇒h0∼3×10−26. - Для fgw=100f_{\rm gw}=100fgw=100 Hz и ϵ=10−6\epsilon=10^{-6}ϵ=10−6: h0∼10−26h_0\sim10^{-26}h0∼10−26. - Ожидаемые правдоподобные ϵ\epsilonϵ: теоретически ϵ∼10−9 − 10−5\epsilon\sim10^{-9}\!-\!10^{-5}ϵ∼10−9−10−5 для обычных нейтронных звёзд (крышка/пластичность/магн. деформации). Экстремальные сценарии (твёрдая кварковая материя, внутренние поля B≳1015 − 1016B\gtrsim10^{15}\!-\!10^{16}B≳1015−1016 G) могут дать ϵ∼10−5 − 10−3\epsilon\sim10^{-5}\!-\!10^{-3}ϵ∼10−5−10−3. 3) Чувствительность детектора и интеграция - Для целенаправленного (targeted) поиска, зная фазу и орбитальные параметры, минимально обнаруживаемая h0h_0h0 при когерентной обработке масштабается как ∼Sn/Tobs\sim\sqrt{S_n/T_{\rm obs}}∼Sn/Tobs (множители зависят от метода, гауссовского шума и числа параметров). Практическая оценка: h0,min∼κSnTobs,
h_{0,\min}\sim \kappa\sqrt{\frac{S_n}{T_{\rm obs}}}, h0,min∼κTobsSn,
где κ∼10\kappa\sim 10κ∼10 (консервативно), SnS_nSn — спектральная плотность шума. Для Tobs∼1T_{\rm obs}\sim 1Tobs∼1 год и современных детекторов (Advanced LIGO/Virgo: Sn∼10−23 − 10−24/Hz \sqrt{S_n}\sim 10^{-23}\!-\!10^{-24}/\sqrt{\mathrm{Hz}}Sn∼10−23−10−24/Hz в лучшем диапазоне) достижимы пределы h0∼10−26 − 10−25h_0\sim10^{-26}\!-\!10^{-25}h0∼10−26−10−25 по некоторым целям. - Следствие: чтобы достать ϵ∼10−6\epsilon\sim10^{-6}ϵ∼10−6 для источника на ∼\sim∼кпк нужен h0∼10−25 − 10−24h_0\sim10^{-25}\!-\!10^{-24}h0∼10−25−10−24 (входящая в возможности текущих таргетированных поисков). Для ϵ∼10−8\epsilon\sim10^{-8}ϵ∼10−8 потребуется чувствительность лучше на 2 порядка — цель для 3‑го поколения детекторов (Einstein Telescope, Cosmic Explorer). 4) Приоритеты в поисковой стратегии (экспериментально) - Targeted searches (наилучшая чувствительность): известные радиопульсары / гамма‑пульсары с точной эпимеридой. Делать полностью когерентные поиски на всё доступное TobsT_{\rm obs}Tobs (месяцы—годы), комбинировать данные нескольких детекторов. - Directed searches: объекты с известным положением (пульсар в SNR, не выявленный в радио), но неопределённая частота — использовать полукохерентные и иерархические методы (StackSlide, Hough, PowerFlux). - Searches for accreting NS (Sco X-1 и другие LMXB): использовать методы учёта орбитальных модулей (sideband, CrossCorr), координировать с X‑ray таймингом. - All‑sky blind searches (неизвестные источники): распределённые вычисления (Einstein@Home), многоступенчатая полукохерентная стратегия — меньше чувствительность, но покрытие большой объём параметров. 5) Какие астрономические наблюдения существенно повышают шанс детекции - Радио‑тайминг (pulsar timing): точная частота, f˙\dot ff˙, орбитальные параметры, наблюдение за глитчами и тайминг‑шумом — критично для когерентной интеграции. Парраллакс/DM → расстояние ddd. - Рентген/Никер (NICER, XMM‑Newton): тайминговые эпимериды для радиотихих или аккрецирующих пульсаров; профиль импульса → ориентация оси (углы инклинации) помогает интерпретировать амплитуду. - Гамма (Fermi LAT): открытие радио‑тихих пульсаров и получение эпимерид для таргетированных поисков. - Оптические/инфрачервонные и SNR‑ассоциации: расстояние и возраст источника, окружение (плотность, магнитные свойства) — важны для оценки возможной ϵ\epsilonϵ. - Измерения магнитного поля (X‑ray спектроскопия, поляриметрия) → оценка возможной магнитной деформации ϵB\epsilon_BϵB. - Непрерывный мониторинг для отслеживания глитчей, активаций, вариабельности (они нарушают фазовую модель и требуют адаптивных анализов). 6) Практические рекомендации (коротко) - Приоритизировать близкие и быстрые пульсары (малые ddd, большие fff) и объекты с высоким спин‑дауном (spin‑down limit даёт верхнюю границу сигналов). - Получать и поддерживать высококачественные эпимериды (радио/гамма/рентген) во время GW‑наблюдений. - Для аккрецирующих систем — синхронизироваться с X‑ray наблюдениями для ограничения орбиты и частоты. - Использовать многодетекторную когерентную обработку, затем иерархические полукохерентные поиски для расширения объёма параметров; глубокие follow‑up кандидатов когерентной интеграцией. - Оценивать результаты в терминах ϵ\epsilonϵ и spin‑down limit: если поставленный предел h0h_0h0 меньше spin‑down limit — это физически интересно. 7) Формулы для перепроверки (полезно) - Обратная формула для ellipticity: ϵ=c44π2Gd h0Ifgw2.
\epsilon=\frac{c^4}{4\pi^2 G}\frac{d\,h_0}{I f_{\rm gw}^2}. ϵ=4π2Gc4Ifgw2dh0.
- Spin‑down limit (максимум, если весь спин‑даун идёт в ГВ): h0sd=1d5GI2c3∣f˙rot∣frot.
h_0^{\rm sd}=\frac{1}{d}\sqrt{\frac{5G I}{2c^3}\frac{|\dot f_{\rm rot}|}{f_{\rm rot}}}. h0sd=d12c35GIfrot∣f˙rot∣. Итого: оптимальная экспериментальная стратегия — таргетированные когерентные поиски по ближайшим/быстрым пульсарам с длительной интеграцией и точными эпимеридными наблюдениями (радио/рентген/гамма), параллельно — направленные и all‑sky полукохерентные поиски. Для доступа к ϵ≲10−7\epsilon\lesssim10^{-7}ϵ≲10−7 потребуется чувствительность уровней h0∼10−26h_0\sim10^{-26}h0∼10−26 и дальше — 3‑го поколения детекторов.
1) Физика и ключевая формула
- Для триаксиальной деформации (излучение на частоте fgw=2frotf_{\rm gw}=2f_{\rm rot}fgw =2frot ) амплитуда волн (на полюсе) даётся стандартно
h0=4π2Gc4I ϵ fgw2d, h_0=\frac{4\pi^2 G}{c^4}\frac{I\,\epsilon\,f_{\rm gw}^2}{d},
h0 =c44π2G dIϵfgw2 , где III — момент инерции (типично I∼1038 kg m2I\sim10^{38}\,\mathrm{kg\,m^2}I∼1038kgm2), ϵ\epsilonϵ — ellipticity, ddd — расстояние.
- Удобный упрощённый масштабный вид (для I=1038 kg m2I=10^{38}\,\mathrm{kg\,m^2}I=1038kgm2, d=1 kpcd=1\,\mathrm{kpc}d=1kpc):
h0≈1.05×10−26(ϵ10−6)(fgw100 Hz)2. h_0\approx 1.05\times10^{-26}\left(\frac{\epsilon}{10^{-6}}\right)\left(\frac{f_{\rm gw}}{100\,\mathrm{Hz}}\right)^2.
h0 ≈1.05×10−26(10−6ϵ )(100Hzfgw )2.
2) Примеры чисел (порядки величин)
- Миллисекундные пульсары: frot∼200 − 700f_{\rm rot}\sim 200\!-\!700frot ∼200−700 Hz ⇒fgw∼400 − 1400\Rightarrow f_{\rm gw}\sim 400\!-\!1400⇒fgw ∼400−1400 Hz. Молодые пульсары: frot∼10 − 100f_{\rm rot}\sim 10\!-\!100frot ∼10−100 Hz ⇒fgw∼20 − 200\Rightarrow f_{\rm gw}\sim 20\!-\!200⇒fgw ∼20−200 Hz. Хорошая чувствительность детекторов — в диапазоне 20 − 150020\!-\!150020−1500 Hz.
- Типичная амплитуда для d=1d=1d=1 kpc:
- ϵ=10−6, fgw=500\epsilon=10^{-6},\,f_{\rm gw}=500ϵ=10−6,fgw =500 Hz ⇒h0∼3×10−25\Rightarrow h_0\sim 3\times10^{-25}⇒h0 ∼3×10−25.
- ϵ=10−7, fgw=500\epsilon=10^{-7},\,f_{\rm gw}=500ϵ=10−7,fgw =500 Hz ⇒h0∼3×10−26\Rightarrow h_0\sim 3\times10^{-26}⇒h0 ∼3×10−26.
- Для fgw=100f_{\rm gw}=100fgw =100 Hz и ϵ=10−6\epsilon=10^{-6}ϵ=10−6: h0∼10−26h_0\sim10^{-26}h0 ∼10−26.
- Ожидаемые правдоподобные ϵ\epsilonϵ: теоретически ϵ∼10−9 − 10−5\epsilon\sim10^{-9}\!-\!10^{-5}ϵ∼10−9−10−5 для обычных нейтронных звёзд (крышка/пластичность/магн. деформации). Экстремальные сценарии (твёрдая кварковая материя, внутренние поля B≳1015 − 1016B\gtrsim10^{15}\!-\!10^{16}B≳1015−1016 G) могут дать ϵ∼10−5 − 10−3\epsilon\sim10^{-5}\!-\!10^{-3}ϵ∼10−5−10−3.
3) Чувствительность детектора и интеграция
- Для целенаправленного (targeted) поиска, зная фазу и орбитальные параметры, минимально обнаруживаемая h0h_0h0 при когерентной обработке масштабается как ∼Sn/Tobs\sim\sqrt{S_n/T_{\rm obs}}∼Sn /Tobs (множители зависят от метода, гауссовского шума и числа параметров). Практическая оценка:
h0,min∼κSnTobs, h_{0,\min}\sim \kappa\sqrt{\frac{S_n}{T_{\rm obs}}},
h0,min ∼κTobs Sn , где κ∼10\kappa\sim 10κ∼10 (консервативно), SnS_nSn — спектральная плотность шума. Для Tobs∼1T_{\rm obs}\sim 1Tobs ∼1 год и современных детекторов (Advanced LIGO/Virgo: Sn∼10−23 − 10−24/Hz \sqrt{S_n}\sim 10^{-23}\!-\!10^{-24}/\sqrt{\mathrm{Hz}}Sn ∼10−23−10−24/Hz в лучшем диапазоне) достижимы пределы h0∼10−26 − 10−25h_0\sim10^{-26}\!-\!10^{-25}h0 ∼10−26−10−25 по некоторым целям.
- Следствие: чтобы достать ϵ∼10−6\epsilon\sim10^{-6}ϵ∼10−6 для источника на ∼\sim∼кпк нужен h0∼10−25 − 10−24h_0\sim10^{-25}\!-\!10^{-24}h0 ∼10−25−10−24 (входящая в возможности текущих таргетированных поисков). Для ϵ∼10−8\epsilon\sim10^{-8}ϵ∼10−8 потребуется чувствительность лучше на 2 порядка — цель для 3‑го поколения детекторов (Einstein Telescope, Cosmic Explorer).
4) Приоритеты в поисковой стратегии (экспериментально)
- Targeted searches (наилучшая чувствительность): известные радиопульсары / гамма‑пульсары с точной эпимеридой. Делать полностью когерентные поиски на всё доступное TobsT_{\rm obs}Tobs (месяцы—годы), комбинировать данные нескольких детекторов.
- Directed searches: объекты с известным положением (пульсар в SNR, не выявленный в радио), но неопределённая частота — использовать полукохерентные и иерархические методы (StackSlide, Hough, PowerFlux).
- Searches for accreting NS (Sco X-1 и другие LMXB): использовать методы учёта орбитальных модулей (sideband, CrossCorr), координировать с X‑ray таймингом.
- All‑sky blind searches (неизвестные источники): распределённые вычисления (Einstein@Home), многоступенчатая полукохерентная стратегия — меньше чувствительность, но покрытие большой объём параметров.
5) Какие астрономические наблюдения существенно повышают шанс детекции
- Радио‑тайминг (pulsar timing): точная частота, f˙\dot ff˙ , орбитальные параметры, наблюдение за глитчами и тайминг‑шумом — критично для когерентной интеграции. Парраллакс/DM → расстояние ddd.
- Рентген/Никер (NICER, XMM‑Newton): тайминговые эпимериды для радиотихих или аккрецирующих пульсаров; профиль импульса → ориентация оси (углы инклинации) помогает интерпретировать амплитуду.
- Гамма (Fermi LAT): открытие радио‑тихих пульсаров и получение эпимерид для таргетированных поисков.
- Оптические/инфрачервонные и SNR‑ассоциации: расстояние и возраст источника, окружение (плотность, магнитные свойства) — важны для оценки возможной ϵ\epsilonϵ.
- Измерения магнитного поля (X‑ray спектроскопия, поляриметрия) → оценка возможной магнитной деформации ϵB\epsilon_BϵB .
- Непрерывный мониторинг для отслеживания глитчей, активаций, вариабельности (они нарушают фазовую модель и требуют адаптивных анализов).
6) Практические рекомендации (коротко)
- Приоритизировать близкие и быстрые пульсары (малые ddd, большие fff) и объекты с высоким спин‑дауном (spin‑down limit даёт верхнюю границу сигналов).
- Получать и поддерживать высококачественные эпимериды (радио/гамма/рентген) во время GW‑наблюдений.
- Для аккрецирующих систем — синхронизироваться с X‑ray наблюдениями для ограничения орбиты и частоты.
- Использовать многодетекторную когерентную обработку, затем иерархические полукохерентные поиски для расширения объёма параметров; глубокие follow‑up кандидатов когерентной интеграцией.
- Оценивать результаты в терминах ϵ\epsilonϵ и spin‑down limit: если поставленный предел h0h_0h0 меньше spin‑down limit — это физически интересно.
7) Формулы для перепроверки (полезно)
- Обратная формула для ellipticity:
ϵ=c44π2Gd h0Ifgw2. \epsilon=\frac{c^4}{4\pi^2 G}\frac{d\,h_0}{I f_{\rm gw}^2}.
ϵ=4π2Gc4 Ifgw2 dh0 . - Spin‑down limit (максимум, если весь спин‑даун идёт в ГВ):
h0sd=1d5GI2c3∣f˙rot∣frot. h_0^{\rm sd}=\frac{1}{d}\sqrt{\frac{5G I}{2c^3}\frac{|\dot f_{\rm rot}|}{f_{\rm rot}}}.
h0sd =d1 2c35GI frot ∣f˙ rot ∣ .
Итого: оптимальная экспериментальная стратегия — таргетированные когерентные поиски по ближайшим/быстрым пульсарам с длительной интеграцией и точными эпимеридными наблюдениями (радио/рентген/гамма), параллельно — направленные и all‑sky полукохерентные поиски. Для доступа к ϵ≲10−7\epsilon\lesssim10^{-7}ϵ≲10−7 потребуется чувствительность уровней h0∼10−26h_0\sim10^{-26}h0 ∼10−26 и дальше — 3‑го поколения детекторов.