Source author record

Craig B. Markwardt

Craig B. Markwardt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2026arXiv

A simple, flexible method for timing cross-calibration of space missions

The timing (cross-)calibration of astronomical instruments is often done by comparing pulsar times-of-arrival (TOAs) to a reference timing model. In high-energy astronomy, the choice of solar system ephemerides and source positions used to barycenter the photon arrival times has a significant impact on the procedure, requiring a full reprocessing of the data each time a new convention is used. Our method, developed as part of the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC), adapts an existing pulsar solution to arbitrary JPL ephemerides and source positions by simulating geocentric TOAs and refitting timing models (implemented with PINT). We validate the procedure and apply it to thousands of observations of the Crab pulsar from 15 missions spanning 1996--2025, demonstrating inter-ephemeris TOA consistency at the $\lesssim5 μ$s level, using the DE200/FK5-based Jodrell Bank Monthly Ephemeris as a common reference. We release the TOAExtractor open-source tool and a TOA database to support future calibration and scientific studies. Instrument timing performance is broadly consistent with mission specifications; the X-ray-to-radio phase offset varies with energy and time at a level that is marginally consistent with the uncertainties of the radio ephemeris, motivating coordinated multiwavelength follow-up.

preprint2022arXiv

The NASA Multi-Messenger Astrophysics Science Support Center (MOSSAIC)

The era of multi-messenger astrophysics has arrived, leading to key new discoveries and revealing a need for coordination, collaboration, and communication between world-wide communities using ground and space-based facilities. To fill these critical needs, NASA's Goddard Space Flight Center and Marshall Space Flight Center are jointly proposing to establish a virtual Multi-Messenger Astrophysics Science Support Center that focuses entirely on community-directed services. In this article, we describe the baseline plan for the virtual Support Center which will position the community and NASA as an Agency to extract maximum science from multi-messenger events, leading to new breakthroughs and fostering increased coordination and collaboration.

preprint2021arXiv

A Detection of Red Noise in PSR J1824$-$2452A and Projections for PSR B1937+21 using NICER X-ray Timing Data

We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations ("red noise") in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar medium. Our Bayesian search methodology yields strong evidence (natural log Bayes factor of $9.634 \pm 0.016$) for red noise in PSR J1824$-$2452A, but is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulate longer and higher precision X-ray datasets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5-year mission with a time-of-arrival (TOA) error of 2 microseconds and an observing cadence of 20 observations per month compared to the 5 microsecond TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that an injected stochastic gravitational wave background (GWB) with an amplitude of $A_{\rm GWB}=2\times10^{-15}$ and spectral index of $γ_{\rm GWB}=13/3$ can be detected in a pulsar with similar TOA precision to PSR B1937+21, but with no additional red noise, in a 10-year mission that observes the pulsar 15 times per month and has an average TOA error of 1 microsecond.

preprint2021arXiv

Timing Calibration of the NuSTAR X-ray Telescope

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission is the first focusing X-ray telescope in the hard X-ray (3-79 keV) band. Among the phenomena that can be studied in this energy band, some require high time resolution and stability: rotation-powered and accreting millisecond pulsars, fast variability from black holes and neutron stars, X-ray bursts, and more. Moreover, a good alignment of the timestamps of X-ray photons to UTC is key for multi-instrument studies of fast astrophysical processes. In this Paper, we describe the timing calibration of the NuSTAR mission. In particular, we present a method to correct the temperature-dependent frequency response of the on-board temperature-compensated crystal oscillator. Together with measurements of the spacecraft clock offsets obtained during downlinks passes, this allows a precise characterization of the behavior of the oscillator. The calibrated NuSTAR event timestamps for a typical observation are shown to be accurate to a precision of ~65 microsec.

preprint2020arXiv

A NICER View of Spectral and Profile Evolution for Three X-ray Emitting Millisecond Pulsars

We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at $\gt 3σ$ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821-24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821-24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission.

preprint2020arXiv

Hard X-ray properties of radio-selected blazars

