Search

LaCrGe$_3$ has attracted attention as a paradigm example of the avoidance of ferromagnetic (FM) quantum criticality in an itinerant magnet. By combining thermodynamic, transport, x-ray and neutron scattering as well as $\mu$SR measurements, we refined the temperature-pressure phase diagram of LaCrGe$_3$. We provide thermodynamic evidence (i) for the first-order character of the FM transition when it is suppressed to low temperatures and (ii) for the formation of new phases at high pressures. From our microscopic data, we infer that short-range FM ordered clusters exist in these high-pressure phases. These results suggest that LaCrGe$_3$ is a rare example, which fills the gap between the two extreme limits of avoided FM quantum criticality in clean and strongly disordered metals., Comment: 7 pages, 4 figures plus Supplemental Information

We present a new optical sample of three Supernova Remnants and 16 Supernova Remnant (SNR) candidates in the Large Magellanic Cloud(LMC). These objects were originally selected using deep H$\alpha$, [SII] and [OIII] narrow-band imaging. Most of the newly found objects are located in less dense regions, near or around the edges of the LMC's main body. Together with previously suggested MCSNR J0541-6659, we confirm the SNR nature for two additional new objects: MCSNR J0522-6740 and MCSNRJ0542-7104. Spectroscopic follow-up observations for 12 of the LMC objects confirm high [SII]/H$\alpha$ a emission-line ratios ranging from 0.5 to 1.1. We consider the candidate J0509-6402 to be a special example of the remnant of a possible Type Ia Supernova which is situated some 2$^\circ$ ($\sim 1.75$kpc) north from the main body of the LMC. We also find that the SNR candidates in our sample are significantly larger in size than the currently known LMC SNRs by a factor of $\sim 2$. This could potentially imply that we are discovering a previously unknown but predicted, older class of large LMC SNRs that are only visible optically. Finally, we suggest that most of these LMC SNRs are residing in a very rarefied environment towards the end of their evolutionary span where they become less visible to radio and X-ray telescopes., Comment: 24 pages, 25 figures, to appear in MNRAS

eROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the primary instrument on the Spectrum-Roentgen-Gamma (SRG) mission, which was successfully launched on July 13, 2019, from the Baikonour cosmodrome. After the commissioning of the instrument and a subsequent calibration and performance verification phase, eROSITA started a survey of the entire sky on December 13, 2019. By the end of 2023, eight complete scans of the celestial sphere will have been performed, each lasting six months. At the end of this program, the eROSITA all-sky survey in the soft X-ray band (0.2--2.3\,keV) will be about 25 times more sensitive than the ROSAT All-Sky Survey, while in the hard band (2.3--8\,keV) it will provide the first ever true imaging survey of the sky. The eROSITA design driving science is the detection of large samples of galaxy clusters up to redshifts $z>1$ in order to study the large-scale structure of the universe and test cosmological models including Dark Energy. In addition, eROSITA is expected to yield a sample of a few million AGNs, including obscured objects, revolutionizing our view of the evolution of supermassive black holes. The survey will also provide new insights into a wide range of astrophysical phenomena, including X-ray binaries, active stars, and diffuse emission within the Galaxy. Results from early observations, some of which are presented here, confirm that the performance of the instrument is able to fulfil its scientific promise. With this paper, we aim to give a concise description of the instrument, its performance as measured on ground, its operation in space, and also the first results from in-orbit measurements., Comment: 16 pages, 19 figures, accepted by Astronomy & Astrophysics

Supersoft X-ray sources (SSS) have been identified as white dwarfs accreting from binary companions and undergoing nuclear-burning of the accreted material on their surface. Although expected to be a relatively numerous population from both binary evolution models and their identification as Type Ia supernova progenitor candidates, given the very soft spectrum of SSSs relatively few are known. Here we report on the X-ray and optical properties of 1RXS J050526.3-684628, a previously unidentified accreting nuclear-burning white dwarf located in the Large Magellanic Cloud (LMC). XMM-Newton observations enabled us to study its X-ray spectrum and measure for the first time short period oscillations of ~170 s. By analysing newly obtained X-ray data by eROSITA, together with Swift observations and archival ROSAT data, we have followed its long-term evolution over the last 3 decades. We identify 1RXS J050526.3-684628 as a slowly-evolving post-nova SSS undergoing residual surface nuclear-burning, which finally reached its peak in 2013 and is now declining. Though long expected on theoretical grounds, such long-lived residual-burning objects had not yet been found. By comparison with existing models, we find that the effective temperature and luminosity evolution are consistent with a ~0.7 $M_{\odot}$ carbon-oxygen white dwarf accreting ~10$^{-9}~\rm{M}_{\odot}$/yr. Our results suggest there may be many more undiscovered SSSs and "missed" novae awaiting dedicated deep X-ray searches in the LMC and elsewhere., Comment: 9 pages, 4 figures, accepted for publication in MNRAS

