Source author record

Laura Chomiuk

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

60works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

60 published item(s)

preprint2022arXiv

1RXH J082623.6-505741: a new long-period cataclysmic variable with an evolved donor and a low mass transfer rate

We report the discovery of 1RXH J082623.6-505741, a 10.4 hr orbital period compact binary. Modeling extensive optical photometry and spectroscopy reveals a $\sim 0.4 M_{\odot}$ K-type secondary transferring mass through a low-state accretion disk to a non-magnetic $\sim 0.8 M_{\odot}$ white dwarf. The secondary is overluminous for its mass and dominates the optical spectra at all epochs, and must be evolved to fill its Roche Lobe at this orbital period. The X-ray luminosity $L_X \sim 1$-$2 \times 10^{32}$ erg s$^{-1}$ derived from both new XMM-Newton and archival observations, although high compared to most CVs, still only requires a modest accretion rate onto the white dwarf of $\dot{M} \sim 3 \times 10^{-11}$ to $3 \times 10^{-10} M_{\odot}$ yr$^{-1}$, lower than expected for a cataclysmic variable with an evolved secondary. No dwarf nova outbursts have yet been observed from the system, consistent with the low derived mass transfer rate. Several other cataclysmic variables with similar orbital periods also show unexpectedly low mass transfer rates, even though selection effects disfavor the discovery of binaries with these properties. This suggests the abundance and evolutionary state of long-period, low mass transfer rate cataclysmic variables is worthy of additional attention.

preprint2022arXiv

4FGL J1120.0-2204: A Unique Gamma-ray Bright Neutron Star Binary with an Extremely Low Mass Proto-White Dwarf

We have discovered a new X-ray emitting compact binary that is the likely counterpart to the unassociated Fermi-LAT GeV $γ$-ray source 4FGL J1120.0-2204, the second brightest Fermi source that still remains formally unidentified. Using optical spectroscopy with the SOAR telescope, we have identified a warm ($T_{\textrm{eff}}\sim8500$ K) companion in a 15.1-hr orbit around an unseen primary, which is likely a yet-undiscovered millisecond pulsar. A precise Gaia parallax shows the binary is nearby, at a distance of only $\sim 820$ pc. Unlike the typical "spider" or white dwarf secondaries in short-period millisecond pulsar binaries, our observations suggest the $\sim 0.17\,M_{\odot}$ companion is in an intermediate stage, contracting on the way to becoming an extremely low-mass helium white dwarf (a "pre-ELM" white dwarf). Although the companion is apparently unique among confirmed or candidate millisecond pulsar binaries, we use binary evolution models to show that in $\sim 2$ Gyr, the properties of the binary will match those of several millisecond pulsar-white dwarf binaries with very short ($< 1$ d) orbital periods. This makes 4FGL J1120.0-2204 the first system discovered in the penultimate phase of the millisecond pulsar recycling process.

preprint2022arXiv

A New Flaring Black Widow Candidate and Demographics of Black Widow Millisecond Pulsars in the Galactic Field

We present the discovery of a new optical/X-ray source likely associated with the Fermi $γ$-ray source 4FGL J1408.6-2917. Its high-amplitude periodic optical variability, large spectroscopic radial velocity semi-amplitude, evidence for optical emission lines and flaring, and X-ray properties together imply the source is probably a new black widow millisecond pulsar binary. We compile the properties of the 41 confirmed and suspected field black widows, finding a median secondary mass of $0.027\pm0.003\,M_{\odot}$. Considered jointly with the more massive redback millisecond pulsar binaries, we find that the "spider" companion mass distribution remains strongly bimodal, with essentially zero systems having companion masses between $\sim0.07-0.1\,M_{\odot}$. X-ray emission from black widows is typically softer and less luminous than in redbacks, consistent with less efficient particle acceleration in the intrabinary shock in black widows, excepting a few systems that appear to have more efficient "redback-like" shocks. Together black widows and redbacks dominate the census of the fastest-spinning field millisecond pulsars in binaries with known companion types, making up $\gtrsim$80% of systems with $P_{\rm{spin}}<2\,\rm{ms}$. Similar to redbacks, the neutron star masses in black widows appear on average significantly larger than the canonical $1.4\,M_{\odot}$, and many of the highest-mass neutron stars claimed to date are black widows with $M_{\rm{NS}}\gtrsim2.1\,M_{\odot}$. Both of these observations are consistent with an evolutionary picture where spider millisecond pulsars emerge from short orbital period progenitors that had a lengthy period of mass transfer initiated while the companion was on the main sequence, leading to fast spins and high masses.

preprint2022arXiv

Testing the Momentum-driven Supernova Feedback Paradigm in M31

Momentum feedback from isolated supernova remnants (SNRs) have been increasingly recognized by modern cosmological simulations as a resolution-independent means to implement the effects of feedback in galaxies, such as turbulence and winds. However, the integrated momentum yield from SNRs is uncertain due to the effects of SN clustering and interstellar medium (ISM) inhomogeneities. In this paper, we use spatially-resolved observations of the prominent 10-kpc star-forming ring of M31 to test models of mass-weighted ISM turbulence driven by momentum feedback from isolated, non-overlapping SNRs. We use a detailed stellar-age distribution (SAD) map from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, observationally-constrained SN delay-time distributions, and maps of the atomic and molecular hydrogen to estimate the mass-weighted velocity dispersion using the Martizzi et al. ISM turbulence model. Our estimates are within a factor of 2 of the observed mass-weighted velocity dispersion in most of the ring, but exceed observations at densities $\lesssim 0.2$ cm$^{-3}$ and SN rates $>2.1\times 10^{-4}$ SN yr$^{-1}$ kpc$^{-2}$, even after accounting for plausible variations in stellar-age distribution models and ISM scale height assumptions. We conclude that at high SN rates the momentum deposited is most likely suppressed by the non-linear effects of SN clustering, while at low densities, SNRs reach pressure equilibrium before the cooling phase. These corrections should be introduced in models of momentum-driven feedback and ISM turbulence.

preprint2022arXiv

The first nova eruption in a novalike variable: YZ Ret as seen in X-rays and gamma-rays

Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT and NuSTAR observations complemented by XMM-Newton X-ray grating spectroscopy to probe the physical parameters of the shocked ejecta and the nova-hosting white dwarf. The XMM-Newton observations revealed a super-soft X-ray emission which is dominated by emission lines of CV, CVI, NVI, NVII, and OVIII rather than a blackbody-like continuum, suggesting CO-composition of the white dwarf in a high-inclination binary system. Fermi/LAT detected YZ Ret for 15 days with the gamma-ray spectrum best described by a power law with an exponential cut-off at 1.9 +/-0.6 GeV. In stark contrast with theoretical predictions and in keeping with previous NuSTAR observations of Fermi-detected classical novae (V5855 Sgr and V906 Car), the 3.5-78 keV X-ray emission is found to be two orders of magnitude fainter than the GeV emission. The X-ray emission observed by NuSTAR is consistent with a single-temperature thermal plasma. We detect no non-thermal tail of the GeV emission expected to extend down to the NuSTAR band. NuSTAR observations continue to challenge theories of high-energy emission from shocks in novae.

preprint2022arXiv

The MAVERIC survey: A catalogue of radio sources in southern globular clusters from the Australia Telescope Compact Array

Radio continuum observations offer a new window on compact objects in globular clusters compared to typical X-ray or optical studies. As part of the MAVERIC survey, we have used the Australia Telescope Compact Array to carry out a deep (median central noise level of approximately 4 microJy per beam) radio continuum survey of 26 southern globular clusters at central frequencies of 5.5 and 9.0 GHz. This paper presents a catalogue of 1285 radio continuum sources in the fields of these 26 clusters. Considering the surface density of background sources, we find significant evidence for a population of radio sources in seven of the 26 clusters, and also identify at least 11 previously known compact objects (6 pulsars and 5 X-ray binaries). While the overall density of radio continuum sources with 7.25-GHz flux densities greater than about 20 microJy in typical globular clusters is relatively low, the survey has already led to the discovery of several exciting compact binaries, including a candidate ultracompact black hole X-ray binary in 47 Tuc. Many of the unclassified radio sources near the centres of the clusters are likely to be true cluster sources, and multi-wavelength follow-up will be necessary to classify these objects and better understand the demographics of accreting compact binaries in globular clusters.

preprint2022arXiv

The MAVERIC Survey: Radio catalogs and source counts from deep Very Large Array imaging of 25 Galactic globular clusters

The MAVERIC survey is the first deep radio continuum imaging survey of Milky Way globular clusters, with a central goal of finding and classifying accreting compact binaries, including stellar-mass black holes. Here we present radio source catalogs for 25 clusters with ultra-deep Karl G. Jansky Very Large Array observations. The median observing time was 10 hr per cluster, resulting in typical rms sensitivities of 2.3 and 2.1 uJy per beam at central frequencies of 5.0 and 7.2 GHz, respectively. We detect nearly 1300 sources in our survey at 5 sigma, and while many of these are likely to be background sources, we also find strong evidence for an excess of radio sources in some clusters. The radio spectral index distribution of sources in the cluster cores differs from the background, and shows a bimodal distribution. We tentatively classify the steep-spectrum sources (those much brighter at 5.0 GHz) as millisecond pulsars and the flat-spectrum sources as compact or other kinds of binaries. These provisional classifications will be solidified with the future addition of X-ray and optical data. The outer regions of our images represent a deep, relatively wide field (~ 0.4/sq. deg) and high resolution C band background survey, and we present source counts calculated for this area. We also release radio continuum images for these 25 clusters to the community.

