Search

Outflows are critical components of many astrophysical systems, including accreting compact binaries and active galactic nuclei (AGN). These outflows can significantly affect a system's evolution and alter its observational appearance by reprocessing the radiation produced by the central engine. Sirocco (Simulating Ionization and Radiation in Outflows Created by Compact Objects - or "the code formerly known as Python") is a Sobolev-based Monte Carlo ionization and radiative transfer code. It is designed to simulate the spectra produced by any system with an azimuthally-symmetric outflow, from spherical stellar winds to rotating, biconical accretion disc winds. Wind models can either be parametrized or imported, e.g. from hydrodynamical simulations. The radiation sources include an optically thick accretion disc and various central sources with flexible spectra and geometries. The code tracks the "photon packets" produced by the sources in any given simulation as they traverse and interact with the wind. The code assumes radiative near-equilibrium, so the thermal and ionization state can be determined iteratively from these interactions. Once the physical properties in the wind have converged, Sirocco can be used to generate synthetic spectra at a series of observer sightlines. Here, we describe the physical assumptions, operation, performance and limitations of the code. We validate it against tardis, cmfgen and cloudy, finding good agreement, and present illustrative synthetic spectra from disc winds in cataclysmic variables, tidal disruption events, AGN and X-ray binaries. Sirocco is publicly available on GitHub, alongside its associated data, documentation and sample input files covering a wide range of astrophysical applications., Comment: 26 pages, 16 figures. Submitted to MNRAS. This is the release paper for the SIROCCO code, which can be found at https://github.com/sirocco-rt/sirocco with links to documentation. Underlying data and figure scripts available at https://github.com/sirocco-rt/release-models -- comments on the code, paper or documentation are welcomed

We have used the Dark Energy Camera (DECam) on the CTIO Blanco 4-m telescope to perform a new emission-line survey of the Large Magellanic Cloud (LMC) using narrow-band H-alpha and [SII] filters in addition to a continuum band for use in creating pure emission-line images. We refer to this new survey as DeMCELS, to distinguish it from the earlier Magellanic Cloud Emission Line Survey (MCELS) that has been in service for nearly 25 years. DeMCELS covers $\sim 54$ degrees$^{2}$, encompassing most of the bright optical disk of the LMC. With DECam's pixel size of only 0.27", our DeMCELS survey provides a seeing-limited improvement of 3-5 times over MCELS and is comparable in depth, with surface brightness limits of 3.3E-17 erg cm$^{-2}$ s$^{-1}$ arcsec$^{-2}$ in H-alpha and 2.9E-17 erg cm$^{-2}$ s$^{-1}$ arcsec$^{-2}$ in H-alpha and [SII], respectively. DeMCELS provides detailed morphological information on nebulae of all scales, from the largest supershells to individual [HII] regions and supernova remnants, to bubbles of emission surrounding individual stars, and even to faint structures in the diffuse ionized gas of the LMC. Many complex regions of emission show significant variations in the ratio of [SII] to H-alpha, a sign of the mixture of shocks from stellar winds and/or supernovae with photoionization by embedded hot, young stars. We present the details of the observing strategy and data processing for this survey, and show selected results in comparison with previous data. A companion project for the Small Magellanic Cloud is in progress and will be reported separately. We are making these new data available to the community at large via the NOIRLab's Data Lab site., Comment: 30 pages, 16 figures Accepted for publication (ApJ)

Tidal disruption events (TDEs) can generate non-spherical, relativistic and optically thick outflows. Simulations show that the radiation we observe is reprocessed by these outflows. According to a unified model suggested by these simulations, the spectral energy distributions (SEDs) of TDEs depend strongly on viewing angle: low [high] optical-to-X-ray ratios (OXRs) correspond to face-on [edge-on] orientations. Post-processing with radiative transfer codes have simulated the emergent spectra, but have so far been carried out only in a quasi-1D framework, with three atomic species (H, He and O). Here, we present 2.5D Monte Carlo radiative transfer simulations which model the emission from a non-spherical outflow, including a more comprehensive set of cosmically abundant species. While the basic trend of OXR increasing with inclination is preserved, the inherently multi-dimensional nature of photon transport through the non-spherical outflow significantly affects the emergent SEDs. Relaxing the quasi-1D approximation allows photons to preferentially escape in (polar) directions of lower optical depth, resulting in a greater variation of bolometric luminosity as a function of inclination. According to our simulations, inclination alone may not fully explain the large dynamic range of observed TDE OXRs. We also find that including metals, other than Oxygen, changes the emergent spectra significantly, resulting in stronger absorption and emission lines in the extreme ultraviolet, as well a greater variation in the OXR as a function of inclination. Whilst our results support previously proposed unified models for TDEs, they also highlight the critical importance of multi-dimensional ionization and radiative transfer., Comment: 17 pages, 15 figures. Submitted to MNRAS

SDSS-V is carrying out a dedicated survey for white dwarfs, single and in binaries, and we report the analysis of the spectroscopy of 504 cataclysmic variables (CVs) and CV candidates obtained during the first 34 months of observations of SDSS-V. We developed a convolutional neural network (CNN) to aid with the identification of CV candidates among the over 2 million SDSS-V spectra obtained with the BOSS spectrograph. The CNN reduced the number of spectra that required visual inspection to $\simeq2$ per cent of the total. We identified 776 CV spectra among the CNN-selected candidates, plus an additional 27 CV spectra that the CNN misclassified, but that were found serendipitously by human inspection of the data. Analysing the SDSS-V spectroscopy and ancillary data of the 504 CVs in our sample, we report 61 new CVs, spectroscopically confirm 248 and refute 13 published CV candidates, and we report 82 new or improved orbital periods. We discuss the completeness and possible selection biases of the machine learning methodology, as well as the effectiveness of targeting CV candidates within SDSS-V. Finally, we re-assess the space density of CVs, and find $1.2\times 10^{-5}\,\mathrm{pc^{-3}}$., Comment: Accepted by MNRAS. 24 pages, 25 figures. Supplementary data substantially updated from previous version and contains details of 504 CVs including spectra and light curves

Many high-state cataclysmic variables (CVs) exhibit blue-shifted absorption features in their ultraviolet (UV) spectra -- a smoking-gun signature of outflows. However, the impact of these outflows on {\em optical} spectra remains much more uncertain. During its recent outburst, the eclipsing dwarf nova V455 And displayed strong optical emission lines whose cores were narrower than expected from a Keplerian disc. Here, we explore whether disc + wind models developed for matching UV observations of CVs can also account for these optical spectra. Importantly, V455~And was extremely bright at outburst maximum: the accretion rate implied by fitting the optical continuum with a standard disc model is $\dot{M}_{\rm acc} \simeq 10^{-7}~{\rm M}_\odot~{\rm yr^{-1}}$. Allowing for continuum reprocessing in the outflow helps to relax this constraint. A disk wind can also broadly reproduce the optical emission lines, but only if the wind is (i) highly mass-loaded, with a mass-loss rate reaching $\dot{M}_{\rm wind} \simeq 0.4 \dot{M}_{\rm acc}$, and/or (ii) clumpy, with a volume filling factor $f_V \simeq 0.1$. The same models can describe the spectral evolution across the outburst, simply by lowering $\dot{M}_{\rm acc}$ and $\dot{M}_{\rm wind}$. Extending these models to lower inclinations and into the UV produces spectra consistent with those observed in face-on high-state CVs. We also find, for the first time in simulations of this type, P-Cygni-like absorption features in the Balmer series, as have been observed in both CVs and X-ray binaries. Overall, dense disc winds provide a promising framework for explaining multiple observational signatures seen in high-state CVs, but theoretical challenges persist., Comment: 13 pages, 9 figures, 3 tables. Accepted for publication in MNRAS

