Search NASA⌕ Search

Engineering topics

Zaven Arzoumanian

Publications and source records attributed to Zaven Arzoumanian.

At least 19 records

Highly Significant Detection of X-Ray Polarization from the Brightest Accreting Neutron Star Sco X-1

The Imaging X-ray Polarimetry Explorer measured with high significance the X-ray polarization of the brightest Z-source, Sco X-1, resulting in the nominal 2–8 keV energy band in a polarization degree of 1.0% ± 0.2% and a polarization angle of 8° ± 6° at a 90% confidence level. This observation was strictly simultaneous with observations performed by NICER, NuSTAR, and Insight-HXMT, which allowed for a precise characterization of its broadband spectrum from soft to hard X-rays. The source has been observed mainly in its soft state, with short periods of flaring. We also observed low-frequency quasiperiodic oscillations. From a spectropolarimetric analysis, we associate a polarization to the accretion disk at <3.2% at 90% confidence level, compatible with expectations for an electron scattering dominated optically thick atmosphere at the Sco X-1 inclination of ∼44°; for the higher-energy Comptonized component, we obtain a polarization of 1.3% ± 0.4%, in agreement with expectations for a slab of Thomson optical depth of ∼7 and an electron temperature of ∼3 keV. A polarization rotation with respect to previous observations by OSO-8 and PolarLight, and also with respect to the radio-jet position angle, is observed. This result may indicate a variation of the polarization with the source state that can be related to relativistic precession or a change in the corona geometry with the accretion flow.

Fabio La Monaca↗

Delayed Development of Cool Plasmas in X-Ray Flares from the Young Sun-like Star κ1 Ceti

The Neutron star Interior Composition Explorer (NICER) X-ray observatory observed two powerful X-ray flares equivalent to superflares from the nearby young solar-like star κ 1 Ceti in 2019. NICER follows each flare from the onset through the early decay, collecting over 30 counts s −1 near the peak, enabling a detailed spectral variation study of the flare rise. The flare in September varies quickly in ∼800 s, while the flare in December has a few times longer timescale. In both flares, the hard-band (2–4 keV) light curves show typical stellar X-ray flare variations with a rapid rise and slow decay, while the soft X-ray light curves, especially of the September flare, have prolonged flat peaks. The time-resolved spectra require two temperature plasma components at kT ∼0.3–1 and ∼2–4 keV. Both components vary similarly, but the cool component lags by ∼200 s with a four to six times smaller emission measure (EM) compared to the hot component. A comparison with hydrodynamic flare loop simulations indicates that the cool component originates from X-ray plasma near the magnetic loop footpoints that mainly cools via thermal conduction. The time lag represents the travel time of the evaporated gas through the entire flare loop. The cool component has a several times smaller EM than its simulated counterpart, suggesting a suppression of conductive cooling, possibly by the expansion of the loop cross-sectional area or turbulent fluctuations. The cool component's time lag and EM ratio provide important constraints on the flare loop geometry.

main sequence stars↗

X-Ray Burst and Persistent Emission Properties of the Magnetar SGR 1830-0645 in Outburst

We report on NICER X-ray monitoring of the magnetar SGR 1830−0645 covering 223 days following its 2020 October outburst, as well as Chandra and radio observations. We present the most accurate spin ephemerides of the source so far: ν = 0.096008680(2) Hz, $\dot{\n𝜈}=-6.2(1) x (10)^(-14) Hz/s, and significant second and third frequency derivative terms indicative of nonnegligible timing noise. The phase-averaged 0.8–7 keV spectrum is well fit with a double-blackbody (BB) model throughout the campaign. The BB temperatures remain constant at 0.46 and 1.2 keV. The areas and flux of each component decreased by a factor of 6, initially through a steep decay trend lasting about 46 days, followed by a shallow long-term one. The pulse shape in the same energy range is initially complex, exhibiting three distinct peaks, yet with clear continuous evolution throughout the outburst toward a simpler, single-pulse shape. The rms pulsed fraction is high and increases from about 40% to 50%. We find no dependence of pulse shape or fraction on energy. These results suggest that multiple hot spots, possibly possessing temperature gradients, emerged at outburst onset and shrank as the outburst decayed. We detect 84 faint bursts with NICER, having a strong preference for occurring close to the surface emission pulse maximum—the first time this phenomenon is detected in such a large burst sample. This likely implies a very low altitude for the burst emission region and a triggering mechanism connected to the surface active zone. Finally, our radio observations at several epochs and multiple frequencies reveal no evidence of pulsed or burst-like radio emission.