preprint2022arXiv

The MAVERIC Survey: The first radio and X-ray limits on the detached black holes in NGC 3201

The Galactic globular cluster NGC 3201 is the first Galactic globular cluster observed to host dynamically-confirmed stellar-mass black holes, containing two confirmed and one candidate black hole. This result indicates that globular clusters can retain black holes, which has important implications for globular cluster evolution. NGC 3201 has been observed as part of the MAVERIC survey of Galactic globular clusters. We use these data to confirm that there is no radio or X-ray detection of the three black holes, and present the first radio and X-ray limits on these sources. These limits indicate that any accretion present is at an extremely low rate and may be extremely inefficient. In particular, for the system ACS ID #21859, by assuming the system is tidally locked and any accretion is through the capture of the companion's winds, we constrain the radiative efficiency of any accretion to $\leq1.5\times10^{-5}$. We also combine the radio and X-ray source catalogues from the MAVERIC survey with the existing MUSE spectroscopic surveys and the HUGS catalogue of NGC 3201 to provide a catalogue of 42 multiwavelength sources in this cluster. We identify a new red straggler source with X-ray emission, and investigate the multiwavelength properties of the sub-subgiant population in the cluster.

preprint2021arXiv

High-Energy Neutrinos and Gamma-Rays from Non-Relativistic Shock-Powered Transients

Shock interaction has been argued to play a role in powering a range of optical transients, including supernovae (particularly the superluminous class), classical novae, stellar mergers, tidal disruption events, and fast blue optical transients. These same shocks can accelerate relativistic ions, generating high-energy neutrino and gamma-ray emission via hadronic pion production. The recent discovery of time-correlated optical and gamma-ray emission in classical novae has revealed the important role of radiative shocks in powering these events, enabling an unprecedented view of the properties of ion acceleration, including its efficiency and energy spectrum, under similar physical conditions to shocks in extragalactic transients. Here we introduce a model for connecting the radiated optical fluence of non-relativistic transients to their maximal neutrino and gamma-ray fluence. We apply this technique to a wide range of extragalactic transient classes in order to place limits on their contributions to the cosmological high-energy gamma-ray and neutrino backgrounds. Based on a simple model for diffusive shock acceleration at radiative shocks, calibrated to novae, we demonstrate that several of the most luminous transients can accelerate protons up to energies $E_{\rm max} \gtrsim 10^{16}$ eV, sufficient to contribute to the IceCube astrophysical background. Furthermore, several of the considered sources$-$particularly hydrogen-poor supernovae$-$may serve as "hidden" gamma-ray sources due to the high gamma-ray opacity of their ejecta, evading constraints imposed by the non-blazar Fermi-LAT background. However, adopting an ion acceleration efficiency $\sim$ 0.3-1$\%$ motivated by nova observations, we find that currently known classes of non-relativistic, potentially shock-powered transients contribute at most a few percent of the total IceCube background.

preprint2021arXiv

Shocks and dust formation in nova V809 Cep

The discovery that many classical novae produce detectable GeV $γ$-ray emission has raised the question of the role of shocks in nova eruptions. Here we use radio observations of nova V809 Cep (Nova Cep 2013) with the Jansky Very Large Array to show that it produced non-thermal emission indicative of particle acceleration in strong shocks for more than a month starting about six weeks into the eruption, quasi-simultaneous with the production of dust. Broadly speaking, the radio emission at late times -- more than a six months or so into the eruption -- is consistent with thermal emission from $10^{-4} M_\odot$ of freely expanding, $10^4$~K ejecta. At 4.6 and 7.4 GHz, however, the radio light-curves display an initial early-time peak 76 days after the discovery of the eruption in the optical ($t_0$). The brightness temperature at 4.6 GHz on day 76 was greater than $10^5 K$, an order of magnitude above what is expected for thermal emission. We argue that the brightness temperature is the result of synchrotron emission due to internal shocks within the ejecta. The evolution of the radio spectrum was consistent with synchrotron emission that peaked at high frequencies before low frequencies, suggesting that the synchrotron from the shock was initially subject to free-free absorption by optically thick ionized material in front of the shock. Dust formation began around day 37, and we suggest that internal shocks in the ejecta were established prior to dust formation and caused the nucleation of dust.

preprint2020arXiv

A New Likely Redback Millisecond Pulsar Binary with a Massive Neutron Star: 4FGL J2333.1--5527

We present the discovery of a likely new redback millisecond pulsar binary associated with the \emph{Fermi} $γ$-ray source 4FGL J2333.1--5527. Using optical photometric and spectroscopic observations from the SOAR telescope, we identify a low-mass, main sequence-like companion in a 6.9-hr, highly inclined orbit around a suspected massive neutron star primary. Archival XMM-Newton X-ray observations show this system has a hard power-law spectrum $Γ= 1.6\pm0.3$ and $L_X \sim 5 \times 10^{31}$ erg s$^{-1}$, consistent with redback millisecond pulsar binaries. Our data suggest that for secondary masses typical of redbacks, the mass of the neutron star is likely well in excess of $\sim1.4\,M_{\odot}$, but future timing of the radio pulsar is necessary to bolster this tentative conclusion. This work shows that a bevy of nearby compact binaries still await discovery, and that unusually massive neutron stars continue to be common in redbacks.

preprint2020arXiv

Direct evidence for shock-powered optical emission in a nova

Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

preprint2020arXiv

Fermi-LAT Observations of V549 Vel 2017: a Sub-Luminous Gamma-Ray Nova?

We report on the Fermi-LAT detection (with ~5.7 sigma significance) as well as the multi-wavelength analysis of the 2017 nova eruption V549 Vel. Unlike the recent shock-powered novae ASASSN-16ma and V906 Car, the optical and gamma-ray light curves of V549 Vel show no correlation, likely implying relatively weak shocks in the eruption. Gaia detected a candidate progenitor of V549 Vel and found a parallax measurement of 1.91+/-0.39 mas, equivalent to a mode distance of d=560 pc (90% credible interval of 380-1050 pc). The progenitor was also observed by the 2MASS and WISE surveys. When adopting the Gaia distance, the spectral energy distribution of the progenitor is close to that of a G-type star. Swift/XRT detected the supersoft X-ray emission of the nova (kT=30-40 keV) since day 236, and the inferred blackbody size is comparable to that of other novae assuming d=560 pc (i.e., R_bb~5*10^8 cm). However, there is also an unknown astrometric excess noise of 3.2 mas found in the Gaia data, and the inferred distance becomes controversial. If the Gaia distance is accurate, the gamma-ray luminosity of V549 Vel will be as low as L_gamma~4*10^33 erg/s, making it the least luminous gamma-ray nova known so far. This may imply that the shock properties responsible for the gamma-ray emission in V549 Vel are different from those of the more luminous events. If the nova is located farther away, it is likely a symbiotic system with a giant companion as the observed progenitor.

preprint2020arXiv

The Flare-Dominated Accretion Mode of a Radio-Bright Candidate Transitional Millisecond Pulsar

We report new simultaneous X-ray and radio continuum observations of 3FGL J0427.9-6704, a candidate member of the enigmatic class of transitional millisecond pulsars. These XMM-Newton and Australia Telescope Compact Array observations of this nearly edge-on, eclipsing low-mass X-ray binary were taken in the sub-luminous disk state at an X-ray luminosity of ~10^33(d/2.3kpc)^2 erg/s. Unlike the few well-studied transitional millisecond pulsars, which spend most of their disk state in a characteristic high or low accretion mode with occasional flares, 3FGL J0427.9-6704 stayed in the flare mode for the entire X-ray observation of ~20 hours, with the brightest flares reaching ~2*10^34 erg/s. The source continuously exhibited flaring activity on time-scales of ~10-100 sec in both the X-ray and optical/UV. No measurable time delay between the X-ray and optical/UV flares is observed, but the optical/UV flares last longer, and the relative amplitudes of the X-ray and optical/UV flares show a large scatter. The X-ray spectrum can be well-fit with a partially-absorbed power-law (Gamma~1.4-1.5), perhaps due to the edge-on viewing angle. Modestly variable radio continuum emission is present at all epochs, and is not eclipsed by the secondary, consistent with the presence of a steady radio outflow or jet. The simultaneous radio/X-ray luminosity ratio of 3FGL J0427.9-6704 is higher than any known transitional millisecond pulsars and comparable to that of stellar-mass black holes of the same X-ray luminosity, providing additional evidence that some neutron stars can be as radio-loud as black holes.

preprint2020arXiv

The MAVERIC survey: A hidden pulsar and a black hole candidate in ATCA radio imaging of the globular cluster NGC 6397