eROSITA, launched on 13 July 2019, will be completing the first all-sky survey in the soft and medium X-ray band in nearly three decades. This 4-year survey, finishing in late 2023, will present a rich legacy for the entire astrophysics community and complement upcoming multi-wavelength surveys (with, e.g. the Large Synoptic Survey Telescope and the Dark Energy Survey). Besides the major scientific aim to study active galactic nuclei (AGN) and galaxy clusters, eROSITA will contribute significantly to X-ray studies of normal (i.e., not AGN) galaxies. Starting from multi-wavelength catalogues, we measure star formation rates and stellar masses for 60 212 galaxies constrained to distances of 50-200 Mpc. We chose this distance range to focus on the relatively unexplored volume outside the local Universe, where galaxies will be largely spatially unresolved and probe a range of X-ray luminosities that overlap with the low luminosity and/or highly obscured AGN population. We use the most recent X-ray scaling relations as well as the on-orbit eROSITA instrument performance to predict the X-ray emission from XRBs and diffuse hot gas and to perform both an analytic prediction and an end-to-end simulation using the mission simulation software, SIXTE. We consider potential contributions from hidden AGN and comment on the impact of normal galaxies on the measurement of the faint end of the AGN luminosity function. We predict that the eROSITA 4-year survey, will detect $\gtrsim$ 15 000 galaxies (3 $\sigma$ significance) at 50-200 Mpc, which is ~100X more normal galaxies than detected in any X-ray survey to date., Comment: 18 pages, 15 figures. Accepted for publication in MNRAS

In late 2019, ASHRAE hosted the Great Energy Predictor III (GEPIII) machine learning competition on the Kaggle platform. This launch marked the third energy prediction competition from ASHRAE and the first since the mid-1990s. In this updated version, the competitors were provided with over 20 million points of training data from 2,380 energy meters collected for 1,448 buildings from 16 sources. This competition's overall objective was to find the most accurate modeling solutions for the prediction of over 41 million private and public test data points. The competition had 4,370 participants, split across 3,614 teams from 94 countries who submitted 39,403 predictions. In addition to the top five winning workflows, the competitors publicly shared 415 reproducible online machine learning workflow examples (notebooks), including over 40 additional, full solutions. This paper gives a high-level overview of the competition preparation and dataset, competitors and their discussions, machine learning workflows and models generated, winners and their submissions, discussion of lessons learned, and competition outputs and next steps. The most popular and accurate machine learning workflows used large ensembles of mostly gradient boosting tree models, such as LightGBM. Similar to the first predictor competition, preprocessing of the data sets emerged as a key differentiator.

We consider a second-order elliptic boundary value problem with strongly monotone and Lipschitz-continuous nonlinearity. We design and study its adaptive numerical approximation interconnecting a finite element discretization, the Banach-Picard linearization, and a contractive linear algebraic solver. We in particular identify stopping criteria for the algebraic solver that on the one hand do not request an overly tight tolerance but on the other hand are sufficient for the inexact (perturbed) Banach-Picard linearization to remain contractive. Similarly, we identify suitable stopping criteria for the Banach-Picard iteration that leave an amount of linearization error that is not harmful for the residual a-posteriori error estimate to steer reliably the adaptive mesh-refinement. For the resulting algorithm, we prove a contraction of the (doubly) inexact iterates after some amount of steps of mesh-refinement / linerization / algebraic solver, leading to its linear convergence. Moreover, for usual mesh-refinement rules, we also prove that the overall error decays at the optimal rate with respect to the number of elements (degrees of freedom) added with respect to the initial mesh. Finally, we prove that our fully adaptive algorithm drives the overall error down with the same optimal rate also with respect to the overall algorithmic cost expressed as the cumulated sum of the number of mesh elements over all mesh-refinement, linearization, and algebraic solver steps. Numerical experiments support these theoretical findings and illustrate the optimal overall algorithmic cost of the fully adaptive algorithm on several test cases.

Confined microexplosions induced in silicon by powerful ultrashort laser pulses can lead to new Si phases. Some of these have not previously been observed via near-equilibrium compression of silicon. In this study, confocal Raman micro-spectroscopy and Raman imaging of arrays of microexplosions have been conducted to search for Raman signatures of these novel allotropes of silicon. A microexplosion is generated at the interface between a thick silicon dioxide confinement layer and underlying silicon. It is characterised by a void at the interface above a region of compressed silicon. Raman data show a rich assembly of silicon phases within the modified silicon. Residual stresses up to 4.5 GPa in the modifications have been determined from the shift in the main diamond-cubic Si Raman peak. The computed Raman spectra for a number of Si allotropes show reasonable agreement with the experimental spectra. Two structurally similar tetragonal phases of silicon (the rhombohedral r8 and the body-centred bc8) phases as well as recently identified bt8-Si are all highly likely to be contained in Raman spectra from many laser-modified sites. Although the st12-Si phase, previously observed in our electron diffraction studies of the highly compressively stressed laser-modified regions, was not reliably identified from Raman data, we suggest this could be due to the possible difference in residual stress level in the sites analysed by the electron diffraction and Raman spectra. Several other unidentified Raman peaks were observed, suggesting the presence of other unknown silicon phases. All of these silicon phases are expected to have attractive semiconducting properties including narrow band gap that open up novel applications., Comment: 18 pages, 6 figures