George Younes↗

A Comprehensive X-ray Report on AT2019wey

The Galactic low-mass X-ray binary AT2019wey (ATLAS19bcxp, SRGAJ043520.9+552226, SRGEJ043523.3+552234, ZTF19acwrvzk) was discovered as a new optical transient in Dec 2019, and independently as an X-ray transient in Mar 2020. In this paper, we present comprehensive NICER, NuSTAR, Chandra, Swift, and MAXI observations of AT2019wey from ~1 year prior to the discovery to the end of September 2020. AT2019wey appeared as a ~ 1 mCrab source and stayed at this flux density for several months, displaying a hard X-ray spectrum that can be modeled as a power-law with photon index Γ ~ 1.8. In June 2020 it started to brighten and reached ~20 mCrab in ~2 months. The inclination of this system can be constrained to I ≲ 30° by modelling the reflection spectrum. Starting from late-August (~ 59082 MJD), AT2019wey entered into the hard-intermediate state (HIMS), and underwent a few week-long timescale outbursts, where the brightening in soft X-rays is correlated with the enhancement of a thermal component. Low-frequency quasi-periodic oscillation (QPO) was observed in the HIMS. We detect no pulsation and in timing analysis of the NICER and NuSTAR data. The X-ray states and power spectra of AT2019wey are discussed against the landscape of low-mass X-ray binaries.

Individual X-rays↗

On the Impact of an Intermediate Duration X-ray Burst on the Accretion Environment in IGR J17062-6143

We report on a spectroscopic analysis of the X-ray emission from IGR J17062−6143 in the after-math of its June 2020 intermediate duration Type I X-ray burst. Using the Neutron Star Interior Composition Explorer, we started observing the source three hours after the burst was detected with MAXI/GSC, and monitored the source for the subsequent twelve days. We observed the tail end of the X-ray burst cooling phase, and find that the X-ray flux is severely depressed relative to its historic value for a three day period directly following the burst. We interpret this intensity dip as the inner accretion disk gradually restoring itself after being perturbed by the burst irradiation. Superimposed on this trend we observed a 1.5 d interval during which the X-ray flux is sharply lower than the wider trend. This drop in flux could be isolated to the non-thermal components in the energy spectrum, suggesting that it may be caused by an evolving corona. Additionally, we detected a 3.4 keV absorption line at 6.3σsignificance in a single 472 s observation while the burst emission was still bright. We tentatively identify the line as a gravitationally redshifted absorption line from burning ashes on the stellar surface, possibly associated with 40Ca or 44Ti.

Peter Bult↗

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars↗

The NANOGrav 12.5 Year Data Set: Monitoring Interstellar Scattering Delays

We extract interstellar scintillation parameters for pulsars observed by the NANOGrav radio pulsar timing program. Dynamic spectra for the observing epochs of each pulsar were used to obtain estimates of scintillation timescales, scintillation bandwidths, and the corresponding scattering delays using a stretching algorithm to account for frequency-dependent scaling. We were able to measure scintillation bandwidths for 28 pulsars at1500 MHz and 15 pulsars at 820 MHz. We examine scaling behavior for 17 pulsars and find power-law indices ranging from−0.7 to−3.6, though these may be biased shallow due to insufficient frequency resolution at lower frequencies. We were also able to measure scintillation timescales for six pulsars at 1500 MHz and seven pulsars at820 MHz. There is fair agreement between our scattering delay measurements and electron-density model predictions for most pulsars. We derive interstellar scattering-based transverse velocities assuming isotropic scattering and a scattering screen halfway between the pulsar and Earth. We also estimate the location of the scattering screens assuming proper motion and interstellar scattering-derived transverse velocities are equal. We find no correlations between variations in scattering delay and either variations in dispersion measure or flux density. For most pulsars for which scattering delays are measurable, we find that time-of-arrival uncertainties for a given epoch are larger than our scattering delay measurements, indicating that variable scattering delays are currently subdominant in our overall noise budget but are important for achieving precisions of tens of nanoseconds or less.

Jacob E Turner↗

