Pulse Profile Research Articles

Phase retrieval from integrated intensity of auto-convolution

Ultra-fast optical pulses are the most ephemeral sensing paradigm ever devised, examining events over incredibly brief timescales with broadband illumination. A consequence of sensing at timescales lower than a picosecond is that pulse characterization cannot be done with traditional analog-to-digital samplers and must be ascertained from integrating intensity sensors. Techniques for pulse characterization have been constructed using combinations of time-invariant and time-variant filter responses to create non-linear but invertible intensity datasets (Walmsley & Dorrer, 2009). In this paper, we develop a novel high-order phase retrieval technique to perform pulse characterization from a single-pixel integrating sensor measuring integrated intensity of auto-convolution (IIAC). We examine gradient descent’s ability to recover signals as a function of signal dimension and measurement count, and we demonstrate the effective use of iterative hard tensor thresholding as an initializer. Finally, we demonstrate IIAC recovery in a laboratory setting to recover the time profile of a complex laser pulse. We assert that the IIAC recovery solution demonstrated here simultaneously provides the optics community with a pulse characterization technique that scales to low-power microscopy systems and provides the optimization community with a physically motivated high-order phase retrieval problem enhanced by low-rank tensor processing.

Detection of pulsed X-ray emission from the isolated neutron star candidate eRASSU J131716.9–402647

The X-ray source eRASSU J131716.9–402647 was recently identified from observations with Spectrum Roentgen Gamma (SRG)/eROSITA as a promising X-ray dim isolated neutron star (XDINS) candidate on the premise of a soft energy distribution, absence of catalogued counterparts, and a high X-ray-to-optical flux ratio. Here, we report the results of a multi-wavelength observational campaign with XMM-Newton, NICER and the FORS2 instrument at the ESO-VLT. We found in both the XMM-Newton and NICER data that the X-ray emission is strongly pulsed at a period of 12.757 s (pulsed fraction pf = (29.1 ± 2.6)% in the 0.2–2 keV band). The pulse profile is double-humped, and the pulsed fraction increases with energy. The XMM-Newton and NICER epochs allow us to derive a 3σ upper limit of Ṗ ≤ 8 × 10−11 s s−1 on the spin-down rate of the neutron star. The source spectrum is well described by a purely thermal continuum, either a blackbody with kT ∼ 95 eV or a magnetised neutron star atmosphere model with kT ∼ 35 eV. Similarly to other thermally emitting isolated neutron stars, we found in either case strong deviations from the continuum, a broad absorption feature at energy ∼260 eV and a narrow one around 590 eV. The FORS2 instrument at ESO-VLT has not detected the optical counterpart (mR > 27.5 mag, 5σ detection limit), implying an X-ray-to-optical flux ratio of 104 at least. The properties of eRASSU J131716.9–402647 strongly resemble those of a highly magnetised isolated neutron star and favour an XDINS or high-B pulsar nature.

Long-term Study of the First Galactic Ultraluminous X-Ray Source Swift J0243.6+6124 Using NICER

We present the results obtained from detailed X-ray timing and spectral studies of X-ray pulsar Swift J0243.6+6124 during its giant and normal X-ray outbursts between 2017 and 2023 observed by the Neutron star Interior Composition Explorer (NICER). We focused on a timing analysis of the normal outbursts. A distinct break is found in the power density spectra of the source. The corresponding break frequency and slopes of the power laws around the break vary with luminosity, indicating a change in the accretion dynamics with the mass accretion rate. Interestingly, we detected quasiperiodic oscillations within a specific luminosity range, providing further insights into the underlying physical processes. We also studied the neutron star spin period evolution and a luminosity variation in the pulse profile during the recent 2023 outburst. The spectral analysis was conducted comprehensively for the giant and all other normal outbursts. We identified a double transition at luminosities of ≈7.5 × 1037 and 2.1 × 1038 erg s−1 in the evolution of continuum parameters like the photon index and cutoff energy with luminosity. This indicates three distinct accretion modes experienced by the source, mainly during the giant X-ray outburst. A soft blackbody component with a temperature of 0.08–0.7 keV is also detected in the spectra. The observed temperature undergoes a discontinuous transition when the pulsar evolves from a sub- to super-Eddington state. Notably, in addition to an evolving 6–7 keV iron line complex, a 1 keV emission line was observed during the super-Eddington state of the source, implying X-ray reflection from the accretion disk or outflow material.

