During a pointed 2018 NuSTAR observation, we detected a flare with a 2.2 hour duration from the magnetar 1RXS J170849.0$-$400910. The flare, which rose in $\sim25$ seconds to a maximum flux 6 times larger than the persistent emission, is highly pulsed with an rms pulsed fraction of $53\%$. The pulse profile shape consists of two peaks separated by half a rotational cycle, with a peak flux ratio of $\sim$2. The flare spectrum is thermal with an average temperature of 2.1 keV. Phase resolved spectroscopy show that the two peaks possess the same temperature, but differ in size. These observational results along with simple light curve modeling indicate that two identical antipodal spots, likely the magnetic poles, are heated simultaneously at the onset of the flare and for its full duration. Hence, the origin of the flare has to be connected to the global dipolar structure of the magnetar. This might best be achieved externally, via twists to closed magnetospheric dipolar field lines seeding bombardment of polar footpoint locales with energetic pairs. Approximately 1.86 hours following the onset of the flare, a short burst with its own 3-minute thermal tail occurred. The burst tail is also pulsating at the spin period of the source and phase-aligned with the flare profile, implying an intimate connection between the two phenomena. The burst may have been caused by a magnetic reconnection event in the same twisted dipolar field lines anchored to the surface hot spots, with subsequent return currents supplying extra heat to these polar caps.

MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $\Gamma{}$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods.

This paper describes the design and implementation of Stingray, a library in Python built to perform time series analysis and related tasks on astronomical light curves. Its core functionality comprises a range of Fourier analysis techniques commonly used in spectral-timing analysis, as well as extensions for analyzing pulsar data, simulating data sets, and statistical modeling. Its modular build allows for easy extensions and incorporation of its methods into data analysis workflows and pipelines. We aim for the library to be a platform for the implementation of future spectral-timing techniques. Here, we describe the overall vision and framework, core functionality, extensions, and connections to high-level command-line and graphical interfaces. The code is well-tested, with a test coverage of currently 95%, and is accompanied by extensive API documentation and a set of step-by-step tutorials.

The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

Selecting a cohort from a set of candidates is a common task within and beyond academia. Admitting students, awarding grants, choosing speakers for a conference are situations where human biases may affect the make-up of the final cohort. We propose a new algorithm, Entrofy, designed to be part of a larger decision making strategy aimed at making cohort selection as just, quantitative, transparent, and accountable as possible. We suggest this algorithm be embedded in a two-step selection procedure. First, all application materials are stripped of markers of identity that could induce conscious or sub-conscious bias. During blind review, the committee selects all applicants, submissions, or other entities that meet their merit-based criteria. This often yields a cohort larger than the admissible number. In the second stage, the target cohort can be chosen from this meritorious pool via a new algorithm and software tool. Entrofy optimizes differences across an assignable set of categories selected by the human committee. Criteria could include gender, academic discipline, experience with certain technologies, or other quantifiable characteristics. The Entrofy algorithm yields the computational maximization of diversity by solving the tie-breaking problem with provable performance guarantees. We show how Entrofy selects cohorts according to pre-determined characteristics in simulated sets of applications and demonstrate its use in a case study. This cohort selection process allows human judgment to prevail when assessing merit, but assigns the assessment of diversity to a computational process less likely to be beset by human bias. Importantly, the stage at which diversity assessments occur is fully transparent and auditable with Entrofy. Splitting merit and diversity considerations into their own assessment stages makes it easier to explain why a given candidate was selected or rejected.

We present an active-learning strategy for undergraduates that applies Bayesian analysis to candy-covered chocolate $\text{m&m's}^\circledR$. The exercise is best suited for small class sizes and tutorial settings, after students have been introduced to the concepts of Bayesian statistics. The exercise takes advantage of the non-uniform distribution of $\text{m&m's}^\circledR~$ colours, and the difference in distributions made at two different factories. In this paper, we provide the intended learning outcomes, lesson plan and step-by-step guide for instruction, and open-source teaching materials. We also suggest an extension to the exercise for the graduate-level, which incorporates hierarchical Bayesian analysis.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