We present the Sloan Digital Sky Survey V (SDSS-V) Local Volume Mapper (LVM). The LVM is an integral-field spectroscopic survey of the Milky Way, Magellanic Clouds, and of a sample of local volume galaxies, connecting resolved pc-scale individual sources of feedback to kpc-scale ionized interstellar medium (ISM) properties. The 4-year survey covers the southern Milky Way disk at spatial resolutions of 0.05 to 1 pc, the Magellanic Clouds at 10 pc resolution, and nearby large galaxies at larger scales totaling $>4300$ square degrees of sky, and more than 55M spectra. It utilizes a new facility of alt-alt mounted siderostats feeding 16 cm refractive telescopes, lenslet-coupled fiber-optics, and spectrographs covering 3600-9800A at R ~ 4000. The ultra-wide field IFU has a diameter of 0.5 degrees with 1801 hexagonally packed fibers of 35.3 arcsec apertures. The siderostats allow for a completely stationary fiber system, avoiding instability of the line spread function seen in traditional fiber feeds. Scientifically, LVM resolves the regions where energy, momentum, and chemical elements are injected into the ISM at the scale of gas clouds, while simultaneously charting where energy is being dissipated (via cooling, shocks, turbulence, bulk flows, etc.) to global scales. This combined local and global view enables us to constrain physical processes regulating how stellar feedback operates and couples to galactic kinematics and disk-scale structures, such as the bar and spiral arms, as well as gas in- and out-flows., Comment: 29 pages, 12 figures, accepted for publication in The Astronomical Journal

Disc winds are a common feature in accreting astrophysical systems on all scales. In active galactic nuclei (AGN) and accreting white dwarfs (AWDs), specifically, radiation pressure mediated by spectral lines is a promising mechanism for driving these outflows. Previous hydrodynamical simulations have largely supported this idea, but relied on highly approximate treatments of ionization and radiative transfer. Given the sensitivity of line driving to the ionization state and radiation field in the outflow, here we present a new method for carrying out 2.5D radiation-hydrodynamic simulations that takes full account of the frequency-dependent radiative transfer through the wind, the corresponding ionization state and the resulting radiative accelerations. Applying our method to AWDs, we find that it is much harder to drive a powerful line-driven outflow when the interaction between matter and radiation is treated self-consistently. This conclusion is robust to changes in the adopted system parameters. The fundamental difficulty is that discs luminous enough to drive such a wind are also hot enough to over-ionize it. As a result, the mass-loss rates in our simulations are much lower than those found in earlier, more approximate calculations. We also show that the ultraviolet spectra produced by our simulations do not match those observed in AWDs. We conclude that, unless the over-ionization problem can be mitigated (e.g. by sub-grid clumping or a softer-than-expected radiation field), line driving may not be a promising mechanism for powering the outflows from AWDs. These conclusions are likely to have significant implications for disc winds in AGN also., Comment: Accepted for publication in MNRAS. 14 pages, 10 figures + 3 figures in Appendix

The nearby irregular galaxy NGC4449 has a star formation rate of about 0.4 solar masses/yr and should host of order 70 SNRs younger than 20,000 years, a typical age for SNRs expanding into to an ISM with unit density to reach the radiative phase. We have carried out an optical imaging and spectroscopic survey in an attempt to identify these SNRs. This task is challenging because diffuse gas with elevated ratios of [SII]:H-alpha is omnipresent in NGC4449, {causing confusion when using this common diagnostic for SNRs. Using narrow-band interference-filter images, we first identified 49 objects that have elevated [SII]:H-alpha ratios compared to nearby HII regions. Using Gemini-N and GMOS, we then obtained high-resolution spectra of 30 of these SNR candidates, 25 of which have [SII]:H-alpha ratios greater than 0.5. Of these, 15 nebulae are almost certainly SNRs, based on a combination of characteristics: higher [OI]:H-alpha ratios and broader line widths than observed from HII regions. The remainder are good candidates as well, but need additional confirmation. Surprisingly, despite having superior imaging and spectroscopic data sets to examine, we are unable to confirm most of the candidates suggested by Leonidaki (2013). While NGC4449 is likely an extreme case because of the high surface brightness and elevated [SII]:H-alpha ratio of diffuse gas, it highlights the need for sensitive high-resolution optical spectroscopy, or high spatial resolution radio or X-ray observations that can ensure accurate SNR identifications in external galaxies.

Blueshifts - or, more accurately, blue asymmetries - in broad emission lines such as CIV $\lambda$1550 are common in luminous quasars and correlate with fundamental properties such as Eddington ratio and broad absorption line (BAL) characteristics. However, the formation of these blueshifts is still not understood, and neither is their physical connection to the BAL phenomenon or accretion disc. In this work, we present Monte Carlo radiative transfer and photoionization simulations using parametrized biconical disc-wind models. We take advantage of the azimuthal symmetry of a quasar and show that we can reproduce CIV blueshifts provided that (i) the disc-midplane is optically thick out to radii beyond the line formation region, so that the receding wind bicone is obscured; and (ii) the system is viewed from relatively low (that is, more face-on) inclinations ($\lesssim40^\circ$). We show that CIV emission line blueshifts and BALs can form in the same wind structure. The velocity profile of the wind has a significant impact on the location of the line formation region and the resulting line profile, suggesting that the shape of the emission lines can be used as a probe of wind-driving physics. While we are successful at producing blueshifts/blue asymmetries in outflows, we struggle to match the detailed shape or skew of the observed emission line profiles. In addition, our models produce redshifted emission-line asymmetries for certain viewing angles. We discuss our work in the context of the CIV $\lambda$1550 emission blueshift versus equivalent-width space and explore the implications for quasar disc wind physics., Comment: 20 pages, 12 figures + 1 appendix. Supplementary material, simulation data and plotting scripts available at: https://github.com/jhmatthews/blueshifts

Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit., Comment: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

In order to better characterize the rich supernova remnant (SNR) population of M83 (NGC 5236), we have obtained high-resolution (about 85 km/s) spectra of 119 of the SNRs and SNR candidates in M83 with Gemini/GMOS, as well as new spectra of the young SNRs B12-174a and SN1957D. Most of the SNRs and SNR candidates have [S II]:H{\alpha} ratios that exceed 0.4. Combining these results with earlier studies we have carried out with MUSE and at lower spectroscopic resolution with GMOS, we have confirmed a total of 238 emission nebulae to be SNRs on the basis of their [S II]:H{\alpha} ratios, about half of which have emission lines that show velocity broadening greater than 100 km/s, providing a kinematic confirmation that they are SNRs and not H II regions. Looking at the entire sample, we find a strong correlation between velocity widths and the line ratios of [O I]{\lambda}6300:H{\alpha}, [N II]{\lambda}6584:H{\alpha} and [S II]{\lambda}{\lambda}6716,6731:H{\alpha}. The density-sensitive [S II]{\lambda}6716:{\lambda}6731 line ratio is strongly correlated with SNR diameter, but not with the velocity width. We discuss these results in the context of previously published shock models., Comment: 39 pages, including 9 figures (some with multiple sub-figures) and 6 tables

We study pre-supernova feedback in a sample of $\sim$ 4700 HII regions in the nearby spiral galaxy M83, identified on their H$\alpha$ emission. We pectroscopically identify Wolf-Rayet (WR) stars populating the star-forming regions. For each HII region, we compute the pressure of ionised gas ($P_{\rm ion}$) and the direct radiation pressure ($P_{\rm dir}$) acting in the region, and investigate how they vary with galactocentric distance, with the physical properties of the region, and with the pressure of the galactic environment ($P_\mathrm{DE}$). For a subset of $\sim$ 500 regions, we also investigate the link between the pressure terms and the properties of the cluster population (age, mass, and LyC flux). We find that $P_{\rm ion}$ dominates over $P_{\rm dir}$ by at least a factor of 10 on average over the disk. Both pressure terms are strongly enhanced and become almost comparable in the central starburst region. In the disk ($R \geq 0.15~R_e$), we observe that $P_{\rm dir}$ stays approximately constant with galactocentric distance. $P_{\rm dir}$ is positively correlated with an increase in radiation field strength (linked to the negative metallicity gradient in the galaxy), while it decreases in low extinction regions. $P_{\rm ion}$ decreases constantly for increasing galactocentric distances. In general, we observe that HII regions near the center are underpressured with respect to their surroundings, whereas regions in the disk are overpressured and hence expanding. We find that regions hosting younger clusters or having more mass in young star clusters have a higher internal pressure, indicating that clustered star formation is likely playing a dominant role in setting the pressure. Finally, we estimate that only 13 % of HII regions hosting young clusters and WR stars have $f_{\rm esc} \geq 0$.[Abridged], Comment: Accepted for publication in A&A

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) on board the Hubble Space Telescope (HST). We compare the catalogue of FUV objects to ten known X-ray sources and find six confident matches of two cataclysmic variables (CVs), one RS CVn, one red giant with strong FUV emission and two sources only detected in the FUV. We also searched for variable sources in our dataset and found a total of seven blue stragglers (BSs), four horizontal branch (HB) stars, five red giant branch stars, 28 main sequence stars and four gap objects that demonstrated variability. One BS star is a known W-UMa contact binary, one of the gap objects is a known CV identified in this work to be a dwarf nova, and the three other gap sources are weak variables. The periods and positions of two of the variable HB stars match them to two previously known RR Lyrae variables of types RRab and RRc., Comment: 13 pages, 12 figures, accepted for publication to MNRAS