The NANOGrav 11yr Data Set: Limits on Supermassive Black Hole Binaries in Galaxies within 500 Mpc

Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At subparsec separations, binaries become strong sources of low-frequency gravitational waves (GWs), targeted by Pulsar Timing Arrays. We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11 yr data set to place constraints on putative SMBHBs in nearby massive galaxies. We compiled a comprehensive catalog of ∼44,000 galaxies in the local universe (up to redshift ∼0.05) and populated them with hypothetical binaries, assuming that the total mass of the binary is equal to the SMBH mass derived from global scaling relations. Assuming circular equal-mass binaries emitting at NANOGrav's most sensitive frequency of 8 nHz, we found that 216 galaxies are within NANOGrav's sensitivity volume. We ranked the potential SMBHBs based on GW detectability by calculating the total signal-to-noise ratio such binaries would induce within the NANOGrav array. We placed constraints on the chirp mass and mass ratio of the 216 hypothetical binaries. For 19 galaxies, only very unequal-mass binaries are allowed, with the mass of the secondary less than 10% that of the primary, roughly comparable to constraints on an SMBHB in the Milky Way. However, we demonstrated that the (typically large) uncertainties in the mass measurements can weaken the upper limits on the chirp mass. Additionally, we were able to exclude binaries delivered by major mergers (mass ratio of at least 1/4) for several of these galaxies. We also derived the first limit on the density of binaries delivered by major mergers purely based on GW data.

Zaven Arzoumanian↗

Constraining the Neutron Star Mass–Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes

We describe the X-ray pulse profile models we use and how we use them to analyze Neutron Star Interior Composition Explorer(NICER)observations of rotation-powered millisecond pulsars to obtain information about the mass–radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we analyzed in Miller et al. and Riley et al. and describe a cross-verification of computations of the pulse profiles of a star with R/M 3, in case stars this compact need to be considered in future analyses. We also present our early cross-verification efforts of the parameter estimation procedures used by Miller et al. and Riley et al. by analyzing two distinct synthetic data sets. Both codes yielded credible regions in the mass–radius plane that are statistically consistent with one another, and both gave posterior distributions for model parameter values consistent with the values that were used to generate the data. We also summarize the additional tests of the parameter estimation procedure of Miller et al. that used synthetic pulse profiles and the NICER pulse profile of PSR J0030+0451. We then illustrate how the precision of mass and radius estimates depends on the pulsar’s spin rate and the size of its hot spot by analyzing four different synthetic pulse profiles. Finally, we assess possible sources of systematic error in the estimates made using this technique, some of which may warrant further investigation.

Slavko Bogdanov↗

Astrophysics Milestones for Pulsar Timing Array Gravitational-wave Detection

The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its12.5 yr pulsar timing array data set, with median characteristic strain amplitude at periods of a year of A(yr) = 1.92(+0.75,-0.55) x 10^(-15). However, evidence for the quadrupolar Hellings & Downs interpulsar correlations, which are characteristic of gravitational-wave signals, was not yet significant. We emulate and extend the NANOGrav data set, injecting a wide range of stochastic gravitational-wave background(GWB)signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones.(I)Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15–17 yr of data, i.e., an additional 2–5 yr from the 12.5 yr data set.(II)At the initial detection, we expect a fractional uncertainty of 40% on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black hole binary origin from some models predicting more exotic origins.(III)Similarly, the measured GWB amplitude will have an uncertainty of 44% upon initial detection, allowing us to arbitrate between some population models of supermassive black hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20 yr of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar timing array data sets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier.

Nihan S. Pol↗

Enhanced X-Ray Emission Coinciding With Giant Radio Pulses From the Crab Pulsar

Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars that last a few microseconds and are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside of radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percentage points during GRPs. We observed the Crab Pulsar simultaneously at x-ray and radio wavelengths, finding enhancement of the x-ray emission by 3.8 ± 0.7% (a 5.4sdetection)coinciding with GRPs. This implies that the total emitted energy from GRPs is tens to hundreds of times higher than previously known. We discuss the implications for the pulsar emission mechanism and extragalactic fast radio bursts.

Teruaki Enoto↗

The RS CVn Type star GT Mus Shows Most Energetic X-ray Flares Throughout the 2010s