We report the analysis of simultaneous XMM-Newton+NuSTAR observations of two low-luminosity Active Galactic Nuclei (LLAGN), NGC 3998 and NGC 4579. We do not detect any significant variability in either source over the ~3 day length of the NuSTAR observations. The broad-band 0.5-60 keV spectrum of NGC 3998 is best fit with a cutoff power-law, while the one for NGC 4579 is best fit with a combination of a hot thermal plasma model, a power-law, and a blend of Gaussians to fit an Fe complex observed between 6 and 7 keV. Our main spectral results are the following: (1) neither source shows any reflection hump with a $3\sigma$ reflection fraction upper-limits $R<0.3$ and $R<0.18$ for NGC 3998 and NGC 4579, respectively; (2) the 6-7 keV line complex in NGC 4579 could either be fit with a narrow Fe K line at 6.4 keV and a moderately broad Fe XXV line, or 3 relatively narrow lines, which includes contribution from Fe XXVI; (3) NGC 4579 flux is 60% brighter than previously detected with XMM-Newton, accompanied by a hardening in the spectrum; (4) we measure a cutoff energy $E_{\rm cut}=107_{-18}^{+27}$ keV in NGC 3998, which represents the lowest and best constrained high-energy cutoff ever measured for an LLAGN; (5) NGC 3998 spectrum is consistent with a Comptonization model with either a sphere ($\tau\approx3\pm1$) or slab ($\tau\approx1.2\pm0.6$) geometry, corresponding to plasma temperatures between 20 and 150 keV. We discuss these results in the context of hard X-ray emission from bright AGN, other LLAGN, and hot accretion flow models.

Relativistic effects dominate the emission of blazar jets complicating our understanding of their intrinsic properties. Although many methods have been proposed to account for them, the variability Doppler factor method has been shown to describe the blazar populations best. We use a Bayesian hierarchical code called {\it Magnetron} to model the light curves of 1029 sources observed by the Owens Valley Radio Observatory's 40-m telescope as a series of flares with an exponential rise and decay, and estimate their variability brightness temperature. Our analysis allows us to place the most stringent constraints on the equipartition brightness temperature i.e., the maximum achieved intrinsic brightness temperature in beamed sources which we found to be $\rm \langle T_{eq}\rangle=2.78\times10^{11}K\pm26\%$. Using our findings we estimated the variability Doppler factor for the largest sample of blazars increasing the number of available estimates in the literature by almost an order of magnitude. Our results clearly show that $\gamma$-ray loud sources have faster and higher amplitude flares than $\gamma$-ray quiet sources. As a consequence they show higher variability brightness temperatures and thus are more relativistically beamed, with all of the above suggesting a strong connection between the radio flaring properties of the jet and $\gamma$-ray emission.

We present a spectral analysis of the lobes and X-ray jets of Cygnus A, using more than 2 Ms of $\textit{Chandra}$ observations. The X-ray jets are misaligned with the radio jets and significantly wider. We detect non-thermal emission components in both lobes and jets. For the eastern lobe and jet, we find 1 keV flux densities of $71_{-10}^{+10}$ nJy and $24_{-4}^{+4}$ nJy, and photon indices of $1.72_{-0.03}^{+0.03}$ and $1.64_{-0.04}^{+0.04}$ respectively. For the western lobe and jet, we find flux densities of $50_{-13}^{+12}$ nJy and $13_{-5}^{+5}$ nJy, and photon indices of $1.97_{-0.10}^{+0.23}$ and $1.86_{-0.12}^{+0.18}$ respectively. Using these results, we modeled the electron energy distributions of the lobes as broken power laws with age breaks. We find that a significant population of non-radiating particles is required to account for the total pressure of the eastern lobe. In the western lobe, no such population is required and the low energy cutoff to the electron distribution there needs to be raised to obtain pressures consistent with observations. This discrepancy is a consequence of the differing X-ray photon indices, which may indicate that the turnover in the inverse-Compton spectrum of the western lobe is at lower energies than in the eastern lobe. We modeled the emission from both jets as inverse-Compton emission. There is a narrow region of parameter space for which the X-ray jet can be a relic of an earlier active phase, although lack of knowledge about the jet's electron distribution and particle content makes the modelling uncertain.

