A Search for a Contribution from Axion-Like Particles to the X-Ray Diffuse Background Utilizing the Earth's Magnetic Field

The Axion Like Particle (ALP) is a hypothetical pseudo-scalar particle beyond the Standard Model, with a compelling possible connection to dark matter and early universe physics. ALPs can be converted into photons via interactions with magnetic fields in the universe, i.e., the so-called inverse Primakoff effect. In this paper, we propose a novel method to explore ALP-induced photons from X-ray data obtained from the {\it Suzaku} satellite, arising from a possible interaction of ALPs with the direction-dependent Earth's magnetic field viewed from the satellite. {\it Suzaku} data is suitable for this purpose because its low-altitude Earth orbit result in intrinsically low cosmic-ray background radiation. We study whether the X-ray diffuse background (XDB) spectra estimated from the four deep fields collected over eight years, vary with the integrated Earth's magnetic strength in the direction of each target field at each observation epoch, which amounts to $10^2$ Tm-a value greater than that achieved by terrestrial experiments due to the large coherent length. From the detailed analysis, we did not find evidence of the XDB confidence level spectra having dependence on the Earth's magnetic strength. We obtained 99 % confidence level upper limit on a possible residual contribution to the cosmic X-ray background (CXB) surface brightness to be $1.6\times 10^{-9}~{\rm ergs~s}^{-1}{\rm cm}^{-2}{\rm sr}^{-1}$ normalized at $10^4$ T${}^2$ m${}^2$ in the 2-6 keV range, which corresponds to 6-15 % of the observed CXB brightness, depending on which model of unresolved point sources are used in the interpretation. It is consistent with 80-90 % of the CXB now being resolved into point sources.


Introduction
Various cosmological observations have provided strong evidence for the existence of dark matter. However, if dark matter is to be an elementary particle, it is a yet-unknown particle beyond the Standard Model. The Axion-Like Particle (ALP) is a hypothetical pseudo-scalar particle beyond the Standard Model, and is a consequence of the quantum field for conserving CP symmetry in the strong interaction [1,2]. ALPs are attractive because they act like cold dark matter (CDM) in the formation of cosmic structure. It is possible that ALPs are created by a decay of other CDM-candidate particles. If ALP mass is too low, a direct experiment with present techniques is unlikely to find it. A possible channel is to observe photons, which are created by ALPs following the inverse Primakoff process in an electromagnetic field. As we show in detail in Section 2, the ALP-photon conversion probability P a→γ is approximately proportional to the squared product of magnetic field strength orthogonal to the ALP momentum direction, B ⊥ , and the length, L, i.e., P a→γ ∝ (B ⊥ L) 2 .
There have been many attempts to detect ALP signals in terrestrial experiments and astronomical observations. One candidate signal is in the direction of galaxies or galaxy clusters, which was proposed to be due to ALP interaction with inter-stellar or galactic magnetic fields [3][4][5], although the results are still under discussion. If ALPs are CDM itself or produced from a decay of CDM at cosmological distances and if those can be observed, the distribution of ALPs or ALP-induced photons via its interaction with magnetic fields in cosmic structures should appear to be isotropic in the sky to the zero-th order approximation, unless we have high-sensitivity and high-angular resolution data to resolve the distribution tracing inhomogeneous cosmic structures in the universe.
In this paper, we propose a novel method to search for ALP-induced photons from the satellite X-ray data, arising from the Primakoff interaction of the ALPs with the Earth's magnetic field. The Earth's magnetic field is known to have a dipole structure around the Earth in a north-south direction and we have a good knowledge of its strength and field configuration from various observations. Therefore, we can expect the ALP-induced photons in X-ray wavelengths, if produced, to vary with the Earth's magnetic field strength integrated along the line-of-sight direction for each observation, even if ALPs arriving on Earth have an isotropic distribution in the sky. To search for such ALP-origin X-ray radiation, we focus on the Suzaku X-ray data in the four deep fields, collected over eight years. These fields had been observed frequently by the Suzaku, but the magnetic field strengths varied with each observation depending on its location in the orbit. For a null hypothesis of ALP-induced photons, the diffuse X-ray background brightness estimated from the same field should not show any dependence on the integrated magnetic field. This is the signal we will search for in this paper. Suzaku data is suitable for our purpose, because the satellite, compared to the XMM-Newton 1 or Chandra 2 , has lower background due to its low-altitude orbit that prevents cosmic rays from penetrating the X-ray detectors [6].
Our study is somewhat similar to Fraser et al., (2014) [7], which claimed a detection of seasonal variations in the XMM-Newton X-ray data. The work claimed that the X-ray flux modulation at a level 4.6 × 10 −12 ergs s −1 cm −2 deg −2 in 2-6 keV might be due to a conversion of solar axions by their interaction with the Earth's magnetic field [also see 8,9]. However, Roncadelli and Tavecchio, (2015) [10] claimed that the XMM-Newton satellite which never points toward the Sun cannot observe such ALP-induced photons originating from solar axions due to momentum conservation.
The structure of this paper is as follows. In Section 2, we briefly review basics of the inverse Primakoff effect and how photons can be induced by the effect from ALPs that are created by dark matter in the expanding universe. In Section 3 we show the main results of this paper using the Suzaku data, when combined with data of the Earth's magnetic field at each orbit of the Suzaku satellite at each observation. Section 4 cotains the discussion and conclusion.