Blind source separation for decomposing X-ray pulsar profiles

Accretion-powered X-ray pulsars offer a unique opportunity to study physics under extreme conditions. To fully exploit this potential, the interrelated problems of modelling radiative transport and the dynamical structure of the accretion flow must, however, be solved. This task is challenging both from a theoretical and observational point of view and is further complicated by a lack of direct correspondence between the properties of emission emerging from the neutron star and observed far away from it. In general, a mixture of emission from both poles of the neutron star viewed from different angles is indeed observed at some or even all phases of the pulse cycle. It is essential, therefore, to reconstruct the contributions of each pole to the observed flux in order to test and refine models describing the formation of the spectra and pulse profiles of X-ray pulsars. In this paper we propose a novel data-driven approach to address this problem using the pulse-to-pulse variability in the observed flux, and demonstrate its application to RXTE observations of the bright persistent X-ray pulsar Cen X-3. We then discuss the comparison of our results with previous work attempting to solve the same problem and how they can be qualitatively interpreted in the framework of a toy model describing emission from the poles of a neutron star.

Silicon x-ray backlighter improvement by targets with spike-like microstructures

In order to accurately probe high energy density matter states, it is vital to create powerful x-ray backlighters. One approach to create such x-ray sources is based on the usage of short-pulse, high-energy lasers, which greatly benefits from an optimization of the laser target coupling. Here, the spectral and temporal x-ray emission profiles of structured silicon targets with micron-sized spikes on the front surface are studied at laser intensities of 1017 W cm−2. The laser pulse length is varied between 1 and 20 ps with an energy of up to 1 kJ. The structured targets show an up to 13× enhancement of silicon Heα emission compared to flat foils with a well-defined, sharp emission pulse profile. Furthermore, the performance of the microstructured targets is compared to targets with a CH shield as well as foils irradiated with a UV prepulse.

Investigation of Nulling and Subpulse-drifting Properties of PSR J1649+2533 with FAST

We report the nulling and subpulse-drifting properties of PSR J1649+2533 with observations of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). The FAST observations reveal that the nulling fraction of this pulsar is about 20.9% ± 0.8% at 1250 MHz. The statistical study of the durations of the nulling and burst states shows that the burst states can persist for extended pulse periods, exceeding 100 periods in some cases, while the shortest lasts only a few pulse periods. The null states have a much shorter duration, with a maximum duration of less than 30 pulse periods. A comparative study between the pulse profiles of the first active pulse (FAP) and the last active pulse (LAP) shows that the pulse profiles of LAPs are stronger and wider than those of FAPs. An analysis of the two-dimensional fluctuation spectrum for the single-pulse stack indicates that the pulsar exhibits clear amplitude modulation and subpulse drifting. The periods are a P 3 = 2.5 ± 0.1 spin period and P 2 = 17.°0 ± 0.°5 at 1250 MHz, respectively. In addition, the multiband investigation shows that P 2 seems to increase with the increase of the observing frequency, i.e., P 2 ∝ ν 0.8±0.1. This seems to be caused by the increasing pulse-profile width with the frequency. The high-sensitivity FAST observations will enable a detailed understanding of the emission of this pulsar and provide important clues for theoretical studies of the radiation mechanism.

BGS Distribution Deformation in BOTDR Induced by Self-Phase Modulation

We report the numerical and experimental studies of the two-dimensional Brillouin gain spectrum (BGS) distribution deformation induced by the self-phase modulation in the Brillouin optical time domain reflectometry (BOTDR) with a 20.6 km sensing distance. The BGS distribution deformation is investigated by analyzing the evolution of the point spread function along the fiber in the two-dimensional model of the BOTDR. In the simulation and experimental results, the specific deformation degree of the BGS distribution induced by the self-phase modulation is related to the pump pulse profile, pump pulse peak power, BGS demodulation method, and detected scattered light component. By comprehensively analyzing the evolution of the point spread function induced by the self-phase modulation and using the image deconvolution, a typical BOTDR sensor with a 25 ns pump pulse reaches the 20 cm spatial resolution over the 20.6 km sensing fiber.