Here we describe a new study of the SNRs and SNR candidates in nearby face-on spiral galaxy M83, based primarily on MUSE integral field spectroscopy. Our revised catalog of SNR candidates in M83 has 366 objects, 81 of which are reported here for the first time. Of these, 229 lie within the MUSE observation region, 160 of which have spectra with [SII]:Halpha ratios exceeding 0.4, the value generally accepted as confirmation that an emission nebula is shock-heated. Combined with 51 SNR candidates outside the MUSE region with high [SII]:Halpha ratios, there are 211 spectroscopically-confirmed SNRs in M83, the largest number of confirmed SNRs in any external galaxy. MUSE's combination of relatively high spectral resolution and broad wavelength coverage has allowed us to explore two other properties of SNRs that could serve as the basis of future SNR searches. Specifically, most of the objects identified as SNRs on the basis of [SII]:Halpha ratios exhibit more velocity broadening and lower ratios of [SIII]:[SII] emission than HII regions. A search for nebulae with the very broad emission lines expected from young, rapidly expanding remnants revealed none, except for the previously identified B12-174a. The SNRs identified in M83 are, with few exceptions, middle-aged ISM-dominated ones. Smaller diameter candidates show a larger range of velocity broadening and a larger range of gas densities than the larger diameter objects, as expected if the SNRs expanding into denser gas brighten and then fade from view at smaller diameters than those expanding into a more tenuous ISM, Comment: 28 page, including 14 figures and 4 tables; accepted for publication in ApJ

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) which were both on board the Hubble Space Telescope (HST). The FUV-UV colour-magnitude diagram (CMD) shows a main sequence (MS) turnoff at FUV $\approx$ 22 mag and FUV-UV $\approx$ 3 mag. The MS extends 4 mag below the turnoff, and a prominent horizontal branch (HB) and blue straggler (BS) sequence can be seen. A total of 1218 MS stars, 185 red giant branch stars, 47 BS stars and 41 HB stars are identified, along with 78 sources blueward of the MS which consist of white dwarfs (WDs) and objects in the gap between the WDs and the MS that include potential cataclysmic variable (CV) candidates. The radial distribution of the BS population is concentrated towards the cluster centre, indicating that mass segregation has occurred. The blue and red sub-populations of the double BS sequence appear mixed in the ultraviolet CMD, and no significant central concentration of CV candidates is seen in this cluster., Comment: 11 pages, 7 figures. Accepted for publication in MNRAS

We present results from deep Chandra observations of the young Type Ia supernova remnant (SNR) 0509-68.7, also known as N103B, located in the Large Magellanic cloud (LMC). The remnant displays an asymmetry in brightness, with the western hemisphere appearing significantly brighter than the eastern half. Previous multi-wavelength observations have attributed the difference to a density gradient and suggested circumstellar material origins, drawing similarities to Kepler's SNR. We apply a clustering technique combined with traditional imaging analysis to spatially locate various emission components within the remnant. We find that O and Mg emission is strongest along the blast wave, and coincides with Spitzer observations of dust emission and optical emission from the non-radiative shocks. The abundances of O and Mg in these regions are enhanced relative to the average LMC abundances and appear as a distinct spatial distribution compared to the ejecta products, supporting the circumstellar medium (CSM) interpretation. We also find that the spatial distribution of Cr is identical to that of Fe in the interior of the remnant, and does not coincide at all with the O and Mg emission., Comment: 23 pages, 17 Figures. Accepted for publication in ApJ

A significant number of tidal disruption events (TDEs) radiate primarily at optical and ultraviolet (UV) wavelengths, with only weak soft X-ray components. One model for this optical excess proposes that thermal X-ray emission from a compact accretion disc is reprocessed to longer wavelengths by an optically thick envelope. Here, we explore this reprocessing scenario in the context of an optically thick accretion disc wind. Using state-of-the-art Monte Carlo radiative transfer and ionization software, we produce synthetic UV and optical spectra for wind and disc-hosting TDEs. Our models are inspired by observations, spanning a realistic range of accretion rates and wind kinematics. We find that such outflows can efficiently reprocess the disc emission and produce the broad Balmer and helium recombination features commonly seen in TDEs and exhibit asymmetric red wings. Moreover, the characteristic colour temperature of the reprocessed spectral energy distribution (SED) is much lower than that of the accretion disc. We show explicitly how changes in black hole mass, accretion rate and wind properties affect the observed broadband SED and line spectrum. In general, slower, denser winds tend to reprocess more radiation and produce stronger Balmer emission. Most of the outflows we consider are too highly ionized to produce UV absorption features, but this is sensitive to the input SED. For example, truncating the inner disc at just 4 $R_{ISCO}$ lowers the wind ionization state sufficiently to produce UV absorption features for sight lines looking into the wind, Comment: Accepted to MNRAS. 19 pages, 11 figures