Hard X-ray properties of beamed AGN have been published in the 105-month Swift/BAT catalog, but there have not been any studies carried out so far on a well-defined, radio-selected sample of low-peaked blazars in the hard X-ray band. Using the statistically complete MOJAVE-1 sample, we aim to determine the hard X-ray properties of radio-selected blazars, including the enigmatic gamma-ray-faint type. Also, we aim to determine the contribution of radio-selected low-peaked blazars to the diffuse CXB. We determined photon indices, fluxes, and luminosities in the range of 20 keV - 100 keV of blazars and other extragalactic jets from the MOJAVE-1 sample, derived from the 105-month Swift/BAT survey. We calculated log N-log S distributions and the luminosity functions. The majority of the MOJAVE-1 blazars are found to be hard X-ray emitters albeit many at low count rates. The log N-log S distribution for the hard X-ray emission of radio-selected blazars is clearly non-Euclidean, in contrast to the radio flux density distribution. Approximately 0.2% of the CXB in the 20 keV - 100 keV band can be resolved into MOJAVE-1 blazars. The peculiar log N-log S distribution disparity might be attributed to different evolutionary paths in the X-ray and radio bands, as tested by luminosity-function modeling. X-ray variability can be ruled out as the dominant contributor. Low-peaked blazars constitute an intrinsically different source population in terms of CXB contribution compared to similar studies of X-ray-selected blazars. The hard X-ray flux and spectral index can serve as a good proxy for the gamma-ray detection probability of individual sources. Future observations combining deep X-ray survey, for example, with eROSITA, and targeted gamma-ray observations with CTA can benefit strongly from the tight connection between these high-energy bands for the different blazar sub-classes.

preprint2016arXiv

Getting NuSTAR on target: predicting mast motion

The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing high energy (3-79 keV) X-ray observatory operating for four years from low Earth orbit. The X-ray detector arrays are located on the spacecraft bus with the optics modules mounted on a flexible mast of 10.14m length. The motion of the telescope optical axis on the detectors during each observation is measured by a laser metrology system and matches the pre-launch predictions of the thermal flexing of the mast as the spacecraft enters and exits the Earths shadow each orbit. However, an additional motion of the telescope field of view was discovered during observatory commissioning that is associated with the spacecraft attitude control system and an additional flexing of the mast correlated with the Solar aspect angle for the observation. We present the methodology developed to predict where any particular target coordinate will fall on the NuSTAR detectors based on the Solar aspect angle at the scheduled time of an observation. This may be applicable to future observatories that employ optics deployed on extendable masts. The automation of the prediction system has greatly improved observatory operations efficiency and the reliability of observation planning.

preprint2015arXiv

Calibration of the NuSTAR High Energy Focusing X-ray Telescope

We present the calibration of the \textit{Nuclear Spectroscopic Telescope Array} (\nustar) X-ray satellite. We used the Crab as the primary effective area calibrator and constructed a piece-wise linear spline function to modify the vignetting response. The achieved residuals for all off-axis angles and energies, compared to the assumed spectrum, are typically better than $\pm 2$\% up to 40\,keV and 5--10\,\% above due to limited counting statistics. An empirical adjustment to the theoretical 2D point spread function (PSF) was found using several strong point sources, and no increase of the PSF half power diameter (HPD) has been observed since the beginning of the mission. We report on the detector gain calibration, good to 60\,eV for all grades, and discuss the timing capabilities of the observatory, which has an absolute timing of $\pm$ 3\,ms. Finally we present cross-calibration results from two campaigns between all the major concurrent X-ray observatories (\textit{Chandra}, \textit{Swift}, \textit{Suzaku} and \textit{XMM-Newton}), conducted in 2012 and 2013 on the sources 3C\,273 and PKS\,2155-304, and show that the differences in measured flux is within $\sim$10\% for all instruments with respect to \nustar.

preprint2014arXiv

NuSTAR Spectroscopy of Multi-Component X-ray Reflection from NGC 1068

