Neutron stars and Pulsar timing arrays as Axion giant gyroscopes
Yiming Liu
7520220161@bit.edu.cnSchool of Physics, Beijing Institute of Technology, Beijing, 100081, China
Jinneng Luo
3120221503@bit.edu.cnSchool of Physics, Beijing Institute of Technology, Beijing, 100081, China
Sichun Sun
sichunssun@gmail.comSchool of Physics, Beijing Institute of Technology, Beijing, 100081, China
Abstract
We consider the three-dimensional rotating motions of neutron stars blown by the “axion wind”. Neutron star precession and spin can change from the magnetic moment coupling to the oscillating axion background field, in analogy to the gyroscope motions with a driving force and the laboratory Nuclear Magnetic Resonance(NMR) detections of the axion. This effect modulates the pulse arrival time of the pulsar timing arrays. It shows up as a signal on the timing residual and two-point correlation function on the recent data of Nanograv and PPTA. The current measurement of PTAs can thus cast constraints on the axion-nucleon coupling as
††preprint: APS/123-QED
I Introduction
Oscillating background fields across the cosmic space, which behave like matter, are interesting dark matter candidates Hu et al. (2000); Hui et al. (2017); Arias et al. (2012); Feng (2010); Gibney (2020). The most common ones are either axionlike pseudoscalars or scalar fields such as moduli Preskill et al. (1983); Abbott and Sikivie (1983); Dine and Fischler (1983); Svrcek and Witten (2006). Since axion has a large parameter space, various axion detection schemes are proposed across the spectrum111https://github.com/cajohare/AxionLimits, from the quantum laboratory and collider studies to astrophysical observations. Inspired by the nucleon spin precession detection of axion particles such as CASPEr Garcon et al. (2017a), here we extend the detection principle using the microscopic nucleon spins to the macroscopic neutron stars, such that the three-dimensional rotating motion of spinning neutron stars blown by the axion wind can yield signals for a large range of parameter space for the axion mass and decay constant. Here in this paper we mainly discuss the case that ultralight axion wind blows on the pulsars in PTA, which yields the common spectrum on the residual with and two-point correlation function curve. We also briefly comment on the modulation of the pulsar period with different axion masses e.g.QCD axion , or the resonant frequency with pulsar period .
Recently several Pulsar Timing Arrays (PTAs) observatories including NANOGrav(North American Nanohertz Observatory for Gravitational Waves) Agazie et al. (2023a), Parkes PTA (PPTA) Reardon et al. (2023), European PTA (EPTA) Antoniadis et al. (2023a, b) and
Chinese PTA (CPTA) Xu et al. (2023),
report strong evidence for gravitational wave (GW) background at the nano-Hertz waveband,
which is consistent with previous claims Arzoumanian et al. (2020); Chen et al. (2021); Goncharov et al. (2021); Antoniadis et al. (2022).
Extensive studies on potential sources of the observed
GW background Cang et al. (2023); Franciolini et al. (2023a, b); Liu et al. (2023); Ellis et al. (2023); Wu et al. (2023a); Li et al. (2023a); Sun et al. (2022); Li et al. (2023b); Battista and De Falco (2021); De Falco and Battista (2023); Konoplya and Zhidenko (2023); Ben-Dayan et al. (2023); Balaji et al. (2023); Kohri and Terada (2021); Inomata et al. (2023a); Vagnozzi (2021); Benetti et al. (2022); Vagnozzi (2023); Guo et al. (2023); Oikonomou (2023a, b); Choudhury (2023); Choudhury et al. (2023); Bhattacharya et al. (2023); Madge et al. (2023); Sun and Zhang (2021),
including supermassive black holes Afzal et al. (2023); Middleton et al. (2021); Agazie et al. (2023b); Antoniadis et al. (2023c),
merging PBHs Depta et al. (2023); Gouttenoire et al. (2023),
phase transitions Bian et al. (2021); Arzoumanian et al. (2021); Xue et al. (2021); Wang (2022a) and
axion topological defects Wang (2022b); Ferreira et al. (2023); Inomata et al. (2023b), etc are conducted.
Nongravitational wave signals, such as the ultralight dark matter can also be observed through the timing residuals Khmelnitsky and Rubakov (2014); Porayko and Postnov (2014); Porayko et al. (2018); Sun et al. (2022); Wu et al. (2023b); Omiya et al. (2023); Nomura et al. (2020).
Most of the pulsars we observed are at the distances of order kpc from the Earth.
The de Broglie wavelength of the ultralight dark matter is around .
Fuzzy dark matter oscillates with the frequency
and the oscillation period .
The field behaves like pressureless cold dark matter on the cosmic scale. Notably, the oscillation frequencies of fuzzy dark matter models fell into the sensitive region of PTA observations.
II Axion background coupling to the pulsar magnetic moment
Taking into account the huge occupation number of the dark matter particle in the universe the axion field in the Galaxy can be viewed as a classical field and be approximated byDuan et al. (2023); Khmelnitsky and Rubakov (2014)
(1)
(2)
(3)
where is the mass of the axion, is the typical velocity of the axion and is the local dark matter density.
The Hamiltonian of the interaction between the axion and the nucleon is given byGraham and Rajendran (2013); Berlin et al. (2024)
(4)
where is the coulping constant and is the nuclear-spin operator.
As it is known that the magnetic moment can be expressed in terms of the nuclear spin
(5)
where is the gyromagnetic ratio of the nuclear spin, related to the property of the nucleon.
Define the effective magnetic field induced by the axionGarcon et al. (2017b):
(6)
therefore we can write Eq.(4) into the form below, similar to the Hamiltonian of the particle possessing the magnetic moment in an external magnetic field:
(7)
Now we generalize the above microscopic interaction to the macroscopic astrophysical object with a huge magnetic moment and the effective axion dark matter magnetic field.
Classically, the pulsar’s magnetic moment and spin angular momentum are given by:
(8)
(9)
where is the magnetic field of the pulsar’s pole, represents the pulsar’s mass, stands for the radius of the pulsar and is pulsar’s angular velocity.
Substitute Eq.(8), (9) into Eq.(5), and we can get the gyromagnetic ratio of the pulsar:
(10)
Then we can arrive at the energy change of the pulsar induced by the motion of the axion
(11)
where we have defined three unit vectors in the respective orientation.
Because the magnetic axis rotates around the spin axis with angular velocity , we can decompose as
(12)
Define:
(13)
(14)
(15)
We then arrive at:
(16)
Assume that all the energy is converted into the rotational energy of the pulsar
(17)
A slight change in the rotation period can be defined,
(18)
Dimensionless quantity then becomes
(19)
Next, we turn to the determination of the pulsar’s timing residuals. Assume that is the time pulsar takes to turn under the action of the axion field, is the time pulsar takes to turn the same angle without the action of the axion field, we then can determine the timing residual of the pulsar per unit timeMaggiore (2018)
(20)
According to the definition of the angular velocity, we can write , as
(21)
can be written as:
(22)
The timing residual of the pulsar is defined to be the integration of , measured with respect to a reference time :
(23)
Notice that when we do the integral , we do not include the case when which has to be computed separately, and the result is
(24)
The first term in (24) has a linear dependence in and the timing residual will tend to infinity as time goes on. The physical meaning of this divergence is that the axion and the pulsar happen to have a resonance around . and the angular velocity of the pulsar is changed permanently which can be searched accordingly. Nonetheless, given the wide spectrum of the axion and different types of pulsars, it seems plausible for the neutron star to resonate with the dark matter.
If the dark matter is the QCD axion with pulsars . Then we have . Then, the axion background will induce a nutation of the magnetic axis of pulsars, which might be observed by closely looking at the neutron star’s movement. We will leave this and the resonant case to a future study.
For the other case, to determine the value of , we choose:
(25)
(26)
(27)
as the initial direction of the three axes.
Then we have
(28)
(29)
(30)
Consider a special case when is parallel to , and is vertical to
(31)
For the ultralight dark matter ,
then we have
(32)
Consider another special case without losing generality:
(33)
The maximum value of the timing residual of the ultralight dark matter is
(34)
This is within the reach of current PTA observation. In Fig.1 we plot the timing residual induced by the axion together with the current Nanograv15 data. We can see that PTA data can already put constraints on the axion-nucleon coupling with some assumptions.
III Timing residual the two-point correlation
The cross-correlation Hellings-Downs curve between the noise-independent data is essential to detect the gravitational wave backgroundHellings and Downs (1983); Armaleo et al. (2021).
The cross-correlation between the timing residuals of the pulsar a and b can be calculated through , where represents a time average or an ensemble average. Having the explicit form of the redshift in hand, we can then calculate,
(35)
To simplify the calculation below, we then choose another frame that contains the pulsar location unit vector which points from the Earth towards the pulsar as the initial direction:
(36)
(37)
Without losing generality we set the initial direction of the magnetic axis of the pulsars a and b as:
(38)
(39)
We also set the orientation of the axion’s velocity and the rotational axes of pulsars a and b:
(40)
(41)
(42)
According to Eq.(15), we can easily calculate of pulsar a and b:
(43)
(44)
The pulsars we observe are distributed all over the sky, we can assume that their relative angle to the direction of the axion’s velocity is random. This means that the average of the cross-correlation signals among many pulsars effectively takes an average over the direction of the axion’s velocity. The rotational axis can also change its direction as a result of the motion of the galaxy, so we also take an average over the rotational axis and perform the integral over :
(45)
This curve is similar to the vector type modes in massive gravity types of studies Lee et al. (2008); Chamberlin and Siemens (2012); Cornish et al. (2018); Hellings and Downs (1983); Omiya et al. (2023) as expected, due to the one-dimensional directional nature of this modulation. We plot this axion-nucleon curve in Fig.2.
IV discussion
We consider a new detection possibility to employ the neutron stars as the axion wind detectors. Especially current PTA network observations can provide interesting sensitivity on the axion-nucleon couplings. This classical macroscopic detection mechanism is complementary to axion microscopic production and capture mechanism Wu et al. (2023c); Beznogov et al. (2018); Dev et al. (2024); Berlin et al. (2024); Chao et al. (2023); Gao et al. (2024); Caputo et al. (2020) to set the limits on the related coupling. Due to the complexity of the neutron star structure models, we expect noises from several different sources to get into the neutron star motion. Nevertheless, a close analysis of neutron star behavior might provide further insight into this type of coupling.
Acknowledgements
We thank Fei Gao for the discussion on the neutron stars. This work is supported by the National Key RD Program of China (grant 2023YFE0117200) and the National Natural Science Foundation of China (Nos. 12105013).
Bhattacharya et al. (2023)G. Bhattacharya, S. Choudhury, K. Dey,
S. Ghosh, A. Karde, and N. S. Mishra, (2023), arXiv:2309.00973
[astro-ph.CO] .
Madge et al. (2023)E. Madge, E. Morgante,
C. Puchades-Ibáñez,
N. Ramberg, W. Ratzinger, S. Schenk, and P. Schwaller, (2023), arXiv:2306.14856
[hep-ph] .
Garcon et al. (2017b)A. Garcon, D. Aybas,
J. W. Blanchard, G. Centers, N. L. Figueroa, P. W. Graham, D. F. J. Kimball, S. Rajendran, M. G. Sendra, A. O. Sushkov, L. Trahms, T. Wang, A. Wickenbrock, T. Wu, and D. Budker, Quantum Science and Technology 3, 014008 (2017b).
Maggiore (2018)M. Maggiore, Gravitational Waves.
Vol. 2: Astrophysics and Cosmology (Oxford
University Press, 2018).