We present results of the first dedicated observation of the young X-ray pulsar SXP 1062 in the broad X-ray energy band obtained during its 2019 outburst with the NuSTAR and XMM-Newton observatories. The analysis of the pulse-phase averaged and phase-resolved spectra in the energy band from 0.5 to 70 keV did not reveal any evidence for the presence of a cyclotron line. The spin period of the pulsar was found to have decreased to 979.48+/-0.06 s implying a ~10% reduction compared to the last measured period during the monitoring campaign conducted about five years ago, and is puzzling considering that the system apparently did not show major outbursts ever since. The switch of the pulsar to the spin-up regime supports the common assumption that torques acting on the accreting neutron star are nearly balanced and thus SXP 1062 likely also spins with a period close to the equilibrium value for this system. The current monitoring of the source revealed also a sharp drop of its soft X-ray flux right after the outburst, which is in drastic contrast to the behaviour during the previous outburst when the pulsar remained observable for years with only a minor flux decrease after the end of the outburst. This unexpected off state of the source lasted for at most 20 days after which SXP 1062 returned to the level observed during previous campaigns. We discuss this and other findings in context of the modern models of accretion onto strongly magnetized neutron stars., Comment: 7 pages, 6 figures; accepted for publication in A&A

We consider adaptive finite element methods for second-order elliptic PDEs, where the arising discrete systems are not solved exactly. For contractive iterative solvers, we formulate an adaptive algorithm which monitors and steers the adaptive mesh-refinement as well as the inexact solution of the arising discrete systems. We prove that the proposed strategy leads to linear convergence with optimal algebraic rates. Unlike prior works, however, we focus on convergence rates with respect to the overall computational costs. In explicit terms, the proposed adaptive strategy thus guarantees quasi-optimal computational time. In particular, our analysis covers linear problems, where the linear systems are solved by an optimally preconditioned CG method as well as nonlinear problems with strongly monotone nonlinearity which are linearized by the so-called Zarantonello iteration.

The application of standard accretion theory to observations of X-ray binaries provides valuable insights into neutron star properties, such as their spin period and magnetic field. However, most studies concentrate on relatively old systems, where the neutron star is in its late propeller, accretor, or nearly spin equilibrium phase. Here we use an analytic model from standard accretion theory to illustrate the evolution of high-mass X-ray binaries early in their life. We show that a young neutron star is unlikely to be an accretor because of the long duration of ejector and propeller phases. We apply the model to the recently discovered ~4000 yr old high-mass X-ray binary XMMU J051342.6-672412 and find that the system's neutron star, with a tentative spin period of 4.4 s, cannot be in the accretor phase and has a magnetic field B > (a few)x10^13 G, which is comparable to the magnetic field of many older high-mass X-ray binaries and is much higher than the spin equilibrium inferred value of (a few)x10^11 G. The observed X-ray luminosity could be the result of thermal emission from a young cooling magnetic neutron star or a small amount of accretion that can occur in the propeller phase., Comment: 7 pages, 4 figures; accepted for publication in MNRAS

Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and strongly affected by high background. In 2018, we obtained a new, deeper (~90 ks) XMM-Newton observation that allowed us to perform a detailed characterization of the ULXs hosted in the galaxy. ULX-1 and ULX-2, the two brightest objects (Lx~(6-10)e39 erg/s), have spectra that can be described by a two-thermal component model as often found in ULXs. ULX-1 shows also one order of magnitude in flux variability on short-term timescales (hundreds to thousand ks). The other sources (ULX-3 and ULX-4) show flux changes of at least an order of magnitude, and these objects may be candidate transient ULXs although longer X-ray monitoring or further studies are required to ascribe them to the ULX population. In addition, we found a previously undetected source that might be a new candidate ULX (labelled as ULX-5) with a luminosity of ~1e39 erg/s and hard power-law spectral shape, whose nature is still unclear and for which a background Active Galactic Nucleus cannot be excluded. We discuss the properties of all the ULXs in NGC 7456 within the framework of super-Eddington accretion onto stellar mass compact objects. Although no pulsations were detected, we cannot exclude that the sources host neutron stars., Comment: Accepted on ApJ; 10 pages, 5 figures, 3 tables

