X-ray emission spectrum for axion-photon conversion in magnetospheres of strongly magnetized neutron stars

Detecting axionic dark matter (DM) could be possible in an X-ray spectrum from strongly magnetized neutron stars (NSs). We examine the possibility of axion-photon conversion in the magnetospheres of strongly magnetized NSs. In the current work, we investigate how the modified Tolman Oppenheimer Volkoff (TOV) system of equations (in the presence of a magnetic field) affects the energy spectrum of axions and axions-converted-photon flux. We have considered the distance-dependent magnetic field in the modified TOV system of equations. We employ three different equations of states (EoSs), namely APR, FPS, and SLY, to solve these equations. We obtain the axions emission rate by including the Cooper-pair-breaking formation process and Bremsstrahlung process in the core of NSs using the NSCool code. We primarily focus on Magnificient seven (M7) star RXJ 1856.5-3754. We further investigate the impact of the magnetic field on the actual observables, such as axion energy spectrum and axion-converted-photon flux at an axion mass in meV range by assuming mass $M_{NS} \sim 1.4M_{\odot}$. We compare our calculated axion-converted-photon flux from all available archival data sets from PN+MOS+Chandra. We also study the variation of the energy spectrum at a fixed energy with varying central magnetic fields. Our predicted axion-converted-photon flux values as a function of axion energy closely follow the experimentally archival data, which allows us to put bounds on the axion mass for the three EoS.