We report that the RS CVn-type star GT Mus (HR4492, HD101379+HD101380) was the most active star in the X-ray sky in the last decade in terms of the scale of recurrent energetic flares. We detected eleven flares from GT Mus in eight years of observations with MAXI (Monitor of All-sky X-ray Image)from 2009 August to 2017 August. The detected flare peak luminosities were 1–4×1033erg s−1in the2.0–20.0 keV band for its distance of 109.6 pc. The durations of the flares ranged from 1–4 days. The flare energies are calculated to be 1–8×1038erg, which are at the upper end of the observed stellar flare. We carried out X-ray follow-up observations for one of the eleven flares with NICER (Neutron star Interior Composition Explorer) on 2017 July 18 and found that the flare cooled quasi-statically. On the basis of a quasi-static cooling model, the flare loop length is derived to be 4×1012cm (or60R). The electron density is derived to be1×1010cm−3, which is consistent with the typical value of solar and stellar flares (1010−11cm−3).

Ryo Sasaki↗

NICER Discovery of Millisecond X-Ray Pulsations and an Ultracompact Orbit in IGR J17494-3030

We report the detection of 376.05 Hz (2.66 ms) coherent X-ray pulsations in NICER observations of a transient outburst of the low-mass X-ray binary IGR J17494−3030 in 2020 October/November. The system is an accreting millisecond X-ray pulsar in a 75-minute ultracompact binary. The mass donor is most likely a ≃0.02 Mꙩ finite-entropy white dwarf composed of He or C/O. The fractional rms pulsed amplitude is 7.4%, and the soft (1–3 keV) X-ray pulse profile contains a significant second harmonic. The pulsed amplitude and pulse phase lag (relative to our mean timing model) are energy dependent, each having a local maximum at 4 and 1.5 keV, respectively. We also recovered the X-ray pulsations in archival 2012 XMM-Newton observations, allowing us to measure a long-term pulsar spin-down rate of ύ= -2.1(7) x (10)^(-14) Hz/s and to infer a pulsar surface dipole magnetic field strength of ≃10^(9)G. We show that the mass transfer in the binary is likely nonconservative, and we discuss various scenarios for mass loss from the system.

Mason Ng↗

The X-ray Bursts of XTE J1739-285: A NICER Sample

In this work, we report on observations with the Neutron Star Interior Composition Explorer of the known neutron star X-ray transient XTE J1739–285. We observed the source in 2020 February and March, finding it in a highly active bursting state. Across a 20 day period, we detected 32 thermonuclear X-ray bursts, with an average burst recurrence time of-+2.0 hr0.30.4. A timing and spectral analysis of the ensemble of X-ray bursts reveals homogeneous burst properties, evidence for short-recurrence time bursts, and the detection of a 386.5 Hz burst oscillation candidate. The latter is especially notable, given that a previous study of this source claimed a 1122 Hz burst oscillation candidate. We did not find any evidence of variability near 1122 Hz and instead find that the 386.5 Hz oscillation is the more prominent signal of the two burst oscillation candidates. Hence, we conclude it is unlikely that XTE J1739–285 has a submillisecond rotation period.

X-ray burst↗

Multimessenger Gravitational-wave Searches with Pulsar Timing Arrays: Application to 3C 66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and emit low-frequency gravitational radiation in the process. In this paper, we consider the galaxy 3C 66B, which was used as the target of the first multimessenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational-wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C 66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C 66B to less than (1.65±0.02)×109Me using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data over “blind” pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

Zaven Arzoumanian↗

Return of the Big Glitcher: NICER Timing and Glitches of PSR J0537-6910

PSR J0537−6910, also known as the Big Glitcher, is the most prolific glitching pulsar known, and its spin-induced pulsations are only detectable in X-ray. We present results from analysis of 2.7 yr of NICER timing observations, from 2017 August to 2020 April. We obtain a rotation phase-connected timing model for the entire time span, which overlaps with the third observing run of LIGO/Virgo, thus enabling the most sensitive gravitational wave searches of this potentially strong gravitational wave-emitting pulsar. We find that the short-term braking index between glitches decreases towards a value of 7 or lower at longer times since the preceding glitch. By combining NICER and RXTE data, we measure a long-term braking index n = −1.25 ± 0.01. Our analysis reveals eight new glitches, the first detected since 2011, near the end of RXTE, with a total NICER and RXTE glitch activity of 8.88 × 10−7 yr−1. The new glitches follow the seemingly unique time-to-next-glitch–glitch-size correlation established previously using RXTE data, with a slope of 5 d μHz−1. For one glitch around which NICER observes 2 d on either side, we search for but do not see clear evidence of spectral nor pulse profile changes that may be associated with the glitch.

Wynn C G Ho↗

Timing the Pulsations of The Accreting Millisecond Pulsar SAX J1808.4-3658 During Its 2019 Outburst

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

Peter Bult↗