We report on observations of NGC1068 with NuSTAR, which provide the best constraints to date on its $>10$~keV spectral shape. We find no strong variability over the past two decades, consistent with its Compton-thick AGN classification. The combined NuSTAR, Chandra, XMM-Newton, and Swift-BAT spectral dataset offers new insights into the complex reflected emission. The critical combination of the high signal-to-noise NuSTAR data and a spatial decomposition with Chandra allow us to break several model degeneracies and greatly aid physical interpretation. When modeled as a monolithic (i.e., a single N_H) reflector, none of the common Compton-reflection models are able to match the neutral fluorescence lines and broad spectral shape of the Compton reflection. A multi-component reflector with three distinct column densities (e.g., N_H~1.5e23, 5e24, and 1e25 cm^{-2}) provides a more reasonable fit to the spectral lines and Compton hump, with near-solar Fe abundances. In this model, the higher N_H components provide the bulk of the Compton hump flux while the lower N_H component produces much of the line emission, effectively decoupling two key features of Compton reflection. We note that ~30% of the neutral Fe Kalpha line flux arises from >2" (~140 pc), implying that a significant fraction of the <10 keV reflected component arises from regions well outside of a parsec-scale torus. These results likely have ramifications for the interpretation of poorer signal-to-noise observations and/or more distant objects [Abridged].

preprint2014arXiv

The Black Hole Binary V4641 Sagitarii: Activity in Quiescence & Improved Mass Determinations

We examine ~10 years of photometric data and find that the black hole X-ray binary V4641 Sgr has two optical states, passive and active, during X-ray quiescence. The passive state is dominated by ellipsoidal variations and is stable in the shape and variability of the light curve. The active state is brighter and more variable. Emission during the active state varies over the course of the orbital period and is redder than the companion star. These optical/infrared states last for weeks or months. V4641 Sgr spends approximately 85% of X-ray quiescence in the passive state and 15% in the active. We analyze passive colors and spectroscopy of V4641 Sgr and show that they are consistent with a reddened B9III star (with E(B-V) = 0.37 +/- 0.19) with little or no contribution from the accretion disk. We use X-ray observations with an updated ephemeris to place an upper limit on the duration of an X-ray eclipse of <8.3 deg in phase (~1.6 hours). High resolution spectroscopy yields a greatly improved measurement of the rotational velocity of the companion star of V_rot_sin(i) = 100.9 +/- 0.8 km s^-1. We fit ellipsoidal models to the passive state data and find an inclination angle of i = 72.3 +/- 4.1 deg, a mass ratio of Q = 2.2 +/- 0.2, and component masses for the system of M_BH = 6.4 +/- 0.6 M_sun and M_2 = 2.9 +/- 0.4 M_sun. Using these values we calculate an updated distance to V4641 Sgr of 6.2 +/- 0.7 kpc.

preprint2013arXiv

NuSTAR discovery of a 3.76-second transient magnetar near Sagittarius A*

We report the discovery of 3.76-s pulsations from a new burst source near Sgr A* observed by the NuSTAR Observatory. The strong signal from SGR J1745-29 presents a complex pulse profile modulated with pulsed fraction 27+/-3 % in the 3-10 keV band. Two observations spaced 9 days apart yield a spin-down rate of Pdot = (6.5+/-1.4)x10^{-12}. This implies a magnetic field B = 1.6x10^14 G, spin-down power Edot = 5x10^33 erg/s, and characteristic age P/2Pdot = 9x10^3 yr, for the rotating dipole model. However, the current Pdot may be erratic, especially during outburst. The flux and modulation remained steady during the observations and the 3-79 keV spectrum is well fitted by a combined blackbody plus power-law model with temperature kT_BB = 0.96+/-0.02 keV and photon index = 1.5+/-0.4, respectively. The neutral hydrogen column density (nH ~ 1.4x10^23 cm^{-2}) measured by NuSTAR and Swift suggests that SGR J1745-29 is located at or near the Galactic Center. The lack of an X-ray counterpart in the published Chandra survey catalog sets a quiescent 2-8 keV luminosity limit of Lx ~< 10^32 erg/s. The bursting, timing, and spectral properties indicate a transient magnetar undergoing an outburst with 2-79 keV luminosity up to 3.5x10^35 erg/s for a distance of 8 kpc. SGR J1745-29 joins a growing subclass of transient magnetars, indicating that many magnetars in quiescence remain undetected in the X-ray band or have been detected as high-B radio pulsars. The peculiar location of SGR J1745-29 has important implications for the formation and dynamics of neutron stars in the Galactic Center region.

preprint2013arXiv

NuSTAR observations of magnetar 1E 1841-045

