Constraints on massive graviton dark matter from precision pulsar timing and astrometry Konstantin POSTNOV (Sternberg Astronomical Institute) Collaborators: - презентация

1 Constraints on massive graviton dark matter from precision pulsar timing and astrometry Konstantin POSTNOV (Sternberg Astronomical Institute) Collaborators: Maxim Pshirkov (PRAO Lebedev Phyical Institute), Artyom Tuntsov (SAI), Aleksandr Polnarev (QMC London UK), Deepak Baskaran (Cardiff UK) QUARKS-2008

2 Plan Pulsars as GW detectors Observational constraints on massive graviton CDM Surfing effect of massive gravitons and limits on their propagation speed

3 Pulsars as GW detectors Gravitational waves (1/2) General Relativity GW propagation velocity in empty space is с: Along axis z: Polarizatrion tensor have two non-zero components Monochromatic transverse GW has two polarizatios (GR) &

4 GW energy density (monochromatic plane): Stochastic isotropic background: Is the critical density Or: Pulsars as GW detectors Gravitational waves (1/2)

5 Pulsars as monochromatic GW detectors Monochromatic GW (1/3) z x y PSR GW changes the observed pulsar frequency (Sazhin (1978), Detweiler (1979)) In GR interaction is independent of distance (if ) – no secular increase ~D. Is the GW polarization vector

6 Pulsars as GW detectors Monochromatic GW (2/3) Variation of the observed frequency results in time residuals in time of arrival (TOA): h Maximum sensitivity at frequencies ~ 1/T obs Longer GWs also contribute to the observed Pulsar period and its derivative 1/T obs 1/T samp 1/T int T obs ~ 10 years T samp ~ 10 days T int ~ 1 hour

7 h In 2003 periodic motions in 3C66b were explained by binary SMBH (Sudou et al., 2003)-80 Mpc, 1.5x10 10 M Timing of PSR B rejected this possibility (Jenet et al., 2004) Pulsars as GW detectors Monochromatic GW (3/3)

8 Pulsars as GW detectors Stochastic GWB (1/3) RMS of TOA residuals is (Detweiler, 1979): RMS of TOA residuals depend on GW energy density For flat GW spectrum of width Δf~f centered at f - the critical density

9 For arbitrary GWB («red noise»): Kaspi, Taylor, Ryba, 1994 Pulsars as GW detectors Stochastic GWB (2/3)

10 Pair correlation of the TOA residuals for 20 pulsars (simulation, R,Manchester, 2007 ) GW noise is the same for all pulsars It is advantageous to observe ensemble of pulsars and correlate rms of TOA residuals between each pair of pulsars Pulsars as GW detectors Stochastic GWB (3/3)

11 Pulsars as GW detectors Present limits and prospects (Manchester, 2007 – arXiv: v2)

12 Tests in Solar systems Doppler tracking (1/2) Estabrook & Wahlquist, 1975, principle similar to pulsar timing Best current limits: Cassini mission, Hz(Armstrong et al. 2003)

13 Solar system tests Doppler tracking (2/2) Reynaud et al Future projects: Search for Anomalous Gravity using Atomic Sensors, SAGAS

14 Astrometric constraints A GW causes «drizzling» of visual position of a source on the sky (e.g, Kaiser&Jaffe, 1997): The observed quantity is the arc length between two sources Ψ: In the presence of a GW sources on the sky would oscillate w.r.t. to their true position with amplitude h. Modern ICRF precision (~100 μas) constrain low- frequency GWB: h

15 Содержание Пульсары как детектор ГВ фона гравитационные волны отклик на монохроматическую волну отклик на стохастический спектр Теории с массивным гравитоном Предпосылки Наблюдаемые проявления Наблюдательная проверка теорий с массивным гравитоном Surfing effect

16 Theories with massive gravitons Massive gravity ( Rubakov 2004, Dubovsky 2004) with spontaneous Lorentz braking (Rubakov & Tinyakov arXiv: for a review) Healthy theory: no ghosts, no vDVZ discontinuity, no strong coupling at low scale Interesting phenomenology: DE-like term in Fridmann equations + possibility to produce massive gravitons in the early Universe copiously enough to explain all of CDM (Dubovsky, Tinyakov & Tkachev 2005) Taking graviton mass < (10 15 cm) -1 (binary PSR constraints) and assuming all galactic CDM due to massive gravitons leads to a strong almost monochromatic (Δf/f~10 -6 ) GW signal with amplitude

17 Содержание Пульсары как детектор ГВ фона гравитационные волны отклик на монохроматическую волну отклик на стохастический спектр Теории с массивным гравитоном Предпосылки Наблюдаемые проявления Наблюдательная проверка теорий с массивным гравитоном Surfing effect