Process of photon emission from ALPs
In this section we describe a mechanism of photon emission from ALPs via the interaction with magnetic fields. To do this, we consider a model in which dark matter, which fills up space of the universe, preferentially decays into ALPs. This model is an example of moduli dark matter model in a string-theory-inspired scenario. When a dark matter particle decays into two ALPs, i.e. DM → 2ALPs, each ALP has a monochromatic energy E a = m φ /2, where m φ is the mass of the dark matter particle. The emissivity of DM → 2ALPs decaying process is given in terms of the energy density of dark matter, ρ φ (r), the decay rate, Γ φ→2a , and m φ as Considering the spatial distribution of dark matter along the line-of-sight direction, the ALP intensity, I a,line [counts s −1 cm −2 sr −1 ], is given as at E a = m φ /2, and S φ is the column density of dark matter in the line-of-sight direction [11], defined as In this case, the converted photon spectrum is a line emission.
1 http://sci.esa.int/xmm-newton/ 2 https://chandra.harvard.edu/ If dark matter is uniformly distributed in the universe, we would observe a continuum spectrum of the ALP intensity because free-streaming ALPs undergo a cosmological redshift in the expanding universe. Assuming light-mass ALPs, i.e. relativistic ALPs, produced by dark matter decay, a superposition of line spectra over different redshifts leads us to observe a continuum spectrum of ALPs [12,13]: where δ D (x) is the Dirac delta function, and the function f (x) is defined as In the above equation z is the redshift at decay, ρ φ 0 is the present energy density, H 0 is present the Hubble constant, Ω m0 and Ω Λ0 are the density parameters of non-relativistic matter and the cosmological constant, respectively. The spectral shape of ALPs is transcribed as a simple power-law function whose number index is a positive value of +1/2. In this case, the converted photon spectrum is also expected as a power-law function with a photon index of +1/2. The ALP-photon conversion probability in a vacuum with a magnetic field via inverse Primakoff effect is given in Ref. [14] as (2.8) Here, g aγγ is an ALP-photon coupling constant, m a and E a are mass and energy scales of ALP, and B ⊥ (x) is the perpendicular component of magnetic field to the ALP momentum direction, denoted as e a . The ALP-photon momentum transfer q is defined as Assuming the B ⊥ (x ) is uniform in the range 0 < x < L, we can write Equation (2.7) as: (2.10) In the limit of light ALP masses compared to the photon energy scale satisfying qL 1, 1 − cos (qL) (qL) 2 /2, and the conversion rate is simply given by under the coherence condition of qL < π → m a < 2πE a L (2.12) The following analysis uses Equation (2.10) to constrain the ALP-photon coupling constant. As shown above, the probability of ALP particles converting to photons proportional to (B ⊥ L) 2 in the light mass limit. Plugging typical values of the strength and coherent length scale of Earth's magnetic field, Equation (2.11) gives P a→γ 2.45 × 10 −21 g aγγ 10 −10 GeV −1 3 Analysis and results: A search for the correlation between residual Suzaku background radiation with the Earth's magnetic strength