The QCD axion and axion-like particles (ALPs) are anticipated to pair derivatively to fermionic matter and electromagnetism, thus enabling their production inside the heated cores of NSs.The effective Lagrangian of axions is given by: The second term in the expression denotes axion coupling to a fermionic field ψ f [17].According to the recent experimental results of the axionic DM experiment (ADMX) [72,[72][73][74] and Large-scale microwave cavity search [75], axions can solve the problem of mysterious DM.Currently, there are two possible axion models, namely the Dean-Fischler-Srednitsky-Zhitnitsky -DFSZ model, which contains an additional coupling of axions with the charged leptons, and the Kim-Shifman-Weinstein-Zakharov (hadronic) -KSVZ model which is related to interaction of axions with photons and hadrons [76,77].The production rate of axions depends on their emissivity, critical temperature profiles of the NSs, and proton and neutron Fermi momenta profiles.There is an extensive literature on the axion mass bounds: plasma haloscopes and MADMAX [78,79] reveals an axion mass between m a ∼ 40-400 µeV, ex-perimental data from HAYSTAC and ADMX reports an axion mass m a ∼ 1-100 µeV, and CASPEr [80][81][82], DM-Radio [83] and ABRACADABRA [60,84,85], confirms that mass of axion cannot be in µeV range (m a ≪ µeV).The experimental data of ADMX (in context to Dean-Fischler-Srednitsky-Zhitnitsky (DFSZ)) model presents an axion mass range m a ∼ 2.7 − 4.2 µeV [53,[72][73][74]77].Nevertheless, the black hole superradiance [86,87] reports mass of QCD axion as m a > 2 × 10 −11 eV.Buschmann et al. [6] (about Kim-Shifman-Weinstein-Zakharov (KSVZ) axion model), restricts the mass of axion at m a 16 meV.For the KSVZ model [50,51] the experimental results from the HAYSTAC [67] and ADMX [75] indicate that the range of axion mass is still unconfirmed within the parameter range of about ten powers of magnitude.Exploring the mass of axionic DM (that is, not in the meV range) in a laboratory is a bit challenging [88][89][90][91][92][93].Recent study [94][95][96], reports the axion mass bound as m a 10 meV.
Various studies describe the potential use of resonant conversion of axions, particularly in the µeV-meV mass range, relevant to QCD (Quantum Chromodynamics) axionic DM.Resonant conversion refers to axion-photon conversion under plasma densities and intense magnetic fields in the magnetosphere of strongly magnetized NSs.The axion mass is crucial in determining the frequency of the associated photon signal.These features have drawn researchers' attention to conduct studies on axionphotons conversion in the NSs magnetosphere [17,97,98].They have been actively characterizing and searching for potential X-ray and radio signals.This area of research holds promise for advancing our understanding of DM and theories related to fundamental particle physics, provided that the predicted signals can be detected and verified through observational data.It has also been revealed in the literature [18,20,21,29,36,[99][100][101] that the distribution of a strong magnetic field affects both internal structure as well as the cooling rates of NSs.Witte et al. [102] performed a calculation to compute the radio signal from the NS magnetosphere by adopting a charge-symmetric Goldreich-Julian model.They accomplish this by employing a cutting-edge autodifferentiation ray tracing technique, which enables them to closely monitor photon absorption, energy exchanges with the surrounding plasma, axion-photon dephasing, and the properties of reflections and refraction. A. Hook et al. [97] reported that axions can be detected through the narrow radio lines from axionic DM conversion in the NS magnetosphere.They solved the axion-photon mixing equations by including magnetized plasma.The conversion can be possible when plasma frequency matches the mass of axions.McDonald et al. [103] derived the essential transport equations in magnetized plasmas.Due to the resonant conversion of axions to photons, these equations describe the creation and propagation of photons in any 3-D medium.Millar et al. [104] shows the axion-photon conversion in a highly magne-tized anisotropic material in 3-D compared to the simplistic 1-D calculation mentioned in the literature.The results show orders of magnitude difference in conversion for the 3-D medium compared to 1-D.The findings will significantly impact the radio signal seen through a telescope.In the current work, we investigate how the modified Tolman Oppenheimer Volkoff (TOV) system of equations (in the presence of magnetic field) [105,106] affects the axions energy spectrum and axion-convertedphoton flux at an axion mass in meV range by assuming mass M N S ∼ 1.4M ⊙ of the non-rotating NSs.Initially, we solved the modified TOV system of equations by employing three EoSs, namely APR, FPS, and SLY, and the distance-dependent magnetic field to build the profiles.We change the NSCool code by incorporating the N-N Bremsstrahlung process and the Cooper-pair-breaking formation process to determine the axion emission rate inside the core of NSs.In particular, our fiducial analysis includes Magnificient seven star, RXJ 1856.5-3754(M7) with the characteristic age of 3.7 × 10 6 yrs [107,108].
Furthermore, we investigate the influence of the strong magnetic field on the fundamental observables, such as the axion-converted-photon flux and energy spectrum of axions.We also compare our computed axionconverted-photon flux with archival available data from PN+MOS+Chandra [17,109] by taking axion-photon coupling in the magnetosphere of NSs.We also obtained the axion mass bounds from the three mentioned EoSs and compared them with the existing bounds in the literature.The paper is organized as follows.In Section I, we begin with a brief Introduction.In Section II, we describe the modified TOV system of equations, Equation of states, Neutron star cooling, and the axion emission rates by Cooper-Breaking and formation and Bremsstrahlung process in the core of NSs, Energy Spectrum of Axions, Conversion Probability and Axionconverted-Photon Flux.In Section III, we explore the impact of magnetic fields on various observables, such as the energy spectrum of axions and axion-convertedphoton flux for different NSs at axion mass in the meV range for three different EoSs.Finally, Section IV summarizes and concludes the work.

II. FORMALISM A. Modified TOV equations for non-rotating neutron star
The present work relies on the modified TOV system of equations (in the presence of a strong magnetic field) [110].As discussed in [22,110,111], this would require us to solve the Einstein's field equations in the stellar interior including the full electromagnetic Energy-Momentum tensor [112].Perturbative solutions to study the neutron star's deformation developed in [113] shall also be valid only for low magnetic fields.In this work, we adopt a more phenomenological approach to replicate the effects of the magnetic field (at least qualitatively).This is done by including a magnetic energy density to the usual fluid energy density in the TOV equations.Further, the conservation of the Energy Momentum tensor in the general relativistic framework [112] motivates us to add a Lorentz force term in the pressure gradient equation.This Lorentz force term, which imprints the effect of the magnetic field is also motivated by the anisotropic term in the pressure gradient equation studied in [114].To examine the effects of magnetic fields, we require the suitable EoS [22,23,115,116], mass, pressure, and baryon density as a function of distance (profiles) and the magnetic field profile.The TOV system of equations for mass, gravitational potential, and hydrostatic equilibrium, respectively, are given by: where Λ(r) corresponds to a Lorentz force term [112].Instead of solving the full Maxwell-Einstein system of equations, we have adopted the phenomenological fit function for Λ(r) from [110].This parametrized fit function depends on the star's mean radius and central magnetic field B c , given as: (5) The detailed formalism is presented in our recent work [96].Various magnetic field profiles are available in the literature to describe the state of super-dense matter of NSs in strong magnetic fields.Chatterjee et al. [110] explored the impact of the magnetic field using an axisymmetric numerical code by changing three parameters: mass, EoSs, and the strength of the magnetic field for non-rotating NSs.According to the azimuthal symmetry, the magnitude of the magnetic field B = ( B. B) 1/2 is given by the expression: As per the numerical results, the magnetic monopole term dominates over all other terms.All the anisotropies (first approximation) are beyond our current work [96].The simplified radial profile for the monopole term of B is given as: here r is the mean radius, and B c corresponds to the value of a central magnetic field for the NSs.Finding a mean radius is difficult in the presence of a magnetic field.So, r value is taken more than the actual radius of the NSs.
We have assumed the mean radius to be r = 13 km and the central magnetic field B c = 10 18 Gauss.Further, it should be noted that in this work, the possible impact of the magnetic field on the composition of NS matter (populations of neutrons, protons, and electrons) has not been considered.

B. Equation of State
We have employed three EoSs, namely APR [117][118][119], FPS [120], and SLY [121] EoSs to study the thermodynamic behavior of the interiors of NSs.The Akmal-Pandharipande-Ravenhall (APR) EoSs describe the interaction potential by isospin asymmetry and baryon density.This EoS depicts the transition from the lowerdensity phase to the higher-density phase.In principle, it has also been observed that phase transition speeds up the shrink rate of NSs.The FPS EoS shows a unified description of the liquid core and inner crust of the NSs.In FPS EoSs, the liquid-crust core transition usually takes place at density ρ edge = 1.6 × 10 14 gm/cm 3 .All the existing phase transitions are first-order, with relative density jumps less than 1%.Considering different nuclear Hamiltonian, FPS EoS is more suitable than any other EoS.The SLY EoS is based on the effective nuclear interaction model in which the nuclei in the ground state remain spherically downward towards the bottom edge of the inner crust of NSs.In this EoS, the transition towards the ionized matter plasma happens at the density ρ = 1.3 × 10 14 gm/cm 3 .The above EoSs in detail are described in our previous work [96].The maximum mass for the three EoS is obtained using [22,36] (corresponding to the ratio of the magnetic and fluid pressure close to unity) is shown in Table I.We have also tabulated the radius of the compact object for the three EoS for a mass of 1.4 M ⊙ .

C. Neutron star cooling
To analyze the state of highly dense matter and the emission properties of the strongly magnetized NSs, we have modified the publicly available NSCool code [122].A Fortran-based cooling computational code that solves the equations related to energy balance and heat transport.This package also contains various pre-built stars, different EoSs, and a TOV integrator to form the stars.Considering the spherically symmetric NSs with interiors as isothermal and solving its energy balance equation, one can reasonably understand the various cooling simulations of such compact objects.The equation for the conservation of energy of a star as per the Newtonian formulation is given by [123]: Here C v specific heat, L ∞ ν is the energy sinks are the total neutrino luminosity, L ∞ a is the luminosity of axions, L ∞ γ is luminosity of photons, T s is the surface temperature and T ∞ b internal temperature of the NS.
, where σ is Stefan-Boltzmann's constant and R is the radius of the NS.T ∞ s = T s 1 − 2GM/c 2 R, here infinity superscript means that the external observer is located at the infinity and measures these terms on the redshifted scale.We have not included any heating mechanisms, H = 0, in the above equation.Specifically, for non-rotating NSs of mass M ∼ 1.4M ⊙ , T ∞ s /T s ∼ 0.7.The NS cooling mechanisms [124,125] is the crucial ingredient to study the hypothesis related to the axionic DM.The structure of heat blanketing envelopes (HBEs) may be affected by neutrino cooling mechanisms [3,26,27,39,126,127].The boundary of the HBEs is ǫ b (ρ b ∼ 10 10 gm/cm 3 ) and magnetic field at HBEs is taken as 10 13 Gauss.The HBEs comprise iron or successive iron, carbon, hydrogen, and helium layers.The formula that determines the radiative opacity for a completely ionized, non-degenerate plasma is made up of electrons and ions and gets modified in the presence of a magnetic field [3,128].Also, the heat blanketing layer of NSs will strongly affect the electron thermal conductivity.In the presence of a magnetic field, the thermal conductivity of electrons considered for current calculations is given in [3,124].A detailed explanation is provided in our previous work [96].

D. Cooper pair Breaking and Formation process and N-N Bremsstrahlung process
The composition of the inner core of NSs is somewhat uncertain, and it depends on the considered EoSs, which is still an open area of research [129].The concept of pairing has opened up new areas for the dominant neutrino emissions.The photon and neutrino emissions are controlled by the stellar outer layers' structure and the state of superdense nuclear matter, respectively.There is a continuous formation and breaking of n − n and p − p cooper pairs, leading to the formation of the superfluid material.This happens lower than the superfluid critical temperature T c .During the breaking and formation process, axions may be emitted.The energy liberated in this process is carried out by a ν ν pairs [130].
Here, X can be protons or neutrons.

Axion emission rate by Cooper pair Breaking and Formation process (PBF)
This section briefly describes the axion emission rate from the Cooper pair breaking and formation (PBF) process [130,131] inside the NS core.To calculate the axion production rates in the cores, the prerequisites are the temperature profiles (from modified TOV equations), metric, critical temperature profiles, and neutron and proton Fermi's momenta profile (which depends on the EoSs).The axion emission rate due to the neutron Swave pairing is given by [130]: The integral I s a is given by: where, ǫ s a is emissivity of axions, f a is axion decay constant, v n (0) density of states at fermi surface and v F (n) is the fermi velocity of the neutron.
with x = ω 2T , where ω is the axion energy.
The ratio between axion and neutrino emissivity is given by: where, The numerical values of the r(z) = z 2 I s a /I s ν are related to integrals corresponding to axion and neutrino emissivity for the different values of z.
Here f a is an axion model-dependent variable.The axion mass m a is related to f a given by the equation: Here, the axion emissivity due to proton S-wave pairing is the same as neutron S-wave pairing provided C n = C p and T cn = T cp , where T cn and T cp are critical temperatures for neutron superfluid and proton superfluid, respectively.We have assumed T cp = T cn = 10 9 K.The axion emissivity due to the P-wave paired neutron superfluid is given by [130,132]: where C n is a constant that depends on the axion model.The integral I p a is expressed as: where z n = ∆ P (T, θ)/T where θ is the angle between the quantization axis and momentum of the neutron.There exist two states (A and B) of P-wave superfluid pairing corresponding to m J = 0 and 2 (representing the projection of the angular momentum of the Cooper-pair on the z-axis).The superfluid gaps are given by [17,130]

Axion emission rate by N-N Bremsstrahlung process
The axions, along with neutrinos, are expected to emit during the nucleon collisions by the bremsstrahlung process [130,[133][134][135][136]. To calculate the production rates in the core, the prerequisites are the processes involving protons, neutrons, and electrons, the purely nucleonic EOS at all possible densities.In the current analysis, the axion production occurs in the fully degenerate nucleonnucleon bremsstrahlung process, N + N → N + N + a, where N can be either a neutron or a proton, and a is the emitted axion.The current process involves the emission of radiation (in this case, axions) during the interaction of charged particles (neutrons and protons in this context) [133,136].The expression for axion emission rate is given by: The axion integral is expressed as: where, Γ σ means the spin rate (nucleon) changes under collisions with the other nucleon.For the degenerate limit: Here p F is the fermi's momentum of the nucleon.We note that we have not included the effect of the magnetic field on the mentioned cooling mechanisms.

E. Energy Spectrum of Axions
Cooper Breaking and Formation process (PBF) As mentioned in the previous work [17], the local energy spectrum of axions emitted from the core by the Cooper Breaking and Formation process follows the modified thermal distribution given as: where, E is local axion energy and F is the axionconverted-photon flux.Both spin−0 S-wave and spin−1 P-wave nucleon superfluids could be possible inside the core of NSs.The energy spectrum of axions due to the spin−0 S-wave process is given as: where, is the normalization constant derived from the expression ∞ 2∆T J s,P BF a dω a = ǫ s,P BF a and 2y∆T is an axion energy.T corresponds to a locally measurable quantity at some radius within the interiors of NSs.The energy spectrum of axions due to the P-wave pairing is given by [17]: is the normalisation constant which is defined as ∞ 2∆ P (T,θ) J p,P BF a dω a = ǫ p,P BF a and 2y∆ P (T, θ) is the axion energy ω a .

N-N Bremsstrahlung process
The energy spectrum of axions due to the N-N bremsstrahlung process is given by: T is the dimensionless quantity and ω a (equivalently x × T ) is the axion energy.

F. Conversion Probability and Axion-converted-Photon Flux
Axions produced in the cores escape the stars because of their weak interaction with the ordinary matter and further get converted to the photons in an extended magnetosphere of strongly magnetized NSs [98,137].The axions mixing with photons is based on the axion-photon coupling term given by the Lagrangian: Here, F & F corresponds to an electromagnetic field strength tensor and corresponding dual tensor, respectively, g aγγ is axion-photon coupling parameter given as [17]: where, C γ is a axion model dependent factor and α = 1/137 is QED fine structure constant.In the magnetosphere of strongly magnetized NSs, the axion-photon mixing term may rotate the initial state of an axion to an electromagnetic wave that is polarised exactly in the same direction as the external magnetic fields.The following expression provides an approximate relationship for conversion probability [17,98,137,138]: where, ω is the axion energy and R N S is radius of NS (in km).B 0 is the strength of the surface magnetic field at the pole, and θ is the magnetic axis polar angle.Finally, we have multiplied the axion-photon conversion probability with the axion energy spectrum to obtain the final axion-converted-photon flux at the different values of axion energy.

III. RESULTS AND DISCUSSIONS
Figure (1a) shows the energy spectrum of axions as a function of axion energy for APR EoS.The local energy spectrum of axions emitted by the Bremsstrahlung plus PBF process remains dominant over those emitted only by the Bremsstrahlung process for the entire range of the axion energies.As the axion energies increase, for a no-magnetic field case, the energy spectrum of axions remains higher than when there is a magnetic field at all values of axion energies.In the Bremsstrahlung plus PBF process, the separation in the energy spectrum of axions is almost constant up to 10 keV axion energy, but after that, separation decreases, and both the curves overlap around 10 2 keV energy.However, if we consider only the Bremsstrahlung process, a significant difference can be observed in the energy spectrum of axions in both the presence and absence of a magnetic field.The magnitude of the axion energy spectrum and axion-convertedphoton flux is quite sensitive to the axion mass, and a significant departure is observed in the presence of the magnetic field.In figure (1b), we have presented the axion-converted-photon flux as a function of axion energies for APR EoS.The axion-converted-photon flux from the Bremsstrahlung plus PBF process dominates from the only Bremsstrahlung process for all ranges of axion energies.Including a magnetic field does not change the qualitative behavior, although a significant difference is observed.As we focus on the Bremsstrahlung plus PBF process, a significant difference in the axion energies spectrum can be observed upto axion energy of ∼ 10 keV.After that, the separation decreases and the curves overlap beyond ∼ 70 keV.We also show the data points from X-ray spectra at an axion mass 16 meV for RXJ 1856.5-3754(M7) star by combining all the available archival data sets from PN+MOS+Chandra [17,109].The axion mass bound is around 16 meV and is obtained by fitting the various axion mass spectral models corresponding to available data sets.On comparing with the observational data, our predicted results show a better agreement with the magnetic field case.The figure (2a) shows the variation of the energy spectrum of axions versus axion ener-   gies for FPS EoS.Here, the local energy spectrum of axions emitted by the Bremsstrahlung plus PBF processes dominates over the only Bremsstrahlung process for the entire range of axions energies.The separation between (with/without) magnetic fields in the energy spectrum of axions is constant up to ∼ 10 keV axion energies for the Bremsstrahlung plus PBF processes.If we include both processes, we find that in the presence of a magnetic field, the energy spectrum of axions in the low axion energies region is lesser than without a magnetic field.The gap becomes negligible at higher values of axion energies, ω a > 95 keV.However, for only the Bremsstrahlung process, a constant difference can be seen in the energy spectrum of axions in the presence/absence of magnetic field for axion energy ∼ 4 keV.Beyond this energy, this difference decreases until both curves overlap between 10 − 20 keV axion energies.The figure (2b) depicts the variation of axion-converted-photon flux as a function of axion energies for FPS EoS.After including the magnetic field, the axion-converted-photon flux when both the Bremsstrahlung and PBF processes are included is more significant than the Bremsstrahlung process.This is true for all values of the axion energies.The separation (in the presence/absence of a magnetic field) for the axion-converted-photon flux emitted by the Bremsstrahlung plus PBF process is almost similar as compared to other mentioned EoSs.The separation decreases from ∼ 30 keV axion energy, and finally, both the curves (with/without magnetic field) overlap beyond ∼ 90 keV axion energies.If we consider axionconverted-photon flux from the Bremsstrahlung process only.In this case, the separation remains constant up to ∼ 3 keV, and after that, curves (with/without magnetic field) overlap beyond the axion energy of ∼ 10 keV.We also show the data points from X-ray spectra for RXJ 1856.5-3754(M7) star from all available archival data sets from PN+MOS+Chandra [17,109].Again, by fitting the several axion mass spectral models corresponding to the available data sets, we obtain the axion mass bound of 10 meV in the case of FPS EoS.Compared to the available experimental data, our calculated results demonstrate better agreement with the magnetic field.For FPS EoS, the magnitude of the axion's energy spectrum and axion-converted-photon flux is quite sensitive to the mass of an axion.One can also observe a significant departure in the presence of the magnetic field.In the figure (3a), we have shown a variation of the energy spectrum of axions versus axion energies for SLY EoS.In the presence of a magnetic field, at all values of the axion energies, the local energy spectrum of axions emitted by the Bremsstrahlung plus PBF process dominates over the case when only the Bremsstrahlung process is included.This behavior is similar to other mentioned EOSs.Considering the Bremsstrahlung plus PBF processes, the separation (with/without magnetic field) is almost similar to the other EoSs but for lesser values of the axion's energy only.One can observe a significant difference in the energy spectrum for the entire range of axion energy.The imprint of the strong magnetic field on the energy spectrum of axions from the Bremsstrahlung plus PBF processes can be seen in this EoS but overall qualitative features remain the same.For the Bremsstrahlung process, a significant departure in the energy spectrum of axions is observed at lower energies before the spectrum declines quickly at the higher values of axion energies.In the figure (3b), we have presented the variation of axion-converted-photon flux as a function of axion energies for SLY EoS.For the Bremsstrahlung plus PBF process, in the presence of a magnetic field, the axionconverted-photon flux consistently deviates as compared to the case of a no-magnetic field.Still, it does not behave qualitatively differently from the other mentioned EoS.Additionally, the axion-converted-photon flux is always higher in the no-magnetic field scenario than in the case of the magnetic field, regardless of all axion energies.After considering both processes, a significant difference can be seen for the entire range of axion energies (∼ 1 − 100 keV).Specifically, the axion-convertedphoton flux is steeply peaked at low axion energies for the Bremsstrahlung process.This rapidly decreases at higher values of axion energies because of energy conservation.Here again, we have presented the X-ray spectra of the RXJ 1856.5-3754(M7) star using the data sets from PN+MOS+Chandra by employing the SLY EoS [17,109].After fitting the various axion mass spectral models with the available data sets, we report an axion mass bound around 7 meV for the SLY EoS.Our anticipated results, on comparsion with the observational data, show better agreements with the magnetic field case than the no-magnetic field situation.
The strong magnetic fields and the corresponding high value of magnetic field energy lead to a drop in the star's thermal energy while the contributions from the overall gravitational energy stay unchanged.This explains the correlation between the luminosity suppression and the magnetic field strength [42].This work has not considered the effect of magnetic fields on the mentioned EoS.It is anticipated that the axion mass bounds would exhibit degeneracies with the mentioned EoS.However, the current analysis suggests that the EoS influences the macroscopic properties (mass and pressure profiles) by modified TOV equations, thereby influencing the bounds on axion mass.
In the Bremsstrahlung plus PBF process, the energy spectrum of axions for the APR, FPS, and SLY EoSs show an ∼ 18.9%, 18.1%, 16.5% departure from without the magnetic field at an axion energy ∼ 5 keV.Similarly, for the Bremsstrahlung process, the energy spectrum of axions for the APR, FPS, and SLY EoSs show an ∼ 6.8%, 4.7%, 5.9% departure from without the magnetic field at an axion energy ∼ 5 keV.
In addition to this, Table II summarises the percentage changes of axion-converted-photon flux for APR, FPS, and SLY EoS at an axion energy of 5 keV in the presence and absence of a magnetic field for RXJ 1856.5-3754(M7) NSs.The effect of the magnetic field can be visualized by studying the variation of energy spectrum at a fixed energy with the central magnetic field.Figure (4a), (4b) and (4c) shows the variations of the energy spectrum of axions at axion energies ω a = 5 keV and ω a = 40 keV with varying central magnetic field for three EoSs.The rate of decrease in the energy spectrum is slightly lesser at a higher axion energy of 40 keV compared to the axion energy of 5 keV from B = 0 G to B = 10 18 G.However, the qualitative behavior remains the same for all three EoSs.

IV. SUMMARY AND CONCLUSIONS
In the present work, we study how the presence of a strong magnetic field affects the energy spectrum of axions and axion-converted-photon flux.Our current work is based on the hypothesis that the axions are produced in the core of NSs by the Cooper pair breaking and formation process and the Bremsstrahlung process, then converted into X-rays due to a strong magnetic field in NS magnetospheres.We have calculated the axion mass bound for the APR, FPS, and SLY EoS.We have also explored the effect of the magnetic field by comparing our numerically computed results with the available archival data of PN+MOS+Chandra [17,109] for Magnificient seven (M7) star RXJ 1856.5-3754(M7).Additionally, the summary of our findings is as follows: • We have not included the influence of magnetic field on the employed EoSs and mechanism related to the production of axions.Initially, the modified Tolman-Oppenheimer-Volkoff (TOV) system of equations was solved for a given EoSs to study the macroscopic properties (mass/pressure) profiles, which are then used to compute luminosities of axions and surface temperature as a function of time using the NSCool code.
• Our work encompassed the energy spectrum of axions and axion-converted-photon flux versus axion energies in the presence and absence of magnetic fields.We find the average best-fit axion mass bound for the mentioned three EOSs; namely APR, FPS, and SLY, as 16 meV, 10 meV, and 7 meV, respectively.These mass bounds lie within the limits of the recently reported axion mass bound of 10 meV and 16 meV in the case of a no-magnetic field.
The results have been calculated with a fixed value of the star's mean radius r = 13 km and a high value of the central magnetic field 10 18 Gauss.
• We emphasize that the axion mass bound slightly varies for APR, FPS, and SLY EoSs.On comparing our results with all available archival data sets of PN+MOS+Chandra, we find that the results for axion-converted-photon flux closely matches with the magnetic field, especially for APR and FPS EOSs.
• Our investigation reports that the influence of a strong magnetic field on the observables depends on the employed EoS.The uncertainties are more pronounced, especially with low temperatures and high-density conditions.Thus, the study of axionphoton conversion remains a subject of active research for these highly magnetized NSs. (2007).

FIG. 1 :
FIG.1:The variation of the energy spectrum of axions (left) with axion energies by Bremsstrahlung plus PBF and only Bremsstrahlung process for RXJ 1856.5-3754(M7) NS at an axion mass 16 meV for APR EoS in the presence of magnetic field.

FIG. 2 :
FIG. 2:The variation of the energy spectrum of axions (left) and axion-converted-photon flux (right) with axion energies by Bremsstrahlung plus PBF and only Bremsstrahlung process for RXJ 1856.5-3754(M7) NS at an axion mass 10 meV for FPS EoS in the presence and absence of magnetic field.
FIG. 3:The variation of the energy spectrum of axions (left) and axion-converted-photon flux (right) with axion energies by Bremsstrahlung plus PBF and only Bremsstrahlung process for RXJ 1856.5-3754(M7) NS at an axion mass 7 meV for SLY EoS in the presence and absence of magnetic field.

FIG. 4 :
FIG. 4: The variation of the energy spectrum of axions with different values of central magnetic field Bc at axion energies ωa = 5 keV and ωa = 40 keV from PBF plus Bremsstrahlung process for RXJ 1856.5-3754(M7) NS and for three EoS.The results for energy spectrum of axions corresponding to axion energy ωa = 40 keV are multiplied by a factor of 100 to compare with the results of ωa = 5 keV.
Figure (5a) and Figure (5b) show the variation of the energy spectrum of axions and axion-converted-photon flux as a function of axion energies from Bremsstrahlung plus PBF and only Bremsstrahlung process for star PSR B0531+21 at an axion mass 16 meV for APR EoS, respectively.In Figure (6a) and Figure (6b), we have presented the variation of the energy spectrum of axions and axionconverted-photon flux as a function of axion energies from Bremsstrahlung plus PBF and only Bremsstrahlung process for star PSR J0538+2817 at an axion mass 16 meV for APR EoS, respectively.

FIG. 5 :
FIG.5:The variation of the energy spectrum of axions (left) and the axion-converted-photon flux (right) with the axion energies by Bremsstrahlung plus PBF and only Bremsstrahlung process at an axion mass 16 meV of star PSR B0531+21 for APR EoS with and without including magnetic field.

FIG. 6 :
FIG.6:The variation of the energy spectrum of axions (left) and the axion-converted-photon flux (right) with the axion energies by Bremsstrahlung plus PBF and only Bremsstrahlung process at an axion mass 16 meV of star PSR J0538+2817 for APR EoS in the presence and absence of magnetic field.

FIG. 7 :
FIG. 7:The variation of the energy spectrum of axions (left) and the axion-converted-photon flux (right) with the axion energies by Bremsstrahlung plus PBF and only Bremsstrahlung process at an axion mass 16 meV of star PSR J0531+21 for SLY EoS in the presence and absence of magnetic field.

FIG. 8 :
FIG.8:The variation of the energy spectrum of axions (left) and the axion-converted-photon flux (right) with the axion energies by Bremsstrahlung plus PBF and only Bremsstrahlung process at an axion mass 16 meV of star PSR J0538+2817 for SLY EoS with and without including the magnetic field.

TABLE I :
Table showing the maximum mass (M⊙) and radius (km) corresponding to the mass M = 1.4M⊙ of NSs for three EoS.

TABLE II :
Table shows the percentage changes of axionconverted-photon flux by Bremsstrahlung plus PBF and only Bremsstrahlung process of RXJ 1856.5-3754(M7) NSs for three EoS at an axion energy of 5 keV.