Showing 1–49 of 49 results for author: Archibald, R

Convergence Analysis for A Stochastic Maximum Principle Based Data Driven Feedback Control Algorithm

Authors: Siming Liang, Hui Sun, Richard Archibald, Feng Bao

Abstract: This paper presents convergence analysis of a novel data-driven feedback control algorithm designed for generating online controls based on partial noisy observational data. The algorithm comprises a particle filter-enabled state estimation component, estimating the controlled system's state via indirect observations, alongside an efficient stochastic maximum principle type optimal control solver.… ▽ More This paper presents convergence analysis of a novel data-driven feedback control algorithm designed for generating online controls based on partial noisy observational data. The algorithm comprises a particle filter-enabled state estimation component, estimating the controlled system's state via indirect observations, alongside an efficient stochastic maximum principle type optimal control solver. By integrating weak convergence techniques for the particle filter with convergence analysis for the stochastic maximum principle control solver, we derive a weak convergence result for the optimization procedure in search of optimal data-driven feedback control. Numerical experiments are performed to validate the theoretical findings. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: arXiv admin note: text overlap with arXiv:2404.05734

Sparse $L^1$-Autoencoders for Scientific Data Compression

Authors: Matthias Chung, Rick Archibald, Paul Atzberger, Jack Michael Solomon

Abstract: Scientific datasets present unique challenges for machine learning-driven compression methods, including more stringent requirements on accuracy and mitigation of potential invalidating artifacts. Drawing on results from compressed sensing and rate-distortion theory, we introduce effective data compression methods by developing autoencoders using high dimensional latent spaces that are $L^1$-regul… ▽ More Scientific datasets present unique challenges for machine learning-driven compression methods, including more stringent requirements on accuracy and mitigation of potential invalidating artifacts. Drawing on results from compressed sensing and rate-distortion theory, we introduce effective data compression methods by developing autoencoders using high dimensional latent spaces that are $L^1$-regularized to obtain sparse low dimensional representations. We show how these information-rich latent spaces can be used to mitigate blurring and other artifacts to obtain highly effective data compression methods for scientific data. We demonstrate our methods for short angle scattering (SAS) datasets showing they can achieve compression ratios around two orders of magnitude and in some cases better. Our compression methods show promise for use in addressing current bottlenecks in transmission, storage, and analysis in high-performance distributed computing environments. This is central to processing the large volume of SAS data being generated at shared experimental facilities around the world to support scientific investigations. Our approaches provide general ways for obtaining specialized compression methods for targeted scientific datasets. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 11 pages, 6 figures

An Online Algorithm for Solving Feedback Optimal Control Problems with Partial Observations

Authors: Siming Liang, Ruoyu Hu, Feng Bao, Richard Archibald, Guannan Zhang

Abstract: This paper presents a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation… ▽ More This paper presents a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation, which is often not feasible in real-world applications. Our approach underscores the significance of utilizing observational data to inform control policy design. Specifically, we introduce a kernel learning backward stochastic differential equation (SDE) filter to enhance data assimilation efficiency and propose a sample-wise stochastic optimization method within the stochastic maximum principle framework. Numerical experiments validate the efficacy and accuracy of our algorithm, showcasing its high efficiency in solving feedback optimal control problems with partial observation. △ Less

Submitted 21 March, 2024; originally announced April 2024.

Comments: arXiv admin note: substantial text overlap with arXiv:2201.10600

MGARD: A multigrid framework for high-performance, error-controlled data compression and refactoring

Authors: Qian Gong, Jieyang Chen, Ben Whitney, Xin Liang, Viktor Reshniak, Tania Banerjee, Jaemoon Lee, Anand Rangarajan, Lipeng Wan, Nicolas Vidal, Qing Liu, Ana Gainaru, Norbert Podhorszki, Richard Archibald, Sanjay Ranka, Scott Klasky

Abstract: We describe MGARD, a software providing MultiGrid Adaptive Reduction for floating-point scientific data on structured and unstructured grids. With exceptional data compression capability and precise error control, MGARD addresses a wide range of requirements, including storage reduction, high-performance I/O, and in-situ data analysis. It features a unified application programming interface (API)… ▽ More We describe MGARD, a software providing MultiGrid Adaptive Reduction for floating-point scientific data on structured and unstructured grids. With exceptional data compression capability and precise error control, MGARD addresses a wide range of requirements, including storage reduction, high-performance I/O, and in-situ data analysis. It features a unified application programming interface (API) that seamlessly operates across diverse computing architectures. MGARD has been optimized with highly-tuned GPU kernels and efficient memory and device management mechanisms, ensuring scalable and rapid operations. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 20 pages, 8 figures

Journal ref: SoftwareX, 24(2023), 101590

Streaming Compression of Scientific Data via weak-SINDy

Authors: Benjamin P. Russo, M. Paul Laiu, Richard Archibald

Abstract: In this paper a streaming weak-SINDy algorithm is developed specifically for compressing streaming scientific data. The production of scientific data, either via simulation or experiments, is undergoing an stage of exponential growth, which makes data compression important and often necessary for storing and utilizing large scientific data sets. As opposed to classical ``offline" compression algor… ▽ More In this paper a streaming weak-SINDy algorithm is developed specifically for compressing streaming scientific data. The production of scientific data, either via simulation or experiments, is undergoing an stage of exponential growth, which makes data compression important and often necessary for storing and utilizing large scientific data sets. As opposed to classical ``offline" compression algorithms that perform compression on a readily available data set, streaming compression algorithms compress data ``online" while the data generated from simulation or experiments is still flowing through the system. This feature makes streaming compression algorithms well-suited for scientific data compression, where storing the full data set offline is often infeasible. This work proposes a new streaming compression algorithm, streaming weak-SINDy, which takes advantage of the underlying data characteristics during compression. The streaming weak-SINDy algorithm constructs feature matrices and target vectors in the online stage via a streaming integration method in a memory efficient manner. The feature matrices and target vectors are then used in the offline stage to build a model through a regression process that aims to recover equations that govern the evolution of the data. For compressing high-dimensional streaming data, we adopt a streaming proper orthogonal decomposition (POD) process to reduce the data dimension and then use the streaming weak-SINDy algorithm to compress the temporal data of the POD expansion. We propose modifications to the streaming weak-SINDy algorithm to accommodate the dynamically updated POD basis. By combining the built model from the streaming weak-SINDy algorithm and a small amount of data samples, the full data flow could be reconstructed accurately at a low memory cost, as shown in the numerical tests. △ Less

Submitted 28 August, 2023; originally announced August 2023.

MSC Class: 37M10; 62J07; 65L60; 41A10; 68T99; 68V99

Convergence Analysis for Training Stochastic Neural Networks via Stochastic Gradient Descent

Authors: Richard Archibald, Feng Bao, Yanzhao Cao, Hui Sun

Abstract: In this paper, we carry out numerical analysis to prove convergence of a novel sample-wise back-propagation method for training a class of stochastic neural networks (SNNs). The structure of the SNN is formulated as discretization of a stochastic differential equation (SDE). A stochastic optimal control framework is introduced to model the training procedure, and a sample-wise approximation scheme… ▽ More In this paper, we carry out numerical analysis to prove convergence of a novel sample-wise back-propagation method for training a class of stochastic neural networks (SNNs). The structure of the SNN is formulated as discretization of a stochastic differential equation (SDE). A stochastic optimal control framework is introduced to model the training procedure, and a sample-wise approximation scheme for the adjoint backward SDE is applied to improve the efficiency of the stochastic optimal control solver, which is equivalent to the back-propagation for training the SNN. The convergence analysis is derived with and without convexity assumption for optimization of the SNN parameters. Especially, our analysis indicates that the number of SNN training steps should be proportional to the square of the number of layers in the convex optimization case. Numerical experiments are carried out to validate the analysis results, and the performance of the sample-wise back-propagation method for training SNNs is examined by benchmark machine learning examples. △ Less

Submitted 17 December, 2022; originally announced December 2022.

A Stochastic Maximum Principle Approach for Reinforcement Learning with Parameterized Environment

Authors: Richard Archibald, Feng Bao, Jiongmin Yong

Abstract: In this work, we introduce a stochastic maximum principle (SMP) approach for solving the reinforcement learning problem with the assumption that the unknowns in the environment can be parameterized based on physics knowledge. For the development of numerical algorithms, we shall apply an effective online parameter estimation method as our exploration technique to estimate the environment parameter… ▽ More In this work, we introduce a stochastic maximum principle (SMP) approach for solving the reinforcement learning problem with the assumption that the unknowns in the environment can be parameterized based on physics knowledge. For the development of numerical algorithms, we shall apply an effective online parameter estimation method as our exploration technique to estimate the environment parameter during the training procedure, and the exploitation for the optimal policy will be achieved by an efficient backward action learning method for policy improvement under the SMP framework. Numerical experiments will be presented to demonstrate that our SMP approach for reinforcement learning can produce reliable control policy, and the gradient descent type optimization in the SMP solver requires less training episodes compared with the standard dynamic programming principle based methods. △ Less

Submitted 6 January, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Modular Light Sources for Microscopy and Beyond (ModLight)

Authors: Graham M Gibson, Robert Archibald, Mark Main, Akhil Kallepalli

Abstract: Delivering light to an object is one of the key steps in any imaging exercise. Tools such as LEDs and lasers are available to achieve this. These components are integrated into systems such as microscopy, medical imaging, remote sensing, and so many more. Motivated by the need for affordable and open access alternatives that are globally relevant, we share the designs and build instructions for mo… ▽ More Delivering light to an object is one of the key steps in any imaging exercise. Tools such as LEDs and lasers are available to achieve this. These components are integrated into systems such as microscopy, medical imaging, remote sensing, and so many more. Motivated by the need for affordable and open access alternatives that are globally relevant, we share the designs and build instructions for modular light source devices that use simple, off-the-shelf components. Light emitted by near-infrared, red, green and blue LEDs are combined with a choice of mirrors or X-Cube prisms to deliver collimated beams of light. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 10 pages, 7 figures