Using a 16.2 hr radio observation by the Australia Telescope Compact Array (ATCA) and archival Chandra data, we found $>5σ$ radio counterparts to 4 known and 3 new X-ray sources within the half-light radius ($r_\mathrm{h}$) of the Galactic globular cluster NGC 6397. The previously suggested millisecond pulsar (MSP) candidate, U18, is a steep-spectrum ($S_ν\propto ν^α$; $α=-2.0^{+0.4}_{-0.5}$) radio source with a 5.5 GHz flux density of $54.7\pm 4.3~\mathrm{μJy}$. We argue that U18 is most likely a "hidden" MSP that is continuously hidden by plasma shocked at the collision between the winds from the pulsar and companion star. The nondetection of radio pulsations so far is probably the result of enhanced scattering in this shocked wind. On the other hand, we observed 5.5 GHz flux of the known MSP PSR J1740-5340 (U12) to decrease by a factor of $>2.8$ during epochs of 1.4 GHz eclipse, indicating that the radio flux is absorbed in its shocked wind. If U18 is indeed a pulsar whose pulsations are scattered, we note the contrast with U12's flux decrease in eclipse, which argues for two different eclipse mechanisms at the same radio frequency. In addition to U12 and U18, we also found radio associations for 5 other Chandra X-ray sources, four of which are likely background galaxies. The last, U97, which shows strong H$α$ variability, is mysterious; it may be either a quiescent black hole low-mass X-ray binary, or something more unusual.

preprint2020arXiv

Thirty Years of Radio Observations of Type Ia SN 1972E and SN 1895B: Constraints on Circumstellar Shells

We have imaged over 35 years of archival Very Large Array (VLA) observations of the nearby (d$_{\rm{L}}$ $=$ 3.15 Mpc) Type Ia supernovae SN\,1972E and SN\,1895B between 9 and 121 years post-explosion. No radio emission is detected, constraining the 8.5 GHz luminosities of SN\,1972E and SN\,1895B to be L$_{ν,8.5\rm{GHz}}$ $<$ 6.0 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 45 years post-explosion and L$_{ν,8.5\rm{GHz}}$ $<$ 8.9 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 121 years post-explosion, respectively. These limits imply a clean circumstellar medium (CSM), with $n$ $<$ 0.9 cm$^{-3}$ out to radii of a few $\times$ 10$^{18}$ cm, if the SN blastwave is expanding into uniform density material. Due to the extensive time coverage of our observations, we also constrain the presence of CSM shells surrounding the progenitor of SN\,1972E. We rule out essentially all medium and thick shells with masses of 0.05$-$0.3 M$_\odot$ at radii between $\sim$10$^{17}$ and 10$^{18}$ cm, and thin shells at specific radii with masses down to $\lesssim$0.01 M$_\odot$. These constraints rule out swaths of parameter space for a range of single and double degenerate progenitor scenarios, including recurrent nova, core-degenerate objects, ultra-prompt explosions and white dwarf (WD) mergers with delays of a few hundred years between the onset of merger and explosion. Allowed progenitors include WD-WD systems with a significant ($>$ 10$^{4}$ years) delay from the last episode of common envelope evolution and single degenerate systems undergoing recurrent nova, provided that the recurrence timescale i short and the system has been in the nova phase for $\gtrsim$10$^{4}$ yr, such that a large ($>$ 10$^{18}$ cm) cavity has been evacuated. Future multi-epoch observations of additional intermediate-aged Type Ia SNe will provide a comprehensive view of the large-scale CSM environments around these explosions.

preprint2020arXiv

X-ray Spectroscopy of the gamma-ray Brightest Nova V906 Car (ASASSN-18fv)

Shocks in gamma-ray emitting classical novae are expected to produce bright thermal and non-thermal X-rays. We test this prediction with simultaneous NuSTAR and Fermi/LAT observations of nova V906 Car, which exhibited the brightest GeV gamma-ray emission to date. The nova is detected in hard X-rays while it is still gamma-ray bright, but contrary to simple theoretical expectations, the detected 3.5-78 keV emission of V906 Car is much weaker than the simultaneously observed >100 MeV emission. No non-thermal X-ray emission is detected, and our deep limits imply that the gamma-rays are likely hadronic. After correcting for substantial absorption (N_H ~ 2 x 10^23 cm^-2), the thermal X-ray luminosity (from a 9 keV optically-thin plasma) is just ~2% of the gamma-ray luminosity. We consider possible explanations for the low thermal X-ray luminosity, including the X-rays being suppressed by corrugated, radiative shock fronts or the X-rays from the gamma-ray producing shock are hidden behind an even larger absorbing column (N_H >10^25 cm^-2). Adding XMM-Newton and Swift/XRT observations to our analysis, we find that the evolution of the intrinsic X-ray absorption requires the nova shell to be expelled 24 days after the outburst onset. The X-ray spectra show that the ejecta are enhanced in nitrogen and oxygen, and the nova occurred on the surface of a CO-type white dwarf. We see no indication of a distinct super-soft phase in the X-ray lightcurve, which, after considering the absorption effects, may point to a low mass of the white dwarf hosting the nova.

preprint2016arXiv

A new gamma-ray loud, eclipsing low-mass X-ray binary

We report the discovery of an eclipsing low-mass X-ray binary at the center of the 3FGL error ellipse of the unassociated Fermi/Large Area Telescope gamma-ray source 3FGL J0427.9-6704. Photometry from OGLE and the SMARTS 1.3-m telescope and spectroscopy from the SOAR telescope have allowed us to classify the system as an eclipsing low-mass X-ray binary (P = 8.8 hr) with a main sequence donor and a neutron star accretor. Broad double-peaked H and He emission lines suggest the ongoing presence of an accretion disk. Remarkably, the system shows shows separate sets of absorption lines associated with the accretion disk and the secondary, and we use their radial velocities to find evidence for a massive (~ 1.8-1.9 M_sun) neutron star primary. In addition to a total X-ray eclipse of duration ~ 2200 s observed with NuSTAR, the X-ray light curve also shows properties similar to those observed among known transitional millisecond pulsars: short-term variability, a hard power-law spectrum (photon index ~ 1.7), and a comparable 0.5-10 keV luminosity (~ 2.4 x 10^33 erg/s). We find tentative evidence for a partial (~ 60%) gamma-ray eclipse at the same phase as the X-ray eclipse, suggesting the gamma-ray emission may not be confined to the immediate region of the compact object. The favorable inclination of this binary is promising for future efforts to determine the origin of gamma-rays among accreting neutron stars.

preprint2016arXiv

An X-ray and Radio Study of the Varying Expansion Velocities in Tycho's Supernova Remnant

We present newly obtained X-ray and radio observations of Tycho's supernova remnant using {\it Chandra} and the Karl G. Jansky Very Large Array in 2015 and 2013/14, respectively. When combined with earlier epoch observations by these instruments, we now have time baselines for expansion measurements of the remnant of 12-15 year in the X-rays and 30 year in the radio. The remnant's large angular size allows for proper motion measurements at many locations around the periphery of the blast wave. We find, consistent with earlier measurements, a clear gradient in the expansion velocity of the remnant, despite its round shape. The proper motions on the western and southwestern sides of the remnant are about a factor of two higher than those in the east and northeast. We showed in an earlier work that this is related to an offset of the explosion site from the geometric center of the remnant due to a density gradient in the ISM, and using our refined measurements reported here, we find that this offset is $\sim 23"$ towards the northeast. An explosion center offset in such a circular remnant has implications for searches for progenitor companions in other remnants.

preprint2016arXiv

ASASSN-16ae: A Powerful White-Light Flare on an Early-L Dwarf

We report the discovery and classification of SDSS~J053341.43+001434.1 (SDSS0533), an early-L dwarf first discovered during a powerful $ΔV < -11$ magnitude flare observed as part of the ASAS-SN survey. Optical and infrared spectroscopy indicate a spectral type of L0 with strong H$α$ emission and a blue NIR spectral slope. Combining the photometric distance, proper motion, and radial velocity of SDSS0533 yields three-dimensional velocities of $(U,V,W)=(14\pm13,-35\pm14,-94\pm22)$~km~s$^{-1}$, indicating that it is most likely part of the thick disk population and probably old. The three detections of SDSS0533 obtained during the flare are consistent with a total $V$-band flare energy of at least $4.9\times10^{33}$~ergs (corresponding to a total thermal energy of at least $E_{\rm tot}>3.7\times10^{34}$~erg), placing it among the strongest detected M dwarf flares. The presence of this powerful flare on an old L0 dwarf may indicate that stellar-type magnetic activity persists down to the end of the main sequence and on older ML transition dwarfs.

preprint2016arXiv

Discovery of a Redback Millisecond Pulsar Candidate: 3FGL J0212.1+5320

We present a multi-wavelength study of the unidentified Fermi object, 3FGL J0212.1+5320. Within the 95% error ellipse, Chandra detects a bright X-ray source (i.e., F(0.5-7keV) = 1.4e-12 erg/cm^2/s), which has a low-mass optical counterpart (M < 0.4 Msun and T ~ 6000 K). A clear ellipsoidal modulation is shown in optical/infrared at 20.87 hours. The gamma-ray properties of 3FGL J0212.1+5320 are all consistent with that of a millisecond pulsar, suggesting that it is a gamma-ray redback millisecond pulsar binary with a low-mass companion filling > 64% of the Roche-lobe. If confirmed, it will be a redback binary with one of the longest orbital periods known. Spectroscopic data taken in 2015 from the Lijiang observatory show no evidence of strong emission lines, revealing that the accretion is currently inactive (the rotation-powered pulsar state). This is consistent with the low X-ray luminosities (Lx ~ 10^32 erg/s) and the possible X-ray modulation seen by Chandra and Swift. Considering that the X-ray luminosity and the high X-ray-to-gamma-ray flux ratio (8%) are both comparable to that of the two known gamma-ray transitional millisecond pulsars, we suspect that 3FGL J0212.1+5320 could be a potential target to search for future transition to the accretion active state.

preprint2016arXiv

Highest Redshift Image of Neutral Hydrogen in Emission: A CHILES Detection of a Starbursting Galaxy at z=0.376