Selection of blank sky observations from Suzaku archival data
To locate ALP-induced photons, we use the Suzaku X-ray data, and search for photons in the detector's field of view (FoV) depending on the integrated magnetic strength along the lineof-sight direction, (B ⊥ L) 2 . Because most X-ray data contains X-ray emission photons from targeted or unresolved sources, we need to study the X-ray diffuse background (XDB) in blank fields, and search for a residual signal in the background that is correlated with the magnetic strengths following the scaling of (B ⊥ L) 2 . The X-ray satellite Suzaku is suitable for this study because of its low instrumental background noise and the low background radiation from cosmic rays (compared to other X-ray satellites) due to its low-altitude Earth orbit; an altitude of ∼ 570 km and an inclination of 31 • from the Earth's equatorial plane, where the Earth's magnetic field prevents cosmic rays from penetrating the satellite's detectors [6]. Figure 1 is a schematic illustration of the Suzaku orbit and the Earth's magnetic field configuration. Even if the satellite observes the same field or the same angular direction-as denoted by the black dotted line-the integrated strength of perpendicular magnetic components along the line-of-sight direction varies with the satellite position. The Suzaku satellite orbits the Earth with a period of approximately 96 minutes and it causes a modulation of the integrated magnetic strength (B ⊥ L) 2 with the orbit, or when the target field is observed. Thus, we expect variations in the ALP-induced photons, if they exist, depending on the strength (B ⊥ L) 2 . We calculated the Earth's magnetic field every 60 seconds for each line-of-sight direction of a given target field using the software, International Geomagnetic Reference Field: the 12th generation (IGRF-12 [15]) for up to 6 times the Earth's radius (R E ), where typically B ∼ 10 −7 T. The right panel in Figure 1 shows a typical case of (B ⊥ L) 2 as a function of the satellite position or equivalently the observation time. It can be found that a typical value of (B ⊥ L) 2 is of order of 10 4 -10 5 T 2 m 2 , which is greater than that of terrestrial experiments such as the CAST experiment 3 . If we apply the non-oscillation condition of qL 1 (Equation (2.12)), the corresponding mass of ALP is limited to be m a ≤ µeV if we assume that the converted photons are in X-ray wavelengths. Note that we considered an oscillation regime of qL ∼ 1 to obtain constraints on the ALP-photon coupling constant.
To estimate the XDB spectrum, we consider blank sky data from four deep fields selected from the Suzaku archives as tabulated in Table 1. The selection criteria are as follows.
1. No bright sources in the FoV of Suzaku X-ray Imaging Spectrometer (XIS) [16], and compact sources in the FoV are already identified and can be masked in our analysis.