2022 Review of Data-Driven Plasma Science

Authors: Rushil Anirudh, Rick Archibald, M. Salman Asif, Markus M. Becker, Sadruddin Benkadda, Peer-Timo Bremer, Rick H. S. Budé, C. S. Chang, Lei Chen, R. M. Churchill, Jonathan Citrin, Jim A Gaffney, Ana Gainaru, Walter Gekelman, Tom Gibbs, Satoshi Hamaguchi, Christian Hill, Kelli Humbird, Sören Jalas, Satoru Kawaguchi, Gon-Ho Kim, Manuel Kirchen, Scott Klasky, John L. Kline, Karl Krushelnick , et al. (38 additional authors not shown)

Abstract: Data science and technology offer transformative tools and methods to science. This review article highlights latest development and progress in the interdisciplinary field of data-driven plasma science (DDPS). A large amount of data and machine learning algorithms go hand in hand. Most plasma data, whether experimental, observational or computational, are generated or collected by machines today.… ▽ More Data science and technology offer transformative tools and methods to science. This review article highlights latest development and progress in the interdisciplinary field of data-driven plasma science (DDPS). A large amount of data and machine learning algorithms go hand in hand. Most plasma data, whether experimental, observational or computational, are generated or collected by machines today. It is now becoming impractical for humans to analyze all the data manually. Therefore, it is imperative to train machines to analyze and interpret (eventually) such data as intelligently as humans but far more efficiently in quantity. Despite the recent impressive progress in applications of data science to plasma science and technology, the emerging field of DDPS is still in its infancy. Fueled by some of the most challenging problems such as fusion energy, plasma processing of materials, and fundamental understanding of the universe through observable plasma phenomena, it is expected that DDPS continues to benefit significantly from the interdisciplinary marriage between plasma science and data science into the foreseeable future. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 112 pages (including 700+ references), 44 figures, submitted to IEEE Transactions on Plasma Science as a part of the IEEE Golden Anniversary Special Issue

Report number: Los Alamos Report number LA-UR-22-24834

Journal ref: IEEE Transactions on Plasma Science 51, 1750 - 1838 (2023)

A PDE-based Adaptive Kernel Method for Solving Optimal Filtering Problems

Authors: Zezhong Zhang, Richard Archibald, Feng Bao

Abstract: In this paper, we introduce an adaptive kernel method for solving the optimal filtering problem. The computational framework that we adopt is the Bayesian filter, in which we recursively generate an optimal estimate for the state of a target stochastic dynamical system based on partial noisy observational data. The mathematical model that we use to formulate the propagation of the state dynamics i… ▽ More In this paper, we introduce an adaptive kernel method for solving the optimal filtering problem. The computational framework that we adopt is the Bayesian filter, in which we recursively generate an optimal estimate for the state of a target stochastic dynamical system based on partial noisy observational data. The mathematical model that we use to formulate the propagation of the state dynamics is the Fokker-Planck equation, and we introduce an operator decomposition method to efficiently solve the Fokker-Planck equation. An adaptive kernel method is introduced to adaptively construct Gaussian kernels to approximate the probability distribution of the target state. Bayesian inference is applied to incorporate the observational data into the state model simulation. Numerical experiments have been carried out to validate the performance of our kernel method. △ Less

Submitted 9 March, 2022; originally announced March 2022.

A Kernel Learning Method for Backward SDE Filter

Authors: Richard Archibald, Feng Bao

Abstract: In this paper, we develop a kernel learning backward SDE filter method to estimate the state of a stochastic dynamical system based on its partial noisy observations. A system of forward backward stochastic differential equations is used to propagate the state of the target dynamical model, and Bayesian inference is applied to incorporate the observational information. To characterize the dynamica… ▽ More In this paper, we develop a kernel learning backward SDE filter method to estimate the state of a stochastic dynamical system based on its partial noisy observations. A system of forward backward stochastic differential equations is used to propagate the state of the target dynamical model, and Bayesian inference is applied to incorporate the observational information. To characterize the dynamical model in the entire state space, we introduce a kernel learning method to learn a continuous global approximation for the conditional probability density function of the target state by using discrete approximated density values as training data. Numerical experiments demonstrate that the kernel learning backward SDE is highly effective and highly efficient. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Network Security Modeling using NetFlow Data: Detecting Botnet attacks in IP Traffic

Authors: Ganesh Subramaniam, Huan Chen, Ravi Varadhan, Robert Archibald

Abstract: Cybersecurity, security monitoring of malicious events in IP traffic, is an important field largely unexplored by statisticians. Computer scientists have made significant contributions in this area using statistical anomaly detection and other supervised learning methods to detect specific malicious events. In this research, we investigate the detection of botnet command and control (C&C) hosts in… ▽ More Cybersecurity, security monitoring of malicious events in IP traffic, is an important field largely unexplored by statisticians. Computer scientists have made significant contributions in this area using statistical anomaly detection and other supervised learning methods to detect specific malicious events. In this research, we investigate the detection of botnet command and control (C&C) hosts in massive IP traffic. We use the NetFlow data, the industry standard for monitoring of IP traffic for exploratory analysis and extracting new features. Using statistical as well as deep learning models, we develop a statistical intrusion detection system (SIDS) to predict traffic traces identified with malicious attacks. Employing interpretative machine learning techniques, botnet traffic signatures are derived. These models successfully detected botnet C&C hosts and compromised devices. The results were validated by matching predictions to existing blacklists of published malicious IP addresses. △ Less

Submitted 19 August, 2021; originally announced August 2021.

A Backward SDE Method for Uncertainty Quantification in Deep Learning

Authors: Richard Archibald, Feng Bao, Yanzhao Cao, He Zhang

Abstract: We develop a probabilistic machine learning method, which formulates a class of stochastic neural networks by a stochastic optimal control problem. An efficient stochastic gradient descent algorithm is introduced under the stochastic maximum principle framework. Numerical experiments for applications of stochastic neural networks are carried out to validate the effectiveness of our methodology. We develop a probabilistic machine learning method, which formulates a class of stochastic neural networks by a stochastic optimal control problem. An efficient stochastic gradient descent algorithm is introduced under the stochastic maximum principle framework. Numerical experiments for applications of stochastic neural networks are carried out to validate the effectiveness of our methodology. △ Less

Submitted 3 April, 2021; v1 submitted 28 November, 2020; originally announced November 2020.

Integrating Deep Learning in Domain Sciences at Exascale

Authors: Rick Archibald, Edmond Chow, Eduardo D'Azevedo, Jack Dongarra, Markus Eisenbach, Rocco Febbo, Florent Lopez, Daniel Nichols, Stanimire Tomov, Kwai Wong, Junqi Yin

Abstract: This paper presents some of the current challenges in designing deep learning artificial intelligence (AI) and integrating it with traditional high-performance computing (HPC) simulations. We evaluate existing packages for their ability to run deep learning models and applications on large-scale HPC systems efficiently, identify challenges, and propose new asynchronous parallelization and optimiza… ▽ More This paper presents some of the current challenges in designing deep learning artificial intelligence (AI) and integrating it with traditional high-performance computing (HPC) simulations. We evaluate existing packages for their ability to run deep learning models and applications on large-scale HPC systems efficiently, identify challenges, and propose new asynchronous parallelization and optimization techniques for current large-scale heterogeneous systems and upcoming exascale systems. These developments, along with existing HPC AI software capabilities, have been integrated into MagmaDNN, an open-source HPC deep learning framework. Many deep learning frameworks are targeted at data scientists and fall short in providing quality integration into existing HPC workflows. This paper discusses the necessities of an HPC deep learning framework and how those needs can be provided (e.g., as in MagmaDNN) through a deep integration with existing HPC libraries, such as MAGMA and its modular memory management, MPI, CuBLAS, CuDNN, MKL, and HIP. Advancements are also illustrated through the use of algorithmic enhancements in reduced- and mixed-precision, as well as asynchronous optimization methods. Finally, we present illustrations and potential solutions for enhancing traditional compute- and data-intensive applications at ORNL and UTK with AI. The approaches and future challenges are illustrated in materials science, imaging, and climate applications. △ Less

Submitted 22 November, 2020; originally announced November 2020.

An efficient numerical algorithm for solving data driven feedback control problems

Authors: Richard Archibald, Feng Bao, Jiongmin Yong, Tao Zhou

Abstract: The goal of this paper is to solve a class of stochastic optimal control problems numerically, in which the state process is governed by an Itô type stochastic differential equation with control process entering both in the drift and the diffusion, and is observed partially. The optimal control of feedback form is determined based on the available observational data. We call this type of control p… ▽ More The goal of this paper is to solve a class of stochastic optimal control problems numerically, in which the state process is governed by an Itô type stochastic differential equation with control process entering both in the drift and the diffusion, and is observed partially. The optimal control of feedback form is determined based on the available observational data. We call this type of control problems the data driven feedback control. The computational framework that we introduce to solve such type of problems aims to find the best estimate for the optimal control as a conditional expectation given the observational information. To make our method feasible in providing timely feedback to the controlled system from data, we develop an efficient stochastic optimization algorithm to implement our computational framework. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Reconstruction of effective potential from statistical analysis of dynamic trajectories