Our current understanding of galaxy evolution still has many uncertainties associated with the details of accretion, processing, and removal of gas across cosmic time. The next generation of radio telescopes will image the neutral hydrogen (HI) in galaxies over large volumes at high redshifts, which will provide key insights into these processes. We are conducting the COSMOS HI Large Extragalactic Survey (CHILES) with the Karl G. Jansky Very Large Array, which is the first survey to simultaneously observe HI from z=0 to z~0.5. Here, we report the highest redshift HI 21-cm detection in emission to date of the luminous infrared galaxy (LIRG) COSMOS J100054.83+023126.2 at z=0.376 with the first 178 hours of CHILES data. The total HI mass is $(2.9\pm1.0)\times10^{10}~M_\odot$, and the spatial distribution is asymmetric and extends beyond the galaxy. While optically the galaxy looks undisturbed, the HI distribution suggests an interaction with candidate a candidate companion. In addition, we present follow-up Large Millimeter Telescope CO observations that show it is rich in molecular hydrogen, with a range of possible masses of $(1.8-9.9)\times10^{10}~M_\odot$. This is the first study of the HI and CO in emission for a single galaxy beyond z~0.2.

preprint2016arXiv

Identifying IGR J14091-6108 as a magnetic CV with a massive white dwarf using X-ray and optical observations

IGR J14091-6108 is a Galactic X-ray source known to have an iron emission line, a hard X-ray spectrum, and an optical counterpart. Here, we report on X-ray observations of the source with XMM-Newton and NuSTAR as well as optical spectroscopy with ESO/VLT and NOAO/SOAR. In the X-rays, this provides data with much better statistical quality than the previous observations, and this is the first report of the optical spectrum. Timing analysis of the XMM data shows a very significant detection of 576.3+/-0.6 s period. The signal has a pulsed fraction of 30%+/-3% in the 0.3-12 keV range and shows a strong drop with energy. The optical spectra show strong emission lines with significant variability in the lines and continuum, indicating that they come from an irradiated accretion disk. Based on these measurements, we identify the source as a magnetic Cataclysmic Variable of Intermediate Polar (IP) type where the white dwarf spin period is 576.3 s. The X-ray spectrum is consistent with the continuum emission mechanism being due to thermal Bremsstrahlung, but partial covering absorption and reflection are also required. In addition, we use the IP mass (IPM) model, which suggests that the white dwarf in this system has a high mass, possibly approaching the Chandrasekhar limit.

preprint2016arXiv

Shock-powered radio emission from V5589 Sagittarii (Nova Sgr 2012 1)

Since the Fermi discovery of $γ$-rays from novae, one of the biggest questions in the field has been how novae generate such high-energy emission. Shocks must be a fundamental ingredient. Six months of radio observations of the 2012 nova V5589 Sgr with the VLA and 15 weeks of X-ray observations with Swift/XRT show that the radio emission consisted of: 1) a shock-powered, non-thermal flare; and 2) weak thermal emission from $10^{-5}$ M$_\odot$ of freely expanding, photoionized ejecta. Absorption features in the optical spectrum and the peak optical brightness suggest that V5589 Sgr lies 4 kpc away (3.2-4.6 kpc). The shock-powered flare dominated the radio light curve at low frequencies before day 100. The spectral evolution of the radio flare, its high radio brightness temperature, the presence of unusually hard ($kT_x > 33$ keV) X-rays, and the ratio of radio to X-ray flux near radio maximum all support the conclusions that the flare was shock-powered and non-thermal. Unlike most other novae with strong shock-powered radio emission, V5589 Sgr is not embedded in the wind of a red-giant companion. Based on the similar inclinations and optical line profiles of V5589 Sgr and V959 Mon, we propose that shocks in V5589 Sgr formed from collisions between a slow flow with an equatorial density enhancement and a subsequent faster flow. We speculate that the relatively high speed and low mass of the ejecta led to the unusual radio emission from V5589 Sgr, and perhaps also to the non-detection of $γ$-rays.

preprint2016arXiv

Supernova Remnants in the Local Group I: A model for the radio luminosity function and visibility times of supernova remnants