Sufficiently distant from regions of high X-ray diffuse emissions such as the North Polar
Spur.
4. Exposure time obtained by standard processing should be more than 200 ksec.
The above criteria are met by the following four fields, also shown in Table 1. First, we use the multiple observation data in the Lockman hole field, which is a famous region with minimum neutral hydrogen column density that was annually observed with Suzaku for calibration. We also use the data in the South Ecliptic Pole (SEP) and North Ecliptic Pole (NEP) fields. Finally, we use the data in the field of high latitude, the neutral hydrogen cloud or the so-called MBM16 field. We use the data reduction pipelines, the Ftools in HEAsoft version 6.16 and XSPEC version 12.8.2, to analyze the X-ray data in the four fields of Table 1, collected from the archive of Suzaku XIS. To avoid a possible contamination from high X-ray background, we removed data during the South Atlantic Anomaly region, Earth occultation, low elevation angle from the Earth's rim, and low cut-off-rigidity (COR; < 8 GV/c) regions. We stacked the X-ray image in the 0.5-7 keV band for each of the four fields, where we removed the point sources whose flux is larger than 1 × 10 −14 ergs s −1 cm −2 , with a radius of 1.5 arcminutes corresponding to encircled power function of 90% for Suzaku's mirror. We then calculated the (B ⊥ L) 2 every 60 seconds for each observation, as a function of the satellite position in orbit and observing line-of-sight direction. Figure 2 shows the distribution of (B ⊥ L) 2 in each of the four fields. We subdivided the data into four to six bins of the (B ⊥ L) 2 values, where the binning was determined so that each bin had almost the same photon statistics, as denoted by the different colored histograms in the figure.

An assessment of non-Xray background contamination
Before presenting the main results, we need to assess the level of non-Xray background (NXB) contamination in the data. Although the NXB of Suzaku is usually low, there could be a residual NXB contamination, up to 16-50% of the observed CXB in the 2-6 keV range. A part of NXB is due to fluorescence lines by materials such as Si, Al, Au, and Ni around the detector. These lines are distinguishable, if their emission lines are identified at their corresponding energy scales, in each XIS spectrum, as identified in Ref. [18]. Not only X-rays but charged particles too can produce pseudo-events in the CCD instrument. Pseudo events produced by particles originating from cosmic rays in orbit were studied by GEANT4 Monte-Carlo simulation. The continuum spectra by pseudo events are reproducible at an accuracy of 20% in its amplitude in each energy bin [19]. The production processes of pseudo events are well understood, but the input cosmic-ray flux varies by time and position of the satellite. The reproduction of the events was studied by using the event database collected during the periods when the FoV was blocked by the night side of the Earth (NTE) [20]. It was found that the intensity of the background could be estimated as a function of COR, and that the spectra were similar. They proposed a background estimation method to make a spectrum from ± 150 days of stacked night-Earth data weighted to reproduce the distribution of COR. They also found that the fluctuation of the background is larger than the simple Poisson statistics. Uncertainty of reproducibility for a typical 50 ksec exposure was reported to be 3.4%, although the expected statistical error by Poisson statistics is 1/10 [20]. This procedure is used as a standard background estimation for Suzaku and adapted as a HEASoft tool. In our analysis, we needed to sort the data by (B ⊥ L) 2 , which correlated with the COR. If the orbital position or (B ⊥ L) 2 is a potential control parameter of the NXB, it will affect our determination of the (B ⊥ L) 2 modulated signal. COR parameters used in the Suzaku analysis (defined as COR2 in the calibration database) are defined by the projected geographic coordinates, and calculated by the geomagnetic model on 2006 Apr. Actual COR would change gradually with time, and the cosmic ray flux is affected by the solar activity. We thus stepped into further NXB analysis of Suzaku, to evaluate possible range of background fluctuation, and to define further data reduction methods if needed.
We evaluated the fluctuation of input cosmic-ray flux by event rates at 12-15 keV of the XIS1. As the effective area of the X-ray mirror dropped rapidly above the Au L-edge below 1%, the event rate above 12 keV is considered to be an indicator of the cosmic ray flux. Due to the back-illumination structure, the background rate of XIS1 is higher than the other front-illuminated CCD, XIS0 and XIS3, and is more sensitive. Apparently, the fluctuation of the background count rate exceeds the Poisson statistics. In [20], the intrinsic fluctuation is evaluated as σ 2 calc − σ 2 Poisson . We evaluate the intrinsic fluctuation as follows: 1. Counting the number of events in 12-15 keV during the 60 sec for each COR range.
2. Calculating an mean of the count every 60 sec bins, denoted as µ, from the distribution of count rate as shown in histogram in Figure 3.
3. Assuming a certain value σ, and simulating a P NXB by Monte-Carlo method as shown in Equation (3.1) (lines in Figure 3).
where k is an observation frequency per interval, λ (k, σ) is an average number of events per interval, σ is an estimated systematic error by the variation of λ (k, σ), and µ is an average of observed NXB events.
4. Comparing the observed and simulated histograms by Pearson's chi-squared test and obtaining 95% upper limit for σ.
It assumes that the variations of the mean value of the count rates follow a Gaussian distribution, and the detected count follows a Poisson distribution. Thus, the observed count rate is expressed by a convolution of these functions. We use Pearson's chi-squared test to set a quantitative upper limit for the short term variability. A sample of these tests is shown in Figure 3. The 95% confidence range for the standard deviation of the Gaussian is obtained as 22-39% of the mean value.  The background rate anomaly of the geographical position is checked as follows. We divided the orbital position projected onto the Earth's surface by every 10 • in longitude and 5 • in latitude as defined as Loc_ID, and sorted them into 4 COR classes. For each observation, the NTE count rate in 2-5.6 keV for ±150 days, which is used by the standard background estimation, was accumulated on every Loc_ID. If the count rate at an Loc_ID is higher than the averaged value over the same COR class by 3σ, the events occurred at that Loc_ID were discarded from the spectral analysis.
After these data were reduced, the count rate in 2-6 keV before and after the standard background subtraction were plotted as a function of (B ⊥ L) 2 . We found that there is an negative correlation between the count rate and the (B ⊥ L) 2 , contrary to the ALP origin signal prediction. We checked the count rate of the upper discriminator (PIN-UD) of the Hard X-ray Detector (HXD) onboard Suzaku, which corresponds to energy deposited by protons approximatelly > 100 MeV [21],and found the same trend. The PIN-UD is affected by the radio-activation of HXD itself, and cannot be used to estimate the XIS background. We evaluated the correlation by a linear function fit, and decided that only those data satisfying (B ⊥ L) 2 ≥ 2 × 10 4 T 2 m 2 would be used for the analysis.