Authors: Ali Yousefzadi Nobakht, Ondrej Dyck, David B. Lingerfelt, Feng Bao, Maxim Ziatdinov, Artem Maksov, Bobby G. Sumpter, Richard Archibald, Stephen Jesse, Sergei V. Kalinin, Kody J. H. Law

Abstract: The broad incorporation of microscopic methods is yielding a wealth of information on atomic and mesoscale dynamics of individual atoms, molecules, and particles on surfaces and in open volumes. Analysis of such data necessitates statistical frameworks to convert observed dynamic behaviors to effective properties of materials. Here we develop a method for stochastic reconstruction of effective act… ▽ More The broad incorporation of microscopic methods is yielding a wealth of information on atomic and mesoscale dynamics of individual atoms, molecules, and particles on surfaces and in open volumes. Analysis of such data necessitates statistical frameworks to convert observed dynamic behaviors to effective properties of materials. Here we develop a method for stochastic reconstruction of effective acting potentials from observed trajectories. Using the Silicon vacancy defect in graphene as a model, we develop a statistical framework to reconstruct the free energy landscape from calculated atomic displacements. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: 12 pages, 5 figures. This manuscript is a part of update to previous work (arXiv:1804.03729v1) authors found some of the analysis in the previous work to be not accurate and this manuscript is a partial update to that work

Two new outbursts and transient hard X-rays from 1E 1048.1$-$5937

Authors: R. F. Archibald, P. Scholz, V. M. Kaspi, S. P. Tendulkar, A. P. Beardmore

Abstract: Since its discovery, 1E 1048.1$-$5937 has been one of the most active magnetars, both in terms of radiative outbursts, and changes to its spin properties. Here we report on a continuing monitoring campaign with the Neil Gehrels Swift Observatory X-ray Telescope in which we observe two new outbursts from this source. The first outburst occurred in 2016 July, and the second in 2017 December, reachin… ▽ More Since its discovery, 1E 1048.1$-$5937 has been one of the most active magnetars, both in terms of radiative outbursts, and changes to its spin properties. Here we report on a continuing monitoring campaign with the Neil Gehrels Swift Observatory X-ray Telescope in which we observe two new outbursts from this source. The first outburst occurred in 2016 July, and the second in 2017 December, reaching peak 0.5-10 keV absorbed fluxes of $3.2^{+0.2}_{-0.3}\times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ and $2.2^{+0.2}_{-0.2}\times10^{-11}$ erg s$^{-1}$ cm$^{-2}$, respectively, factors of $\sim$5 and $\sim 4$ above the quiescent flux. Both new outbursts were accompanied by spin-up glitches with amplitudes of $Δν= 4.47(6)\times10^{-7}$ Hz and $Δν= 4.32(5)\times10^{-7}$ Hz, respectively. Following the 2016 July outburst, we observe, as for past outbursts, a period of delayed torque fluctuations, which reach a peak spin-down of $1.73\pm0.01$ times the quiescent rate, and which dominates the spin evolution compared to the spin-up glitches. We also report an observation near the peak of the first of these outbursts with NuSTAR in which hard X-ray emission is detected from the source. This emission is well characterized by an absorbed blackbody plus a broken power law, with a power-law index above $13.4\pm0.6$ keV of $0.5_{-0.2}^{+0.3}$, similar to those observed in both persistent and transient magnetars. The hard X-ray results are broadly consistent with models of electron/positron cooling in twisted magnetic field bundles in the outer magnetosphere. However the repeated outbursts and associated torque fluctuations in this source remain puzzling. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 11 pages, 5 figures, accepted in ApJ

Super-resolution energy spectra from neutron direct-geometry spectrometers

Authors: Fahima Islam, Jiao Y. Y. Lin, Richard Archibald, Douglas L. Abernathy, Iyad Al-Qasir, Anne A. Campbell, Matthew B. Stone, Garrett E. Granroth

Abstract: Neutron direct-geometry time-of-flight chopper spectroscopy is instrumental in studying fundamental excitations of vibrational and/or magnetic origin. We report here that techniques in super-resolution optical imagery (which is in real-space) can be adapted to enhance resolution and reduce noise for a neutron spectroscopy (an instrument for mapping excitations in reciprocal space). The procedure t… ▽ More Neutron direct-geometry time-of-flight chopper spectroscopy is instrumental in studying fundamental excitations of vibrational and/or magnetic origin. We report here that techniques in super-resolution optical imagery (which is in real-space) can be adapted to enhance resolution and reduce noise for a neutron spectroscopy (an instrument for mapping excitations in reciprocal space). The procedure to reconstruct super-resolution energy spectra of phonon density of states relies on a realization of multi-frame registration, accurate determination of the energy-dependent point spread function, asymmetric nature of instrument resolution broadening, and iterative reconstructions. Applying these methods to phonon density of states data for a graphite sample demonstrates contrast enhancement, noise reduction, and ~5-fold improvement over nominal energy resolution. The data were collected at three different incident energies measured at the Wide Angular-Range Chopper Spectrometer at the Spallation Neutron Source. △ Less

Submitted 22 June, 2019; originally announced June 2019.

X-ray timing studies of the low-field magnetar CXOU J164710.2-455216

Authors: Hongjun An, Robert Archibald