Supernova remnants (SNRs) in Local Group galaxies offer unique insights into the origin of different types of supernovae. In order to take full advantage of these insights, one must understand the intrinsic and environmental diversity of SNRs in the context of their host galaxies. We introduce a semi-analytic model that reproduces the statistical properties of a radio continuum-selected SNR population, taking into account the detection limits of radio surveys, the range of SN kinetic energies, the measured ISM and stellar mass distribution in the host galaxy from multi-wavelength images and the current understanding of electron acceleration and field amplification in SNR shocks from first-principle kinetic simulations. Applying our model to the SNR population in M33, we reproduce the SNR radio luminosity function with a median SN rate of $\sim 3.1 \times 10^{-3}$ per year and an electron acceleration efficiency, $ε_{\rm{e}} \sim 4.2 \times 10^{-3}$. We predict that the radio visibility times of $\sim 70\%$ of M33 SNRs will be determined by their Sedov-Taylor lifetimes, and correlated with the measured ISM column density, $N_H$ ($t_{\rm{vis}} \propto N_H^{-a}$, with $a \sim 0.33$) while the remaining will have visibility times determined by the detection limit of the radio survey. These observational constraints on the visibility time of SNRs will allow us to use SNR catalogs as `SN surveys' to calculate SN rates and delay time distributions in the Local Group.

preprint2015arXiv

1FGL J1417.7-4407: A likely gamma-ray bright binary with a massive neutron star and a giant secondary

We present multiwavelength observations of the persistent Fermi-LAT unidentified gamma-ray source 1FGL J1417.7-4407, showing it is likely to be associated with a newly discovered X-ray binary containing a massive neutron star (nearly 2 M_sun) and a ~ 0.35 M_sun giant secondary with a 5.4 day period. SOAR optical spectroscopy at a range of orbital phases reveals variable double-peaked H-alpha emission, consistent with the presence of an accretion disk. The lack of radio emission and evidence for a disk suggests the gamma-ray emission is unlikely to originate in a pulsar magnetosphere, but could instead be associated with a pulsar wind, relativistic jet, or could be due to synchrotron self-Compton at the disk--magnetosphere boundary. Assuming a wind or jet, the high ratio of gamma-ray to X-ray luminosity (~ 20) suggests efficient production of gamma-rays, perhaps due to the giant companion. The system appears to be a low-mass X-ray binary that has not yet completed the pulsar recycling process. This system is a good candidate to monitor for a future transition between accretion-powered and rotational-powered states, but in the context of a giant secondary.

preprint2015arXiv

A Deep Search for Prompt Radio Emission from Thermonuclear Supernovae with the Very Large Array

Searches for circumstellar material around Type Ia supernovae (SNe Ia) are one of the most powerful tests of the nature of SN Ia progenitors, and radio observations provide a particularly sensitive probe of this material. Here we report radio observations for SNe Ia and their lower-luminosity thermonuclear cousins. We present the largest, most sensitive, and spectroscopically diverse study of prompt (delta t <~ 1 yr) radio observations of 85 thermonuclear SNe, including 25 obtained by our team with the unprecedented depth of the Karl G. Jansky Very Large Array. With these observations, SN 2012cg joins SN 2011fe and SN 2014J as a SN Ia with remarkably deep radio limits and excellent temporal coverage (six epochs, spanning 5--216 days after explosion, yielding Mdot/v_w <~ 5 x 10^-9 M_sun/yr / (100 km/s), assuming epsilon_B = 0.1 and epsilon_e = 0.1). All observations yield non-detections, placing strong constraints on the presence of circumstellar material. We present analytical models for the temporal and spectral evolution of prompt radio emission from thermonuclear SNe as expected from interaction with either wind-stratified or uniform density media. These models allow us to constrain the progenitor mass loss rates, with limits ranging from Mdot <~ 10^-9--10^-4 M_sun/yr, assuming a wind velocity v_w=100 km/s. We compare our radio constraints with measurements of Galactic symbiotic binaries to conclude that <~10% of thermonuclear SNe have red giant companions.

preprint2015arXiv

A Missing-Link in the Supernova-GRB Connection: The Case of SN 2012ap

Gamma Ray Bursts (GRBs) are characterized by ultra-relativistic outflows, while supernovae are generally characterized by non-relativistic ejecta. GRB afterglows decelerate rapidly usually within days, because their low-mass ejecta rapidly sweep up a comparatively larger mass of circumstellar material. However supernovae, with heavy ejecta, can be in nearly free expansion for centuries. Supernovae were thought to have non-relativistic outflows except for few relativistic ones accompanied by GRBs. This clear division was blurred by SN 2009bb, the first supernova with a relativistic outflow without an observed GRB. Yet the ejecta from SN 2009bb was baryon loaded, and in nearly-free expansion for a year, unlike GRBs. We report the first supernova discovered without a GRB, but with rapidly decelerating mildly relativistic ejecta, SN 2012ap. We discovered a bright and rapidly evolving radio counterpart driven by the circumstellar interaction of the relativistic ejecta. However, we did not find any coincident GRB with an isotropic fluence of more than a sixth of the fluence from GRB 980425. This shows for the first time that central engines in type Ic supernovae, even without an observed GRB, can produce both relativistic and rapidly decelerating outflows like GRBs.

preprint2015arXiv

Distance and Reddening of the Enigmatic Gamma-ray-Detected Nova V1324 Sco

It has recently been discovered that some, if not all, classical novae emit GeV gamma-rays during outburst. Despite using an unreliable method to determine its distance, previous work showed that nova V1324 Sco was the most gamma-ray luminous of all gamma-ray-detected novae. We present here a different, more robust, method to determine the reddening and distance to V1324 Sco using high-resolution optical spectroscopy. Using two independent methods we derived a reddening of E(B-V) = 1.16 +/- 0.12 and a distance rD > 6.5 kpc. This distance is >40% greater than previously estimated, meaning that V1324 Sco has an even higher gamma-ray luminosity than previously calculated. We also use periodic modulations in the brightness, interpreted as the orbital period, in conjunction with pre-outburst photometric limits to show that a main-sequence companion is strongly favored.

preprint2015arXiv

Gamma-ray novae as probes of relativistic particle acceleration at non-relativistic shocks

The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into optical emission, similar to Type IIn (interacting) supernovae. Gamma-rays are produced by collisions between relativistic protons with the nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic particle energy into LAT gamma-rays is at most a few tens of per cent. The ratio of gamma-ray and optical luminosities, L_gam/L_opt, thus sets a lower limit on the fraction of the shock power used to accelerate relativistic particles, e_nth. The measured values of L_gam/L_opt for two classical novae, V1324 Sco and V339 Del, constrains e_nth > 1e-2 and > 1e-3, respectively. Inverse Compton models for the gamma-ray emission are disfavored given the low electron acceleration efficiency, e_nth ~ 1e-4-1e-3, inferred from observations of Galactic cosmic rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(0.01/e_nth) and > 10(0.01/e_nth) per cent of the optical luminosity is powered by shocks in V1324 Sco and V339 Del, respectively. Such high fractions challenge standard models that instead attribute all nova optical emission to the direct outwards transport of thermal energy released near the white dwarf surface.

preprint2015arXiv

Limits on thermal variations in a dozen quiescent neutron stars over a decade

In quiescent low-mass X-ray binaries (qLMXBs) containing neutron stars, the origin of the thermal X-ray component may be either release of heat from the core of the neutron star, or continuing low-level accretion. In general, heat from the core should be stable on timescales $<10^4$ years, while continuing accretion may produce variations on a range of timescales. While some quiescent neutron stars (e.g. Cen X-4, Aql X-1) have shown variations in their thermal components on a range of timescales, several others, particularly those in globular clusters with no detectable nonthermal hard X-rays (fit with a powerlaw), have shown no measurable variations. Here, we constrain the spectral variations of 12 low mass X-ray binaries in 3 globular clusters over $\sim10$ years. We find no evidence of variations in 10 cases, with limits on temperature variations below 11% for the 7 qLMXBs without powerlaw components, and limits on variations below 20% for 3 other qLMXBs that do show non-thermal emission. However, in 2 qLMXBs showing powerlaw components in their spectra (NGC 6440 CX 1 & Terzan 5 CX 12) we find marginal evidence for a 10% decline in temperature, suggesting the presence of continuing low-level accretion. This work adds to the evidence that the thermal X-ray component in quiescent neutron stars without powerlaw components can be explained by heat deposited in the core during outbursts. Finally, we also investigate the correlation between hydrogen column density (N$_H$) and optical extinction (A$_V$) using our sample and current models of interstellar X-ray absorption, finding $N_H ({\rm cm}^{-2}) = (2.81\pm0.13)\times10^{21} A_V$.

preprint2015arXiv

Next Generation Very Large Array Memo No. 9 Science Working Group 4: Time Domain, Fundamental Physics, and Cosmology

We report here on key science topics for the Next Generation Very Large Array in the areas of time domain, fundamental physics, and cosmology. Key science cases considered are pulsars in orbit around the Galactic Center massive black hole, Sagittarius A*, electromagnetic counterparts to gravitational waves, and astrometric cosmology. These areas all have the potential for ground-breaking and transformative discovery. Numerous other topics were discussed during the preparation of this report and some of those discussions are summarized here, as well. There is no doubt that further investigation of the science case will reveal rich and compelling opportunities.

preprint2015arXiv

Non-Thermal Radio Emission from Colliding Flows in Classical Nova V1723 Aql

The importance of shocks in nova explosions has been highlighted by Fermi's discovery of γ-ray producing novae. Over three years of multi-band VLA radio observations of the 2010 nova V1723 Aql show that shocks between fast and slow flows within the ejecta led to the acceleration of particles and the production of synchrotron radiation. Soon after the start of the eruption, shocks in the ejecta produced an unexpected radio flare, resulting in a multi-peaked radio light curve. The emission eventually became consistent with an expanding thermal remnant with mass $2 \times 10^{-4} M_\odot$ and temperature $10^4$ K. However, during the first two months, the $\gtrsim 10^6$ K brightness temperature at low frequencies was too high to be due to thermal emission from the small amount of X-ray producing shock-heated gas. Radio imaging showed structures with velocities of 400 km s$^{-1}$ (d/6 kpc) in the plane of the sky, perpendicular to a more elongated 1500 km s$^{-1}$ (d/6 kpc) flow. The morpho-kinematic structure of the ejecta from V1723 Aql appears similar to nova V959 Mon, where collisions between a slow torus and a faster flow collimated the fast flow and gave rise to γ-ray producing shocks. Optical spectroscopy and X-ray observations of V1723 Aql during the radio flare are consistent with this picture. Our observations support the idea that shocks in novae occur when a fast flow collides with a slow collimating torus. Such shocks could be responsible for hard X-ray emission, γ-ray production, and double-peaked radio light curves from some classical novae.

preprint2015arXiv

Optical spectroscopy of the high-mass gamma-ray binary 1FGL J1018.6-5856: A probable neutron star primary

We present medium-resolution optical spectroscopy with the SOAR telescope of the O star secondary of the high-mass gamma-ray binary 1FGL J1018.6-5856 to help determine whether the primary is a neutron star or black hole. We find that the secondary has a low radial velocity semi-amplitude of 11-12 km/s, with consistent values obtained for H and He absorption lines. This low value strongly favors a neutron star primary: while a black hole cannot be excluded if the system is close to face on, such inclinations are disallowed by the observed rotation of the secondary. We also find the high-energy (X-ray and gamma-ray) flux maxima occur when the star is behind the compact object along our line of sight, inconsistent with a simple model of anisotropic inverse Compton scattering for the gamma-ray photons.

preprint2015arXiv

VLBI detection of internal shocks in nova V959 Mon

V959 Mon is a classical nova detected at GeV gamma-ray wavelengths on 2012 June 19. While classical novae are now routinely detected in gamma-rays, the origin of the shocks that produce relativistic particles has remained unknown. We carried out electronic European VLBI Network (e-EVN) observations that revealed a pair of compact synchrotron emission features in V959 Mon on 2012 Sep 18. Since synchrotron emission requires strong shocks as well, we identify these features as the location where the gamma rays were produced. We also detected the extended ejecta in the follow-up EVN observations. They expanded much faster in East-West direction than the compact knots detected in the aforementioned e-EVN measurements. By comparing the VLBI results with lower resolution images obtained using e-MERLIN and the VLA - as reported by Chomiuk et al. (2014) - it appears that 1) influenced by the binary orbit, the nova ejecta was highly asymmetric with a dense and slow outflow in the equatorial plane and low-density and faster ejecta along the poles; and 2) the VLBI knots were related to shocks formed in the interaction region of these outflows.

preprint2014arXiv

1FGL J0523.5-2529: A New Probable Gamma-ray Pulsar Binary

We report optical photometric and SOAR spectroscopic observations of an X-ray source found within the localization error of the Fermi-LAT unidentified gamma-ray source J0523.5-2529. The optical data show periodic flux modulation and radial velocity variations indicative of a binary with a 16.5-hr period. The data suggest a massive non-degenerate secondary (~> 0.8 M_sun), and we argue the source is likely a pulsar binary. The radial velocities have good phase coverage and show evidence for a measurable eccentricity (e=0.04). There is no clear sign of irradiation of the secondary in either photometry or spectroscopy. The spatial location out of the Galactic plane and gamma-ray luminosity of the source are more consistent with classification as a recycled millisecond pulsar than as a young pulsar. Future radio timing observations can confirm the identity of the primary and further characterize this interesting system.

preprint2014arXiv

A Mid-life crisis? Sudden Changes in Radio and X-Ray Emission from Supernova 1970G

Supernovae provide a backdrop from which we can probe the end state of stellar evolution in the final years before the progenitor star explodes. As the shock from the supernova expands, the timespan of mass loss history we are able to probe also extends, providing insight to rapid time-scale processes that govern the end state of massive stars. While supernovae transition into remnants on timescales of decades to centuries, observations of this phase are currently limited. Here we present observations of SN 1970G, serendipitously observed during the monitoring campaign of SN 2011fe that shares the same host galaxy. Utilizing the new Jansky Very Large Array upgrade and a deep X-ray exposure taken by the Chandra Space Telescope, we are able to recover this middle-aged supernova and distinctly resolve it from the HII cloud with which it is associated. We find that the flux density of SN 1970G has changed significantly since it was last observed - the X-ray luminosity has increased by a factor of ~3, while we observe a significantly lower radio flux of only 27.5 micro-Jy at 6.75 GHz, a level only detectable through the upgrades now in operation at the Jansky Very Large Array. These changes suggest that SN 1970G has entered a new stage of evolution towards a supernova remnant, and we may be detecting the turn-on of the pulsar wind nebula. Deep radio observations of additional middle-aged supernovae with the improved radio facilities will provide a statistical census of the delicate transition period between supernova and remnant.

preprint2014arXiv

Binary orbits as the driver of gamma-ray emission and mass ejection in classical novae

Classical novae are the most common astrophysical thermonuclear explosions, occurring on the surfaces of white dwarf stars accreting gas from companions in binary star systems. Novae typically expel ~10^(-4) solar masses of material at velocities exceeding 1,000 kilometres per second. However, the mechanism of mass ejection in novae is poorly understood, and could be dominated by the impulsive flash of thermonuclear energy, prolonged optically thick winds, or binary interaction with the nova envelope. Classical novae are now routinely detected in gigaelectronvolt gamma-ray wavelengths, suggesting that relativistic particles are accelerated by strong shocks in the ejecta. Here we report high-resolution radio imaging of the gamma-ray-emitting nova V959 Mon. We find that its ejecta were shaped by the motion of the binary system: some gas was expelled rapidly along the poles as a wind from the white dwarf, while denser material drifted out along the equatorial plane, propelled by orbital motion. At the interface between the equatorial and polar regions, we observe synchrotron emission indicative of shocks and relativistic particle acceleration, thereby pinpointing the location of gamma-ray production. Binary shaping of the nova ejecta and associated internal shocks are expected to be widespread among novae, explaining why many novae are gamma-ray emitters.

preprint2014arXiv

Hunting The Most Distant Stars in the Milky Way: Methods and Initial Results

We present a new catalog of 404 M giant candidates found in the UKIRT Infrared Deep Sky Survey (UKIDSS). The 2,400 deg$^2$ available in the UKIDSS Large Area Survey Data Release 8 resolve M giants through a volume four times larger than that of the entire Two Micron All Sky Survey. Combining near-infrared photometry with optical photometry and proper motions from the Sloan Digital Sky Survey yields an M giant candidate catalog with less M dwarf and quasar contamination than previous searches for similarly distant M giants. Extensive follow-up spectroscopy of this sample will yield the first map of our Galaxy's outermost reaches over a large area of sky. Our initial spectroscopic follow-up of $\sim$ 30 bright candidates yielded the positive identification of five M giants at distances $\sim 20-90$ kpc. Each of these confirmed M giants have positions and velocities consistent with the Sagittarius stream. The fainter M giant candidates in our sample have estimated photometric distances $\sim 200$ kpc (assuming $[Fe/H]$ = 0.0), but require further spectroscopic verification. The photometric distance estimates extend beyond the Milky Way's virial radius, and increase by $\sim 50\%$ for each 0.5 dex decrease in assumed $[Fe/H]$. Given the number of M giant candidates, initial selection efficiency, and volume surveyed, we loosely estimate that at least one additional Sagittarius-like accretion event could have contributed to the hierarchical build-up of the Milky Way's outer halo.

preprint2014arXiv

Shocks in nova outflows. I. Thermal emission

Evidence for shocks in nova outflows include (1) multiple velocity components in the optical spectra; (2) keV X-ray emission weeks to months after the outburst; (3) early radio flare on timescales of months, in excess of that predicted from the freely expanding photo-ionized gas; and (4) ~ GeV gamma-rays. We present a 1D model for the shock interaction between the fast nova outflow and a dense external shell (DES) and its associated thermal X-ray, optical, and radio emission. The forward shock is radiative initially when the density of shocked gas is highest, at which times radio emission originates from the dense cooling layer immediately downstream of the shock. The radio light curve is characterized by sharper rises to maximum and later peak times at progressively lower frequencies, with a peak brightness temperature that is approximately independent of frequency. We apply our model to the recent gamma-ray classical nova V1324 Sco, obtaining an adequate fit to the early radio maximum for reasonable assumptions about the fast nova outflow and assuming the DES possesses a velocity ~1e3 km/s and mass ~ 2e-4 M_sun; the former is consistent with the velocities of narrow line absorption systems observed previously in nova spectra, while the total ejecta mass of the DES and fast outflow is consistent with that inferred independently by modeling the late radio peak. Rapid evolution of the early radio light curves require the DES possess a steep outer density profile, which may indicate that the onset of mass loss from the white dwarf was rapid, providing indirect evidence that the DES was expelled by the thermonuclear runaway event. Reprocessed X-rays from the shock absorbed by the DES at early times may contribute significantly to the optical/UV emission, which we speculate is responsible for the previously unexplained `plateaus' and secondary maxima in nova optical light curves.