The trigger for the short bursts observed in $\gamma$-rays from many magnetar sources remains unknown. One particular open question in this context is the localization of burst emission to a singular active region or a larger area across the neutron star. While several observational studies have attempted to investigate this question by looking at the phase dependence of burst properties, results have been mixed. At the same time, it is not obvious a priori that bursts from a localized active region would actually give rise to a detectable phase-dependence, taking into account issues such as geometry, relativistic effects, and intrinsic burst properties such brightness and duration. In this paper, we build a simple theoretical model to investigate the circumstances under which the latter effects could affect detectability of a dependence of burst emission on rotational phase. We find that even for strongly phase-dependent emission, inferred burst properties may not show a rotational phase dependence depending on the geometry of the system and the observer. Furthermore, the observed properties of bursts with durations short as 10-20% of the spin period can vary strongly depending on the rotational phase at which the burst was emitted. We also show that detectability of a rotational phase dependence depends strongly on the minimum number of bursts observed, and find that existing burst samples may simply be too small to rule out a phase dependence.

Dead time affects many of the instruments used in X-ray astronomy, by producing a strong distortion in power density spectra. This can make it difficult to model the aperiodic variability of the source or look for quasi-periodic oscillations. Whereas in some instruments a simple a priori correction for dead-time-affected power spectra is possible, this is not the case for others such as NuSTAR, where the dead time is non-constant and long (~2.5 ms). Bachetti et al. 2015 suggested the cospectrum obtained from light curves of independent detectors within the same instrument as a possible way out, but this solution has always only been a partial one: the measured rms was still affected by dead time, because the width of the power distribution of the cospectrum was modulated by dead time in a frequency-dependent way. In this Letter we suggest a new, powerful method to normalize cospectra and, with some caveats, even power density spectra. Our approach uses the difference of the Fourier amplitudes from two independent detectors to characterize and filter out the effect of dead time. This method is crucially important for the accurate modelling of periodograms derived from instruments affected by dead time on board current missions like NuSTAR and ASTROSAT, but also future missions such as IXPE

We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparable to the population of Solar System C/D asteroids, Trojans, or comets. We find no evidence of any cometary activity at a heliocentric distance of 1.46 au, approximately 1.5 months after 1I's closest approach distance to the Sun. Significant brightness variability was seen in the $r$ observations, with the object becoming notably brighter towards the end of the run. By combining our APO photometric time series data with the Discovery Channel Telescope (DCT) data of Knight et al. (2017), taken 20 h later on 2017 October 30, we construct an almost complete light curve with a most probable lightcurve period of $P \simeq 4~{\rm h}$. Our results imply a double peaked rotation period of 8.1 $\pm$ 0.02 h, with a peak-to-peak amplitude of 1.5 - 2.1 mags. Assuming that 1I's shape can be approximated by an ellipsoid, the amplitude constraint implies that 1I has an axial ratio of 3.5 to 10.3, which is strikingly elongated. Assuming that 1I is rotating above its critical break up limit, our results are compatible with 1I having having modest cohesive strength and may have obtained its elongated shape during a tidal disruption event before being ejected from its home system. Astrometry useful for constraining 1I's orbit was also obtained and published in Weaver et al. (2017).

In recent years, the cross spectrum has received considerable attention as a means of characterising the variability of astronomical sources as a function of wavelength. While much has been written about the statistics of time and phase lags, the cospectrum has only recently been understood as means of mitigating instrumental effects dependent on temporal frequency in astronomical detectors, as well as a method of characterizing the coherent variability in two wavelength ranges on different time scales. In this paper, we lay out the statistical foundations of the cospectrum, starting with the simplest case of detecting a periodic signal in the presence of white noise. This case is especially relevant for detecting faint X-ray pulsars in detectors heavily affected by instrumental effects, including NuSTAR, Astrosat and IXPE. We show that the statistical distributions of both single and averaged cospectra differ considerably from those for standard periodograms. While a single cospectrum follows a Laplace distribution exactly, averaged cospectra are approximated by a Gaussian distribution only for more than ~30 averaged segments, dependent on the number of trials. We provide an instructive example of a quasi-periodic oscillation in NuSTAR and show that applying standard periodogram statistics leads to underestimated tail probabilities for period detection. We also demonstrate the application of these distributions to a NuSTAR observation of the X-ray pulsar Hercules X-1.