Abstract: We report results of X-ray timing analyses for the low-field magnetar CXOU~J164710.2$-$455216 which exhibited multiple outbursts. We use data taken with NICER, NuSTAR, Chandra, and Neil-Gehrels-Swift telescopes between 2017 and 2018 when the source was in an active state. We perform semi-phase-coherent timing analyses to measure the spin parameters and a spin-inferred magnetic-field strength (… ▽ More We report results of X-ray timing analyses for the low-field magnetar CXOU~J164710.2$-$455216 which exhibited multiple outbursts. We use data taken with NICER, NuSTAR, Chandra, and Neil-Gehrels-Swift telescopes between 2017 and 2018 when the source was in an active state. We perform semi-phase-coherent timing analyses to measure the spin parameters and a spin-inferred magnetic-field strength ($B_s$) of the magnetar. Using a semi-phase-coherent method, we infer the magnetic field strengths to be $3-4\times 10^{13}\rm \ G$ at the observation period ($\sim$MJD 58000), and by comparing with previous frequency measurements (MJD 54000) a long-term average value of $B_s$ is estimated to be $\approx4\times 10^{13}\rm \ G$. So this analysis may add CXOU~J164710.2$-$455216 to the ranks of low-field magnetars. The inferred characteristic age ($τ_c$) is 1--2 Myr which is smaller than the age of Westerlund~1, so the magnetar's association with the star cluster is still secure. For the low dipole field and the large age, recent multiple outbursts observed from the source are hard to explain unless it has strong magnetic multipole components. We also find timing anomalies around outburst epochs, which suggests that there may be spin-down torque applied to the magnetar near the epochs as was proposed in magnetar models. △ Less

Submitted 7 May, 2019; originally announced May 2019.

Comments: 7 pages, 2 figures, accepted for publication in ApJL

Edge-adaptive l2 regularization image reconstruction from non-uniform Fourier data

Authors: Victor Churchill, Rick Archibald, Anne Gelb

Abstract: Total variation regularization based on the l1 norm is ubiquitous in image reconstruction. However, the resulting reconstructions are not always as sparse in the edge domain as desired. Iteratively reweighted methods provide some improvement in accuracy, but at the cost of extended runtime. In this paper we examine these methods for the case of data acquired as non-uniform Fourier samples. We then… ▽ More Total variation regularization based on the l1 norm is ubiquitous in image reconstruction. However, the resulting reconstructions are not always as sparse in the edge domain as desired. Iteratively reweighted methods provide some improvement in accuracy, but at the cost of extended runtime. In this paper we examine these methods for the case of data acquired as non-uniform Fourier samples. We then develop a non-iterative weighted regularization method that uses a pre-processing edge detection to find exactly where the sparsity should be in the edge domain. We show that its performance in terms of both accuracy and speed has the potential to outperform reweighted TV regularization methods. △ Less

Submitted 20 November, 2018; originally announced November 2018.

MSC Class: 68U10; 65F22

Journal ref: Inverse Problems & Imaging, 13(5), 931-958 (2019)

The 2016 outburst of PSR J1119-6127: cooling & a spin-down dominated glitch

Authors: R. F. Archibald, V. M. Kaspi, S. P. Tendulkar, P. Scholz

Abstract: We report on the aftermath of a magnetar outburst from the young, high-magnetic-field radio pulsar PSR J1119-6127 that occurred on 2016 July 27. We present the results of a monitoring campaign using the Neil Gehrels Swift X-ray Telescope, NuSTAR, and XMM-Newton. After reaching a peak luminosity of ~300 times the quiescent luminosity, the pulsar's X-ray flux declined by factor of ~50 on a time scal… ▽ More We report on the aftermath of a magnetar outburst from the young, high-magnetic-field radio pulsar PSR J1119-6127 that occurred on 2016 July 27. We present the results of a monitoring campaign using the Neil Gehrels Swift X-ray Telescope, NuSTAR, and XMM-Newton. After reaching a peak luminosity of ~300 times the quiescent luminosity, the pulsar's X-ray flux declined by factor of ~50 on a time scale of several months. The X-ray spectra are well described by a blackbody and a hard power-law tail. After an initial rapid decline during the first day of the outburst, we observe the blackbody temperature rising from kT = 0.9 keV to 1.05 keV during the first two weeks of the outburst, before cooling to 0.9 keV. During this time, the blackbody radius decreases monotonically by a factor of ~4 over a span of nearly 200 days. We also report a heretofore unseen highly pulsed hard X-ray emission component, which fades on a similar timescale to the soft X-ray flux, as predicted by models of relaxation of magnetospheric current twists. The previously reported spin-up glitch which accompanied this outburst was followed by a period of enhanced and erratic torque, leading to a net spin-down of $\sim3.5\times10^{-4}$ Hz, a factor of ~24 over-recovery. We suggest that this and other radiatively loud magnetar-type glitch recoveries are dominated by magnetospheric processes, in contrast to conventional radio pulsar glitch recoveries which are dominated by internal physics. △ Less

Submitted 4 June, 2018; originally announced June 2018.

Comments: Submitted to ApJ

Backward SDE Filter for Jump Diffusion Processes and Its Applications in Material Sciences

Authors: Richard Archibald, Feng Bao, Peter Maksymovych

Abstract: The connection between forward backward doubly stochastic differential equations and the optimal filtering problem is established without using the Zakai's equation. The solutions of forward backward doubly stochastic differential equations are expressed in terms of conditional law of a partially observed Markov diffusion process. It then follows that the adjoint time-inverse forward backward doub… ▽ More The connection between forward backward doubly stochastic differential equations and the optimal filtering problem is established without using the Zakai's equation. The solutions of forward backward doubly stochastic differential equations are expressed in terms of conditional law of a partially observed Markov diffusion process. It then follows that the adjoint time-inverse forward backward doubly stochastic differential equations governs the evolution of the unnormalized filtering density in the optimal filtering problem. △ Less

Submitted 28 May, 2018; originally announced May 2018.

Single atom force measurements: mapping potential energy landscapes via electron beam induced single atom dynamics

Authors: Ondrej Dyck, Feng Bao, Maxim Ziatdinov, Ali Yousefzadi Nobakht, Seungha Shin, Kody Law, Artem Maksov, Bobby G. Sumpter, Richard Archibald, Stephen Jesse, Sergei V. Kalinin

Abstract: In the last decade, the atomically focused beam of a scanning transmission electron microscope (STEM) was shown to induce a broad set of transformations of material structure, open pathways for probing atomic-scale reactions and atom-by-atom matter assembly. However, the mechanisms of beam-induced transformations remain largely unknown, due to an extreme mismatch between the energy and time scales… ▽ More In the last decade, the atomically focused beam of a scanning transmission electron microscope (STEM) was shown to induce a broad set of transformations of material structure, open pathways for probing atomic-scale reactions and atom-by-atom matter assembly. However, the mechanisms of beam-induced transformations remain largely unknown, due to an extreme mismatch between the energy and time scales of electron passage through solids and atomic and molecular motion. Here, we demonstrate that a single dopant Si atom in the graphene lattice can be used as an atomic scale force sensor, providing information on the random force exerted by the beam on chemically-relevant time scales. Using stochastic reconstruction of molecular dynamic simulations, we recover the potential energy landscape of the atom and use it to determine the beam-induced effects in the thermal (i.e. white noise) approximation. We further demonstrate that the moving atom under beam excitation can be used to map potential energy along step edges, providing information about atomic-scale potentials in solids. These studies open the pathway for quantitative studies of beam-induced atomic dynamics, elementary mechanisms of solid-state transformations, and predictive atom-by-atom fabrication. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Revival of the magnetar PSR J1622-4950: observations with MeerKAT, Parkes, XMM-Newton, Swift, Chandra, and NuSTAR

Authors: F. Camilo, P. Scholz, M. Serylak, S. Buchner, M. Merryfield, V. M. Kaspi, R. F. Archibald, M. Bailes, A. Jameson, W. van Straten, J. Sarkissian, J. E. Reynolds, S. Johnston, G. Hobbs, T. D. Abbott, R. M. Adam, G. B. Adams, T. Alberts, R. Andreas, K. M. B. Asad, D. E. Baker, T. Baloyi, E. F. Bauermeister, T. Baxana, T. G. H. Bennett , et al. (183 additional authors not shown)

Abstract: New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x la… ▽ More New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation. △ Less

Submitted 5 April, 2018; originally announced April 2018.

Comments: Published in ApJ (2018 April 5); 13 pages, 4 figures

Journal ref: ApJ 856 (2018) 180

Magnetar-like X-ray bursts suppress pulsar radio emission

Authors: R. F. Archibald, M. Burgay, M. Lyutikov, V. M. Kaspi, P. Esposito, G. Israel, M. Kerr, A. Possenti, N. Rea, J. Sarkissian, P. Scholz, S. P. Tendulkar

Abstract: Rotation-powered pulsars and magnetars are two different observational manifestations of neutron stars: rotation powered pulsars are rapidly spinning objects that are mostly observed as pulsating radio sources, while magnetars, neutron stars with the highest known magnetic fields, often emit short-duration X-ray bursts. Here we report simultaneous observations of the high-magnetic-field radio puls… ▽ More Rotation-powered pulsars and magnetars are two different observational manifestations of neutron stars: rotation powered pulsars are rapidly spinning objects that are mostly observed as pulsating radio sources, while magnetars, neutron stars with the highest known magnetic fields, often emit short-duration X-ray bursts. Here we report simultaneous observations of the high-magnetic-field radio pulsar PSR J1119-6127 at X-ray, with XMM-Newton & NuSTAR, and at radio energies with Parkes radio telescope, during a period of magnetar-like bursts. The rotationally powered radio emission shuts off coincident with the occurrence of multiple X-ray bursts, and recovers on a time scale of ~70 seconds. These observations of related radio and X-ray phenomena further solidify the connection between radio pulsars and magnetars, and suggest that the pair plasma produced in bursts can disrupt the acceleration mechanism of radio emitting particles. △ Less

Submitted 10 October, 2017; originally announced October 2017.

Comments: 8 pages, 4 figures. Submitted to ApJL

The glitches and rotational history of the highly energetic young pulsar PSR J0537$-$6910

Authors: Robert D. Ferdman, Robert F. Archibald, Konstantinos N. Gourgouliatos, Victoria M. Kaspi

Abstract: We present a timing and glitch analysis of the young X-ray pulsar PSR J0537$-$6910, located within the Large Magellanic Cloud, using 13 years of data from the now decommissioned Rossi X-ray Timing Explorer. Rotating with a spin period of 16 ms, PSR J0537$-$6910 is the fastest spinning and most energetic young pulsar known. It also displays the highest glitch activity of any known pulsar. We have f… ▽ More We present a timing and glitch analysis of the young X-ray pulsar PSR J0537$-$6910, located within the Large Magellanic Cloud, using 13 years of data from the now decommissioned Rossi X-ray Timing Explorer. Rotating with a spin period of 16 ms, PSR J0537$-$6910 is the fastest spinning and most energetic young pulsar known. It also displays the highest glitch activity of any known pulsar. We have found 42 glitches over the data span, corresponding to a glitch rate of 3.2 yr$^{-1}$, with an overall glitch activity rate of $8.8\times 10^{-7}\,$yr$^{-1}$. The high glitch frequency has allowed us to study the glitch behavior in ways that are inaccessible in other pulsars. We observe a strong linear correlation between spin frequency glitch magnitude and wait time to the following glitch. We also find that the post-glitch spin-down recovery is well described by a single two-component model fit to all glitches for which we have adequate input data. This consists of an exponential amplitude $A = (7.6 \pm 1.0)\times 10^{-14}\,$s$^{-2}$ and decay timescale $τ= 27_{-6}^{+7}\,$d, and linear slope $m = (4.1\pm 0.4)\times 10^{-16}\,$s$^{-2}\,$d$^{-1}$. The latter slope corresponds to a second frequency derivative $\ddotν = (4.7\pm 0.5) \times 10^{-22}\,$s$^{-3}$, from which we find an implied braking index $n=7.4 \pm 0.8$. We also present a maximum-likelihood technique for searching for periods in event-time data, which we used to both confirm previously published values and determine rotation frequencies in later observations. We discuss the implied constraints on glitch models from the observed behavior of this system, which we argue cannot be fully explained in the context of existing theories. △ Less

Submitted 11 December, 2017; v1 submitted 29 August, 2017; originally announced August 2017.

Comments: 17 pages, 2 tables, 8 figures. Accepted for publication in the Astrophysical Journal

Accelerating optimization-based computed tomography via sparse matrix approximations

Authors: Richard C. Barnard, Rick Archibald

Abstract: Variational formulations of reconstruction in computed tomography have the notable drawback of requiring repeated evaluations of both the forward Radon transform and either its adjoint or an approximate inverse transform which are relatively expensive. We look at two methods for reducing the effect of this resulting computational bottleneck via approximating the transform evaluation with sparse ma… ▽ More Variational formulations of reconstruction in computed tomography have the notable drawback of requiring repeated evaluations of both the forward Radon transform and either its adjoint or an approximate inverse transform which are relatively expensive. We look at two methods for reducing the effect of this resulting computational bottleneck via approximating the transform evaluation with sparse matrix multiplications. The first method is applicable for general iterative optimization algorithms. The second is applicable in error-forgetting algorithms such as split Bregman. We demonstrate these approximations significantly reduce the needed computational time needed for the iterative algorithms needed to solve the reconstruction problem while still providing good reconstructions. △ Less

Submitted 21 May, 2017; originally announced May 2017.

MSC Class: 44A12; 49N45

Recovery guarantees for compressed sensing with unknown errors

Authors: Simone Brugiapaglia, Ben Adcock, Richard K. Archibald

Abstract: From a numerical analysis perspective, assessing the robustness of l1-minimization is a fundamental issue in compressed sensing and sparse regularization. Yet, the recovery guarantees available in the literature usually depend on a priori estimates of the noise, which can be very hard to obtain in practice, especially when the noise term also includes unknown discrepancies between the finite model… ▽ More From a numerical analysis perspective, assessing the robustness of l1-minimization is a fundamental issue in compressed sensing and sparse regularization. Yet, the recovery guarantees available in the literature usually depend on a priori estimates of the noise, which can be very hard to obtain in practice, especially when the noise term also includes unknown discrepancies between the finite model and data. In this work, we study the performance of l1-minimization when these estimates are not available, providing robust recovery guarantees for quadratically constrained basis pursuit and random sampling in bounded orthonormal systems. Several applications of this work are approximation of high-dimensional functions, infinite-dimensional sparse regularization for inverse problems, and fast algorithms for non-Cartesian Magnetic Resonance Imaging. △ Less

Submitted 9 May, 2017; v1 submitted 14 February, 2017; originally announced February 2017.

Comments: Final manuscript for "Sampling Theory and Applications" 12th International Conference (SampTA 2017)

Swift observations of two outbursts from the magnetar 4U 0142+61

Authors: R. F. Archibald, V. M. Kaspi, P. Scholz, A. P. Beardmore, N. Gehrels, J. A. Kennea

Abstract: 4U 0142+61 is one of a small class of persistently bright magnetars. Here we report on a monitoring campaign of 4U 0142+61 from 2011 July 26 - 2016 June 12 using the Swift X-ray Telescope, continuing a 16 year timing campaign with the Rossi X-ray Timing Explorer. We show that 4U 0142+61 had two radiatively loud timing events, on 2011 July 29 and 2015 February 28, both with short soft gamma-ray bur… ▽ More 4U 0142+61 is one of a small class of persistently bright magnetars. Here we report on a monitoring campaign of 4U 0142+61 from 2011 July 26 - 2016 June 12 using the Swift X-ray Telescope, continuing a 16 year timing campaign with the Rossi X-ray Timing Explorer. We show that 4U 0142+61 had two radiatively loud timing events, on 2011 July 29 and 2015 February 28, both with short soft gamma-ray bursts, and a long-lived flux decay associated with each case. We show that the 2015 timing event resulted in a net spin-down of the pulsar due to over-recovery of a glitch. We compare this timing event to previous such events in other high-magnetic-field pulsars, and discuss net spin-down glitches now seen in several young, high-B pulsars. △ Less

Submitted 29 November, 2016; originally announced November 2016.

Comments: 9 pages, 5 figures, Submitted to ApJ

Near Infrared Counterpart of 2E 1613.5-5053, the Central Source in Supernova Remnant RCW 103

Authors: Shriharsh P. Tendulkar, Victoria M. Kaspi, Robert F. Archibald, Paul Scholz

Abstract: On 2016 June 22, 2E 1613.5-5053, the puzzling central compact object in supernova remnant RCW 103, emitted a magnetar-like burst. Using Director's Discretionary Time, we observed 2E 1613.5-5053 with the Hubble Space Telescope (WFC3/IR) and we report here on the detection of a previously unseen infrared counterpart. In observations taken on 2016 July 4 and August 11, we detect a new source (… ▽ More On 2016 June 22, 2E 1613.5-5053, the puzzling central compact object in supernova remnant RCW 103, emitted a magnetar-like burst. Using Director's Discretionary Time, we observed 2E 1613.5-5053 with the Hubble Space Telescope (WFC3/IR) and we report here on the detection of a previously unseen infrared counterpart. In observations taken on 2016 July 4 and August 11, we detect a new source ($m_\mathrm{F110W} = 26.3$ AB mag and $m_\mathrm{F160W} = 24.2$ AB mag) at the Chandra position of 2E 1613.5-5053 which was not detected in HST/NICMOS images from 2002 August 15 and October 8 to a depth of 24.5 AB mag (F110W) and 25.5 AB mag (F160W). We show that these deep IR observations rule out the possibility of an accreting binary but mimic IR emission properties of magnetars and isolated neutron stars. The presence or absence of a low-mass fallback disk cannot be confirmed from our observations. △ Less

Submitted 7 October, 2016; originally announced October 2016.

Comments: 7 pages, 3 figures. Submitted to ApJ

A Magnetar-like Outburst from a High-B Radio Pulsar

Authors: R. F. Archibald, V. M. Kaspi, S. P. Tendulkar, P. Scholz

Abstract: Radio pulsars are believed to have their emission powered by the loss of rotational kinetic energy. By contrast, magnetars show intense X-ray and gamma-ray radiation whose luminosity greatly exceeds that due to spin-down and is believed to be powered by intense internal magnetic fields. A basic prediction of this picture is that radio pulsars of high magnetic field should show magnetar-like emissi… ▽ More Radio pulsars are believed to have their emission powered by the loss of rotational kinetic energy. By contrast, magnetars show intense X-ray and gamma-ray radiation whose luminosity greatly exceeds that due to spin-down and is believed to be powered by intense internal magnetic fields. A basic prediction of this picture is that radio pulsars of high magnetic field should show magnetar-like emission. Here we report on a magnetar-like X-ray outburst from the radio pulsar PSR J1119-6127, heralded by two short bright X-ray bursts on 2016 July 27 and 28 (Kennea et al. 2016; Younes et al. 2016). Using Target-of-Opportunity data from the Swift X-ray Telescope and NuSTAR, we show that this pulsar's flux has brightened by a factor of > 160 in the 0.5-10 keV band, and its previously soft X-ray spectrum has undergone a strong hardening, with strong pulsations appearing for the first time above 2.5 keV, with phase-averaged emission detectable up to 25 keV. By comparing Swift-XRT and NuSTAR timing data with a pre-outburst ephemeris derived from Fermi Large Area Telescope data, we find that the source has contemporaneously undergone a large spin-up glitch of amplitude df/f = 5.74(8) E-6. The collection of phenomena observed thus far in this outburst strongly mirrors those in most magnetar outbursts and provides an unambiguous connection between the radio pulsar and magnetar populations. △ Less

Submitted 31 August, 2016; v1 submitted 2 August, 2016; originally announced August 2016.

Comments: 7 pages, 3 figures, accepted to ApJL

A High Braking Index for a Pulsar

Authors: R. F. Archibald, E. V. Gotthelf, R. D. Ferdman, V. M. Kaspi, S. Guillot, F. A. Harrison, E. F. Keane, M. J. Pivovaroff, D. Stern, S. P. Tendulkar, J. A. Tomsick

Abstract: We present a phase-coherent timing solution for PSR J1640-4631, a young 206 ms pulsar using X-ray timing observations taken with NuSTAR. Over this timing campaign, we have measured the braking index of PSR J1640-4631 to be n = 3.15+/-0.03. Using a series of simulations, we argue that this unusually high braking index is not due to timing noise, but is intrinsic to the pulsar's spin-down. We cannot… ▽ More We present a phase-coherent timing solution for PSR J1640-4631, a young 206 ms pulsar using X-ray timing observations taken with NuSTAR. Over this timing campaign, we have measured the braking index of PSR J1640-4631 to be n = 3.15+/-0.03. Using a series of simulations, we argue that this unusually high braking index is not due to timing noise, but is intrinsic to the pulsar's spin-down. We cannot, however, rule out contamination due to an unseen glitch recovery, although the recovery timescale would have to be longer than most yet observed. If this braking index is eventually proven to be stable, it demonstrates that pulsar braking indices greater than 3 are allowed in nature, hence other physical mechanisms such as mass or magnetic quadrupoles are important in pulsar spin-down. We also present a 3-sigma upper limit on the pulsed flux at 1.4 GHz of 0.018 mJy. △ Less

Submitted 1 March, 2016; originally announced March 2016.

Comments: 4 pages, 2 figures. Accepted at ApJL

The NuSTAR view of the non-thermal emission from PSR J0437-4715

Authors: Sebastien Guillot, Victoria M. Kaspi, Robert F. Archibald, Matteo Bachetti, Chris Flynn, Fabian Jankowski, Matthew Bailes, Steven Boggs, Finn E. Christensen, William W. Craig, Charles A. Hailey, Fiona A. Harrison, Daniel Stern, Willam W. Zhang

Abstract: We present a hard X-ray NuSTAR observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period ~5.76 ms are detected at energies up to 20 keV. We measure a photon index $Γ= 1.65\pm0.24$ (90% confidence) for the power law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton sp… ▽ More We present a hard X-ray NuSTAR observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period ~5.76 ms are detected at energies up to 20 keV. We measure a photon index $Γ= 1.65\pm0.24$ (90% confidence) for the power law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase-resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsar's phase-folded lightcurve with the pulsar's well-defined mass and distance from radio timing observations. △ Less

Submitted 14 September, 2016; v1 submitted 12 December, 2015; originally announced December 2015.

Comments: Accepted for publication in MNRAS. 12 pages, 8 figures, 3 table

NuSTAR Observations of Magnetar 1E 1048.1-5937

Authors: C. Yang, R. F. Archibald, J. K. Vogel, H. An, V. M. Kaspi, S. Guillot, A. M. Beloborodov, M. Pivovaroff

Abstract: We report on simultaneous NuSTAR and XMM-Newton observations of the magnetar 1E 1048.1$-$5937 obtained in July 2013, along with Rossi X-ray Timing Explorer (RXTE) data for the same source obtained from December 2002 through March 2012. The NuSTAR data provide a clear detection of this magnetar's persistent emission up to 20 keV. We detect a previously unreported small secondary peak in the average… ▽ More We report on simultaneous NuSTAR and XMM-Newton observations of the magnetar 1E 1048.1$-$5937 obtained in July 2013, along with Rossi X-ray Timing Explorer (RXTE) data for the same source obtained from December 2002 through March 2012. The NuSTAR data provide a clear detection of this magnetar's persistent emission up to 20 keV. We detect a previously unreported small secondary peak in the average pulse profile in the 7-10 keV band, which grows to an amplitude comparable to that of the main peak in the 10--20 keV band. We show using RXTE data that this secondary peak is likely transient. We find that the pulsed fraction increases with energy from a value of $\sim$0.55 at $\sim$2 keV to a value of $\sim$0.75 near 8 keV but shows evidence for decreasing at higher energies. After filtering out multiple bright X-ray bursts during the observation, we find that the phase-averaged spectrum from combined NuSTAR and XMM data is well described by an absorbed double blackbody plus power-law model, with no evidence for the spectral turn-up near $\sim$10 keV as has been seen in some other magnetars. Our data allow us to rule out a spectral turn-up similar to those seen in magnetars 4U 0142+61 and 1E 2259+586 of $ΔΓ>= 2$, where $ΔΓ$ is the difference between the soft-band and hard-band photon indexes. The absence of a significant spectral turn-up is consistent with what has been observed from a particularly active subset of magnetars having high spin-inferred magnetic fields, and with previously reported trends suggesting the degree of spectral turn-up is correlated with spin-down rate and/or spin-inferred magnetic field strength. △ Less

Submitted 31 August, 2016; v1 submitted 27 September, 2015; originally announced September 2015.

Comments: 14 pages, 8 figures. Accepted to APJ

On the braking index of the unusual high-B rotation-powered pulsar PSR J1846-0258

Authors: R. F. Archibald, V. M. Kaspi, A. P. Beardmore, N. Gehrels, J. A. Kennea

Abstract: PSR J1846-0258 is an object which straddles the boundary between magnetars and rotation powered pulsars. Though behaving for many years as a rotation-powered pulsar, in 2006, it exhibited distinctly magnetar-like behavior - emitting several short hard X-ray bursts, and a flux increase. Here we report on 7 years of post-outburst timing observations of PSR J1846-0258 using the Rossi X-ray Timing Exp… ▽ More PSR J1846-0258 is an object which straddles the boundary between magnetars and rotation powered pulsars. Though behaving for many years as a rotation-powered pulsar, in 2006, it exhibited distinctly magnetar-like behavior - emitting several short hard X-ray bursts, and a flux increase. Here we report on 7 years of post-outburst timing observations of PSR J1846-0258 using the Rossi X-ray Timing Explorer and the Swift X-ray Telescope. We measure the braking index over the post-magnetar outburst period to be $n=2.19\pm0.03$. This represents a change of $Δn=-0.46\pm0.03$ or a 14.5$\;σ$ difference from the pre-outburst braking index of $n=2.65\pm0.01$, which itself was measured over a span of 6.5 yr. So large and long-lived a change to a pulsar braking index is unprecedented and poses a significant challenge to models of pulsar spin-down. △ Less

Submitted 14 August, 2015; v1 submitted 19 June, 2015; originally announced June 2015.

Comments: 7 pages, 5 figures. APJ, accepted

Long-term timing and emission behavior of the young Crab-like pulsar PSR B0540-69

Authors: Robert D. Ferdman, Robert F. Archibald, Victoria M. Kaspi

Abstract: We present timing solutions and spin properties of the young pulsar PSR B0540-69 from analysis of 15.8 yr of data from the Rossi X-Ray Timing Explorer. We perform a partially phase-coherent timing analysis in order to mitigate the pronounced effects of timing noise in this pulsar. We also perform fully coherent timing over large subsets of the data set in order to arrive at a more precise solution… ▽ More We present timing solutions and spin properties of the young pulsar PSR B0540-69 from analysis of 15.8 yr of data from the Rossi X-Ray Timing Explorer. We perform a partially phase-coherent timing analysis in order to mitigate the pronounced effects of timing noise in this pulsar. We also perform fully coherent timing over large subsets of the data set in order to arrive at a more precise solution. In addition to the previously reported first glitch undergone by this pulsar, we find a second glitch, which occurred at MJD 52927 $\pm$ 4, with fractional changes in spin frequency $Δν/ν= (1.64 \pm 0.05) \times 10^{-9}$ and spin-down rate $Δ\dotν/\dotν = (0.930 \pm 0.011) \times 10^{-4}$ (taken from our fully coherent analysis). We measure a braking index that is consistent over the entire data span, with a mean value $n = 2.129 \pm 0.012$, from our partially coherent timing analysis. We also investigated the emission behavior of this pulsar, and have found no evidence for significant flux changes, flares, burst-type activity, or pulse profile shape variations. While there is strong evidence for the much-touted similarity of PSR B0540-69 to the Crab pulsar, they nevertheless differ in several aspects, including glitch activity, where PSR B0540-69 can be said to resemble certain other very young pulsars. It seems clear that the specific processes governing the formation, evolution, and interiors of this population of recently born neutron stars can vary significantly, as reflected in their observed properties. △ Less

Submitted 2 October, 2015; v1 submitted 30 May, 2015; originally announced June 2015.

Comments: 10 pages, 1 table, 5 figures. Accepted for publication in ApJ

Deep NuSTAR and Swift Monitoring Observations of the Magnetar 1E 1841-045

Authors: Hongjun An, Robert F. Archibald, Romain Hascoet, Victoria M. Kaspi, Andrei M. Beloborodov, Anne M. Archibald, Andy Beardmore, Steven E. Boggs, Finn E. Christensen, William W. Craig, Niel Gehrels, Charles J. Hailey, Fiona A. Harrison, Jamie Kennea, Chryssa Kouveliotou, Daniel Stern, George Younes, William W. Zhang

Abstract: We report on a 350-ks NuSTAR observation of the magnetar 1E 1841-045 taken in 2013 September. During the observation, NuSTAR detected six bursts of short duration, with $T_{90}<1$ s. An elevated level of emission tail is detected after the brightest burst, persisting for $\sim$1 ks. The emission showed a power-law decay with a temporal index of 0.5 before returning to the persistent emission level… ▽ More We report on a 350-ks NuSTAR observation of the magnetar 1E 1841-045 taken in 2013 September. During the observation, NuSTAR detected six bursts of short duration, with $T_{90}<1$ s. An elevated level of emission tail is detected after the brightest burst, persisting for $\sim$1 ks. The emission showed a power-law decay with a temporal index of 0.5 before returning to the persistent emission level. The long observation also provided detailed phase-resolved spectra of the persistent X-ray emission of the source. By comparing the persistent spectrum with that previously reported, we find that the source hard-band emission has been stable over approximately 10 years. The persistent hard X-ray emission is well fitted by a coronal outflow model, where $e^{+/-}$ pairs in the magnetosphere upscatter thermal X-rays. Our fit of phase-resolved spectra allowed us to estimate the angle between the rotational and magnetic dipole axes of the magnetar, $α_{mag}=0.25$, the twisted magnetic flux, $2.5\times10^{26}\rm \ G\ cm^2$, and the power released in the twisted magnetosphere, $L_j=6\times10^{36}\rm \ erg\ s^{-1}$. Assuming this model for the hard X-ray spectrum, the soft X-ray component is well fit by a two-blackbody model, with the hotter blackbody consistent with the footprint of the twisted magnetic field lines on the star. We also report on the 3-year Swift monitoring observations obtained since 2011 July. The soft X-ray spectrum remained stable during this period, and the timing behavior was noisy, with large timing residuals. △ Less

Submitted 19 May, 2015; v1 submitted 13 May, 2015; originally announced May 2015.

Comments: 18 pages, 12 figures. Accepted for publication in ApJ

Repeated, Delayed Torque Variations Following X-ray Flux Enhancements in the Magnetar 1E$\;$1048.1$-$5937

Authors: R. F. Archibald, V. M. Kaspi, C. -Y. Ng, P. Scholz, A. P. Beardmore, N. Gehrels, J. A. Kennea

Abstract: We report on two years of flux and spin evolution monitoring of 1E 1048.1$-$5937, a 6.5-s X-ray pulsar identified as a magnetar. Using {\it Swift} XRT data, we observed an X-ray outburst consisting of an increase in the persistent 1--10 keV flux by a factor of 6.3$\pm$0.2, beginning on 2011 December 31 (MJD 55926). Following a delay of $\sim100$ days, the magnetar entered a period of large torque… ▽ More We report on two years of flux and spin evolution monitoring of 1E 1048.1$-$5937, a 6.5-s X-ray pulsar identified as a magnetar. Using {\it Swift} XRT data, we observed an X-ray outburst consisting of an increase in the persistent 1--10 keV flux by a factor of 6.3$\pm$0.2, beginning on 2011 December 31 (MJD 55926). Following a delay of $\sim100$ days, the magnetar entered a period of large torque variability, with $\dotν$ reaching a factor of $4.55\pm0.05$ times the nominal value, before decaying in an oscillatory manner over a time scale of months. We show by comparing to previous outbursts from the source that this pattern of behavior may repeat itself with a quasi-period of $\sim1800$ days. We compare this phenomenology to periodic torque variations in radio pulsars, finding some similarities which suggest a magnetospheric origin for the behavior of 1E 1048.1$-$5937. △ Less

Submitted 8 December, 2014; originally announced December 2014.

Comments: Submitted to ApJ. 10 pages, 7 figures

Green Bank Telescope and Swift X-ray Telescope Observations of the Galactic Center Radio Magnetar SGR J1745-2900

Authors: Ryan S. Lynch, Robert F. Archibald, Victoria M. Kaspi, Paul Scholz

Abstract: We present results from eight months of Green Bank Telescope 8.7-GHz observations and nearly 18 months Swift X-ray telescope observations of the radio magnetar SGR J1745-2900, which is located 2.4" from Sgr A*. We tracked the magnetar's radio flux density, polarization properties, pulse profile evolution, rotation, and single-pulse behavior. We identified two main periods of activity in SGR J1745-… ▽ More We present results from eight months of Green Bank Telescope 8.7-GHz observations and nearly 18 months Swift X-ray telescope observations of the radio magnetar SGR J1745-2900, which is located 2.4" from Sgr A*. We tracked the magnetar's radio flux density, polarization properties, pulse profile evolution, rotation, and single-pulse behavior. We identified two main periods of activity in SGR J1745-2900. The first is characterized by approximately 5.5 months of relatively stable evolution in radio flux density, rotation, and profile shape, while in the second these properties varied substantially. Specifically, a third profile component emerged and the radio flux increased on average, but also became more variable. Bright single pulses are visible and are well described by a log-normal energy distribution at low to moderate energies, but with an excess at high energies. The 2-10 keV flux has decayed steadily since the initial X-ray outburst, in contrast with the radio flux. Our timing analysis includes Green Bank Telescope, Swift, and NuSTAR data. When we include the X-ray data in our analyses, we find that SGR J1745-2900 exhibited a level of timing noise unprecedented in a radio magnetar, though an analysis of only the radio data indicates timing noise at a level similar to that observed in other radio magnetars. We conclude that, while SGR J1745-2900 is similar to other radio magnetars in many regards, it differs by having experienced a period of relative stability in the radio band that now appears to have ended, while the X-ray properties have evolved independently. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: Submitted to ApJ; 10 pages (emulateapj style), 6 figures, 1 table

NuSTAR Observations of the Magnetar 1E 2259+586

Authors: Julia K. Vogel, Romain Hascoet, Victoria M. Kaspi, Hongjun An, Robert Archibald, Andrei M. Beloborodov, Steven E. Boggs, Finn E. Christensen, William W. Craig, Eric V. Gotthelf, Brian W. Grefenstette, Charles J. Hailey, Fiona A. Harrison, Jamie A. Kennea, Kristin K. Madsen, Michael J. Pivovaroff, Daniel Stern, William W. Zhang

Abstract: We report on new broad band spectral and temporal observations of the magnetar 1E 2259+586, which is located in the supernova remnant CTB 109. Our data were obtained simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift, and cover the energy range from 0.5-79 keV. We present pulse profiles in various energy bands and compare them to previous RXTE results. The NuSTAR data… ▽ More We report on new broad band spectral and temporal observations of the magnetar 1E 2259+586, which is located in the supernova remnant CTB 109. Our data were obtained simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift, and cover the energy range from 0.5-79 keV. We present pulse profiles in various energy bands and compare them to previous RXTE results. The NuSTAR data show pulsations above 20 keV for the first time and we report evidence that one of the pulses in the double-peaked pulse profile shifts position with energy. The pulsed fraction of the magnetar is shown to increase strongly with energy. Our spectral analysis reveals that the soft X-ray spectrum is well characterized by an absorbed double-blackbody or blackbody plus power-law model in agreement with previous reports. Our new hard X-ray data, however, suggests that an additional component, such as a power-law, is needed to describe the NuSTAR and Swift spectrum. We also fit the data with the recently developed coronal outflow model by Beloborodov for hard X-ray emission from magnetars. The outflow from a ring on the magnetar surface is statistically preferred over outflow from a polar cap. △ Less

Submitted 4 August, 2014; originally announced August 2014.

Comments: 37 pages, 9 figures, corresponding author, vogel9@llnl.gov

Journal ref: The Astrophysical Journal, 789, 75 (11pp), 2014 July 1

NuSTAR observations of X-ray bursts from the magnetar 1E 1048.1-5937

Authors: Hongjun An, Victoria M. Kaspi, Andrei M. Beloborodov, Chryssa Kouveliotou, Robert F. Archibald, Steven E. Boggs, Finn E. Christensen, William W. Craig, Eric V. Gotthelf, Brian W. Grefenstette, Charles J. Hailey, Fiona A. Harrison, Kristin K. Madsen, Kaya Mori, Daniel Stern, William W. Zhang

Abstract: We report the detection of eight bright X-ray bursts from the 6.5-s magnetar 1E 1048.1-5937, during a 2013 July observation campaign with the Nuclear Spectroscopic Telescope Array (NuSTAR). We study the morphological and spectral properties of these bursts and their evolution with time. The bursts resulted in count rate increases by orders of magnitude, sometimes limited by the detector dead time,… ▽ More We report the detection of eight bright X-ray bursts from the 6.5-s magnetar 1E 1048.1-5937, during a 2013 July observation campaign with the Nuclear Spectroscopic Telescope Array (NuSTAR). We study the morphological and spectral properties of these bursts and their evolution with time. The bursts resulted in count rate increases by orders of magnitude, sometimes limited by the detector dead time, and showed blackbody spectra with kT=6-8 keV in the T90 duration of 1-4 s, similar to earlier bursts detected from the source. We find that the spectra during the tail of the bursts can be modeled with an absorbed blackbody with temperature decreasing with flux. The bursts flux decays followed a power-law of index 0.8-0.9. In the burst tail spectra, we detect a ~13 keV emission feature, similar to those reported in previous bursts from this source as well as from other magnetars observed with the Rossi X-ray Timing Explorer (RXTE). We explore possible origins of the spectral feature such as proton cyclotron emission, which implies a magnetic field strength of B~2X10^15 G in the emission region. However, the consistency of the energy of the feature in different objects requires further explanation. △ Less

Submitted 12 June, 2014; originally announced June 2014.

Comments: 10 pages, 6 figures, accepted for publication in ApJ

Timing and Flux Evolution of the Galactic Center Magnetar SGR J1745-2900

Authors: V. M. Kaspi, R. F. Archibald, V. Bhalerao, F. Dufour, E. V. Gotthelf, H. An, M. Bachetti, A. M. Beloborodov, S. E. Boggs, F. E. Christensen, W. W. Craig, B. W. Grefenstette, C. J. Hailey, F. A. Harrison, J. A. Kennea, C. Kouveliotou, K. K. Madsen, K. Mori, C. B. Markwardt, D. Stern, J. K. Vogel, W. W. Zhang

Abstract: We present the X-ray timing and spectral evolution of the Galactic Center magnetar SGR J1745-2900 for the first ~4 months post-discovery using data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR)} and Swift observatories. Our timing analysis reveals a large increase in the magnetar spin-down rate by a factor of 2.60 +/- 0.07 over our data span. We further show that the change in s… ▽ More We present the X-ray timing and spectral evolution of the Galactic Center magnetar SGR J1745-2900 for the first ~4 months post-discovery using data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR)} and Swift observatories. Our timing analysis reveals a large increase in the magnetar spin-down rate by a factor of 2.60 +/- 0.07 over our data span. We further show that the change in spin evolution was likely coincident with a bright X-ray burst observed in 2013 June by Swift, and if so, there was no accompanying discontinuity in the frequency. We find that the source 3-10 keV flux has declined monotonically by a factor of ~2 over an 80-day period post-outburst accompanied by a ~20% decrease in the source's blackbody temperature, although there is evidence for both flux and kT having levelled off. We argue that the torque variations are likely to be magnetospheric in nature and will dominate over any dynamical signatures of orbital motion around Sgr A*. △ Less