preprint2014arXiv

The 2011 Outburst of Recurrent Nova T Pyx: Radio Observations Reveal the Ejecta Mass and Hint at Complex Mass Loss

Despite being the prototype of its class, T Pyx is arguably the most unusual and poorly understood recurrent nova. Here, we use radio observations from the Karl G. Jansky Very Large Array to trace the evolution of the ejecta over the course of the 2011 outburst of T Pyx. The radio emission is broadly consistent with thermal emission from the nova ejecta. However, the radio flux began rising surprisingly late in the outburst, indicating that the bulk of the radio-emitting material was either very cold, or expanding very slowly, for the first ~50 days of the outburst. Considering a plausible range of volume filling factors and geometries for the ejecta, we find that the high peak flux densities of the radio emission require a massive ejection of 1-30 x 10^{-5} solar masses. This ejecta mass is much higher than the values normally associated with recurrent novae, and is more consistent with a nova on a white dwarf well below the Chandrasekhar limit.

preprint2014arXiv

The 2011 Outburst of Recurrent Nova T Pyx: X-ray Observations Expose the White Dwarf Mass and Ejection Dynamics

The recurrent nova T Pyx underwent its sixth historical outburst in 2011, and became the subject of an intensive multi-wavelength observational campaign. We analyze data from the Swift and Suzaku satellites to produce a detailed X-ray light curve augmented by epochs of spectral information. X-ray observations yield mostly non-detections in the first four months of outburst, but both a super-soft and hard X-ray component rise rapidly after Day 115. The super-soft X-ray component, attributable to the photosphere of the nuclear-burning white dwarf, is relatively cool (~45 eV) and implies that the white dwarf in T Pyx is significantly below the Chandrasekhar mass (~1 M_sun). The late turn-on time of the super-soft component yields a large nova ejecta mass (>~10^-5 M_sun), consistent with estimates at other wavelengths. The hard X-ray component is well fit by a ~1 keV thermal plasma, and is attributed to shocks internal to the 2011 nova ejecta. The presence of a strong oxygen line in this thermal plasma on Day 194 requires a significantly super-solar abundance of oxygen and implies that the ejecta are polluted by white dwarf material. The X-ray light curve can be explained by a dual-phase ejection, with a significant delay between the first and second ejection phases, and the second ejection finally released two months after outburst. A delayed ejection is consistent with optical and radio observations of T Pyx, but the physical mechanism producing such a delay remains a mystery.

preprint2014arXiv

Towards Non-spherical Radio Models

Radio observations of novae in outburst are of particular interest due to the physical parameters that may be retrieved from fitting the radio light curves. Most models that have fitted previous data assumed spherical symmetry however, it is becoming more and more clear that this is not the case. We explore morpho-kinematical techniques to retrieve the free-free radio light curves of non-spherical models and explore the effects of a non-spherical outburst on the physical parameters. In particular, we find that we may have been over estimating the ejected masses in the outburst of non-spherical novae.

preprint2013arXiv

A Pilot for a VLA HI Deep Field

High-resolution 21-cm HI deep fields provide spatially and kinematically resolved neutral gas maps at different redshifts, which are key to understanding galaxy evolution across cosmic time and testing predictions of cosmological simulations. Here we present results from a pilot for the COSMOS HI Large Extragalactic Survey (CHILES) done with the Karl G. Jansky Very Large Array (VLA). We take advantage of the newly expanded capabilities of the telescope to probe the redshift interval 0<z<0.193 in one observation. We observe the COSMOS field for 50 hours, which contains 413 galaxies with optical spectroscopic redshifts in the imaged field of view of 34' x 34' and the observed redshift interval. We have detected neutral hydrogen gas in 33 galaxies in different environments spanning the probed redshift range, including three without a previously known spectroscopic redshift. The detections have a range of HI and stellar masses, indicating the diversity of galaxies we are probing. We discuss the observations, data reduction, results and highlight interesting detections. We find that the VLA's B-array is the ideal configuration for HI deep fields since its long spacings mitigate RFI. This pilot shows that the VLA is ready to carry out such a survey, and serves as a test for future HI deep fields planned with other SKA pathfinders.

preprint2013arXiv

A Radio-Selected Black Hole X-ray Binary Candidate in the Milky Way Globular Cluster M62

We report the discovery of a candidate stellar-mass black hole in the Milky Way globular cluster M62. We detected the black hole candidate, which we term M62-VLA1, in the core of the cluster using deep radio continuum imaging from the Karl G. Jansky Very Large Array. M62-VLA1 is a faint source, with a flux density of 18.7 +/- 1.9 microJy at 6.2 GHz and a flat radio spectrum (alpha=-0.24 +/- 0.42, for S_nu = nu^alpha). M62 is the second Milky Way cluster with a candidate stellar-mass black hole; unlike the two candidate black holes previously found in the cluster M22, M62-VLA1 is associated with a Chandra X-ray source, supporting its identification as a black hole X-ray binary. Measurements of its radio and X-ray luminosity, while not simultaneous, place M62-VLA1 squarely on the well-established radio--X-ray correlation for stellar-mass black holes. In archival Hubble Space Telescope imaging, M62-VLA1 is coincident with a star near the lower red giant branch. This possible optical counterpart shows a blue excess, H alpha emission, and optical variability. The radio, X-ray, and optical properties of M62-VLA1 are very similar to those for V404 Cyg, one of the best-studied quiescent stellar-mass black holes. We cannot yet rule out alternative scenarios for the radio source, such as a flaring neutron star or background galaxy; future observations are necessary to determine whether M62-VLA1 is indeed an accreting stellar-mass black hole.