Here we present a radio continuum study based on new and archival data from the Australia Telescope Compact Array towards N 103B, a young (<=1000 yrs) spectroscopically confirmed type Ia SNR in the Large Magellanic Cloud (LMC) and proposed to have originated from a single degenerate progenitor. The radio morphology of this SNR is asymmetrical with two bright regions towards the north-west and south-west of the central location as defined by radio emission. N 103B identified features include: a radio spectral index of -0.75+-0.01 (consistent with other young type Ia SNRs in the Galaxy); a bulk SNR expansion rate as in X-rays; morphology and polarised electrical field vector measurements where we note radial polarisation peak towards the north-west of the remnant at both 5500 and 9000 MHz. The spectrum is concave-up and the most likely reason is the non-linear diffusive shock acceleration effects or presence of two different populations of ultra-relativistic electrons. We also note unpolarized clumps near the south-west region which is in agreement with this above scenario. We derive a typical magnetic field strength for N 103B, of 16.4 microG for an average rotation measurement of 200 rad m^-2. However, we estimate the equipartition field to be of the order of ~235 microG with an estimated minimum energy of Emin=6.3*10^48 erg. The close (~0.5 degree) proximity of N 103B to the LMC mid-plane indicates that an early encounter with dense interstellar medium may have set an important constrain on SNR evolution. Finally, we compare features of N 103B, to six other young type Ia SNRs in the LMC and Galaxy, with a range of proposed degeneracy scenarios to highlight potential differences due to a different models. We suggest that the single degenerate scenario might point to morphologically asymmetric type Ia supernova explosions.

Aims: We use new deep XMM-Newton observations of the northern disk of M 31 to trace the hot interstellar medium (ISM) in unprecedented detail and to characterise the physical properties of the X-ray emitting plasmas. Methods: We used all XMM-Newton data up to and including our new observations to produce the most detailed image yet of the hot ISM plasma in a grand design spiral galaxy such as our own. We compared the X-ray morphology to multi-wavelength studies in the literature to set it in the context of the multi-phase ISM. We performed spectral analyses on the extended emission using our new observations as they offer sufficient depth and count statistics to constrain the plasma properties. Data from the Panchromatic Hubble Andromeda Treasury were used to estimate the energy injected by massive stars and their supernovae. We compared these results to the hot gas properties. Results: The brightest emission regions were found to be correlated with populations of massive stars, notably in the 10 kpc star-forming ring. The plasma temperatures in the ring regions are ~0.2 keV up to ~0.6 keV. We suggest this emission is hot ISM heated in massive stellar clusters and superbubbles. We derived X-ray luminosities, densities, and pressures for the gas in each region. We also found large extended emission filling low density gaps in the dust morphology of the northern disk, notably between the 5 kpc and 10 kpc star-forming rings. We propose that the hot gas was heated and expelled into the gaps by the populations of massive stars in the rings. Conclusions. It is clear that the massive stellar populations are responsible for heating the ISM to X-ray emitting temperatures, filling their surroundings, and possibly driving the hot gas into the low density regions. Overall, the morphology and spectra of the hot gas in the northern disk of M 31 is similar to other galaxy disks., Comment: 23 pages, 14 figures, accepted for publication in Astronomy & Astrophysics

The nearby galaxy NGC 300 is hosting two luminous transient supersoft X-ray sources with bolometric luminosities above 3x10^38 erg/s, assuming simple black-body spectra with temperatures around 60-70 eV. For one of these, SSS1, a periodic modulation of 5.4h was observed in an XMM-Newton observation from 1st of January 2001 lasting 47 ks, but not visible 6 days earlier when the luminosity was higher. We report here the detection of a new outburst from this source, which occurred during two more recent XMM-Newton observations performed on 17 to 20 December 2016 lasting for 310 ks. The luminosity was similar as in December 2000, and the 0.2-2.0 keV light curve revealed again a periodic modulation, with a period of 4.68+-0.26h, significant only in the first of the two observations. Taking into account the large uncertainties (the 2001 period was re-estimated at 5.7+-1.1h), the two values could be marginally compatible, and maybe associated with an orbital period, although the signal strength is highly variable. Thanks to the new long exposures, an additional absorption feature is now visible in the spectra, that we modelled with an absorption edge. This component decreases the bolometric luminosity below 3x10^38 erg/s and would therefore allow the presence of a white dwarf with a mass close to the Chandrasekhar limit. The system was found in outburst in 1992, 2000, 2008, and 2016 suggesting a possible recurrence period of about 8 years. We discuss viable models involving white dwarfs, neutron stars or black holes., Comment: 7 pages, 7 figures