18 Observational constraints: PTP08 Pulsar timing (1/2) Isotropic GW background affects pulsar timing GW amplitude can be constrained from rms residuals of TOA of even one pulsar Strong monochromatic signal (e.g. if all of galactic DM is due to massive gravitons, as in Dubovsky et al 2005) will manifest itself at frequencies < 1/T int (PSR integration time ~ 1-2 hrs) (PTP08):, Limit on the GW amplitude from the existing rms residuals of TOA of pulsars: 2008arXiv : Pshirkov, Tuntsov, Postnov

19 Constraints using existing rms TOA residuals (Manchester, 2007), PSR B Observational constraints: PTP08 Pulsar timing (2/2)

20 Observational constraints: BPPP08 «Surfing effect» (1/4) Unlike in GR, massive gravitons propagate with velocity less than c : Mass of the graviton is expressed through phenomenological parameter ε : Pulsar frequency change by massive GW:: arXiv: arXiv: : Baskaran, Polnarev, Pshirkov, Postnov

21 Observational constraints: BPPP08 «Surfing effect» (2/4) TOA residuals: Unlike GR, residuals seculary increase with distance to the source D ! Above results for a monochromatic GW can be generalized to stochastic GWB:

22 Observational constraints: BPPP08 «Surfing effect» (3/4) Response to any harmonics is known: The observed TOA residuals will be expressed through this «transfer function»:

23 Observational constraints: BPPP08 «Surfing effect» (4/4) R(k) depends on ε (term ) I. II. For example, power-law spectrum:

24 Observational constraints: BPPP08 «Surfing effect»: limits (1/5) Depending on ε, PSR timing put bounds on energy density of GWB: Or some combination of GW energy density and ε :

25 Observational constraints: BPPP08 «Surfing effect»: limits (2/5) For known GW amplitude, the parameter ε can be constrained: For theoretically motivated GWB from SMBH: or

26 Observational constraints: BPPP08 «Surfing effect»: limits (3/5)

27 Observational constraints: BPPP08 «Surfing effect»: limits (5/5) In terms of the graviton mass: From modern pulsar timing (Manchester 2007), which is by 3 orders of magnitude better than from Solar system bounds can be increased by one order with increasing observational time comparable to the future LISA constraints.

28 CONCLUSIONS Precise astronomical observations, especially pulsar timing, put strong bounds on massive graviton parameters: Cold massive gravitons cannot constitute all of the galactic dark matter

29 Спасибо за внимание!

30 Theories with massive gravitons (Тиняков 2007)

31 Теории с массивным гравитоном Модель (2/4) (Тиняков 2007)

32 Теории с массивным гравитоном Модель (4/4) (Тиняков 2007)

33 Теории с массивным гравитоном Модель (3/4) (Тиняков 2007)

34 Теории с массивным гравитоном Наблюдаемые проявления (1/4) (Тиняков 2007)

35 Теории с массивным гравитоном Наблюдаемые проявления (2/4) (Тиняков 2007) (H i – параметр Хаббла в инфл. эпоху)

36 Теории с массивным гравитоном Наблюдаемые проявления (3/4) (Тиняков 2007)

37 Принципы тайминга Одиночные пульсары(1/4) J J B J Stairs, 2003

38 Принципы тайминга Одиночные пульсары(2/4) Радиотелескоп РТ-64 КРАО (ТНА-1500 ОКБ МЭИ)

39 Принципы тайминга Одиночные пульсары(3/4) N-ый импульс от пульсара приходит на РТ в момент времени N. Редукция в барицентр Солнечной системы. Момент прихода в барицентр СС: Считается, что пульсар вращается по известным законам. Момент прихода N-го импульса связан с его номером, частотой вращения и её производными и может быть предсказан. В действительности, между наблюдаемыми моментами прихода N-го импульса и предсказанными значениями всегда существует разница-остаточные уклонения:

40 Принципы тайминга Одиночные пульсары(4/4) Уточнение параметров происходит по МНК. Минимизируются остаточные уклонения: -поправки к принятым значениям

41 Принципы тайминга Остаточные уклонения После процедуры остаются остаточные уклонения моментов прихода импульсов Остаточные уклонения пульсаров B и B ( , Аресибо), Kaspi, Taylor&Ryba(1994)

42 Принципы тайминга Двойные пульсары Движение в двойной системе описывается стандартными кеплеровскими параметрами: Период обращения:P b Проекция большой полуоси: Эксцентриситет:e Долгота периастра: ω Эпоха периастра: T 0 В сильных гравитационных полях появляются ПК-параметры ( и т.д. ) Все эти параметры могут быть найдены из тайминга (аналогично, МНК-методом)

43 Принципы тайминга Алгоритм 1.Наблюдения, вычисление моментов прихода импульсов пульсаров (МПИ) в барицентре Солнечной системы. 2.Вычисление теоретических значений МПИ с использованием модели хронометрирования. 3.Определение отклонения значений теоретических МПИ от наблюдаемых (расчет остаточных уклонений – ОУ МПИ). 4.Уточнение параметров модели хронометрирования (далее к п.3 до сходимости модели).