Submitted 20 March, 2014; originally announced March 2014.

Comments: 26 pre-print pages, 6 figures, Accepted for publication in ApJ

NuSTAR results and future plans for magnetar and rotation-powered pulsar observations

Authors: Hongjun An, Victoria M. Kaspi, Robert Archibald, Matteo Bachetti, Varun Bhalerao, Eric C. Bellm, Andrei M. Beloborodov, Steven E. Boggs, Deepto Chakrabarty, Finn E. Christensen, William W. Craig, Francois Dufour, Karl Forster, Eric V. Gotthelf, Brian W. Grefenstette, Charles J. Hailey, Fiona A. Harrison, Romain Hascoet, Takao Kitaguchi, Chryssa Kouveliotou, Kristin K. Madsen, Kaya Mori, Michael J. Pivovaroff, Vikram R. Rana, Daniel Stern , et al. (5 additional authors not shown)

Abstract: The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray mission in orbit and operates in the 3-79 keV range. NuSTAR's sensitivity is roughly two orders of magnitude better than previous missions in this energy band thanks to its superb angular resolution. Since its launch in 2012 June, NuSTAR has performed excellently and observed many interesting sources including four… ▽ More The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray mission in orbit and operates in the 3-79 keV range. NuSTAR's sensitivity is roughly two orders of magnitude better than previous missions in this energy band thanks to its superb angular resolution. Since its launch in 2012 June, NuSTAR has performed excellently and observed many interesting sources including four magnetars, two rotation-powered pulsars and the cataclysmic variable AE Aquarii. NuSTAR also discovered 3.76-s pulsations from the transient source SGR J1745-29 recently found by Swift very close to the Galactic Center, clearly identifying the source as a transient magnetar. For magnetar 1E 1841-045, we show that the spectrum is well fit by an absorbed blackbody plus broken power-law model with a hard power-law photon index of ~1.3. This is consistent with previous results by INTEGRAL and RXTE. We also find an interesting double-peaked pulse profile in the 25-35 keV band. For AE Aquarii, we show that the spectrum can be described by a multi-temperature thermal model or a thermal plus non-thermal model; a multi-temperature thermal model without a non-thermal component cannot be ruled out. Furthermore, we do not see a spiky pulse profile in the hard X-ray band, as previously reported based on Suzaku observations. For other magnetars and rotation-powered pulsars observed with NuSTAR, data analysis results will be soon available. △ Less