Across almost all scientific disciplines, the instruments that record our experimental data and the methods required for storage and data analysis are rapidly increasing in complexity. This gives rise to the need for scientific communities to adapt on shorter time scales than traditional university curricula allow for, and therefore requires new modes of knowledge transfer. The universal applicability of data science tools to a broad range of problems has generated new opportunities to foster exchange of ideas and computational workflows across disciplines. In recent years, hack weeks have emerged as an effective tool for fostering these exchanges by providing training in modern data analysis workflows. While there are variations in hack week implementation, all events consist of a common core of three components: tutorials in state-of-the-art methodology, peer-learning and project work in a collaborative environment. In this paper, we present the concept of a hack week in the larger context of scientific meetings and point out similarities and differences to traditional conferences. We motivate the need for such an event and present in detail its strengths and challenges. We find that hack weeks are successful at cultivating collaboration and the exchange of knowledge. Participants self-report that these events help them both in their day-to-day research as well as their careers. Based on our results, we conclude that hack weeks present an effective, easy-to-implement, fairly low-cost tool to positively impact data analysis literacy in academic disciplines, foster collaboration and cultivate best practices.

High energy ($> 250$ keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares. It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper we undertake a detailed examination of this model. We investigate the properties of the required scatterers, and whether the mechanism is consistent with the degree of pulsed emission observed in the tail of the giant flare. We conclude that the mechanism is consistent with current data, although the origin of the scattering population remains an open question. We propose an alternative picture in which the emission is closer to that star and is dominated by synchrotron radiation. The $RHESSI$ observations of the December 2004 flare modestly favor this latter picture. We assess the prospects for the Fermi Gamma-Ray Space Telescope to detect and characterize a similar high energy component in a future giant flare. Such a detection should help to resolve some of the outstanding issues.

We analyzed broad-band X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a BB+PL or 2BB model during all three outbursts. NuSTAR observations revealed a hard X-ray tail, $\Gamma=0.9$, extending up to 79 keV, with flux larger than the one detected <10 keV. Imaging analysis of Chandra data did not reveal small-scale extended emission around the source. Following the outbursts, the total 0.5-10 keV flux from SGR J1935+2154 increased in concordance to its bursting activity, with the flux at activation onset increasing by a factor of $\sim7$ following its strongest June 2016 outburst. A Swift/XRT observation taken 1.5 days prior to the onset of this outburst showed a flux level consistent with quiescence. We show that the flux increase is due to the PL or hot BB component, which increased by a factor of $25$ compared to quiescence, while the cold BB component $kT=0.47$ keV remained more or less constant. The 2014 and 2015 outbursts decayed quasi-exponentially with time-scales of $\sim40$ days, while the stronger May and June 2016 outbursts showed a quick short-term decay with time-scales of $\sim4$ days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 $\mu$Jy for the 4.6 GHz observations and 7 $\mu$Jy for the 1.4 GHz observations. We discuss these results in the framework of the current magnetar theoretical models.

Among the population of known galactic black hole X-ray binaries, GRS 1915+105 stands out in multiple ways. It has been in continuous outburst since 1992, and has shown a wide range of different states that can be distinguished by their timing and spectral properties. These states, also observed in IGR J17091-3624, have in the past been linked to accretion dynamics. Here, we present the first comprehensive study into the long-term evolution of GRS 1915+105, using the entire data set observed with RXTE over its sixteen-year lifetime. We develop a set of descriptive features allowing for automatic separation of states, and show that supervised machine learning in the form of logistic regression and random forests can be used to efficiently classify the entire data set. For the first time, we explore the duty cycle and time evolution of states over the entire sixteen-year time span, and find that the temporal distribution of states has significantly changed over the span of the observations. We connect the machine classification with physical interpretations of the phenomenology in terms of chaotic and stochastic processes.

In June 2015, the black hole X-ray binary (BHXRB) V404 Cygni went into outburst for the first time since 1989. Here, we present a comprehensive search for quasi-periodic oscillations (QPOs) of V404 Cygni during its recent outburst, utilizing data from six instruments on board five different X-ray missions: Swift/XRT, Fermi/GBM, Chandra/ACIS, INTEGRAL's IBIS/ISGRI and JEM-X, and NuSTAR. We report the detection of a QPO at 18 mHz simultaneously with both Fermi/GBM and Swift/XRT, another example of a rare but slowly growing new class of mHz-QPOs in BHXRBs linked to sources with a high orbital inclination. Additionally, we find a duo of QPOs in a Chandra/ACIS observation at 73 mHz and 1.03 Hz, as well as a QPO at 136 mHz in a single Swift/XRT observation that can be interpreted as standard Type-C QPOs. Aside from the detected QPOs, there is significant structure in the broadband power, with a strong feature observable in the Chandra observations between 0.1 and 1 Hz. We discuss our results in the context of current models for QPO formation.

We have compiled the most comprehensive burst sample from magnetar 4U 0142+61, comprising 27 bursts from its three burst-active episodes in 2011, 2012 and the latest one in 2015 observed with Swift/BAT and Fermi/GBM. Bursts from 4U 0142+61 morphologically resemble typical short bursts from other magnetars. However, 4U 0142+61 bursts are less energetic compared to the bulk of magnetar bursts. We uncovered an extended tail emission following a burst on 2015 February 28, with a thermal nature, cooling over a time-scale of several minutes. During this tail emission, we also uncovered pulse peak phase aligned X-ray bursts , which could originate from the same underlying mechanism as that of the extended burst tail, or an associated and spatially coincident but different mechanism.

Two energetic hard X-ray bursts have recently triggered the Fermi and Swift space observatories from the rotation powered pulsar, PSR J1119-6127. We have performed in depth spectral and temporal analyses of these two events. Our extensive searches in both observatory data for lower luminosity bursts uncovered 10 additional events from the source. We report here on the timing and energetics of the 12 bursts from PSR J1119-6127 during its burst active phase of 2016 July 26 and 28. We also found a spectral softer X-ray flux enhancement in a post burst episode, which shows evidence of cooling. We discuss here the implications of these results on the nature of this unusual high-field radio pulsar, which firmly place it within the typical magnetar population.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

There have recently been several reports of apparently periodic variations in the light curves of quasars, e.g. PG 1302-102 by Graham et al. (2015a). Any quasar showing periodic oscillations in brightness would be a strong candidate to be a close binary supermassive black hole and, in turn, a candidate for gravitational wave studies. However, normal quasars -- powered by accretion onto a single, supermassive black hole -- usually show stochastic variability over a wide range of timescales. It is therefore important to carefully assess the methods for identifying periodic candidates from among a population dominated by stochastic variability. Using a Bayesian analysis of the light curve of PG 1302-102, we find that a simple stochastic process is preferred over a sinusoidal variations. We then discuss some of the problems one encounters when searching for rare, strictly periodic signals among a large number of irregularly sampled, stochastic time series, and use simulations of quasar light curves to illustrate these points. From a few thousand simulations of steep spectrum (`red noise') stochastic processes, we find many simulations that display few-cycle periodicity like that seen in PG 1302-102. We emphasise the importance of calibrating the false positive rate when the number of targets in a search is very large.

We report on the analysis of two deep XMM-Newton observations of the magnetar Swift J1834.9-0846 and its surrounding extended emission taken in March 2014 and October 2014, 2.5 and 3.1 years after the source went into outburst. The magnetar is only weakly detected in the first observation with an absorption corrected flux $F_{\rm 0.5-10 keV}\approx4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and a $3\sigma$ upper limit during the second observation of about $3\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. This flux level is more than 3 orders of magnitude lower than the flux measured at the outburst onset on September 2011. The extended emission, centered at the magnetar position and elongated towards the south-west, is clearly seen in both observations; it is best fit by a highly absorbed power-law (PL), with a hydrogen column density of $N_{\rm H}=8.0\times10^{22}$ cm$^{-2}$ and PL photon index $\Gamma=2.2\pm0.2$. Its flux is constant between the two observations at $F_{\rm 0.5-10 keV}=1.3\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$. We find no statistically significant changes in the spectral shape or the flux of this extended emission over a period of 9 years from 2005 to 2014. These new results strongly support the extended emission nature as a wind nebula and firmly establish Swift J1834.9-0846 as the first magnetar to show a surrounding wind nebula. Further, our results imply that such nebulae are no longer exclusive to rotation-powered pulsars and narrow the gap between these two sub-populations of isolated neutron stars. The size and spectrum of the nebula are compatible with those of pulsar-wind nebulae but its radiative efficiency $\eta_{\rm X}=L_{\rm X}/\dot{E}\approx0.1$ is markedly high, possibly pointing to an additional wind component in Swift J1834.9-0846.

Neutron stars are a prime laboratory for testing physical processes under conditions of strong gravity, high density, and extreme magnetic fields. Among the zoo of neutron star phenomena, magnetars stand out for their bursting behaviour, ranging from extremely bright, rare giant flares to numerous, less energetic recurrent bursts. The exact trigger and emission mechanisms for these bursts are not known; favoured models involve either a crust fracture and subsequent energy release into the magnetosphere, or explosive reconnection of magnetic field lines. In the absence of a predictive model, understanding the physical processes responsible for magnetar burst variability is difficult. Here, we develop an empirical model that decomposes magnetar bursts into a superposition of small spike-like features with a simple functional form, where the number of model components is itself part of the inference problem. The cascades of spikes that we model might be formed by avalanches of reconnection, or crust rupture aftershocks. Using Markov Chain Monte Carlo (MCMC) sampling augmented with reversible jumps between models with different numbers of parameters, we characterise the posterior distributions of the model parameters and the number of components per burst. We relate these model parameters to physical quantities in the system, and show for the first time that the variability within a burst does not conform to predictions from ideas of self-organised criticality. We also examine how well the properties of the spikes fit the predictions of simplified cascade models for the different trigger mechanisms.

Since launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has detected many hundreds of bursts from magnetar sources. While the vast majority of these bursts have been attributed to several known magnetars, there is also a small sample of magnetar-like bursts of unknown origin. Here we present the Fermi/GBM magnetar catalog, giving the results of the temporal and spectral analyses of 440 magnetar bursts with high temporal and spectral resolution. This catalog covers the first five years of GBM magnetar observations, from July 2008 to June 2013. We provide durations, spectral parameters for various models, fluences and peak fluxes for all the bursts, as well as a detailed temporal analysis for SGR J1550-5418 bursts. Finally, we suggest that some of the bursts of unknown origin are associated with the newly discovered magnetar 3XMM J185246.6+0033.7.

Quasi-periodic oscillations (QPOs) observed in the giant flares of magnetars are of particular interest due to their potential to open up a window into the neutron star interior via neutron star asteroseismology. However, only three giant flares have been observed. We therefore make use of the much larger data set of shorter, less energetic recurrent bursts. Here, we report on a search for QPOs in a large data set of bursts from the two most burst-active magnetars, SGR 1806-20 and SGR 1900+14, observed with the Rossi X-ray Timing Explorer (RXTE). We find a single detection in an averaged periodogram comprising 30 bursts from SGR 1806-20, with a frequency of 57 Hz and a width of 5 Hz, remarkably similar to a giant flare QPO observed from SGR 1900+14. This QPO fits naturally within the framework of global magneto-elastic torsional oscillations employed to explain the giant flare QPOs. Additionally, we uncover a limit on the applicability of Fourier analysis for light curves with low background count rates and strong variability on short timescales. In this regime, standard Fourier methodology and more sophisticated Fourier analyses fail in equal parts by yielding an unacceptably large number of false positive detections. This problem is not straightforward to solve in the Fourier domain. Instead, we show how simulations of light curves can offer a viable solution for QPO searches in these light curves.

Quasi-periodic oscillations (QPOs) detected in the 2004 giant flare from SGR 1806-20 are often interpreted as global magneto-elastic oscillations of the neutron star. There is, however, a large discrepancy between theoretical models, which predict that the highest frequency oscillations should die out rapidly, and the observations, which suggested that the highest-frequency signals persisted for ~100s in X-ray data from two different spacecraft. This discrepancy is particularly important for the high-frequency QPO at ~625 Hz. However, previous analyses did not systematically test whether the signal could also be there in much shorter data segments, more consistent with the theoretical predictions. Here, we test for the presence of the high-frequency QPO at 625 Hz in data from both the Rossi X-ray Timing Explorer (RXTE) and the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) systematically both in individual rotational cycles of the neutron star, as well as averaged over multiple successive rotational cycles at the same phase. We find that the QPO in the RXTE data is consistent with being only present in a single cycle, for a short duration of ~0.5s, whereas the RHESSI data are as consistent with a short-lived signal that appears and disappears as with a long-lived QPO. Taken together, this data provides evidence for strong magnetic interaction between the crust and the core.

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. The scarcity of giant flares makes a search for QPOs in the shorter, far more numerous bursts from Soft Gamma Repeaters (SGRs) desirable. In Huppenkothen et al (2013), we developed a Bayesian method for searching for QPOs in short magnetar bursts, taking into account the effects of the complicated burst structure, and have shown its feasibility on a small sample of bursts. Here, we apply the same method to a much larger sample from a burst storm of 286 bursts from SGR J1550-5418. We report a candidate signal at 260 Hz in a search of the individual bursts, which is fairly broad. We also find two QPOs at 93 Hz and one at 127 Hz, when averaging periodograms from a number of bursts in individual triggers, at frequencies close to QPOs previously observed in magnetar giant flares. Finally, for the first time, we explore the overall burst variability in the sample, and report a weak anti-correlation between the power-law index of the broadband model characterising aperiodic burst variability, and the burst duration: shorter bursts have steeper power law indices than longer bursts. This indicates that longer bursts vary over a broader range of time scales, and are not simply longer versions of the short bursts.

We report on time-resolved spectroscopy of the 63 brightest bursts of SGR J1550-5418, detected with Fermi/Gamma-ray Burst Monitor during its 2008-2009 intense bursting episode. We performed spectral analysis down to 4 ms time-scales, to characterize the spectral evolution of the bursts. Using a Comptonized model, we find that the peak energy, E_peak, anti-correlates with flux, while the low-energy photon index remains constant at -0.8 up to a flux limit F~10^-5 erg s-1 cm-2. Above this flux value the E_peak-flux correlation changes sign, and the index positively correlates with flux reaching 1 at the highest fluxes. Using a two black-body model, we find that the areas and fluxes of the two emitting regions correlate positively. Further, we study here for the first time, the evolution of the temperatures and areas as a function of flux. We find that the area-kT relation follows lines of constant luminosity at the lowest fluxes, R^2 \propto kT^-4, with a break at higher fluxes ($F>10^-5.5 erg s-1 cm-2). The area of the high-kT component increases with flux while its temperature decreases, which we interpret as due to an adiabatic cooling process. The area of the low-kT component, on the other hand, appears to saturate at the highest fluxes, towards R_max~30 km. Assuming that crust quakes are responsible for SGR bursts and considering R_max as the maximum radius of the emitting photon-pair plasma fireball, we relate this saturation radius to a minimum excitation radius of the magnetosphere, and put a lower limit on the internal magnetic field of SGR J1550-5418, B_int>~4.5x10^15 G.

In order to discern the physical nature of many gamma-ray sources in the sky, we must look not only in spectral and spatial dimensions, but also understand their temporal variability. However, timing analysis of sources with a highly transient nature, such as magnetar bursts, is difficult: standard Fourier techniques developed for long-term variability generally observed, for example, from AGN often do not apply. Here, we present newly developed timing methods applicable to transient events of all kinds, and show their successful application to magnetar bursts observed with Fermi/GBM. Magnetars are a prime subject for timing studies, thanks to the detection of quasi-periodicities in magnetar Giant Flares and their potential to help shed light on the structure of neutron stars. Using state-of-the art statistical techniques, we search for quasi-periodicities (QPOs) in a sample of bursts from Soft Gamma Repeater SGR J0501+4516 observed with Fermi/GBM and provide upper limits for potential QPO detections. Additionally, for the first time, we characterise the broadband variability behaviour of magnetar bursts and highlight how this new information could provide us with another way to probe these mysterious objects.

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. However, with only three giant flares ever recorded, and only two with data of sufficient quality to search for QPOs, such analysis is seriously data limited. We set out a procedure for doing QPO searches in the far more numerous, short, less energetic magnetar bursts. The short, transient nature of these bursts requires the implementation of sophisticated statistical techniques to make reliable inferences. Using Bayesian statistics, we model the periodogram as a combination of red noise at low frequencies and white noise at high frequencies, which we show is a conservative approach to the problem. We use empirical models to make inferences about the potential signature of periodic and quasi-periodic oscillations at these frequencies. We compare our method with previously used techniques and find that although it is on the whole more conservative, it is also more reliable in ruling out false positives. We illustrate our Bayesian method by applying it to a sample of 27 bursts from the magnetar SGR J0501+4516 observed by the Fermi Gamma-ray Burst Monitor, and we find no evidence for the presence of QPOs in any of the bursts in the unbinned spectra, but do find a candidate detection in the binned spectra of one burst. However, whether this signal is due to a genuine quasi-periodic process, or can be attributed to unmodeled effects in the noise is at this point a matter of interpretation.

Jets around low- and intermediate-mass young stellar objects (YSOs) contain a fossil record of the recent accretion and outflow activity of their parent star-forming systems. We aim to understand whether the accretion/ejection process is similar across the entire stellar mass range of the parent YSOs. To this end we have obtained VLT/X-shooter spectra of HH 1042 and HH 1043, two newly discovered jets in the massive star-forming region RCW 36. HH 1042 is associated with the intermediate-mass YSO 08576nr292. Over 90 emission lines are detected in the spectra. High-velocity (up to 220 km/s) blue- and redshifted emission from a bipolar flow is observed in typical shock tracers. Low-velocity emission from the background cloud is detected in nebular tracers, including lines from high ionization species. We applied combined optical and infrared spectral diagnostic tools in order to derive the physical conditions (density, temperature, and ionization) in the jets. The measured mass outflow rates are Mjet ~ 10^-7 Msun/yr. We measure a high accretion rate for HH 1042 (Macc ~ 10^-6 Msun/yr) and Mjet/Macc ~ 0.1, comparable to low-mass sources and consistent with models for magneto-centrifugal jet launching. The knotted structure and velocity spread in both jets are interpreted as fossil signatures of a variable outflow rate. The mean velocities in both lobes of the jets are comparable, but the variations in Mjet and velocity in the two lobes are not symmetric, suggesting that the launching mechanism on either side of the accretion disk is not synchronized. For HH 1042, we have constructed an interpretative physical model with a stochastic or periodic outflow rate and a description of a ballistic flow as its constituents. The knotted structure and velocity spread can be reproduced qualitatively with the model, indicating that the outflow velocity varies on timescales on the order of 100 yr.

In early October 2008, the Soft Gamma Repeater SGRJ1550 - 5418 (1E 1547.0 - 5408, AXJ155052 - 5418, PSR J1550 - 5418) became active, emitting a series of bursts which triggered the Fermi Gamma-ray Burst Monitor (GBM) after which a second especially intense activity period commenced in 2009 January and a third, less active period was detected in 2009 March-April. Here we analyze the GBM data all the bursts from the first and last active episodes. We performed temporal and spectral analysis for all events and found that their temporal characteristics are very similar to the ones of other SGR bursts, as well the ones reported for the bursts of the main episode (average burst durations \sim 170 ms). In addition, we used our sample of bursts to quantify the systematic uncertainties of the GBM location algorithm for soft gamma-ray transients to < 8 deg. Our spectral analysis indicates significant spectral evolution between the first and last set of events. Although the 2008 October events are best fit with a single blackbody function, for the 2009 bursts an Optically Thin Thermal Bremsstrahlung (OTTB) is clearly preferred. We attribute this evolution to changes in the magnetic field topology of the source, possibly due to effects following the very energetic main bursting episode.

Supernova remnants (SNRs) are a complex class of sources, and their heterogeneous nature has hindered the characterization of their general observational properties. To overcome this challenge, we use statistical tools to analyze the Chandra X-ray images of Galactic and Large Magellanic Cloud SNRs. We apply two techniques, a power-ratio method (a multipole expansion) and wavelet-transform analysis, to measure the global and local morphological properties of the X-ray line and thermal emission in twenty-four SNRs. We find that Type Ia SNRs have statistically more spherical and mirror symmetric thermal X-ray emission than core-collapse (CC) SNRs. The ability to type SNRs based on thermal emission morphology alone enables, for the first time, the typing of SNRs with weak X-ray lines or with low-resolution spectra. We identify one source, SNR G344.7-0.1, as originating from a CC explosion that was previously unknown, and we confirm the tentative Type Ia classifications of G337.2-0.7 and G272.2-3.2. Although the global morphology is indicative of the explosion type, the relative morphology of the X-ray line emission within SNRs is not: all sources in our sample have well-mixed ejecta, irrespective of stellar origin. In particular, we find that 90% of the bright metal-line emitting substructures are spatially coincident and have similar scales, even if the metals arise from different burning processes. Moreover, the overall X-ray line morphologies within each SNR are the same, with <6% differences. These findings reinforce that hydrodynamical instabilities can efficiently mix ejecta in Type Ia and CC SNRs. The only exception is W49B, which is likely from its jet-driven/bipolar SN explosion. Finally, we describe observational constraints that can be used to test hydrodynamical models of SNR evolution; notably, the filling factor of X-ray emission decreases with SNR age.

We present a new observational method to type the explosions of young supernova remnants (SNRs). By measuring the morphology of the Chandra X-ray line emission in seventeen Galactic and Large Magellanic Cloud SNRs with a multipole expansion analysis (using power ratios), we find that the core-collapse SNRs are statistically more asymmetric than the Type Ia SNRs. We show that the two classes of supernovae can be separated naturally using this technique because X-ray line morphologies reflect the distinct explosion mechanisms and structure of the circumstellar material. These findings are consistent with recent spectropolarimetry results showing that core-collapse SNe are intrinsically more asymmetric.