Particle acceleration to suprathermal energies in strong astrophysical shock waves is a widespread phenomenon, generally explained by diffusive shock acceleration. Such shocks can also amplify upstream magnetic field considerably beyond simple compression. The complex plasma physics processes involved are often parameterized by assuming that shocks put some fraction $\epsilon_e$ of their energy into fast particles, and another fraction $\epsilon_B$ into magnetic field. Modelers of shocks in supernovae, supernova remnants, and gamma-ray bursters, among other locations, often assume typical values for these fractions, presumed to remain constant in time. However, it is rare that enough properties of a source are independently constrained that values of the epsilons can be inferred directly. Supernova remnants (SNRs) can provide such circumstances. Here we summarize results from global fits to spatially integrated emission in six young SNRs, finding $10^{-4} \le \epsilon_e \le 0.05$ and $0.001 \le \epsilon_B \le 0.1$. These large variations might be put down to the differing ages and environments of these SNRs, so we conduct a detailed analysis of a single remnant, that of Kepler's supernova. Both epsilons can be determined at seven different locations around the shock, and we find even larger ranges for both epsilons, as well as for their ratio (thus independent of the shock energy itself). We conclude that unknown factors have a large influence on the efficiency of both processes. Shock obliquity, upstream neutral fraction, or other possibilities need to be explored, while calculations assuming fixed values of the epsilons should be regarded as provisional., Comment: 16 pages, 7 figures; accepted for publication in ApJ

We constrained the progenitor masses for 169 supernova remnants, 8 historically observed supernovae, and the black hole formation candidate in NGC 6946, finding that they are consistent with originating from a standard initial mass function. Additionally, there were 16 remnants that showed no sign of nearby star formation consistent with a core-collapse supernova, making them good Type Ia candidates. Using $Hubble$ $Space$ $Telescope$ broadband imaging, we measured stellar photometry of ACS/WFC fields in F435W, F555W, F606W, and F814W filters as well as WFC3/UVIS fields in F438W, F606W, and F814W. We then fitted this photometry with stellar evolutionary models to determine the ages of the young populations present at the positions of the SNRs and SNe. We then infer a progenitor mass probability distribution from the fitted age distribution. For 37 SNRs we tested how different filter combinations affected the inferred masses. We find that filters sensitive to H$\alpha$, [N II], and [S II] gas emission can bias mass estimates for remnants that rely on our technique. Using a KS-test analysis on our most reliable measurements, we find the progenitor mass distribution is well-matched by a power-law index of $-2.6^{+0.5}_{-0.6}$, which is consistent with a standard initial mass function., Comment: 28 pages, 13 figures, 4 tables

The spectra arising from the disks of nova-like variables show many of the features seen in stellar atmospheres. They are typically modelled either from an appropriated weighted set of stellar atmospheres or a disk atmosphere with energy is dissipated near the disk plane, with the effective temperature distribution expected from a steady state accretion disk. However these models generally over-predict the depth of the Balmer jump and the slope of the spectrum in the ultraviolet. The problem is likely due to energy dissipation in the disk atmosphere, which produces a flatter vertical temperature profile than is observed in stars. Here, we provide validation for this hypothesis in the form of spectra generated using the stellar atmosphere code TLUSTY using a parametric prescription for energy dissipation as a function of depth and closely match the spectrum of the nova-like IX Vel over the wavelength range 1150-6000 \AA., Comment: 12 pages, 7 figures

Using a combination of ground-based and HST imaging, we have constructed a catalog of 179 supernova remnants (SNRs) and SNR candidates in the nearby spiral galaxy M51. Follow-up spectroscopy of 66 of the candidates confirms 61 of these as SNRs, and suggests that the vast majority of the unobserved objects are SNRs as well. A total of 55 of the candidates are coincident with (mostly soft) X-ray sources identified in deep Chandra observations of M51; searching the positions of other soft X-ray sources resulted in several additional possible optical candidates. There are 16 objects in the catalog coincident with known radio sources. None of the sources with spectra shows the high velocities (>500 km/s) characteristic of young, ejecta-dominated SNRs like Cas A; instead, most if not all appear to be middle-aged SNRs. The general properties of the SNRs, size distribution and spectral characteristics, resemble those in other nearby spiral galaxies, notably M33, M83, and NGC6946, where similar samples exist. However, the spectroscopically observed [N II]:H{\alpha} ratios appear to be significantly higher than in any of these other galaxies. Although we have explored various ideas to explain the high ratios in M51, none of the explanations appears to be satisfactory., Comment: 44 pages, 17 figures, and 7 tables To be published in The Astrophysical Journal

We carry out a comparative analysis of the metallicities from the stellar, neutral-gas, and ionized-gas components in the metal-rich spiral galaxy M83. We analyze spectroscopic observations taken with the Hubble Space Telescope (HST), the Large Binocular Telescope (LBT) and the Very Large Telescope (VLT). We detect a clear depletion of the HI gas, as observed from the HI column densities in the nuclear region of this spiral galaxy. We find column densities of log[$N$(HI) cm$^{-2}$] $<$ 20.0 at galactocentric distances of $<$ 0.18 kpc, in contrast to column densities of log[$N$(HI) cm$^{-2}$] $\sim$ 21.0 in the galactic disk, a trend observed in other nearby spiral galaxies. We measure a metallicity gradient of $-$0.03 $\pm$ 0.01 dex kpc$^{-1}$ for the ionized gas, comparable to the metallicity gradient of a local benchmark of 49 nearby star-forming galaxies of $-$0.026 $\pm$ 0.002 dex kpc$^{-1}$. Our co-spatial metallicity comparison of the multi-phase gas and stellar populations shows excellent agreement outside of the nucleus of the galaxy hinting at a scenario where the mixing of newly synthesized metals from the most massive stars in the star clusters takes longer than their lifetimes ($\sim$10 Myr). Finally, our work shows that caution must be taken when studying the metallicity gradient of the neutral-gas component in star-forming galaxies, since this can be strongly biased, as these environments can be dominated by molecular gas. In these regions the typical metallicity tracers can provide inaccurate abundances as they may trace both the neutral- and molecular-gas components., Comment: 29 pages, 14 tables, 9 figures, Accepted for Publication in ApJ

We present HST/WFC3 multiband imagery of N103B, the remnant of a Type Ia supernova in the Large Magellanic Cloud, as well as HST/COS ultraviolet spectroscopy of the brightest radiatively shocked region. The images show a wide range of morphology and relative emission-line intensities, from smooth Balmer-line dominated collisionless shocks due to the primary blast wave, to clumpy radiative shock filaments due to secondary shocks in density enhancements. The COS data show strong FUV line emission despite a moderately high extinction along this line of sight. We use the COS data with previous optical spectra to constrain the shock conditions and refine the abundance analysis, finding abundances typical of the local interstellar medium within the uncertainties. Under an assumption that the material being shocked was shed from the pre-supernova system, this finding places constraints on any significant enrichment in that material, and thus on the non-degenerate star in what was presumably a single-degenerate Type Ia supernova., Comment: 25 pages, 6 figures, 2 tables, accepted for Astrophysical Journal, August 27, 2020