We report the discovery of a very young high-mass X-ray binary (HMXB) system associated with the supernova remnant (SNR) MCSNRJ0513-6724 in the Large Magellanic Cloud (LMC), using XMM-Newton X-ray observations. The HMXB is located at the geometrical centre of extended soft X-ray emission, which we confirm as an SNR. The HMXB spectrum is consistent with an absorbed power law with spectral index ~1.6 and a luminosity of 7x10^{33} ergs/s (0.2--12 keV). Tentative X-ray pulsations are observed with a periodicity of 4.4 s and the OGLE I-band light curve of the optical counterpart from more than 17.5 years reveals a period of 2.2324\pm0.0003 d, which we interpret as the orbital period of the binary system. The X-ray spectrum of the SNR is consistent with non-equilibrium shock models as expected for young/less evolved SNRs. From the derived ionisation time scale we estimate the age of the SNR to be <6 kyr. The association of the HMXB with the SNR makes it the youngest HMXB, in the earliest evolutionary stage known to date. A HMXB as young as this can switch on as an accreting pulsar only when the spin period has reached a critical value. Under this assumption, we obtain an upper limit to the magnetic field of < 5x10^{11} G. This implies several interesting possibilities including magnetic field burial, possibly by an episode of post-supernova hyper-critical accretion. Since these fields are expected to diffuse out on a timescale of 10^{3}-10^{4} years, the discovery of a very young HMXB can provide us the unique opportunity to observe the evolution of the observable magnetic field for the first time in X-ray binaries., Comment: 9 pages, 8 figures, accepted for publication in MNRAS

We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$\sigma$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$\sigma$ limiting luminosity in the 4$-$25 keV band of $\sim$ $10^{35}$ erg s$^{-1}$, allowing us to probe fainter X-ray binary (XRB) populations than has been possible with other extragalactic NuSTAR surveys. We used hard X-ray colors and luminosities to constrain the compact-object type, exploiting the spectral differences between accreting black holes and neutron stars at $E>10$ keV. Several of our sources demonstrate variability consistent with previously observed behavior. We confirmed pulsations for seven pulsars in our 3$\sigma$ sample. We present the first detection of pulsations from a Be-XRB, SXP305 (CXO J005215.4$-$73191), with an X-ray pulse period of $305.69\pm0.16$ seconds and a likely orbital period of $\sim$1160-1180 days. Bright sources ($\gtrsim 5\times 10^{36}$ erg s$^{-1}$) in our sample have compact-object classifications consistent with their previously reported types in the literature. Lower luminosity sources ($\lesssim 5\times 10^{36}$ erg s$^{-1}$) have X-ray colors and luminosities consistent with multiple classifications. We raise questions about possible spectral differences at low luminosity between SMC pulsars and the Galactic pulsars used to create the diagnostic diagrams., Comment: 8 tables, 18 figures

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies., Comment: 18 pages, 23 figues, 9 tables, Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Aims: We present a comprehensive study of the supernova remnant (SNR) population of the Small Magellanic Cloud (SMC). We measure multiwavelength properties of the SMC SNRs and compare them to those of the Large Magellanic Cloud (LMC) population. Methods: This study combines the large dataset of XMM-Newton observations of the SMC, archival and recent radio continuum observations, an optical line emission survey, and new optical spectroscopic observations. We can thus build a complete and clean sample of 19 confirmed and 4 candidate SNRs. The homogeneous X-ray spectral analysis allows to search for SN ejecta and Fe K line emission, and to measure interstellar medium (ISM) abundances. We estimate the ratio of core-collapse to type Ia supernova rates of the SMC based on the X-ray properties and the local stellar environment of each SNR. Results : After the removal of unconfirmed or misclassified objects, and the addition of two newly confirmed SNRs based on multi-wavelength features, we present a final list of 21 confirmed SNRs and 2 candidates. While no Fe K line is detected even for the brightest and youngest SNR, we find X-ray evidence of SN ejecta in 11 SNRs. We estimate a ratio of 4.7$_{-1.9} ^{+0.6}$ core-collapse supernova to every type Ia SN, three times higher than in the LMC. The difference can be ascribed to the absence of the enhanced star formation episode in the SMC, which occurred in the LMC 0.5-1.5 Gyr ago. The hot-gas abundances of O, Ne, Mg, and Fe are 0.1-0.2 times solar. Their ratios with respect to SMC stellar abundances reflect the effects of dust depletion and partial dust destruction in SNR shocks. We find evidence that the ambient medium probed by SMC SNRs is less disturbed and less dense on average than in the LMC, consistent with the different morphologies of the two galaxies., Comment: Submitted to Astronomy & Astrophysics. 29 pages, 8 figures, 5 tables; 3 additional tables and 9 figures in Appendices