We report new spectral and temporal observations of the magnetar 1E 1841-045 in the Kes 73 supernova remnant obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR). Combined with new Swift and archival XMM-Newton and Chandra observations, the phase-averaged spectrum is well characterized by a blackbody plus double power-law model, in agreement with previous, multi-mission X-ray results. However, we are unable to reproduce the spectral results reported using Suzaku observations. The pulsed fraction of the source is found to increase with photon energy. The measured rms pulsed fraction is ~12% and ~17% at ~20 keV and ~50 keV, respectively. We detect a new feature in the 24--35 keV band pulse profile that is uniquely double-peaked. This feature may be associated with a possible absorption or emission feature in the phase-resolved spectrum. We fit the X-ray data using the recently developed electron-positron outflow model of Beloborodov (2013) for the hard X-ray emission from magnetars. This produces a satisfactory fit allowing a constraint on the angle between the rotation and magnetic axes of the neutron star of ~20 degrees and on the angle between the rotation axis and line-of-sight of ~50 degrees. In this model, the soft X-ray component is inconsistent with a single blackbody; adding a second blackbody or a power-law component fits the data. The two-blackbody interpretation suggests a hot spot of temperature kT~0.9 keV occupying ~1% of the stellar surface.

preprint2013arXiv

The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

preprint2013arXiv

The Swift/BAT Hard X-ray Transient Monitor

The Swift/Burst Alert Telescope (BAT) hard X-ray transient monitor provides near real-time coverage of the X-ray sky in the energy range 15-50 keV. The BAT observes 88% of the sky each day with a detection sensitivity of 5.3 mCrab for a full-day observation and a time resolution as fine as 64 seconds. The three main purposes of the monitor are (1) the discovery of new transient X-ray sources, (2) the detection of outbursts or other changes in the flux of known X-ray sources, and (3) the generation of light curves of more than 900 sources spanning over eight years. The primary interface for the BAT transient monitor is a public web page. Between 2005 February 12 and 2013 April 30, 245 sources have been detected in the monitor, 146 of them persistent and 99 detected only in outburst. Among these sources, 17 were previously unknown and were discovered in the transient monitor. In this paper, we discuss the methodology and the data processing and filtering for the BAT transient monitor and review its sensitivity and exposure. We provide a summary of the source detections and classify them according to the variability of their light curves. Finally, we review all new BAT monitor discoveries; for the new sources that are previously unpublished, we present basic data analysis and interpretations.

preprint2012arXiv

Fast X-ray transients in the Galactic bulge with the Rossi X-ray Timing Explorer

In X-ray binaries, rapid variability in X-ray flux of greater than an order of magnitude on time-scales of a day or less appears to be a signature of wind accretion from a supergiant companion. When the variability takes the form of rare, brief, bright outbursts with only faint emission between them, the systems are called Supergiant Fast X-ray Transients (SFXTs). We present data from twice-weekly scans of the Galactic bulge by the Rossi X-ray Timing Explorer (RXTE) that allow us to compare the behaviour of known SFXTs and possible SFXT candidates with the persistently bright supergiant X-ray binary 4U 1700-377. We independently confirm the orbital periods reported by other groups for SFXTs SAX J1818.6-1703 and IGR J17544-2619. The new data do not independently reproduce the orbital period reported for XTE J1739-302, but slightly improve the significance of the original result when the data are combined. The bulge source XTE J1743-363 shows a combination of fast variability and a long-term decline in activity, the latter behaviour not being characteristic of supergiant X-ray binaries. A far-red spectrum of the companion suggests that it is a symbiotic neutron star binary rather than a high-mass binary, and the reddest known of this class: the spectral type is approximately M8 III.

preprint2011arXiv

Support for temporally varying behavior of the Pioneer anomaly from the extended Pioneer 10 and 11 Doppler data sets

The Pioneer anomaly is a small sunward anomalous acceleration found in the trajectory analysis of the Pioneer 10 and 11 spacecraft. As part of the investigation of the effect, analysis of recently recovered Doppler data for both spacecraft has been completed. The presence of a small anomalous acceleration is confirmed using data spans more than twice as long as those that were previously analyzed. We examine the constancy and direction of the Pioneer anomaly, and conclude that: i) the data favor a temporally decaying anomalous acceleration (~2\times 10^{-11} m/s^2/yr) with an over 10% improvement in the residuals compared to a constant acceleration model; ii) although the direction of the acceleration remains imprecisely determined, we find no support in favor of a Sun-pointing direction over the Earth-pointing or along the spin-axis directions, and iii) support for an early "onset" of the acceleration remains weak in the pre-Saturn Pioneer 11 tracking data. We present these new findings and discuss their implications for the nature of the Pioneer anomaly.