NGC6946 is a high star formation rate face-on spiral galaxy that has hosted ten supernovae since 1917. Not surprisingly, a large number of supernova remnants and candidates have been identified either as optical nebulae with high [S II]:H$\alpha$ line ratios (147) or as compact non-thermal radio sources (35). However, there are only seven overlaps between these two samples. Here, we apply [Fe II] 1.644 $\mu$m emission as a new diagnostic to search for supernova remnants in an attempt to resolve this discrepancy. [Fe II] is expected to be relatively strong in the radiative shocks of supernova remnants and almost absent in HII regions. It is less susceptible to the effects of absorption along the line of sight than the optical lines normally used to identify remnants. Using data from the WFC3 camera on HST}, we identify 132 [Fe II] emission nebulae in NGC6946 as likely supernova remnants. Of these, 54 align with previously known optical supernova remnants. The remaining 78 objects are new; of these 44 are visible in new HST imagery in H$\alpha$ and [S II]. This brings the total number of supernova remnant candidates (from optical and/or IR data) in NGC6946 to 225. A total of 14 coincidences with radio supernova remnant candidates (out of 30 in our search area) are found in this expanded list. The identification of so many new remnant candidates validates the use of [Fe II] imagery for finding remnants, and suggests that previous remnant searches in other galaxies may be far from complete., Comment: 43 pages, 9 figures, 2 tables, accepted for publication in ApJ

We present a new catalogue of radio sources in the face-on spiral galaxy M83. Radio observations taken in 2011, 2015, and 2017 with the Australia Telescope Compact Array (ATCA) at 5.5 and 9 GHz have detected 270 radio sources. Although a small number of these sources are background extragalactic sources, most are either H II regions or supernova remnants (SNRs) within M83 itself. Three of the six historical supernovae are detected, as is the very young remnant that had been identified in a recent study, which is likely the result of a supernova that exploded in the last ~100 years but was missed. All of these objects are generally fading with time. Confusion limits our ability to measure the radio emission from a number of the SNRs in M83, but 64 were detected in unconfused regions, and these have the approximate power-law luminosity function which has been observed in other galaxies. The SNRs in M83 are systematically smaller in diameter and brighter than those that have been detected at radio wavelengths in M33. A number of the radio sources are coincident with X-ray sources in M83; most of these coincident sources turn out to be supernova remnants. Our dual frequency observations are among the most sensitive to date for a spiral galaxy outside the Local Group; despite this we were not able to place realistic constraints on the spectral indices, and as a result, it was not possible to search for supernova remnants based on their radio properties alone., Comment: 23 pages, 17 figures, 8 tables. Accepted for publication in MNRAS. The full machine-readable catalogues are provided in the source (MRT) folder

Some tidal disruption events (TDEs) exhibit blueshifted broad absorption lines (BALs) in their rest-frame ultraviolet (UV) spectra, while others display broad emission lines (BELs). Similar phenomenology is observed in quasars and accreting white dwarfs, where it can be interpreted as an orientation effect associated with line formation in an accretion disc wind.We propose and explore a similar unification scheme for TDEs. We present synthetic UV spectra for disc and wind-hosting TDEs, produced by a state-of-the-art Monte Carlo ionization and radiative transfer code. Our models cover a wide range of disc wind geometries and kinematics. Such winds naturally reproduce both BALs and BELs. In general, sight lines looking into the wind cone preferentially produce BALs, while other orientations preferentially produce BELs. We also study the effect of wind clumping and CNO-processed abundances on the observed spectra. Clumpy winds tend to produce stronger UV emission and absorption lines, because clumping increases both the emission measure and the abundances of the relevant ionic species, the latter by reducing the ionization state of the outflow. The main effect of adopting CNO-processed abundances is a weakening of C~{\sc iv}~1550~\AA~ and an enhancement of N \textsc{v}~1240~\AA~ in the spectra. We conclude that line formation in an accretion disc wind is a promising mechanism for explaining the diverse UV spectra of TDEs. If this is correct, the relative number of BAL and BEL TDEs can be used to estimate the covering factor of the outflow. The models in this work are publicly available online and upon request., Comment: 16 pages, 9 figures

We present a study of the eclipses in the accreting white dwarf EX Dra during TESS Cycles 14 and 15. During both of the two outbursts present in this dataset, the eclipses undergo a hysteretic loop in eclipse-depth/out-of-eclipse-flux space. In each case, the direction in which the loops are executed strongly suggests an outburst which is triggered near the inner edge of the accretion disk and propagates outwards. This in turn suggests that the outbursts in EX Dra are 'Inside Out' outbursts; events predicted by previous hydrodynamic studies of dwarf nova accretion disks and confirmed spectroscopically in a number of other accreting white dwarf systems. We therefore propose that the direction of the loop executed in eclipse-depth/out-of-eclipse flux space be used as a test to phenomenologically distinguish between 'inside out' and 'outside in' outbursts in other eclipsing dwarf novae; a reliable and purely photometric test to differentiate between these phenomena., Comment: 10 pages, 6 figures

X-ray signatures of outflowing gas have been detected in several accreting black-hole binaries, always in the soft state. A key question raised by these observations is whether these winds might also exist in the hard state. Here, we carry out the first full-frequency radiation hydrodynamic simulations of luminous ($\rm{L = 0.5 \, L_{\mathrm{Edd}}}$) black-hole X-ray binary systems in both the hard and the soft state, with realistic spectral energy distributions (SEDs). Our simulations are designed to describe X-ray transients near the peak of their outburst, just before and after the hard-to-soft state transition. At these luminosities, it is essential to include radiation driving, and we include not only electron scattering, but also photoelectric and line interactions. We find powerful outflows with $\rm{\dot{M}_{wind} \simeq 2 \,\dot{M}_{acc}}$ are driven by thermal and radiation pressure in both hard and soft states. The hard-state wind is significantly faster and carries approximately 20 times as much kinetic energy as the soft-state wind. However, in the hard state the wind is more ionized, and so weaker X-ray absorption lines are seen over a narrower range of viewing angles. Nevertheless, for inclinations $\gtrsim 80^{\circ}$, blue-shifted wind-formed Fe XXV and Fe XXVI features should be observable even in the hard state. Given that the data required to detect these lines currently exist for only a single system in a {\em luminous} hard state -- the peculiar GRS~1915+105 -- we urge the acquisition of new observations to test this prediction. The new generation of X-ray spectrometers should be able to resolve the velocity structure., Comment: Accepted for publication by MNRAS