We examine the new Galactic supernova remnant (SNR) candidate, G23.11+0.18, as seen by the Murchison Widefield Array (MWA) radio telescope. We describe the morphology of the candidate and find a spectral index of -0.63+/-0.05 in the 70-170MHz domain. A coincident TeV gamma-ray detection in High-Energy Stereoscopic System (HESS) data supports the SNR nature of G23.11+0.18 and suggests that G23.11+0.18 is accelerating particles beyond TeV energies, thus making this object a promising new cosmic ray hadron source candidate. The remnant cannot be seen in current optical, infrared and X-ray data-sets. We do find, however, a dip in CO-traced molecular gas at a line-of-sight velocity of ~85 km/s, suggesting the existence of a G23.11+0.18 progenitor wind-blown bubble. Furthermore, the discovery of molecular gas clumps at a neighbouring velocity towards HESS J1832-085 adheres to the notion that a hadronic gamma-ray production mechanism is plausible towards the north of the remnant. Based on these morphological arguments, we propose an interstellar medium association for G23.11+0.18 at a kinematic distance of 4.6+/-0.8 kpc., Comment: 22 pages, 12 figures, Accepted for publication in The Astrophysical Journal

We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$)., Comment: Accepted for publication in The Astrophysical Journal

NGC300 ULX1 is an ultra-luminous X-ray pulsar, showing an unprecedented spin evolution, from about 126 s to less than 20 s in only 4 years, consistent with steady mass accretion rate. Following its discovery we have been monitoring the system with Swift/XRT and NICER to further study its properties. We found that even though the observed flux of the system dropped by a factor of $\gtrsim$20, the spin-up rate remained almost constant. A possible explanation is that the decrease in the observed flux is a result of increased absorption of obscuring material due to outflows or a precessing accretion disk., Comment: 7 pages, 3 figures, MNRAS in press

A clamp pressure cell for neutron scattering experiments at low temperatures and in external magnetic fields under pressure up to 2~GPa has been fabricated and tested. The cell provides optical access to the sample space that allows instantaneous pressure determination during sample loading, cooling and measuring using ruby and/or samarium doped strontium tetraborate fluorescence monitoring. A new calibration curve of the pressure-induced shift of the $^7D_0 - ^5F_0$ (0-0) line in the fluorescent spectrum of SrB$_4$O$_7$:Sm$^{2+}$ for moderate pressures, $P\leqslant{2}$ GPa, is given., Comment: 11 pages, 6 figures

We study several trace operators and spaces that are related to the bi-Laplacian. They are motivated by the development of ultraweak formulations for the bi-Laplace equation with homogeneous Dirichlet condition, but are also relevant to describe conformity of mixed approximations. Our aim is to have well-posed (ultraweak) formulations that assume low regularity, under the condition of an $L_2$ right-hand side function. We pursue two ways of defining traces and corresponding integration-by-parts formulas. In one case one obtains a non-closed space. This can be fixed by switching to the Kirchhoff-Love traces from [F\"uhrer, Heuer, Niemi, An ultraweak formulation of the Kirchhoff-Love plate bending model and DPG approximation, Math. Comp., 88 (2019)]. Using different combinations of trace operators we obtain two well-posed formulations. For both of them we report on numerical experiments with the DPG method and optimal test functions. In this paper we consider two and three space dimensions. However, with the exception of a given counterexample in an appendix (related to the non-closedness of a trace space), our analysis applies to any space dimension larger than or equal to two., Comment: 26 pages, 3 figures

In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$., Comment: Accepted for publication in A&A

In response to the comment posted by Nakashima et al. (arXiv:1903.1176), regarding prior claims for the features that we referred to as the Galactic Center Chimneys (2019, Nature, 567, 347), we point out the following: 1) The Nakashima et al. 2019 paper appeared in the arXiv on March 8th (1903.02571), after our paper was in the final stage of printing (accepted on January 30th). It is however interesting to see that the morphology of the brightest portions of the two results are in broad agreement (compare their Fig. 1 to our Extended Data Figs. 1 and 2). 2) Nakashima et al. 2013 ApJ 773, 20 claim the discovery of a blob of recombining plasma ~1deg south of Sgr A*, implying peculiar abundances. Again, their image (Fig. 1) agrees with the brightest portions of our images, although it does not show any direct connection between the plasma blob and the central parsec (e.g., such as the quasi-continuous chimney that we reported), nor evidence for an outflow from the center. We apologize for overlooking an appropriate citation to this contribution by Nakashima et al. 3) We fitted the XMM-Newton and Chandra data at the same position of the claimed recombining plasma and we did not find any clear-cut evidence for the presence of either an over-ionised plasma or peculiar abundances. Future X-ray calorimetric observations will presumably clarify this disagreement. 4) The continuity of the Chimney features, their quasi-symmetrical placement relative to Sgr A*, and their relatively sharp and well-defined edges are the essential features of our data that have led us to propose that the Chimneys are a unified columnar structure that represents a channel for the outflow of energy from the central region, possibly contributing to the stocking of the relativistic particle population manifested in the Fermi Bubbles.