preprint2011arXiv

The 5 hr pulse period and broadband spectrum of the Symbiotic X-ray Binary 3A 1954+319

We present an analysis of the highly variable accreting X-ray pulsar 3A 1954+319 using 2005-2009 monitoring data obtained with INTEGRAL and Swift. This considerably extends the pulse period history and covers flaring episodes in 2005 and 2008. In 2006 the source was identified as one of only a few known symbiotic X-ray binaries (SyXBs), i.e., systems composed of a neutron star accreting from the inhomogeneous medium around an M-giant star. The extremely long pulse period of 5.3 hr is directly visible in the 2008 INTEGRAL-ISGRI outburst light curve. The pulse profile is double peaked and generally not significantly energy dependent although there is an indication of possible softening during the main pulse. During the outburst a strong spin-up of -1.8 10^(-4) hr hr^(-1) occurred. Between 2005 and 2008 a long-term spin-down trend of 2.1 10^-5 hr hr^(-1) was observed for the first time for this source. The 3-80 keV pulse peak spectrum of 3A 1954+319 during the 2008 flare could be well described by a thermal Comptonization model. We interpret the results within the framework of a recently developed quasi-spherical accretion model for SyXBs.

preprint2007arXiv

The long-term evolution of the spin, pulse shape, and orbit of the accretion-powered millisecond pulsar SAX J1808.4-3658

We present a 7 yr timing study of the 2.5 ms X-ray pulsar SAX J1808.4-3658, an X-ray transient with a recurrence time of ~2 yr, using data from the Rossi X-ray Timing Explorer covering 4 transient outbursts (1998-2005). We verify that the 401 Hz pulsation traces the spin frequency fundamental and not a harmonic. Substantial pulse shape variability, both stochastic and systematic, was observed during each outburst. Analysis of the systematic pulse shape changes suggests that, as an outburst dims, the X-ray "hot spot" on the pulsar surface drifts longitudinally and a second hot spot may appear. The overall pulse shape variability limits the ability to measure spin frequency evolution within a given X-ray outburst (and calls previous nudot measurements of this source into question), with typical upper limits of |nudot| < 2.5x10^{-14} Hz/s (2 sigma). However, combining data from all the outbursts shows with high (6 sigma) significance that the pulsar is undergoing long-term spin down at a rate nudot = (-5.6+/-2.0)x10^{-16} Hz/s, with most of the spin evolution occurring during X-ray quiescence. We discuss the possible contributions of magnetic propeller torques, magnetic dipole radiation, and gravitational radiation to the measured spin down, setting an upper limit of B < 1.5x10^8 G for the pulsar's surface dipole magnetic field and and Q/I < 5x10^{-9} for the fractional mass quadrupole moment. We also measured an orbital period derivative of Pdot = (3.5+/-0.2)x10^{-12} s/s. This surprising large Pdot is reminiscent of the large and quasi-cyclic orbital period variation observed in the so-called "black widow" millisecond radio pulsars, supporting speculation that SAX J1808.4-3658 may turn on as a radio pulsar during quiescence. In an appendix we derive an improved (0.15 arcsec) source position from optical data.

preprint2006arXiv

RXTE Observations of 1A 1744-361: Correlated Spectral and Timing Behavior

We analyze Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) data of the transient low mass X-ray binary (LMXB) system 1A 1744-361. We explore the X-ray intensity and spectral evolution of the source, perform timing analysis, and find that 1A 1744-361 shows `atoll' behavior during the outbursts. The color-color diagram indicates that this LMXB was observed in a low intensity spectrally hard (low-hard) state and in a high intensity `banana' state. The low-hard state shows a horizontal pattern in the color-color diagram, and the previously reported `dipper QPO' appears only during this state. We also perform energy spectral analyses, and report the first detection of broad iron emission line and iron absorption edge from 1A 1744-361.