Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

A candidate proton cyclotron feature in the ultraluminous X-ray source NGC 4656 ULX-1

Nelson Cruz-Sanchez, Enzo A. Saavedra, Federico A. Fogantini, Federico García, Jorge A. Combi, Matteo Bachetti, Matteo Imbrogno, Lara Sidoli, Alessio Marino

Abstract

Ultraluminous X-ray sources represent extreme super-Eddington accretion regimes, and a subset is now known to host highly magnetized neutron stars. However, direct observational probes of their surface magnetic fields remain scarce. In this Letter, we report the detection of a narrow X-ray absorption feature at $3.29\pm0.02$ keV in the XMM$-$Newton/EPIC-pn spectrum of NGC 4656 ULX-1. The source exhibits a hard-ultraluminous state, while our timing analysis reveals a candidate pulsation at $\sim$0.9736 Hz, with a local significance of $5.5σ$ and a pulsed fraction of $\sim11\%$. The feature is robust against changes in continuum modeling and data-selection criteria, retaining a statistical significance of $\gtrsim3σ$ in Monte Carlo simulations. Interpreting the absorption as a proton cyclotron resonant scattering feature implies a local magnetic field of $B\sim(6-7)\times10^{14}$ G in the line-forming region. This value is consistent with strong magnetic fields anchored near the neutron star surface, even if the large-scale dipole is substantially weaker. Although we discuss electron cyclotron features and atomic transitions as possible alternatives, they appear less consistent with the observed phenomenology.

A candidate proton cyclotron feature in the ultraluminous X-ray source NGC 4656 ULX-1

Abstract

Ultraluminous X-ray sources represent extreme super-Eddington accretion regimes, and a subset is now known to host highly magnetized neutron stars. However, direct observational probes of their surface magnetic fields remain scarce. In this Letter, we report the detection of a narrow X-ray absorption feature at keV in the XMMNewton/EPIC-pn spectrum of NGC 4656 ULX-1. The source exhibits a hard-ultraluminous state, while our timing analysis reveals a candidate pulsation at 0.9736 Hz, with a local significance of and a pulsed fraction of . The feature is robust against changes in continuum modeling and data-selection criteria, retaining a statistical significance of in Monte Carlo simulations. Interpreting the absorption as a proton cyclotron resonant scattering feature implies a local magnetic field of G in the line-forming region. This value is consistent with strong magnetic fields anchored near the neutron star surface, even if the large-scale dipole is substantially weaker. Although we discuss electron cyclotron features and atomic transitions as possible alternatives, they appear less consistent with the observed phenomenology.
Paper Structure (16 sections, 7 figures, 2 tables)

This paper contains 16 sections, 7 figures, 2 tables.

Table of Contents

  1. Introduction
  2. Observation and data reduction
  3. Data analysis
  4. Timing analysis
  5. Spectral analysis
  6. Discussion and conclusions
  7. XMM--Newton observation
  8. Data reduction
  9. Timing analysis
  10. Phase-resolved spectroscopy
  11. NuSTAR observation
  12. Data reduction
  13. Timing analysis
  14. Spectral analysis
  15. Chandra observation
  16. ...and 1 more sections

Figures (7)

  • Figure 1: Broadband XMM--Newton and NuSTAR spectra of NGC 4656 ULX-1. EPIC-pn data are shown in black, MOS2 in orange, FPMA in red, and FPMB in blue. Top panel: Unfolded source spectra over 0.3--25 keV fit with the diskbb+powerlaw continuum. The light green vertical band marks the region around $E_{\rm line}\simeq3.3$ keV. Middle panel: Residuals of the diskbb+powerlaw fit without the absorption component ($\chi^2$/d.o.f. = 608.20/582). Bottom panel: Residuals after including the gabs component ($\chi^2$/d.o.f. = 587.74/579). Error bars correspond to 1$\sigma$ uncertainties.
  • Figure 2: Map of the $Z^2_2$ statistic in the frequency ($f$) versus frequency derivative ($\dot{f}$) plane around the candidate signal. The color scale represents the signal power, with the maximum ($Z^2_2 \approx 42$) marked by the green cross.
  • Figure 3: Left: Folded pulse profile at the candidate pulsation frequency, normalized by subtracting the mean count rate and dividing by the standard deviation. The re-binned profile (15 phase bins) is shown as black points with uncertainties. The red line represents a smoothed version of the folded light curve. The dashed horizontal line indicates the zero-mean level. Right: Residuals ($\Delta\chi$) of the phase-resolved spectra fit with an absorbed power-law model for the pulse-off (top) and pulse-on (bottom) intervals. The shaded green region indicates the expected energy range of the $\sim$3.3 keV absorption feature reported in the phase-averaged analysis.
  • Figure 4: EPIC-pn (left) and MOS2 (right) 0.3--10 keV images of the field around NGC 4656 ULX-1. The source extraction regions (solid circles, 20) are centered on the ULX position, while the dashed circles mark the background regions (40), chosen on the same CCD and free of contaminating sources. The colorbar scale indicates surface brightness in units of counts per pixel.
  • Figure 5: Distribution of the simulated $\Delta\chi^2$ values for the EPIC-pn spectra, obtained from $10^5$ Likelihood Ratio Test (LRT) simulations in XSPEC. Left panel: simulations performed under the null model ( powerlaw), comparing powerlaw versus gabs*powerlaw. Right panel: simulations performed under the null model ( ThComp(diskbb)), comparing ThComp(diskbb) versus gabs*ThComp(diskbb). In both panels, the red dashed lines mark the observed $\Delta\chi^2$ obtained from the real data when including the multiplicative Gaussian absorption component at $\sim$3.3 keV. The shaded histograms show the simulated $\Delta\chi^2$ distributions under the null hypothesis, while the vertical gray dashed lines indicate the nominal 1--3$\sigma$ thresholds for reference.
  • ...and 2 more figures