Submitted 5 February, 2014; originally announced February 2014.

Comments: 5 pages, 1 figure. Accepted for publication in Astronomische Nachrichten (AN), Proceedings of the XMM-Newton 2013 Science Workshop "The Fast and the Furious: Energetic Phenomena in Isolated Neutron Stars, Pulsar Wind Nebulae and Supernova Remnants"

On the X-ray variability of magnetar 1RXS J170849.0-400910

Authors: P. Scholz, R. F. Archibald, V. M. Kaspi, C. -Y. Ng, A. P. Beardmore, N. Gehrels, J. A. Kennea

Abstract: We present a long-term X-ray flux and spectral analysis for 1RXS J170849.0-400910 using Swift/XRT spanning over 8 years from 2005-2013. We also analyze two observations from Chandra and XMM in the period from 2003-2004. In this 10-yr period, 1RXS J170849.0-400910 displayed several rotational glitches. Previous studies have claimed variations in the X-ray emission associated with some of the glitch… ▽ More We present a long-term X-ray flux and spectral analysis for 1RXS J170849.0-400910 using Swift/XRT spanning over 8 years from 2005-2013. We also analyze two observations from Chandra and XMM in the period from 2003-2004. In this 10-yr period, 1RXS J170849.0-400910 displayed several rotational glitches. Previous studies have claimed variations in the X-ray emission associated with some of the glitches. From our analysis we find no evidence for significant X-ray flux variations and evidence for only low-level spectral variations. We also present an updated timing solution for 1RXS J170849.0-400910, from RXTE and Swift observations, which includes a previously unreported glitch at MJD 56019. We discuss the frequency and implications of radiatively quiet glitches in magnetars. △ Less

