sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Williams, Peter K. G.

Normalized to: Williams, P.

175 article(s) in total. 2698 co-authors, from 1 to 55 common article(s). Median position in authors list is 6,0.

[1]  oai:arXiv.org:2005.10275  [pdf] - 2126816
Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array
Comments: 22 pages, 12 figures (including sub-figures). Published in PhRvD. Abstract may be slightly abridged compared to the actual manuscript due to length limitations on arXiv
Submitted: 2020-05-20, last modified: 2020-07-02
Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.
[2]  oai:arXiv.org:2006.13428  [pdf] - 2121549
EarthFinder Probe Mission Concept Study: Characterizing nearby stellar exoplanet systems with Earth-mass analogs for future direct imaging
Comments: NASA Probe Mission concept white paper for 2020 Astrophysics National Academies Decadal Survey
Submitted: 2020-06-23
EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to the challenges of PRV measurements at the cm/s level. This new paradigm takes advantage of: 1) broad wavelength coverage from the UV to NIR which is only possible from space to minimize the effects of stellar activity; 2) extremely compact, highly stable, highly efficient spectrometers (R>150,000) which require the diffraction-limited imaging possible only from space over a broad wavelength range; 3) the revolution in laser-based wavelength standards to ensure cm/s precision over many years; 4) a high cadence observing program which minimizes sampling-induced period aliases; 5) exploiting the absolute flux stability from space for continuum normalization for unprecedented line-by-line analysis not possible from the ground; and 6) focusing on the bright stars which will be the targets of future imaging missions so that EarthFinder can use a ~1.5 m telescope.
[3]  oai:arXiv.org:2006.10066  [pdf] - 2117698
TDCOSMO II: 6 new time delays in lensed quasars from high-cadence monitoring at the MPIA 2.2m telescope
Comments: 14 pages, 9 figures, 4 Tables, submitted to A&A
Submitted: 2020-06-17
We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the MPIA 2.2m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measure time delays between all pairs of multiple images with only one or two seasons of monitoring at the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate is obtained for the delay between image A and image B of DES 0407-5006, with $\tau_{AB} = -128.4^{+3.5}_{-3.8}$ days (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combine our high-cadence data with measurements from decade-long light-curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 day on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.
[4]  oai:arXiv.org:2006.08695  [pdf] - 2115454
Revisiting the Impact of Dust Production from Carbon-Rich Wolf-Rayet Binaries
Comments: 34 pages, 12 figures, 7 tables, Accepted for publication in ApJ
Submitted: 2020-06-15
We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $\chi_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $\mu$m) dust grain size composition is favored for efficient dustars where $\chi_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $\chi_C$, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a co-eval $10^6$ M$_\odot$ stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and co-eval models show that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of $2-3$, whereas at solar metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr).
[5]  oai:arXiv.org:2006.07434  [pdf] - 2114314
A Late-time Radio Survey of Short GRBs at $z<0.5$: New Constraints on the Remnants of Neutron Star Mergers
Comments: 20 pages, 7 figures, submitted to ApJ
Submitted: 2020-06-12
Massive, rapidly-spinning magnetar remnants produced as a result of binary neutron star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. Here, we present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (GRBs; $z<0.5$) on rest-frame timescales of $\approx 2.4-13.9$ yr following the bursts. We place $3\sigma$ limits on radio continuum emission of $F_{\nu}\lesssim 6-20\,\mu$Jy at the burst positions, or $L_{\nu} \lesssim (0.6-8.3) \times10^{28}$ erg s$^{-1}$ Hz$^{-1}$. Comparing these limits with new light curve modeling which properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of $E_{\rm ej}\lesssim (0.6-6.7)\times 10^{52}$ erg ($E_{\rm ej}\lesssim (1.7-14.1) \times 10^{52}$ erg) for an ejecta mass of $0.03\,M_{\odot}$ ($0.1\,M_{\odot}$). We present a uniform re-analysis of 27 short GRBs with 5.5-6.0 GHz observations, and find that $\gtrsim 50\%$ of short GRBs did not form stable magnetar remnants in their mergers. Assuming short GRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction $f_{\rm GW190425}$ of cases, we place constraints on the maximum mass of a non-rotating neutron star (NS) ($M_{\rm TOV}$), finding $M_{\rm TOV} \lesssim 2.4\,M_{\odot}$ for $f_{\rm GW190425}=0.4$. The detection (or lack thereof) of radio remnants in untargeted surveys such as the VLA Sky Survey (VLASS) could provide more stringent constraints on the fraction of mergers that produce stable remnants. If $\gtrsim 100-1000$ radio remnants are discovered in VLASS, this suggests that short GRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce.
[6]  oai:arXiv.org:2005.12174  [pdf] - 2129379
Measuring HERA's primary beam in-situ: methodology and first results
Comments: 22 pages, 22 figures, accepted for publication in ApJ
Submitted: 2020-05-25
The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.
[7]  oai:arXiv.org:2001.02688  [pdf] - 2109955
Wandering Massive Black Holes or Analogs of the First Repeating Fast Radio Burst?
Comments: 11 pages, 6 figures; replaced with referee revisions including minor changes to figures and text
Submitted: 2020-01-08, last modified: 2020-05-14
The discovery of a persistent radio source coincident with the first repeating fast radio burst, FRB 121102, and offset from the center of its dwarf host galaxy has been used as evidence for a link with young millisecond magnetars born in superluminous supernovae (SLSNe) or long-duration gamma-ray bursts (LGRBs). A prediction of this scenario is that compact radio sources offset from the centers of dwarf galaxies may serve as signposts for at least some FRBs. Recently, Reines et al. 2019 presented the discovery of 20 such radio sources in nearby ($z\lesssim 0.055$) dwarf galaxies, and argued that these cannot be explained by emission from HII regions, normal supernova remnants, or normal radio supernovae. Instead, they attribute the emission to accreting wandering massive black holes. Here, we explore the alternative possibility that these sources are analogs of FRB 121102. We compare their properties -- radio luminosities, spectral energy distributions, light curves, ratios of radio-to-optical flux, and spatial offsets -- to FRB 121102, a few other well-localized FRBs, and potentially related systems, and find that these are all consistent as arising from the same population. We further compare their properties to the magnetar nebula model used to explain FRB 121102, as well as to theoretical off-axis LGRB light curves, and find overall consistency. Finally, we find a consistent occurrence rate relative to repeating FRBs and LGRBs. We outline key follow-up observations to further test these possible connections.
[8]  oai:arXiv.org:2004.06041  [pdf] - 2097404
Foreground modelling via Gaussian process regression: an application to HERA data
Comments: 15 pages, 15 figures, 1 table, Accepted to MNRAS
Submitted: 2020-04-13, last modified: 2020-05-12
The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.
[9]  oai:arXiv.org:2005.00531  [pdf] - 2093346
Two Wolf-Rayet stars at the heart of colliding-wind binary Apep
Comments: Accepted for publication in MNRAS, 10 pages, 9 figures, 1 table
Submitted: 2020-05-01
Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the difficulties in interpreting the dynamics of Apep is that the spectral composition of the stars in the system was unclear. Here we present visual to near-infrared spectra that demonstrate that the central component of Apep is composed of two classical Wolf-Rayet stars of carbon- (WC8) and nitrogen-sequence (WN4-6b) subtypes. We argue that such an assignment represents the strongest case of a classical WR+WR binary system in the Milky Way. The terminal line-of-sight wind velocities of the WC8 and WN4-6b stars are measured to be $2100 \pm 200$ and $3500 \pm 100$ km s$^{-1}$, respectively. If the mass-loss rate of the two stars are typical for their spectral class, the momentum ratio of the colliding winds is expected to be $\approx$ 0.4. Since the expansion velocity of the dust plume is significantly smaller than either of the measured terminal velocities, we explore the suggestion that one of the Wolf-Rayet winds is anisotropic. We can recover a shock-compressed wind velocity consistent with the observed dust expansion velocity if the WC8 star produces a significantly slow equatorial wind with a velocity of $\approx$530 km s$^{-1}$. Such slow wind speeds can be driven by near-critical rotation of a Wolf-Rayet star.
[10]  oai:arXiv.org:2004.10579  [pdf] - 2093318
$K2$ Ultracool Dwarfs Survey. VI. White light superflares observed on an L5 dwarf and flare rates of L dwarfs
Comments: 12 pages, 4 figures, accepted for publication in MNRAS
Submitted: 2020-04-22
Kepler K2 long cadence data are used to study white light flares in a sample of 45 L dwarfs. We identified 11 flares on 9 L dwarfs with equivalent durations of (1.3 - 198) hr and total (UV/optical/IR) energies of $\geq$0.9 $\times$ 10$^{32}$ erg. Two superflares with energies of $>$10$^{33}$ erg were detected on an L5 dwarf: this is the coolest object so far on which flares have been identified. The larger superflare on this L5 dwarf has an energy of 4.6$\times$ 10$^{34}$ ergs and an amplitude of $>$300 times the photospheric level: so far, this is the largest amplitude flare detected by the $Kepler/K2$ mission. The next coolest star on which we identified a flare was an L2 dwarf: 2MASS J08585891+1804463. Combining the energies of all the flares which we have identified on 9 L dwarfs with the total observation time which was dedicated by $Kepler$ to all 45 L dwarfs, we construct a composite flare frequency distribution (FFD). The FFD slope is quite shallow (-0.51$\pm$0.17), consistent with earlier results reported by Paudel et al. (2018) for one particular L0 dwarf, for which the FFD slope was found to be -0.34. Using the composite FFD, we predict that, in early and mid-L dwarfs, a superflare of energy 10$^{33}$ erg occurs every 2.4 years and a superflare of energy 10$^{34}$ erg occurs every 7.9 years. Analysis of our L dwarf flares suggests that magnetic fields of $\geq$0.13-1.3 kG are present on the stellar surface: such fields could suppress Type II radio bursts.
[11]  oai:arXiv.org:1908.08913  [pdf] - 2076637
A Galaxy-Targeted Search for the Optical Counterpart of the Candidate NS-BH Merger S190814bv with Magellan
Comments: 10 pages, 1 table, 4 figures, accepted to ApJL
Submitted: 2019-08-23, last modified: 2020-04-11
On 2019 August 14 the Laser Interferometer Gravitational Wave Observatory (LIGO) and the Virgo gravitational wave interferometer announced the detection of a binary merger, S190814bv, with a low false alarm rate (FAR) of about 1 in $1.6\times 10^{25}$ years, a distance of $267\pm 52$ Mpc, a 90\% (50\%) localization region of about 23 (5) deg$^2$, and a probability of being a neutron star--black hole (NS-BH) merger of $>99\%$. The LIGO/Virgo Collaboration (LVC) defines NS-BH such that the lighter binary member has a mass of $<3$ M$_\odot$ and the more massive one has $>5$ M$_\odot$, and this classification is in principle consistent with a BH-BH merger depending on the actual upper mass cut-off for neutron stars. Additionally, the LVC designated a probability that the merger led to matter outside the final BH remnant of $<1\%$, suggesting that an electromagnetic (EM) counterpart is unlikely. Here we report our optical follow-up observations of S190814bv using the Magellan Baade 6.5 m telescope to target all 96 galaxies in the GLADE catalog within the 50\% localization volume (representing about 70\% of the integrated luminosity within this region). No counterpart was identified to a median $3\sigma$ limiting magnitude of $i=22.2$ ($M_i\approx -14.9$ mag), comparable to the brightness of the optical counterpart of the binary neutron star merger GW170817 at the distance of S190814bv; similarly, we can rule out an on-axis jet typical of short GRBs. However, we cannot rule out other realistic models, such as a kilonova with only $\sim 0.01$ M$_\odot$ of lanthanide-rich material, or an off-axis jet with a viewing angle of $\theta_{\rm obs}\gtrsim 15^\circ$.
[12]  oai:arXiv.org:2003.08399  [pdf] - 2067038
Redundant-Baseline Calibration of the Hydrogen Epoch of Reionization Array
Comments: 23 Pages, 19 Figures
Submitted: 2020-03-18
In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.
[13]  oai:arXiv.org:1909.01395  [pdf] - 2042190
Modelling the AGN broad line region using single-epoch spectra II. Nearby AGN
Comments: Accepted for publication in MNRAS
Submitted: 2019-09-03, last modified: 2020-01-28
The structure of the broad line region (BLR) is an essential ingredient in the determination of active galactic nuclei (AGN) virial black hole masses, which in turn are important to study the role of black holes in galaxy evolution. Constraints on the BLR geometry and dynamics can be obtained from velocity-resolved studies using reverberation mapping data (i.e. monitoring data). However, monitoring data are observationally expensive and only available for a limited sample of AGN, mostly confined to the local Universe. Here we explore a new version of a Bayesian inference, physical model of the BLR which uses an individual spectrum and prior information on the BLR size from the radius-luminosity relation, to model the AGN BLR geometry and dynamics. We apply our model to a sample of 11 AGN, which have been previously modelled using monitoring data. Our single-epoch BLR model is able to constrain some of the BLR parameters with inferred parameter values that agree within the uncertainties with those determined from the modelling of monitoring data. We find that our model is able to derive stronger constraints on the BLR for AGN with broad emission lines that qualitatively have more substructure and more asymmetry, presumably as they contain more information to constrain the physical model. The performance of this model makes it a practical and cost-effective tool to determine some of the BLR properties of a large sample of low and high redshift AGN, for which monitoring data are not available.
[14]  oai:arXiv.org:2001.06551  [pdf] - 2036338
Optical follow-up of gravitational wave triggers with DECam during the first two LIGO/VIRGO observing runs
Comments: 15 pages, 7 figures, submitted to Astronomy and Computing
Submitted: 2020-01-17, last modified: 2020-01-22
Gravitational wave (GW) events detectable by LIGO and Virgo have several possible progenitors, including black hole mergers, neutron star mergers, black hole--neutron star mergers, supernovae, and cosmic string cusps. A subset of GW events are expected to produce electromagnetic (EM) emission that, once detected, will provide complementary information about their astrophysical context. To that end, the LIGO--Virgo Collaboration (LVC) sends GW candidate alerts to the astronomical community so that searches for their EM counterparts can be pursued. The DESGW group, consisting of members of the Dark Energy Survey (DES), the LVC, and other members of the astronomical community, uses the Dark Energy Camera (DECam) to perform a search and discovery program for optical signatures of LVC GW events. DESGW aims to use a sample of GW events as standard sirens for cosmology. Due to the short decay timescale of the expected EM counterparts and the need to quickly eliminate survey areas with no counterpart candidates, it is critical to complete the initial analysis of each night's images as quickly as possible. We discuss our search area determination, imaging pipeline, and candidate selection processes. We review results from the DESGW program during the first two LIGO--Virgo observing campaigns and introduce other science applications that our pipeline enables.
[15]  oai:arXiv.org:1910.12943  [pdf] - 2054054
Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum
Comments: Accepted to ApJ
Submitted: 2019-10-28, last modified: 2020-01-04
We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.
[16]  oai:arXiv.org:1907.01981  [pdf] - 2042157
The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science case and survey design
Comments: 34 pages, accepted by PASP (modified from prior version to address referee's and coauthor comments). (v2) Minor fixes to author list
Submitted: 2019-07-03, last modified: 2019-12-30
The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.
[17]  oai:arXiv.org:1911.07688  [pdf] - 1999764
emcee v3: A Python ensemble sampling toolkit for affine-invariant MCMC
Comments: Published in the Journal for Open Source Software
Submitted: 2019-11-18
emcee is a Python library implementing a class of affine-invariant ensemble samplers for Markov chain Monte Carlo (MCMC). This package has been widely applied to probabilistic modeling problems in astrophysics where it was originally published, with some applications in other fields. When it was first released in 2012, the interface implemented in emcee was fundamentally different from the MCMC libraries that were popular at the time, such as PyMC, because it was specifically designed to work with "black box" models instead of structured graphical models. This has been a popular interface for applications in astrophysics because it is often non-trivial to implement realistic physics within the modeling frameworks required by other libraries. Since emcee's release, other libraries have been developed with similar interfaces, such as dynesty (Speagle 2019). The version 3.0 release of emcee is the first major release of the library in about 6 years and it includes a full re-write of the computational backend, several commonly requested features, and a set of new "move" implementations.
[18]  oai:arXiv.org:1909.06393  [pdf] - 2025811
Two years of non-thermal emission from the binary neutron star merger GW170817: rapid fading of the jet afterglow and first constraints on the kilonova fastest ejecta
Comments: version accepted for publication in ApJL, 13 pages, 6 figures
Submitted: 2019-09-13, last modified: 2019-11-08
We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_{\nu}\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_{\nu}\propto \nu^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $\theta_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $\Gamma_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $\theta_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $\epsilon_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $\epsilon_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5\sigma$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (\Gamma\beta)^{-\alpha}$ into the environment, finding that shallow stratification indexes $\alpha\le6$ are disfavored.
[19]  oai:arXiv.org:1909.11733  [pdf] - 2030542
Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array
Comments: Submitted to ApJ
Submitted: 2019-09-25, last modified: 2019-10-29
We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.
[20]  oai:arXiv.org:1908.08046  [pdf] - 1966983
The Optical Afterglow of GW170817: An Off-axis Structured Jet and Deep Constraints on a Globular Cluster Origin
Comments: 11 pages, 6 figures, 2 tables. Accepted to ApJ Letters
Submitted: 2019-08-21
We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 110-362$ days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously indicated by radio and X-ray data. Combined with contemporaneous radio and X-ray observations, we find no spectral evolution, with a weighted average spectral index of $\langle \beta \rangle = -0.583 \pm 0.013$, demonstrating that no synchrotron break frequencies evolve between the radio and X-ray bands over these timescales. We find that an extrapolation of the post-peak temporal slope of GW170817 to the luminosities of cosmological short GRBs matches their observed jet break times, suggesting that their explosion properties are similar, and that the primary difference in GW170817 is viewing angle. Additionally, we place a deep limit on the luminosity and mass of an underlying globular cluster of $L \lesssim 6.7 \times 10^{3}\,L_{\odot}$, or $M \lesssim 1.3 \times 10^{4}\,M_{\odot}$, at least 4 standard deviations below the peak of the globular cluster mass function of the host galaxy, NGC4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a globular cluster. As highlighted here, HST (and soon JWST) enables critical observations of the optical emission from neutron star merger jets and outflows.
[21]  oai:arXiv.org:1908.02383  [pdf] - 1931014
Electrical and electromagnetic co-simulations of the HERA Phase I receiver system including the effects of mutual coupling, and impact on the EoR window
Comments: 25 pages, 29 figures - Submitted to MNRAS
Submitted: 2019-08-06
The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.
[22]  oai:arXiv.org:1905.02186  [pdf] - 1919413
Follow-up of the Neutron Star Bearing Gravitational Wave Candidate Events S190425z and S190426c with MMT and SOAR
Comments: The first six authors contributed equally to this work
Submitted: 2019-05-06, last modified: 2019-07-18
On 2019 April 25.346 and 26.640 UT the LIGO and Virgo gravitational wave (GW) observatories announced the detection of the first candidate events in Observing Run 3 that contain at least one neutron star. S190425z is a likely binary neutron star (BNS) merger at $d_L = 156 \pm 41$ Mpc, while S190426c is possibly the first NS-BH merger ever detected, at $d_L = 377 \pm 100$ Mpc, although with marginal statistical significance. Here we report our optical follow-up observations for both events using the MMT 6.5-m telescope, as well as our spectroscopic follow-up of candidate counterparts (which turned out to be unrelated) with the 4.1-m SOAR telescope. We compare to publicly reported searches, explore the overall areal coverage and depth, and evaluate those in relation to the optical/NIR kilonova emission from the BNS merger GW170817, to theoretical kilonova models, and to short GRB afterglows. We find that for a GW170817-like kilonova, the partial volume covered spans up to about 40% for S190425z and 60% for S190426c. For an on-axis jet typical of short GRBs, the search effective volume is larger, but such a configuration is expected in at most a few percent of mergers. We further find that wide-field $\gamma$-ray and X-ray limits rule out luminous on-axis SGRBs, for a large fraction of the localization regions, although these searches are not sufficiently deep in the context of the $\gamma$-ray emission from GW170817 or off-axis SGRB afterglows. The results indicate that some optical follow-up searches are sufficiently deep for counterpart identification to about 300 Mpc, but that localizations better than 1000 deg$^2$ are likely essential.
[23]  oai:arXiv.org:1907.05383  [pdf] - 1914780
IDEAS: Immersive Dome Experiences for Accelerating Science
Comments: Astro2020 White Paper submission, 10 pages, 2 figures
Submitted: 2019-07-11
Astrophysics lies at the crossroads of big datasets (such as the Large Synoptic Survey Telescope and Gaia), open source software to visualize and interpret high dimensional datasets (such as Glue, WorldWide Telescope, and OpenSpace), and uniquely skilled software engineers who bridge data science and research fields. At the same time, more than 4,000 planetariums across the globe immerse millions of visitors in scientific data. We have identified the potential for critical synergy across data, software, hardware, locations, and content that -- if prioritized over the next decade -- will drive discovery in astronomical research. Planetariums can and should be used for the advancement of scientific research. Current facilities such as the Hayden Planetarium in New York City, Adler Planetarium in Chicago, Morrison Planetarium in San Francisco, the Iziko Planetarium and Digital Dome Research Consortium in Cape Town, and Visualization Center C in Norrkoping are already developing software which ingests catalogs of astronomical and multi-disciplinary data critical for exploration research primarily for the purpose of creating scientific storylines for the general public. We propose a transformative model whereby scientists become the audience and explorers in planetariums, utilizing software for their own investigative purposes. In this manner, research benefits from the authentic and unique experience of data immersion contained in an environment bathed in context and equipped for collaboration. Consequently, in this white paper we argue that over the next decade the research astronomy community should partner with planetariums to create visualization-based research opportunities for the field. Realizing this vision will require new investments in software and human capital.
[24]  oai:arXiv.org:1906.11595  [pdf] - 1912976
Variable dust emission by WC type Wolf-Rayet stars observed in the NEOWISE-R survey
Comments: accepted to MNRAS
Submitted: 2019-06-27
Photometry at 3.4 and 4.6 micron of 128 Population~I WC type Wolf-Rayet stars in the Galaxy and 12 in the Large Magellanic Cloud (LMC) observed in the WISE NEOWISE-R survey was searched for evidence of circumstellar dust emission and its variation. Infrared spectral energy distributions (SEDs) were assembled, making use of archival r, i, Z and Y photometry to determine reddening and stellar wind levels for the WC stars found in recent IR surveys and lacking optical photometry. From their SEDs, ten apparently non-variable stars were newly identified as dust makers, including three, WR 102-22, WR 110-10 and WR 124-10, having subtype earlier than WC8-9, the first such stars to show this phenomenon. The 11 stars found to show variable dust emission include six new episodic dust-makers, WR 47c, WR 75-11, WR 91-1, WR 122-14 and WR 125-1 in the Galaxy and HD 38030 in the LMC. Of previously known dust makers, NEOWISE-R photometry of WR 19 captured its rise to maximum in 2018, confirming the 10.1-y period, that of WR 125 the beginning of a new episode of dust formation suggesting a period near 28.3~y., while that of HD 36402 covered almost a whole period and forced revision of it to 5.1~y.
[25]  oai:arXiv.org:1906.11272  [pdf] - 1907615
Do Reverberation Mapping Analyses Provide an Accurate Picture of the Broad Line Region?
Comments: 22 pages, 24 figures
Submitted: 2019-06-26
Reverberation mapping (RM) is a powerful approach for determining the nature of the broad-line region (BLR) in active galactic nuclei. However, inferring physical BLR properties from an observed spectroscopic time series is a difficult inverse problem. Here, we present a blind test of two widely used RM methods: MEMEcho (developed by Horne) and CARAMEL (developed by Pancoast and collaborators). The test data are simulated spectroscopic time series that track the H$\alpha$ emission line response to an empirical continuum light curve. The underlying BLR model is a rotating, biconical accretion disc wind, and the synthetic spectra are generated via self-consistent ionization and radiative transfer simulations. We generate two mock data sets, representing Seyfert galaxies and QSOs. The Seyfert model produces a largely *negative* response, which neither method can recover. However, both fail $``gracefully''$, neither generating spurious results. For the QSO model both CARAMEL and expert interpretation of MEMEcho's output both capture the broadly annular, rotation-dominated nature of the line-forming region, though MEMEcho analysis overestimates its size by 50%, but CARAMEL is unable to distinguish between additional inflow and outflow components. Despite fitting individual spectra well, the CARAMEL velocity-delay maps and RMS line profiles are strongly inconsistent with the input data. Finally, since the H$\alpha$ line-forming region is rotation dominated, neither method recovers the disc wind nature of the underlying BLR model. Thus considerable care is required when interpreting the results of RM analyses in terms of physical models.
[26]  oai:arXiv.org:1802.04151  [pdf] - 1958034
The HERA-19 Commissioning Array: Direction Dependent Effects
Comments: 21 pages, 12 figures, submitted to ApJ
Submitted: 2018-02-07, last modified: 2019-06-10
Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.
[27]  oai:arXiv.org:1905.08260  [pdf] - 1966767
COSMOGRAIL XVIII: time delays of the quadruply lensed quasar WFI2033-4723
Comments: Submitted to Astronomy and Astrophysics
Submitted: 2019-05-20
We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.
[28]  oai:arXiv.org:1904.04320  [pdf] - 1864136
X-ray Studies of Exoplanets: A 2020 Decadal Survey White Paper
Comments: 7 pages with 2 figures. A 2020 Decadal Survey White Paper
Submitted: 2019-04-08
Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray imaging and spectroscopy provide powerful and unique windows into the high energy flux that an exoplanet experiences, and X-ray photons also serve as proxies for potentially transfigurative coronal mass ejections. Finally, if the host star is a bright enough X-ray source, transit measurements akin to those in the optical and infrared are possible and allow for direct characterization of the upper atmospheres of exoplanets. In this brief white paper, we discuss contributions to the study of exoplanets and their environs which can be made by X-ray data of increasingly high quality that are achievable in the next 10--15 years.
[29]  oai:arXiv.org:1903.12338  [pdf] - 1858405
High-Energy Photon and Particle Effects onExoplanet Atmospheres and Habitability
Comments: Astro2020 Decadal Survey Science White Paper
Submitted: 2019-03-28
It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections.
[30]  oai:arXiv.org:1903.11586  [pdf] - 1857308
Magnetism in the Brown Dwarf Regime
Comments: Decadal 2020 science white paper
Submitted: 2019-03-27
A suite of discoveries in the last two decades demonstrate that we are now at a point where incorporating magnetic behavior is key for advancing our ability to characterize substellar and planetary systems. The next decade heralds the exciting maturation of the now-burgeoning field of brown dwarf magnetism, and investing now in brown dwarf magnetism will provide a key platform for exploring exoplanetary magnetism and habitability beyond the solar system. We anticipate significant discoveries including: the nature of substellar and planetary magnetic dynamos, the characterization of exo-aurora physics and brown dwarf magnetospheric environments, and the role of satellites in manifestations of substellar magnetic activity. These efforts will require significant new observational capabilities at radio and near infrared wavelengths, dedicated long-term monitoring programs, and committed support for the theoretical modeling efforts underpinning the physical processes of the magnetic phenomena
[31]  oai:arXiv.org:1901.01540  [pdf] - 1905656
First measurement of the Hubble constant from a dark standard siren using the Dark Energy Survey galaxies and the LIGO/Virgo binary-black-hole merger GW170814
The DES Collaboration; the LIGO Scientific Collaboration; the Virgo Collaboration; Soares-Santos, M.; Palmese, A.; Hartley, W.; Annis, J.; Garcia-Bellido, J.; Lahav, O.; Doctor, Z.; Fishbach, M.; Holz, D. E.; Lin, H.; Pereira, M. E. S.; Garcia, A.; Herner, K.; Kessler, R.; Peiris, H. V.; Sako, M.; Allam, S.; Brout, D.; Rosell, A. Carnero; Chen, H. Y.; Conselice, C.; deRose, J.; deVicente, J.; Diehl, H. T.; Gill, M. S. S.; Gschwend, J.; Sevilla-Noarbe, I.; Tucker, D. L.; Wechsler, R.; Berger, E.; Cowperthwaite, P. S.; Metzger, B. D.; Williams, P. K. G.; Abbott, T. M. C.; Abdalla, F. B.; Avila, S.; Bechtol, K.; Bertin, E.; Brooks, D.; Buckley-Geer, E.; Burke, D. L.; Kind, M. Carrasco; Carretero, J.; Castander, F. J.; Crocce, M.; Cunha, C. E.; D'Andrea, C. B.; da Costa, L. N.; Davis, C.; Desai, S.; Doel, P.; Drlica-Wagner, A.; Eifler, T. F.; Evrard, A. E.; Flaugher, B.; Fosalba, P.; Frieman, J.; Gaztanaga, E.; Gerdes, D. W.; Gruen, D.; Gruendl, R. A.; Gutierrez, G.; Hollowood, D. L.; Hoyle, B.; James, D. J.; Jeltema, T.; Kuehn, K.; Kuropatkin, N.; Li, T. S.; Lima, M.; Maia, M. A. G.; Marshall, J. L.; Menanteau, F.; Miquel, R.; Neilsen, E.; Ogando, R. L. C.; Plazas, A. A.; Romer, A. K.; Roodman, A.; Sanchez, E.; Scarpine, V.; Schindler, R.; Schubnell, M.; Serrano, S.; Smith, M.; Smith, R. C.; Sobreira, F.; Suchyta, E.; Swanson, M. E. C.; Tarle, G.; Thomas, R. C.; Walker, A. R.; Wester, W.; Zuntz, J.; Abbott, B. P.; Abbott, R.; Abbott, T. D.; Abraham, S.; Acernese, F.; Ackley, K.; Adams, C.; Adhikari, R. X.; Adya, V. B.; Affeldt, C.; Agathos, M.; Agatsuma, K.; Aggarwal, N.; Aguiar, O. D.; Aiello, L.; Ain, A.; Ajith, P.; Allen, G.; Allocca, A.; Aloy, M. A.; Altin, P. A.; Amato, A.; Ananyeva, A.; Anderson, S. B.; Anderson, W. G.; Angelova, S. V.; Appert, S.; Arai, K.; Araya, M. C.; Areeda, J. S.; Arene, M.; Ascenzi, S.; Ashton, G.; Aston, S. M.; Astone, P.; Aubin, F.; Aufmuth, P.; AultONeal, K.; Austin, C.; Avendano, V.; Avila-Alvarez, A.; Babak, S.; Bacon, P.; Badaracco, F.; Bader, M. K. M.; Bae, S.; Baker, P. T.; Baldaccini, F.; Ballardin, G.; Ballmer, S. W.; Banagiri, S.; Barayoga, J. C.; Barclay, S. E.; Barish, B. C.; Barker, D.; Barkett, K.; Barnum, S.; Barone, F.; Barr, B.; Barsotti, L.; Barsuglia, M.; Barta, D.; Bartlett, J.; Bartos, I.; Bassiri, R.; Basti, A.; Bawaj, M.; Bayley, J. C.; Bazzan, M.; Becsy, B.; Bejger, M.; Bell, A. S.; Beniwal, D.; Bergmann, G.; Bernuzzi, S.; Bero, J. J.; Berry, C. P. L.; Bersanetti, D.; Bertolini, A.; Betzwieser, J.; Bhandare, R.; Bidler, J.; Bilenko, I. A.; Bilgili, S. A.; Billingsley, G.; Birch, J.; Birney, R.; Birnholtz, O.; Biscans, S.; Biscoveanu, S.; Bisht, A.; Bitossi, M.; Blackburn, J. K.; Blair, C. D.; Blair, D. G.; Blair, R. M.; Bloemen, S.; Bode, N.; Boer, M.; Boetzel, Y.; Bogaert, G.; Bondu, F.; Bonilla, E.; Bonnand, R.; Booker, P.; Boom, B. A.; Booth, C. D.; Bork, R.; Boschi, V.; Bose, S.; Bossie, K.; Bossilkov, V.; Bosveld, J.; Bouffanais, Y.; Bozzi, A.; Bradaschia, C.; Brady, P. R.; Bramley, A.; Branchesi, M.; Brau, J. E.; Briant, T.; Briggs, J. H.; Brighenti, F.; Brillet, A.; Brinkmann, M.; Brockill, P.; Brooks, A. F.; Brown, D. D.; Brunett, S.; Buikema, A.; Bulik, T.; Bulten, H. J.; Buonanno, A.; Buskulic, D.; Buy, C.; Byer, R. L.; Cabero, M.; Cadonati, L.; Cagnoli, G.; Cahillane, C.; Bustillo, J. Calderon; Callister, T. A.; Calloni, E.; Camp, J. B.; Campbell, W. A.; Cannon, K. C.; Cao, H.; Cao, J.; Capocasa, E.; Carbognani, F.; Caride, S.; Carney, M. F.; Carullo, G.; Diaz, J. Casanueva; Casentini, C.; Caudill, S.; Cavaglia, M.; Cavalieri, R.; Cella, G.; Cerda-Duran, P.; Cerretani, G.; Cesarini, E.; Chaibi, O.; Chakravarti, K.; Chamberlin, S. J.; Chan, M.; Chao, S.; Charlton, P.; Chase, E. A.; Chassande-Mottin, E.; Chatterjee, D.; Chaturvedi, M.; Chatziioannou, K.; Cheeseboro, B. D.; Chen, H. Y.; Chen, X.; Chen, Y.; Cheng, H. -P.; Cheong, C. K.; Chia, H. Y.; Chincarini, A.; Chiummo, A.; Cho, G.; Cho, H. S.; Cho, M.; Christensen, N.; Chu, Q.; Chua, S.; Chung, K. W.; Chung, S.; Ciani, G.; Ciobanu, A. A.; Ciolfi, R.; Cipriano, F.; Cirone, A.; Clara, F.; Clark, J. A.; Clearwater, P.; Cleva, F.; Cocchieri, C.; Coccia, E.; Cohadon, P. -F.; Colgan, R.; Colleoni, M.; Collette, C. G.; Collins, C.; Cominsky, L. R.; Constancio, M.; Conti, L.; Cooper, S. J.; Corban, P.; Corbitt, T. R.; Cordero-Carrion, I.; Corley, K. R.; Cornish, N.; Corsi, A.; Cortese, S.; Costa, C. A.; Cotesta, R.; Coughlin, M. W.; Coughlin, S. B.; Coulon, J. -P.; Countryman, S. T.; Couvares, P.; Covas, P. B.; Cowan, E. E.; Coward, D. M.; Cowart, M. J.; Coyne, D. C.; Coyne, R.; Creighton, J. D. E.; Creighton, T. D.; Cripe, J.; Croquette, M.; Crowder, S. G.; Cullen, T. J.; Cumming, A.; Cunningham, L.; Cuoco, E.; Canton, T. Dal; Dalya, G.; Danilishin, S. L.; D'Antonio, S.; Danzmann, K.; Dasgupta, A.; Costa, C. F. Da Silva; Datrier, L. E. H.; Dattilo, V.; Dave, I.; Davis, D.; Daw, E. J.; DeBra, D.; Deenadayalan, M.; Degallaix, J.; De Laurentis, M.; Deleglise, S.; Del Pozzo, W.; DeMarchi, L. M.; Demos, N.; Dent, T.; De Pietri, R.; Derby, J.; De Rosa, R.; De Rossi, C.; DeSalvo, R.; de Varona, O.; Dhurandhar, S.; Diaz, M. C.; Dietrich, T.; Di Fiore, L.; Di Giovanni, M.; Di Girolamo, T.; Di Lieto, A.; Ding, B.; Di Pace, S.; Di Palma, I.; Di Renzo, F.; Dmitriev, A.; Donovan, F.; Dooley, K. L.; Doravari, S.; Dorrington, I.; Downes, T. P.; Drago, M.; Driggers, J. C.; Du, Z.; Dupej, P.; Dwyer, S. E.; Easter, P. J.; Edo, T. B.; Edwards, M. C.; Effler, A.; Ehrens, P.; Eichholz, J.; Eikenberry, S. S.; Eisenmann, M.; Eisenstein, R. A.; Estelles, H.; Estevez, D.; Etienne, Z. B.; Etzel, T.; Evans, M.; Evans, T. M.; Fafone, V.; Fair, H.; Fairhurst, S.; Fan, X.; Farinon, S.; Farr, B.; Farr, W. M.; Fauchon-Jones, E. J.; Favata, M.; Fays, M.; Fazio, M.; Fee, C.; Feicht, J.; Fejer, M. M.; Feng, F.; Fernandez-Galiana, A.; Ferrante, I.; Ferreira, E. C.; Ferreira, T. A.; Ferrini, F.; Fidecaro, F.; Fiori, I.; Fiorucci, D.; Fisher, R. P.; Fishner, J. M.; Fitz-Axen, M.; Flaminio, R.; Fletcher, M.; Flynn, E.; Fong, H.; Font, J. A.; Forsyth, P. W. F.; Fournier, J. -D.; Frasca, S.; Frasconi, F.; Frei, Z.; Freise, A.; Frey, R.; Fritschel, P.; Frolov, V. V.; Fulda, P.; Fyffe, M.; Gabbard, H. A.; Gadre, B. U.; Gaebel, S. M.; Gair, J. R.; Gammaitoni, L.; Ganija, M. R.; Gaonkar, S. G.; Garcia, A.; Garcia-Quiros, C.; Garufi, F.; Gateley, B.; Gaudio, S.; Gaur, G.; Gayathri, V.; Gemme, G.; Genin, E.; Gennai, A.; George, D.; George, J.; Gergely, L.; Germain, V.; Ghonge, S.; Ghosh, Abhirup; Ghosh, Archisman; Ghosh, S.; Giacomazzo, B.; Giaime, J. A.; Giardina, K. D.; Giazotto, A.; Gill, K.; Giordano, G.; Glover, L.; Godwin, P.; Goetz, E.; Goetz, R.; Goncharov, B.; Gonzalez, G.; Castro, J. M. Gonzalez; Gopakumar, A.; Gorodetsky, M. L.; Gossan, S. E.; Gosselin, M.; Gouaty, R.; Grado, A.; Graef, C.; Granata, M.; Grant, A.; Gras, S.; Grassia, P.; Gray, C.; Gray, R.; Greco, G.; Green, A. C.; Green, R.; Gretarsson, E. M.; Groot, P.; Grote, H.; Grunewald, S.; Guidi, G. M.; Gulati, H. K.; Guo, Y.; Gupta, A.; Gupta, M. K.; Gustafson, E. K.; Gustafson, R.; Haegel, L.; Halim, O.; Hall, B. R.; Hall, E. D.; Hamilton, E. Z.; Hammond, G.; Haney, M.; Hanke, M. M.; Hanks, J.; Hanna, C.; Hannuksela, O. A.; Hanson, J.; Hardwick, T.; Haris, K.; Harms, J.; Harry, G. M.; Harry, I. W.; Haster, C. -J.; Haughian, K.; Hayes, F. J.; Healy, J.; Heidmann, A.; Heintze, M. C.; Heitmann, H.; Hemming, G.; Hendry, M.; Heng, I. S.; Hennig, J.; Heptonstall, A. W.; Vivanco, Francisco Hernandez; Heurs, M.; Hild, S.; Hinderer, T.; Hoak, D.; Hochheim, S.; Hofman, D.; Holgado, A. M.; Holland, N. A.; Holt, K.; Hopkins, P.; Horst, C.; Hough, J.; Howell, E. J.; Hoy, C. G.; Hreibi, A.; Huerta, E. A.; Hughey, B.; Hulko, M.; Husa, S.; Huttner, S. H.; Huynh-Dinh, T.; Idzkowski, B.; Iess, A.; Ingram, C.; Inta, R.; Intini, G.; Irwin, B.; Isa, H. N.; Isac, J. -M.; Isi, M.; Iyer, B. R.; Izumi, K.; Jacqmin, T.; Jadhav, S. J.; Jani, K.; Janthalur, N. N.; Jaranowski, P.; Jenkins, A. C.; Jiang, J.; Johnson, D. S.; Jones, A. W.; Jones, D. I.; Jones, R.; Jonker, R. J. G.; Ju, L.; Junker, J.; Kalaghatgi, C. V.; Kalogera, V.; Kamai, B.; Kandhasamy, S.; Kang, G.; Kanner, J. B.; Kapadia, S. J.; Karki, S.; Karvinen, K. S.; Kashyap, R.; Kasprzack, M.; Katsanevas, S.; Katsavounidis, E.; Katzman, W.; Kaufer, S.; Kawabe, K.; Keerthana, N. V.; Kefelian, F.; Keitel, D.; Kennedy, R.; Key, J. S.; Khalili, F. Y.; Khan, H.; Khan, I.; Khan, S.; Khan, Z.; Khazanov, E. A.; Khursheed, M.; Kijbunchoo, N.; Kim, Chunglee; Kim, J. C.; Kim, K.; Kim, W.; Kim, W. S.; Kim, Y. -M.; Kimball, C.; King, E. J.; King, P. J.; Kinley-Hanlon, M.; Kirchhoff, R.; Kissel, J. S.; Kleybolte, L.; Klika, J. H.; Klimenko, S.; Knowles, T. D.; Koch, P.; Koehlenbeck, S. M.; Koekoek, G.; Koley, S.; Kondrashov, V.; Kontos, A.; Koper, N.; Korobko, M.; Korth, W. Z.; Kowalska, I.; Kozak, D. B.; Kringel, V.; Krishnendu, N.; Krolak, A.; Kuehn, G.; Kumar, A.; Kumar, P.; Kumar, R.; Kumar, S.; Kuo, L.; Kutynia, A.; Kwang, S.; Lackey, B. D.; Lai, K. H.; Lam, T. L.; Landry, M.; Lane, B. B.; Lang, R. N.; Lange, J.; Lantz, B.; Lanza, R. K.; Lasky, P. D.; Laxen, M.; Lazzarini, A.; Lazzaro, C.; Leaci, P.; Leavey, S.; Lecoeuche, Y. K.; Lee, C. H.; Lee, H. K.; Lee, H. M.; Lee, H. W.; Lee, J.; Lee, K.; Lehmann, J.; Lenon, A.; Letendre, N.; Levin, Y.; Li, J.; Li, K. J. L.; Li, T. G. F.; Li, X.; Lin, F.; Linde, F.; Linker, S. D.; Littenberg, T. B.; Liu, J.; Liu, X.; Lo, R. K. L.; Lockerbie, N. A.; London, L. T.; Longo, A.; Lorenzini, M.; Loriette, V.; Lormand, M.; Losurdo, G.; Lough, J. D.; Lousto, C. O.; Lovelace, G.; Lower, M. E.; Luck, H.; Lumaca, D.; Lundgren, A. P.; Lynch, R.; Ma, Y.; Macas, R.; Macfoy, S.; MacInnis, M.; Macleod, D. M.; Macquet, A.; Hernandez, I. Magana; Magana-Sandoval, F.; Zertuche, L. Magana; Magee, R. M.; Majorana, E.; Maksimovic, I.; Malik, A.; Man, N.; Mandic, V.; Mangano, V.; Mansell, G. L.; Manske, M.; Mantovani, M.; Marchesoni, F.; Marion, F.; Marka, S.; Marka, Z.; Markakis, C.; Markosyan, A. S.; Markowitz, A.; Maros, E.; Marquina, A.; Marsat, S.; Martelli, F.; Martin, I. W.; Martin, R. M.; Martynov, D. V.; Mason, K.; Massera, E.; Masserot, A.; Massinger, T. J.; Masso-Reid, M.; Mastrogiovanni, S.; Matas, A.; Matichard, F.; Matone, L.; Mavalvala, N.; Mazumder, N.; McCann, J. J.; McCarthy, R.; McClelland, D. E.; McCormick, S.; McCuller, L.; McGuire, S. C.; McIver, J.; McManus, D. J.; McRae, T.; McWilliams, S. T.; Meacher, D.; Meadors, G. D.; Mehmet, M.; Mehta, A. K.; Meidam, J.; Melatos, A.; Mendell, G.; Mercer, R. A.; Mereni, L.; Merilh, E. L.; Merzougui, M.; Meshkov, S.; Messenger, C.; Messick, C.; Metzdorff, R.; Meyers, P. M.; Miao, H.; Michel, C.; Middleton, H.; Mikhailov, E. E.; Milano, L.; Miller, A. L.; Miller, A.; Millhouse, M.; Mills, J. C.; Milovich-Goff, M. C.; Minazzoli, O.; Minenkov, Y.; Mishkin, A.; Mishra, C.; Mistry, T.; Mitra, S.; Mitrofanov, V. P.; Mitselmakher, G.; Mittleman, R.; Mo, G.; Moffa, D.; Mogushi, K.; Mohapatra, S. R. P.; Montani, M.; Moore, C. J.; Moraru, D.; Moreno, G.; Morisaki, S.; Mours, B.; Mow-Lowry, C. M.; Mukherjee, Arunava; Mukherjee, D.; Mukherjee, S.; Mukund, N.; Mullavey, A.; Munch, J.; Muniz, E. A.; Muratore, M.; Murray, P. G.; Nardecchia, I.; Naticchioni, L.; Nayak, R. K.; Neilson, J.; Nelemans, G.; Nelson, T. J. N.; Nery, M.; Neunzert, A.; Ng, K. Y.; Ng, S.; Nguyen, P.; Nichols, D.; Nissanke, S.; Nocera, F.; North, C.; Nuttall, L. K.; Obergaulinger, M.; Oberling, J.; O'Brien, B. D.; O'Dea, G. D.; Ogin, G. H.; Oh, J. J.; Oh, S. H.; Ohme, F.; Ohta, H.; Okada, M. A.; Oliver, M.; Oppermann, P.; Oram, Richard J.; O'Reilly, B.; Ormiston, R. G.; Ortega, L. F.; O'Shaughnessy, R.; Ossokine, S.; Ottaway, D. J.; Overmier, H.; Owen, B. J.; Pace, A. E.; Pagano, G.; Page, M. A.; Pai, A.; Pai, S. A.; Palamos, J. R.; Palashov, O.; Palomba, C.; Pal-Singh, A.; Pan, Huang-Wei; Pang, B.; Pang, P. T. H.; Pankow, C.; Pannarale, F.; Pant, B. C.; Paoletti, F.; Paoli, A.; Parida, A.; Parker, W.; Pascucci, D.; Pasqualetti, A.; Passaquieti, R.; Passuello, D.; Patil, M.; Patricelli, B.; Pearlstone, B. L.; Pedersen, C.; Pedraza, M.; Pedurand, R.; Pele, A.; Penn, S.; Perez, C. J.; Perreca, A.; Pfeiffer, H. P.; Phelps, M.; Phukon, K. S.; Piccinni, O. J.; Pichot, M.; Piergiovanni, F.; Pillant, G.; Pinard, L.; Pirello, M.; Pitkin, M.; Poggiani, R.; Pong, D. Y. T.; Ponrathnam, S.; Popolizio, P.; Porter, E. K.; Powell, J.; Prajapati, A. K.; Prasad, J.; Prasai, K.; Prasanna, R.; Pratten, G.; Prestegard, T.; Privitera, S.; Prodi, G. A.; Prokhorov, L. G.; Puncken, O.; Punturo, M.; Puppo, P.; Purrer, M.; Qi, H.; Quetschke, V.; Quinonez, P. J.; Quintero, E. A.; Quitzow-James, R.; Radkins, H.; Radulescu, N.; Raffai, P.; Raja, S.; Rajan, C.; Rajbhandari, B.; Rakhmanov, M.; Ramirez, K. E.; Ramos-Buades, A.; Rana, Javed; Rao, K.; Rapagnani, P.; Raymond, V.; Razzano, M.; Read, J.; Regimbau, T.; Rei, L.; Reid, S.; Reitze, D. H.; Ren, W.; Ricci, F.; Richardson, C. J.; Richardson, J. W.; Ricker, P. M.; Riles, K.; Rizzo, M.; Robertson, N. A.; Robie, R.; Rocchi, A.; Rolland, L.; Rollins, J. G.; Roma, V. J.; Romanelli, M.; Romano, R.; Romel, C. L.; Romie, J. H.; Rose, K.; Rosinska, D.; Rosofsky, S. G.; Ross, M. P.; Rowan, S.; Rudiger, A.; Ruggi, P.; Rutins, G.; Ryan, K.; Sachdev, S.; Sadecki, T.; Sakellariadou, M.; Salconi, L.; Saleem, M.; Samajdar, A.; Sammut, L.; Sanchez, E. J.; Sanchez, L. E.; Sanchis-Gual, N.; Sandberg, V.; Sanders, J. R.; Santiago, K. A.; Sarin, N.; Sassolas, B.; Saulson, P. R.; Sauter, O.; Savage, R. L.; Schale, P.; Scheel, M.; Scheuer, J.; Schmidt, P.; Schnabel, R.; Schofield, R. M. S.; Schonbeck, A.; Schreiber, E.; Schulte, B. W.; Schutz, B. F.; Schwalbe, S. G.; Scott, J.; Scott, S. M.; Seidel, E.; Sellers, D.; Sengupta, A. S.; Sennett, N.; Sentenac, D.; Sequino, V.; Sergeev, A.; Shaddock, D. A.; Shaffer, T.; Shahriar, M. S.; Shaner, M. B.; Shao, L.; Sharma, P.; Shawhan, P.; Shen, H.; Shink, R.; Shoemaker, D. H.; Shoemaker, D. M.; ShyamSundar, S.; Siellez, K.; Sieniawska, M.; Sigg, D.; Silva, A. D.; Singer, L. P.; Singh, N.; Singhal, A.; Sintes, A. M.; Sitmukhambetov, S.; Skliris, V.; Slagmolen, B. J. J.; Slaven-Blair, T. J.; Smith, J. R.; Smith, R. J. E.; Somala, S.; Son, E. J.; Sorazu, B.; Sorrentino, F.; Souradeep, T.; Sowell, E.; Spencer, A. P.; Srivastava, A. K.; Srivastava, V.; Staats, K.; Stachie, C.; Standke, M.; Steer, D. A.; Steinke, M.; Steinlechner, J.; Steinlechner, S.; Steinmeyer, D.; Stevenson, S. P.; Stocks, D.; Stone, R.; Stops, D. J.; Strain, K. A.; Stratta, G.; Strigin, S. E.; Strunk, A.; Sturani, R.; Stuver, A. L.; Sudhir, V.; Summerscales, T. Z.; Sun, L.; Sunil, S.; Sur, A.; Suresh, J.; Sutton, P. J.; Swinkels, B. L.; Szczepanczyk, M. J.; Tacca, M.; Tait, S. C.; Talbot, C.; Talukder, D.; Tanner, D. B.; Tapai, M.; Taracchini, A.; Tasson, J. D.; Taylor, R.; Thies, F.; Thomas, M.; Thomas, P.; Thondapu, S. R.; Thorne, K. A.; Thrane, E.; Tiwari, Shubhanshu; Tiwari, Srishti; Tiwari, V.; Toland, K.; Tonelli, M.; Tornasi, Z.; Torres-Forne, A.; Torrie, C. I.; Toyra, D.; Travasso, F.; Traylor, G.; Tringali, M. C.; Trovato, A.; Trozzo, L.; Trudeau, R.; Tsang, K. W.; Tse, M.; Tso, R.; Tsukada, L.; Tsuna, D.; Tuyenbayev, D.; Ueno, K.; Ugolini, D.; Unnikrishnan, C. S.; Urban, A. L.; Usman, S. A.; Vahlbruch, H.; Vajente, G.; Valdes, G.; van Bakel, N.; van Beuzekom, M.; Brand, J. F. J. van den; Broeck, C. Van Den; Vander-Hyde, D. C.; van Heijningen, J. V.; van der Schaaf, L.; van Veggel, A. A.; Vardaro, M.; Varma, V.; Vass, S.; Vasuth, M.; Vecchio, A.; Vedovato, G.; Veitch, J.; Veitch, P. J.; Venkateswara, K.; Venugopalan, G.; Verkindt, D.; Vetrano, F.; Vicere, A.; Viets, A. D.; Vine, D. J.; Vinet, J. -Y.; Vitale, S.; Vo, T.; Vocca, H.; Vorvick, C.; Vyatchanin, S. P.; Wade, A. R.; Wade, L. E.; Wade, M.; Walet, R.; Walker, M.; Wallace, L.; Walsh, S.; Wang, G.; Wang, H.; Wang, J. Z.; Wang, W. H.; Wang, Y. F.; Ward, R. L.; Warden, Z. A.; Warner, J.; Was, M.; Watchi, J.; Weaver, B.; Wei, L. -W.; Weinert, M.; Weinstein, A. J.; Weiss, R.; Wellmann, F.; Wen, L.; Wessel, E. K.; Wessels, P.; Westhouse, J. W.; Wette, K.; Whelan, J. T.; Whiting, B. F.; Whittle, C.; Wilken, D. M.; Williams, D.; Williamson, A. R.; Willis, J. L.; Willke, B.; Wimmer, M. H.; Winkler, W.; Wipf, C. C.; Wittel, H.; Woan, G.; Woehler, J.; Wofford, J. K.; Worden, J.; Wright, J. L.; Wu, D. S.; Wysocki, D. M.; Xiao, L.; Yamamoto, H.; Yancey, C. C.; Yang, L.; Yap, M. J.; Yazback, M.; Yeeles, D. W.; Yu, Hang; Yu, Haocun; Yuen, S. H. R.; Yvert, M.; Zadrozny, A. K.; Zanolin, M.; Zelenova, T.; Zendri, J. -P.; Zevin, M.; Zhang, J.; Zhang, L.; Zhang, T.; Zhao, C.; Zhou, M.; Zhou, Z.; Zhu, X. J.; Zimmerman, A.; Zucker, M. E.; Zweizig, J.
Comments: 18 pages, 2 Figures, submitted to ApJL
Submitted: 2019-01-06, last modified: 2019-03-22
We present a multi-messenger measurement of the Hubble constant H_0 using the binary-black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the LIGO/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DES Year 3 data. Black-hole mergers such as GW170814 are expected to lack bright electromagnetic emission to uniquely identify their host galaxies and build an object-by-object Hubble diagram. However, they are suitable for a statistical measurement, provided that a galaxy catalog of adequate depth and redshift completion is available. Here we present the first Hubble parameter measurement using a black-hole merger. Our analysis results in $H_0 = 75.2^{+39.5}_{-32.4}~{\rm km~s^{-1}~Mpc^{-1}}$, which is consistent with both SN Ia and CMB measurements of the Hubble constant. The quoted 68% credible region comprises 60% of the uniform prior range [20,140] ${\rm km~s^{-1}~Mpc^{-1}}$, and it depends on the assumed prior range. If we take a broader prior of [10,220] ${\rm km~s^{-1}~Mpc^{-1}}$, we find $H_0 = 78^{+ 96}_{-24}~{\rm km~s^{-1}~Mpc^{-1}}$ ($57\%$ of the prior range). Although a weak constraint on the Hubble constant from a single event is expected using the dark siren method, a multifold increase in the LVC event rate is anticipated in the coming years and combinations of many sirens will lead to improved constraints on $H_0$.
[32]  oai:arXiv.org:1903.04553  [pdf] - 1847028
Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s
Foley, R. J.; Alexander, K. D.; Andreoni, I.; Arcavi, I.; Auchettl, K.; Barnes, J.; Baym, G.; Bellm, E. C.; Beloborodov, A. M.; Blagorodnova, N.; Blakeslee, J. P.; Brady, P. R.; Branchesi, M.; Brown, J. S.; Butler, N.; Cantiello, M.; Chornock, R.; Cook, D. O.; Cooke, J.; Coppejans, D. L.; Corsi, A.; Couch, S. M.; Coughlin, M. W.; Coulter, D. A.; Cowperthwaite, P. S.; Dietrich, T.; Dimitriadis, G.; Drout, M. R.; Elias, J. H.; Farr, B.; Fernandez, R.; Filippenko, A. V.; Fong, W.; Fragos, T.; Frail, D. A.; Freedman, W. L.; Fryer, C. L.; Golkhou, V. Z.; Hiramatsu, D.; Hjorth, J.; Horesh, A.; Hosseinzadeh, G.; Hotokezaka, K.; Howell, D. A.; Hung, T.; Jones, D. O.; Kalogera, V.; Kasen, D.; Kerzendorf, W. E.; Kilpatrick, C. D.; Kirshner, R. P.; Krisciunas, K.; Lattimer, J. M.; Lazzati, D.; Levan, A. J.; MacFadyen, A. I.; Maeda, K.; Mandel, I.; Mandel, K. S.; Margalit, B.; Margutti, R.; McIver, J.; Metzger, B. D.; Mooley, K.; Moriya, T.; Murguia-Berthier, A.; Narayan, G.; Nicholl, M.; Nissanke, S.; Nomoto, K.; O'Meara, J. M.; O'Shaughnessy, R.; O'Connor, E.; Palmese, A.; Pan, Y. -C.; Pankow, C.; Paterson, K.; Perley, D. A.; Perna, R.; Piro, A. L.; Pritchard, T. A.; Quataert, E.; Radice, D.; Ramirez-Ruiz, E.; Reddy, S.; Rest, A.; Riess, A. G.; Rodriguez, C. L.; Rojas-Bravo, C.; Rossi, E. M.; Rosswog, S.; Ruiz, M.; Shapiro, S. L.; Shoemaker, D. H.; Siebert, M. R.; Siegel, D. M.; Siellez, K.; Smith, N.; Soares-Santos, M.; Suntzeff, N. B.; Surman, R.; Tanaka, M.; Tanvir, N. R.; Terreran, G.; Valenti, S.; Villar, V. A.; Wang, L.; Webb, S. A.; Wheeler, J. C.; Williams, P. K. G.; Woosley, S.; Zaldarriaga, M.; Zevin, M.
Comments: submitted to Astro2020
Submitted: 2019-03-11
As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.
[33]  oai:arXiv.org:1903.04691  [pdf] - 1847063
Radio Time-Domain Signatures of Magnetar Birth
Comments: Science white paper for Astro2020
Submitted: 2019-03-11
The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques.
[34]  oai:arXiv.org:1902.08244  [pdf] - 1916817
Optimizing Sparse RFI Prediction using Deep Learning
Comments: 11 pages, 7 figures
Submitted: 2019-02-21
Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations.
[35]  oai:arXiv.org:1901.10479  [pdf] - 1878722
A Radio Source Coincident with the Superluminous Supernova PTF10hgi: Evidence for a Central Engine and an Analogue of the Repeating FRB121102?
Comments: 10 pages, 5 figures; submitted to ApJL
Submitted: 2019-01-29
We present the detection of an unresolved radio source coincident with the position of the Type I superluminous supernova (SLSN) PTF10hgi ($z=0.098$) about 7.5 years post-explosion, with a luminosity of $L_\nu(6\,\,{\rm GHz})\approx 1.1\times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$. This represents the first detection of radio emission coincident with a SLSN on any timescale. We investigate various scenarios for the origin of the radio emission: star formation activity, an active galactic nucleus, an off-axis jet, and a non-relativistic supernova blastwave. While any of these would be quite novel if confirmed, none appear likely when taken in context of the other properties of the host galaxy, previous radio observations of SLSNe, the sample of long gamma-ray bursts (LGRBs), and the general population of hydrogen-poor SNe. Instead, the radio emission is reminiscent of the quiescent radio source associated with the repeating FRB 121102, which has been argued to be powered by a magnetar born in a SLSN or LGRB explosion several decades ago. We show that such a central engine powered nebula is consistent with the age and luminosity of the radio source. Our directed search for FRBs from the location of PTF10hgi using 40 min of VLA phased-array data reveals no detections to a limit of $22$ mJy ($7\sigma$; 10 ms duration). We outline several follow-up observations that can conclusively establish the origin of the radio emission.
[36]  oai:arXiv.org:1901.05485  [pdf] - 1818897
Signs of accretion in the white dwarf + brown dwarf binary NLTT5306
Comments: 10 pages, 5 figures
Submitted: 2019-01-16
We present new XSHOOTER spectra of NLTT5306, a 0.44 $\pm$ 0.04\msun white dwarf in a short period (101\,min) binary system with a brown dwarf companion that is likely to have previously undergone common envelope evolution. We have confirmed the presence of H$\alpha$ emission and discovered Na I absorption associated with the white dwarf. These observations are indicative of accretion. Accretion is typically evidenced by high energy emission in the UV and X-ray regime. However our \textit{Swift} observations covering the full orbital period in three wavebands (uvw1, uvm2, uvw2) revealed no UV excess or modulation. We used the X-ray non-detection to put an upper limit on the accretion rate of 2$\times$10$^{-15}$\msun yr$^{-1}$. We compare NLTT5306 to similar accreting binaries with brown dwarf donors and suggest the inferred accretion rate could be from wind accretion or accretion from a debris/dust disk. The lack of evidence for a disk implies NLTT5306 is magnetically funnelling a weak wind from a potentially low gravity brown dwarf. The upper limit on the accretion rate suggests a magnetic field as low as 0.45\,kG would be sufficient to achieve this. If confirmed this would constitute the first detection of a brown dwarf wind and could provide useful constraints on mass loss rates.
[37]  oai:arXiv.org:1806.08017  [pdf] - 1814931
An Unexpectedly Small Emission Region Size Inferred from Strong High-Frequency Diffractive Scintillation in GRB 161219B
Comments: 14 pages, 7 figures, 1 table. Accepted to ApJ. The data used to create Figures 2-6 is available with the published article, or by request to kate.alexander@northwestern.edu
Submitted: 2018-06-20, last modified: 2019-01-09
We present Karl G. Jansky Very Large Array radio observations of the long gamma-ray burst GRB 161219B ($z=0.147$) spanning $1-37$ GHz. The data exhibit unusual behavior, including sharp spectral peaks and minutes-timescale large-amplitude variability centered at $20$ GHz and spanning the full frequency range. We attribute this behavior to scattering of the radio emission by the turbulent ionized Galactic interstellar medium (ISM), including both diffractive and refractive scintillation. However, the scintillation is much stronger than predicted by a model of the Galactic electron density distribution (NE2001); from the measured variability timescale and decorrelation bandwidth we infer a scattering measure of $SM\approx {(8-70)\times 10^{-4}}$ kpc m$^{-20/3}$ (up to ${25}$ times larger than predicted in NE2001) and a scattering screen distance of $d_{\rm scr}\approx {0.2-3}$ kpc. We infer an emission region size of $\theta_s \approx {0.9-4}$ $\mu$as ($\approx {(1-4)}\times 10^{16}$ cm) at $\approx4$ days, and find that prior to 8 days the source size is an order of magnitude smaller than model predictions for a uniformly illuminated disk or limb-brightened ring, indicating a slightly off-axis viewing angle or significant substructure in the emission region. Simultaneous multi-hour broadband radio observations of future GRB afterglows will allow us to characterize the scintillation more completely, and hence to probe the observer viewing angle, the evolution of the jet Lorentz factor, the structure of the afterglow emission regions, and ISM turbulence at high Galactic latitudes.
[38]  oai:arXiv.org:1812.07631  [pdf] - 1863928
K2 Ultracool Dwarfs Survey. V. High superflare rates on rapidly rotating late-M dwarfs
Comments: 10 pages, 10 figures, submitted to MNRAS
Submitted: 2018-12-18
We observed strong superflares (defined as flares with energy in excess of 10^33 erg) on three late-M dwarfs: 2MASS J08315742+2042213 (hereafter 2M0831+2042; M7 V), 2MASS J08371832+2050349 (hereafter 2M0837+2050; M8 V) and 2MASS J08312608+2244586 (hereafter 2M0831+2244; M9 V). 2M0831+2042 and 2M0837+2050 are members of the young (~700 Myr) open cluster Praesepe. The strong superflare on 2M0831+2042 has an equivalent duration (ED) of 13.7 hr and an estimated energy of 1.3 X 10^35 erg. We observed five superflares on 2M0837+2050, on which the strongest superflare has an ED of 46.4 hr and an estimated energy of 3.5 X 10^35 erg. This energy is larger by 2.7 orders of magnitude than the largest flare observed on the older (7.6 Gyr) planet-hosting M8 dwarf TRAPPIST-1. Furthermore, we also observed five superflares on 2M0831+2244 which is probably a field star. The estimated energy of the strongest superflare on 2M0831+2244 is 6.1 X 10^34 erg. 2M0831+2042, 2M0837+2050 and 2MASS J0831+2244 have rotation periods of 0.556\pm0.002, 0.193\pm0.000 and 0.292\pm0.001 d respectively, which are measured by using K2 light curves. We compare the flares of younger targets with those of TRAPPIST-1 and discuss the possible impacts of such flares on planets in the habitable zone of late-M dwarfs.
[39]  oai:arXiv.org:1811.06985  [pdf] - 1787063
Anisotropic winds in Wolf-Rayet binary identify potential gamma-ray burst progenitor
Comments: Accepted to Nature Astronomy on 2018 Oct 1. 27 pages, 9 figures, and 4 tables
Submitted: 2018-11-16
The massive evolved Wolf-Rayet stars sometimes occur in colliding-wind binary systems in which dust plumes are formed as a result of the collision of stellar winds. These structures are known to encode the parameters of the binary orbit and winds. Here, we report observations of a previously undiscovered Wolf-Rayet system, 2XMM J160050.7-514245, with a spectroscopically determined wind speed of $\approx$3400 km s$^{-1}$. In the thermal infrared, the system is adorned with a prominent $\approx$12$''$ spiral dust plume, revealed by proper motion studies to be expanding at only $\approx$570 km s$^{-1}$. As the dust and gas appear coeval, these observations are inconsistent with existing models of the dynamics of such colliding wind systems. We propose that this contradiction can be resolved if the system is capable of launching extremely anisotropic winds. Near-critical stellar rotation is known to drive such winds, suggesting this Wolf-Rayet system as a potential Galactic progenitor system for long-duration gamma-ray bursts.
[40]  oai:arXiv.org:1811.03616  [pdf] - 1842409
Meter- to Millimeter Emission from Cool Stellar Systems: Latest Results, Synergies Across the Spectrum, and Outlook for the Next Decade
Comments: Splinter session summary, to appear in the proceedings of the 20th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun (ed. S. J. Wolk)
Submitted: 2018-11-08
Radio observations of cool stellar systems provide unique information on their magnetic fields, high-energy processes, and chemistry. Buoyed by powerful new instruments (e.g. ALMA, JVLA, LOFAR), advances in related fields (e.g., the Gaia astrometric revolution), and above all a renewed interest in the relevant stellar astrophysics, stellar radio astronomy is experiencing a renaissance. In this splinter session, participants took stock of the present state of stellar radio astronomy to chart a course for the field's future.
[41]  oai:arXiv.org:1802.09525  [pdf] - 1754934
Associating Fast Radio Bursts with Extragalactic Radio Sources: General Methodology and a Search for a Counterpart to FRB 170107
Comments: Updated to match published ApJ version (minor modifications to text)
Submitted: 2018-02-26, last modified: 2018-09-23
The discovery of a repeating fast radio burst has led to the first precise localization, an association with a dwarf galaxy, and the identification of a coincident persistent radio source. However, further localizations are required to determine the nature of FRBs, the sources powering them, and the possibility of multiple populations. Here we investigate the use of associated persistent radio sources to establish FRB counterparts, taking into account the localization area and the persistent source flux density. Due to the lower areal number density of radio sources compared to faint optical sources, robust associations can be achieved for less precise localizations as compared to direct optical host galaxy associations. For generally larger localizations which preclude robust associations, the number of candidate hosts can be reduced based on the ratio of radio-to-optical brightness. We find that confident associations with $\sim 0.01-$1 mJy sources, comparable to the luminosity of the persistent source associated with FRB 121102 over the redshift range $z \approx 0.1 - 1$, require FRB localizations of $\lesssim 20''$. In the absence of a robust association, constraints can be placed on the luminosity of an associated radio source as a function of localization and DM. For DM $\approx 1000 \rm \ pc \ cm^{-3}$, an upper limit comparable to the luminosity of the FRB 121102 persistent source can be placed if the localization is $\lesssim 10''$. We apply our analysis to the case of the ASKAP FRB 170107, using optical and radio observations of the localization region. We identify two candidate hosts based on a ratio of radio-to-optical brightness of $\gtrsim 100$. We find that if one of these is associated with FRB 170107, the resulting radio luminosity ($1 \times 10^{29} - 4 \times 10^{30} \ \rm erg \ s^{-1} \ Hz^{-1}$) is comparable to the luminosity of the FRB 121102 persistent source.
[42]  oai:arXiv.org:1809.05113  [pdf] - 1771690
The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad-Line Region
Comments: 26 pages, 14 figures. Accepted for publication in ApJ
Submitted: 2018-09-13
We present models of the H$\beta$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ionizing source. The dynamics in our sample range from near-circular elliptical orbits to inflowing or outflowing trajectories. We measure black hole masses of $\log_{10}(M_{\rm BH}/M_\odot) = 6.48^{+0.21}_{-0.18}$ for PG 1310$-$108, $7.50^{+0.25}_{-0.18}$ for Mrk 50, $7.46^{+0.15}_{-0.21}$ for Mrk 141, $7.58^{+0.08}_{-0.08}$ for Mrk 279, $7.11^{+0.20}_{-0.17}$ for Mrk 1511, $6.65^{+0.27}_{-0.15}$ for NGC 4593, and $6.94^{+0.14}_{-0.14}$ for Zw 229$-$015. We use these black hole mass measurements along with cross-correlation time lags and line widths to recover the scale factor $f$ used in traditional reverberation mapping measurements. Combining our results with other studies that use this modeling technique, bringing our sample size to 16, we calculate a scale factor that can be used for measuring black hole masses in other reverberation mapping campaigns. When using the root-mean-square (rms) spectrum and using the line dispersion to measure the line width, we find $\log_{10}(f_{{\rm rms},\sigma})_{\rm pred} = 0.57 \pm 0.19$. Finally, we search for correlations between $f$ and other AGN and BLR parameters and find marginal evidence that $f$ is correlated with $M_{\rm BH}$ and the BLR inclination angle, but no significant evidence of a correlation with the AGN luminosity or Eddington ratio.
[43]  oai:arXiv.org:1808.04838  [pdf] - 1745972
The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2016 follow-up campaign. I. Overview and classification of candidates selected by two techniques
Comments: (MNRAS in press; color scheme of Figures 1-4 improved with respect to v1)
Submitted: 2018-08-14, last modified: 2018-09-05
The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the Outlier Selection Technique and a morphological algorithm, and presenting lens models of a system, which could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars (NIQs), which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range 6-35\%. This rate is comparable to that of previous searches like SQLS even though the parent dataset of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.
[44]  oai:arXiv.org:1808.09476  [pdf] - 1741322
First ALMA Light Curve Constrains Refreshed Reverse Shocks and Jet Magnetization in GRB 161219B
Comments: 23 pages, 17 figures, 9 tables. Published in ApJ
Submitted: 2018-08-28
We present detailed multi-wavelength observations of GRB 161219B at $z=0.1475$, spanning the radio to X-ray regimes, and the first ALMA light curve of a GRB afterglow. The cm- and mm-band observations before $8.5$ d require emission in excess of that produced by the afterglow forward shock (FS). These data are consistent with radiation from a refreshed reverse shock (RS) produced by the injection of energy into the FS, signatures of which are also present in the X-ray and optical light curves. We infer a constant-density circumburst environment with an extremely low density, $n_0\approx 3\times10^{-4}$ cm$^{-3}$ and show that this is a characteristic of all strong RS detections to date. The VLA observations exhibit unexpected rapid variability on $\sim$ minute timescales, indicative of strong interstellar scintillation. The X-ray, ALMA, and VLA observations together constrain the jet break time, $t_{\rm jet}\approx32$ day, yielding a wide jet opening angle of $\theta_{\rm jet}\approx13^{\circ}$, implying beaming corrected $\gamma$-ray and kinetic energies of $E_{\gamma}\approx4.9\times10^{48}$ erg and $E_{\rm K}\approx1.3\times10^{50}$ erg, respectively. Comparing the RS and FS emission, we show that the ejecta are only weakly magnetized, with relative magnetization, $R_{\rm B}\approx1$, compared to the FS. These direct, multi-frequency measurements of a refreshed RS spanning the optical to radio bands highlight the impact of radio and millimeter data in probing the production and nature of GRB jets.
[45]  oai:arXiv.org:1808.08597  [pdf] - 1739200
Design and Characterization of a Balloon-Borne Diffraction-Limited Submillimeter Telescope Platform for BLAST-TNG
Comments: 19 pages, 11 figures, 2018 SPIE Astronomical Telescopes & Instrumentation Conference
Submitted: 2018-08-26
The Next Generation Balloon-borne Large Aperture Submillimeter Telescope (BLAST-TNG) is a submillimeter mapping experiment planned for a 28 day long-duration balloon (LDB) flight from McMurdo Station, Antarctica during the 2018-2019 season. BLAST-TNG will detect submillimeter polarized interstellar dust emission, tracing magnetic fields in galactic molecular clouds. BLAST-TNG will be the first polarimeter with the sensitivity and resolution to probe the $\sim$0.1 parsec-scale features that are critical to understanding the origin of structures in the interstellar medium. With three detector arrays operating at 250, 350, and 500 $\mu$m (1200, 857, and 600 GHz), BLAST-TNG will obtain diffraction-limited resolution at each waveband of 30, 41, and 59 arcseconds respectively. To achieve the submillimeter resolution necessary for its science goals, the BLAST-TNG telescope features a 2.5 m aperture carbon fiber composite primary mirror, one of the largest mirrors flown on a balloon platform. Successful performance of such a large telescope on a balloon-borne platform requires stiff, lightweight optical components and mounting structures. Through a combination of optical metrology and finite element modeling of thermal and mechanical stresses on both the telescope optics and mounting structures, we expect diffraction-limited resolution at all our wavebands. We expect pointing errors due to deformation of the telescope mount to be negligible. We have developed a detailed thermal model of the sun shielding, gondola, and optical components to optimize our observing strategy and increase the stability of the telescope over the flight. We present preflight characterization of the telescope and its platform.
[46]  oai:arXiv.org:1808.08489  [pdf] - 1739192
Preflight Characterization of the BLAST-TNG Receiver and Detector Arrays
Comments: 15 pages, 7 figures, SPIE Astronomical Telescopes and Instrumentation Conference
Submitted: 2018-08-25
The Next Generation Balloon-borne Large Aperture Submillimeter Telescope (BLAST-TNG) is a submillimeter mapping experiment planned for a 28 day long-duration balloon (LDB) flight from McMurdo Station, Antarctica during the 2018-2019 season. BLAST-TNG will detect submillimeter polarized interstellar dust emission, tracing magnetic fields in galactic molecular clouds. BLAST-TNG will be the first polarimeter with the sensitivity and resolution to probe the $\sim$0.1 parsec-scale features that are critical to understanding the origin of structures in the interstellar medium. BLAST-TNG features three detector arrays operating at wavelengths of 250, 350, and 500 $\mu$m (1200, 857, and 600 GHz) comprised of 918, 469, and 272 dual-polarization pixels, respectively. Each pixel is made up of two crossed microwave kinetic inductance detectors (MKIDs). These arrays are cooled to 275 mK in a cryogenic receiver. Each MKID has a different resonant frequency, allowing hundreds of resonators to be read out on a single transmission line. This inherent ability to be frequency-domain multiplexed simplifies the cryogenic readout hardware, but requires careful optical testing to map out the physical location of each resonator on the focal plane. Receiver-level optical testing was carried out using both a cryogenic source mounted to a movable xy-stage with a shutter, and a beam-filling, heated blackbody source able to provide a 10-50 $^\circ$C temperature chop. The focal plane array noise properties, responsivity, polarization efficiency, instrumental polarization were measured. We present the preflight characterization of the BLAST-TNG cryogenic system and array-level optical testing of the MKID detector arrays in the flight receiver.
[47]  oai:arXiv.org:1805.02870  [pdf] - 1736377
A Decline in the X-ray through Radio Emission from GW170817 Continues to Support an Off-Axis Structured Jet
Comments: Accepted to ApJL. Updated version includes new VLA observations extending through 2018 June 2
Submitted: 2018-05-08, last modified: 2018-08-13
We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-290$ days post-merger, at radio (Karl G. Jansky Very Large Array; VLA), X-ray (Chandra X-ray Observatory) and optical (Hubble Space Telescope; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio-to-X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously published models of a successful structured jet expanding into a low-density circumbinary medium, but pure cocoon models with a choked jet cannot be ruled out. If future observations continue to track our predictions, we expect that the radio and X-ray emission will remain detectable until $\sim 1000$ days post-merger.
[48]  oai:arXiv.org:1807.03304  [pdf] - 1721210
The stellar rotation--activity relationship in fully convective M dwarfs
Comments: 12 pages, 5 figures, accepted for publication in MNRAS
Submitted: 2018-07-09
The coronal activity-rotation relationship is considered to be a proxy for the underlying stellar dynamo responsible for magnetic activity in solar and late-type stars. While this has been studied in considerable detail for partly-convective stars that are believed to operate an interface dynamo, it is poorly unconstrained in fully-convective stars that lack the necessary shear layer between radiative core and the convective envelope. We present new X-ray observations of 19 slowly-rotating fully-convective stars with rotation periods from the MEarth Project. We use these to calculate X-ray luminosities (or upper limits for undetected sources) and combine these with existing measurements from \citet{wrig16b}. We confirm the existence of fully-convective stars in the X-ray unsaturated regime and find that these objects follow the same rotation-activity relationship seen for partly-convective stars. We measure the power-law slope of the relationship between Rossby number (the ratio of the rotation period to the convective turnover time) and the fractional X-ray luminosity for X-ray unsaturated fully-convective stars for the first time, and find it to be consistent with that found for partly-convective stars. We discuss this implications of this result for our understanding of stellar magnetic dynamos in fully- and partly-convective stars. Finally, we also use this data to improve empirical estimates of the convective turnover time for fully-convective stars.
[49]  oai:arXiv.org:1806.06313  [pdf] - 1701057
Magnetic fields in forming stars with the ngVLA
Comments: 14 pages, 2 figures. An ngVLA Science Book chapter
Submitted: 2018-06-16
The magnetic field plays an important role in every stage of the star-formation process from the collapse of the initial protostellar core to the star's arrival on the main sequence. Consequently, the goal of this science case is to explore a wide range of magnetic phenomena that can be investigated using the polarization capabilities of the Next Generation Very Large Array (ngVLA). These include (1) magnetic fields in protostellar cores via polarized emission from aligned dust grains, including in regions optically thick at wavelengths observable by the Atacama Large Millimeter/submillimeter Array (ALMA); (2) magnetic fields in both protostellar cores and molecular outflows via spectral-line polarization from the Zeeman and Goldreich-Kylafis effects; (3) magnetic fields in protostellar jets via polarized synchrotron emission; and (4) gyrosynchrotron emission from magnetospheres around low-mass stars.
[50]  oai:arXiv.org:1805.12151  [pdf] - 1732700
The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2016 follow-up campaign. II. New quasar lenses from double component fitting
Comments: 13 pages, 8 figures. MNRAS submitted (05/24). Approved by the DES Collaboration Wide Review
Submitted: 2018-05-30
We report upon the follow up of 34 candidate lensed quasars found in the Dark Energy Survey using NTT-EFOSC, Magellan-IMACS, KECK-ESI and SOAR-SAMI. These candidates were selected by a combination of double component fitting, morphological assessment and color analysis. Most systems followed up are indeed composed of at least one quasar image and 13 with two or more quasar images: two lenses, four projected binaries and seven Nearly Identical Quasar Pairs (NIQs). The two systems confirmed as genuine gravitationally lensed quasars are one quadruple at $z_s=1.713$ and one double at $z_s=1.515$. Lens modeling of these two systems reveals that both systems require very little contribution from the environment to reproduce the image configuration. Nevertheless, small flux anomalies can be observed in one of the images of the quad. Further observations of 9 inconclusive systems (including 7 NIQs) will allow to confirm (or not) their gravitational lens nature.
[51]  oai:arXiv.org:1805.11185  [pdf] - 1721145
K2 Ultracool Dwarfs Survey. IV. Monster flares observed on the young brown dwarf CFHT-BD-Tau 4
Comments: 10 pages, 4 figures, Accepted for publication in ApJ
Submitted: 2018-05-28
We present photometric measurements of two superflares observed on a very young brown dwarf CFHT-BD-Tau 4, observed during Campaign 13 of the \textit{Kepler K2} mission. The stronger of the two superflares brightened by a factor of $\sim$48 relative to the quiescent photospheric level, with an increase in \textit{Kepler} magnitude $\Delta \tilde{K_{p}}$ = -4.20. It has an equivalent duration of $\sim$107 hour, a flare duration of 1.7 day, and an estimated total bolometric (ultraviolet/optical/infrared) energy up to 2.1 $\times$ 10$^{38}$ erg. The weaker of the two superflares is a complex (multipeaked) flare with an estimated total bolometric (UV/optical/IR) energy up to 4.7 $\times$ 10$^{36}$ erg. They are the strongest flares observed on any brown dwarf so far. The flare energies are strongly dependent on the value of visual extinction parameter $A_{V}$ used for extinction correction. If we apply a solar flare-model to interpret the two superflares, we find that the magnetic fields are required to be stronger by as much as an order of magnitude than previous reports of field measurements in CFHT-BD-Tau 4 by Reiners et al. (2009b). On the other hand, if we interpret our data in terms of accretion, we find that the requisite rate of accretion for the stronger superflare exceeds the rates which have been reported for other young brown dwarfs.
[52]  oai:arXiv.org:1805.08192  [pdf] - 1728995
Spitzer Space Telescope Infrared Observations of the Binary Neutron Star Merger GW170817
Comments: 6 pages, 2 figures, submitted to ApJL
Submitted: 2018-05-21
We present Spitzer Space Telescope 3.6 and 4.5 micron observations of the binary neutron star merger GW170817 at 43, 74, and 264 days post-merger. Using the final observation as a template, we uncover a source at the position of GW170817 at 4.5 micron with a brightness of 22.9+/-0.3 AB mag at 43 days and 23.8+/-0.3 AB mag at 74 days (the uncertainty is dominated by systematics from the image subtraction); no obvious source is detected at 3.6 micron to a 3-sigma limit of >23.3 AB mag in both epochs. The measured brightness is dimmer by a factor of about 2-3 times compared to our previously published kilonova model, which is based on UV, optical, and near-IR data at <30 days. However, the observed fading rate and color (m_{3.6}-m_{4.5}> 0 AB mag) are consistent with our model. We suggest that the discrepancy is likely due to a transition to the nebular phase, or a reduced thermalization efficiency at such late time. Using the Spitzer data as a guide, we briefly discuss the prospects of observing future binary neutron star mergers with Spitzer (in LIGO/Virgo Observing Run 3) and the James Webb Space Telescope (in LIGO/Virgo Observing Run 4 and beyond).
[53]  oai:arXiv.org:1805.00953  [pdf] - 1771536
HI 21cm Cosmology and the Bi-spectrum: Closure Diagnostics in Massively Redundant Interferometric Arrays
Comments: 32 pages. 11 figures. Accepted to Radio Science
Submitted: 2018-05-02
New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.
[54]  oai:arXiv.org:1804.09183  [pdf] - 1747897
COSMOGRAIL XVII: Time delays for the quadruply imaged quasar PG 1115+080
Comments: 15 pages, 6 figures, submitted to A&A. Comments welcome
Submitted: 2018-04-24
We present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Dt(AB) = 8.3+1.5-1.6 days (18.7% precision), Dt(AC) = 9.9+1.1-1.1 days (11.1%) and Dt(BC) = 18.8+1.6-1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope.
[55]  oai:arXiv.org:1803.07708  [pdf] - 1694059
K2 Ultracool Dwarfs Survey. III. White Light Flares are Ubiquitous in M6-L0 Dwarfs
Comments: 19 pages, 12 figures, Accepted for publication in ApJ, no any major changes from the previous version
Submitted: 2018-03-20, last modified: 2018-03-28
We report the white light flare rates for 10 ultracool dwarfs (UCDs) using \textit{Kepler K2} short cadence data. Among our sample stars, two have spectral type M6, three are M7, three are M8 and two are L0. Most of our targets are old low mass stars. We identify a total of 283 flares in all of the stars in our sample, with \textit{Kepler} energies in the range log \textit{E$_{Kp}$} $\sim$(29 - 33.5) erg. Using the maximum-likelihood method of line fitting, we find that the flare frequency distribution (FFD) for each star in our sample follows a power law with slope -$\alpha$ in range -(1.3-2.0). We find that cooler objects tend to have shallower slopes. For some of our targets, the FFD follows either a broken power law, or a power law with an exponential cutoff. For the L0 dwarf 2MASS J12321827-0951502, we find a very shallow slope (-$\alpha$ $=$ -1.3) in the \textit{Kepler} energy range (0.82-130)$\times$10$^{30}$ erg: this L0 dwarf has flare rates which are comparable to the rates of high energy flares in stars of earlier spectral types. In addition, we report photometry of two superflares: one on the L0 dwarf 2MASS J12321827-0951502 and another on the M7 dwarf 2MASS J08352366+1029318. In case of 2MASS J12321827-0951502, we report a flare brightening by a factor of $\sim$144 relative to the quiescent photospheric level. Likewise, for 2MASS J08352366+1029318, we report a flare brightening by a factor of $\sim$60 relative to the quiescent photospheric level. These two superflares have bolometric (UV/optical/infrared) energies 3.6 $\times$ 10$^{33}$ erg and 8.9 $\times$ 10$^{33}$ erg respectively, while the FWHM time scales are very short, $\sim$2 minutes. We find that the M8 star TRAPPIST-1 is more active than the M8.5 dwarf: 2M03264453+1919309, but less active than another M8 dwarf (2M12215066-0843197).
[56]  oai:arXiv.org:1803.06094  [pdf] - 1651022
Turbulence, Accretion Braking Torques and Efficient Jets Without Magnetocentrifugal Acceleration: Core Concepts
Comments: 37 pages, 9 figures, submitted to ApJ
Submitted: 2018-03-16
I discuss three mutually-supportive notions or assumptions regarding jets and accretion. The first is magnetocentrifugal acceleration (MCA), the overwhelmingly favored mechanism for the production of jets in most steady accreting systems. The second is the zero-torque inner boundary condition. The third is that effective viscous dissipation is like real dissipation, leading directly to heating. All three assumptions fit nicely together in a manner that is simple, persuasive, and mutually-consistent. All, I argue, are incorrect. For concreteness I focus on protostars. Magnetohydrodynamic (MHD) turbulence in accretion is not a sink of energy, but a reservoir, capable of doing mechanical work directly and therefore efficiently, rather than solely through ohmic ("viscous") heating. Advection of turbulence energy reduces the effective radiative efficiency, and may help solve the missing boundary-layer emission problem. The angular momentum problem, whereby accretion spins up a protostar to breakup, is resolved by allowing direct viscous coupling to the protostar, permitting substantially greater energy to be deposited into the accretion flow than otherwise possible. This goes not into heat, but into a turbulent, tangled, buoyant toroidal magnetic field. I argue that there is neither an angular momentum problem nor an efficiency problem that MCA is needed to solve. Moreover, the turbulent magnetic field has ample ability not just to collimate but to accelerate gas, first radially inwards through tension forces and then vertically through pressure forces, without any MCA mechanism. I suggest then that jets, particularly the most powerful and well-collimated protostellar jets, are not magnetocentrifugally driven.
[57]  oai:arXiv.org:1801.06080  [pdf] - 1648816
A Precise Distance to the Host Galaxy of the Binary Neutron Star Merger GW170817 Using Surface Brightness Fluctuations
Comments: ApJ Letters in press
Submitted: 2018-01-18, last modified: 2018-02-27
The joint detection of gravitational waves and electromagnetic radiation from the binary neutron star (BNS) merger GW170817 has provided unprecedented insight into a wide range of physical processes: heavy element synthesis via the $r$-process; the production of relativistic ejecta; the equation of state of neutron stars and the nature of the merger remnant; the binary coalescence timescale; and a measurement of the Hubble constant via the "standard siren" technique. In detail, all of these results depend on the distance to the host galaxy of the merger event, NGC4993. In this paper we measure the surface brightness fluctuation (SBF) distance to NGC4993 in the F110W and F160W passbands of the Wide Field Camera 3 Infrared Channel on the Hubble Space Telescope (HST). For the preferred F110W passband we derive a distance modulus of $m{-}M=33.05\pm0.08\pm0.10$ mag, or a linear distance $d=40.7\pm1.4\pm1.9$ Mpc (random and systematic errors, respectively); a virtually identical result is obtained from the F160W data. This is the most precise distance to NGC4993 available to date. Combining our distance measurement with the corrected recession velocity of NGC4993 implies a Hubble constant $H_0=71.9\pm 7.1~km~s^{-1}~Mpc^{-1}$. A comparison of our result to the GW-inferred value of $H_0$ indicates a binary orbital inclination of $i\,{\gtrsim}\,137~\deg$. The SBF technique can be applied to early-type host galaxies of BNS mergers to ${\sim\,}100$ Mpc with HST and possibly as far as ${\sim\,}300$ Mpc with the James Webb Space Telescope, thereby helping to break the inherent distance-inclination degeneracy of the GW signals at distances where many future BNS mergers are likely to be detected.
[58]  oai:arXiv.org:1801.03531  [pdf] - 1663278
The Binary Neutron Star event LIGO/VIRGO GW170817 a hundred and sixty days after merger: synchrotron emission across the electromagnetic spectrum
Comments: Updated with the latest VLA and Chandra data
Submitted: 2018-01-10, last modified: 2018-02-25
We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d after merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency $\nu_c$ is above the X-ray band and the synchrotron frequency $\nu_m$ is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index $p$ of the distribution of non-thermal relativistic electrons $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{not} constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, \emph{and} emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.
[59]  oai:arXiv.org:1710.02145  [pdf] - 1686642
MOSFiT: Modular Open-Source Fitter for Transients
Comments: 18 pages, 6 figures, accepted to ApJS. Revision cleans up descriptions and adds additional code examples. Documentation available at https://mosfit.readthedocs.io
Submitted: 2017-10-05, last modified: 2018-02-20
Much of the progress made in time-domain astronomy is accomplished by relating observational multi-wavelength time series data to models derived from our understanding of physical laws. This goal is typically accomplished by dividing the task in two: collecting data (observing), and constructing models to represent that data (theorizing). Owing to the natural tendency for specialization, a disconnect can develop between the best available theories and the best available data, potentially delaying advances in our understanding new classes of transients. We introduce MOSFiT: the Modular Open-Source Fitter for Transients, a Python-based package that downloads transient datasets from open online catalogs (e.g., the Open Supernova Catalog), generates Monte Carlo ensembles of semi-analytical light curve fits to those datasets and their associated Bayesian parameter posteriors, and optionally delivers the fitting results back to those same catalogs to make them available to the rest of the community. MOSFiT is designed to help bridge the gap between observations and theory in time-domain astronomy; in addition to making the application of existing models and creation of new models as simple as possible, MOSFiT yields statistically robust predictions for transient characteristics, with a standard output format that includes all the setup information necessary to reproduce a given result. As large-scale surveys such as LSST discover entirely new classes of transients, tools such as MOSFiT will be critical for enabling rapid comparison of models against data in statistically consistent, reproducible, and scientifically beneficial ways.
[60]  oai:arXiv.org:1710.11576  [pdf] - 1605053
The Combined Ultraviolet, Optical, and Near-Infrared Light Curves of the Kilonova Associated with the Binary Neutron Star Merger GW170817: Unified Data Set, Analytic Models, and Physical Implications
Comments: 35 pages; 5 figures; 3 tables. Significant changes made to Table 3. Accepted to ApJL
Submitted: 2017-10-31, last modified: 2017-11-28
We present the first effort to aggregate, homogenize, and uniformly model the combined ultraviolet, optical, and near-infrared dataset for the electromagnetic counterpart of the binary neutron star merger GW170817. By assembling all of the available data from 18 different papers and 46 different instruments, we are able to identify and mitigate systematic offsets between individual datasets, and to identify clear outlying measurements, with the resulting pruned and adjusted dataset offering an opportunity to expand the study of the kilonova. The unified dataset includes 647 individual flux measurements, spanning 0.45 to 29.4 days post-merger, and thus has greater constraining power for physical models than any single dataset. We test a number of semi-analytical models and find that the data are well modeled with a three-component kilonova model: a "blue" lanthanide-poor component with Mej~0.020 Msol and vej~0.27c; an intermediate opacity "purple" component with Mej~0.047 Msol and vej~0.15c; and a "red" lanthanide-rich component with Mej~0.011 Msol and vej~0.14c. We further explore the possibility of ejecta asymmetry and its impact on the estimated parameters. From the inferred parameters we draw conclusions about the physical mechanisms responsible for the various ejecta components, the properties of the neutron stars, and, combined with an up-to-date merger rate, the implications for r-process enrichment via this channel. To facilitate future studies of this keystone event we make the unified dataset and our modeling code public.
[61]  oai:arXiv.org:1710.06426  [pdf] - 1608518
Improved constraints on H0 from a combined analysis of gravitational-wave and electromagnetic emission from GW170817
Comments: Submitted to ApJL. Comments Welcome. Revised uncertainties in v2
Submitted: 2017-10-17, last modified: 2017-10-19
The luminosity distance measurement of GW170817 derived from GW analysis in Abbott et al. 2017 (here, A17:H0) is highly correlated with the measured inclination of the NS-NS system. To improve the precision of the distance measurement, we attempt to constrain the inclination by modeling the broad-band X-ray-to-radio emission from GW170817, which is dominated by the interaction of the jet with the environment. We update our previous analysis and we consider the radio and X-ray data obtained at $t<40$ days since merger. We find that the afterglow emission from GW170817 is consistent with an off-axis relativistic jet with energy $10^{48}\,\rm{erg}<E_{k}\le 3\times 10^{50} \,\rm{erg}$ propagating into an environment with density $n\sim10^{-2}-10^{-4} \,\rm{cm^{-3}}$, with preference for wider jets (opening angle $\theta_j=15$ deg). For these jets, our modeling indicates an off-axis angle $\theta_{\rm obs}\sim25-50$ deg. We combine our constraints on $\theta_{\rm obs}$ with the joint distance-inclination constraint from LIGO. Using the same $\sim 170$ km/sec peculiar velocity uncertainty assumed in A17:H0 but with an inclination constraint from the afterglow data, we get a value of $H_0=$$74.0 \pm \frac{11.5}{7.5}$ $\mbox{km/s/Mpc}$, which is higher than the value of $H_0=$$70.0 \pm \frac{12.0}{8.0}$ $\mbox{km/s/Mpc}$ found in A17:H0. Further, using a more realistic peculiar velocity uncertainty of 250 km/sec derived from previous work, we find $H_0=$$75.5 \pm \frac{11.6}{9.6}$ km/s/Mpc for H0 from this system. We note that this is in modestly better agreement with the local distance ladder than the Planck CMB, though a significant such discrimination will require $\sim 50$ such events. Future measurements at $t>100$ days of the X-ray and radio emission will lead to tighter constraints.
[62]  oai:arXiv.org:1710.05454  [pdf] - 1589770
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. IV. Detection of Near-infrared Signatures of r-process Nucleosynthesis with Gemini-South
Comments: 9 pages, 4 figures, 1 table, ApJL, in press
Submitted: 2017-10-16
We present a near-infrared spectral sequence of the electromagnetic counterpart to the binary neutron star merger GW170817 detected by Advanced LIGO/Virgo. Our dataset comprises seven epochs of J+H spectra taken with FLAMINGOS-2 on Gemini-South between 1.5 and 10.5 days after the merger. In the initial epoch, the spectrum is dominated by a smooth blue continuum due to a high-velocity, lanthanide-poor blue kilonova component. Starting the following night, all subsequent spectra instead show features that are similar to those predicted in model spectra of material with a high concentration of lanthanides, including spectral peaks near 1.07 and 1.55 microns. Our fiducial model with 0.04 M_sun of ejecta, an ejection velocity of v=0.1c, and a lanthanide concentration of X_lan=1e-2 provides a good match to the spectra taken in the first five days, although it over-predicts the late-time fluxes. We also explore models with multiple fitting components, in each case finding that a significant abundance of lanthanide elements is necessary to match the broad spectral peaks that we observe starting at 2.5 d after the merger. These data provide direct evidence that binary neutron star mergers are significant production sites of even the heaviest r-process elements.
[63]  oai:arXiv.org:1710.05431  [pdf] - 1589749
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. V. Rising X-ray Emission from an Off-Axis Jet
Comments: 7 Pages, 4 Figures, ApJL, In Press. Keywords: GW170817, LVC
Submitted: 2017-10-16
We report the discovery of rising X-ray emission from the binary neutron star (BNS) merger event GW170817. This is the first detection of X-ray emission from a gravitational-wave source. Observations acquired with the Chandra X-ray Observatory (CXO) at t~2.3 days post merger reveal no significant emission, with L_x<=3.2x10^38 erg/s (isotropic-equivalent). Continued monitoring revealed the presence of an X-ray source that brightened with time, reaching L_x\sim 9x10^39 erg/s at ~15.1 days post merger. We interpret these findings in the context of isotropic and collimated relativistic outflows (both on- and off-axis). We find that the broad-band X-ray to radio observations are consistent with emission from a relativistic jet with kinetic energy E_k~10^49-10^50 erg, viewed off-axis with theta_obs~ 20-40 deg. Our models favor a circumbinary density n~ 0.0001-0.01 cm-3, depending on the value of the microphysical parameter epsilon_B=10^{-4}-10^{-2}. A central-engine origin of the X-ray emission is unlikely. Future X-ray observations at $t\gtrsim 100$ days, when the target will be observable again with the CXO, will provide additional constraints to solve the model degeneracies and test our predictions. Our inferences on theta_obs are testable with gravitational wave information on GW170817 from Advanced LIGO/Virgo on the binary inclination.
[64]  oai:arXiv.org:1710.05438  [pdf] - 1589756
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VIII. A Comparison to Cosmological Short-duration Gamma-ray Bursts
Comments: 11 Pages, 4 Figures, 1 Table. ApJL, In Press. Keywords: GW170817, LVC
Submitted: 2017-10-16
We present a comprehensive comparison of the properties of the radio through X-ray counterpart of GW170817 and the properties of short-duration gamma-ray bursts (GRBs). For this effort, we utilize a sample of 36 short GRBs spanning a redshift range of $z \approx 0.12-2.6$ discovered over 2004-2017. We find that the counterpart to GW170817 has an isotropic-equivalent luminosity that is $\approx 3000$ times less than the median value of on-axis short GRB X-ray afterglows, and $\gtrsim10^{4}$ times less than that for detected short GRB radio afterglows. Moreover, the allowed jet energies and particle densities inferred from the radio and X-ray counterparts to GW170817 and on-axis short GRB afterglows are remarkably similar, suggesting that viewing angle effects are the dominant, and perhaps only, difference in their observed radio and X-ray behavior. From comparison to previous claimed kilonovae following short GRBs, we find that the optical and near-IR counterpart to GW170817 is comparatively under-luminous by a factor of $\approx 3-5$, indicating a range of kilonova luminosities and timescales. A comparison of the optical limits following short GRBs on $\lesssim 1$ day timescales also rules out a "blue" kilonova of comparable optical isotropic-equivalent luminosity in one previous short GRB. Finally, we investigate the host galaxy of GW170817, NGC4993, in the context of short GRB host galaxy stellar population properties. We find that NGC4993 is superlative in terms of its large luminosity, old stellar population age, and low star formation rate compared to previous short GRB hosts. Additional events within the Advanced LIGO/VIRGO volume will be crucial in delineating the properties of the host galaxies of NS-NS mergers, and connecting them to their cosmological counterparts.
[65]  oai:arXiv.org:1710.05456  [pdf] - 1589772
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. III. Optical and UV Spectra of a Blue Kilonova From Fast Polar Ejecta
Comments: ApJL, in press (GW170817, LVC)
Submitted: 2017-10-16
We present optical and ultraviolet spectra of the first electromagnetic counterpart to a gravitational wave (GW) source, the binary neutron star merger GW170817. Spectra were obtained nightly between 1.5 and 9.5 days post-merger, using the SOAR and Magellan telescopes; the UV spectrum was obtained with the \textit{Hubble Space Telescope} at 5.5 days. Our data reveal a rapidly-fading blue component ($T\approx5500$ K at 1.5 days) that quickly reddens; spectra later than $\gtrsim 4.5$ days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at $\sim 7900$ \AA\ at $t\lesssim 4.5$ days. The colours, rapid evolution and featureless spectrum are consistent with a "blue" kilonova from polar ejecta comprised mainly of light $r$-process nuclei with atomic mass number $A\lesssim 140$. This indicates a sight-line within $\theta_{\rm obs}\lesssim 45^{\circ}$ of the orbital axis. Comparison to models suggests $\sim0.03$ M$_\odot$ of blue ejecta, with a velocity of $\sim 0.3c$. The required lanthanide fraction is $\sim 10^{-4}$, but this drops to $<10^{-5}$ in the outermost ejecta. The large velocities point to a dynamical origin, rather than a disk wind, for this blue component, suggesting that both binary constituents are neutron stars (as opposed to a binary consisting of a neutron star and a black hole). For dynamical ejecta, the high mass favors a small neutron star radius of $\lesssim 12$ km. This mass also supports the idea that neutron star mergers are a major contributor to $r$-process nucleosynthesis.
[66]  oai:arXiv.org:1710.05457  [pdf] - 1589773
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VI. Radio Constraints on a Relativistic Jet and Predictions for Late-Time Emission from the Kilonova Ejecta
Comments: 8 pages, 4 figures, 1 table. ApJL, in press. Keywords: GW170817, LVC
Submitted: 2017-10-16
We present Very Large Array (VLA) and Atacama Large Millimeter/sub-millimeter Array ALMA radio observations of GW\,170817, the first Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo gravitational wave (GW) event from a binary neutron star merger and the first GW event with an electromagnetic (EM) counterpart. Our data include the first observations following the discovery of the optical transient at both the centimeter ($13.7$ hours post merger) and millimeter ($2.41$ days post merger) bands. We detect faint emission at 6 GHz at 19.47 and 39.23 days after the merger, but not in an earlier observation at 2.46 d. We do not detect cm/mm emission at the position of the optical counterpart at frequencies of 10-97.5 GHz at times ranging from 0.6 to 30 days post merger, ruling out an on-axis short gamma-ray burst (SGRB) for energies $\gtrsim 10^{48}$ erg. For fiducial SGRB parameters, our limits require an observer viewer angle of $\gtrsim 20^{\circ}$. The radio and X-ray data can be jointly explained as the afterglow emission from an SGRB with a jet energy of $\sim 10^{49}-10^{50}$ erg that exploded in a uniform density environment with $n\sim 10^{-4}-10^{-2}$ cm$^{-3}$, viewed at an angle of $\sim 20^{\circ}-40^{\circ}$ from the jet axis. Using the results of our light curve and spectral modeling, in conjunction with the inference of the circumbinary density, we predict the emergence of late-time radio emission from the deceleration of the kilonova (KN) ejecta on a timescale of $\sim 5-10$ years that will remain detectable for decades with next-generation radio facilities, making GW\,170817 a compelling target for long-term radio monitoring.
[67]  oai:arXiv.org:1710.05458  [pdf] - 1589774
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VII. Properties of the Host Galaxy and Constraints on the Merger Timescale
Comments: 9 Pages, 3 Figures, 2 Tables, ApJL, In Press. Keywords: GW170817, LVC
Submitted: 2017-10-16
We present the properties of NGC 4993, the host galaxy of GW170817, the first gravitational wave (GW) event from the merger of a binary neutron star (BNS) system and the first with an electromagnetic (EM) counterpart. We use both archival photometry and new optical/near-IR imaging and spectroscopy, together with stellar population synthesis models to infer the global properties of the host galaxy. We infer a star formation history peaked at $\gtrsim 10$ Gyr ago, with subsequent exponential decline leading to a low current star formation rate of 0.01 M$_{\odot}$ yr$^{-1}$, which we convert into a binary merger timescale probability distribution. We find a median merger timescale of $11.2^{+0.7}_{-1.4}$ Gyr, with a 90% confidence range of $6.8-13.6$ Gyr. This in turn indicates an initial binary separation of $\approx 4.5$ R$_{\odot}$, comparable to the inferred values for Galactic BNS systems. We also use new and archival $Hubble$ $Space$ $Telescope$ images to measure a projected offset of the optical counterpart of $2.1$ kpc (0.64$r_{e}$) from the center of NGC 4993 and to place a limit of $M_{r} \gtrsim -7.2$ mag on any pre-existing emission, which rules out the brighter half of the globular cluster luminosity function. Finally, the age and offset of the system indicates it experienced a modest natal kick with an upper limit of $\sim 200$ km s$^{-1}$. Future GW$-$EM observations of BNS mergers will enable measurement of their population delay time distribution, which will directly inform their viability as the dominant source of $r$-process enrichment in the Universe.
[68]  oai:arXiv.org:1710.05459  [pdf] - 1589775
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. I. Dark Energy Camera Discovery of the Optical Counterpart
Soares-Santos, M.; Holz, D. E.; Annis, J.; Chornock, R.; Herner, K.; Berger, E.; Brout, D.; Chen, H.; Kessler, R.; Sako, M.; Allam, S.; Tucker, D. L.; Butler, R. E.; Palmese, A.; Doctor, Z.; Diehl, H. T.; Frieman, J.; Yanny, B.; Lin, H.; Scolnic, D.; Cowperthwaite, P.; Neilsen, E.; Marriner, J.; Kuropatkin, N.; Hartley, W. G.; Paz-Chinchón, F.; Alexander, K. D.; Balbinot, E.; Blanchard, P.; Brown, D. A.; Carlin, J. L.; Conselice, C.; Cook, E. R.; Drlica-Wagner, A.; Drout, M. R.; Durret, F.; Eftekhari, T.; Farr, B.; Finley, D. A.; Foley, R. J.; Fong, W.; Fryer, C. L.; García-Bellido, J.; Gill, M. S . S.; Gruendl, R. A.; Hanna, C.; Kasen, D.; Li, T. S.; Lopes, P. A. A.; Lourenço, A. C. C.; Margutti, R.; Marshall, J. L.; Matheson, T.; Medina, G. E.; Metzger, B. D.; Muñoz, R. R.; Muir, J.; Nicholl, M.; Quataert, E.; Rest, A.; Sauseda, M.; Schlegel, D. J.; Secco, L. F.; Sobreira, F.; Stebbins, A.; Villar, V. A.; Walker, A. R.; Wester, W.; Williams, P. K. G.; Zenteno, A.; Zhang, Y.; Abbott, T. M. C.; Abdalla, F. B.; Banerji, M.; Bechtol, K.; Benoit-Lévy, A.; Bertin, E.; Brooks, D.; Buckley-Geer, E.; Burke, D. L.; Rosell, A. Carnero; Kind, M. Carrasco; Carretero, J.; Castander, F. J.; Crocce, M.; Cunha, C. E.; D'Andrea, C. B.; da Costa, L. N.; Davis, C.; Desai, S.; Dietrich, J. P.; Doel, P.; Eifler, T. F.; Fernandez, E.; Flaugher, B.; Fosalba, P.; Gaztanaga, E.; Gerdes, D. W.; Giannantonio, T.; Goldstein, D. A.; Gruen, D.; Gschwend, J.; Gutierrez, G.; Honscheid, K.; Jain, B.; James, D. J.; Jeltema, T.; Johnson, M. W. G.; Johnson, M. D.; Kent, S.; Krause, E.; Kron, R.; Kuehn, K.; Kuhlmann, S.; Lahav, O.; Lima, M.; Maia, M. A. G.; March, M.; McMahon, R. G.; Menanteau, F.; Miquel, R.; Mohr, J. J.; Nichol, R. C.; Nord, B.; Ogando, R. L. C.; Petravick, D.; Plazas, A. A.; Romer, A. K.; Roodman, A.; Rykoff, E. S.; Sanchez, E.; Scarpine, V.; Schubnell, M.; Sevilla-Noarbe, I.; Smith, M.; Smith, R. C.; Suchyta, E.; Swanson, M. E. C.; Tarle, G.; Thomas, D.; Thomas, R. C.; Troxel, M. A.; Vikram, V.; Wechsler, R. H.; Weller, J.
Comments: 9 Pages, 3 Figures, 2 Tables. ApJL, In Press. Keywords: GW170817, LVC
Submitted: 2017-10-16
We present the Dark Energy Camera (DECam) discovery of the optical counterpart of the first binary neutron star merger detected through gravitational wave emission, GW170817. Our observations commenced 10.5 hours post-merger, as soon as the localization region became accessible from Chile. We imaged 70 deg$^2$ in the $i$ and $z$ bands, covering 93\% of the initial integrated localization probability, to a depth necessary to identify likely optical counterparts (e.g., a kilonova). At 11.4 hours post-merger we detected a bright optical transient located $10.6''$ from the nucleus of NGC\,4993 at redshift $z=0.0098$, consistent (for $H_0 = 70$\, km s$^{-1}$ Mpc$^{-1}$) with the distance of $40 \pm 8$\, Mpc reported by the LIGO Scientific Collaboration and the Virgo Collaboration (LVC). At detection the transient had magnitudes $i\approx 17.30$ and $z\approx 17.45$, and thus an absolute magnitude of $M_i = -15.7$, in the luminosity range expected for a kilonova. We identified 1,500 potential transient candidates. Applying simple selection criteria aimed at rejecting background events such as supernovae, we find the transient associated with NGC\,4993 as the only remaining plausible counterpart, and reject chance coincidence at the 99.5\% confidence level. We therefore conclude that the optical counterpart we have identified near NGC\,4993 is associated with GW170817. This discovery ushers in the era of multi-messenger astronomy with gravitational waves, and demonstrates the power of DECam to identify the optical counterparts of gravitational-wave sources.
[69]  oai:arXiv.org:1710.05835  [pdf] - 1589824
A gravitational-wave standard siren measurement of the Hubble constant
Abbott, B. P.; Abbott, R.; Abbott, T. D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R. X.; Adya, V. B.; Affeldt, C.; Afrough, M.; Agarwal, B.; Agathos, M.; Agatsuma, K.; Aggarwal, N.; Aguiar, O. D.; Aiello, L.; Ain, A.; Ajith, P.; Allen, B.; Allen, G.; Allocca, A.; Altin, P. A.; Amato, A.; Ananyeva, A.; Anderson, S. B.; Anderson, W. G.; Angelova, S. V.; Antier, S.; Appert, S.; Arai, K.; Araya, M. C.; Areeda, J. S.; Arnaud, N.; Arun, K. G.; Ascenzi, S.; Ashton, G.; Ast, M.; Aston, S. M.; Astone, P.; Atallah, D. V.; Aufmuth, P.; Aulbert, C.; AultONeal, K.; Austin, C.; Avila-Alvarez, A.; Babak, S.; Bacon, P.; Bader, M. K. M.; Bae, S.; Baker, P. T.; Baldaccini, F.; Ballardin, G.; Ballmer, S. W.; Banagiri, S.; Barayoga, J. C.; Barclay, S. E.; Barish, B. C.; Barker, D.; Barkett, K.; Barone, F.; Barr, B.; Barsotti, L.; Barsuglia, M.; Barta, D.; Bartlett, J.; Bartos, I.; Bassiri, R.; Basti, A.; Batch, J. C.; Bawaj, M.; Bayley, J. C.; Bazzan, M.; Bécsy, B.; Beer, C.; Bejger, M.; Belahcene, I.; Bell, A. S.; Berger, B. K.; Bergmann, G.; Bero, J. J.; Berry, C. P. L.; Bersanetti, D.; Bertolini, A.; Betzwieser, J.; Bhagwat, S.; Bhandare, R.; Bilenko, I. A.; Billingsley, G.; Billman, C. R.; Birch, J.; Birney, R.; Birnholtz, O.; Biscans, S.; Biscoveanu, S.; Bisht, A.; Bitossi, M.; Biwer, C.; Bizouard, M. A.; Blackburn, J. K.; Blackman, J.; Blair, C. D.; Blair, D. G.; Blair, R. M.; Bloemen, S.; Bock, O.; Bode, N.; Boer, M.; Bogaert, G.; Bohe, A.; Bondu, F.; Bonilla, E.; Bonnand, R.; Boom, B. A.; Bork, R.; Boschi, V.; Bose, S.; Bossie, K.; Bouffanais, Y.; Bozzi, A.; Bradaschia, C.; Brady, P. R.; Branchesi, M.; Brau, J. E.; Briant, T.; Brillet, A.; Brinkmann, M.; Brisson, V.; Brockill, P.; Broida, J. E.; Brooks, A. F.; Brown, D. A.; Brown, D. D.; Brunett, S.; Buchanan, C. C.; Buikema, A.; Bulik, T.; Bulten, H. J.; Buonanno, A.; Buskulic, D.; Buy, C.; Byer, R. L.; Cabero, M.; Cadonati, L.; Cagnoli, G.; Cahillane, C.; Bustillo, J. Calderón; Callister, T. A.; Calloni, E.; Camp, J. B.; Canepa, M.; Canizares, P.; Cannon, K. C.; Cao, H.; Cao, J.; Capano, C. D.; Capocasa, E.; Carbognani, F.; Caride, S.; Carney, M. F.; Diaz, J. Casanueva; Casentini, C.; Caudill, S.; Cavaglià, M.; Cavalier, F.; Cavalieri, R.; Cella, G.; Cepeda, C. B.; Cerdá-Durán, P.; Cerretani, G.; Cesarini, E.; Chamberlin, S. J.; Chan, M.; Chao, S.; Charlton, P.; Chase, E.; Chassande-Mottin, E.; Chatterjee, D.; Chatziioannou, K.; Cheeseboro, B. D.; Chen, H. Y.; Chen, X.; Chen, Y.; Cheng, H. -P.; Chia, H.; Chincarini, A.; Chiummo, A.; Chmiel, T.; Cho, H. S.; Cho, M.; Chow, J. H.; Christensen, N.; Chu, Q.; Chua, A. J. K.; Chua, S.; Chung, A. K. W.; Chung, S.; Ciani, G.; Ciolfi, R.; Cirelli, C. E.; Cirone, A.; Clara, F.; Clark, J. A.; Clearwater, P.; Cleva, F.; Cocchieri, C.; Coccia, E.; Cohadon, P. -F.; Cohen, D.; Colla, A.; Collette, C. G.; Cominsky, L. R.; Constancio, M.; Conti, L.; Cooper, S. J.; Corban, P.; Corbitt, T. R.; Cordero-Carrión, I.; Corley, K. R.; Cornish, N.; Corsi, A.; Cortese, S.; Costa, C. A.; Coughlin, M. W.; Coughlin, S. B.; Coulon, J. -P.; Countryman, S. T.; Couvares, P.; Covas, P. B.; Cowan, E. E.; Coward, D. M.; Cowart, M. J.; Coyne, D. C.; Coyne, R.; Creighton, J. D. E.; Creighton, T. D.; Cripe, J.; Crowder, S. G.; Cullen, T. J.; Cumming, A.; Cunningham, L.; Cuoco, E.; Canton, T. Dal; Dálya, G.; Danilishin, S. L.; D'Antonio, S.; Danzmann, K.; Dasgupta, A.; Costa, C. F. Da Silva; Datrier, L. E. H.; Dattilo, V.; Dave, I.; Davier, M.; Davis, D.; Daw, E. J.; Day, B.; De, S.; DeBra, D.; Degallaix, J.; De Laurentis, M.; Deléglise, S.; Del Pozzo, W.; Demos, N.; Denker, T.; Dent, T.; De Pietri, R.; Dergachev, V.; De Rosa, R.; DeRosa, R. T.; De Rossi, C.; DeSalvo, R.; de Varona, O.; Devenson, J.; Dhurandhar, S.; Díaz, M. C.; Di Fiore, L.; Di Giovanni, M.; Di Girolamo, T.; Di Lieto, A.; Di Pace, S.; Di Palma, I.; Di Renzo, F.; Doctor, Z.; Dolique, V.; Donovan, F.; Dooley, K. L.; Doravari, S.; Dorrington, I.; Douglas, R.; Álvarez, M. Dovale; Downes, T. P.; Drago, M.; Dreissigacker, C.; Driggers, J. C.; Du, Z.; Ducrot, M.; Dupej, P.; Dwyer, S. E.; Edo, T. B.; Edwards, M. C.; Effler, A.; Eggenstein, H. -B.; Ehrens, P.; Eichholz, J.; Eikenberry, S. S.; Eisenstein, R. A.; Essick, R. C.; Estevez, D.; Etienne, Z. B.; Etzel, T.; Evans, M.; Evans, T. M.; Factourovich, M.; Fafone, V.; Fair, H.; Fairhurst, S.; Fan, X.; Farinon, S.; Farr, B.; Farr, W. M.; Fauchon-Jones, E. J.; Favata, M.; Fays, M.; Fee, C.; Fehrmann, H.; Feicht, J.; Fejer, M. M.; Fernandez-Galiana, A.; Ferrante, I.; Ferreira, E. C.; Ferrini, F.; Fidecaro, F.; Finstad, D.; Fiori, I.; Fiorucci, D.; Fishbach, M.; Fisher, R. P.; Fitz-Axen, M.; Flaminio, R.; Fletcher, M.; Fong, H.; Font, J. A.; Forsyth, P. W. F.; Forsyth, S. S.; Fournier, J. -D.; Frasca, S.; Frasconi, F.; Frei, Z.; Freise, A.; Frey, R.; Frey, V.; Fries, E. M.; Fritschel, P.; Frolov, V. V.; Fulda, P.; Fyffe, M.; Gabbard, H.; Gadre, B. U.; Gaebel, S. M.; Gair, J. R.; Gammaitoni, L.; Ganija, M. R.; Gaonkar, S. G.; Garcia-Quiros, C.; Garufi, F.; Gateley, B.; Gaudio, S.; Gaur, G.; Gayathri, V.; Gehrels, N.; Gemme, G.; Genin, E.; Gennai, A.; George, D.; George, J.; Gergely, L.; Germain, V.; Ghonge, S.; Ghosh, Abhirup; Ghosh, Archisman; Ghosh, S.; Giaime, J. A.; Giardina, K. D.; Giazotto, A.; Gill, K.; Glover, L.; Goetz, E.; Goetz, R.; Gomes, S.; Goncharov, B.; González, G.; Castro, J. M. Gonzalez; Gopakumar, A.; Gorodetsky, M. L.; Gossan, S. E.; Gosselin, M.; Gouaty, R.; Grado, A.; Graef, C.; Granata, M.; Grant, A.; Gras, S.; Gray, C.; Greco, G.; Green, A. C.; Gretarsson, E. M.; Groot, P.; Grote, H.; Grunewald, S.; Gruning, P.; Guidi, G. M.; Guo, X.; Gupta, A.; Gupta, M. K.; Gushwa, K. E.; Gustafson, E. K.; Gustafson, R.; Halim, O.; Hall, B. R.; Hall, E. D.; Hamilton, E. Z.; Hammond, G.; Haney, M.; Hanke, M. M.; Hanks, J.; Hanna, C.; Hannam, M. D.; Hannuksela, O. A.; Hanson, J.; Hardwick, T.; Harms, J.; Harry, G. M.; Harry, I. W.; Hart, M. J.; Haster, C. -J.; Haughian, K.; Healy, J.; Heidmann, A.; Heintze, M. C.; Heitmann, H.; Hello, P.; Hemming, G.; Hendry, M.; Heng, I. S.; Hennig, J.; Heptonstall, A. W.; Heurs, M.; Hild, S.; Hinderer, T.; Hoak, D.; Hofman, D.; Holt, K.; Holz, D. E.; Hopkins, P.; Horst, C.; Hough, J.; Houston, E. A.; Howell, E. J.; Hreibi, A.; Hu, Y. M.; Huerta, E. A.; Huet, D.; Hughey, B.; Husa, S.; Huttner, S. H.; Huynh-Dinh, T.; Indik, N.; Inta, R.; Intini, G.; Isa, H. N.; Isac, J. -M.; Isi, M.; Iyer, B. R.; Izumi, K.; Jacqmin, T.; Jani, K.; Jaranowski, P.; Jawahar, S.; Jiménez-Forteza, F.; Johnson, W. W.; Jones, D. I.; Jones, R.; Jonker, R. J. G.; Ju, L.; Junker, J.; Kalaghatgi, C. V.; Kalogera, V.; Kamai, B.; Kandhasamy, S.; Kang, G.; Kanner, J. B.; Kapadia, S. J.; Karki, S.; Karvinen, K. S.; Kasprzack, M.; Katolik, M.; Katsavounidis, E.; Katzman, W.; Kaufer, S.; Kawabe, K.; Kéfélian, F.; Keitel, D.; Kemball, A. J.; Kennedy, R.; Kent, C.; Key, J. S.; Khalili, F. Y.; Khan, I.; Khan, S.; Khan, Z.; Khazanov, E. A.; Kijbunchoo, N.; Kim, Chunglee; Kim, J. C.; Kim, K.; Kim, W.; Kim, W. S.; Kim, Y. -M.; Kimbrell, S. J.; King, E. J.; King, P. J.; Kinley-Hanlon, M.; Kirchhoff, R.; Kissel, J. S.; Kleybolte, L.; Klimenko, S.; Knowles, T. D.; Koch, P.; Koehlenbeck, S. M.; Koley, S.; Kondrashov, V.; Kontos, A.; Korobko, M.; Korth, W. Z.; Kowalska, I.; Kozak, D. B.; Krämer, C.; Kringel, V.; Krishnan, B.; Królak, A.; Kuehn, G.; Kumar, P.; Kumar, R.; Kumar, S.; Kuo, L.; Kutynia, A.; Kwang, S.; Lackey, B. D.; Lai, K. H.; Landry, M.; Lang, R. N.; Lange, J.; Lantz, B.; Lanza, R. K.; Lartaux-Vollard, A.; Lasky, P. D.; Laxen, M.; Lazzarini, A.; Lazzaro, C.; Leaci, P.; Leavey, S.; Lee, C. H.; Lee, H. K.; Lee, H. M.; Lee, H. W.; Lee, K.; Lehmann, J.; Lenon, A.; Leonardi, M.; Leroy, N.; Letendre, N.; Levin, Y.; Li, T. G. F.; Linker, S. D.; Littenberg, T. B.; Liu, J.; Liu, X.; Lo, R. K. L.; Lockerbie, N. A.; London, L. T.; Lord, J. E.; Lorenzini, M.; Loriette, V.; Lormand, M.; Losurdo, G.; Lough, J. D.; Lousto, C. O.; Lovelace, G.; Lück, H.; Lumaca, D.; Lundgren, A. P.; Lynch, R.; Ma, Y.; Macas, R.; Macfoy, S.; Machenschalk, B.; MacInnis, M.; Macleod, D. M.; Hernandez, I. Magaña; Magaña-Sandoval, F.; Zertuche, L. Magaña; Magee, R. M.; Majorana, E.; Maksimovic, I.; Man, N.; Mandic, V.; Mangano, V.; Mansell, G. L.; Manske, M.; Mantovani, M.; Marchesoni, F.; Marion, F.; Márka, S.; Márka, Z.; Markakis, C.; Markosyan, A. S.; Markowitz, A.; Maros, E.; Marquina, A.; Martelli, F.; Martellini, L.; Martin, I. W.; Martin, R. M.; Martynov, D. V.; Mason, K.; Massera, E.; Masserot, A.; Massinger, T. J.; Masso-Reid, M.; Mastrogiovanni, S.; Matas, A.; Matichard, F.; Matone, L.; Mavalvala, N.; Mazumder, N.; McCarthy, R.; McClelland, D. E.; McCormick, S.; McCuller, L.; McGuire, S. C.; McIntyre, G.; McIver, J.; McManus, D. J.; McNeill, L.; McRae, T.; McWilliams, S. T.; Meacher, D.; Meadors, G. D.; Mehmet, M.; Meidam, J.; Mejuto-Villa, E.; Melatos, A.; Mendell, G.; Mercer, R. A.; Merilh, E. L.; Merzougui, M.; Meshkov, S.; Messenger, C.; Messick, C.; Metzdorff, R.; Meyers, P. M.; Miao, H.; Michel, C.; Middleton, H.; Mikhailov, E. E.; Milano, L.; Miller, A. L.; Miller, B. B.; Miller, J.; Millhouse, M.; Milovich-Goff, M. C.; Minazzoli, O.; Minenkov, Y.; Ming, J.; Mishra, C.; Mitra, S.; Mitrofanov, V. P.; Mitselmakher, G.; Mittleman, R.; Moffa, D.; Moggi, A.; Mogushi, K.; Mohan, M.; Mohapatra, S. R. P.; Montani, M.; Moore, C. J.; Moraru, D.; Moreno, G.; Morriss, S. R.; Mours, B.; Mow-Lowry, C. M.; Mueller, G.; Muir, A. W.; Mukherjee, Arunava; Mukherjee, D.; Mukherjee, S.; Mukund, N.; Mullavey, A.; Munch, J.; Muñiz, E. A.; Muratore, M.; Murray, P. G.; Napier, K.; Nardecchia, I.; Naticchioni, L.; Nayak, R. K.; Neilson, J.; Nelemans, G.; Nelson, T. J. N.; Nery, M.; Neunzert, A.; Nevin, L.; Newport, J. M.; Newton, G.; Ng, K. K. Y.; Nguyen, T. T.; Nichols, D.; Nielsen, A. B.; Nissanke, S.; Nitz, A.; Noack, A.; Nocera, F.; Nolting, D.; North, C.; Nuttall, L. K.; Oberling, J.; O'Dea, G. D.; Ogin, G. H.; Oh, J. J.; Oh, S. H.; Ohme, F.; Okada, M. A.; Oliver, M.; Oppermann, P.; Oram, Richard J.; O'Reilly, B.; Ormiston, R.; Ortega, L. F.; O'Shaughnessy, R.; Ossokine, S.; Ottaway, D. J.; Overmier, H.; Owen, B. J.; Pace, A. E.; Page, J.; Page, M. A.; Pai, A.; Pai, S. A.; Palamos, J. R.; Palashov, O.; Palomba, C.; Pal-Singh, A.; Pan, Howard; Pan, Huang-Wei; Pang, B.; Pang, P. T. H.; Pankow, C.; Pannarale, F.; Pant, B. C.; Paoletti, F.; Paoli, A.; Papa, M. A.; Parida, A.; Parker, W.; Pascucci, D.; Pasqualetti, A.; Passaquieti, R.; Passuello, D.; Patil, M.; Patricelli, B.; Pearlstone, B. L.; Pedraza, M.; Pedurand, R.; Pekowsky, L.; Pele, A.; Penn, S.; Perez, C. J.; Perreca, A.; Perri, L. M.; Pfeiffer, H. P.; Phelps, M.; Piccinni, O. J.; Pichot, M.; Piergiovanni, F.; Pierro, V.; Pillant, G.; Pinard, L.; Pinto, I. M.; Pirello, M.; Pitkin, M.; Poe, M.; Poggiani, R.; Popolizio, P.; Porter, E. K.; Post, A.; Powell, J.; Prasad, J.; Pratt, J. W. W.; Pratten, G.; Predoi, V.; Prestegard, T.; Prijatelj, M.; Principe, M.; Privitera, S.; Prodi, G. A.; Prokhorov, L. G.; Puncken, O.; Punturo, M.; Puppo, P.; Pürrer, M.; Qi, H.; Quetschke, V.; Quintero, E. A.; Quitzow-James, R.; Raab, F. J.; Rabeling, D. S.; Radkins, H.; Raffai, P.; Raja, S.; Rajan, C.; Rajbhandari, B.; Rakhmanov, M.; Ramirez, K. E.; Ramos-Buades, A.; Rapagnani, P.; Raymond, V.; Razzano, M.; Read, J.; Regimbau, T.; Rei, L.; Reid, S.; Reitze, D. H.; Ren, W.; Reyes, S. D.; Ricci, F.; Ricker, P. M.; Rieger, S.; Riles, K.; Rizzo, M.; Robertson, N. A.; Robie, R.; Robinet, F.; Rocchi, A.; Rolland, L.; Rollins, J. G.; Roma, V. J.; Romano, J. D.; Romano, R.; Romel, C. L.; Romie, J. H.; Rosińska, D.; Ross, M. P.; Rowan, S.; Rüdiger, A.; Ruggi, P.; Rutins, G.; Ryan, K.; Sachdev, S.; Sadecki, T.; Sadeghian, L.; Sakellariadou, M.; Salconi, L.; Saleem, M.; Salemi, F.; Samajdar, A.; Sammut, L.; Sampson, L. M.; Sanchez, E. J.; Sanchez, L. E.; Sanchis-Gual, N.; Sandberg, V.; Sanders, J. R.; Sassolas, B.; Sathyaprakash, B. S.; Saulson, P. R.; Sauter, O.; Savage, R. L.; Sawadsky, A.; Schale, P.; Scheel, M.; Scheuer, J.; Schmidt, J.; Schmidt, P.; Schnabel, R.; Schofield, R. M. S.; Schönbeck, A.; Schreiber, E.; Schuette, D.; Schulte, B. W.; Schutz, B. F.; Schwalbe, S. G.; Scott, J.; Scott, S. M.; Seidel, E.; Sellers, D.; Sengupta, A. S.; Sentenac, D.; Sequino, V.; Sergeev, A.; Shaddock, D. A.; Shaffer, T. J.; Shah, A. A.; Shahriar, M. S.; Shaner, M. B.; Shao, L.; Shapiro, B.; Shawhan, P.; Sheperd, A.; Shoemaker, D. H.; Shoemaker, D. M.; Siellez, K.; Siemens, X.; Sieniawska, M.; Sigg, D.; Silva, A. D.; Singer, L. P.; Singh, A.; Singhal, A.; Sintes, A. M.; Slagmolen, B. J. J.; Smith, B.; Smith, J. R.; Smith, R. J. E.; Somala, S.; Son, E. J.; Sonnenberg, J. A.; Sorazu, B.; Sorrentino, F.; Souradeep, T.; Spencer, A. P.; Srivastava, A. K.; Staats, K.; Staley, A.; Steer, D.; Steinke, M.; Steinlechner, J.; Steinlechner, S.; Steinmeyer, D.; Stevenson, S. P.; Stone, R.; Stops, D. J.; Strain, K. A.; Stratta, G.; Strigin, S. E.; Strunk, A.; Sturani, R.; Stuver, A. L.; Summerscales, T. Z.; Sun, L.; Sunil, S.; Suresh, J.; Sutton, P. J.; Swinkels, B. L.; Szczepańczyk, M. J.; Tacca, M.; Tait, S. C.; Talbot, C.; Talukder, D.; Tanner, D. B.; Tápai, M.; Taracchini, A.; Tasson, J. D.; Taylor, J. A.; Taylor, R.; Tewari, S. V.; Theeg, T.; Thies, F.; Thomas, E. G.; Thomas, M.; Thomas, P.; Thorne, K. A.; Thrane, E.; Tiwari, S.; Tiwari, V.; Tokmakov, K. V.; Toland, K.; Tonelli, M.; Tornasi, Z.; Torres-Forné, A.; Torrie, C. I.; Töyrä, D.; Travasso, F.; Traylor, G.; Trinastic, J.; Tringali, M. C.; Trozzo, L.; Tsang, K. W.; Tse, M.; Tso, R.; Tsukada, L.; Tsuna, D.; Tuyenbayev, D.; Ueno, K.; Ugolini, D.; Unnikrishnan, C. S.; Urban, A. L.; Usman, S. A.; Vahlbruch, H.; Vajente, G.; Valdes, G.; van Bakel, N.; van Beuzekom, M.; Brand, J. F. J. van den; Broeck, C. Van Den; Vander-Hyde, D. C.; van der Schaaf, L.; van Heijningen, J. V.; van Veggel, A. A.; Vardaro, M.; Varma, V.; Vass, S.; Vasúth, M.; Vecchio, A.; Vedovato, G.; Veitch, J.; Veitch, P. J.; Venkateswara, K.; Venugopalan, G.; Verkindt, D.; Vetrano, F.; Viceré, A.; Viets, A. D.; Vinciguerra, S.; Vine, D. J.; Vinet, J. -Y.; Vitale, S.; Vo, T.; Vocca, H.; Vorvick, C.; Vyatchanin, S. P.; Wade, A. R.; Wade, L. E.; Wade, M.; Walet, R.; Walker, M.; Wallace, L.; Walsh, S.; Wang, G.; Wang, H.; Wang, J. Z.; Wang, W. H.; Wang, Y. F.; Ward, R. L.; Warner, J.; Was, M.; Watchi, J.; Weaver, B.; Wei, L. -W.; Weinert, M.; Weinstein, A. J.; Weiss, R.; Wen, L.; Wessel, E. K.; Weßels, P.; Westerweck, J.; Westphal, T.; Wette, K.; Whelan, J. T.; Whitcomb, S. E.; Whiting, B. F.; Whittle, C.; Wilken, D.; Williams, D.; Williams, R. D.; Williamson, A. R.; Willis, J. L.; Willke, B.; Wimmer, M. H.; Winkler, W.; Wipf, C. C.; Wittel, H.; Woan, G.; Woehler, J.; Wofford, J.; Wong, K. W. K.; Worden, J.; Wright, J. L.; Wu, D. S.; Wysocki, D. M.; Xiao, S.; Yamamoto, H.; Yancey, C. C.; Yang, L.; Yap, M. J.; Yazback, M.; Yu, Hang; Yu, Haocun; Yvert, M.; Zadrożny, A.; Zanolin, M.; Zelenova, T.; Zendri, J. -P.; Zevin, M.; Zhang, L.; Zhang, M.; Zhang, T.; Zhang, Y. -H.; Zhao, C.; Zhou, M.; Zhou, Z.; Zhu, S. J.; Zhu, X. J.; Zimmerman, A. B.; Zucker, M. E.; Zweizig, J.; Foley, R. J.; Coulter, D. A.; Drout, M. R.; Kasen, D.; Kilpatrick, C. D.; Madore, B. F.; Murguia-Berthier, A.; Pan, Y. -C.; Piro, A. L.; Prochaska, J. X.; Ramirez-Ruiz, E.; Rest, A.; Rojas-Bravo, C.; Shappee, B. J.; Siebert, M. R.; Simon, J. D.; Ulloa, N.; Annis, J.; Soares-Santos, M.; Brout, D.; Scolnic, D.; Diehl, H. T.; Frieman, J.; Berger, E.; Alexander, K. D.; Allam, S.; Balbinot, E.; Blanchard, P.; Butler, R. E.; Chornock, R.; Cook, E. R.; Cowperthwaite, P.; Drlica-Wagner, A.; Drout, M. R.; Durret, F.; Eftekhari, T.; Finley, D. A.; Fong, W.; Fryer, C. L.; García-Bellido, J.; Gill, M. S . S.; Gruendl, R. A.; Hanna, C.; Hartley, W.; Herner, K.; Huterer, D.; Kasen, D.; Kessler, R.; Li, T. S.; Lin, H.; Lopes, P. A. A.; Lourenço, A. C. C.; Margutti, R.; Marriner, J.; Marshall, J. L.; Matheson, T.; Medina, G. E.; Metzger, B. D.; Muñoz, R. R.; Muir, J.; Nicholl, M.; Nugent, P.; Palmese, A.; Paz-Chinchón, F.; Quataert, E.; Sako, M.; Sauseda, M.; Schlegel, D. J.; Secco, L. F.; Smith, N.; Sobreira, F.; Stebbins, A.; Villar, V. A.; Vivas, A. K.; Wester, W.; Williams, P. K. G.; Yanny, B.; Zenteno, A.; Abbott, T. M. C.; Abdalla, F. B.; Bechtol, K.; Benoit-Lévy, A.; Bertin, E.; Bridle, S. L.; Brooks, D.; Buckley-Geer, E.; Burke, D. L.; Rosell, A. Carnero; Kind, M. Carrasco; Carretero, J.; Castander, F. J.; Cunha, C. E.; D'Andrea, C. B.; da Costa, L. N.; Davis, C.; DePoy, D. L.; Desai, S.; Dietrich, J. P.; Estrada, J.; Fernandez, E.; Flaugher, B.; Fosalba, P.; Gaztanaga, E.; Gerdes, D. W.; Giannantonio, T.; Goldstein, D. A.; Gruen, D.; Gutierrez, G.; Hartley, W. G.; Honscheid, K.; Jain, B.; James, D. J.; Jeltema, T.; Johnson, M. W. G.; Kent, S.; Krause, E.; Kron, R.; Kuehn, K.; Kuhlmann, S.; Kuropatkin, N.; Lahav, O.; Lima, M.; Maia, M. A. G.; March, M.; Miller, C. J.; Miquel, R.; Neilsen, E.; Nord, B.; Ogando, R. L. C.; Plazas, A. A.; Romer, A. K.; Roodman, A.; Rykoff, E. S.; Sanchez, E.; Scarpine, V.; Schubnell, M.; Sevilla-Noarbe, I.; Smith, M.; Smith, R. C.; Suchyta, E.; Tarle, G.; Thomas, D.; Thomas, R. C.; Troxel, M. A.; Tucker, D. L.; Vikram, V.; Walker, A. R.; Weller, J.; Zhang, Y.; Haislip, J. B.; Kouprianov, V. V.; Reichart, D. E.; Tartaglia, L.; Sand, D. J.; Valenti, S.; Yang, S.; Arcavi, Iair; Hosseinzadeh, Griffin; Howell, D. Andrew; McCully, Curtis; Poznanski, Dovi; Vasylyev, Sergiy; Tanvir, N. R.; Levan, A. J.; Hjorth, J.; Cano, Z.; Copperwheat, C.; de Ugarte-Postigo, A.; Evans, P. A.; Fynbo, J. P. U.; González-Fernández, C.; Greiner, J.; Irwin, M.; Lyman, J.; Mandel, I.; McMahon, R.; Milvang-Jensen, B.; O'Brien, P.; Osborne, J. P.; Perley, D. A.; Pian, E.; Palazzi, E.; Rol, E.; Rosetti, S.; Rosswog, S.; Rowlinson, A.; Schulze, S.; Steeghs, D. T. H.; Thöne, C. C.; Ulaczyk, K.; Watson, D.; Wiersema, K.; Lipunov, V. M.; Gorbovskoy, E.; Kornilov, V. G.; Tyurina, N .; Balanutsa, P.; Vlasenko, D.; Gorbunov, I.; Podesta, R.; Levato, H.; Saffe, C.; Buckley, D. A. H.; Budnev, N. M.; Gress, O.; Yurkov, V.; Rebolo, R.; Serra-Ricart, M.
Comments: 26 pages, 5 figures, Nature in press. For more information see https://dcc.ligo.org/LIGO-P1700296/public
Submitted: 2017-10-16
The detection of GW170817 in both gravitational waves and electromagnetic waves heralds the age of gravitational-wave multi-messenger astronomy. On 17 August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong signal from the merger of a binary neutron-star system. Less than 2 seconds after the merger, a gamma-ray burst (GRB 170817A) was detected within a region of the sky consistent with the LIGO-Virgo-derived location of the gravitational-wave source. This sky region was subsequently observed by optical astronomy facilities, resulting in the identification of an optical transient signal within $\sim 10$ arcsec of the galaxy NGC 4993. These multi-messenger observations allow us to use GW170817 as a standard siren, the gravitational-wave analog of an astronomical standard candle, to measure the Hubble constant. This quantity, which represents the local expansion rate of the Universe, sets the overall scale of the Universe and is of fundamental importance to cosmology. Our measurement combines the distance to the source inferred purely from the gravitational-wave signal with the recession velocity inferred from measurements of the redshift using electromagnetic data. This approach does not require any form of cosmic "distance ladder;" the gravitational wave analysis can be used to estimate the luminosity distance out to cosmological scales directly, without the use of intermediate astronomical distance measurements. We determine the Hubble constant to be $70.0^{+12.0}_{-8.0} \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ (maximum a posteriori and 68% credible interval). This is consistent with existing measurements, while being completely independent of them. Additional standard-siren measurements from future gravitational-wave sources will provide precision constraints of this important cosmological parameter.
[70]  oai:arXiv.org:1710.05840  [pdf] - 1589828
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. II. UV, Optical, and Near-IR Light Curves and Comparison to Kilonova Models
Cowperthwaite, P. S.; Berger, E.; Villar, V. A.; Metzger, B. D.; Nicholl, M.; Chornock, R.; Blanchard, P. K.; Fong, W.; Margutti, R.; Soares-Santos, M.; Alexander, K. D.; Allam, S.; Annis, J.; Brout, D.; Brown, D. A.; Butler, R. E.; Chen, H. -Y.; Diehl, H. T.; Doctor, Z.; Drout, M. R.; Eftekhari, T.; Farr, B.; Finley, D. A.; Foley, R. J.; Frieman, J. A.; Fryer, C. L.; García-Bellido, J.; Gill, M. S. S.; Guillochon, J.; Herner, K.; Holz, D. E.; Kasen, D.; Kessler, R.; Marriner, J.; Matheson, T.; Neilsen,, E. H.; Quataert, E.; Palmese, A.; Rest, A.; Sako, M.; Scolnic, D. M.; Smith, N.; Tucker, D. L.; Williams, P. K. G.; Balbinot, E.; Carlin, J. L.; Cook, E. R.; Durret, F.; Li, T. S.; Lopes, P. A. A.; Lourenço, A. C. C.; Marshall, J. L.; Medina, G. E.; Muir, J.; Muñoz, R. R.; Sauseda, M.; Schlegel, D. J.; Secco, L. F.; Vivas, A. K.; Wester, W.; Zenteno, A.; Zhang, Y.; Abbott, T. M. C.; Banerji, M.; Bechtol, K.; Benoit-Lévy, A.; Bertin, E.; Buckley-Geer, E.; Burke, D. L.; Capozzi, D.; Rosell, A. Carnero; Kind, M. Carrasco; Castander, F. J.; Crocce, M.; Cunha, C. E.; D'Andrea, C. B.; da Costa, L. N.; Davis, C.; DePoy, D. L.; Desai, S.; Dietrich, J. P.; Drlica-Wagner, A.; Eifler, T. F.; Evrard, A. E.; Fernandez, E.; Flaugher, B.; Fosalba, P.; Gaztanaga, E.; Gerdes, D. W.; Giannantonio, T.; Goldstein, D. A.; Gruen, D.; Gruendl, R. A.; Gutierrez, G.; Honscheid, K.; Jain, B.; James, D. J.; Jeltema, T.; Johnson, M. W. G.; Johnson, M. D.; Kent, S.; Krause, E.; Kron, R.; Kuehn, K.; Kuropatkin, N.; Lahav, O.; Lima, M.; Lin, H.; Maia, M. A. G.; March, M.; Martini, P.; McMahon, R. G.; Menanteau, F.; Miller, C. J.; Miquel, R.; Mohr, J. J.; Neilsen, E.; Nichol, R. C.; Ogando, R. L. C.; Plazas, A. A.; Roe, N.; Romer, A. K.; Roodman, A.; Rykoff, E. S.; Sanchez, E.; Scarpine, V.; Schindler, R.; Schubnell, M.; Sevilla-Noarbe, I.; Smith, M.; Smith, R. C.; Sobreira, F.; Suchyta, E.; Swanson, M. E. C.; Tarle, G.; Thomas, D.; Thomas, R. C.; Troxel, M. A.; Vikram, V.; Walker, A. R.; Wechsler, R. H.; Weller, J.; Yanny, B.; Zuntz, J.
Comments: 13 Pages, 3 Figures, 2 Tables. ApJL, In Press. Keywords: GW170817, LVC
Submitted: 2017-10-16
We present UV, optical, and NIR photometry of the first electromagnetic counterpart to a gravitational wave source from Advanced LIGO/Virgo, the binary neutron star merger GW170817. Our data set extends from the discovery of the optical counterpart at $0.47$ days to $18.5$ days post-merger, and includes observations with the Dark Energy Camera (DECam), Gemini-South/FLAMINGOS-2 (GS/F2), and the {\it Hubble Space Telescope} ({\it HST}). The spectral energy distribution (SED) inferred from this photometry at $0.6$ days is well described by a blackbody model with $T\approx 8300$ K, a radius of $R\approx 4.5\times 10^{14}$ cm (corresponding to an expansion velocity of $v\approx 0.3c$), and a bolometric luminosity of $L_{\rm bol}\approx 5\times10^{41}$ erg s$^{-1}$. At $1.5$ days we find a multi-component SED across the optical and NIR, and subsequently we observe rapid fading in the UV and blue optical bands and significant reddening of the optical/NIR colors. Modeling the entire data set we find that models with heating from radioactive decay of $^{56}$Ni, or those with only a single component of opacity from $r$-process elements, fail to capture the rapid optical decline and red optical/NIR colors. Instead, models with two components consistent with lanthanide-poor and lanthanide-rich ejecta provide a good fit to the data, the resulting "blue" component has $M_\mathrm{ej}^\mathrm{blue}\approx 0.01$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{blue}\approx 0.3$c, and the "red" component has $M_\mathrm{ej}^\mathrm{red}\approx 0.04$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{red}\approx 0.1$c. These ejecta masses are broadly consistent with the estimated $r$-process production rate required to explain the Milky Way $r$-process abundances, providing the first evidence that BNS mergers can be a dominant site of $r$-process enrichment.
[71]  oai:arXiv.org:1706.09424  [pdf] - 1614751
COSMOGRAIL XVI: Time delays for the quadruply imaged quasar DES J0408-5354 with high-cadence photometric monitoring
Comments: 9 pages, 5 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2017-06-28, last modified: 2017-10-10
We present time-delay measurements for the new quadruply imaged quasar DES J0408-5354, the first quadruply imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image. This data quality allows us to catch small photometric variations (a few mmag rms) of the quasar, acting on temporal scales much shorter than microlensing, hence making the time delay measurement very robust against microlensing. In only 7 months we measure very accurately one of the time delays in DES J0408-5354: Dt(AB) = -112.1 +- 2.1 days (1.8%) using only the MPIA 2.2m data. In combination with data taken with the 1.2m Euler Swiss telescope, we also measure two delays involving the D component of the system Dt(AD) = -155.5 +- 12.8 days (8.2%) and Dt(BD) = -42.4 +- 17.6 days (41%), where all the error bars include systematics. Turning these time delays into cosmological constraints will require deep HST imaging or ground-based Adaptive Optics (AO), and information on the velocity field of the lensing galaxy.
[72]  oai:arXiv.org:1710.02144  [pdf] - 1682478
An Empirical Study of Contamination in Deep, Rapid, and Wide-Field Optical Follow-Up of Gravitational Wave Events
Comments: 19 pages, 8 figures, 5 tables, Submitted to ApJ (Replacement 2017-10-08: Updated Author List)
Submitted: 2017-10-05, last modified: 2017-10-08
We present an empirical study of contamination in deep, rapid, and wide-field optical follow-up searches of GW sources from aLIGO. We utilize dedicated observations during four nights of imaging with DECam. Our search covered $\sim56$ deg$^2$, with two visits per night separated by $\approx 3$~hours, in $i$- and $z$-band, followed by an additional set of $griz$ images three weeks later to serve as reference images for subtraction, and for the purpose of identifying galaxy and stellar counterparts for any transient sources. We achieve $5\sigma$ point-source limiting magnitudes of $i \approx 23.5$ and $z \approx 22.4$ mag in the coadded single-epoch images. We conduct a search for transient objects that can mimic the $i-z$ color behavior of both red ($i-z > 0.5$~mag) and blue ($i-z < 0$~mag) kilonova emission, finding 11 and 10 contaminants, respectively. Independent of color, we identify 48 transients of interest. Additionally, we leverage the rapid cadence of our observations to search for sources with characteristic timescales of $\approx1$ day and $\approx3$ hours, finding no potential contaminants. We assess the efficiency of our pipeline and search methodology with injected point sources, finding that we are 90\% (60\%) efficient when searching for red (blue) kilonova-like sources to a limiting magnitude of $i \lesssim 22.5$ mag. Applying these efficiencies, we derive sky rates for kilonova contaminants in the red and blue regimes of $\mathcal{R}_{\rm red} \approx 0.16$ deg$^{-2}$ and $\mathcal{R}_{\rm blue} \approx 0.80$ deg$^{-2}$. The total contamination rate, independent of color, is $\mathcal{R}_{\rm all} \approx 1.79$ deg$^{-2}$. We compare our derived results to optical follow-up searches of the GW events GW150914 and GW151226 and comment on the outlook for GW follow-up searches as additional GW detectors (e.g., KAGRA, LIGO India) come online in the next decade.
[73]  oai:arXiv.org:1705.08455  [pdf] - 1583792
A Reverse Shock and Unusual Radio Properties in GRB 160625B
Comments: 17 pages, 8 figures, 3 tables. Accepted for publication in ApJ
Submitted: 2017-05-23, last modified: 2017-09-06
We present multi-wavelength observations and modeling of the exceptionally bright long $\gamma$-ray burst GRB 160625B. The optical and X-ray data are well-fit by synchrotron emission from a collimated blastwave with an opening angle of $\theta_j\approx 3.6^\circ$ and kinetic energy of $E_K\approx 2\times10^{51}$ erg, propagating into a low density ($n\approx 5\times10^{-5}$ cm$^{-3}$) medium with a uniform profile. The forward shock is sub-dominant in the radio band; instead, the radio emission is dominated by two additional components. The first component is consistent with emission from a reverse shock, indicating an initial Lorentz factor of $\Gamma_0\gtrsim 100$ and an ejecta magnetization of $R_B\approx 1-100$. The second component exhibits peculiar spectral and temporal evolution and is most likely the result of scattering of the radio emission by the turbulent Milky Way interstellar medium (ISM). Such scattering is expected in any sufficiently compact extragalactic source and has been seen in GRBs before, but the large amplitude and long duration of the variability seen here are qualitatively more similar to extreme scattering events previously observed in quasars, rather than normal interstellar scintillation effects. High-cadence, broadband radio observations of future GRBs are needed to fully characterize such effects, which can sensitively probe the properties of the ISM and must be taken into account before variability intrinsic to the GRB can be interpreted correctly.
[74]  oai:arXiv.org:1707.05784  [pdf] - 1697908
A VLA Study of High-redshift GRBs II - The Complex Radio Afterglow of GRB 140304A: Shell Collisions and Two Reverse Shocks
Comments: 24 pages, 17 figures, 7 tables; submitted to ApJ
Submitted: 2017-07-18
We present detailed multi-frequency, multi-epoch radio observations of GRB 140304A at $z=5.283$ from 1 to 86 GHz and 0.45 d to 89 d. The radio and mm data exhibit unusual multiple spectral components, which cannot be simply explained by standard forward and reverse shock scenarios. Through detailed multi-wavelength analysis spanning radio to X-rays, we constrain the forward shock parameters to $E_{\rm K, iso}\approx4.9\times10^{54}\,$erg, $A_* \approx 2.6\times10^{-2}$, $\epsilon_{\rm e}\approx2.5\times10^{-2}$, $\epsilon_{\rm B}\approx5.9\times10^{-2}$, $p\approx2.6$, and $\theta_{\rm jet}\approx1.1^{\circ}$, yielding a beaming corrected $\gamma$-ray and kinetic energy, $E_{\gamma}\approx2.3\times10^{49}\,$erg and $E_{\rm K}\approx9.5\times10^{50}\,$erg, respectively. We model the excess radio emission as due to a combination of a late-time reverse shock (RS) launched by a shell collision, which also produces a re-brightening in the X-rays at $\approx0.26\,$d, and either a standard RS or diffractive interstellar scintillation. Under the standard RS interpretation, we invoke consistency arguments between the forward and reverse shocks to derive a deceleration time, $t_{\rm dec}\approx100\,$s, the ejecta Lorentz factor, $\Gamma(t_{\rm dec})\approx300$, and a low RS magnetization, $R_{\rm B}\approx0.6$. Our observations highlight both the power of radio observations in capturing RS emission and thus constraining the properties of GRB ejecta and central engines, and the challenge presented by interstellar scintillation in conclusively identifying RS emission in GRB radio afterglows.
[75]  oai:arXiv.org:1707.04264  [pdf] - 1790554
Radio Emission from Ultra-Cool Dwarfs
Comments: 19 pages; submitted chapter to the Handbook of Exoplanets, eds. Hans J. Deeg and Juan Antonio Belmonte (Springer-Verlag)
Submitted: 2017-07-13
The 2001 discovery of radio emission from ultra-cool dwarfs (UCDs), the very low-mass stars and brown dwarfs with spectral types of ~M7 and later, revealed that these objects can generate and dissipate powerful magnetic fields. Radio observations provide unparalleled insight into UCD magnetism: detections extend to brown dwarfs with temperatures <1000 K, where no other observational probes are effective. The data reveal that UCDs can generate strong (kG) fields, sometimes with a stable dipolar structure; that they can produce and retain nonthermal plasmas with electron acceleration extending to MeV energies; and that they can drive auroral current systems resulting in significant atmospheric energy deposition and powerful, coherent radio bursts. Still to be understood are the underlying dynamo processes, the precise means by which particles are accelerated around these objects, the observed diversity of magnetic phenomenologies, and how all of these factors change as the mass of the central object approaches that of Jupiter. The answers to these questions are doubly important because UCDs are both potential exoplanet hosts, as in the TRAPPIST-1 system, and analogues of extrasolar giant planets themselves.
[76]  oai:arXiv.org:1704.00022  [pdf] - 1582283
Empirical constraints on the origin of fast radio bursts: volumetric rates and host galaxy demographics as a test of millisecond magnetar connection
Comments:
Submitted: 2017-03-31, last modified: 2017-07-11
The localization of the repeating FRB 121102 to a low-metallicity dwarf galaxy at $z=0.193$, and its association with a quiescent radio source, suggests the possibility that FRBs originate from magnetars, formed by the unusual supernovae in such galaxies. We investigate this via a comparison of magnetar birth rates, the FRB volumetric rate, and host galaxy demographics. We calculate average volumetric rates of possible millisecond magnetar production channels such as superluminous supernovae (SLSNe), long and short gamma-ray bursts (GRBs), and general magnetar production via core-collapse supernovae. For each channel we also explore the expected host galaxy demographics using their known properties. We determine for the first time the number density of FRB emitters (the product of their volumetric birthrate and lifetime), $R_{\rm FRB}\tau\approx 10^4$Gpc$^{-3}$, assuming that FRBs are predominantly emitted from repetitive sources similar to FRB 121102 and adopting a beaming factor of 0.1. By comparing rates we find that production via rare channels (SLSNe, GRBs) implies a typical FRB lifetime of $\approx$30-300 yr, in good agreement with other lines of argument. The total energy emitted over this time is consistent with the available energy stored in the magnetic field. On the other hand, any relation to magnetars produced via normal core-collapse supernovae leads to a very short lifetime of $\approx$0.5yr, in conflict with both theory and observation. We demonstrate that due to the diverse host galaxy distributions of the different progenitor channels, many possible sources of FRB birth can be ruled out with $\lesssim 10$ host galaxy identifications. Conversely, targeted searches of galaxies that have previously hosted decades-old SLSNe and GRBs may be a fruitful strategy for discovering new FRBs and related quiescent radio sources, and determining the nature of their progenitors.
[77]  oai:arXiv.org:1707.02107  [pdf] - 1585643
Radio Exploration of Planetary Habitability: Conference Summary
Comments: Five pages; conference Web site: http://aas.org/meetings/aastcs/radiohab
Submitted: 2017-07-07
Radio Exploration of Planetary Habitability was the fifth in the series of American Astronomical Society's Topical Conference Series. Notable aspects of the conference included the interdisciplinary nature of both the topics and the intellectual breadth of the participants, the diversity of approaches to studying this topic presented by recent discoveries and of the participants themselves, the expanding meaning of the topic of "star-planet interactions," and the expectation of an increasingly statistical approach to the topic. Potential areas of future research include the actual extent to which planetary magnetic fields shield planetary atmospheres; the planetary dynamo process itself, particularly once multiple extrasolar planetary magnetic fields are confirmed; and "planet-star interactions." A major major topic of the conference concerned observational opportunities, highlighted by a number of new or upcoming, specialized observatories to observe exoplanets especially at radio wavelengths. This article summarizes these main points of the conference and expands briefly upon these potential avenues for future investigation. A future meeting on this topic, given the variety of data sets being generated over the next few years, is warranted.
[78]  oai:arXiv.org:1703.08745  [pdf] - 1582130
K2 Ultracool Dwarfs Survey II: The White Light Flare Rate of Young Brown Dwarfs
Comments: Accepted version, The Astrophysical Journal
Submitted: 2017-03-25, last modified: 2017-07-04
We use Kepler K2 Campaign 4 short-cadence (one-minute) photometry to measure white light flares in the young, moving group brown dwarfs 2MASS J03350208+2342356 (2M0335+23) and 2MASS J03552337+1133437 (2M0355+11), and report on long-cadence (thirty-minute) photometry of a superflare in the Pleiades M8 brown dwarf CFHT-PL-17. The rotation period (5.24 hr) and projected rotational velocity ($45$ km s$^{-1}$) confirm 2M0335+23 is inflated ($R \ge 0.20 R_\odot$) as predicted for a $0.06M_\odot$, 26-Myr old brown dwarf $\beta $Pic moving group member. We detect 22 white light flares on 2M0335+23. The flare frequency distribution follows a power-law distribution with slope $-\alpha = -1.8 \pm 0.2$ over the range $10^{31}$ to $10^{33}$ erg. This slope is similar to that observed in the Sun and warmer flare stars, and is consistent with lower energy flares in previous work on M6-M8 very-low-mass stars; taken the two datasets together, the flare frequency distribution for ultracool dwarfs is a power law over 4.3 orders of magnitude. The superflare ($2.6\times10^{34}$ erg) on CFHT-PL-17 shows higher energy flares are possible. We detect no flares down to a limit of $2 \times 10^{30}$ erg in the nearby L$5\gamma$ AB Dor Moving Group brown dwarf 2M0355+11, consistent with the view that fast magnetic reconnection is suppressed in cool atmospheres. We discuss two multi-peaked flares observed in 2M0335+23, and argue that these complex flares can be understood as sympathetic flares, in which a fast-mode MHD waves similar to EUV waves in the Sun trigger magnetic reconnection in different active regions.
[79]  oai:arXiv.org:1706.01506  [pdf] - 1656087
Discovery of three strongly lensed quasars in the Sloan Digital Sky Survey
Comments: 9 pages, 3 figures, 4 tables; submitted to MNRAS
Submitted: 2017-06-05
We present the discovery of 3 quasar lenses in the Sloan Digital Sky Survey (SDSS), selected using two novel photometry-based selection techniques. The J0941+0518 system, with two point sources separated by 5.46" on either side of a galaxy, has source and lens redshifts $z_s = 1.54$ and $z_l = 0.343$. The AO-assisted images of J2211+1929 show two point sources separated by 1.04", corresponding to the same quasar at $z_s = 1.07,$ besides the lens galaxy and Einstein ring. Images of J2257+2349 show two point sources separated by 1.67" on either side of an E/S0 galaxy. The extracted spectra show two images of the same quasar at redshift $z_s = 2.10$. In total, the two selection techniques identified 309 lens candidates, including 47 known lenses, and 6 previously ruled out candidates. 55 of the remaining candidates were observed using NIRC2 and ESI at Keck Observatory, EFOSC2 at the ESO-NTT (La Silla), and SAM and the Goodman spectrograph at SOAR. Of the candidates observed, 3 were confirmed as lenses, 36 were ruled out, and 16 remain inconclusive. Taking into account that we recovered known lenses, this gives us a success rate of at least 50/309 (16%). This initial campaign demonstrates the power of purely photometric selection techniques in finding lensed quasars. Developing and refining these techniques is essential for efficient identification of these rare lenses in ongoing and future photometric surveys.
[80]  oai:arXiv.org:1705.07553  [pdf] - 1598048
A Multi-telescope Campaign on FRB 121102: Implications for the FRB Population
Comments: 17 pages, 7 figures. Submitted to AAS Journals
Submitted: 2017-05-22
We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.
[81]  oai:arXiv.org:1704.05720  [pdf] - 1565396
Observation of a Deep Visual "Eclipse" in the WC9-Type Wolf-Rayet Star, WR 76
Comments: 9 pages, 2 figures, accepted for publication in the Journal of the AAVSO
Submitted: 2017-04-19
The WC9-Type Wolf-Rayet star WR 76 is one of the most prolific dust makers identified from its infrared emission. WR 76 experienced a deep fading eclipse in 2016. The ~3.1 magnitude depth of the eclipse exceeds fadings in similar eclipses observed in WR stars thus far. Conclusions from recent and earlier analyses of eclipses observed suggests that WR 76 may be a prolific eclipser.
[82]  oai:arXiv.org:1701.03209  [pdf] - 1652278
The Hydrogen Epoch of Reionization Array Dish III: Measuring Chromaticity of Prototype Element with Reflectometry
Comments: 13 pages, 10 figures, Submitted to APJ
Submitted: 2017-01-11, last modified: 2017-03-14
The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna.
[83]  oai:arXiv.org:1703.04692  [pdf] - 1547792
Enabling New ALMA Science with Improved Support for Time-Domain Observations
Comments: 9 pages; whitepaper submitted to the ALMA Science Advisory Council; corresponding author P. K. G. Williams (pwilliams@cfa.harvard.edu)
Submitted: 2017-03-14
While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here as "time-domain" observations. In this whitepaper we identify some of these policies, describe the scientific applications they impair, and suggest changes that would increase ALMA's science impact in Cycle 6 and beyond. Parties interested in time-domain science with ALMA are encouraged to join the ALMA Time-domain Special Interest Group (ATSIG) by signing up for the ATSIG mailing list at https://groups.google.com/group/alma-td-sig .
[84]  oai:arXiv.org:1611.07080  [pdf] - 1553934
K2 Ultracool Dwarfs Survey I: Photometry of an L Dwarf Superflare
Comments: Accepted to ApJ
Submitted: 2016-11-21, last modified: 2017-02-18
We report on K2 Campaign 8 measurements of a huge white light flare on the L1 dwarf SDSSp J005406.55-003101.8 (EPIC 220186653). The source is a typical L1 dwarf at a distance of $\sim50$ pc, probably an old hydrogen-burning star rather than a young brown dwarf. In the long (30-minute) cadence photometry, the flare peak is 21 times the flux of the stellar photosphere in the broad optical Kepler filter, which we estimate corresponds to $\Delta V \approx -7.1$. The total equivalent duration of the flare is 15.4 hr. We estimate the total bolometric energy of the flare was $4 \times 10^{33}$ erg, more powerful than the previously reported Kepler white light flares for the L1 dwarf WISEP J190648.47+401106.8, but weaker than the $\Delta V = -11$ L0 dwarf superflare ASASSN-16ae. The initial (impulsive) cooling phase is too rapid to resolve with our 30-minute cadence data, but after one hour the gradual cooling phase has an exponential time constant of 1.8 hours. We use template fitting to estimate that the full-time-width-at-half-amplitude of the light curve is $<10$ minutes and that the true flare maximum reached $\sim70$ times the stellar photosphere, or $\Delta V \approx -8$. This flare is comparable to the most powerful Kepler flares observed on the active M4 dwarf GJ 1243.
[85]  oai:arXiv.org:1612.03821  [pdf] - 1533279
Population mixtures and searches of lensed and extended quasars across photometric surveys
Comments: 17 pages, 10 figures, 3 tables; Accepted for publication in MNRAS
Submitted: 2016-12-12
Wide-field photometric surveys enable searches of rare yet interesting objects, such as strongly lensed quasars or quasars with a bright host galaxy. Past searches for lensed quasars based on their optical and near infrared properties have relied on photometric cuts and spectroscopic pre-selection (as in the Sloan Quasar Lens Search), or neural networks applied to photometric samples. These methods rely on cuts in morphology and colours, with the risk of losing many interesting objects due to scatter in their population properties, restrictive training sets, systematic uncertainties in catalog-based magnitudes, and survey-to-survey photometric variations. Here, we explore the performance of a Gaussian Mixture Model to separate point-like quasars, quasars with an extended host, and strongly lensed quasars using griz psf and model magnitudes and WISE W1, W2. The choice of optical magnitudes is due to their presence in all current and upcoming releases of wide-field surveys, whereas UV information is not always available. We then assess the contamination from blue galaxies and the role of additional features such as W3 magnitudes or psf-model terms as morphological information. As a demonstration, we conduct a search in a random 10% of the SDSS footprint, and we provide the catalog of the 43 SDSS object with the highest `lens' score in our selection that survive visual inspection, and are spectroscopically confirmed to host active nuclei. We inspect archival data and find images of 5/43 objects in the Hubble Legacy Archive, including 2 known lenses. The code and materials are available to facilitate follow-up.
[86]  oai:arXiv.org:1610.04587  [pdf] - 1532067
On the Progenitor of the Type IIb Supernova 2016gkg
Comments: 8 pages, 5 figures, submitted to MNRAS
Submitted: 2016-10-14
We present a detection in pre-explosion Hubble Space Telescope (HST) imaging of a point source consistent with being the progenitor star of the Type IIb supernova (SN IIb) 2016gkg. Post-explosion imaging from the Keck Adaptive Optics system was used to perform relative astrometry between the Keck and HST imaging. We identify a single point source in the HST images coincident with the SN position to 0.89-sigma. The HST photometry is consistent with the progenitor star being an A0Ia star with T=9500 K and log (L/Lsun)=5.15. We find that the SN 2016gkg progenitor star appears more consistent with binary than single-star evolutionary models. In addition, early-time light curve data from SN 2016gkg revealed a rapid rise in luminosity within ~0.4 days of non-detection limits, consistent with models of the cooling phase after shock break-out. We use these data to determine an explosion date of 20.15 September 2016 and progenitor star radius of log (R/Rsun)=2.41, which agrees with photometry from the progenitor star. Our findings are also consistent with detections of other SNe IIb progenitor stars, although more luminous and bluer than most other examples.
[87]  oai:arXiv.org:1610.03159  [pdf] - 1494915
The Astropy Problem
Muna, Demitri; Alexander, Michael; Allen, Alice; Ashley, Richard; Asmus, Daniel; Azzollini, Ruyman; Bannister, Michele; Beaton, Rachael; Benson, Andrew; Berriman, G. Bruce; Bilicki, Maciej; Boyce, Peter; Bridge, Joanna; Cami, Jan; Cangi, Eryn; Chen, Xian; Christiny, Nicholas; Clark, Christopher; Collins, Michelle; Comparat, Johan; Cook, Neil; Croton, Darren; Davids, Isak Delberth; Depagne, Éric; Donor, John; Santos, Leonardo A. dos; Douglas, Stephanie; Du, Alan; Durbin, Meredith; Erb, Dawn; Faes, Daniel; Fernández-Trincado, J. G.; Foley, Anthony; Fotopoulou, Sotiria; Frimann, Søren; Frinchaboy, Peter; Garcia-Dias, Rafael; Gawryszczak, Artur; George, Elizabeth; Gonzalez, Sebastian; Gordon, Karl; Gorgone, Nicholas; Gosmeyer, Catherine; Grasha, Katie; Greenfield, Perry; Grellmann, Rebekka; Guillochon, James; Gurwell, Mark; Haas, Marcel; Hagen, Alex; Haggard, Daryl; Haines, Tim; Hall, Patrick; Hellwing, Wojciech; Herenz, Edmund Christian; Hinton, Samuel; Hlozek, Renee; Hoffman, John; Holman, Derek; Holwerda, Benne Willem; Horton, Anthony; Hummels, Cameron; Jacobs, Daniel; Jensen, Jens Juel; Jones, David; Karick, Arna; Kelley, Luke; Kenworthy, Matthew; Kitchener, Ben; Klaes, Dominik; Kohn, Saul; Konorski, Piotr; Krawczyk, Coleman; Kuehn, Kyler; Kuutma, Teet; Lam, Michael T.; Lane, Richard; Liske, Jochen; Lopez-Camara, Diego; Mack, Katherine; Mangham, Sam; Mao, Qingqing; Marsh, David J. E.; Mateu, Cecilia; Maurin, Loïc; McCormac, James; Momcheva, Ivelina; Monteiro, Hektor; Mueller, Michael; Munoz, Roberto; Naidu, Rohan; Nelson, Nicholas; Nitschelm, Christian; North, Chris; Nunez-Iglesias, Juan; Ogaz, Sara; Owen, Russell; Parejko, John; Patrício, Vera; Pepper, Joshua; Perrin, Marshall; Pickering, Timothy; Piscionere, Jennifer; Pogge, Richard; Poleski, Radek; Pourtsidou, Alkistis; Price-Whelan, Adrian M.; Rawls, Meredith L.; Read, Shaun; Rees, Glen; Rein, Hanno; Rice, Thomas; Riemer-Sørensen, Signe; Rusomarov, Naum; Sanchez, Sebastian F.; Santander-García, Miguel; Sarid, Gal; Schoenell, William; Scholz, Aleks; Schuhmann, Robert L.; Schuster, William; Scicluna, Peter; Seidel, Marja; Shao, Lijing; Sharma, Pranav; Shulevski, Aleksandar; Shupe, David; Sifón, Cristóbal; Simmons, Brooke; Sinha, Manodeep; Skillen, Ian; Soergel, Bjoern; Spriggs, Thomas; Srinivasan, Sundar; Stevens, Abigail; Streicher, Ole; Suchyta, Eric; Tan, Joshua; Telford, O. Grace; Thomas, Romain; Tonini, Chiara; Tremblay, Grant; Tuttle, Sarah; Urrutia, Tanya; Vaughan, Sam; Verdugo, Miguel; Wagner, Alexander; Walawender, Josh; Wetzel, Andrew; Willett, Kyle; Williams, Peter K. G.; Yang, Guang; Zhu, Guangtun; Zonca, Andrea
Comments:
Submitted: 2016-10-10
The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.
[88]  oai:arXiv.org:1606.07473  [pdf] - 1542684
Hydrogen Epoch of Reionization Array (HERA)
Comments: 26 pages, 24 figures, 2 tables
Submitted: 2016-06-23, last modified: 2016-09-12
The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology.
[89]  oai:arXiv.org:1606.04538  [pdf] - 1477626
A DECam Search for an Optical Counterpart to the LIGO Gravitational Wave Event GW151226
Comments: 7 Pages, 2 Figures, 1 Table. Accepted to ApJL, updated to better match published version
Submitted: 2016-06-14, last modified: 2016-09-08
We report the results of a Dark Energy Camera (DECam) optical follow-up of the gravitational wave (GW) event GW151226, discovered by the Advanced LIGO detectors. Our observations cover 28.8 deg$^2$ of the localization region in the $i$ and $z$ bands (containing 3% of the BAYESTAR localization probability), starting 10 hours after the event was announced and spanning four epochs at $2-24$ days after the GW detection. We achieve $5\sigma$ point-source limiting magnitudes of $i\approx21.7$ and $z\approx21.5$, with a scatter of $0.4$ mag, in our difference images. Given the two day delay, we search this area for a rapidly declining optical counterpart with $\gtrsim 3\sigma$ significance steady decline between the first and final observations. We recover four sources that pass our selection criteria, of which three are cataloged AGN. The fourth source is offset by $5.8$ arcsec from the center of a galaxy at a distance of 187 Mpc, exhibits a rapid decline by $0.5$ mag over $4$ days, and has a red color of $i-z\approx 0.3$ mag. These properties roughly match the expectations for a kilonova. However, this source was detected several times, starting $94$ days prior to GW151226, in the Pan-STARRS Survey for Transients (dubbed as PS15cdi) and is therefore unrelated to the GW event. Given its long-term behavior, PS15cdi is likely a Type IIP supernova that transitioned out of its plateau phase during our observations, mimicking a kilonova-like behavior. We comment on the implications of this detection for contamination in future optical follow-up observations.
[90]  oai:arXiv.org:1608.05456  [pdf] - 1459468
Instrumental performance and results from testing of the BLAST-TNG receiver, submillimeter optics, and MKID arrays
Comments: Presented at SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII, June 29th, 2016
Submitted: 2016-08-18
Polarized thermal emission from interstellar dust grains can be used to map magnetic fields in star forming molecular clouds and the diffuse interstellar medium (ISM). The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) flew from Antarctica in 2010 and 2012 and produced degree-scale polarization maps of several nearby molecular clouds with arcminute resolution. The success of BLASTPol has motivated a next-generation instrument, BLAST-TNG, which will use more than 3000 linear polarization sensitive microwave kinetic inductance detectors (MKIDs) combined with a 2.5m diameter carbon fiber primary mirror to make diffraction-limited observations at 250, 350, and 500 $\mu$m. With 16 times the mapping speed of BLASTPol, sub-arcminute resolution, and a longer flight time, BLAST-TNG will be able to examine nearby molecular clouds and the diffuse galactic dust polarization spectrum in unprecedented detail. The 250 $\mu$m detector array has been integrated into the new cryogenic receiver, and is undergoing testing to establish the optical and polarization characteristics of the instrument. BLAST-TNG will demonstrate the effectiveness of kilo-pixel MKID arrays for applications in submillimeter astronomy. BLAST-TNG is scheduled to fly from Antarctica in December 2017 for 28 days and will be the first balloon-borne telescope to offer a quarter of the flight for "shared risk" observing by the community.
[91]  oai:arXiv.org:1608.04390  [pdf] - 1531273
Variable and polarized radio emission from the T6 brown dwarf WISEP J112254.73+255021.5
Comments: 9 pages, 5 figures, AAS submitted
Submitted: 2016-08-15
Route & Wolszczan (2016) recently detected five radio bursts from the T6 dwarf WISEP J112254.73+255021.5 and used the timing of these events to propose that this object rotates with an ultra-short period of ~17.3 minutes. We conducted follow-up observations with the Very Large Array and Gemini-North but found no evidence for this periodicity. We do, however, observe variable, highly circularly polarized radio emission possibly with a period of 116 minutes, although our observation lasted only 162 minutes and so more data are needed to confirm it. Our proposed periodicity is typical of other radio-active ultracool dwarfs. The handedness of the circular polarization alternates with time and there is no evidence for any unpolarized emission component, the first time such a phenomenology has been observed in radio studies of very low-mass stars and brown dwarfs. We suggest that the object's magnetic dipole axis may be highly misaligned relative to its rotation axis.
[92]  oai:arXiv.org:1608.00978  [pdf] - 1531171
M Dwarf Activity in the Pan-STARRS 1 Medium-Deep Survey: First Catalog and Rotation Periods
Comments: 20 pages, 11 figures. Submitted to ApJ
Submitted: 2016-08-02
We report on an ongoing project to investigate activity in the M dwarf stellar population observed by the Pan-STARRS 1 Medium Deep Survey (PS1-MDS). Using a custom-built pipeline, we refine an initial sample of $\approx$ 4 million sources in PS1-MDS to a sample of 184,148 candidate cool stars using color cuts. Motivated by the well-known relationship between rotation and stellar activity, we use a multi-band periodogram analysis and visual vetting to identify 271 sources that are likely rotating M dwarfs. We derive a new set of polynomials relating M dwarf PS1 colors to fundamental stellar parameters and use them to estimate the masses, distances, effective temperatures, and bolometric luminosities of our sample. We present a catalog containing these values, our measured rotation periods, and cross-matches to other surveys. Our final sample spans periods of $\lesssim$1-130 days in stars with estimated effective temperatures of $\approx$ 2700-4000 K. Twenty-two of our sources have X-ray cross-matches, and they are found to be relatively X-ray bright as would be expected from selection effects. Our data set provides evidence that Kepler-based searches have not been sensitive to very slowly-rotating stars ($P_{\rm rot} \gtrsim 70$ d), implying that the observed emergence of very slow rotators in studies of low-mass stars may be a systematic effect. We also see a lack of low-amplitude ($<$2%) variability in objects with intermediate (10-40 d) rotation periods, which, considered in conjunction with other observational results, may be a signpost of a loss of magnetic complexity associated with a phase of rapid spin-down in intermediate-age M dwarfs. This work represents just a first step in exploring stellar variability in data from the PS1-MDS and, in the farther future, LSST.
[93]  oai:arXiv.org:1607.00943  [pdf] - 1501979
WISEP J060738.65+242953.4: A Nearby. Pole-On L8 Brown Dwarf with Radio Emission
Comments: Accepted to A.J
Submitted: 2016-07-04
We present a simultaneous, multi-wavelength campaign targeting the nearby (7.2 pc) L8/L9 (optical/near-infrared) dwarf WISEP J060738.65+242953.4 in the mid-infrared, radio, and optical. Spitzer Space Telescope observations show no variability at the 0.2% level over 10 hours each in the 3.6 and 4.5 micron bands. Kepler K2 monitoring over 36 days in Campaign 0 rules out stable periodic signals in the optical with amplitudes great than 1.5% and periods between 1.5 hours and 2 days. Non-simultaneous Gemini optical spectroscopy detects lithium, constraining this L dwarf to be less than ~2 Gyr old, but no Balmer emission is observed. The low measured projected rotation velocity (v sin i < 6 km/s) and lack of variability are very unusual compared to other brown dwarfs, and we argue that this substellar object is likely viewed pole-on. We detect quiescent (non-bursting) radio emission with the VLA. Amongst radio detected L and T dwarfs, it has the lowest observed L_nu and the lowest v sin i. We discuss the implications of a pole-on detection for various proposed radio emission scenarios.
[94]  oai:arXiv.org:1606.09586  [pdf] - 1475401
The CHARA Array resolves the long-period Wolf-Rayet binaries WR 137 and WR 138
Comments: 11 pages, 4 tables, 7 figures, accepted to MNRAS
Submitted: 2016-06-30
We report on interferometric observations with the CHARA Array of two classical Wolf-Rayet stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H-band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems ($f_{\rm WR, 137} = 0.59\pm0.04$; $f_{\rm WR, 138} = 0.67\pm0.01$), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modeling of each system. Our spectral modeling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edge-on. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane.
[95]  oai:arXiv.org:1606.08873  [pdf] - 1530905
A Reverse Shock in GRB 160509A
Comments: 8 pages, 2 tables, 4 figures. Submitted to ApJ
Submitted: 2016-06-28
We present the second multi-frequency radio detection of a reverse shock in a $\gamma$-ray burst. By combining our extensive radio observations of the Fermi-LAT GRB 160509A at $z = 1.17$ up to $20$ days after the burst with Swift X-ray observations and ground-based optical and near-infrared data, we show that the afterglow emission comprises distinct reverse shock and forward shock contributions: the reverse shock emission dominates in the radio band at $\lesssim10~$days, while the forward shock emission dominates in the X-ray, optical, and near-infrared bands. Through multi-wavelength modeling, we determine a circumburst density of $n_0\approx10^{-3}~$cm$^{-3}$, supporting our previous suggestion that a low-density circumburst environment is conducive to the production of long-lasting reverse shock radiation in the radio band. We infer the presence of a large excess X-ray absorption column, $N_{\rm H} \approx 1.5\times10^{22}~$cm$^{-2}$, and a high rest-frame optical extinction, $A_{\rm V}\approx3.4~$mag. We identify a jet break in the X-ray light curve at $t_{\rm jet}\approx6~$d, and thus derive a jet opening angle of $\theta_{\rm jet}\approx4~$deg, yielding a beaming-corrected kinetic energy and radiated $\gamma$-ray energy of $E_{\rm K}\approx4\times10^{50}~$erg and $E_{\gamma}\approx1.3\times10^{51}~$erg ($1$-$10^4~$keV, rest frame), respectively. Consistency arguments connecting the forward and reverse shocks suggest a deceleration time of $t_{\rm dec} \approx 460~$s$~\approx T_{90}$, a Lorentz factor of $\Gamma(t_{\rm dec})\approx330$, and a reverse shock to forward shock fractional magnetic energy density ratio of $R_{\rm B}\equiv\epsilon_{\rm B,RS}/\epsilon_{\rm B,FS}\approx8$.
[96]  oai:arXiv.org:1602.08434  [pdf] - 1396901
No precise localization for FRB 150418: claimed radio transient is AGN variability
Comments: ApJL in press; 7 pages, 3 figures, 2 tables. v2: resubmitted version; v1 was an early arxiv-only note
Submitted: 2016-02-26, last modified: 2016-03-25
Keane et al. have recently claimed to have obtained the first precise localization for a Fast Radio Burst (FRB) thanks to the identification of a contemporaneous fading slow (~week-timescale) radio transient. They use this localization to pinpoint the FRB to a galaxy at z ~ 0.49 that exhibits no discernable star formation activity. We argue that the transient is not genuine and that the host candidate, WISE J071634.59-190039.2, is instead a radio variable: the available data did not exclude this possibility; a random radio variable consistent with the observations is not unlikely to have a redshift compatible with the FRB dispersion measure; and the proposed transient light curve is better explained as a scintillating steady source, perhaps also showing an active galactic nucleus (AGN) flare, than a synchrotron-emitting blastwave. The radio luminosity of the host candidate implies that it is an AGN and we present new late-time Very Large Array observations showing that the galaxy is indeed variable at a level consistent with the claimed transient. Therefore the claimed precise localization and redshift determination for FRB 150418 cannot be justified.
[97]  oai:arXiv.org:1510.01226  [pdf] - 1374120
Discovery of an outflow from radio observations of the tidal disruption event ASASSN-14li
Comments: 10 pages, 6 figures, 2 tables; published in ApJ Letters. Fixes an error in our density calculation and expands discussion, major conclusions unchanged from previous version
Submitted: 2015-10-05, last modified: 2016-03-07
We report the discovery of transient radio emission from the nearby optically-discovered TDE ASASSN-14li (distance of 90 Mpc), making it the first typical TDE detected in the radio, and unambiguously pointing to the formation of a non-relativistic outflow with a kinetic energy of $\approx (4-10)\times10^{47}$ erg, a velocity of $\approx 12,000-36,000$ km s$^{-1}$, and a mass of $\approx 3\times10^{-5}-7\times10^{-4}$ M$_{\odot}$. We show that the outflow was ejected on 2014 August 11--25, in agreement with an independent estimate of the timing of super-Eddington accretion based on the optical, UV, and X-ray observations, and that the ejected mass corresponds to about $1-10\%$ of the mass accreted in the super-Eddington phase. The temporal evolution of the radio emission also uncovers the circumnuclear density profile, $\rho(R)\propto R^{-2.5}$ on a scale of about 0.01 pc, a scale that cannot be probed via direct measurements even in the nearest SMBHs. Our discovery of radio emission from the nearest well-studied TDE to date, with a radio luminosity lower than all previous limits, indicates that non-relativistic outflows are ubiquitous in TDEs, and that future, more sensitive, radio surveys will uncover similar events.
[98]  oai:arXiv.org:1511.05559  [pdf] - 1325415
The First Millimeter Detection of a Non-Accreting Ultracool Dwarf
Comments: 7 pages, 3 figures, 1 table; ApJ in press
Submitted: 2015-11-17
The well-studied M9 dwarf TVLM 513-46546 is a rapid rotator (P_rot ~ 2 hr) hosting a stable, dipolar magnetic field of ~3 kG surface strength. Here we report its detection with ALMA at 95 GHz at a mean flux density of $56 \pm 12$ uJy, making it the first ultracool dwarf detected in the millimeter band, excluding young, disk-bearing objects. We also report flux density measurements from unpublished archival VLA data and new optical monitoring data from the Liverpool Telescope. The ALMA data are consistent with a power-law radio spectrum that extends continuously between centimeter and millimeter wavelengths. We argue that the emission is due to the synchrotron process, excluding thermal, free-free, and electron cyclotron maser emission as possible sources. During the interval of the ALMA observation that phases with the maximum of the object's optical variability, the flux density is higher at a ~1.8 sigma significance level. These early results show how ALMA opens a new window for studying the magnetic activity of ultracool dwarfs, particularly shedding light on the particle acceleration mechanism operating in their immediate surroundings.
[99]  oai:arXiv.org:1509.08479  [pdf] - 1308321
Suzaku monitoring of the Wolf-Rayet binary WR140 around periastron passage: An approach for quantifying the wind parameters
Comments: 16 pages, 16 figures, accepted for publication in PASJ
Submitted: 2015-09-28
Suzaku observations of the Wolf-Rayet binary WR 140 (WC7pd+O5.5fc) were made at four different times around periastron passage in 2009 January. The spectra changed in shape and flux with the phase. As periastron approached, the column density of the low-energy absorption increased, which indicates that the emission from the wind-wind collision plasma was absorbed by the dense W-R wind. The spectra can be mostly fitted with two different components: a warm component with kT=0.3--0.6 keV and a dominant hot component with kT~3 keV. The emission measure of the dominant, hot component is not inversely proportional to the distance between the two stars. This can be explained by the O star wind colliding before it has reached its terminal velocity, leading to a reduction in its wind momentum flux. At phases closer to periastron, we discovered a cool plasma component in a recombining phase, which is less absorbed. This component may be a relic of the wind-wind collision plasma, which was cooled down by radiation, and may represent a transitional stage in dust formation.
[100]  oai:arXiv.org:1509.07186  [pdf] - 1308311
Kepler Monitoring of an L Dwarf II. Clouds with Multiyear Lifetimes
Comments: Accepted to ApJ
Submitted: 2015-09-23
We present Kepler, Spitzer Space Telescope, Gemini-North, MMT, and Kitt Peak observations of the L1 dwarf WISEP J190648.47+401106.8. We find that the Kepler optical light curve is consistent in phase and amplitude over the nearly two years of monitoring with a peak-to-peak amplitude of 1.4%. Spitzer Infrared Array Camera 3.6 micron observations are in phase with Kepler with similar light curve shape and peak-to-peak amplitude 1.1%, but at 4.5 micron, the variability has amplitude $<$0.1%. Chromospheric H$\alpha$ emission is variable but not synced with the stable Kepler light curve. A single dark spot can reproduce the light curve but is not a unique solution. An inhomogeneous cloud deck, specifically a region of thick cloud cover, can explain the multi-wavelength data of this ultracool dwarf and need not be coupled with the asynchronous magnetic emission variations. The long life of the cloud is in contrast with weather changes seen in cooler brown dwarfs on the timescale of hours and days.
[101]  oai:arXiv.org:1503.02632  [pdf] - 1262935
Search for Gamma-Ray Emission from DES Dwarf Spheroidal Galaxy Candidates with Fermi-LAT Data
Collaboration, The Fermi-LAT; Collaboration, The DES; :; Drlica-Wagner, A.; Albert, A.; Bechtol, K.; Wood, M.; Strigari, L.; Sanchez-Conde, M.; Baldini, L.; Essig, R.; Cohen-Tanugi, J.; Anderson, B.; Bellazzini, R.; Bloom, E. D.; Caputo, R.; Cecchi, C.; Charles, E.; Chiang, J.; de Angelis, A.; Funk, S.; Fusco, P.; Gargano, F.; Giglietto, N.; Giordano, F.; Guiriec, S.; Gustafsson, M.; Kuss, M.; Loparco, F.; Lubrano, P.; Mirabal, N.; Mizuno, T.; Morselli, A.; Ohsugi, T.; Orlando, E.; Persic, M.; Raino, S.; Sehgal, N.; Spada, F.; Suson, D. J.; Zaharijas, G.; Zimmer, S.; Abbott, T.; Allam, S.; Balbinot, E.; Bauer, A. H.; Benoit-Levy, A.; Bernstein, R. A.; Bernstein, G. M.; Bertin, E.; Brooks, D.; Buckley-Geer, E.; Burke, D. L.; Rosell, A. Carnero; Castander, F. J.; Covarrubias, R.; D'Andrea, C. B.; da Costa, L. N.; DePoy, D. L.; Desai, S.; Diehl, H. T.; Cunha, C. E; Eifler, T. F.; Estrada, J.; Evrard, A. E.; Neto, A. Fausti; Fernandez, E.; Finley, D. A.; Flaugher, B.; Frieman, J.; Gaztanaga, E.; Gerdes, D.; Gruen, D.; Gruendl, R. A.; Gutierrez, G.; Honscheid, K.; Jain, B.; James, D.; Jeltema, T.; Kent, S.; Kron, R.; Kuehn, K.; Kuropatkin, N.; Lahav, O.; Li, T. S.; Luque, E.; Maia, M. A. G.; Makler, M.; March, M.; Marshall, J.; Martini, P.; Merritt, K. W.; Miller, C.; Miquel, R.; Mohr, J.; Neilsen, E.; Nord, B.; Ogando, R.; Peoples, J.; Petravick, D.; Pieres, A.; Plazas, A. A.; Queiroz, A.; Romer, A. K.; Roodman, A.; Rykoff, E. S.; Sako, M.; Sanchez, E.; Santiago, B.; Scarpine, V.; Schubnell, M.; Sevilla, I.; Smith, R. C.; Soares-Santos, M.; Sobreira, F.; Suchyta, E.; Swanson, M. E. C.; Tarle, G.; Thaler, J.; Thomas, D.; Tucker, D.; Walker, A.; Wechsler, R. H.; Wester, W.; Williams, P.; Yanny, B.; Zuntz, J.
Comments: 9 pages, 4 figures, 1 table. Updated to published version. Readers may be interested in the related work by Ackermann et al. arXiv:1503.02641 (The Fermi-LAT Collaboration), Bechtol, Drlica-Wagner, et al. arXiv:1503.02584 (The DES Collaboration), and Koposov, Belokurov, Torrealba, & Evans arXiv:1503.02079
Submitted: 2015-03-09, last modified: 2015-08-16
Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.
[102]  oai:arXiv.org:1508.00724  [pdf] - 1255358
The colliding-wind WC9+OB system WR 65 and dust formation by WR stars
Comments: International Workshop on Wolf-Rayet Stars, Potsdam, Germany, 1 - 5 June 2015
Submitted: 2015-08-04
Observations of the WC9+OB system WR 65 in the infrared show variations of its dust emission consistent with a period near 4.8~yr, suggesting formation in a colliding-wind binary (CWB) having an elliptical orbit. If we adopt the IR maximum as zero phase, the times of X-ray maximum count and minimum extinction to the hard component measured by Oskinova & Hamann fall at phases 0.4--0.5, when the separation of the WC9 and OB stars is greatest. We consider WR 65 in the context of other WC8-9+OB stars showing dust emission.
[103]  oai:arXiv.org:1503.02584  [pdf] - 1269275
Eight New Milky Way Companions Discovered in First-Year Dark Energy Survey Data
Comments: 33 pages, 12 figures, 3 tables. Accepted for publication in ApJ. Readers may be interested in the concurrent work by Koposov, Belokurov, Torrealba, & Evans (http://arxiv.org/abs/1503.02079). Indirect dark matter search results are presented in Drlica-Wagner, Albert, Bechtol, Wood, Strigari, et al. (The LAT and DES Collaborations, http://arxiv.org/abs/1503.02632)
Submitted: 2015-03-09, last modified: 2015-05-28
We report the discovery of eight new Milky Way companions in ~1,800 deg^2 of optical imaging data collected during the first year of the Dark Energy Survey (DES). Each system is identified as a statistically significant over-density of individual stars consistent with the expected isochrone and luminosity function of an old and metal-poor stellar population. The objects span a wide range of absolute magnitudes (M_V from -2.2 mag to -7.4 mag), physical sizes (10 pc to 170 pc), and heliocentric distances (30 kpc to 330 kpc). Based on the low surface brightnesses, large physical sizes, and/or large Galactocentric distances of these objects, several are likely to be new ultra-faint satellite galaxies of the Milky Way and/or Magellanic Clouds. We introduce a likelihood-based algorithm to search for and characterize stellar over-densities, as well as identify stars with high satellite membership probabilities. We also present completeness estimates for detecting ultra-faint galaxies of varying luminosities, sizes, and heliocentric distances in the first-year DES data.
[104]  oai:arXiv.org:1502.06610  [pdf] - 1254867
The Rotation Period and Magnetic Field of the T Dwarf 2MASSI J1047539+212423 Measured From Periodic Radio Bursts
Comments: ApJ submitted; 8 emulateapj pages, 4 figures
Submitted: 2015-02-23
Periodic radio bursts from very low mass stars and brown dwarfs simultaneously probe their magnetic and rotational properties. The brown dwarf 2MASSI J1047539+212423 (2M 1047+21) is currently the only T dwarf (T6.5) detected at radio wavelengths. Previous observations of this source with the Arecibo observatory revealed intermittent, 100%-polarized radio pulses similar to those detected from other brown dwarfs, but were unable to constrain a pulse periodicity; previous VLA observations detected quiescent emission a factor of ~100 times fainter than the Arecibo pulses but no additional events. Here we present 14 hours of Very Large Array observations of this object that reveal a series of pulses at ~6 GHz with highly variable profiles, showing that the pulsing behavior evolves on time scales that are both long and short compared to the rotation period. We measure a periodicity of ~1.77 hr and identify it with the rotation period. This is just the sixth rotation period measurement in a late T dwarf, and the first obtained in the radio. We detect a pulse at 10 GHz as well, suggesting that the magnetic field strength of 2M 1047+21 reaches at least 3.6 kG. Although this object is the coolest and most rapidly-rotating radio-detected brown dwarf to date, its properties appear continuous with those of other such objects, suggesting that the generation of strong magnetic fields and radio emission may continue to even cooler objects. Further studies of this kind will help to clarify the relationships between mass, age, rotation, and magnetic activity at and beyond the end of the main sequence, where both theories and observational data are currently scarce. A video of a short talk discussing this work is at http://youtu.be/OqCtFxcBB2s, and an experimental Web-native version of the paper may be viewed at http://purl.org/net/pkgwpub/goreham.
[105]  oai:arXiv.org:1502.01350  [pdf] - 931480
Extragalactic Transients in the Era of Wide-Field Radio Surveys. I. Detection Rates and Light Curve Characteristics
Comments: 25 pages, 5 figures, 1 appendix, submitted to ApJ
Submitted: 2015-02-04
The impending era of wide-field radio surveys has the potential to revolutionize our understanding of astrophysical transients. Here we evaluate the prospects of a wide range of planned and hypothetical radio surveys using the properties and volumetric rates of known and hypothetical classes of extragalactic synchrotron radio transients (e.g., on- and off-axis gamma-ray bursts [GRB], supernovae, tidal disruption events [TDE], compact object mergers). Utilizing these sources and physically motivated considerations we assess the allowed phase-space of radio luminosity and peak timescale for extragalactic transients. We also include for the first time effects such as redshift evolution of the rates, K-corrections, and non-Euclidean luminosity distance, which affect the detection rates of the most sensitive surveys. The number of detected events is calculated by means of a Monte Carlo method, using the various survey properties (depth, cadence, area) and realistic detection criteria that include a cut on the minimum variability of the transients during the survey and an assessment of host galaxy contamination. Near-term GHz frequency surveys (ASKAP/VAST, Very Large Array Sky Survey) will detect few events: <~30-50 on- and off-axis long GRBs and off-axis tidal disruption events, and ~10-20 neutron star binary mergers if ~1% of the mergers result in a stable millisecond magnetar. Low-frequency surveys (e.g., LOFAR) are unlikely to detect any transients, while a hypothetical large-scale mm survey may detect ~40 on-axis long GRBs. On the other hand, SKA surveys at ~0.1-1 GHz have the potential to uncover thousands of transients, mainly on- and off-axis long GRBs, on-axis short GRBs, off-axis TDEs, and neutron star binary mergers with magnetar remnants.
[106]  oai:arXiv.org:1409.4411  [pdf] - 1216951
Simultaneous Multiwavelength Observations of Magnetic Activity in Ultracool Dwarfs. IV. The Active, Young Binary NLTT 33370 AB (=2MASS J13142039+1320011)
Comments: emulateapj, 22 pages, 15 figures, ApJ in press; v2: fixes low-impact error in Figure 15; v3: now in-press
Submitted: 2014-09-15, last modified: 2014-11-25
We present multi-epoch simultaneous radio, optical, H{\alpha}, UV, and X-ray observations of the active, young, low-mass binary NLTT 33370 AB (blended spectral type M7e). This system is remarkable for its extreme levels of magnetic activity: it is the most radio-luminous ultracool dwarf (UCD) known, and here we show that it is also one of the most X-ray luminous UCDs known. We detect the system in all bands and find a complex phenomenology of both flaring and periodic variability. Analysis of the optical light curve reveals the simultaneous presence of two periodicities, 3.7859 $\pm$ 0.0001 and 3.7130 $\pm$ 0.0002 hr. While these differ by only ~2%, studies of differential rotation in the UCD regime suggest that it cannot be responsible for the two signals. The system's radio emission consists of at least three components: rapid 100% polarized flares, bright emission modulating periodically in phase with the optical emission, and an additional periodic component that appears only in the 2013 observational campaign. We interpret the last of these as a gyrosynchrotron feature associated with large-scale magnetic fields and a cool, equatorial plasma torus. However, the persistent rapid flares at all rotational phases imply that small-scale magnetic loops are also present and reconnect nearly continuously. We present an SED of the blended system spanning more than 9 orders of magnitude in wavelength. The significant magnetism present in NLTT 33370 AB will affect its fundamental parameters, with the components' radii and temperatures potentially altered by ~+20% and ~-10%, respectively. Finally, we suggest spatially resolved observations that could clarify many aspects of this system's nature.
[107]  oai:arXiv.org:1408.6759  [pdf] - 1216562
Eclipses and dust formation by WC9 type Wolf-Rayet stars
Comments: 9 page, 12 figures
Submitted: 2014-08-28
Visual photometry of 16 WC8-9 dust-making Wolf-Rayet (WR) stars during 2001--2009 was extracted from the All Sky Automated Survey All Star Catalogue (ASAS-3) to search for eclipses attributable to extinction by dust formed in clumps in our line of sight. Data for a comparable number of dust-free WC6-9 stars were also examined to help characterise the dataset. Frequent eclipses were observed from WR 104, and several from WR 106, extending the 1994-2001 studies by Kato et al. (2002a,b), but not supporting their phasing the variations in WR 104 with its `pinwheel' rotation period. Only four other stars showed eclipses, WR 50 (one of the dust-free stars), WR 69, WR 95 and WR 117, and there may have been an eclipse by WR 121, which had shown two eclipses in the past. No dust eclipses were shown by the `historic' eclipsers WR 103 and WR 113. The atmospheric eclipses of the latter were observed but the suggestion by David-Uraz et al. (2012) that dust may be partly responsible for these is not supported. Despite its frequent eclipses, there is no evidence in the infrared images of WR 104 for dust made in its eclipses, demonstrating that any dust formed in this process is not a significant contributor to its circumstellar dust cloud and suggesting that the same applies to the other stars showing fewer eclipses.
[108]  oai:arXiv.org:1407.8265  [pdf] - 1216036
Characterization of the Reflectivity of Various Black Materials
Comments:
Submitted: 2014-07-30
We present total and specular reflectance measurements of various materials that are commonly (and uncommonly) used to provide baffling and/or to minimize the effect of stray light in optical systems. More specifically, we investigate the advantage of using certain black surfaces and their role in suppressing stray light on detectors in optical systems. We measure the total reflectance of the samples over a broad wavelength range (250 < lambda < 2500 nm) that is of interest to astronomical instruments in the ultraviolet, visible, and near-infrared regimes. Additionally, we use a helium-neon laser to measure the specular reflectance of the samples at various angles. Finally, we compare these two measurements and derive the specular fraction for each sample.
[109]  oai:arXiv.org:1407.7609  [pdf] - 1215971
Optical BVRI Photometry of Common Proper Motion F/G/K+M Wide Separation Binaries
Comments: 26 pages, 8 figures, accepted for publication in AJ
Submitted: 2014-07-28
We present optical (BVRI) photometric measurements of a sample of 76 common proper motion wide separation main sequence binary pairs. The pairs are composed of a F-, G-, or K-type primary star and an M-type secondary. The sample is selected from the revised NLTT catalog and the LSPM catalog. The photometry is generally precise to 0.03 mag in all bands. We separate our sample into two groups, dwarf candidates and subdwarf candidates, using the reduced proper motion (RPM) diagram constructed with our improved photometry. The M subdwarf candidates in general have larger $V-R$ colors than the M dwarf candidates at a given $V-I$ color. This is consistent with an average metallicity difference between the two groups, as predicted by the PHOENIX/BT-Settl models. The improved photometry will be used as input into a technique to determine the metallicities of the M-type stars.
[110]  oai:arXiv.org:1310.6757  [pdf] - 800703
Trends in Ultracool Dwarf Magnetism. I. X-Ray Suppression and Radio Enhancement
Comments: 24 pages, 12 figures; ApJ in press; v2: accepted version, including new figure
Submitted: 2013-10-24, last modified: 2014-02-20
Although ultracool dwarfs (UCDs) are now known to generate and dissipate strong magnetic fields, a clear understanding of the underlying dynamo is still lacking. We have performed X-ray and radio observations of seven UCDs in a narrow range of spectral type (M6.5-M9.5) but spanning a wide range of projected rotational velocities (v sin i ~ 3-40 km/s). We have also analyzed unpublished archival Chandra observations of four additional objects. All of the newly-observed targets are detected in the X-ray, while only one is detected in the radio, with the remainder having sensitive upper limits. We present a database of UCDs with both radio and X-ray measurements and consider the data in light of the so-called G\"udel-Benz relation (GBR) between magnetic activity in these bands. Some UCDs have very bright radio emission and faint X-ray emission compared to what would be expected for rapid rotators, while others show opposite behavior. We show that UCDs would still be radio-over-luminous relative to the GBR even if their X-ray emission were at standard rapid-rotator "saturation" levels. Recent results from Zeeman-Doppler imaging and geodynamo simulations suggest that rapidly-rotating UCDs may harbor a bistable dynamo that supports either a stronger, axisymmetric magnetic field or a weaker, non-axisymmetric field. We suggest that the data can be explained in a scenario in which strong-field objects obey the GBR while weak-field objects are radio-over-luminous and X-ray-under-luminous, possibly because of a population of gyrosynchrotron-emitting coronal electrons that is continuously replenished by low-energy reconnection events.
[111]  oai:arXiv.org:1310.6758  [pdf] - 800704
Trends in Ultracool Dwarf Magnetism. II. The Inverse Correlation Between X-ray Activity and Rotation as Evidence for a Bimodal Dynamo
Comments: 11 pages, 5 figures, 3 tables; ApJ in press
Submitted: 2013-10-24, last modified: 2014-02-20
Observations of magnetic activity indicators in solar-type stars exhibit a relationship with rotation with an increase until a "saturation" level and a moderate decrease in activity in the very fastest rotators ("supersaturation"). While X-ray data have suggested that this relationship is strongly violated in ultracool dwarfs (UCDs; spectral type >M7), the limited number of X-ray detections has prevented firm conclusions. In this paper, we analyze the X-ray activity-rotation relation in 38 ultracool dwarfs. Our sample represents the largest catalog of X-ray active ultracool dwarfs to date, including seven new and four previously-unpublished Chandra observations presented in a companion paper. We identify a substantial number of rapidly-rotating UCDs with X-ray activity extending two orders of magnitude below the expected saturation level and measure a "supersaturation"-type anticorrelation between rotation and X-ray activity. The scatter in UCD X-ray activity at a fixed rotation is ~3 times larger than that in earlier-type stars. We discuss several mechanisms that have been proposed to explain the data, including centrifugal stripping of the corona, and find them to be inconsistent with the observed trends. Instead, we suggest that an additional parameter correlated with both X-ray activity and rotation is responsible for the observed effects. Building on the results of Zeeman-Doppler imaging of UCD magnetic fields and our companion study of radio/X-ray flux ratios, we argue that this parameter is the magnetic field topology, and that the large scatter in UCD X-ray fluxes reflects the presence of two dynamo modes that produce distinct topologies.
[112]  oai:arXiv.org:1401.1221  [pdf] - 767190
Radio Supernovae in the Local Universe
Comments: Submitted in response to NRAO's call for white papers on the potential VLA Sky Survey. 10 pages, 2 figures
Submitted: 2014-01-06
In the last three decades, about 50 radio supernovae have been detected as a result of targeted searches of optically discovered supernovae in the local universe. Despite this relatively small number some diversity among them has already been identified which is an indication of the underlying richness of radio supernovae waiting to be discovered. For example, comparison of star formation and supernova discovery rate imply that as many as half of the supernovae remain undetected in the traditional optical searches, either because of intrinsic dimness or due to dust obscuration. This has far reaching consequences to the models of stellar and galaxy evolution. A radio sky survey would be ideal to uncover larger supernova population. Transient radio sky would benefit significantly from such a survey. With the advent of advanced gravitational wave detectors a new window is set to open on the local Universe. Localization of these gravitational detectors is poor to identify electromagnetic counterparts of the gravitational wave sources. However, the longer lasting radio emission accompanied in these sources could be effectively identified in a radio sky survey. We advocate a medium area ($\approx$ a few thousand $\rm deg^{2}$) radio sky survey at C-band. Alternatively, a survey at S-band has advantage of larger sky coverage without serious loss of the science case as presented here. Understanding the background in radio sky will be of paramount importance for the upcoming sensitive radio facilities including the Square Kilometer Array.
[113]  oai:arXiv.org:1310.5940  [pdf] - 1180139
Kepler Monitoring of an L Dwarf I. The Photometric Period and White Light Flares
Comments: accepted to ApJ
Submitted: 2013-10-22
We report on the results of fifteen months of monitoring the nearby field L1 dwarf WISEP J190648.47+401106.8 (W1906+40) with the Kepler mission. Supporting observations with the Karl G. Jansky Very Large Array and Gemini North telescope reveal that the L dwarf is magnetically active, with quiescent radio and variable H$\alpha$ emission. A preliminary trigonometric parallax shows that W1906+40 is at a distance of $16.35^{+0.36}_{-0.34}$ pc, and all observations are consistent with W1906+40 being an old disk star just above the hydrogen-burning limit. The star shows photometric variability with a period of 8.9 hours and an amplitude of 1.5%, with a consistent phase throughout the year. We infer a radius of $0.92 \pm 0.07 R_J$ and $\sin i > 0.57$ from the observed period, luminosity ($10^{-3.67 \pm 0.03} L_\odot$), effective temperature ($2300 \pm 75$K), and $v \sin i$ ($11.2 \pm 2.2$ km/s). The light curve may be modeled with a single large, high latitude dark spot. Unlike many L-type brown dwarfs, there is no evidence of other variations at the $\gtrsim 2%$ level, either non-periodic or transient periodic, that mask the underlying rotation period. We suggest that the long-lived surface features may be due to starspots, but the possibility of cloud variations cannot be ruled out without further multi-wavelength observations. During the Gemini spectroscopy, we observed the most powerful flare ever seen on an L dwarf, with an estimated energy of $\sim 1.6 \times 10^{32}$ ergs in white light emission. Using the Kepler data, we identify similar flares and estimate that white light flares with optical/ultraviolet energies of $10^{31}$ ergs or more occur on W1906+40 as often as 1-2 times per month.
[114]  oai:arXiv.org:1301.2321  [pdf] - 648360
Quasi-Quiescent Radio Emission from the First Radio-Emitting T Dwarf
Comments: 7 pages, 3 figures; ApJ Letters in press
Submitted: 2013-01-10, last modified: 2013-03-13
Radio detections of ultracool dwarfs provide insight into their magnetic fields and the dynamos that maintain them, especially at the very bottom of the main sequence, where other activity indicators dramatically weaken. Until recently, radio emission was only detected in the M and L dwarf regimes, but this has changed with the Arecibo detection of rapid polarized flares from the T6.5 dwarf 2MASS J10475385+2124234. Here, we report the detection of quasi-quiescent radio emission from this source at 5.8 GHz using the Karl G. Jansky Very Large Array. The luminosity is {\nu}L_{\nu} = (1.3 +- 0.4) x 10^22 erg/s, a factor of ~100 times fainter than the Arecibo flares. Our detection is the lowest-luminosity yet achieved for an ultracool dwarf. Although the emission is fully consistent with being steady, unpolarized, and broadband, we find tantalizing hints for variability. We exclude the presence of short-duration flares as seen by Arecibo, although this is not unexpected given estimates of the duty cycle. Follow-up observations of this object will offer the potential to constrain its rotation period, electron density, and the strength and configuration of the magnetic field. Equally important, follow-up will address the question of whether the electron cyclotron maser instability, which is thought to produce the flares seen by Arecibo, also operates in the very different parameter regime of the emission we detect, or whether instead this ultracool dwarf exhibits both maser and gyrosynchrotron radiation, potentially originating from substantially different locations.
[115]  oai:arXiv.org:1302.5720  [pdf] - 630672
DECal: A Spectrophotometric Calibration System For DECam
Comments: 10 pages, 5 figures. To appear in conference proceedings: "Calibration and Standardization of Large Surveys and Missions in Astronomy and Astrophysics" (Fermilab, April 16-19, 2012)
Submitted: 2013-02-22
DECal is a new calibration system for the CTIO 4 m Blanco telescope. It is currently being installed as part of the Dark Energy Survey and will provide both broadband flat fields and narrowband (about 1 nm bandwidth) spectrophotometric calibration for the new Dark Energy Camera (DECam). Both of these systems share a new Lambertian flat field screen. The broadband flat field system uses LEDs to illuminate each photometric filter. The spectrophotometric calibration system consists of a monochromator-based tunable light source that is projected onto the flat field screen using a custom line-to-spot fiber bundle and an engineered diffuser. Several calibrated photodiodes positioned along the beam monitor the telescope throughput as a function of wavelength. This system will measure the wavelength-dependent instrumental response function of the total telescope+instrument system in the range 300 < lambda < 1100nm. The spectrophotometric calibration will be performed regularly (roughly once per month) to determine the spectral response of the DECam system and to monitor changes in instrumental throughput during the five year Dark Energy Survey project.
[116]  oai:arXiv.org:1302.2002  [pdf] - 1164487
Variable dust formation by the colliding-wind Wolf-Rayet system HD 36402 in the Large Magellanic Cloud
Comments: 8 pages, 4 figures, to be published in MNRAS
Submitted: 2013-02-08
Infrared photometry of the probable triple WC4(+O?)+O8I: Wolf-Rayet system HD 36402 (= BAT99-38) in the Large Magellanic Cloud (LMC) shows emission characteristic of heated dust. The dust emission is variable on a time-scale of years, with a period near 4.7 yr, possibly associated with orbital motion of the O8 supergiant and the inner P ~ 3.03-d WC4+O binary. The phase of maximum dust emission is close to that of the X-ray minimum, consistent with both processes being tied to colliding wind effects in an elliptical binary orbit. It is evident that Wolf-Rayet dust formation occurs also in metal-poor environments.
[117]  oai:arXiv.org:1211.1042  [pdf] - 606331
ASGARD: A Large Survey for Slow Galactic Radio Transients. I. Overview and First Results
Comments: 35 pages, 18 figures, ApJ in press (v2: accepted version)
Submitted: 2012-11-05, last modified: 2012-11-13
Searches for slow radio transients and variables have generally focused on extragalactic populations, and the basic parameters of Galactic populations remain poorly characterized. We present a large 3 GHz survey performed with the Allen Telescope Array (ATA) that aims to improve this situation: ASGARD, the ATA Survey of Galactic Radio Dynamism. ASGARD observations spanned 2 years with weekly visits to 23 deg^2 in two fields in the Galactic Plane, totaling 900 hr of integration time on science fields and making it significantly larger than previous efforts. The typical blind unresolved source detection limit was 10 mJy. We describe the observations and data analysis techniques in detail, demonstrating our ability to create accurate wide-field images while effectively modeling and subtracting large-scale radio emission, allowing standard transient-and-variability analysis techniques to be used. We present early results from the analysis of two pointings: one centered on the microquasar Cygnus X-3 and one overlapping the Kepler field of view (l = 76{\deg}, b = +13.5{\deg}). Our results include images, catalog statistics, completeness functions, variability measurements, and a transient search. Out of 134 sources detected in these pointings, the only compellingly variable one is Cygnus X-3, and no transients are detected. We estimate number counts for potential Galactic radio transients and compare our current limits to previous work and our projection for the fully-analyzed ASGARD dataset.
[118]  oai:arXiv.org:1210.8246  [pdf] - 1678414
Primary Beam and Dish Surface Characterization at the Allen Telescope Array by Radio Holography
Comments: 19 pages, 23 figures, 3 tables, Authors indicated by an double dagger ({\ddag}) are affiliated with the SETI Institute, Mountain View, CA 95070. Authors indicated by a section break ({\S}) are affiliated with the Hat Creek Radio Observatory and/or the Radio Astronomy Laboratory, both affiliated with the University of California Berkeley, Berkeley CA
Submitted: 2012-10-31
The Allen Telescope Array (ATA) is a cm-wave interferometer in California, comprising 42 antenna elements with 6-m diameter dishes. We characterize the antenna optical accuracy using two-antenna interferometry and radio holography. The distortion of each telescope relative to the average is small, with RMS differences of 1 percent of beam peak value. Holography provides images of dish illumination pattern, allowing characterization of as-built mirror surfaces. The ATA dishes can experience mm-scale distortions across -2 meter lengths due to mounting stresses or solar radiation. Experimental RMS errors are 0.7 mm at night and 3 mm under worst case solar illumination. For frequencies 4, 10, and 15 GHz, the nighttime values indicate sensitivity losses of 1, 10 and 20 percent, respectively. The ATA.s exceptional wide-bandwidth permits observations over a continuous range 0.5 to 11.2 GHz, and future retrofits may increase this range to 15 GHz. Beam patterns show a slowly varying focus frequency dependence. We probe the antenna optical gain and beam pattern stability as a function of focus and observation frequency, concluding that ATA can produce high fidelity images over a decade of simultaneous observation frequencies. In the day, the antenna sensitivity and pointing accuracy are affected. We find that at frequencies greater than 5 GHz, daytime observations greater than 5 GHz will suffer some sensitivity loss and it may be necessary to make antenna pointing corrections on a 1 to 2 hourly basis.
[119]  oai:arXiv.org:1207.1528  [pdf] - 1124630
VAST: An ASKAP Survey for Variables and Slow Transients
Comments: 29 pages, 8 figures. Submitted for publication in Pub. Astron. Soc. Australia
Submitted: 2012-07-06
The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities of ASKAP to enable the discovery and investigation of variable and transient phenomena from the local to the cosmological, including flare stars, intermittent pulsars, X-ray binaries, magnetars, extreme scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. In addition, it will allow us to probe unexplored regions of parameter space where new classes of transient sources may be detected. In this paper we review the known radio transient and variable populations and the current results from blind radio surveys. We outline a comprehensive program based on a multi-tiered survey strategy to characterise the radio transient sky through detection and monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and greater. We also present an analysis of the expected source populations that we will be able to detect with VAST.
[120]  oai:arXiv.org:1203.0330  [pdf] - 530976
Rapid Development of Interferometric Software Using MIRIAD and Python
Comments: Submitted to PASP; 13 pages, 4 figures
Submitted: 2012-03-01
New and upgraded radio interferometers produce data at massive rates and will require significant improvements in analysis techniques to reach their promised levels of performance in a routine manner. Until these techniques are fully developed, productivity and accessibility in scientific programming environments will be key bottlenecks in the pipeline leading from data-taking to research results. We present an open-source software package, miriad-python, that allows access to the MIRIAD interferometric reduction system in the Python programming language. The modular design of MIRIAD and the high productivity and accessibility of Python provide an excellent foundation for rapid development of interferometric software. Several other projects with similar goals exist and we describe them and compare miriad-python to them in detail. Along with an overview of the package design, we present sample code and applications, including the detection of millisecond astrophysical transients, determination and application of nonstandard calibration parameters, interactive data visualization, and a reduction pipeline using a directed acyclic graph dependency model analogous to that of the traditional Unix tool "make". The key aspects of the miriad-python software project are documented. We find that miriad-python provides an extremely effective environment for prototyping new interferometric software, though certain existing packages provide far more infrastructure for some applications. While equivalent software written in compiled languages can be much faster than Python, there are many situations in which execution time is profitably exchanged for speed of development, code readability, accessibility to nonexpert programmers, quick interlinking with foreign software packages, and other virtues of the Python language.
[121]  oai:arXiv.org:1201.5413  [pdf] - 469090
The Allen Telescope Array Commensal Observing System
Comments: 22 pages, 1 figure
Submitted: 2012-01-25
This memo describes the system used to conduct commensal correlator and beamformer observations at the Allen Telescope Array (ATA). This system was deployed for ~2 years until the ATA hibernation in 2011 and was responsible for collecting >5 TB of data during thousands of hours of observations. The general system design is presented and the implementation is discussed in detail. I emphasize the rationale for various design decisions and attempt to document a few aspects of ATA operations that might not be obvious to non-insiders. I close with some recommendations from my experience developing the software infrastructure and managing the correlator observations. These include: reuse existing systems; solve, don't avoid, tensions between projects, and share infrastructure; plan to make standalone observations to complement the commensal ones; and be considerate of observatory staff when deploying new and unusual observing modes. The structure of the software codebase is documented.
[122]  oai:arXiv.org:1111.5194  [pdf] - 1091835
Recurrent dust formation by WR 48a on a 30-year timescale
Comments: 15 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2011-11-22
We present infrared photometry of the WC8 Wolf-Rayet system WR 48a observed with telescopes at ESO, the SAAO and the AAT between 1982 and 2011 which show a slow decline in dust emission from the previously reported outburst in 1978--79 until about 1997, when significant dust emission was still evident. This was followed by a slow rise, accelerating to reach and overtake the first (1978) photometry, demonstrating that the outburst observed in 1978--79 was not an isolated event, but that they recur at intervals of 32+ years. This suggests that WR 48a is a long-period dust maker and colliding-wind binary (CWB). The locus of WR 48a in the (H-L), K colour-magnitude diagram implies that the rate of dust formation fell between 1979 and about 1997 and then increased steadily until 2011. Superimposed on the long-term variation are secondary (`mini') eruptions in (at least) 1990, 1994, 1997, 1999 and 2004, characteristic of relatively brief episodes of additional dust formation. Spectra show evidence for an Oe or Be companion to the WC8 star, supporting the suggestion that WR 48a is a binary system and indicating a system luminosity consistent with the association of WR 48a and the young star clusters Danks 1 and Danks 2. The range of dust formation suggests that these stars are in an elliptical orbit having e ~ 0.6. The size of the orbit implied by the minimum period, together with the WC wind velocity and likely mass-loss rate, implies that the post-shock WC wind is adiabatic throughout the orbit -- at odds with the observed dust formation. A similar conflict is observed in the `pinwheel' dust-maker WR 112.
[123]  oai:arXiv.org:1111.1266  [pdf] - 1091428
First visual orbit for the prototypical colliding-wind binary WR 140
Comments: Complete OIFITS dataset included via Data Conservancy Project
Submitted: 2011-11-04
Wolf-Rayet stars represent one of the final stages of massive stellar evolution. Relatively little is known about this short-lived phase and we currently lack reliable mass, distance, and binarity determinations for a representative sample. Here we report the first visual orbit for WR 140(=HD193793), a WC7+O5 binary system known for its periodic dust production episodes triggered by intense colliding winds near periastron passage. The IOTA and CHARA interferometers resolved the pair of stars in each year from 2003--2009, covering most of the highly-eccentric, 7.9 year orbit. Combining our results with the recent improved double-line spectroscopic orbit of Fahed et al. (2011), we find the WR 140 system is located at a distance of 1.67 +/- 0.03 kpc, composed of a WR star with M_WR = 14.9 +/- 0.5 Msun and an O star with M_O = 35.9 +/- 1.3 Msun. Our precision orbit yields key parameters with uncertainties times 6 smaller than previous work and paves the way for detailed modeling of the system. Our newly measured flux ratios at the near-infrared H and Ks bands allow an SED decomposition and analysis of the component evolutionary states.
[124]  oai:arXiv.org:1109.1014  [pdf] - 406945
Astronomy Career Profiles from the AAS Newsletter Archives
Comments: Original articles edited by Travis Metcalfe, Tony Beasley and Liam McDaid
Submitted: 2011-09-05
This is a collection of articles that were originally published in the Newsletter of the American Astronomical Society (AAS) between May 2008 and September 2011 by the Committee on Employment. Authors representing a wide range of career paths tell their stories and provide insight and advice that is relevant to success in various job sectors. Although all of these articles are available individually from the AAS archives, we are posting the complete collection here to make them more accessible as a resource for the astronomy community. The collection includes the following articles: (1) Changing Priorities: the Hard Money Wild Card, (2) Beyond Ivory Towers, (3) Astronomers Working in Public Outreach, (4) Bush-Whacking a Career Trail, (5) Science Communication as a Press Officer, (6) Jobs in Industry, (7) Back to School: A Ph.D. Enters the Classroom, (8) Working at a Soft-Money Institute, (9) Balancing Research and Service at NOAO, (10) Succeeding in a Large Research Collaboration, and (11) Preparing for the College Teaching Job Market. The final published versions of the articles can be found at http://aas.org/career/
[125]  oai:arXiv.org:1106.4876  [pdf] - 1077519
Millisecond Imaging of Radio Transients with the Pocket Correlator
Comments: 13 pages, accepted to ApJ
Submitted: 2011-06-23, last modified: 2011-08-10
We demonstrate a signal processing concept for imaging the sky at millisecond rates with radio interferometers. The "Pocket Correlator" (PoCo) correlates the signals from multiple elements of a radio interferometer fast enough to image brief, dispersed pulses. By the nature of interferometry, a millisecond correlator functions like a large, single-dish telescope, but with improved survey speed, spatial localization, calibration, and interference rejection. To test the concept, we installed PoCo at the Allen Telescope Array (ATA) to search for dispersed pulses from the Crab pulsar, B0329+54, and M31 using total-power, visibility-based, and image-plane techniques. In 1.7 hours of observing, PoCo detected 191 giant pulses from the Crab pulsar brighter than a typical 5 sigma sensitivity limit of 60 Jy over pulse widths of 3 milliseconds. Roughly 40% of pulses from pulsar B0329+54 were detected by using novel visibility-based techniques. Observations of M31 constrain the rate of pulses brighter than 190 Jy in a three degree region surrounding the galaxy to <4.3/hr. We calculate the computational demand of various visibility-based pulse search algorithms and demonstrate how compute clusters can help meet this demand. Larger implementations of the fast imaging concept will conduct blind searches for millisecond pulses in our Galaxy and beyond, providing a valuable probe of the interstellar/intergalactic media, discovering new kinds of radio transients, and localizing them to constrain models of their origin.
[126]  oai:arXiv.org:1107.1517  [pdf] - 1077843
The Allen Telescope Array Pi GHz Sky Survey II. Daily and Monthly Monitoring for Transients and Variability in the Bootes Field
Comments: Accepted for publication in ApJ
Submitted: 2011-07-07
We present results from daily radio continuum observations of the Bootes field as part of the Pi GHz Sky Survey (PiGSS). These results are part of a systematic and unbiased campaign to characterize variable and transient sources in the radio sky. The observations include 78 individual epochs distributed over 5 months at a radio frequency of 3.1 GHz with a median RMS image noise in each epoch of 2.8 mJy. We produce 5 monthly images with a median RMS of 0.6 mJy. No transient radio sources are detected in the daily or monthly images. At 15 mJy, we set an upper limit (2 sigma) to the surface density of 1-day radio transients at 0.025 deg^-2. At 5 mJy, we set an upper limit (2\sigma) to the surface density of 1-month radio transients at 0.18 deg^-2. We also produce light curves for 425 sources and explore the variability properties of these sources. Approximately 20% of the sources exhibit some variability on daily and monthly time scales. The maximum RMS fractional modulations on the 1 day and 1 month time scales for sources brighter than 10 mJy are 2 and 0.5, respectively. The probability of a daily fluctuation for all sources and all epochs by a factor of 10 is less than 10^-4. We compare the radio to mid-infrared variability for sources in the field and find no correlation. Finally, we apply the statistics of transient and variable populations to constrain models for a variety of source classes.
[127]  oai:arXiv.org:1104.5248  [pdf] - 355598
The 2010 May Flaring Episode of Cygnus X-3 in Radio, X-Rays, and {\gamma}-Rays
Comments: 6 pages, 3 figures, accepted ApJL
Submitted: 2011-04-27
In 2009, Cygnus X-3 (Cyg X-3) became the first microquasar to be detected in the GeV {\gamma}-ray regime, via the satellites Fermi and AGILE. The addition of this new band to the observational toolbox holds promise for building a more detailed understanding of the relativistic jets of this and other systems. We present a rich dataset of radio, hard and soft X-ray, and {\gamma}-ray observations of Cyg X-3 made during a flaring episode in 2010 May. We detect a ~3-d softening and recovery of the X-ray emission, followed almost immediately by a ~1-Jy radio flare at 15 GHz, followed by a 4.3{\sigma} {\gamma}-ray flare (E > 100 MeV) ~1.5 d later. The radio sampling is sparse, but we use archival data to argue that it is unlikely the {\gamma}-ray flare was followed by any significant unobserved radio flares. In this case, the sequencing of the observed events is difficult to explain in a model in which the {\gamma}-ray emission is due to inverse Compton scattering of the companion star's radiation field. Our observations suggest that other mechanisms may also be responsible for {\gamma}-ray emission from Cyg X-3.
[128]  oai:arXiv.org:1103.1696  [pdf] - 1052588
Comparisons of Supergranule Characteristics During the Solar Minima of Cycles 22/23 and 23/24
Comments: 22 pages, 5 figures. Solar Physics, in press
Submitted: 2011-03-08
Supergranulation is a component of solar convection that manifests itself on the photosphere as a cellular network of around 35 Mm across, with a turnover lifetime of 1-2 days. It is strongly linked to the structure of the magnetic field. The horizontal, divergent flows within supergranule cells carry local field lines to the cell boundaries, while the rotational properties of supergranule upflows may contribute to the restoration of the poloidal field as part of the dynamo mechanism that controls the solar cycle. The solar minimum at the transition from cycle 23 to 24 was notable for its low level of activity and its extended length. It is of interest to study whether the convective phenomena that influences the solar magnetic field during this time differed in character to periods of previous minima. This study investigates three characteristics (velocity components, sizes and lifetimes) of solar supergranulation. Comparisons of these characteristics are made between the minima of cycles 22/23 and 23/24 using MDI Doppler data from 1996 and 2008, respectively. It is found that whereas the lifetimes are equal during both epochs (around 18 h), the sizes are larger in 1996 (35.9 +/- 0.3 Mm) than in 2008 (35.0 +/- 0.3 Mm), while the dominant horizontal velocity flows are weaker (139 +/- 1 m/s in 1996; 141 +/- 1 m/s in 2008). Although numerical differences are seen, they are not conclusive proof of the most recent minimum being inherently unusual.
[129]  oai:arXiv.org:1102.2227  [pdf] - 1052030
The Allen Telescope Array Twenty-centimeter Survey -- A 700-Square-Degree, Multi-Epoch Radio Dataset -- II: Individual Epoch Transient Statistics
Comments: 28 pages, 12 figures, ApJ accepted
Submitted: 2011-02-10
We present our second paper on the Allen Telescope Array Twenty-centimeter Survey (ATATS), a multi-epoch, ~700 sq. deg. radio image and catalog at 1.4 GHz. The survey is designed to detect rare, bright transients as well as to commission the ATA's wide-field survey capabilities. ATATS explores the challenges of multi-epoch transient and variable source surveys in the domain of dynamic range limits and changing (u,v) coverage. Here we present images made using data from the individual epochs, as well as a revised image combining data from all ATATS epochs. The combined image has RMS noise 3.96 mJy / beam, with a circular beam of 150 arcsec FWHM. The catalog, generated using a false detection rate algorithm, contains 4984 sources, and is >90% complete to 37.9 mJy. The catalogs generated from snapshot images of the individual epochs contain between 1170 and 2019 sources over the 564 sq. deg. area in common to all epochs. The 90% completeness limits of the single epoch catalogs range from 98.6 to 232 mJy. We compare the catalog generated from the combined image to those from individual epochs, and from the NRAO VLA Sky Survey (NVSS), a legacy survey at the same frequency. We are able to place new constraints on the transient population: fewer than 6e-4 transients / sq. deg., for transients brighter than 350 mJy with characteristic timescales of minutes to days. This strongly rules out an astronomical origin for the ~1 Jy sources reported by Matsumura et al. (2009), based on their stated rate of 3.1e-3 / sq. deg.
[130]  oai:arXiv.org:1101.1430  [pdf] - 955892
The WR 140 periastron passage 2009: first results from MONS and other optical sources
Comments: 5 pages, 5 figures, Published in proceedings of "Stellar Winds in Interaction", editors T. Eversberg and J.H. Knapen. Full proceedings volume is available on http://www.stsci.de/pdf/arrabida.pdf
Submitted: 2011-01-07
We present the results from the spectroscopic follow-up of WR140 (WC7 + O4-5) during its last periastron passage in January 2009. This object is known as the archetype of colliding wind binaries and has a relatively large period (~ 8 years) and eccentricity (~ 0.89). We provide updated values for the orbital parameters, new estimates for the WR and O star masses and new constraints on the mass-loss rates.
[131]  oai:arXiv.org:1101.1046  [pdf] - 287109
WR 140 in the Infrared
Comments: 11 pages, 16 figures, Published in proceedings of "Stellar Winds in Interaction", editors T. Eversberg and J.H. Knapen. Full proceedings volume is available on http://www.stsci.de/pdf/arrabida.pdf
Submitted: 2011-01-05
Observations in the infrared have played a significant role in the development of our understanding of WR140. Two sets of observations are described here: the changing profile of the 1.083-{\mu}m He I line observed for the 2009 campaign to study evolution of the wind-collision region near periastron passage, and multiwavelength photometry and imaging of the dust made by WR140 in previous periastron passages.
[132]  oai:arXiv.org:1012.0945  [pdf] - 1042380
Spectropolarimetry with the Allen Telescope Array: Faraday Rotation toward Bright Polarized Radio Galaxies
Comments: Accepted to ApJ, 83 (preprint format) pages, 49 figures
Submitted: 2010-12-04
We have observed 37 bright, polarized radio sources with the Allen Telescope Array (ATA) to present a novel analysis of their Faraday rotation properties. Each source was observed during the commissioning phase with 2 to 4 100-MHz bands at frequencies ranging from 1 to 2 GHz. These observations demonstrate how the continuous frequency coverage of the ATA's log-periodic receiver can be applied to the study of Faraday rotation measures (RMs). We use RM synthesis to show that wide-bandwidth data can find multiple RM components toward a single source. Roughly a quarter of the sources studied have extra RM components with high confidence (brighter than ~40 mJy), when observing with a RM resolution of roughly 100 rad/m2. These extra components contribute 10%-70% of the total polarized flux. This is the first time multiple RM components have been identified in a large sample of point sources. For our observing configuration, these extra RM components bias the measurement of the peak RM by 10-15 rad/m2 ; more generally, the peak RM cannot be determined more precisely than the RM beam size. Comparing our 1-2 GHz RM spectra to VLBA polarimetric maps shows both techniques can identify complicated Faraday structures in the sources. However, the RM values and fractional polarization are generally smaller at lower frequencies than in the higher-frequency VLBA maps. With a few exceptions, the RMs from this work are consistent with that of earlier, narrow-bandwidth, all-sky surveys. This work also describes the polarimetry calibration procedure and that on-axis ATA observations of linear polarization can be calibrated to an accuracy of 0.2% of Stokes I. Future research directions include studying the time-dependent RM structure in Active Galactic Nuclei (AGNs) and enabling accurate, wide-area RM surveys to test models of Galactic and extragalactic magnetic fields.
[133]  oai:arXiv.org:1011.0777  [pdf] - 253304
Cyg OB2 #5: When three stars are just not enough
Comments: 5 pages, 4 figs. To appear in "The Multi-wavelength view of hot, massive stars - 39th Liege Int. Astroph. Coll., 12-16 July 2010
Submitted: 2010-11-02
Archival observations from the Very Large Array (VLA) at frequencies between 1.4 GHz and 43 GHz of the 6.6-day O6.5-7+O5.5-6 binary Cyg OB2 #5 over 20 years are re-examined. The aim is to determine the location and character of its known variable radio emission. The radio emission consists of a primary component associated with the binary, and a non-thermal source (NE), $0.8" to the NE. This work reveals that NE shows no evidence of variation demonstrating that the variable emission arises in the primary component. With NE constant, the radio flux from the primary can now be well determined for the first time, especially in observations that do not resolve both the primary and NE components. The variable radio emission from the primary has a period of 6.7+/-0.3 years which is described by a simple model of a non-thermal source orbiting within the stellar wind envelope of the binary. Such a model implies the presence of a third, unresolved stellar companion (Star C) orbiting the 6.6-day binary with a period of 6.7 years. The non-thermal emission arises from either a WCR between Star C and the binary system, or possibly from Star C directly. Examination of radial velocity observations suggests reflex motion of the binary due to Star C, for which a mass of 23^+22_-14 solar mass is deduced. Together with the star associated with NE, this implies that Cyg OB2 #5 is a quadruple system.
[134]  oai:arXiv.org:1010.1567  [pdf] - 239199
Commensal observing with the Allen Telescope array: software command and control
Comments: SPIE Conference Proceedings, Software and Cyberinfrastructure for Astronomy, Nicole M. Radziwill; Alan Bridger, Editors, 77400Z, Vol 7740
Submitted: 2010-10-07, last modified: 2010-10-11
The Allen Telescope Array (ATA) is a Large-Number-Small-Diameter radio telescope array currently with 42 individual antennas and 5 independent back-end science systems (2 imaging FX correlators and 3 time domain beam formers) located at the Hat Creek Radio Observatory (HCRO). The goal of the ATA is to run multiple back-ends simultaneously, supporting multiple science projects commensally. The primary software control systems are based on a combination of Java, JRuby and Ruby on Rails. The primary control API is simplified to provide easy integration with new back-end systems while the lower layers of the software stack are handled by a master observing system. Scheduling observations for the ATA is based on finding a union between the science needs of multiple projects and automatically determining an efficient path to operating the various sub-components to meet those needs. When completed, the ATA is expected to be a world-class radio telescope, combining dedicated SETI projects with numerous radio astronomy science projects.
[135]  oai:arXiv.org:1010.1064  [pdf] - 1041094
Primary Beam Shape Calibration from Mosaicked, Interferometric Observations
Comments: Accepted by PASP. 8 pages, 8 figures
Submitted: 2010-10-06
Image quality in mosaicked observations from interferometric radio telescopes is strongly dependent on the accuracy with which the antenna primary beam is calibrated. The next generation of radio telescope arrays such as the Allen Telescope Array (ATA) and the Square Kilometer Array (SKA) have key science goals that involve making large mosaicked observations filled with bright point sources. We present a new method for calibrating the shape of the telescope's mean primary beam that uses the multiple redundant observations of these bright sources in the mosaic. The method has an analytical solution for simple Gaussian beam shapes but can also be applied to more complex beam shapes through $\chi^2$ minimization. One major benefit of this simple, conceptually clean method is that it makes use of the science data for calibration purposes, thus saving telescope time and improving accuracy through simultaneous calibration and observation. We apply the method both to 1.43 GHz data taken during the ATA Twenty Centimeter Survey (ATATS) and to 3.14 GHz data taken during the ATA's Pi Gigahertz Sky Survey (PiGSS). We find that the beam's calculated full width at half maximum (FWHM) values are consistent with the theoretical values, the values measured by several independent methods, and the values from the simulation we use to demonstrate the effectiveness of our method on data from future telescopes such as the expanded ATA and the SKA. These results are preliminary, and can be expanded upon by fitting more complex beam shapes. We also investigate, by way of a simulation, the dependence of the accuracy of the telescope's FWHM on antenna number. We find that the uncertainty returned by our fitting method is inversely proportional to the number of antennas in the array.
[136]  oai:arXiv.org:1009.4443  [pdf] - 1040771
The Allen Telescope Array Pi GHz Sky Survey I. Survey Description and Static Catalog Results for the Bootes Field
Comments: Accepted for publication in ApJ; revision submitted with extraneous figure removed
Submitted: 2010-09-22, last modified: 2010-09-27
The Pi GHz Sky Survey (PiGSS) is a key project of the Allen Telescope Array. PiGSS is a 3.1 GHz survey of radio continuum emission in the extragalactic sky with an emphasis on synoptic observations that measure the static and time-variable properties of the sky. During the 2.5-year campaign, PiGSS will twice observe ~250,000 radio sources in the 10,000 deg^2 region of the sky with b > 30 deg to an rms sensitivity of ~1 mJy. Additionally, sub-regions of the sky will be observed multiple times to characterize variability on time scales of days to years. We present here observations of a 10 deg^2 region in the Bootes constellation overlapping the NOAO Deep Wide Field Survey field. The PiGSS image was constructed from 75 daily observations distributed over a 4-month period and has an rms flux density between 200 and 250 microJy. This represents a deeper image by a factor of 4 to 8 than we will achieve over the entire 10,000 deg^2. We provide flux densities, source sizes, and spectral indices for the 425 sources detected in the image. We identify ~100$ new flat spectrum radio sources; we project that when completed PiGSS will identify 10^4 flat spectrum sources. We identify one source that is a possible transient radio source. This survey provides new limits on faint radio transients and variables with characteristic durations of months.
[137]  oai:arXiv.org:1009.4627  [pdf] - 1040790
Variability monitoring of OB stars during the Mons campaign
Comments: 2 pages, 1 figure. To appear in proceedings of IAUS272 'Active OB Stars: Structure, Evolution, Mass Loss and Critical Limits'
Submitted: 2010-09-23
We present preliminary results of a 3-month campaign carried out in the framework of the Mons project, where time-resolved Halpha observations are used to study the wind and circumstellar properties of a number of OB stars.
[138]  oai:arXiv.org:1008.4385  [pdf] - 1034504
The Advection of Supergranules by the Sun's Axisymmetric Flows
Comments: 15 pages, 8 figures
Submitted: 2010-08-25
We show that the motions of supergranules are consistent with a model in which they are simply advected by the axisymmetric flows in the Sun's surface shear layer. We produce a 10-day series of simulated Doppler images at a 15-minute cadence that reproduces most spatial and temporal characteristics seen in the SOHO/MDI Doppler data. Our simulated data have a spectrum of cellular flows with just two components -- a granule component that peaks at spherical wavenumbers of about 4000 and a supergranule component that peaks at wavenumbers of about 110. We include the advection of these cellular components by the axisymmetric flows -- differential rotation and meridional flow -- whose variations with latitude and depth (wavenumber) are consistent with observations. We mimic the evolution of the cellular pattern by introducing random variations to the phases of the spectral components at rates that reproduce the levels of cross-correlation as functions of time and latitude. Our simulated data do not include any wave-like characteristics for the supergranules yet can reproduce the rotation characteristics previously attributed to wave-like behavior. We find rotation rates which appear faster than the actual rotation rates and attribute this to the projection effects. We find that the measured meridional flow does accurately represent the actual flow and that the observations indicate poleward flow to $65\degr-70\degr$ latitude with equatorward counter cells in the polar regions.
[139]  oai:arXiv.org:1007.3645  [pdf] - 442239
Heated dust around the LMC Wolf-Rayet system HD 36402 (BAT99-38)
Comments: 4 pages, 4 figures, contribution to 39th Liege Astrophysical Colloquium, 'The multi-wavelength view of hot, massive stars'
Submitted: 2010-07-21
Infrared photometry of the probable triple WC4 (+O?) +O8I: system HD 36402 (BAT99-38) in the LMC shows emission characteristic of heated dust. Although HD 36402 is close to two luminous YSOs, it is possible to distinguish its emission at wavelengths less than 10 microns. Simple modelling indicates a dust temperature near 800 K and mass of about 1.5e-7 Msol amorphous carbon grains. The dust emission appears to be variable. It is apparent that Wolf Rayet dust formation occurs also in metal-poor environments.
[140]  oai:arXiv.org:1006.2003  [pdf] - 1033022
The Allen Telescope Array Twenty-centimeter Survey - A 690-Square-Degree, 12-Epoch Radio Dataset - I: Catalog and Long-Duration Transient Statistics
Comments: 41 pages, 19 figures, ApJ accepted; corrected minor typo in Table 1
Submitted: 2010-06-10, last modified: 2010-06-11
We present the Allen Telescope Array Twenty-centimeter Survey (ATATS), a multi-epoch (12 visits), 690 square degree radio image and catalog at 1.4GHz. The survey is designed to detect rare, very bright transients as well as to verify the capabilities of the ATA to form large mosaics. The combined image using data from all 12 ATATS epochs has RMS noise sigma = 3.94mJy / beam and dynamic range 180, with a circular beam of 150 arcsec FWHM. It contains 4408 sources to a limiting sensitivity of S = 20 mJy / beam. We compare the catalog generated from this 12-epoch combined image to the NRAO VLA Sky Survey (NVSS), a legacy survey at the same frequency, and find that we can measure source positions to better than ~20 arcsec. For sources above the ATATS completeness limit, the median flux density is 97% of the median value for matched NVSS sources, indicative of an accurate overall flux calibration. We examine the effects of source confusion due to the effects of differing resolution between ATATS and NVSS on our ability to compare flux densities. We detect no transients at flux densities greater than 40 mJy in comparison with NVSS, and place a 2-sigma upper limit on the transient rate for such sources of 0.004 per square degree. These results suggest that the > 1 Jy transients reported by Matsumura et al. (2009) may not be true transients, but rather variable sources at their flux density threshold.
[141]  oai:arXiv.org:1004.0936  [pdf] - 143193
Photospheric Manifestations of Supergranules During the Last Two Solar Minima
Comments: 4 pages, 2 figures
Submitted: 2010-04-06
Solar supergranulation plays an important role in generating and structuring the solar magnetic field and as a mechanism responsible for the 11-year solar cycle. It is clearly detected within SOHO/MDI Dopplergrams, from which a variety of properties may be derived. Techniques that extract spatial, temporal and kinematic characteristics and provide comparisons for the two most recent solar minima are described. Although supergranule lifetimes are comparable between these minima, their sizes maybe slightly smaller during the recent minimum.
[142]  oai:arXiv.org:1003.4210  [pdf] - 134560
The Advection of Supergranules by Large-Scale Flows
Comments: GONG 2008/SOHO 21 Meeting (Solar-stellar dynamos as revealed by helio- and asteroseismology, Boulder, CO, USA, August 11-15, 2008), 6 pages, 2 figures
Submitted: 2010-03-22
We produce a 10-day series of simulated Doppler images at a 15-minute cadence that reproduces the spatial and temporal characteristics seen in the SOHO/MDI Doppler data. Our simulated data contains a spectrum of cellular flows with but two necessary components --- a granule component that peaks at wavenumbers of about 4000 and a supergranule component that peaks at wavenumbers of about 110. We include the advection of these cellular components by a differential rotation profile that depends on latitude and wavenumber (depth). We further mimic the evolution of the cellular pattern by introducing random variations to the amplitudes and phases of the spectral components at rates that reproduce the level of cross-correlation as a function of time and latitude. Our simulated data do not include any wave-like characteristics for the supergranules yet can accurately reproduce the rotation characteristics previously attributed to wave-like characteristics.
[143]  oai:arXiv.org:0911.5674  [pdf] - 902580
Modelling the radio emission from Cyg OB2 #5: a quadruple system?
Comments: Accepted for publication in ApJ Dec 1, 2009. 13 pages, 11 figures, 4 tables
Submitted: 2009-11-30
Fifty archival radio observations of the supergiant binary Cyg OB2 #5 using the Very Large Array over 20 years are re-examined to determine the location and character of the previously detected variable radio emission. The radio emission from the system consists of a primary component that is associated with the binary, and a non-thermal source (NE) that has been ascribed to a wind-collision region (WCR) between the stellar winds of the binary and that of a B-type star (Star D) to the NE. NE shows no evidence of variation in 23 epochs where it is resolved separately from the primary radio component, demonstrating that the variable emission arises in the primary. Since NE is non-variable, the radio flux from the primary can now be well determined for the first time, most especially in observations that do not resolve both the primary and NE components. The variable radio emission from the primary has a period of 6.7+/-0.3 years which is described by a simple model of a non-thermal source orbiting within the stellar wind envelope of the binary. Such a model implies the presence of a third, unresolved stellar companion (Star C) orbiting the 6.6-day binary with a period of 6.7 years. The variable non-thermal emission arises from either a WCR between Star C and the binary system, or possibly from Star C directly. The model gives a mass-loss rate of 3.4 x 10^{-5} solar mass/yr for Cyg OB2 #5, unusually high for an Of supergiant and comparable to that of WR stars, and consistent with an unusually strong He I 1.083-micron emission line, also redolent of WR stars. An examination of radial velocity observations suggests reflex motion of the binary due to Star C. The natures of NE and Star D are also examined. (abridged)
[144]  oai:arXiv.org:0912.0014  [pdf] - 902587
Evaluating the Calorimeter Model with Broadband, Continuous Spectra of Starburst Galaxies Observed with the Allen Telescope Array
Comments: 44 pages, 15 figures, ApJ accepted
Submitted: 2009-11-30
Although the relationship between the far-infrared and cm-wave radio luminosities of normal galaxies is one of the most striking correlations in astronomy, a solid understanding of its physical basis is lacking. In one interpretation, the "calorimeter model," rapid synchrotron cooling of cosmic ray electrons is essential in reproducing the observed linear relationship. Observed radio spectra, however, are shallower than what is expected of cooled synchrotron emission. In Thompson et al. (2006), a simple parameterized model is presented to explain how relatively shallow observed spectra might arise even in the presence of rapid synchrotron cooling by accounting for ionization losses and other cooling mechanisms. During the commissioning of the 42-element Allen Telescope Array, we observed the starburst galaxies M82, NGC 253, and Arp 220 at frequencies ranging from 1 to 7 GHz, obtaining unprecedented broadband continuous radio spectra of these sources. We combine our observations with high-frequency data from the literature to separate the spectra into thermal and nonthermal components. The nonthermal components all steepen in the cm-wave regime and cannot be well-modeled as simple power laws. The model of Thompson et al. is consistent with our M82 results when plausible parameters are chosen, and our results in fact significantly shrink the space of allowed model parameters. The model is only marginally consistent with our NGC 253 data. Assuming the Thompson et al. model, a steep electron energy injection index of p = -2.5 is ruled out in M82 and NGC 253 to >99% confidence. We describe in detail the observing procedures, calibration methods, analysis, and consistency checks used for broadband spectral observations with the Allen Telescope Array.
[145]  oai:arXiv.org:0905.3073  [pdf] - 1002280
Archiving multi-epoch data and the discovery of variables in the near infrared
Comments: 30 pages, 31 figures, MNRAS, in press Minor changes from previous version due to refereeing and proof-reading
Submitted: 2009-05-19, last modified: 2009-09-16
We present a description of the design and usage of a new synoptic pipeline and database model for time series photometry in the VISTA Data Flow System (VDFS). All UKIRT-WFCAM data and most of the VISTA main survey data will be processed and archived by the VDFS. Much of these data are multi-epoch, useful for finding moving and variable objects. Our new database design allows the users to easily find rare objects of these types amongst the huge volume of data being produced by modern survey telescopes. Its effectiveness is demonstrated through examples using Data Release 5 of the UKIDSS Deep Extragalactic Survey (DXS) and the WFCAM standard star data. The synoptic pipeline provides additional quality control and calibration to these data in the process of generating accurate light-curves. We find that 0.6+-0.1% of stars and 2.3+-0.6% of galaxies in the UKIDSS-DXS with K<15 mag are variable with amplitudes \Delta K>0.015 mag
[146]  oai:arXiv.org:0904.2786  [pdf] - 1001952
Old, Rich, and Eccentric: Two Jovian Planets Orbiting Evolved Metal-Rich Stars
Comments: 7 pages, 3 figures
Submitted: 2009-04-17, last modified: 2009-06-20
We present radial velocity measurements of two stars observed as part of the Lick Subgiants Planet Search and the Keck N2K survey. Variations in the radial velocities of both stars reveal the presence of Jupiter-mass exoplanets in highly eccentric orbits. HD 16175 is a G0 subgiant from the Lick Subgiants Planet Search, orbited by a planet having a minimum mass of 4.4 M_Jup, in an eccentric (e = 0.59), 2.71 yr orbit. HD 96167 is a G5 subgiant from the N2K ("Next 2000") program at Keck Observatory, orbited by a planet having a minimum mass of 0.68 M_Jup, in an eccentric (e = 0.71), 1.366 yr orbit. Both stars are relatively massive (M_star = 1.3 M_sun) and are very metal rich ([Fe/H] > +0.3). We describe our methods for measuring the stars' radial velocity variations and photometric stability.
[147]  oai:arXiv.org:0904.2571  [pdf] - 23402
Training the Next Generation of Astronomers
Comments: 9 pp + cover page; position paper submitted to the Astro2010 Decadal Survey Subcommittee on the State of the Profession; small typographical corrections made; submitted version at http://www8.nationalacademies.org/astro2010/DetailFileDisplay.aspx?id=460
Submitted: 2009-04-16
While both society and astronomy have evolved greatly over the past fifty years, the academic institutions and incentives that shape our field have remained largely stagnant. As a result, the astronomical community is faced with several major challenges, including: (1) the training that we provide does not align with the skills that future astronomers will need, (2) the postdoctoral phase is becoming increasingly demanding and demoralizing, and (3) our jobs are increasingly unfriendly to families with children. Solving these problems will require conscious engineering of our profession. Fortunately, this Decadal Review offers the opportunity to revise outmoded practices to be more effective and equitable. The highest priority of the Subcommittee on the State of the Profession should be to recommend specific, funded activities that will ensure the field meets the challenges we describe.
[148]  oai:arXiv.org:0904.0762  [pdf] - 596261
The Allen Telescope Array: The First Widefield, Panchromatic, Snapshot Radio Camera for Radio Astronomy and SETI
Comments: Special Issue of Proceedings of the IEEE: "Advances in Radio Telescopes", Baars,J. Thompson,R., D'Addario, L., eds, 2009, in press
Submitted: 2009-04-06
The first 42 elements of the Allen Telescope Array (ATA-42) are beginning to deliver data at the Hat Creek Radio Observatory in Northern California. Scientists and engineers are actively exploiting all of the flexibility designed into this innovative instrument for simultaneously conducting surveys of the astrophysical sky and conducting searches for distant technological civilizations. This paper summarizes the design elements of the ATA, the cost savings made possible by the use of COTS components, and the cost/performance trades that eventually enabled this first snapshot radio camera. The fundamental scientific program of this new telescope is varied and exciting; some of the first astronomical results will be discussed.
[149]  oai:arXiv.org:0902.4540  [pdf] - 900633
Orbitally modulated dust formation by the WC7+O5 colliding-wind binary WR140
Comments: 21 pages, 10 figures, accepted for MNRAS. A version with higher resolution figures is available at ftp://ftp.roe.ac.uk/pub/pmw/wr140dust.ps.gz
Submitted: 2009-02-26
We present high-resolution infrared (2--18 micron) images of the archetypal periodic dust-making Wolf-Rayet binary system WR140 (HD 193793) taken between 2001 and 2005, and multi-colour (J -- [19.5]) photometry observed between 1989 and 2001. The images resolve the dust cloud formed by WR140 in 2001, allowing us to track its expansion and cooling, while the photometry allows tracking the average temperature and total mass of the dust. The combination of the two datasets constrains the optical properties of the dust. The most persistent dust features, two concentrations at the ends of a `bar' of emission to the south of the star, were observed to move with constant proper motions of 324+/-8 and 243+/-7 mas/y. Longer wavelength (4.68-micron and 12.5-micron) images shows dust emission from the corresponding features from the previous (1993) periastron passage and dust-formation episode. A third persistent dust concentration to the east of the binary (the `arm') was found to have a proper motion ~ 320 mas/y. Extrapolation of the motions of the concentrations back to the binary suggests that the eastern `arm' began expansion 4--5 months earlier than those in the southern `bar', consistent with the projected rotation of the binary axis and wind-collision region (WCR) on the sky. Comparison of model dust images and the observations constrain the intervals when the WCR was producing sufficiently compressed wind for dust nucleation in the WCR, and suggests that the distribution of this material was not uniform about the axis of the WCR, but more abundant in the following edge in the orbital plane.
[150]  oai:arXiv.org:0902.4105  [pdf] - 900629
Dust formation by the colliding-wind WC5+O9 binary WR19 at periastron passage
Comments: 5 pages and 4 figures, accepted for publication in MNRAS
Submitted: 2009-02-24
We present infrared photometry of the episodic dust-making Wolf-Rayet system WR19 (LS3), tracking its fading from a third observed dust-formation episode in 2007 and strengthening the view that these episodes are periodic (P = 10.1+/-0.1 y). Radial velocities of the O9 component observed between 2001 and 2008 show RV variations consistent with WC19 being a spectroscopic binary of high eccentricity (e=0.8), having periastron passage in 2007.14, shortly before the phase of dust formation. In this respect, WR19 resembles the archetypical episodic dust-making colliding-wind binary system WR140.
[151]  oai:arXiv.org:0806.1734  [pdf] - 13431
Measurement of the Spin-Orbit Angle of Exoplanet HAT-P-1b
Comments: 28 total pages, 4 figures, 4 tables, ApJ Accepted
Submitted: 2008-06-10
We present new spectroscopic and photometric observations of the HAT-P-1 planetary system. Spectra obtained during three transits exhibit the Rossiter-McLaughlin effect, allowing us to measure the angle between the sky projections of the stellar spin axis and orbit normal, \lambda = 3.7 +/- 2.1 degrees. The small value of \lambda for this and other systems suggests that the dominant planet migration mechanism preserves spin-orbit alignment. Using two new transit light curves, we refine the transit ephemeris and reduce the uncertainty in the orbital period by an order of magnitude. We find a upper limit on the orbital eccentricity of 0.067, with 99% confidence, by combining our new radial-velocity measurements with those obtained previously.
[152]  oai:arXiv.org:0711.4367  [pdf] - 7457
Retired A Stars and Their Companions II: Jovian planets orbiting kappa Coronae Borealis and HD167042
Comments: 8 pages, 3 figures, 4 tables, ApJ Accepted
Submitted: 2007-11-28
We report precise Doppler measurements of two evolved stars, kappa CrB (HD142091) and HD 167042, obtained at Lick Observatory as part of our search for planets orbiting intermediate-mass subgiants. Periodic variations in the radial velocities of both stars reveal the presence of substellar orbital companions. These two stars are notably massive with stellar masses of 1.80 Msun and 1.64 Msun, indicating that they are former A-type dwarfs that have evolved off of the main sequence and are now K-type subgiants. The planet orbiting kappa CrB has a minimum mass Msini = 1.8 Mjup, eccentricity e = 0.146 and a 1208 day period, corresponding to a semimajor axis of 2.7 AU. The planet around HD167042 has a minimum mass Msini = 1.7 Mjup and a 412.6 day orbit, corresponding to a semimajor axis of 1.3 AU. The eccentricity of HD167042b is consistent with circular (e = 0.027+/-0.04), adding to the rare class of known exoplanets in long-period, circular orbits similar to the Solar System gas giants. Like all of the planets previously discovered around evolved A stars, kappa CrBb and HD167042b orbit beyond 0.8 AU.
[153]  oai:arXiv.org:0711.3593  [pdf] - 7305
The WFCAM Science Archive
Comments: 28 pages, 18 figures; accepted for publication in MNRAS (2007 November 8)
Submitted: 2007-11-22
We describe the WFCAM Science Archive (WSA), which is the primary point of access for users of data from the wide-field infrared camera WFCAM on the United Kingdom Infrared Telescope (UKIRT), especially science catalogue products from the UKIRT Infrared Deep Sky Survey (UKIDSS). We describe the database design with emphasis on those aspects of the system that enable users to fully exploit the survey datasets in a variety of different ways. We give details of the database-driven curation applications that take data from the standard nightly pipeline-processed and calibrated files for the production of science-ready survey datasets. We describe the fundamentals of querying relational databases with a set of astronomy usage examples, and illustrate the results.
[154]  oai:arXiv.org:0707.1908  [pdf] - 3042
The Transit Light Curve Project. VII. The Not-So-Bloated Exoplanet HAT-P-1b
Comments: To appear in AJ [19pg, 3 figures]. New co-author added. Minor revisions to match published version
Submitted: 2007-07-12, last modified: 2007-08-20
We present photometry of the G0 star HAT-P-1 during six transits of its close-in giant planet, and we refine the estimates of the system parameters. Relative to Jupiter's properties, HAT-P-1b is 1.20 +/- 0.05 times larger and its surface gravity is 2.7 +/- 0.2 times weaker. Although it remains the case that HAT-P-1b is among the least dense of the known sample of transiting exoplanets, its properties are in accord with previously published models of strongly irradiated, coreless, solar-composition giant planets. The times of the transits have a typical accuracy of 1 min and do not depart significantly from a constant period.
[155]  oai:arXiv.org:astro-ph/0607282  [pdf] - 83485
Local hydrodynamics and operator stability of Keplerian accretion disks
Comments: Old (October 2000) manuscript submitted to ApJL. See discussion in postscript. Posted for archival purposes only
Submitted: 2006-07-13
We discuss non-self-gravitating hydrodynamic disks in the thin disk limit. These systems are stable according to the Rayleigh criterion, and yet there is some evidence that the dissipative and transport processes in these disks are hydrodynamic in nature at least some of the time. We draw on recent work on the hydrodynamics of laboratory shear flows. Such flows are often experimentally unstable even in the absence of a linear instability. The transition to turbulence in these systems, as well as the large linear transient amplification of initial disturbances, may depend upon the non-self-adjoint nature of the differential operator that describes the dynamics of perturbations to the background state. We find that small initial perturbations can produce large growth in accretion disks in the shearing sheet approximation with shearing box boundary conditions, despite the absence of any linear instability. Furthermore, the differential operator that propagates initial conditions forward in time is asymptotically close (as a function of Reynolds number) to possessing growing eigenmodes. The similarity to the dynamics of laboratory shear flows is suggestive that accretion disks might be hydrodynamically unstable despite the lack of any known instability mechanism.
[156]  oai:arXiv.org:astro-ph/0601092  [pdf] - 78950
Resolving the Surfaces of Extrasolar Planets With Secondary Eclipse Light Curves
Comments: AASTeX 5.2, 24 pages, 5 figures, accepted for publication in ApJ; v2: clarifications, updated to version accepted by ApJ; v3: try to reduce spacing
Submitted: 2006-01-04, last modified: 2006-06-28
We present a method that employs the secondary eclipse light curves of transiting extrasolar planets to probe the spatial variation of their thermal emission. This technique permits an observer to resolve the surface of the planet without the need to spatially resolve its central star. We evaluate the feasibility of this technique for the HD 209458 system [..]. We consider two representations of the planetary thermal emission; a simple model parameterized by a sinusoidal dependence on longitude and latitude, as well as the results of a three-dimensional dynamical simulation of the planetary atmosphere previously published by Cooper & Showman. We find that observations of the secondary eclipse light curve are most sensitive to a longitudinal offset in the geometric and photometric centroids of the hemisphere of the planet visible near opposition. To quantify this signal, we define a new parameter, the ``uniform time offset,'' which measures the time lag between the observed secondary eclipse and that predicted by a planet with a uniform surface flux distribution. We compare the predicted amplitude of this parameter for HD 209458 with the precision with which it could be measured with IRAC. We find that IRAC observations at 3.6um a single secondary eclipse should permit sufficient precision to confirm or reject the Cooper & Showman model of the surface flux distribution for this planet. We quantify the signal-to-noise ratio for this offset in the remaining IRAC bands (4.5um, 5.8um, and 8.0um), and find that a modest improvement in photometric precision (as might be realized through observations of several eclipse events) should permit a similarly robust detection.
[157]  oai:arXiv.org:astro-ph/0603038  [pdf] - 80247
Comment Regarding the Functional Form of the Schmidt Law
Comments: 14 pages, 3 figures, in press, New Astronomy. Figs provided in original (png) format as well as ps format for ps/pdf generation
Submitted: 2006-03-02
Star formation rates on the galactic scale are described phenomenologically by two distinct relationships, as emphasized recently by Elmegreen (2002). The first of these is the Schmidt law, which is a power-law relation between the star formation rate and the column density. The other relationship is that there is a cutoff in the gas density below which star formation shuts off. The purpose of this paper is to argue that 1) these two relationships can be accommodated by a single functional form of the Schmidt law, and 2) this functional form is motivated by the hypothesis that star formation is a critical phenomenon, and that as a corollary, 3) the existence of a sharp cutoff may thus be an emergent property of galaxies, as was argued by Seiden (1983), as opposed to the classical view that this cutoff is due to an instability criterion.
[158]  oai:arXiv.org:astro-ph/0602502  [pdf] - 80078
Turbulent Elasticity of the Solar Convective Zone and the Taylor Number Puzzle
Comments: Six pages, no figures, in Solar MHD: Theory and Observations, eds: J. W. Leibacher, H. Uitenbroek, & R. F. Stein, Astron. Soc. Pac. Conf. Ser
Submitted: 2006-02-23
Previous work on the angular momentum balance and meridional circulation of the solar convective zone (SCZ) generally consists either of semi-analytic approaches in which a simple turbulence model is adopted, or full direct numerical simulation (DNS) of the hydrodynamics or magnetohydrodynamics (MHD). In both instances the inclusion of magnetic fields has been troublesome. Also, both approaches have had difficulty reproducing the known angular velocity profile of the SCZ; this is the Taylor Number Puzzle. I discuss preliminary work in which I incorporate magnetic fields into a viscoelastic turbulence model for the SCZ and apply this to the problem of meridional circulation and angular momentum balance. I suggest that such an approach may help solve the Taylor Number Puzzle of the SCZ and bring theoretical predictions for the large-scale motion of the SCZ in line with observations.
[159]  oai:arXiv.org:astro-ph/0509231  [pdf] - 75813
The Colliding Winds of WR146: Seeing the Works
Comments: 4 pages, 3 figures, contribution to "Massive Stars and High-Energy Emission in OB Associations"; JENAM 2005, held in Liege (Belgium)
Submitted: 2005-09-08
WR146 is a WC6+O8 colliding-wind binary (CWB) system with thermal emission from the stellar winds of the two stars, and bright non-thermal emission from the wind-collision region (WCR) where the winds collide. We present high resolution radio observations from 1.4 to 43 GHz that give one of the best quality radio spectra of any CWB to date. Observations at 22 GHz now span 8 years, and reveal the proper motion of the system, allowing comparison of multi-epoch data. VLBI observations show the location of the WCR relative to the stellar components, from which the wind momentum ratio can be shown to be 0.06+/-0.15. The radio spectrum and the spatial distribution of emission are modelled, and we determine the contribution of both stellar winds and the WCR to the observed emission. We show that our current models fail to account for the high frequency spectrum of WR146, and also produce too much emission far from the stagnation point of the wind collision.
[160]  oai:arXiv.org:astro-ph/0508157  [pdf] - 75046
Are WC9 Wolf-Rayet stars in colliding-wind binaries?
Comments: 4 pages, 3 figures, contribution to "Massive Stars and High-Energy Emission in OB Associations"; JENAM 2005, held in Liege (Belgium)
Submitted: 2005-08-05
We present results from a spectroscopic search for massive companions to dust-making Galactic WC9 stars as a step to testing the paradigm that dust formation in these systems requires colliding winds to produce over densities. We find evidence for OB companions to the WC9 stars WR 59 and WR 65, but not WR 121 or WR 117. We identify lines of N III-V and possibly N II in the spectrum of WR 88, one of the few Galactic WC9 stars which do not make circumstellar dust, and suggest that WR 88 is a transitional WN-WC9 object and less evolved than the other WC9 stars. On the other hand, the possible identification of a strong emission line at 4176A in the spectrum of WR 117 with Ne I suggests that this star is more evolved than other WC9 stars studied.
[161]  oai:arXiv.org:astro-ph/0506184  [pdf] - 73584
Three routes to jet collimation by the Balbus-Hawley magnetorotational instability
Comments: 13 pages; in press, MNRAS
Submitted: 2005-06-09
Three completely different lines of work have recently led to the conclusion that the magnetorotational instability (MRI) may create a hoop-stress that collimates jets. One argument is based upon consideration that magnetohydrodynamic turbulence, in general, and turbulence driven by the MRI, in particular, is more nearly viscoelastic than it is viscous. Another argument is based upon the dispersion relation for the MRI in the context of 1D simulations of core collapse. Yet a third argument rests in the results of direct numerical magnetohydrodynamic simulations of collapsars and thick accretion flows. I elaborate on my previous work regarding the first argument above and I briefly discuss how these three sets of results are all related. I also discuss the different roles played by the magnetic tension and the magnetic pressure within the context of this work. I point out that this leads to consideration of the normal stress difference between the hoop stress and the radial stress, in preference to a focus on just the hoop stress itself. Additionally, I argue that simulations of thick accretion flows and collapsars are not self-consistent if they include a phenomenological model for an MRI-induced viscous stress but disregard these other MRI-induced stress components. I comment briefly on the RHESSI observation of polarization in the gamma-ray burst GRB0212206. I argue that this polarization is consistent with a tangled field, and does not require a large-scale organized field. Finally I suggest that the role of magnetic fields in creating jets, as described here, should be understood not to work within the confines of magnetocentrifugal models of jets, but rather as an alternative to them.
[162]  oai:arXiv.org:astro-ph/0411440  [pdf] - 69039
Mapping WR140 from infrared spectroscopy and imaging
Comments: 5 pages, 6 figures, contribution to "Massive Stars in Interacting Binaries", eds. N. St-Louis and A. F. J. Moffat, ASP Conference Series
Submitted: 2004-11-16
Observations of the 1.083-micron He I line in WR140 (HD 193793) show P-Cygni profiles which varied as the binary system went through periastron passage. A sub-peak appeared on the normally flat-topped emission component and then moved across the profile consistent with its formation in the wind-collision region. Variation of the absorption component provided constraints on the opening angle (theta) of the wind-collision region. Infrared (2-10 micron) images observed with a variety of instruments in 2001-04 resolve the dust cloud formed in 2001, and show it to be expanding at a constant rate. Owing to the high eccentricity of the binary orbit, the dust is spread around the orbital plane in a 'splash' and we compare the dust images with the orientation of the orbit.
[163]  oai:arXiv.org:astro-ph/0409085  [pdf] - 67166
Multipixel characterization of imaging CZT detectors for hard X-ray imaging and spectroscopy
Comments: To appears in the SPIE 2004 proceedings (5540: Hard X-ray and gamma-ray detector physics V)
Submitted: 2004-09-03
We report our in-depth study of Cd-Zn-Te (CZT) crystals to determine an optimum pixel and guard band configuration for Hard X-ray imaging and spectroscopy. We tested 20x20x5mm crystals with 8x8 pixels on a 2.46mm pitch. We have studied different types of cathode / anode contacts and different pixel pad sizes. We present the measurements of leakage current as well as spectral response for each pixel. Our I-V measurement setup is custom designed to allow automated measurements of the I-V curves sequentially for all 64 pixels, whereas the radiation properties measurement setup allows for interchangeable crystals with the same XAIM3.2 ASIC readout from IDEAS. We have tested multiple crystals of each type, and each crystal in different positions to measure the variation between individual crystals and variation among the ASIC channels. We also compare the same crystals with and without a grounded guard band deposited on the crystal side walls vs. a floating guard band and compare results to simulations. This study was carried out to find the optimum CZT crystal configuration for prototype detectors for the proposed Black-Hole Finder mission, EXIST.
[164]  oai:arXiv.org:astro-ph/0404551  [pdf] - 64459
WR20a: a massive cornerstone binary system comprising two extreme early-type stars
Comments: 5 pages, 3 figures, accepted by A&A Letters
Submitted: 2004-04-28
We analyse spectroscopic observations of WR20a revealing that this star is a massive early-type binary system with a most probable orbital period of \sim 3.675 days. Our spectra indicate that both components are most likely of WN6ha or O3If^*/WN6ha spectral type. The orbital solution for a period of 3.675 days yields extremely large minimum masses of 70.7 \pm 4.0 and 68.8 \pm 3.8 M_{\odot} for the two stars. These properties make WR20a a cornerstone system for the study of massive star evolution.
[165]  oai:arXiv.org:astro-ph/0403717  [pdf] - 63908
A Preliminary Observational Search for Circumbinary Disks Around Cataclysmic Variables
Comments: accepted for publication in the Astronomical Journal (July 2004)
Submitted: 2004-03-31
Circumbinary (CB) disks have been proposed as a mechanism to extract orbital angular momentum from cataclysmic variables (CVs) during their evolution. As proposed by Taam & Spruit, these disks extend outwards to several a.u. and should be detected observationally via their infrared flux or by absorption lines in the ultraviolet spectra of the CV. We have made use of archival HST/STIS spectra as well as our own near-IR imaging to search for observational evidence of such CB disks in seven CVs. Based on the null result, we place an upper limit on the column density of the disk of N_H~10^17 cm^-2.
[166]  oai:arXiv.org:astro-ph/0403639  [pdf] - 63830
Near IR Spectroscopic monitoring of WR 140 during the 2001 periastron passage
Comments: Accepted for publication in MNRAS - 13 Pages, 10 figures
Submitted: 2004-03-27
We present new spectra of WR 140 (HD 193793) in the JHK bands with some covering the 1.083-micron He I emission line at higher resolution, observed between 2000 October and 2003 May to cover its 2001 periastron passage. The WC7 + O4-5 spectroscopic binary WR 140 is the prototype of colliding-wind, episodic dust-making Wolf-Rayet systems which also show strong variations in radio and X-ray emission. The JHK spectra showed changes in continuum and in the equivalent widths of the WC emission lines, consistent with formation of dust starting between 2001 January 3 and March 26 (orbital phases 0.989 and 0.017) and its subsequent fading and cooling. The 1.083-micron He I line has a P-Cygni profile which showed variations in both absorption and emission components as WR 140 went through periastron passage. The variation of the absorption component yielded tight constraints on the geometry of the wind-collision region, giving theta = 50 +/- 8 degrees for the opening semi-angle of the interaction `cone', indicating a wind-momentum ratio of the O to the WR star=0.1, about three times larger than previously believed. As the system approached periastron, the emission component showed the appearance of a significant sub-peak, movement of which across the profile was seen to be consistent with its formation in wind material flowing along the contact discontinuity between the two stellar winds and the changing orientation of the colliding wind region. The flux carried in the sub-peak exceeded the X-ray fluxes measured at previous periastron passages. This additional source of radiative cooling of the shock-heated gas probably causes it to depart from being adiabatic around periastron passage, thereby accounting for the departure of the X-ray flux from its previously expected $1/d$-dependency.
[167]  oai:arXiv.org:astro-ph/0312513  [pdf] - 61726
Turbulent Magnetohydrodynamic Jet Collimation and Thermal Driving
Comments: To appear in proceedings of: 3-D Signatures in Stellar Explosions (A workshop honoring J. Craig Wheeler's 60th birthday, June 10-13, 2003, Austin, TX, USA)
Submitted: 2003-12-19, last modified: 2003-12-20
We have argued that magnetohydrodynamic (MHD) turbulence in an accretion disk naturally produces hoop-stresses, and that in a geometrically-thick flow these stresses could both drive and collimate an outflow. We based this argument on an analogy of turbulent MHD fluids to viscoelastic fluids, in which azimuthal shear flow creates hoop-stresses that cause a variety of flow phenomena, including the Weissenberg effect in which a fluid climbs a spinning rod. One of the more important differences between the Weissenberg effect and astrophysical jets is the source of power. In our previous analysis, we only considered the power due to the spin-down torque on the central object, and thus found that we could only drive an outflow if the central object were maximally rotating. Here we take into account the energy that is liberated by the accreting matter, and describe a scenario in which this energy couples to the outflow to create a thermodynamic engine.
[168]  oai:arXiv.org:astro-ph/0307186  [pdf] - 57899
Radio emission models of Colliding-Wind Binary Systems
Comments: 18 pages, 18 figures; Accepted by Astronomy and Astrophysics, July 8, 2003
Submitted: 2003-07-09
We present calculations of the spatial and spectral distribution of the radio emission from a wide WR+OB colliding-wind binary system based on high-resolution hydrodynamical simulations and solutions to the radiative transfer equation. We account for both thermal and synchrotron radio emission, free-free absorption in both the unshocked stellar wind envelopes and the shocked gas, synchrotron self-absorption, and the Razin effect. The applicability of these calculations to modelling radio images and spectra of colliding-wind systems is demonstrated with models of the radio emission from the wide WR+OB binary WR147. Its synchrotron spectrum follows a power-law between 5 and 15 GHz but turns down to below this at lower and higher frequencies. We find that while free-free opacity from the circum-binary stellar winds can potentially account for the low-frequency turnover, models that also include a combination of synchrotron self-absorption and Razin effect are favoured. We argue that the high-frequency turn down is a consequence of inverse-Compton cooling. We present our resulting spectra and intensity distributions, along with simulated MERLIN observations of these intensity distributions. From these we argue that the inclination of the WR147 system to the plane of the sky is low. We summarise by considering extensions of the current model that are important for models of the emission from closer colliding wind binaries, in particular the dramatically varying radio emission of WR140.
[169]  oai:arXiv.org:astro-ph/0212556  [pdf] - 53946
Turbulent Magnetohydrodynamic Elasticity: I. Boussinesq-like Approximations for Steady Shear
Comments: (to be submitted to New Astronomy)
Submitted: 2002-12-27
We re-examine the Boussinesq hypothesis of an effective turbulent viscosity within the context of simple closure considerations for models of strong magnetohydrodynamic turbulence. Reynolds-stress and turbulent Maxwell-stress closure models will necessarily introduce a suite of transport coefficients, all of which are to some degree model-dependent. One of the most important coefficients is the relaxation time for the turbulent Maxwell stress, which until recently has been relatively ignored. We discuss this relaxation within the context of magnetohydrodynamic turbulence in steady high Reynolds-number, high magnetic-Reynolds-number shearing flows. The relaxation time for the turbulent Maxwell stress is not limited by the shear time scale, in contrast with Reynolds-stress closure models for purely hydrodynamical turbulence. The anisotropy of the turbulent stress tensor for magnetohydrodynamic turbulence, even for the case of zero mean-field considered here, can therefore not be neglected. This shear-generated anisotropy can be interpreted as being due to an effective turbulent elasticity, in analogy to the Boussinesq turbulent viscosity. We claim that this turbulent elasticity should be important for any astrophysical problem in which the turbulent stress in quasi-steady shear has been treated phenomenologically with an effective viscosity.
[170]  oai:arXiv.org:astro-ph/0206230  [pdf] - 49861
Protostellar Jets Driven by a Disorganized Magnetic Field
Comments: 5 pages. Talk at the IAOC Workshop, "Galactic Star Formation Across the Stellar Mass Spectrum," La Serena, Chile, March 2002; to be published in the ASP Conference Series, edited by J. M. De Buizer
Submitted: 2002-06-13
We have proposed a viscoelastic model of the Maxwell stresses due to the disorganized magnetic field in MRI-driven MHD turbulence. Viscoelastic fluids in the laboratory are known to produce jet-like structures under the action of a rotating sphere. Here we argue that a similar mechanism may help explain jets in protostellar systems. Such jets would be driven not by large-scale organized magnetic fields, but by the mean-field stresses of small-scale tangled magnetic fields.
[171]  oai:arXiv.org:astro-ph/0204212  [pdf] - 48757
High-Resolution X-ray Imaging of the Colliding Wind Shock in WR147
Comments: 13 pages, 6 figures, accepted by A&A
Submitted: 2002-04-12
We analyze new high-resolution Chandra X-ray images of the Wolf-Rayet binary system WR147. This system contains a WN8 star with an early-type companion located 0.6'' to its north, and is the only known early-type binary with a separation on the sky large enough for the wind-wind collision between the stars to currently be resolved at X-ray energies. The 5 ksec Chandra HRC-I image provides the first direct evidence for spatially extended X-ray emission in an early-type binary system. The X-ray emission peaks close to the position of the radio bow shock and north of the WN8 star. A deeper X-ray image is needed to accurately determine the degree of spatial extension, to exactly align the X-ray and optical/radio frames, and to determine whether part of the detected X-ray emission arises in the individual stellar winds. Simulated X-ray images of the wind-wind collision have a FWHM consistent with the data, and maximum likelihood fits suggest that a deeper observation may also constrain the inclination and wind momentum ratio of this system. However, as the WR wind dominates the colliding wind X-ray emission it appears unlikely that the mass-loss rate and the terminal velocity of the companion wind can be separately determined from X-ray observations. We also note an inconsistency between numerical and analytical estimates of the X-ray luminosity ratio of the stronger and weaker wind components, and conclude that the analytical results are in error.
[172]  oai:arXiv.org:astro-ph/0111603  [pdf] - 46369
Viscoelastic Analogy for the Acceleration and Collimation of Astrophysical Jets
Comments: 8 pages, submitted to ApJL 9/25/01, revised 11/30/01
Submitted: 2001-11-30
Jets are ubiquitous in astronomy. It has been conjectured that the existence of jets is intimately connected with the spin of the central object and the viscous angular momentum transport of the inner disk. Bipolar jet-like structures propelled by the viscous torque on a spinning central object are also known in a completely different context, namely the flow in the laboratory of a viscoelastic fluid. On the basis of an analogy of the tangled magnetic field lines of magnetohydrodynamic (MHD) turbulence to the tangled polymers of viscoelastic polymer solutions, we propose a viscoelastic description of the dynamics of highly turbulent conductive fluid. We argue that the same mechanism that forms jets in viscoelastic fluids in the laboratory may be responsible for collimating and powering astrophysical jets by the angular momentum of the central object.
[173]  oai:arXiv.org:astro-ph/0108286  [pdf] - 1339309
The SuperCOSMOS Sky Survey. Paper I: Introduction and Description
Comments: 18 pages, 10 figures; accepted for publication in MNRAS
Submitted: 2001-08-17
In this, the first in a series of three papers concerning the SuperCOSMOS Sky Survey (SSS), we give an introduction and user guide to the survey programme. We briefly describe other wide-field surveys and compare with our own. We give examples of the data, and make a comparison of the accuracies of the various image parameters available with those from the other surveys providing similar data; we show that the SSS database and interface offer advantages over these surveys. Some science applications of the data are also described and some limitations discussed. The series of three papers constitutes a comprehensive description and user guide for the SSS.
[174]  oai:arXiv.org:astro-ph/0104209  [pdf] - 41948
Intermittent Jolts of Galactic UV Radiation: Mutagenetic Effects
Comments: 10 pages, 1 figure, to be published in "Frontiers of Life; 12th Rencontres de Blois", ed. L. M. Celnikier
Submitted: 2001-04-11
We estimate the frequency of intermittent hypermutation events and disruptions of planetary/satellite photochemistry due to ultraviolet radiation from core collapse supernova explosions. Calculations are presented for planetary systems in the local Milky Way, including the important moderating effects of vertical Galactic structure and UV absorption by interstellar dust. The events are particularly frequent for satellites of giant gas planets at \gtrsim 5-10 AU distance from solar-type parent stars, or in the conventional habitable zones for planets orbiting spectral type K and M parent stars, with rates of significant jolts about 10^3 - 10^4 per Gyr. The steep source spectra and existing data on UVA and longer-wavelength radiation damage in terrestrial organisms suggest that the mutational effects may operate even on planets with ozone shields. We argue that the mutation doubling dose for UV radiation should be much smaller than the mean lethal dose, using terrestrial prokaryotic organisms as our model, and that jolts may lead to important real-time evolutionary episodes if the jolt durations are longer than about a week, corresponding to several hundred generation times, or much less if the equivalent of mutator genes exist in extraterrestrial organisms. Longer-term phylogenetic effects are likely if atmospheric photochemical disturbances lead to niche creation or destruction in relevant habitats.
[175]  oai:arXiv.org:astro-ph/0003048  [pdf] - 34922
WR146 - observing the OB-type companion
Comments: 9 pages, 5 figures, ftp://fto.drao.nrc.ca/pub/smd/wr146/accepted.ps.gz To be published in Monthly Notices of the Royal Astronomical Society
Submitted: 2000-03-03
We present new radio and optical observations of the colliding-wind system WR146 aimed at understanding the nature of the companion to the Wolf-Rayet star and the collision of their winds. The radio observations reveal emission from three components: the WR stellar wind, the non-thermal wind-wind interaction region and, for the first time, the stellar wind of the OB companion. This provides the unique possibility of determining the mass-loss rate and terminal wind velocity ratios of the two winds, independent of distance. Respectively, these ratios are determined to be 0.20+/-0.06 and 0.56+/-0.17 for the OB-companion star relative to the WR star. A new optical spectrum indicates that the system is more luminous than had been believed previously. We deduce that the ``companion'' cannot be a single, low luminosity O8 star as previously suggested, but is either a high luminosity O8 star, or possibly an O8+WC binary system.