Identifying reliable periods in 2MASS J09213414–5939068, IGR J16167–4957, and V667 Pup

Detailed timing analyses of three cataclysmic variables, namely 2MASS J09213414−5939068, IGR J16167−4957, and V667 Pup are carried out using the long-baseline and high-cadence optical photometric data from the Transiting Exoplanet Survey Satellite (TESS). Periods of 908.12 ± 0.05 s and 990.10 ± 0.06 s are observed in the optical variation of 2MASS J09213414−5939068 that were not found in earlier studies and appear to be probable spin and beat periods of the system, respectively. The presence of multiple periods at spin, beat, and other sidebands indicates that 2MASS J09213414−5939068 likely belongs to an intermediate polar class of magnetic cataclysmic variables that seems to be accreted via a disc-overflow mechanism. Clear evidence of a period of 582.45 ± 0.04 s is found during the TESS observations of IGR J16167−4957, which can be interpreted as the spin period of the system. Strong modulation at this frequency supports its classification as an intermediate polar, where accretion may primarily be governed by a disc. The dominance of the spin pulse unveils the disc-fed dominance accretion in V667 Pup, but the detection of the previously unknown beat period of 525.77 ± 0.03 s suggests that a portion of the material is also accreted through a stream. Moreover, the double-peaked structure observed in the optical spin pulse profile of V667 Pup suggests the possibility of a two-pole accretion geometry, where each pole accretes at a different rate and is separated by 180°.

ERASSU J060839.5–704014: A double degenerate ultra-compact binary in the direction of the LMC

Context. During four all-sky surveys (eRASS1–4), eROSITA, the soft X-ray instrument aboard Spektrum-Roentgen-Gamma (SRG) detected a new supersoft X-ray source, eRASSU J060839.5−704014, in the direction of the Large Magellanic Cloud (LMC). Aims. We arranged follow-up observations in the X-ray and optical wavelengths and further searched in archival observations to reveal the nature of the object. Methods. Using X-ray observations with XMM-Newton we investigated the temporal and spectral behaviour of the source. Results. We discover pulsations at 374 s with a pulse profile consistent with 100% modulation. We identify two other periodicities in the eROSITA data, which we establish as aliases due to the sampling of the eROSITA light curve. We identify a multi-wavelength counterpart to the X-ray source in UVW1 and g, r, i, and z images obtained by the optical/UV monitor on XMM-Newton and the Dark Energy Camera at the Cerro Tololo Inter-American Observatory. The timing and spectral characteristics of the source are consistent with a double degenerate ultra-compact binary system in the foreground of the LMC. eRASSU J060839.5−704014 belongs to a rare class of AM CVns, which are important to study in the context of progenitors of SN Ia and for persistent gravitational wave detection. Conclusions. We identify eRASSU J060839.5−704014 as a new double degenerate ultra-compact binary located in the foreground of the LMC.

A Transition Discovered in the Subcritical Regime of 1A 0535+262

We present NICER observations of the accreting X-ray pulsar 1A 0535+262 during its faint state (≲6 × 1036 erg s−1), observed in several type I and type II outbursts. We discovered a transition of temporal and spectral properties around the luminosity L t = 3.3 × 1035 erg s−1, below which spectra are relatively soft and the pulse profiles have only a narrow peak. The spectra are harder and a secondary hump gradually appears in the pulse profiles when L ≳ L t. We discuss possible physical mechanisms for this transition, including different Comptonization seed photons, the disappearance of gas shocks on the neutron star surface, and the combination of plasma and vacuum polarization effects.

The AO327 Drift Survey Catalog and Data Release of Pulsar Detections