Evidence has increasingly mounted in recent decades that outflows of matter and energy from the central parsecs of our Galaxy have shaped the observed structure of the Milky Way on a variety of larger scales. On scales of ~15 pc, the Galactic centre has bipolar lobes that can be seen in both X-rays and radio, indicating broadly collimated outflows from the centre, directed perpendicular to the Galactic plane. On far larger scales approaching the size of the Galaxy itself, gamma-ray observations have identified the so-called Fermi Bubble features, implying that our Galactic centre has, or has recently had, a period of active energy release leading to a production of relativistic particles that now populate huge cavities on both sides of the Galactic plane. The X-ray maps from the ROSAT all-sky survey show that the edges of these cavities close to the Galactic plane are bright in X-rays. At intermediate scales (~150 pc), radio astronomers have found the Galactic Centre Lobe, an apparent bubble of emission seen only at positive Galactic latitudes, but again indicative of energy injection from near the Galactic centre. Here we report the discovery of prominent X-ray structures on these intermediate (hundred-parsec) scales above and below the plane, which appear to connect the Galactic centre region to the Fermi bubbles. We propose that these newly-discovered structures, which we term the Galactic Centre Chimneys, constitute a channel through which energy and mass, injected by a quasi-continuous train of episodic events at the Galactic centre, are transported from the central parsecs to the base of the Fermi bubbles., Comment: Submitted version of paper published in Nature on March 20, 2019

What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates., Comment: 9 pages, 2 figures. Submitted to the Astro2020 Decadal Survey

We report the discovery of a new Small Magellanic Cloud Pulsar Wind Nebula (PWN) at the edge of the Supernova Remnant (SNR)-DEM S5. The pulsar powered object has a cometary morphology similar to the Galactic PWN analogs PSR B1951+32 and 'the mouse'. It is travelling supersonically through the interstellar medium. We estimate the Pulsar kick velocity to be in the range of 700-2000 km/s for an age between 28-10 kyr. The radio spectral index for this SNR PWN pulsar system is flat (-0.29 $\pm$ 0.01) consistent with other similar objects. We infer that the putative pulsar has a radio spectral index of -1.8, which is typical for Galactic pulsars. We searched for dispersion measures (DMs) up to 1000 cm/pc^3 but found no convincing candidates with a S/N greater than 8. We produce a polarisation map for this PWN at 5500 MHz and find a mean fractional polarisation of P $\sim 23$ percent. The X-ray power-law spectrum (Gamma $\sim 2$) is indicative of non-thermal synchrotron emission as is expected from PWN-pulsar system. Finally, we detect DEM S5 in Infrared (IR) bands. Our IR photometric measurements strongly indicate the presence of shocked gas which is expected for SNRs. However, it is unusual to detect such IR emission in a SNR with a supersonic bow-shock PWN. We also find a low-velocity HI cloud of $\sim 107$ km/s which is possibly interacting with DEM S5. SNR DEM S5 is the first confirmed detection of a pulsar-powered bow shock nebula found outside the Galaxy., Comment: 20 pages, 15 figures, accepted for publication in MNRAS

Previous XMM-Newton observations of the thermally emitting isolated neutron star RX J1605.3+3249 provided a candidate for a shallow periodic signal and evidence of a fast spin down, which suggested a high dipolar magnetic field and an evolution from a magnetar. We obtained a large programme with XMM-Newton to confirm its candidate timing solution, understand the energy-dependent amplitude of the modulation, and investigate the spectral features of the source. We performed extensive high-resolution and broadband periodicity searches in the new observations, using the combined photons of the three EPIC cameras and allowing for moderate changes of pulsed fraction and the optimal energy range for detection. A deep $4\sigma$ upper limit of $1.33(6)\%$ for modulations in the relevant frequency range conservatively rules out the candidate period previously reported. Blind searches revealed no other periodic signal above the $1.5\%$ level $(3\sigma$) in any of the four new observations. While theoretical models fall short at physically describing the complex energy distribution of the source, best-fit X-ray spectral parameters are obtained for a fully or partially ionized neutron star hydrogen atmosphere model with $B=10^{13}$ G, modified by a broad Gaussian absorption line at $385\pm10$ eV. The deep limits from the timing analysis disfavour equally well-fit double temperature blackbody models where both the star surface and small hotspots contribute to the X-ray flux of the source. We identified a low significance ($1\sigma$) temporal trend on the parameters of the source in the analysis of RGS data dating back to 2002, which may be explained by unaccounted calibration issues and spectral model uncertainties. The new dataset also shows no evidence of the previously reported narrow absorption feature at $\sim570$ eV, whose possible transient nature disfavours an atmospheric origin., Comment: 15 pages, 6 figures, accepted for publication in A&A