preprint2013arXiv

A Wolf in Sheepskin: Extraordinary Supernova 2012au Veiled Behind Ordinary Radio Emission

We present extensive radio and X-ray observations of SN\, 2012au, the energetic radio luminous supernova of type Ib that may be a link between subsets of hydrogen-poor superluminous and normal core-collapse supernovae. The observations closely follow models of synchrotron emission from shock heated circum-burst medium that has a wind density profile ($ρ\propto r^{-2}$). We infer a sub-relativistic velocity for the shock wave $v \approx 0.2\,c$ and a radius of $r \approx 1.4 \times 10^{16} \rm cm$ at 25 days after the estimated date of explosion. For a constant wind velocity of 1000 km/s we determine the constant mass loss rate of the progenitor to be $\dot{M} = 3.6 \times 10^{-6} \rm M_{\odot} yr^{-1}$, consistent with the estimates from X-ray observations. We estimate the total internal energy of the radio emitting material to be $E \approx 10^{47} \rm erg$, which is intermediate to SN\,1998bw and SN\,2002ap. Evolution of the radio light curves of SN\,2012au is consistent with interaction with a smoothly distributed circum-burst medium and absence of stellar shells ejected from previous outbursts out to $r \approx 10^{17} \rm cm$ from the supernova site. Based on this we conclude that the evolution of the SN\,2012au progenitor star was relatively quiet during the final years preceding explosion. We find that the bright radio emission from SN2012au was not dissimilar from other core collapse supernovae despite it's extraordinary optical properties. We speculate that it was the nature of the explosion that led to the unusual demise of the SN2012au progenitor star.

preprint2013arXiv

On the Source of the Dust Extinction in Type Ia Supernovae and the Discovery of Anomalously Strong Na I Absorption

High-dispersion observations of the Na I D 5890, 5896 and K I 7665, 7699 interstellar lines, and the diffuse interstellar band at 5780 Angstroms in the spectra of 32 Type Ia supernovae are used as an independent means of probing dust extinction. We show that the dust extinction of the objects where the diffuse interstellar band at 5780 Angstroms is detected is consistent with the visual extinction derived from the supernova colors. This strongly suggests that the dust producing the extinction is predominantly located in the interstellar medium of the host galaxies and not in circumstellar material associated with the progenitor system. One quarter of the supernovae display anomalously large Na I column densities in comparison to the amount of dust extinction derived from their colors. Remarkably, all of the cases of unusually strong Na I D absorption correspond to "Blueshifted" profiles in the classification scheme of Sternberg et al. (2011). This coincidence suggests that outflowing circumstellar gas is responsible for at least some of the cases of anomalously large Na I column densities. Two supernovae with unusually strong Na I D absorption showed essentially normal K I column densities for the dust extinction implied by their colors, but this does not appear to be a universal characteristic. Overall, we find the most accurate predictor of individual supernova extinction to be the equivalent width of the diffuse interstellar band at 5780 Angstroms, and provide an empirical relation for its use. Finally, we identify ways of producing significant enhancements of the Na abundance of circumstellar material in both the single-degenerate and double-degenerate scenarios for the progenitor system.

preprint2013arXiv

Shocks and Ejecta Mass: Radio Observations of Nova V1723 Aql

The radio light curves of novae rise and fall over the course of months to years, allowing for detailed observations of the evolution of the nova shell. However, the main parameter determined by radio models of nova explosions - the mass of the ejecta - often seems to exceed theoretical expectations by an order of magnitude. With the recent technological improvements on the Karl G. Jansky Very Large Array (VLA), new observations can test the assumptions upon which ejecta mass estimates are based. Early observations of the classical nova V1723 Aql showed an unexpectedly rapid rise in radio flux density and a distinct bump in the radio light curve on the rise to radio maximum, which is inconsistent with the simple model of spherical ejecta expelled in a single discrete event. This initial bump appears to indicate the presence of shocked material in the outer region of the ejected shell, with the emission from the shocks fading over time. We explore possible origins for this emission and its relation to the mass loss history of the nova. The evolution of the radio spectrum also reveals the density profile, the mass of the ejected shell, and other properties of the ejecta. These observations comprise one of the most complete, longterm set of multi-wavelength radio observations for any classical nova to date.

preprint2013arXiv

SN 2011fe: A Laboratory for Testing Models of Type Ia Supernovae

SN 2011fe is the nearest supernova of Type Ia (SN Ia) discovered in the modern multi-wavelength telescope era, and it also represents the earliest discovery of a SN Ia to date. As a normal SN Ia, SN 2011fe provides an excellent opportunity to decipher long-standing puzzles about the nature of SNe Ia. In this review, we summarize the extensive suite of panchromatic data on SN 2011fe, and gather interpretations of these data to answer four key questions: 1) What explodes in a SN Ia? 2) How does it explode? 3) What is the progenitor of SN 2011fe? and 4) How accurate are SNe Ia as standardizeable candles? Most aspects of SN 2011fe are consistent with the canonical picture of a massive CO white dwarf undergoing a deflagration-to-detonation transition. However, there is minimal evidence for a non-degenerate companion star, so SN 2011fe may have marked the merger of two white dwarfs.

preprint2012arXiv

EVLA Observations Constrain the Environment and Progenitor System of Type Ia Supernova 2011fe

We report unique EVLA observations of SN 2011fe representing the most sensitive radio study of a Type Ia supernova to date. Our data place direct constraints on the density of the surrounding medium at radii ~10^15-10^16 cm, implying an upper limit on the mass loss rate from the progenitor system of Mdot <~ 6 x 10^-10 Msol/yr (assuming a wind speed of 100 km/s), or expansion into a uniform medium with density n_CSM <~ 6 cm^-3. Drawing from the observed properties of non-conservative mass transfer among accreting white dwarfs, we use these limits on the density of the immediate environs to exclude a phase space of possible progenitors systems for SN 2011fe. We rule out a symbiotic progenitor system and also a system characterized by high accretion rate onto the white dwarf that is expected to give rise to optically-thick accretion winds. Assuming that a small fraction, 1%, of the mass accreted is lost from the progenitor system, we also eliminate much of the potential progenitor parameter space for white dwarfs hosting recurrent novae or undergoing stable nuclear burning. Therefore, we rule out the most popular single degenerate progenitor models for SN 2011fe, leaving a limited phase space inhabited by some double degenerate systems and exotic progenitor scenarios.

preprint2012arXiv

No Evidence for Intermediate-Mass Black Holes in Globular Clusters: Strong Constraints from the JVLA

With a goal of searching for accreting intermediate-mass black holes (IMBHs), we report the results of ultra-deep Jansky VLA radio continuum observations of the cores of three Galactic globular clusters: M15, M19, and M22. We reach rms noise levels of 1.5-2.1 uJy/beam at an average frequency of 6 GHz. No sources are observed at the center of any of the clusters. For a conservative set of assumptions about the properties of the accretion, we set 3-sigma upper limits on IMBHs from 360-980 M_sun. These limits are among the most stringent obtained for any globular cluster. They add to a growing body of work that suggests either (a) IMBHs ~> 1000 M_sun are rare in globular clusters, or (b) when present, IMBHs accrete in an extraordinarily inefficient manner.

preprint2012arXiv

The Radio Light Curve of the Gamma-Ray Nova in V407 Cyg: Thermal Emission from the Ionized Symbiotic Envelope, Devoured from Within by the Nova Blast

We present multi-frequency radio observations of the 2010 nova event in the symbiotic binary V407 Cygni, obtained with the Karl G. Jansky Very Large Array and spanning 1-45 GHz and 17-770 days following discovery. This nova---the first ever detected in gamma rays---shows a radio light curve dominated by the wind of the Mira giant companion, rather than the nova ejecta themselves. The radio luminosity grew as the wind became increasingly ionized by the nova outburst, and faded as the wind was violently heated from within by the nova shock. This study marks the first time that this physical mechanism has been shown to dominate the radio light curve of an astrophysical transient. We do not observe a thermal signature from the nova ejecta or synchrotron emission from the shock, due to the fact that these components were hidden behind the absorbing screen of the Mira wind. We estimate a mass loss rate for the Mira wind of Mdot_w ~ 10^-6 M_sun/yr. We also present the only radio detection of V407 Cyg before the 2010 nova, gleaned from unpublished 1993 archival VLA data, which shows that the radio luminosity of the Mira wind varies by a factor of >~20 even in quiescence. Although V407 Cyg likely hosts a massive accreting white dwarf, making it a candidate progenitor system for a Type Ia supernova, the dense and radially continuous circumbinary material surrounding V407 Cyg is inconsistent with observational constraints on the environments of most Type Ia supernovae.

preprint2012arXiv

Two stellar-mass black holes in the globular cluster M22