The origin, geometry and kinematics of the broad line region (BLR) gas in quasars and active galactic nuclei (AGN) are uncertain. We demonstrate that clumpy biconical disc winds illuminated by an AGN continuum can produce BLR-like spectra. We first use a simple toy model to illustrate that disc winds make quite good BLR candidates, because they are self-shielded flows and can cover a large portion of the ionizing flux-density ($\phi_H$-$n_H$) plane. We then conduct Monte Carlo radiative transfer and photoionization calculations, which fully account for self-shielding and multiple scattering in a non-spherical geometry. The emergent model spectra show broad emission lines with equivalent widths and line ratios comparable to those observed in AGN, provided that the wind has a volume filling factor of $f_V\lesssim0.1$. Similar emission line spectra are produced for a variety of wind geometries (polar or equatorial) and for launch radii that differ by an order of magnitude. The line emission arises almost exclusively from plasma travelling below the escape velocity, implying that `failed winds' are important BLR candidates. The behaviour of a line-emitting wind (and possibly any `smooth flow' BLR model) is similar to that of the locally optimally-emitting cloud (LOC) model originally proposed by Baldwin et al (1995), except that the gradients in ionization state and temperature are large-scale and continuous, rather than within or between distinct clouds. Our models also produce UV absorption lines and X-ray absorption features, and the stratified ionization structure can partially explain the different classes of broad absorption line quasars., Comment: Accepted to MNRAS. 20 pages, 13 figures plus appendix. Supplementary plots, some simulation data and demo notebooks can be found at https://github.com/jhmatthews/windy-blr-2020

Microquasars are neutron star or black hole X-ray binaries with jets. These jets can create shock-ionized bubbles of hot plasma that can masquerade as peculiar supernova remnants (SNRs) in extragalactic surveys. To see if this is the case in the well-studied spiral galaxy M83, where one microquasar candidate (M83-MQ1) has already been identified, we studied the properties of nine SNR candidates, selected because of their elongated or irregular morphology, from the set of previously identified SNRs in that galaxy. Using multiwavelength data from Chandra, the Hubble Space Telescope, Gemini, and the Australia Telescope Compact Array, we found that at least six of our nine sources are best interpreted as SNRs. For one source, we do not have enough observational data to explain its linear morphology. Another source shows a nebular optical spectrum dominated by photo-ionization by O stars, but its excess [Fe II] and radio luminosity suggest a possible hidden SNR. Finally, one source (S2) shows an elongated structure of ionized gas, two radio sources along that line, and an accretion-powered X-ray source in between them (the Chandra source L14-139). While S2 could be a chance alignment of multiple SNRs and one X-ray binary, it seems more likely that it is a single physical structure powered by the jet from the accreting compact object. In the latter case, the size and luminosity of S2 suggest a kinetic power of around 10^{40} erg/s, in the same class as the most powerful microquasars in the local universe (e.g., S26 in NGC7793 and SS433 in our own Galaxy)., Comment: Accepted by ApJ. Size = 8 MB

Many more supernova remnants (SNRs) are now known in external galaxies than in the Milky Way. Most of these SNRs have been identified using narrow-band imaging, separating SNRs from HII regions on the basis of [SII]:H-alpha ratios that are elevated compared to HII regions. However, the boundary between SNRs and HII regions is not always distinct, especially at low surface brightness. Here we explore velocity structure as a possible criterion for separating SNRs from HII regions, using a sample of well-studied SNRs in the Large Magellanic Cloud (LMC) as well as a small number of SNRs in the galaxy M83. We find, perhaps not surprisingly, that even at large diameters, SNRs exhibit velocity broadening sufficient to readily distinguish them from HII regions. We thus suggest that the purity of most extragalactic samples would be greatly improved through spectroscopic observations with a velocity resolution of order 50~km/s$., Comment: 19 pages, 9 figures. To be published in ApJ

We determine the ages of the young, resolved stellar populations at the locations of 237 optically-identified supernova remnants in M83. These age distributions put constraints on the progenitor masses of the supernovae that produced 199 of the remnants. The other 38 show no evidence for having a young progenitor and are therefore good Type Ia SNR candidates. Starting from Hubble Space Telescope broadband imaging, we measured resolved stellar photometry of seven archival WFC3/UVIS fields in F336W, F438W, and F814W. We generate color-magnitude diagrams of the stars within 50 pc of each SNR and fit them with stellar evolution models to obtain the population ages. From these ages we infer the progenitor mass that corresponds to the lifetime of the most prominent age that is $<$50 Myr. In this sample, there are 47 SNRs with best-fit progenitor masses $>$15 M$_{\odot}$, and 5 of these are $>$15 M$_{\odot}$ at 84% confidence. This is the largest collection of high-mass progenitors to date, including our highest-mass progenitor inference found so far, with a constraint of $<$8 Myr. Overall, the distribution of progenitor masses has a power-law index of $-3.0^{+0.2}_{-0.7}$, steeper than Salpeter initial mass function ($-2.35$). It remains unclear whether the reason for the low number of high-mass progenitors is due to the difficulty of finding and measuring such objects or because only a fraction of very massive stars produce supernovae., Comment: Accepted for publication in ApJ, 33 pages, 7 figures, 3 tables

We have performed new 1.4 GHz and 5 GHz observations of the Local Group galaxy M33 with the Jansky Very Large Array. Our survey has a limiting sensitivity of 20 uJy (4-sigma) and a resolution of 5.9 arcsec (FWHM), corresponding to a spatial resolution of 24 pc at 817 kpc. Using a new multi-resolution algorithm, we have created a catalog of 2875 sources, including 675 with well-determined spectral indices. We detect sources at the position of 319 of the X-ray sources in the Tuellmann et al. (2011) Chandra survey of M33, the majority of which are likely to be background galaxies. The radio source coincident with M33 X-8, the nuclear source, appears to be extended. Along with numerous H II regions or portions of H II region complexes, we detect 155 of the 217 optical supernova remnants included in the lists of Long et al. (2010) and Lee & Lee (2014), making this by far the largest sample of remnants at known distances with multiwavelength coverage. The remnants show a large dispersion in the ratio of radio to X-ray luminosity at a given diameter, a result that challenges the current generation of models for synchrotron radiation evolution in supernova remnants. See http://sundog.stsci.edu/m33 for access to catalogs and images., Comment: 37 pages, 24 figures; Accepted for publication in the Astrophysical Journal; see http://sundog.stsci.edu/m33 for access to catalogs and images

We present extensive optical photometric and spectroscopic observations, from 4 to 482 days after explosion, of the Type II-plateau (II-P) supernova (SN) 2017eaw in NGC 6946. SN 2017eaw is a normal SN II-P intermediate in properties between, for example, SN 1999em and SN 2012aw and the more luminous SN 2004et, also in NGC 6946. We have determined that the extinction to SN 2017eaw is primarily due to the Galactic foreground and that the SN site metallicity is likely subsolar. We have also independently confirmed a tip-of-the-red-giant-branch (TRGB) distance to NGC 6946 of 7.73+/-0.78 Mpc. The distances to the SN that we have also estimated via both the standardized candle method and expanding photosphere method corroborate the TRGB distance. We confirm the SN progenitor identity in pre-explosion archival Hubble Space Telescope (HST) and Spitzer Space Telescope images, via imaging of the SN through our HST Target of Opportunity program. Detailed modeling of the progenitor's spectral energy distribution indicates that the star was a dusty, luminous red supergiant consistent with an initial mass of ~15 Msuns., Comment: Appearing in ApJ