The AO327 drift survey for radio pulsars and transients used the Arecibo telescope from 2010 until its collapse in 2020. AO327 collected ∼3100 hr of data at 327 MHz with a time resolution of 82 μs and a frequency resolution of 24 kHz. While the main motivation for such surveys is the discovery of new pulsars and new, even unforeseen, types of radio transients, they also serendipitously collect a wealth of data on known pulsars. We present an electronic catalog of data and data products of 206 pulsars whose periodic emission was detected by AO327 and are listed in the Australia Telescope National Facility catalog of all published pulsars. The AO327 data products include dedispersed time series at full time resolution, average (“folded”) pulse profiles, Gaussian pulse profile templates, and an absolute phase reference that allows phase aligning the AO327 pulse profiles in a physically meaningful manner with profiles from data taken with other instruments. We also provide machine-readable tables with uncalibrated flux measurements at 327 MHz and pulse widths at 50% and 10% of the pulse peak determined from the fitted Gaussian profile templates. The AO327 catalog data set can be used in applications like population analysis of radio pulsars, pulse profile evolution studies in time and frequency, cone and core emission of the pulsar beam, scintillation, pulse intensity distributions, and others. It also constitutes a ready-made resource for teaching signal-processing and pulsar astronomy techniques.

Broad-band noise and quasi-periodic oscillation characteristics of the X-ray pulsar RX J0440.9+4431

ABSTRACT We present a comprehensive timing analysis on the Be/X-ray binary pulsar RX J0440.9+4431 using observations from NICER and Insight-HXMT during the 2022–2023 outburst. The power density spectrum (PDS) of RX J0440.9+4431 exhibits typical aperiodic variability in X-ray flux across a wide frequency range. During a supercritical accretion state, we detect quasi-periodic oscillations (QPOs) at 0.2–0.5 Hz in the light curves of five pulses for RX J0440.9+4431. The observed QPOs manifest during flares, while the flares appear at the peaks of the pulse profiles on a time-scale of seconds and are primarily caused by an increase in hard photons. These flares can be explained by increased material ingestion in the accretion column at a fixed phase, primarily generating hard photons. Alternatively, an increase in accretion rate, independent of phase, may result in highly beamed hard photons within the accretion column, causing the flares. We argue the origin of QPOs to instabilities within the accretion flow. Additionally, we find that the break frequencies in the noise power spectra align well with $\propto L_{\mathrm{x}}^{3 / 7}$ across three orders of magnitude in the luminosity, which points to a relatively strong magnetic field in RX J0440.9+4431, estimated to be ~1013 G.

The Thousand-Pulsar-Array programme on MeerKAT – XII. Discovery of long-term pulse profile evolution in seven young pulsars

ABSTRACT A number of pulsars are known to have profile evolution on time-scales of months, often correlated with spin-down rate changes. Here, we present the first result from 3 yr of monitoring observations from MeerKAT as part of the Thousand Pulsar Array programme. This programme obtains high-fidelity pulse profiles for ∼ 500 pulsars, which enabled the detection of subtle changes in seven sources not previously known to exhibit long-term profile evolution. A 2D Gaussian convolution is used to highlight correlated emission variability in both the pulse phase and observing epoch direction. Simulations show that for one additional source the observed profile variability is likely to originate from stochastic single-pulse shape variability (jitter). We find that it is common for long-term profile variability to be associated with changes in polarization fractions, but not with polarization position angle (PA) changes. PA changes are expected if emission height changes or precession is responsible for the profile variability. PSR J1741−3927 is the only pulsar in our sample that shows correlated PA variability, and this is associated with orthogonal polarization mode activity. For the six other pulsars limits on possible emission height changes and impact angle changes are derived. These limits are consistent with the small changes in the total intensity profile shape. None of the sources show detectable spin-down variability correlated with the emission changes, which are thought to be driven by magnetospheric current fluctuations. Therefore, the absence of correlated spin-down rate variability allows upper limits to be placed on changes in the magnetospheric charge density.

Studying the nulling, subpulse drifting, and moding in PSR J1921+1948 with the FAST

ABSTRACT We reported the observations of nulling, subpulse drifting, and moding of radio radiation in pulsar (PSR) B1918+19 at 1250 MHz with the Five-hundred-meter Aperture Spherical radio Telescope. The nulling fraction is 2.6 ± 0.1 per cent and no periodicity of nulling can be seen. We confirmed the existence of three different drift modes (A, B, C) and a disordered mode (N) at 1250 MHz. Unlike the other three modes where the second component of the average pulse profile dominates, the first component of the average pulse profile in mode C is brighter than the other components, the second component shifts forward and the fourth component shifts backward. The bidrifting phenomenon is observed in PSR B1918+19, the drifting of the first and second components is positive, and the drift direction of the fourth component is negative. The drifting rate of the drift bands composed of the first and second components has an obvious evolutionary trend. Mode B has a short duration and no clear trend can be seen. But for mode A, the drift rate of each drift band starts relatively fast, then slows down to a steady state, and finally increases slowly until it enters either null or mode N. Further analysis shows that the emergence and significant variation in the drifting period of multidrifting subpulse emission modes for PSR B1918+19 may be due to the aliasing effect. The interesting subpulse emission phenomenon of PSR B1918+19 at different frequencies provides a unique opportunity to understand the switching mechanism of the multidrift mode of the pulsars.