Submitted 20 January, 2014; originally announced January 2014.

Comments: 9 pages, 2 figures, accepted for publication in ApJ

NuSTAR observations of magnetar 1E 1841-045

Authors: Hongjun An, Romain Hascoet, Victoria M. Kaspi, Andrei M. Beloborodov, Francois Dufour, Eric V. Gotthelf, Robert Archibald, Matteo Bachetti, Steven E. Boggs, Finn E. Christensen, William W. Craig, Brian W. Grefenstette, Charles J. Hailey, Fiona A. Harrison, Takao Kitaguchi, Chryssa Kouveliotou, Kristin K. Madsen, Craig B. Markwardt, Daniel Stern, Julia K. Vogel, William W. Zhang

Abstract: 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.… ▽ More 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. △ Less

Submitted 19 November, 2013; v1 submitted 23 October, 2013; originally announced October 2013.

Comments: 9 pages, 4 figures and 4 tables, accepted for publication in ApJ

An Anti-Glitch in a Magnetar

Authors: R. F. Archibald, V. M. Kaspi, C. -Y. Ng, K. N. Gourgouliatos, D. Tsang, P. Scholz, A. P. Beardmore, N. Gehrels, J. A. Kennea

Abstract: Magnetars are neutron stars showing dramatic X-ray and soft $γ$-ray outbursting behaviour that is thought to be powered by intense internal magnetic fields. Like conventional young neutron stars in the form of radio pulsars, magnetars exhibit "glitches" during which angular momentum is believed to be transferred between the solid outer crust and the superfluid component of the inner crust. Hithert… ▽ More Magnetars are neutron stars showing dramatic X-ray and soft $γ$-ray outbursting behaviour that is thought to be powered by intense internal magnetic fields. Like conventional young neutron stars in the form of radio pulsars, magnetars exhibit "glitches" during which angular momentum is believed to be transferred between the solid outer crust and the superfluid component of the inner crust. Hitherto, the several hundred observed glitches in radio pulsars and magnetars have involved a sudden spin-up of the star, due presumably to the interior superfluid rotating faster than the crust. Here we report on X-ray timing observations of the magnetar 1E 2259+586 which we show exhibited a clear "anti-glitch" -- a sudden spin down. We show that this event, like some previous magnetar spin-up glitches, was accompanied by multiple X-ray radiative changes and a significant spin-down rate change. This event, if of origin internal to the star, is unpredicted in models of neutron star spin-down and is suggestive of differential rotation in the neutron star, further supporting the need for a rethinking of glitch theory for all neutron stars. △ Less

