BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Peter K. G. Williams

Peter K. G. Williams appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.IMastro-ph.COastro-ph.HEastro-ph.SRastro-ph.GAastro-ph.EPDigital Libraries

Catalog footprint

What is connected

35works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

35 published item(s)

preprint2022arXiv

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Figure and Figure Caption Extraction for Mixed Raster and Vector PDFs: Digitization of Astronomical Literature with OCR Features

Scientific articles published prior to the "age of digitization" in the late 1990s contain figures which are "trapped" within their scanned pages. While progress to extract figures and their captions has been made, there is currently no robust method for this process. We present a YOLO-based method for use on scanned pages, post-Optical Character Recognition (OCR), which uses both grayscale and OCR-features. When applied to the astrophysics literature holdings of the Astrophysics Data System (ADS), we find F1 scores of 90.9% (92.2%) for figures (figure captions) with the intersection-over-union (IOU) cut-off of 0.9 which is a significant improvement over other state-of-the-art methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The First Short GRB Millimeter Afterglow: The Wide-Angled Jet of the Extremely Energetic SGRB 211106A

We present the discovery of the first millimeter afterglow of a short-duration $γ$-ray burst (SGRB) and the first confirmed afterglow of an SGRB localized by the GUANO system on Swift. Our Atacama Large Millimeter/Sub-millimeter Array (ALMA) detection of SGRB 211106A establishes an origin in a faint host galaxy detected in Hubble Space Telescope (HST) imaging at $0.7\lesssim z\lesssim1.4$. From the lack of a detectable optical afterglow, coupled with the bright millimeter counterpart, we infer a high extinction, $A_{\rm V}\gtrsim2.6$ mag along the line of sight, making this the one of the most highly dust-extincted SGRBs known to date. The millimeter-band light curve captures the passage of the synchrotron peak from the afterglow forward shock and reveals a jet break at $t_{\rm jet}=29.2^{+4.5}_{-4.0}$~days. For a presumed redshift of $z=1$, we infer an opening angle, $θ_{\rm jet}=(15.5\pm1.4)$~degrees, and beaming-corrected kinetic energy of $\log(E_{\rm K}/{\rm erg})=51.8\pm0.3$, making this one of the widest and most energetic SGRB jets known to date. Combining all published millimeter-band upper limits in conjunction with the energetics for a large sample of SGRBs, we find that energetic outflows in high density environments are more likely to have detectable millimeter counterparts. Concerted afterglow searches with ALMA should yield detection fractions of 24-40% on timescales of $\gtrsim2$~days at rates $\approx0.8$-1.6 per year, outpacing the historical discovery rate of SGRB centimeter-band afterglows.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum

We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Galactic Faraday rotation sky 2020

This work gives an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. Observations of extra-Galactic sources in recent years have, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This has culminated in an all-sky data set of 55,190 data points, which is a significant expansion on the 41,330 used in previous works, hence making an updated separation of the Galactic component a promising venture. The increased source density allows us to present our results in a resolution of about $1.3\cdot 10^{-2}\, \mathrm{deg}^2$ ($46.8\,\mathrm{arcmin}^2$), which is a twofold increase compared to previous works. As for previous Faraday rotation sky reconstructions, this work is based on information field theory, a Bayesian inference scheme for field-like quantities which handles noisy and incomplete data. In contrast to previous reconstructions, we find a significantly thinner and pronounced Galactic disc with small-scale structures exceeding values of several thousand $\mathrm{rad}\,\mathrm{m}^{-2}$. The improvements can mainly be attributed to the new catalog of Faraday data, but are also supported by advances in correlation structure modeling within numerical information field theory. We furthermore give a detailed discussion on statistical properties of the Faraday rotation sky and investigate correlations to other data sets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

$K2$ Ultracool Dwarfs Survey. VI. White light superflares observed on an L5 dwarf and flare rates of L dwarfs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array

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σ$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $κ_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $κ_\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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Foreground modelling via Gaussian process regression: an application to HERA data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

IDEAS: Immersive Dome Experiences for Accelerating Science

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Reverse Shock in GRB 160509A

We present the second multi-frequency radio detection of a reverse shock in a $γ$-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 $θ_{\rm jet}\approx4~$deg, yielding a beaming-corrected kinetic energy and radiated $γ$-ray energy of $E_{\rm K}\approx4\times10^{50}~$erg and $E_γ\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 $Γ(t_{\rm dec})\approx330$, and a reverse shock to forward shock fractional magnetic energy density ratio of $R_{\rm B}\equivε_{\rm B,RS}/ε_{\rm B,FS}\approx8$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discovery of an outflow from radio observations of the tidal disruption event ASASSN-14li

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, $ρ(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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

M Dwarf Activity in the Pan-STARRS 1 Medium-Deep Survey: First Catalog and Rotation Periods

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Astropy Problem

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

WISEP J060738.65+242953.4: A Nearby. Pole-On L8 Brown Dwarf with Radio Emission

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Extragalactic Transients in the Era of Wide-Field Radio Surveys. I. Detection Rates and Light Curve Characteristics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Kepler Monitoring of an L Dwarf II. Clouds with Multiyear Lifetimes

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$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kepler Monitoring of an L Dwarf I. The Photometric Period and White Light Flares

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$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quasi-Quiescent Radio Emission from the First Radio-Emitting T Dwarf

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 νL_ν = (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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

ASGARD: A Large Survey for Slow Galactic Radio Transients. I. Overview and First Results

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°, b = +13.5°). 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rapid Development of Interferometric Software Using MIRIAD and Python

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Allen Telescope Array Commensal Observing System

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

VAST: An ASKAP Survey for Variables and Slow Transients

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The 2010 May Flaring Episode of Cygnus X-3 in Radio, X-Rays, and γ-Rays

In 2009, Cygnus X-3 (Cyg X-3) became the first microquasar to be detected in the GeV γ-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 γ-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σ γ-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 γ-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 γ-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 γ-ray emission from Cyg X-3.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Allen Telescope Array Pi GHz Sky Survey II. Daily and Monthly Monitoring for Transients and Variability in the Bootes Field

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σ) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Allen Telescope Array Twenty-centimeter Survey -- A 700-Square-Degree, Multi-Epoch Radio Dataset -- II: Individual Epoch Transient Statistics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Commensal observing with the Allen Telescope array: software command and control

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Primary Beam Shape Calibration from Mosaicked, Interferometric Observations

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 $χ^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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Allen Telescope Array Pi GHz Sky Survey I. Survey Description and Static Catalog Results for the Bootes Field

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Allen Telescope Array Twenty-centimeter Survey - A 690-Square-Degree, 12-Epoch Radio Dataset - I: Catalog and Long-Duration Transient Statistics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Evaluating the Calorimeter Model with Broadband, Continuous Spectra of Starburst Galaxies Observed with the Allen Telescope Array

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Old, Rich, and Eccentric: Two Jovian Planets Orbiting Evolved Metal-Rich Stars

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.

0 citations0 reviewsNo code linkedOpen access