Observations of interstellar scattering of six pulsars using Polish LOFAR station PL611

ABSTRACT We present the preliminary results of 4 yr of observations of the scattering of pulsar radiation in the interstellar medium using the Low Frequency Array (LOFAR) PL611 station located in Lazy near Krakow. In this work we show the initial results for six pulsars from our observing campaign. We used the HBA antennas of the station, with a central frequency of 154 MHz and a 72 MHz bandwidth and we were able to detect the frequency-dependent change in the pulse profiles. Splitting the bandwidth into a number of separate sub-bands we obtained independent scatter time estimates, which allowed for the estimation of the scattering frequency scaling slope for each individual observation. Our project has been accumulating data for more than 4 yr and as a result we are in a position to study the time variability of the scattering parameters over this period. We detected significant changes in the scatter time and its frequency scaling for at least two of pulsars. The average values of the frequency scaling index for the six pulsars selected for this work are below the range of values predicted by the simple thin screen model of interstellar scattering. This is in accordance with previous results shown for both higher dispersion measure pulsars observed in the past, as well as the more recent LOFAR observations and other projects similar in the observing frequency. We also discuss the advantages of using individual LOFAR stations (or similar instruments) for this kind of research.

Multi-GeV wakefield acceleration in a plasma-modulated plasma accelerator.

We investigate the accelerator stage of a plasma-modulated plasma accelerator (P-MoPA) [Jakobsson et al., Phys. Rev. Lett. 127, 184801 (2021)0031-900710.1103/PhysRevLett.127.184801] using both the paraxial wave equationand particle-in-cell (PIC) simulations. We show that adjusting the laser and plasma parameters of the modulator stage of a P-MoPA allows the temporal profile of pulses within the pulse train to be controlled, which in turn allows the wake amplitude in the accelerator stage to be as much as 72% larger than that generated by a plasma beat-wave accelerator with the same total drive laser energy. Our analysis shows that Rosenbluth-Liu detuning is unimportant in a P-MoPA if the number of pulses in the train is less than ∼30, and that this detuning is also partially counteracted by increased red-shifting, and hence increased pulse spacing, towards the back of the train. An analysis of transverse mode oscillations of the driving pulse train is found to be in good agreement with 2D (Cartesian) PIC simulations. PIC simulations demonstrating energy gains of ∼1.5GeV (∼2.5GeV) for drive pulse energies of 2.4J (5.0J) are presented. Our results suggest that P-MoPAs driven by few-joule, picosecond pulses, such as those provided by high-repetition-rate thin-disk lasers, could accelerate electron bunches to multi-GeV energies at pulse repetition rates in the kilohertz range.

Discovery and timing of pulsar J2016+3711 in supernova remnant CTB 87 with FAST

ABSTRACT We report on our discovery of the radio pulsar, PSR J2016+3711, in supernova remnant (SNR) CTB 87, with a ∼10.8σ significance of pulses, which confirms the compact nature of the X-ray point source in CTB 87. It is the first pulsar discovered in SNRs using Five-hundred-meter Aperture Spherical radio Telescope (FAST). Its integrated radio pulse profile can be well described by a single component, with a width at 50 per cent of the peak flux density of about 28.1○ and an effective width of about 32.2○. The mean flux density at 1.25 GHz is estimated to be about 15.5 $\mu$Jy. Combined with the non-detection of the radio pulse at lower frequencies, the radio spectral index of the pulsar is constrained to be ≲2.3. We also present the timing solution based on 28 follow-up FAST observations. Our results reveal a period of 50.81 ms, period derivative of 7.2 × 10−14 s s−1, and dispersion measure of 428 pc cm−3. The strength of the equatorial surface magnetic dipole magnetic field is inferred to be about 1.9 × 1012 G. Using the ephemeris obtained from the radio observations, we searched Fermi-LAT data for gamma-ray pulsations but detected no pulsed signal. We also searched for radio pulses with FAST towards the X-ray counterpart of the gamma-ray binary HESS J1832−093 proximate to SNR G22.7−00.2 but found no signal.