Spectral analysis for (B ⊥ L) 2 sorted data
In the spectral analysis, we assumed that celestial diffuse emission of each blank field is expressed by the sum of Cosmic X-ray Background (CXB), Milky Way Halo (MWH) emission, Solar Wind Charge eXchange (SWCX), Local Hot Bubble (LHB), and unknown High Temperature Component (HTC) as studied by previous works [11,22,23]. These are collectively called XDB. The surface brightness and spectral parameters for the celestial emission can be varied by the FoV in a reasonable range. The ALP signal has a power-law spectral shape with a photon index of +1/2, and with intensities proportional to (B ⊥ L) 2 . The NXB for each observation can be estimated by the standard background estimation method [20], but the intensities can also be varied within the fluctuation studied in the previous subsection.
Steps of spectral analysis for one observing direction are as follows; 0. Apply standard data reduction for XIS 0,1,3 of each observation ID (a unit of archival data, events from continuous pointing for the same observation direction), point source removal, Loc_ID selection and (B ⊥ L) 2 cut. Response matrices [24] and template NXB by standard method [20] are also prepared.
1. Accumulate the energy spectra in the 0.7-7.0 keV range for each XIS 0,3 and the 0.7-5.0 keV range for XIS 1 subtract the standard NXB, and fit them simultaneously by an empirical X-ray background model, obtained the best-fit values and errors with χ 2 statistics and C-statistics [25] in Xspec, and evaluate the validity of parameters.
2. Divide the energy spectra by (B ⊥ L) 2 values and fit them again simultaneously with C-statistics because of low photon statistics in each range. Check the consistency of spectral parameters obtained in step 1.
3. Add ALP emission model as a power-law function with a photon index of +1/2 and with a surface brightness proportional to the (B ⊥ L) 2 , and treat the background as a spectral model whose intensities can be tuned.
The fitting model describing the diffuse X-ray emission is similar to that used in [11]; it is shown by apec SWCX+LHB + phabs(apec MWH + power-law CXB + apec HTC ).
The APEC (Astrophysical Plasma Emission Code) [26] 4 is an emission model from collisional equilibrium and optically thin plasma installed in Xspec and applied to estimate the SWCX and LHB blend, MWH, and HTC. The temperature of apec in the SWCX and LHB blend was fixed to kT = 0.1 keV [27]. The typical temperature of the MWH is kT = 0.15-0.35 keV [11,22], a part of the blank sky spectra requires a HTC with kT = 0.6-0.9 keV to describe emission of approximately 0.9 keV [11]. The CXB was represented by a power-law emission model with a photon index of ∼ 1.4. The solar abundance table of apec model was given by [28]. Phabs describes the absorption by the Galactic interstellar medium, whose column density is fixed from the LAB(Leiden/Argentine/Bonn) survey [29] database. Steps 0 -1 are the standard spectral fitting procedure for Suzaku, and the parameters obtained in Step 1 were consistent with each other within 90% error, and with previous works like Sekiya et al. (2016) [11].
In step 2, we divided the data by (B ⊥ L) 2 . For example, in the case of Lockman hole, there were 11 observations, sorted by (B ⊥ L) 2 into 3 classes, and 3 CCDs, thus 99 spectra were fitted simultaneously with the same emission parameter. The number of the energy spectra for Lockman hole, MBM 16, SEP, and NEP, are 99, 36, 24, and 36, respectively.
The degrees of freedom in the spectral fit also increased, and the number of photons in each energy bin decreased. We applied C-statistics, which assumes that the data follows a Poisson distribution and uses the likelihood ratio to be minimized, and confirmed that the obtained parameters are consistent with Step 1. Because we would divide the spectra in later analysis, complex structure (mainly of Oxygen lines below 0.7 keV) are not well resolved. We only used the data in from the 0.7-7.0 keV range. Some components whose intensities were consistent with null were ignored by setting the intensities to 0.
In step 3, we added the ALP component whose surface brightness is proportional to (B ⊥ L) 2 . In usual spectral fitting like Steps 1-2, we used the spectra after subtraction of the estimated background. Here, we treated the NXB as one of the input models with a normalization factor, which can be variable in each observation ID, CCD, and (B ⊥ L) 2 class. In contrast, the parameters for celestial emission and the normalization of the ALP component at a fixed (B ⊥ L) 2 = 10 4 T 2 m 2 are common for the same FoV. The final fitting results are summarized in Table 2. We assumed the flux of the ALP both in negative and positive mathematically to evaluate proper error ranges, as shown in Figure 4.  Because the normalization of apec allows 0 within the error range, the temperature is not determined.
In all four observational directions, the surface brightness for the ALP components is consistent with 0 within a 90% confidence level. We also checked that the normalizations of the NXB model were within ± 40%, or the fluctuation studied in section 3.2. We made contour plots by surface brightness of ALP and CXB components between 2-6 keV, as shown in Figure 4. As the surface brightness varies with the index and normalization of an assumed power-law part, the contour is not smooth, owing to the steps in the parameter search. The limit obtained from the MBM16 observation is the lowest among these four fields and gives the tightest upper limit on the ALP flux: 1.6 × 10 −9 ergs s −1 cm −2 sr −1 normalized at 10 4 T 2 m 2 . An accumulated spectrum of all fitted data is shown in Figure 6 with averaged XDB and NXB model and the obtained upper limit for ALP. In Table 3, we also tabulated the center values of the CXB surface brightness and the upper limits of the ratio of ALP emission hidden in the CXB.   Table 2  To show the degeneracy among ALP and NXB normalization, we made a contour plot with the Lockman hole observation at one Obs. ID, one BL class, and one XIS, as shown in Figure 5. In the case of the Lockman hole, we have independently determined NXB normalization for all 11Obs. ID, 3 BL classes, and 3 XISs.