Hundreds of stellar-mass black holes likely form in a typical globular star cluster, with all but one predicted to be ejected through dynamical interactions. Some observational support for this idea is provided by the lack of X-ray-emitting binary stars comprising one black hole and one other star ("black-hole/X-ray binaries") in Milky Way globular clusters, even though many neutron-star/X-ray binaries are known. Although a few black holes have been seen in globular clusters around other galaxies, the masses of these cannot be determined, and some may be intermediate-mass black holes that form through exotic mechanisms. Here we report the presence of two flat-spectrum radio sources in the Milky Way globular cluster M22, and we argue that these objects are black holes of stellar mass (each ~ 10-20 times more massive than the Sun) that are accreting matter. We find a high ratio of radio-to-X-ray flux for these black holes, consistent with the larger predicted masses of black holes in globular clusters compared to those outside. The identification of two black holes in one cluster shows that the ejection of black holes is not as efficient as predicted by most models, and we argue that M22 may contain a total population of ~ 5-100 black holes. The large core radius of M22 could arise from heating produced by the black holes.

preprint2012arXiv

X-ray Emission from an Asymmetric Blast Wave and a Massive White Dwarf in the Gamma-ray Emitting Nova V407 Cyg

Classical nova events in symbiotic stars, although rare, offer a unique opportunity to probe the interaction between ejecta and a dense environment in stellar explosions. In this work, we use X-ray data obtained with Swift and Suzaku during the recent classical nova outburst in V407 Cyg to explore such an interaction. We find evidence of both equilibrium and non-equilibrium ionization plasmas at the time of peak X-ray brightness, indicating a strong asymmetry in the density of the emitting region. Comparing a simple model to the data, we find that the X-ray evolution is broadly consistent with nova ejecta driving a forward shock into the dense wind of the Mira companion. We detect a highly absorbed soft X-ray component in the spectrum during the first 50 days of the outburst that is consistent with supersoft emission from the nuclear burning white dwarf. The high temperature and short turn off time of this emission component, in addition to the observed breaks in the optical and UV lightcurves, indicate that the white dwarf in the binary is extremely massive. Finally, we explore the connections between the X-ray and GeV gamma-ray evolution, and propose that the gamma ray turn-off is due to the stalling of the forward shock as the ejecta reach the red giant surface.

preprint2011arXiv

EVLA Nova Project Observations of the Classical Nova V1723 Aquilae

We present radio light curves and spectra of the classical nova V1723 Aql obtained with the Expanded Very Large Array (EVLA). This is the first paper to showcase results from the EVLA Nova Project, which comprises a team of observers and theorists utilizing the greatly enhanced sensitivity and frequency coverage of EVLA radio observations, along with observations at other wavelengths, to reach a deeper understanding of the energetics, morphology, and temporal characteristics of nova explosions. Our observations of V1723 Aql span 1-37 GHz in frequency, and we report on data from 14-175 days following the time of the nova explosion. The broad frequency coverage and frequent monitoring show that the radio behavior of V1723 Aql does not follow the classic Hubble-flow model of homologous spherically expanding thermal ejecta. The spectra are always at least partially optically thin, and the flux rises on faster timescales than can be reproduced with linear expansion. Therefore, any description of the underlying physical processes must go beyond this simple picture. The unusual spectral properties and light curve evolution might be explained by multiple emitting regions or shocked material. Indeed, X-ray observations from Swift reveal that shocks are likely present.

preprint2011arXiv

Microwave Observations of Edge-on Protoplanetary Disks: Program Overview and First Results

We are undertaking a multi-frequency Expanded Very Large Array (EVLA) survey of edge-on protoplanetary disks to probe the growth of solids in each disk, sedimentation of such material into the disk midplane, and the connection of these phenomena to the planet formation process. The projection of edge-on disk systems along our line of sight enables a study of the vertical stratification of large grains with fewer model dependencies than would be required for disks that are more face-on. Robust studies of the spatial distribution of grains up to ~1 cm in size are possible with the wavelength range and sensitivity of the EVLA. In this contribution we describe target selection and observational strategies. First results concerning the Class 0 source IRAS04368+2557 (L1527 IRS) are presented, including a study of this source's 8.46 GHz continuum variability over short and long time baselines and an indication that its protoplanetary disk may have a dearth of pebble-sized grains.

preprint2011arXiv

The Fossil Record of Two-phase Galaxy Assembly: Kinematics and Metallicities in the Nearest S0 Galaxy

We present a global analysis of kinematics and metallicity in the nearest S0 galaxy, NGC 3115, along with implications for its assembly history. The data include high-quality wide-field imaging from Suprime-Cam on the Subaru telescope, and multi-slit spectra of the field stars and globular clusters (GCs) obtained using Keck-DEIMOS/LRIS and Magellan-IMACS. Within two effective radii, the bulge (as traced by the stars and metal-rich GCs) is flattened and rotates rapidly (v/sigma > 1.5). At larger radii, the rotation declines dramatically to v/sigma ~ 0.7, but remains well-aligned with the inner regions. The radial decrease in characteristic metallicity of both the metal-rich and metal-poor GC subpopulations produces strong gradients with power law slopes of -0.17 +/- 0.04 and -0.38 +/- 0.06 dex per dex, respectively. We argue that this pattern is not naturally explained by a binary major merger, but instead by a two-phase assembly process where the inner regions have formed in an early violent, dissipative phase, followed by the protracted growth of the outer parts via minor mergers with typical mass ratios of ~ 15-20:1.

preprint2011arXiv

Toward a Unification of Star Formation Rate Determinations in the Milky Way and Other Galaxies

The star formation rate (SFR) of the Milky Way remains poorly known, with often-quoted values ranging from 1 to 10 solar masses per year. This situation persists despite the potential for the Milky Way to serve as the ultimate SFR calibrator for external galaxies. We show that various estimates for the Galactic SFR are consistent with one another once they have been normalized to the same initial mass function (IMF) and massive star models, converging to 1.9 +/- 0.4 M_sun/yr. However, standard SFR diagnostics are vulnerable to systematics founded in the use of indirect observational tracers sensitive only to high-mass stars. We find that absolute SFRs measured using resolved low/intermediate-mass stellar populations in Galactic H II regions are systematically higher by factors of ~2-3 as compared with calibrations for SFRs measured from mid-IR and radio emission. We discuss some potential explanations for this discrepancy and conclude that it could be allayed if (1) the power-law slope of the IMF for intermediate-mass (1.5 M_sun < m < 5 M_sun) stars were steeper than the Salpeter slope, or (2) a correction factor was applied to the extragalactic 24 micron SFR calibrations to account for the duration of star formation in individual mid-IR-bright H II regions relative to the lifetimes of O stars. Finally, we present some approaches for testing if a Galactic SFR of ~2 M_sun/yr is consistent with what we would measure if we could view the Milky Way as external observers. Using luminous radio supernova remnants and X-ray point sources, we find that the Milky Way deviates from expectations at the 1-3 sigma level, hinting that perhaps the Galactic SFR is overestimated or extragalactic SFRs need to be revised upwards.

preprint2011arXiv

Unveiling Extragalactic Star Formation Using Radio Recombination Lines: An EVLA Pilot Study with NGC 253

Radio recombination lines (RRLs) are powerful, extinction-free diagnostics of the ionized gas in young, star-forming regions. Unfortunately, these lines are difficult to detect in external galaxies. We present the results of EVLA observations of the RRL and radio continuum emission at 33 GHz from NGC 253, a nearby nuclear starburst galaxy. We detect the previously unobserved H58a and H59a RRLs and make simultaneous sensitive measurements of the continuum. We measure integrated line fluxes of $44.3 \pm 0.7$ W m$^{-2}$ and $39.9 \pm 0.8$ W m$^{-2}$ for the H58a and H59a lines, respectively. The thermal gas in NGC 253 is kinematically complex with multiple velocity components. We constrain the density of the thermal gas to $1.4 - 4 \times 10^4$ cm$^{-3}$ and estimate an ionizing photon flux of $1 \times 10^{53}$ s$^{-1}$. We use the RRL kinematics and the derived ionizing photon flux to show that the nuclear region of NGC 253 is not gravitationally bound, which is consistent with the outflow of gas inferred from the X-ray and Halpha measurements. The line profiles, fluxes, and kinematics of the H58a and H59a lines agree with those of RRLs at different frequencies confirming the accuracy of the previous, more difficult, high frequency observations. We find that the EVLA is an order of magnitude more efficient for extragalactic RRL observations than the VLA. These observations demonstrate both the power of the EVLA and the future potential of extragalactic RRL studies with the EVLA.

preprint2009arXiv

A Universal Luminosity Function for Radio Supernova Remnants

We compile radio supernova remnant (SNR) samples from the literature for 19 nearby galaxies ranging from the SMC to Arp 220, and use this data to constrain the SNR luminosity function (LF) at 20 cm. We find that radio SNR populations are strikingly similar across galaxies. The LF can be described as a power law with constant index and scaling proportional to a galaxy's star formation rate (SFR). Unlike previous authors, we do not find any dependence of SNR luminosity on a galaxy's global ISM density. The observed correlation between the luminosity of a galaxy's brightest SNR and a galaxy's SFR can be completely explained by statistical effects, wherein galaxies with higher SFR more thoroughly sample the high-luminosity end of the SNR LF. The LF is well fit by a model of SNR synchrotron emission which includes diffusive shock acceleration and magnetic field amplification, if we assume that all remnants are undergoing adiabatic expansion, the densities of star-forming regions are similar across galaxies, and the efficiency of CR production is constant.