The relatively nearby spiral galaxy NGC~6946 is one of the most actively star forming galaxies in the local Universe. Ten supernovae (SNe) have been observed since 1917, and hence NGC6946 surely contains a large number of supernova remnants (SNRs). Here we report a new optical search for these SNRs using narrow-band images obtained with the WIYN telescope. We identify 147 emission nebulae as likely SNRs, based on elevated [SII]:Halpha ratios compared to HII regions. We have obtained spectra of 102 of these nebulae with Gemini North-GMOS; of these, 89 have [SII]:Halpha ratios greater than 0.4, the canonical optical criterion for identifying SNRs. There is very little overlap between our sample and the SNR candidates identified by Lacey et al. (2001) from radio data. Also, very few of our SNR candidates are known X-ray sources, unlike the situation in some other galaxies such as M33 and M83. The emission line ratios, e.g., [NII]:Halpha, of the candidates in NGC6946 are typical of those observed in SNR samples from other galaxies with comparable metallicity. None of the candidates observed in our low-resolution spectra show evidence of anomalous abundances or significant velocity broadening. A search for emission at the sites of all the historical SNe in NGC6946 resulted in detections for only two: SN1980K and SN2004et. Spectra of both show very broad, asymmetric line profiles, consistent with the interaction between SN ejecta and the progenitor star's circumstellar material, as seen in late spectra from other core-collapse SNe of similar age., Comment: 36 pages, 11 figures, accepted for publication in Astrophysical Journal

We have carried out radiation-hydrodynamic simulations of thermally-driven accretion disc winds in low-mass X-ray binaries. Our main goal is to study the luminosity dependence of these outflows and compare with observations. The simulations span the range $\rm{0.04 \leq L_{acc}/L_{Edd} \leq 1.0}$ and therefore cover most of the parameter space in which disc winds have been observed. Using a detailed Monte Carlo treatment of ionization and radiative transfer, we confirm two key results found in earlier simulations that were carried out in the optically thin limit: (i) the wind velocity -- and hence the maximum blueshift seen in wind-formed absorption lines -- increases with luminosity; (ii) the large-scale wind geometry is quasi-spherical, but observable absorption features are preferentially produced along high-column equatorial sightlines. In addition, we find that (iii) the wind efficiency always remains approximately constant at $\rm{\dot{M}_{wind}/\dot{M}_{acc} \simeq 2}$, a behaviour that is consistent with observations. We also present synthetic Fe XXV and Fe XXVI absorption line profiles for our simulated disc winds in order to illustrate the observational implications of our results., Comment: Accepted for publication in MNRAS

Stellar populations are powerful tools for investigating the evolution of extragalactic environments. We present the first UV integrated-light spectroscopic observations for 15 young star clusters in the starburst M83 with a special focus on metallicity measurements. The data were obtained with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope. We analyse the data applying an abundance technique previously used to study an optical set of star clusters. We estimate a central metallicity of [Z] = $+$0.20 $\pm$ 0.15 dex in agreement with those obtained through independent methods, i.e. $J$-band and blue supergiants. We estimate a UV metallicity gradient of $-$0.041 $\pm$ 0.022 dex kpc$^{-1}$ consistent with the optical metallicity gradient of $-$0.040 $\pm$ 0.032 dex kpc$^{-1}$ for $R/R_{25}<0.5$. Combining our stellar metallicities, UV and optical, with those from HII regions (strong-line abundances based on empirical calibrations) we identify two possible breaks in the gradient of M83 at galactocentric distances of $R\sim0.5$ and $1.0\:R_{25}$. If the abundance breaks are genuine, the metallicity gradient of this galaxy follows a steep-shallow-steep trend, a scenario predicted by three-dimensional (3D) numerical simulations of disc galaxies. The first break is located near the corotation radius. This first steep gradient may have originated by recent star formation episodes and a relatively young bar ($<$1 Gyr). In the numerical simulations the shallow gradient is created by the effects of dilution by outflow where low-metallicity material is mixed with enriched gas. And finally, the second break and last steep gradient mark the farthest galactocentric distances where the outward flow has penetrated., Comment: 18 pages, 7 figures, Accepted for Publication in ApJ

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

We present disc wind models aimed at reproducing the main features of the strong Na I resonance line P-Cygni profiles in the rapidly-accreting pre-main sequence FU Ori objects. We conducted Monte Carlo radiative transfer simulations for a standard magnetocentrifugally driven wind (MHD) model and our own "Genwind" models, which allows for a more flexible wind parameterisation. We find that the fiducial MHD wind and similar Genwind models, which have flows emerging outward from the inner disc edge, and thus have polar cavities with no absorbing gas, cannot reproduce the deep, wide Na I absorption lines in FU Ori objects viewed at low inclination. We find that it is necessary to include an "inner wind" to fill this polar cavity to reproduce observations. In addition, our models assuming pure scattering source functions in the Sobolev approximation at intermediate viewing angles ($30^{\circ} \lesssim i \lesssim 60^{\circ}$) do not yield sufficiently deep line profiles. Assuming complete absorption yields better agreement with observations, but simple estimates strongly suggest that pure scattering should be a much better approximation. The discrepancy may indicate that the Sobolev approximation is not applicable, possibly due to turbulence or non-monotonic velocity fields; there is some observational evidence for the latter. Our results provide guidance for future attempts to constrain FU Ori wind properties using full MHD wind simulations, by pointing to the importance of the boundary conditions necessary to give rise to an inner wind, and by suggesting that the winds must be turbulent to produce sufficiently deep line profiles., Comment: 12 pages, 17 figures, accepted for publication in MNRAS

We present new Hubble Space Telescope (HST) ultraviolet and ground-based optical observations of the hot, metal-rich white dwarf GD 394. Extreme-ultraviolet (EUV) observations in 1992-1996 revealed a 1.15d periodicity with a 25 percent amplitude, hypothesised to be due to metals in a surface accretion spot. We obtained phase-resolved HST/Space Telescope Imaging Spectrograph (STIS) high-resolution far-ultraviolet (FUV) spectra of GD 394 that sample the entire period, along with a large body of supplementary data. We find no evidence for an accretion spot, with the flux, accretion rate and radial velocity of GD 394 constant over the observed timescales at ultraviolet and optical wavelengths. We speculate that the spot may have no longer been present when our observations were obtained, or that the EUV variability is being caused by an otherwise undetected evaporating planet. The atmospheric parameters obtained from separate fits to optical and ultraviolet spectra are inconsistent, as is found for multiple hot white dwarfs. We also detect non-photospheric, high-excitation absorption lines of multiple volatile elements, which could be evidence for a hot plasma cocoon surrounding the white dwarf., Comment: Accepted for publication in MNRAS

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$., Comment: 19 pages, 18 figures, 6 tables. Accepted for publication in PASJ