Discussion and Conclusions
We assumed that the cosmologically distributed ALPs would make a power-law with a photon index of +0.5 (dN/dE ∝ E +0.5 ) emission by the Earth's magnetosphere in proportion to the Ratio of NXB normalization Figure 5. The confidence contour between ratio of NXB normalization and surface brightness of ALP calculated from the photon index Γ ALP and normalization S ALP as shown in Table 2 obtained for Lockman hole at one Obs. ID, one (B ⊥ L) 2 class, and one XIS. 3 confidence levels: 68% (black), 90% (red) and 99% (green). Data all modlels MWH CXB ALP 99% upper limit NXB Figure 6. Accumulated spectrum used in the spectral fit for MBM16 direction. The spectrum is the sum of all Obs.ID, (B ⊥ L) 2 classes, and XISs. A response is weighted by the number of photons, and the NXB model is weighted by the exposure time after applying the normalization constant. Note that actual fitting was done with a set of energy spectra simultaneously, and no residuals are shown.
integrated (B ⊥ L) 2 in the FoV, and analyzed the data with Suzaku for four different directions. We did not detect any possible continuous emission from ALPs reported by previous similar studies [7]. We obtained the 99% upper limit of the X-ray surface brightness and flux  originating from ALPs in the 2.0-6.0 keV energy range as 1.6 × 10 −9 ergs s −1 cm −2 sr −1 , at (B ⊥ L) 2 = 10 4 T 2 m 2 , as shown in Table 3. It corresponds to 6-15% of the apparent CXB surface brightness in the 2-6 keV band, and is consistent with the idea that 80-90% of the CXB in the 2-8 keV band are resolved into point sources [30,31]. In other words, it could not be denied that 10-20% of the unresolved CXB could originate from the ALP converted to X-ray by the Earth atmosphere at Suzaku orbit. If we assume the dark matter density and decay rate, we can limit the ALP-photon coupling constant. By combining Equations (2.5), (2.6) and (2.13), the ALP-photon coupling constant, g aγγ , was constrained in the ALP mass range of m a < 2πE a /L ∼ 3.3 × 10 −6 eV to be g aγγ < 3.3 × 10 −7 GeV −1 m φ 10 keV as shown in Figure 7. Here, we assume a standard cosmology model with a mass density and a Hubble constant H 0 . The decay rate of the dark matter to the ALPs Γ φ→2a = 1/τ φ < 1/t 0 , where t 0 is the Hubble time. The factor f is defined by the Equation (2.6). Note that we neglect the reduction and anisotropy of the ALP flux due to interstellar and intergalactic magnetic fields. For the line emission search in the X-ray band, Sekiya et al. (2016) collected the longest exposure of 12 Msec from 10 years of Suzaku archival data, and obtained a 3σ upper limit for a narrow line emission between 1 and 7 keV to be 0.021 photons s −1 cm −2 sr −1 [11]. The ALP-photon conversion rate, P a→γ ∝ (B ⊥ L) 2 , was also computed by using IGRF-12 model every 60 seconds, and the averaged value was obtained to be (B ⊥ L) = 140 T m. It is larger than the value of 84 T m by CAST [32]. This value gives the upper limits of I a,line ·( g aγγ 10 −10 GeV −1 ) 2 < 4.4×10 14 axions s −1 cm −2 sr −1 in the 1.0 − 7.0 keV band. (4.2) A g aγγ can be also constrained as when the ALP density is connected with dark matter density around our galaxy. It is also shown in Figure 7 as galactic monochromatic ALP. In the plot, we consider the oscillation effect by Equation (2.10). This restriction of a physical parameter of ALPs is less strict than other experiments (e.g. CAST, ADMX), which assume a different axion and ALP model than this research. Nevertheless, it is important to note that we found these restrictions by using a new independent method from X-ray observations.   [33,34]