Cyclotron line evolution revealed with pulse-to-pulse analysis in the 2020 outburst of 1A 0535+262

ABSTRACT We present a detailed analysis of the X-ray luminosity (LX) dependence of the cyclotron absorption line energy (Ecyc) for the X-ray binary pulsar 1A 0535+262 during its 2020 giant outburst based on pulse-to-pulse analysis. By applying this technique to high cadence observations of Insight-Hard X-ray Modulation Telescope, we reveal the most comprehensive Ecyc–LX correlation across a broad luminosity range of ∼(0.03–1.3) × 1038 erg s−1. Apart from the positive and negative correlations between cyclotron line energy and luminosity at LX ∼ (1–3) × 1037 erg s−1 and ∼(7–13) × 1037 erg s−1, which are expected from the typical subcritical and supercritical accretion regimes, respectively, a plateau in the correlation is also detected at ∼(3–7) × 1037 erg s−1. Moreover, at the lowest luminosity level (LX ≲ 1037 erg s−1), the positive Ecyc–LX correlation seems to be broken, and the pulse profile also occurs a significant transition. These discoveries provide the first complete view on the correlation between luminosity and the centriod energy of the cyclotron line, and therefore are relevant for understanding how accretion on to magnetized neutron stars depends on luminosity.

Effects of scattering in the accretion funnel on the pulse profiles of accreting millisecond pulsars

The hotspot emission of accreting millisecond pulsars (AMPs) undergoes scattering in the accretion flow between the disk inner radius and the neutron star surface. The scattering optical depth of the flow depends on the photon emission angle, which is a function of the pulse phase, and reaches its maximum when the hotspot is closest to the observer. At sufficiently large optical depths the observed pulse profile should develop a secondary minimum, the depth of which depends on the accretion rate and the emission geometry. Such a dip evolving with the accretion rate might explain the phase shift and pulse profile evolution observed in AMPs during outbursts. Accounting for scattering is important for accurate modeling of the AMP pulse profiles in order to improve the accuracy of determination of the neutron star parameters, such as their masses and radii. In this paper we present a simplified analytical model for the Thomson optical depth of the accretion funnel, and apply it to simulating the pulse profiles. We show that scattering in the accretion funnel has a significant effect on the pulse profiles at accretion rates of M ≳ 10−10 M⊙yr−1. Our model predicts a gradual evolution of the pulse profile with the accretion rate that appears to be consistent with the observations.

The ongoing spin-down episode of 4U 1626-67

Abstract We report the X-ray characteristics of the persistent X-ray pulsar 4U 1626-67 using simultaneous NuSTAR and NICER observations. The X-ray pulsar 4U 1626-67 has just encountered a torque reversal in 2023 and is presently in the spin-down state. We have examined the temporal and spectral characteristics of the source during its ongoing spin-down episode. The pulse profiles of the source are characterized by multiple substructures at lower energies and a wide asymmetric single-peaked structure at higher energies. The pulse fraction follows an overall increasing trend with energy. We confirm the existence of mHz quasi-periodic oscillation (QPO) exclusively during the current spin-down phase in all the observations. The source is spinning down at 0.00045(4) s yr−1. The broadband spectrum during this phase is described by empirical NPEX model and a soft blackbody component with kT ∼ 0.25 keV. In addition to the iron emission line, we also confirm the presence of cyclotron line at ∼ 36 keV. The source flux continues to decrease during the current spin-down phase, and the corresponding luminosity ∼ (3.3–4.9) × 1036 ergs s−1 lies in the intermediate range of accreting X-ray pulsars that may be associated with a hybrid accretion geometry. The magnetic field strengths estimated using the cyclotron line measurements and QPO frequency are consistent. The evolution of the spectral parameters relative to the pulsed phase is examined using phase-resolved spectroscopy.