We present a second epoch of {\it Chandra} observations of the Type Ia LMC SNR 0509-68.7 (N103B) obtained in 2017. When combined with the earlier observations from 1999, we have a 17.4-year baseline with which we can search for evidence of the remnant's expansion. Although the lack of strong point source detections makes absolute image alignment at the necessary accuracy impossible, we can measure the change in the diameter and the area of the remnant, and find that it has expanded by an average velocity of 4170 (2860, 5450) km s$^{-1}$. This supports the picture of this being a young remnant; this expansion velocity corresponds to an undecelerated age of 850 yr, making the real age somewhat younger, consistent with results from light echo studies. Previous infrared observations have revealed high densities in the western half of the remnant, likely from circumstellar material, so it is likely that the real expansion velocity is lower on that side of the remnant and higher on the eastern side. A similar scenario is seen in Kepler's SNR. N103B joins the rare class of Magellanic Cloud SNRs with measured proper motions., Comment: Accepted for publication in ApJ Letters

The X-ray source populations within galaxies are typically difficult to identify and classify from X-ray data alone. We are able to break through this barrier by combining deep new Chandra ACIS-I observations with extensive Hubble Space Telescope imaging from the PHAT survey of the M31 disk. We detect 373 X-ray sources down to 0.35-8.0 keV flux of 10$^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 0.4 square degrees, 170 of which are reported for the first time. We identify optical counterpart candidates for 188 of the 373 sources, after using the HST data to correct the absolute astrometry of our Chandra imaging to 0.1$"$. While 58 of these 188 are associated with point sources potentially in M31, over half (107) of the counterpart candidates are extended background galaxies, 5 are star clusters, 12 are foreground stars, and 6 are supernova remnants. Sources with no clear counterpart candidate are most likely to be undetected background galaxies and low-mass X-ray binaries in M31. The 58 point sources that are not consistent with foreground stars are bright enough that they could be high mass stars in M31; however, all but 8 have optical colors inconsistent with single stars, suggesting that many could be background galaxies or binary counterparts. For point-like counterparts, we examine the star formation history of the surrounding stellar populations to look for a young component that could be associated with a high mass X-ray binary. For the 40 point-like counterpart candidates associated with young populations, we find that their age distribution has two peaks at 15-20 Myr and 40-50 Myr. If we only consider the 8 counterpart candidates with typical high-mass main sequence optical star colors, their age distribution peaks mimic those of the sample of 40. Finally, we find that intrinsic faintness, and not extinction, is the main limitation for finding further counterpart candidates., Comment: 24 pages, 12 figures, 7 tables, accepted for publication in ApJS

Essentially all low-mass X-ray binaries (LMXBs) in the soft state appear to drive powerful equatorial disc winds. A simple mechanism for driving such outflows involves X-ray heating of the top of the disc atmosphere to the Compton temperature. Beyond the Compton radius, the thermal speed exceeds the escape velocity, and mass loss is inevitable. Here, we present the first coupled radiation-hydrodynamic simulation of such thermally-driven disc winds. The main advance over previous modelling efforts is that the frequency-dependent attenuation of the irradiating SED is taken into account. We can therefore relax the approximation that the wind is optically thin throughout which is unlikely to hold in the crucial acceleration zone of the flow. The main remaining limitations of our simulations are connected to our treatment of optically thick regions. Adopting parameters representative of the wind-driving LMXB GRO~J1655-40, our radiation-hydrodynamic model yields a mass-loss rate that is $\simeq5\times$ lower than that suggested by pure hydrodynamic, optically thin models. This outflow rate still represents more than twice the accretion rate and agrees well with the mass-loss rate inferred from Chandra/HETG observations of GRO~J1655-40 at a time when the system had a similar luminosity to that adopted in our simulations. The Fe XXV and Fe XXVI Lyman $\rm{\alpha}~$ absorption line profiles observed in this state are slightly stronger than those predicted by our simulations but the qualitative agreement between observed and simulated outflow properties means that thermal driving is a viable mechanism for powering the disc winds seen in soft-state LMXBs.

If some or all Type Ia SNe arise from accretion onto a massive WD from a companion, then the companion will remain in some form after the SN explosion. Tycho-B is an unusual, relatively hot star along the line of sight to Tycho's SNR -- conclusively shown to be a Type Ia -- and has been suggested as such a companion. If the interior of Tycho's SNR contains unshocked Fe, and if Tycho-B is either within the SNR shell or in the background, then one might hope to see evidence of this in the UV spectrum. Such is the case for SN\,1006, where spectra of the background Schweizer-Middleditch star, as well as two AGNs, show broad absorption lines of Fe II. To test this idea, we have used STIS on HST to obtain a UV spectrum of Tycho-B. The observed spectrum, however, shows no evidence of Fe II absorption. Furthermore, a luminosity distance estimate using UV and optical spectra of Tycho-B suggests that the star is consistent with a foreground interloper. We conclude either that Tycho B is nearer than Tycho's SNR, or that all of the Fe in the interior of Tycho's SNR is more highly ionized., Comment: 9 pages, 7 figures, submitted to MNRAS and comments (particularly constructive criticisms) are always welcome

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter., Comment: 16 pages, 8 figures, 4 tables. Accepted for publication in PASJ

To date, over 220 emission nebulae in M33 have been identified as supernova remnants (SNRs) or SNR candidates, principally through [SII]:Halpha line ratios that are elevated compared to those in H II regions. In many cases, the determination of a high [SII]:Halpha line ratio was made using narrow-band interference filter images and has not been confirmed spectroscopically. Here we present MMT 6.5 m optical spectra that we use to measure [SII]:Halpha and other line ratios in an attempt to determine the nature of these suggested candidates. Of the 197 objects in our sample, 120 have no previously published spectroscopic observations. We confirm that the majority of candidate SNRs have emission line ratios characteristic of SNRs. While no candidates show Doppler-broadened lines expected from young, ejecta-dominated SNRs (> 1000 km/s), a substantial number do exhibit lines that are broader than H II regions. We argue that the majority of the objects with high [SII]:Halpha line ratios (>0.4) are indeed SNRs, but at low surface brightness the distinction between H II regions and SNRs becomes less obvious, and additional criteria, such as X-ray detection, are needed. We discuss the properties of the sample as a whole and compare it with similar samples in other nearby galaxies., Comment: 43 pages, including 2 tables and 14 figures. Accepted for publication in ApJ

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band., Comment: 29 pages, 19 figures, 9 tables. Accepted for publication in PASJ

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission., Comment: 46 pages, 25 figures, 11 tables. Accepted for publication in PASJ

Supernova remnants (SNRs) arise from the interaction between the ejecta of a supernova (SN) explosion and the surrounding circumstellar and interstellar medium. Some SNRs, mostly nearby SNRs, can be studied in great detail. However, to understand SNRs as a whole, large samples of SNRs must be assembled and studied. Here, we describe the radio, optical, and X-ray techniques which have been used to identify and characterize almost 300 Galactic SNRs and more than 1200 extragalactic SNRs. We then discuss which types of SNRs are being found and which are not. We examine the degree to which the luminosity functions, surface-brightness distributions and multi-wavelength comparisons of the samples can be interpreted to determine the class properties of SNRs and describe efforts to establish the type of SN explosion associated with a SNR. We conclude that in order to better understand the class properties of SNRs, it is more important to study (and obtain additional data on) the SNRs in galaxies with extant samples at multiple wavelength bands than it is to obtain samples of SNRs in other galaxies, Comment: Final 2016 draft of a chapter in "Handbook of Supernovae" edited by Athem W. Alsabti and Paul Murdin. Final version available at https://doi.org/10.1007/978-3-319-20794-0_90-1