The accurate computation of the electrostatic capacity of three dimensional objects is a fascinating benchmark problem with a long and rich history. In particular, the capacity of the unit cube has widely been studied, and recent advances allow to compute its capacity to more than ten digits of accuracy. However, the accurate computation of the capacity for general three dimensional polyhedra is still an open problem. In this paper, we propose a new algorithm based on a combination of ZZ-type a posteriori error estimation and effective operator preconditioned boundary integral formulations to easily compute the capacity of complex three dimensional polyhedra to 5 digits and more. While this paper focuses on the capacity as a benchmark problem, it also discusses implementational issues of adaptive boundary element solvers, and we provide codes based on the boundary element package Bempp to make the underlying techniques accessible to a wide range of practical problems.

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively., Comment: 9 pages, 4 figures, submitted to ApJ, comments welcome

Nearby star-forming galaxies offer a unique environment to study the populations of young ($<$100 Myr) accreting binaries. These systems are tracers of past populations of massive stars that heavily affect their immediate environment and parent galaxies. Using a Chandra X-ray Visionary program, we investigate the young neutron-star binary population in the low metallicity of the Small Magellanic Cloud (SMC) by reaching quiescent X-ray luminosity levels ($\sim$few times $10^{32}$ erg/s). We present the first measurement of the formation efficiency of high-mass X-ray binaries (HMXBs) as a function of the age of their parent stellar populations by using 3 indicators: the number ratio of HMXBs to OB stars, to the SFR, and to the stellar mass produced during the specific star-formation burst they are associated with. In all cases, we find that the HMXB formation efficiency increases as a function of time up to $\sim$40-60 Myr, and then gradually decreases., Comment: 5 pages, 4 figures. To appear in Proc. IAUS 346: High-mass X-ray binaries: illuminating the passage from massive binaries to merging compact objects

Cyclotron lines, also called cyclotron resonant scattering features (CRSF) are spectral features, generally appearing in absorption, in the X-ray spectra of objects containing highly magnetized neutron stars, allowing the direct measurement of the magnetic field strength in these objects. Cyclotron features are thought to be due to resonant scattering of photons by electrons in the strong magnetic fields. The main content of this contribution focusses on electron cyclotron lines as found in accreting X-ray binary pulsars (XRBP) with magnetic fields on the order of several 1012 Gauss. Also, possible proton cyclotron lines from single neutron stars with even stronger magnetic fields are briefly discussed.With regard to electron cyclotron lines, we present an updated list of XRBPs that show evidence of such absorption lines. The first such line was discovered in a 1976 balloon observation of the accreting binary pulsar Hercules X-1, it is considered to be the first direct measurement of the magnetic field of a neutron star. As of today (mid 2018), we list 36 XRBPs showing evidence of one ore more electron cyclotron absorption line(s). A few have been measured only once and must be confirmed (several more objects are listed as candidates). In addition to the Tables of objects, we summarize the evidence of variability of the cyclotron line as a function of various parameters (especially pulse phase, luminosity and time), and add a discussion of the different observed phenomena and associated attempts of theoretical modeling. We also discuss our understanding of the underlying physics of accretion onto highly magnetized neutron stars. For proton cyclotron lines, we present tables with seven neutron stars and discuss their nature and the physics in these objects., Comment: 32 pages, 15 figures, 8 Tables, accepted by A&A 2018

NGC300 X-1 and IC10 X-1 are currently the only two robust extragalactic candidates for being Wolf-Rayet/black hole X-ray binaries, the Galactic analogue being Cyg X-3. These systems are believed to be a late product of high-mass X-ray binary evolution and direct progenitors of black hole mergers. From the analysis of Swift data, the orbital period of NGC 300 X-1 was found to be 32.8h. We here merge the full set of existing data of NGC300 X-1, using XMM-Newton, Chandra and Swift observations to derive a more precise value of the orbital period of 32.7932+-0.0029h above a confidence level of 99.99%. This allows us to phase connect the X-ray light curve of the source with radial velocity measurements of He II lines performed in 2010. We show that, as for IC10 X-1 and Cyg X-3, the X-ray eclipse corresponds to maximum of the blueshift of the He II lines, instead of the expected zero velocity. This indicates that for NGC300 X-1 as well, the wind of the WR star is completely ionised by the black hole radiation and that the emission lines come from the region of the WR star that is in the shadow. We also present for the first time the light curve of two recent very long XMM-Newton observations of the source, performed on the 16th to 20th of December 2016., Comment: 7 pages, 4 figures, from IAU Symposium 346: High-mass X-ray binaries: illuminating the passage from massive binaries to merging compact objects