Submitted 29 May, 2013; originally announced May 2013.

Journal ref: Nature, 497, 2013, p591-593

Spectral and timing properties of the magnetar CXOU~J164710.2-455216

Authors: Hongjun An, Victoria M. Kaspi, Robert Archibald, Andrew Cumming

Abstract: We report on spectral and timing properties of the magnetar CXOU J164710.2-455216 in the massive star cluster Westerlund 1. Using 11 archival observations obtained with Chandra and XMM-Newton over approximately 1000 days after the source's 2006 outburst, we study the flux and spectral evolution of the source. We show that the hardness of the source, as quantified by hardness ratio, blackbody tempe… ▽ More We report on spectral and timing properties of the magnetar CXOU J164710.2-455216 in the massive star cluster Westerlund 1. Using 11 archival observations obtained with Chandra and XMM-Newton over approximately 1000 days after the source's 2006 outburst, we study the flux and spectral evolution of the source. We show that the hardness of the source, as quantified by hardness ratio, blackbody temperature or power-law photon index, shows a clear correlation with the 2--10 keV absorption-corrected flux and that the power-law component flux decayed faster than the blackbody component for the first ~100 days. We also measure the timing properties of the source by analyzing data spanning approximately 2500 days. The measured period and period derivative are 10.610644(17) s (MJD 53999.06) and <4 X 10^{-13} s s^-1 (90% confidence) which imply that the spin-inferred dipolar magnetic field of the source is less than 7 X 10^{13} G. This is significantly smaller than was suggested previously. We find evidence for a second flux increase, suggesting a second outburst between MJDs 55068 and 55832. Finally, based on a crustal cooling model, we find that the source's cooling curve can be reproduced if we assume that the energy was deposited in the outer crust and that the temperature profile of the star right after the 2006 outburst was relatively independent of density. △ Less

Submitted 1 December, 2012; originally announced December 2012.

Comments: 9 pages, 5 figures, 2 tables, Accepted for publication in ApJ

The New Magnetar Swift J1822.3-1606

Authors: P. Scholz, C. -Y. Ng, M. A. Livingstone, V. M. Kaspi, A. Cumming, R. Archibald

Abstract: On 2011 July 14, a transient X-ray source, Swift J1822.3-1606, was detected by Swift BAT via its burst activities. It was subsequently identified as a new magnetar upon the detection of a pulse period of 8.4 s. Using follow-up RXTE, Swift, and Chandra observations, we have determined a spin-down rate of $\dot{P}\sim3\times10^{-13}$, implying a dipole magnetic field of $\sim5\times10^{13}$\,G, seco… ▽ More On 2011 July 14, a transient X-ray source, Swift J1822.3-1606, was detected by Swift BAT via its burst activities. It was subsequently identified as a new magnetar upon the detection of a pulse period of 8.4 s. Using follow-up RXTE, Swift, and Chandra observations, we have determined a spin-down rate of $\dot{P}\sim3\times10^{-13}$, implying a dipole magnetic field of $\sim5\times10^{13}$\,G, second lowest among known magnetars, although our timing solution is contaminated by timing noise. The post-outburst flux evolution is well modelled by surface cooling resulting from heat injection in the outer crust, although we cannot rule out other models. We measure an absorption column density similar to that of the open cluster M17 at 10\arcmin\ away, arguing for a comparable distance of $\sim$1.6\,kpc. If confirmed, this could be the nearest known magnetar. △ Less

Submitted 29 October, 2012; v1 submitted 26 October, 2012; originally announced October 2012.

Comments: Proceedings of IAUS 291 "Neutron Stars and Pulsars: Challenges and Opportunities after 80 years", J. van Leeuwen (ed.); 3 pages, 2 figures

Post-outburst X-ray flux and timing evolution of Swift J1822.3-1606

Authors: P. Scholz, C. -Y. Ng, M. A. Livingstone, V. M. Kaspi, A. Cumming, R. Archibald

Abstract: Swift J1822.3-1606 was discovered on 2011 July 14 by the Swift Burst Alert Telescope following the detection of several bursts. The source was found to have a period of 8.4377 s and was identified as a magnetar. Here we present a phase-connected timing analysis and the evolution of the flux and spectral properties using RXTE, Swift, and Chandra observations. We measure a spin frequency of 0.118515… ▽ More Swift J1822.3-1606 was discovered on 2011 July 14 by the Swift Burst Alert Telescope following the detection of several bursts. The source was found to have a period of 8.4377 s and was identified as a magnetar. Here we present a phase-connected timing analysis and the evolution of the flux and spectral properties using RXTE, Swift, and Chandra observations. We measure a spin frequency of 0.1185154343(8) s$^{-1}$ and a frequency derivative of $-4.3\pm0.3\times10^{-15}$ at MJD 55761.0, in a timing analysis that include significant non-zero second and third frequency derivatives that we attribute to timing noise. This corresponds to an estimated spin-down inferred dipole magnetic field of $B\sim5\times10^{13}$ G, consistent with previous estimates though still possibly affected by unmodelled noise. We find that the post-outburst 1--10 keV flux evolution can be characterized by a double-exponential decay with decay timescales of $15.5\pm0.5$ and $177\pm14$ days. We also fit the light curve with a crustal cooling model which suggests that the cooling results from heat injection into the outer crust. We find that the hardness-flux correlation observed in magnetar outbursts also characterizes the outburst of Swift J1822.3-1606. We compare the properties of Swift J1822.3-1606 with those of other magnetars and their outbursts. △ Less

Submitted 26 October, 2012; v1 submitted 4 April, 2012; originally announced April 2012.

Comments: 12 pages, 8 figures, accepted for publication in ApJ