📻 DRAFT RESEARCH PROPOSAL · DEEP PULSAR TIMING · MEERKAT → SKA-MID

MeerKAT Long-Baseline Pulsar Timing: Definitive Gravitational Potential Map of Omega Centauri

A 5-year high-precision timing program targeting 20+ OC millisecond pulsars to constrain the central mass distribution at ~10⁻⁹ m/s² acceleration sensitivity — delivering Bayesian model discrimination between single IMBH and dark cluster hypotheses before LISA (~2035) · Working draft · April 2026

Relation to companion proposal: The MeerKAT/SKA 3-year PTA proposal is the discovery-phase program; this 5-year program is the precision-phase follow-up with expanded science goals. They can be submitted together as a staged program or independently.

1. Scientific Rationale

1.1 Quantitative Science Case

MeerKAT's timing precision for OC MSPs (flux densities ~10–100 µJy; integration ~1 hr each):

Timing precision (ToA): ~1–3 µs per 1-hr session (bright MSPs) Acceleration sensitivity from the quadratic-fit spin-down term [corrected 2026-09-25, P0-15.6: the previous line, δa = c·δσ_ToA/(P·T²), carried a spurious division by the spin period P]: δa_los = c · σ_ToA · k / T², with k = 2√(180/N) = √(720/N) (quadratic-fit coefficient uncertainty for N uniform ToAs; the linear-fit factor √(12/N) does NOT apply, because spin-down is quadratic) Bi-weekly cadence over T = 5 yr → N ≈ 130 ToAs per pulsar → k ≈ 2.35: δa_los ≈ (3×10⁸ m/s × 2×10⁻⁶ s × 2.35) / (1.575×10⁸ s)² ≈ 5.7×10⁻¹⁴ m/s² per pulsar (ToA-noise-limited only) But ToA noise is not the limiting term. The measured ṗ is the sum of intrinsic spin-down, cluster acceleration, and Shklovskii drift, and the intrinsic term is unknown per pulsar at the ~10⁻²⁰ s/s level, which maps to an acceleration uncertainty c×10⁻²⁰/P ≈ 10⁻⁹ m/s² (P = 3 ms). The usable sensitivity is therefore ~10⁻⁹–10⁻¹⁰ m/s², four orders of magnitude above the ToA-noise floor. Expected acceleration from 30,000 M☉ IMBH at r = 0.3 pc: a_los ≈ GM/r² × cos(i) ≈ 4×10⁻⁸ m/s² (mean projection) Significance, stated honestly: against the ToA-noise floor alone the formal single-pulsar significance would be ~7×10⁵σ, a meaningless number. Against the intrinsic-Ṗ cap the usable ratio is at most ~40–130× per pulsar (P = 3–10 ms), and Paper H's own analysis (verified verbatim: "The pulsar accelerations are nearly uninformative.", omega-centauri-mass-tension.html) shows the set is a weak mass discriminator even at that sensitivity, because the intrinsic terms are unknown and the bounds are one-sided. The ensemble radial-profile test, not any single-pulsar detection, is the program's real deliverable.

1.2 Model Discrimination Power

ObservableSingle IMBH (≥8,200 M☉)Dark cluster (~2–3×10⁵ M☉)
Radial acceleration profileSteep gradient, consistent with 1/r²Shallower, more extended profile
Central acceleration amplitudeHigh: ~10⁻⁷–10⁻⁶ m/s² for inner MSPsLower: ~10⁻⁸–10⁻⁷ m/s²
Higher-order derivatives (jerk)Detectable after 4+ years for inner MSPsDetectable but extended spatial distribution
Bayesian evidence ratioExpected B₁₂ ≫ 10 (decisive) after 5 yearsExpected B₁₂ ≫ 10 favoring extended model

1.3 Extended Science: Jerk and Spin Constraints

2. Observation Strategy

ParameterValue
CadenceBi-weekly (26 sessions/year); sessions ~5 hr each (cluster transit)
Pulsars targetedAll 18+ confirmed OC MSPs; priority on <0.5 pc projected separation
Integration per pulsar1 hr per session via sub-arraying: the ~20+ targets are observed in parallel sub-arrays within each 5-hr cluster transit. Assumption and cost (stated 2026-09-25, P0-15.6): splitting MeerKAT's 64 dishes into ~4 sub-arrays degrades each sub-array's SEFD ~4× relative to the full array, so per-ToA precision for a given pulsar is ~4× worse than a full-array hour at equal integration; the 26 × 1 hr/yr cadence per pulsar partially compensates through baseline, and the ToA figures quoted above (1–3 µs) are for these sub-array integrations
Frequency bandsL-band (1.28 GHz) + UHF (816 MHz); dual-band for DM corrections
New MSP discoveriesCommensal deep search expected to find 5–10 additional MSPs
SKA-MID transitionProgram continues on SKA-MID (~2030+) for 5× timing precision improvement
Total MeerKAT time (5 yr)~650 hours (26 sessions/yr × 5 hr × 5 yr; sub-arraying is already folded into the 5-hr/session figure above, not an additional reduction)

3. Work Plan

YearMilestoneDeliverable
1MeerKAT time allocation; timing pipeline; initial ephemeridesYear-1 timing baseline
2First acceleration measurements; preliminary potential mapConference results
3Radial acceleration profile; initial IMBH vs dark cluster Bayesian testPaper #1 (early constraints)
4Jerk measurements; strong model discriminationPaper #2 (model comparison)
5Full Bayesian reconstruction; SKA-MID handoff preparationPaper #3 (definitive constraints); public dataset

4. Budget

ItemCost (USD)
PI (0.2 FTE × 5 yr)400,000
Co-I × 2 (0.15 FTE × 5 yr)600,000
Postdoc × 2 (1.0 FTE × 5 yr)950,000
Graduate student (1.0 FTE × 5 yr)275,000
MeerKAT time + operational costs250,000
Computing (Bayesian inference, TEMPO2/enterprise)200,000
Travel + publications + fringe225,000
Total~$2.9M (5 yr)

5. References

  1. Bañares-Hernández, A., et al. (2025). New constraints on OC central mass via MSP timing. A&A, 693, A104.
  2. Häberle, M., et al. (2024). Fast-moving stars in ω Cen. Nature, 631, 285.
  3. González Prieto, A., et al. (2025). Growing the IMBH in ω Cen. ApJL, 990, L69.
  4. Chen, W., et al. (2023). MeerKAT discovery of 13 new pulsars in ω Cen, MNRAS, 520, 3847; 18 confirmed (19 after TRAPUM 2026).
  5. Prager, B. J., et al. (2017). Using pulsars to probe the potential well of NGC 6752. ApJ, 845, 148.
  6. Freire, P. C. C., et al. (2008). Timing MSPs in globular clusters as IMBH probes. ApJ, 675, 670.
  7. Babak, S., et al. (2017). LISA EMRI science. Phys. Rev. D, 95, 103012.
Working draft · April 2026 · The 5-year program requires sustained MeerKAT access beyond standard TAC allocations; a large program or partnership proposal with SARAO is recommended. ← Return to omegacentauri.me

Relevant tools

Pulsar Timing Array
MSP timing residuals as mass probe
Velocity Dispersion
Stellar kinematics mass estimator
IMBH Evidence Dashboard
Live multi-constraint overview
Spin-up Timeline
IMBH spin history via accretion